Search
2024 Volume 4
Article Contents
ARTICLE   Open Access    

Genetic selection for growth, wood quality and resin traits of potential Slash pine for multiple industrial uses

  • # Authors contributed equally: Xianyin Ding, Yini Zhang, Jiaming Sun

More Information
  • Received: 29 February 2024
    Revised: 11 May 2024
    Accepted: 04 June 2024
    Published online: 20 June 2024
    Forestry Research  4 Article number: e023 (2024)  |  Cite this article
  • This study aims to understand the genetic basis of key industrial traits in Slash pine (Pinus elliottii Engelm. var. elliottii) to enhance improvement efficiency. Detailed analyses were conducted on inter-family differences, genetic parameters, correlations, and breeding values (BVs) for growth, wood properties, and resin traits of Slash pine planted in Changle Forest Farm of Hangzhou, leading to the identification of elite families. It indicates that growth traits are primarily influenced by environmental effects, while wood properties exhibit a significant impact of genetic effects. The variation in resin traits arises from both genetic and environmental effects. Notably, Beta-pinene exhibits the highest variability and genetic gains among the traits analyzed. The family heritability ranges for growth, wood properties, and resin traits are 0.543−0.794, 0.870−0.885, and 0.285−0.695, respectively. Significant positive correlations are evident between growth and resin traits, while a negative correlation is observed between growth and wood properties. Elite families identified through single-trait and multi-trait combined selection are 8−126 for growth traits, 2−325 and 0−373 for wood properties, and 8−131 for resin traits. The average genetic gains for these elite families are 7.44%, 7.17%, and 8.84%, respectively. These findings provide valuable insights for high-generation breeding of Slash pine and lay a genetic foundation for formulating effective breeding strategies for conifers.
  • 加载中
  • Supplemental Table S1 The information of phenotype data.
    Supplemental Table S2 Estimated breeding value for growth traits.
    Supplemental Table S3 Estimated breeding value for wood properties.
    Supplemental Table S4 Estimated breeding value for resin traits.
  • [1]

    Fengel D, Wegener G. 1984. Wood: chemistry, ultrastructure, reactions. Berlin, Germany: Walter de Gruyter. xiii, 613 pp

    [2]

    Hamilton MG, Freeman JS, Blackburn DP, Downes GM, Pilbeam DJ, et al. 2017. Independent lines of evidence of a genetic relationship between acoustic wave velocity and kraft pulp yield in Eucalyptus globulus. Annals of Forest Science 74:17

    doi: 10.1007/s13595-017-0617-2

    CrossRef   Google Scholar

    [3]

    Silva JC, Borralho NMG, Araújo JA, Vaillancourt RE, Potts BM. 2009. Genetic parameters for growth, wood density and pulp yield in Eucalyptus globulus. Tree Genetics & Genomes 5:291−305

    doi: 10.1007/s11295-008-0174-9

    CrossRef   Google Scholar

    [4]

    Santos GA, Nunes ACP, Resende MDV, Silva LD, Higa A, et al. 2016. Genetic control and genotype-by-environment interaction of wood weight in Eucalyptus clones in the state of Rio Grande do Sul, Brazil. Revista Árvore 40:867−76

    doi: 10.1590/0100-67622016000500010

    CrossRef   Google Scholar

    [5]

    Gao S, Wang X, Wiemann MC, Brashaw BK, Ross RJ, et al. 2017. A critical analysis of methods for rapid and nondestructive determination of wood density in standing trees. Annals of Forest Science 74:27

    doi: 10.1007/s13595-017-0623-4

    CrossRef   Google Scholar

    [6]

    Viana LC, Trugilho PF, Hein PRG, Moreira da Silva JR, Lima JT. 2010. Modelos de calibração e a espectroscopia no infravermelho próximo para predição das propriedades químicas e da densidade básica da madeira de Eucalyptus. Ciência Florestal 20:367−76

    doi: 10.5902/198050981859

    CrossRef   Google Scholar

    [7]

    Zulak KG, Bohlmann J. 2010. Terpenoid biosynthesis and specialized vascular cells of conifer defense. Journal of Integrative Plant Biology 52:86−97

    doi: 10.1111/j.1744-7909.2010.00910.x

    CrossRef   Google Scholar

    [8]

    Kelkar VM, Geils BW, Becker DR, Overby ST, Neary DG. 2006. How to recover more value from small pine trees: essential oils and resins. Biomass and Bioenergy 30:316−20

    doi: 10.1016/j.biombioe.2005.07.009

    CrossRef   Google Scholar

    [9]

    Neis FA, de Costa F, de Araújo AT Jr, Fett JP, Fett-Neto AG. 2019. Multiple industrial uses of non-wood pine products. Industrial Crops and Products 130:248−58

    doi: 10.1016/j.indcrop.2018.12.088

    CrossRef   Google Scholar

    [10]

    Luan Q, Tao X, Diao S, Ding X, Jiang J. 2021. Methods, characteristics, variance, and genetics of pine oleoresin components, and their potential for renewable and sustainable energy. Studies in Natural Products Chemistry 68:221−53.

    doi: 10.1016/B978-0-12-819485-0.00018-9

    CrossRef   Google Scholar

    [11]

    Luan Q, Diao S, Sun H, Ding X, Jiang J. 2022. Prediction and comparisons of turpentine content in Slash pine at different slope positions using near-infrared spectroscopy. Plants 11:914

    doi: 10.3390/plants11070914

    CrossRef   Google Scholar

    [12]

    Ding X, Li Y, Zhang Y, Diao S, Luan Q, et al. 2023. Genetic analysis and elite tree selection of the main resin components of slash pine. Frontiers in Plant Science 14:1079952

    doi: 10.3389/fpls.2023.1079952

    CrossRef   Google Scholar

    [13]

    Diao S, Zhang Y, Luan Q, Ding X, Sun J, et al. 2022. Identification of TPS-d subfamily genes and functional characterization of three monoterpene synthases in Slash pine. Industrial Crops and Products 188:115609

    doi: 10.1016/j.indcrop.2022.115609

    CrossRef   Google Scholar

    [14]

    Serra S, Fuganti C, Brenna E. 2005. Biocatalytic preparation of natural flavours and fragrances. Trends in Biotechnology 23:193−98

    doi: 10.1016/j.tibtech.2005.02.003

    CrossRef   Google Scholar

    [15]

    Cunnighan A. 2012. Pine resin: biology, chemistry and applications. Pine Tapping Recent Advances. Bosto: PCA. pp. 1−8

    Google Scholar

    [16]

    Peralta-Yahya PP, Zhang F, Del Cardayre SB, Keasling JD. 2012. Microbial engineering for the production of advanced biofuels. Nature 488:320−28

    doi: 10.1038/nature11478

    CrossRef   Google Scholar

    [17]

    da Silva Rodrigues-Corrêa KC, de Lima JC, Fett-Neto AG. 2013.Oleoresins from pine: production and industrial uses. In Natural Products, eds Ramawat K, Mérillon JM. Heidelberg: Springer, Berlin. pp. 4037−60. https://doi.org/10.1007/978-3-642-22144-6_175

    [18]

    Lai M, Dong L, Yi M, Sun S, Zhang Y, et al. 2017. Genetic variation, heritability and genotype × environment interactions of resin yield, growth traits and morphologic traits for Pinus elliottii at three progeny trials. Forests 8:409

    doi: 10.3390/f8110409

    CrossRef   Google Scholar

    [19]

    Lai M, Zhang L, Lei L, Liu S, Jia T, et al. 2020. Inheritance of resin yield and main resin components in Pinus elliottii Engelm. at three locations in southern China. Industrial Crops and Products 144:112065

    doi: 10.1016/j.indcrop.2019.112065

    CrossRef   Google Scholar

    [20]

    Williams CG, Megraw RA. 1994. Juvenile-mature relationships for wood density in Pinustaeda. Canadian Journal of Forest Research 24:714−22

    doi: 10.1139/x94-095

    CrossRef   Google Scholar

    [21]

    Mao T. 2007. Genetic analysis and combined selection of the growth and material of superior tree progeny of Masson pine. Thesis. Nanjing Forestry University, China.

    [22]

    Fries A. 1986. Volume growth and wood density of plus tree progenies of Pinus contorta in two Swedish field trials. Scandinavian Journal of Forest Research 1:403−19

    doi: 10.1080/02827588609382433

    CrossRef   Google Scholar

    [23]

    Belonger PJ, McKeand SE, Jett JB. 1996. Genetic and environmental effects on biomass production and wood density in loblolly pine. In Tree Improvement for Sustainable Tropical Forestry, Proceedings of the QFRI-IUFRO Conference, Caloundra, Queensland, Australia. pp. 307–10.

    [24]

    Sun X, Zhang S, Li S, Hou Y. 2005. Combined selection of fine families with multiple traits for pulp wood of Larch japonica. Scientia Silvae Sinicae 41:48−54

    Google Scholar

    [25]

    Wu HX, Powell MB, Yang JL, Ivković M, McRae TA. 2007. Efficiency of early selection for rotation-aged wood quality traits in radiata pine. Annals of Forest Science 64:1−9

    doi: 10.1051/forest:2006082

    CrossRef   Google Scholar

    [26]

    Apiolaza L, Chauhan S, Hayes M, Nakada R, Sharma M, et al. 2013. Selection and breeding for wood quality: a new approach. New Zealand Journal of Forestry 58:33−37

    Google Scholar

    [27]

    O'hehir JF, Nambiar EKS. 2010. Productivity of three successive rotations of Pinus radiata plantations in South Australia over a century. Forest Ecology and Management 259:1857−69

    doi: 10.1016/j.foreco.2009.12.004

    CrossRef   Google Scholar

    [28]

    Li Y, Luan Q, Shen D, Chen B, Jiang J, et al. 2012. Study on genetic variation of resin components among open-pollinated families of slash pine. Forest Research 25:773−79

    doi: 10.13275/j.cnki.lykxyj.2012.06.007

    CrossRef   Google Scholar

    [29]

    Haapanen M. 1996. Impact of family-by-trial interaction on the utility of progeny testing methods for Scots pine. Silvae Genetica 45:130−35

    Google Scholar

    [30]

    Costa e Silva J, Potts BM, Dutkowski GW. 2006. Genotype by environment interaction for growth of Eucalyptus globulus in Australia. Tree Genetics & Genomes 2:61−75

    doi: 10.1007/s11295-005-0025-x

    CrossRef   Google Scholar

    [31]

    Campbell RK. 1992. Genotype × environment interaction: a case study for Douglas-fir in western Oregon. Research Paper. Res. Pap. PNW-RP-455. U.S. Department of Agriculture, Forest Service, Pacific Northwest Research Station, Portland, OR. 21 pp. https://doi.org/10.2737/PNW-RP-455

    [32]

    Lundströmer J, Karlsson B, Berlin M. 2020. Strategies for deployment of reproductive material under supply limitations - a case study of Norway spruce seed sources in Sweden. Scandinavian Journal of Forest Research 35:495−505

    doi: 10.1080/02827581.2020.1833979

    CrossRef   Google Scholar

    [33]

    Sixto H, Salvia J, Barrio M, Ciria MP, Cañellas I. 2011. Genetic variation and genotype-environment interactions in short rotation Populus plantations in southern Europe. New Forests 42:163−77

    doi: 10.1007/s11056-010-9244-6

    CrossRef   Google Scholar

    [34]

    Alexandru AM, Mihai G, Stoica E, Curtu AL. 2023. Multi-trait selection and stability in Norway Spruce (Picea abies) provenance trials in Romania. Forests 14:456

    doi: 10.3390/f14030456

    CrossRef   Google Scholar

    [35]

    Yan W, Frégeau-Reid J. 2008. Breeding line selection based on multiple traits. Crop Science 48:417−23

    doi: 10.2135/cropsci2007.05.0254

    CrossRef   Google Scholar

    [36]

    Smith HF. 1936. A discriminant function for plant selection. Annals of Eugenics 7:240−50

    doi: 10.1111/j.1469-1809.1936.tb02143.x

    CrossRef   Google Scholar

    [37]

    BuenoJ, Vencovsky R. 2000. Efficiency of combined selection over sequential selection in forest tree progeny trials. Silvae Genetica 49:169−73

    Google Scholar

    [38]

    Dong P, Chang J, Huang L, Li C. 2015. Application of multi-factor comprehensive appraisal in maize breeding. Agricultural Science & Technology 16:1614−1616, 1638

    doi: 10.3969/j.issn.1009-4229.2015.08.012

    CrossRef   Google Scholar

    [39]

    Jia Q, Zhang H, Zhang L. 2016. Variation analysis of hybrid larch families and superior famlies selection. Journal of Northeast Forestry University 44:1−7

    doi: 10.3969/j.issn.1000-5382.2016.04.001

    CrossRef   Google Scholar

    [40]

    Dong L, Li T, Huang W, Wang B, Xu L, et al. 2021. Screening and comprehensive evaluation of superior strains of Zhejiang safflower camellia oleifera. Journal of Central South Forestry University 41:35−45

    Google Scholar

    [41]

    Ding X, Diao S, Luan Q, Wu HX, Zhang Y, et al. 2022. A transcriptome-based association study of growth, wood quality, and oleoresin traits in a slash pine breeding population. PLoS Genetics 18:e1010017

    doi: 10.1371/journal.pgen.1010017

    CrossRef   Google Scholar

    [42]

    Ding X, Tao X, Diao S, Luan Q, Jiang J. 2020. Estimation of wood basic density in a Pinus elliottii stand using Pilodyn and Resistograph measurements. Journal of Nanjing Forestry University (Natural Sciences Edition) 44:142−48

    doi: 10.3969/j.issn.1000-2006.201906026

    CrossRef   Google Scholar

    [43]

    Tham VTH, Inagaki T, Tsuchikawa S. 2018. A novel combined application of capacitive method and near-infrared spectroscopy for predicting the density and moisture content of solid wood. Wood Science and Technology 52:115−29

    doi: 10.1007/s00226-017-0974-x

    CrossRef   Google Scholar

    [44]

    Zhang S, Jiang J, Xu Y, Luan Q. 2017. Study on the modulus of elasticity non-destructive evaluation technique of slash pine standing tree. Forestry Scientific Research 30:75−80

    doi: 10.13275/j.cnki.lykxyj.2017.01.011

    CrossRef   Google Scholar

    [45]

    Li Y, Jiang J, Luan Q. 2012. Determination and genetic analysis of lipid-producing power, resin density and oleresin content of slash pine. Journal of Beijing Forestry University 34:48−51

    doi: 10.13332/j.1000-1522.2012.04.022

    CrossRef   Google Scholar

    [46]

    Cuevas J, Montesinos-López O, Juliana P, Guzmán C, Pérez-Rodríguez P, et al. 2019. Deep kernel for genomic and near infrared predictions in multi-environment breeding trials. G3 Genes|Genomes|Genetics 9:2913−24

    doi: 10.1534/g3.119.400493

    CrossRef   Google Scholar

    [47]

    R Core Team. 2014. R: a language and environment for statistical computing.

    [48]

    Vazquez AI, Bates DM, Rosa GJM, Gianola D, Weigel KA. 2013. Technical note: an R package for fitting generalized linear mixed models in animal breeding. Journal of Animal Science 88:497−504

    doi: 10.2527/jas.2009-1952

    CrossRef   Google Scholar

    [49]

    Covarrubias-Pazaran G. 2016. Genome-assisted prediction of quantitative traits using the R package sommer. PLoS One 11:e0156744

    doi: 10.1371/journal.pone.0156744

    CrossRef   Google Scholar

    [50]

    Revelle W. 2020. How to: use the psych package for factor analysis and data reduction. Northwestern University, Department of Psychology, Evanston, IL, USA. 95 pp. https://personality-project.org/r/psych/HowTo/factor.pdf

    [51]

    Wickham H. 2011. ggplot2. WIREs Computational Statistics 3:180−85

    doi: 10.1002/wics.147

    CrossRef   Google Scholar

    [52]

    Yuan C, Zhang Z, Jin G, Zheng Y, Zhou Z, et al. 2021. Genetic parameters and genotype by environment interactions influencing growth and productivity in Masson pine in east and central China. Forest Ecology and Management 487:118991

    doi: 10.1016/j.foreco.2021.118991

    CrossRef   Google Scholar

    [53]

    Liu Q, Zhou Z, Fan H, Liu Y. 2013. Genetic variation and correlation among resin yield, growth, and morphologic traits of Pinus massoniana. Silvae Genetica 62:38−43

    doi: 10.1515/sg-2013-0005

    CrossRef   Google Scholar

    [54]

    Jia Q, Liu G, Zhao J, Li K, Sun W. 2022. Variation analyses of growth traits in half-sib families of Korean pine and superior families selection. Journal of Nanjing Forestry University (Natural Sciences Edition) 46:109−16

    doi: 10.12302/j.issn.1000-2006.202107040

    CrossRef   Google Scholar

    [55]

    Baltunis BS, Wu HX, Powell MB. 2007. Inheritance of density, microfibril angle, and modulus of elasticity in juvenile wood of Pinus radiata at two locations in Australia. Canadian Journal of Forest Research 37:2164−74

    doi: 10.1139/X07-061

    CrossRef   Google Scholar

    [56]

    Lei L, Pan X, Zhang L, Ai Q, Li B, et al. 2015. Genetic variation and comprehensive selection of turpentine composition in high-yielding slash pine (Pinus elliottii). Forest Research 28:804−09

    doi: 10.3969/j.issn.1001-1498.2015.06.007

    CrossRef   Google Scholar

    [57]

    Li Y, Luan Q, Jiang J. 2012. Study on genetic variation of rosin components in open-pollinated slash pine families. Forestry Science Research 25:773−79

    Google Scholar

    [58]

    Li Y, Ding X, Jiang J, Luan Q. 2020. Inheritance and correlation analysis of pulpwood properties, wood density, and growth traits of slash pine. Forests 11:493

    doi: 10.3390/f11050493

    CrossRef   Google Scholar

    [59]

    Lenz PRN, Nadeau S, Mottet MJ, Perron M, Isabel N, et al. 2020. Multi-trait genomic selection for weevil resistance, growth, and wood quality in Norway spruce. Evolutionary Applications 13:76−94

    doi: 10.1111/eva.12823

    CrossRef   Google Scholar

    [60]

    Hannrup B, Cahalan C, Chantre G, Grabner M, Karlsson B, et al. 2004. Genetic parameters of growth and wood quality traits in Picea abies. Scandinavian Journal of Forest Research 19:14−29

    doi: 10.1080/02827580310019536

    CrossRef   Google Scholar

    [61]

    Johnson LPV. 1942. Studies on the relation of growth rate to wood quality in Populus hybrids. Canadian Journal of Research 20:28−40

    doi: 10.1139/cjr42c-003

    CrossRef   Google Scholar

    [62]

    Park YS, Weng Y, Mansfield SD. 2012. Genetic effects on wood quality traits of plantation-grown white spruce (Picea glauca) and their relationships with growth. Tree Genetics & Genomes 8:303−11

    doi: 10.1007/s11295-011-0441-z

    CrossRef   Google Scholar

    [63]

    Gaspar MJ, Lousada JL, Rodrigues JC, Aguiar A, Almeida MH. 2009. Does selecting for improved growth affect wood quality of Pinus pinaster in Portugal. Forest Ecology and Management 258:115−21

    doi: 10.1016/j.foreco.2009.03.046

    CrossRef   Google Scholar

    [64]

    Hood S, Sala A. 2015. Ponderosa pine resin defenses and growth: metrics matter. Tree Physiology 35:1223−35

    doi: 10.1093/treephys/tpv098

    CrossRef   Google Scholar

    [65]

    López-Álvarez Ó, Zas R, Marey-Perez M. 2023. Resin tapping: a review of the main factors modulating pine resin yield. Industrial Crops and Products 202:117105

    doi: 10.1016/j.indcrop.2023.117105

    CrossRef   Google Scholar

    [66]

    Hazel LN. 1943. The genetic basis for constructing selection indexes. Genetics 28:476−90

    doi: 10.1093/genetics/28.6.476

    CrossRef   Google Scholar

    [67]

    Walsh B, Lynch M. 2018. Evolution and selection of quantitative traits. Oxford: Oxford University Press. https://doi.org/10.1093/oso/9780198830870.001.0001

    [68]

    Amaya A, Martínez R, Cerón-Muñoz M. 2021. Selection indexes using principal component analysis for reproductive, beef and milk traits in Simmental cattle. Tropical Animal Health and Production 53:378

    doi: 10.1007/s11250-021-02815-y

    CrossRef   Google Scholar

    [69]

    Fries A. 2012. Genetic parameters, genetic gain and correlated responses in growth, fibre dimensions and wood density in a Scots pine breeding population. Annals of Forest Science 69:783−94

    doi: 10.1007/s13595-012-0202-7

    CrossRef   Google Scholar

  • Cite this article

    Ding X, Zhang Y, Sun J, Tan Z, Huang Q, et al. 2024. Genetic selection for growth, wood quality and resin traits of potential Slash pine for multiple industrial uses. Forestry Research 4: e023 doi: 10.48130/forres-0024-0020
    Ding X, Zhang Y, Sun J, Tan Z, Huang Q, et al. 2024. Genetic selection for growth, wood quality and resin traits of potential Slash pine for multiple industrial uses. Forestry Research 4: e023 doi: 10.48130/forres-0024-0020

Figures(5)  /  Tables(5)

Article Metrics

Article views(1655) PDF downloads(326)

ARTICLE   Open Access    

Genetic selection for growth, wood quality and resin traits of potential Slash pine for multiple industrial uses

Forestry Research  4 Article number: e023  (2024)  |  Cite this article

Abstract: This study aims to understand the genetic basis of key industrial traits in Slash pine (Pinus elliottii Engelm. var. elliottii) to enhance improvement efficiency. Detailed analyses were conducted on inter-family differences, genetic parameters, correlations, and breeding values (BVs) for growth, wood properties, and resin traits of Slash pine planted in Changle Forest Farm of Hangzhou, leading to the identification of elite families. It indicates that growth traits are primarily influenced by environmental effects, while wood properties exhibit a significant impact of genetic effects. The variation in resin traits arises from both genetic and environmental effects. Notably, Beta-pinene exhibits the highest variability and genetic gains among the traits analyzed. The family heritability ranges for growth, wood properties, and resin traits are 0.543−0.794, 0.870−0.885, and 0.285−0.695, respectively. Significant positive correlations are evident between growth and resin traits, while a negative correlation is observed between growth and wood properties. Elite families identified through single-trait and multi-trait combined selection are 8−126 for growth traits, 2−325 and 0−373 for wood properties, and 8−131 for resin traits. The average genetic gains for these elite families are 7.44%, 7.17%, and 8.84%, respectively. These findings provide valuable insights for high-generation breeding of Slash pine and lay a genetic foundation for formulating effective breeding strategies for conifers.

    • Wood, known for its renewable and eco-friendly properties, serves diverse purposes in construction, furniture, and pulp production[1]. Key growth parameters like diameter at breast height (DBH) and tree height (Ht), alongside wood properties such as wood density (WD) and modulus of elasticity (MOE), are essential for breeding programs targeting industrial production[2,3]. Traditional assessment methods for these traits are costly, time-intensive, and destructive, posing challenges in forestry research[4]. However, non-destructive wood quality assessment technologies like Pilodyn and Resistograph offer a promising solution, enabling rapid and accurate target detection in a cost-effective and high-throughput manner[5,6]. Resin, a viscous fluid mixture of secondary metabolites comprises turpentine rich in monoterpenes and sesquiterpenes, and rosin predominantly composed of diterpenes[7]. Turpentine finds applications in fragrance, pharmaceuticals, and food additives, while rosin is extensively used in agrochemicals and adhesives[8,9]. Resin derivatives have diversified applications in biodegradable batteries, green plastics, and petroleum alternatives[9,10]. Each resin component serves distinct purposes in industrial applications[11,12]. For example, pinene, the predominant turpentine, is utilized in natural insecticides, pharmaceuticals, and food additives[1316]. Among pinene isomers, Alpha-pinene is the most abundant, while Beta-pinene holds significant value in the chemical industry[17]. Slash pine (Pinus elliottii Engelm. var. elliottii) resin, characterized by its high Beta-pinene content is of considerable industrial value[18,19].

      Genetic variation is crucial for genetic improvement, providing insights into trait breeding potential[12]. In forestry genetic improvement, breeders often aim to enhance multiple traits simultaneously. Therefore, investigating the relationship between growth traits and wood properties is crucial. Studies have indicated a positive genetic correlation between growth traits and wood properties in P. taeda[20] and P. massoniana[21]. However, the correlation between growth traits and WD in P. contorta[22] is not significant, contrasting with other research suggesting a negative correlation between growth trait and wood properties in P. taeda[23], Larix kaempferi[24], and P. radiata[25]. Consequently, achieving rapid wood maturation through genetic improvement for growth traits may potentially have adverse effects on wood quality[26,27]. Additionally, research on resin traits has revealed significant negative correlations between Alpha- and Beta-pinene. Furthermore, resin yield and resin basic density exhibit a highly significant positive correlation, while turpentine content shows a highly significant negative correlation with resin basic density[28]. Understanding these connections is crucial for implementing multi-trait combined selection, optimizing breeding strategies, and achieving balanced improvements in growth-wood-resin traits[18].

      However, economically valuable traits in tree breeding programs typically exhibit polygenic variation and are strongly influenced by environmental factors. Almost all commercially important tree species have reported significant Genotype-by-Environment interactions (G × E), including Scots pine[29], Eucalyptus[30], Douglas fir[31], Norway spruce[32], and Poplar[33]. Multi-trait combined selection has the potential to mitigate the impact of the environment on genotype stability and expedite the breeding process[34]. Nonetheless, addressing linkage disequilibrium (LD) while simultaneously selecting for multiple desirable traits has led to various methods for multi-trait combined selection[35], such as index selection[36], sequential selection[37], and multi-factor comprehensive evaluation[38]. Principal component analysis (PCA) is one such method, revealing the intrinsic characteristics of test materials by extracting factors from the multifactorial relationships between components and performance, which has demonstrated efficacy in Larix kaempferi[39] and Camellia chekiangoleosa[40].

      Slash pine plays a crucial role in providing wood and resin material, with genetic improvement efforts primarily targeting growth rates, wood quality, and oleresin yield (OY)[41]. This study aims to: (1) investigate the varied G × E effects on growth, wood properties, and resin traits; (2) estimate the heritability, correlations, and breeding values of these traits; and (3) employ multi-trait combined selection to identify elite families. These findings may enhance the efficiency of genetic improvement for key industrial traits in Slash pine and accelerate its high-generation breeding process.

    • The Slash pine clonal seed orchard (located at 30°20' N, 119°50' E) comprises 300 elite clone genotypes selected from Slash pine provenance test forests in seven provinces: Guangdong, Fujian, Jiangxi, Jiangsu, Hunan, Hubei, and Zhejiang. These clones, all exhibiting elite growth performance and unrelated to each other, were propagated in April 1976 by grafting robust branches from the crown of selected elite trees onto 2-year-old Slash pine rootstocks using the method of cambial layer grafting.

      In 1992, seeds from the clonal families of the orchard were collected, and seedlings were raised in 1993. In 1994, a clonal stand test forest was established at coordinates 30°27' N, 119°49' E. The test forest comprises 33 families, arranged in a completely randomized design with six offspring in single-row plots and six replications (blocks), with trees spaced 2 m apart within rows and rows spaced 3 m apart. The local annual average precipitation is 1,480 mm, and the annual average temperature is 17.0 °C. The planting site is characterized by gently sloping hills with yellow-red soil of moderate fertility. Thinning was carried out in the test forest after 10 years of establishment. The test forest had never been tapped for resin before the commencement of this experiment.

      In July 2019, phenotypic surveys were conducted on four blocks of the test forest. The three trees with the largest DBH were measured in each plot, and trees with obvious damage were excluded. This resulted in a total of 240 individual phenotype data collected (Supplemental Table S1).

    • In this study, the phenotypic traits and survey methods for the growth traits of Slash pine are as follows:

      (1) DBH: Measured with a caliper at a height of 1.3 m from the ground on the trunk.

      (2) Ht and height under the crown (Huc): Measured using a height measuring pole.

      (3) Crown width: Diameter of the crown measured in both east-west and north-south directions, with the average result taken and the area calculated.

      (4) Average annual ring width (ARW): ARW was measured using the Resistograph-IML-RESIPD500 (IML-RESI GmbH, Taufkirchen, Germany). The drill was inserted vertically at the tree's DBH, recording depth and amplitude for each 0.1 mm. Peaks and valleys correspond to annual growth cycles. Data selection excluded bark-induced irregularities, focusing on distinct resistance peaks. The ARW and resistance value (Mountain peak, MP) were calculated. Due to needle height, the first-year ring is often unclear, starting ring counting from the outermost. Each tree's data covered 24 years for RW and yearly MP. ARW and AMP were recorded, followed by data analysis using PDToolsPro software. The detailed operational method has been previously described in our earlier research[42].

    • (1) WD: Using a growth increment borer with a diameter of 5 mm, cores were extracted from the trees at breast height in the north-south direction, ensuring that the borehole penetrated through the pith of the tree. Basic density ($ \rho $) was determined using the saturated moisture content method[43]:

      $ \rho =\dfrac{1}{\dfrac{{\mathrm{M}}_{\mathrm{w}}}{{\mathrm{M}}_{\mathrm{d}}}-0.3464} $ (1)

      In Eqn (1), $ \rho $ is the basic density of wood (g·cm−3), $ {\mathrm{M}}_{\mathrm{w}} $ is the mass of the wood core when saturated with water, $ {\mathrm{M}}_{\mathrm{d}} $ is the mass of wood core when completely dry.

      (2) Determination of stress wave velocity $ v $ and MOE: The Hitman ST300 (Fibre-gen, Christchurch, New Zealand) was utilized to measure the $ v $ of the samples. The specific method is detailed in the research findings of Zhang et al.[44]. Subsequently, the MOE of the samples was calculated using Eqn (2) based on the $ \rho $ and $ v $.

      $ MOE=\rho {v}^{2} $ (2)

      (3) Pilodyn (Pd) detection: Using Pilodyn (6J, PROCEQ, Switzerland), cores were extracted from the lower 2 cm of the stem of individual trees in both the south and north directions. Two Pd resistance values (Ps and Pn) were measured, and only values with a difference of no more than 2 mm were considered valid. Otherwise, the two values with the smallest difference were selected. The average resistance values in the south and north directions were denoted as Pd.

      (4) The AMP data acquisition method is as described in the section 'Growth traits phenotyping'.

    • (1) OY: A specialized plastic tube with a diameter of 1.8 mm and a volume of 15 mL is fixed in the borehole on the sunny side of the trunk. After waiting for 24 h, the resin is collected, and then its yield is measured using a balance. The detailed installation method is described in Li et al.[45].

      (2) Relative amount of resin components: Gas chromatography experiments were carried out with a GC 6890 gas chromatograph (Agilent 5975B, Santa Clara, CA, USA) equipped with a DB-5MS capillary column cross-linked with 5% Ph Me silicone (60 m/0.25 mm/0.25 mm) and coupled with a Hewlett Packard GC 5975B mass spectrometer. The qualitative and quantitative analysis of resin composition with the chromatographic conditions was as follows: GC: 0.05 g of pine resin was dissolved in 0.5 mL of ethyl alcohol containing 50 μL tetramethylammonium hydroxide. The initial column temperature was 60 °C, held for 2 min, increased at 8 °C·min−1 to 80 °C, and reached a maximum of 280 °C at a rate of 2 °C·min−1 for 5 min. The helium gas flow was set at 1 mL·min−1. The temperature of the injector was 260 °C, and the volume was 1 μL with a 1/50 split ratio. Mass spectra were recorded under electron impact ionization at an electron energy of 70 eV in the range from m/z 30 to 600 along with solvent delay for 3 min.

      Resin compositions were identified by matching experimental fragmentation patterns in mass spectra with the NIST08 database through the data processing system of Agilent Chem Station and then compared with the relevant literature. Monoterpene content was determined by isobutylbenzene content, and diterpene content was determined by heptadecanoic acid content. The resin component contents were calculated by comparing their peak areas. We then selected four components believed to have broad industrial utility for subsequent analysis, including two monoterpenes (Alpha- and Beta-pinene) and two diterpenes (Abietic and Levoprimaric acid).

    • The linear model for each observation of the resin component trait yijk in each tree was shown as follows:

      $ \mathit{y} _{ \mathit{ijk} } = \mu+{f} _{ \mathit{i} } \mathit{+b} _{ \mathit{j} } \mathit{+fb} _{ \mathit{ij} } \mathit{+e} _{ \mathit{ijk} } $ (3)

      Where the observed values (yijk) for individual plants within each family (i), block (j), and specific plant (k) are considered, The μ represents the average value across all observations. Family effects are denoted by fi. Similarly, bj refers to block effects. The interaction effect between family and block is represented by fbij. Environmental error effects (eijk) encapsulate random variations and unforeseen factors affecting individual plant observations.

      The genetic variation analysis employed the Restricted Maximum Likelihood (REML) method to fit the Generalized Linear Mixed Model (GLMM). Cuevas et al. provided a detailed description of this approach, wherein the model equations are derived by stacking the vectors for all individual plants[46]:

      $ \mathrm{y}=\mathrm{X}\mathrm{m}+{\mathrm{Z}}_{1}\mathrm{f}+{\mathrm{Z}}_{2}\mathrm{b}+{\mathrm{Z}}_{3}\mathrm{f}\mathrm{b}+\mathrm{e} $ (4)

      Where the vector y represents the observed values of the overall phenotype. The vector m denotes the population mean values, while b, f, and e represent the vectors of block effects, family effects, family and block interaction effects, and random error effects, respectively. The design matrices X, Z1, Z2 and Z3 are corresponding correlation matrices used to link the observed values with their respective effects. We define the expected value vector (E) and the discrete matrix (Var) as:

      $ E\left[\mathrm{y}\right]=\mathrm{ }\mathrm{X}\mathrm{m} $ (5)
      $ Var\left[\mathrm{f}\right]={\mathrm{Z}}_{1}\otimes {\mathrm{F}}_{0} $ (6)
      $ Var\left[\mathrm{b}\right]={\mathrm{Z}}_{2}\otimes {\mathrm{B}}_{0} $ (7)
      $ {Var}\left[\mathrm{f}\mathrm{b}\right]={{\mathrm{Z}}_{1}\mathrm{Z}}_{2}\otimes {{\mathrm{F}}_{0}\mathrm{B}}_{0} $ (8)
      $ Var\left[\mathrm{e}\right]=\mathrm{Z}\oplus{\mathrm{R}}_{0} $ (9)

      and

      $ {\mathrm{F}}_{0}=\left[\begin{array}{cc}{\sigma }_{{f}_{1}}^{2}& {\sigma }_{{f}_{1}{f}_{2}}\\ {\sigma }_{{f}_{2}{f}_{1}}& {\sigma }_{{f}_{2}}^{2}\end{array}\right] $ (10)
      $ {\mathrm{B}}_{0}=\left[\begin{array}{cc}{\sigma }_{{b}_{1}}^{2}& {\sigma }_{{b}_{1}{b}_{2}}\\ {\sigma }_{{b}_{2}{b}_{1}}& {\sigma }_{{b}_{2}}^{2}\end{array}\right] $ (11)
      $ {{\mathrm{F}}_{0}\mathrm{B}}_{0}=\left[\begin{array}{cc}{\sigma }_{{f}_{1}{b}_{2}}^{2}& {\sigma }_{{f}_{1}{b}_{2}}\\ {\sigma }_{{f}_{2}{b}_{1}}& {\sigma }_{{f}_{2}{b}_{2}}^{2}\end{array}\right] $ (12)
      $ {\mathrm{R}}_{0}=\left[\begin{array}{cc}{\sigma }_{{e}_{1}}^{2}& {\sigma }_{{e}_{1}{e}_{2}}\\ {\sigma }_{{e}_{2}{e}_{1}}& {\sigma }_{{e}_{2}}^{2}\end{array}\right] $ (13)

      Where the $ \otimes $ and $ \oplus $ are vector product and vector addition respectively; $ {\sigma }_{{f}_{i}}^{2} $, $ {\sigma }_{{b}_{i}}^{2} $, $ {\sigma }_{{f}_{i}{b}_{i}}^{2} $ and $ {\sigma }_{{e}_{i}}^{2} $ represent the variances of the block effect, family effect, family × block interaction effect, and environmental error effect for trait i, respectively; $ {\sigma }_{{b}_{i}{b}_{j}} $, $ {\sigma }_{{f}_{i}{f}_{j}} $ and $ {\sigma }_{{e}_{i}{e}_{j}} $ are the covariances effects of block, family and environmental errors for traits i and j. The variance component of the model was used to calculate the family mean heritability (H2) and individual narrow-sense heritability (h2):

      $ {H}_{i}^{2}=\dfrac{{\sigma }_{{f}_{i}}^{2}+{\sigma }_{{b}_{i}}^{2}}{{\sigma }_{{f}_{i}}^{2}+{\sigma }_{{b}_{i}}^{2}+{\sigma }_{{f}_{i}{b}_{i}}^{2}/{n}_{b}+{\sigma }_{{e}_{i}}^{2}/{n}_{b}{n}_{k}} $ (14)
      $ {h}_{i}^{2}=\dfrac{4({\sigma }_{{f}_{i}}^{2}+{\sigma }_{{b}_{i}}^{2})}{{\sigma }_{{f}_{i}}^{2}+{\sigma }_{{b}_{i}}^{2}+{\sigma }_{{f}_{i}{b}_{i}}^{2}+{\sigma }_{{e}_{i}}^{2}} $ (15)

      Phenotypic correlation $ {r}_{{p}_{ij}} $ and genetic correlation $ {r}_{{g}_{ij}} $ of traits i and j:

      $ {r}_{{p}_{ij}}=\dfrac{{\sigma }_{{f}_{ij}}+{{\sigma }_{{b}_{ij}}}+\sigma _{{e}_{ij}}}{\sqrt{\left({\sigma }_{{f}_{i}}^{2}+{{\sigma }_{{b}_{i}}^{2}+{\sigma }_{{f}_{i}{b}_{i}}^{2}}+\sigma _{{e}_{i}}^{2}\right)\left({\sigma }_{{f}_{j}}^{2}+{\sigma }_{{b}_{j}}^{2}+{\sigma }_{{f}_{j}{b}_{j}}^{2}+{\sigma }_{{e}_{j}}^{2}\right)}} $ (16)
      $ {r}_{{g}_{ij}}=\dfrac{{\sigma }_{{f}_{ij}}{+\sigma }_{{b}_{ij}}}{\sqrt{\left({\sigma }_{{f}_{i}}^{2}+{\sigma }_{{b}_{i}}^{2}\right)\left({\sigma }_{{f}_{j}}^{2}+{\sigma }_{{b}_{j}}^{2}\right)}} $ (17)

      Among them, $ {\sigma }_{{f}_{ij}} $, $ {\sigma }_{{b}_{ij}} $, $ {\sigma }_{{e}_{ij}} $ represent the covariance of family effect, block effect and environmental errors effect for traits i and j; $ {\sigma }_{{f}_{i}}^{2} $ and $ {\sigma }_{{f}_{j}}^{2} $ represent the family variance estimate of traits i and j; $ {\sigma }_{{b}_{i}}^{2} $ and $ {\sigma }_{{b}_{j}}^{2} $ represent the block variance estimate of traits i and j; $ {\sigma }_{{e}_{i}}^{2} $ and $ {\sigma }_{{e}_{j}}^{2} $ represent the environmental errors variance estimate of traits i and j. $ {\sigma }_{{b}_{i}{f}_{i}}^{2} $ and $ {\sigma }_{{b}_{j}{f}_{j}}^{2} $ represent the variances of family × block interaction effect; The $ {n}_{b} $ and $ {n}_{k} $ are block numbers and tree numbers per family, respectively.

      Coefficient of variation (CV/%):

      $ {\rm{CV}}=100\overline{x}/\sigma $ (18)

      Coefficient of phenotypic variation (Ip) and Coefficient of genetic Variation (Ig) :

      $ {I}_{p}=\sqrt{{\sigma }_{p}^{2}}\Big/\overline{x} $ (19)
      $ {I}_{g}=\sqrt{{\sigma }_{g}^{2}}\Big/\overline{x} $ (20)

      Among them, $ \sqrt{{\sigma }_{p}^{2}} $ is the square root of the phenotypic variance component; $ \sqrt{{\sigma }_{g}^{2}} $ is the square root of the genetic variance component; $ \overline{x} $ is the trait mean value.

      The genetic gain ($ \Delta G $) is estimated by equation (17) :

      $ \Delta G=i\cdot H\cdot {\sigma }_{g}/t $ (21)

      Among them, i is the selection intensity; $ H $ is represented by the square root of the heritability; $ {\sigma }_{g} $ represents the square root of the additive genetic variance; t represents the breeding period. In this study, we employed varying selection intensities, specifically 1.40, 1.75, and 2.06, to simulate diverse selective pressures.

      Breeding value (BV) is estimated by the best linear unbiased prediction (BLUP) model.

      PCA was employed in the process of conducting multi-trait combined selection of elite families. The principal components (PCs) for growth traits, wood properties, and resin traits were derived from the breeding values associated with each trait. Subsequently, we computed the loadings of each trait across various PCs, along with eigenvalues, eigenvectors, contribution rates, and cumulative contribution rates for each PC. Elite families were selected based on the comprehensive scores obtained from the PC factors (Ftotal) of each family.

    • All data analyses were based on R software[47]. The analysis of variance was calculated by the function aov, the pedigree relationship matrix was calculated by the R language 'pedigreemm' software package[48], and the genetic parameters and breeding values were estimated using the R language 'sommer' software package[49]. PCA analysis was performed by the function prcomp and the software package 'psych'[50]. All result visualizations were implemented using the R language 'ggplot2' software package[51].

    • Table 1 presents statistical results for growth, wood properties, and resin traits in the progeny test stand of Slash pine half-sib families. Growth traits, including DBH, Ht, Huc, Crown, and ARW, had average values of 19.766 cm, 16.773 m, 9.439 m, 3.469 m2, and 4.254 mm, respectively, with CV ranging from 15.724% to 26.931%. For wood properties, Pd, MOE, and AMP had average values of 19.836 mm, 3.006 GPa, and 0.415%, with CVs ranging from 10.335% to 29.488%. Regarding resin traits, average values for OY, Alpha-pinene, Beta-pinene, Abietic acid, and Levopimaric acid were 1.838 g, 14.790 μg·g−1, 10.661 μg·g−1, 8.183 μg·g−1, and 21.631 μg·g−1, respectively. Beta-pinene exhibited the highest CV at 42.126%, while levopimaric acid had the lowest at 14.868%. Significant differences (p < 0.05) were observed in all traits among families. Moreover, except for DBH and ARW, all other traits were significantly influenced by the block effect. Additionally, Ht, Crown, AMP, Alpha-pinene, and Beta-pinene exhibited effects from family, block, and their interaction.

      Table 1.  Results of descriptive statistics and variance analysis.

      TraitsMeanSDCV/%F-statistic
      FamilyBlockBlock × family
      DBH (cm)19.7664.02620.3701.518*0.6160.894
      Ht (m)16.7732.63715.7242.901**7.738**1.902*
      Huc (m)9.4392.33424.7301.393*5.102**0.875
      Crown (m2)3.4690.93426.9311.529*3.362*1.463*
      ARW (mm)4.2540.81819.2281.477*0.4270.622
      Pd (mm)19.8362.05010.3352.637**4.351**1.007
      MOE (Gpa)3.0060.36812.2441.750*3.842*1.239
      AMP (%)0.4150.12229.4882.731**2.902**1.446*
      OY (g)1.8380.71038.6221.618*1.213*0.972
      Alpha_pinene (μg·g−1)14.7904.30129.0781.184*0.649*0.771*
      Beta_pinene (μg·g−1)10.6614.49142.1261.274*1.241*0.929*
      Abietic_acid (μg·g−1)8.1832.99436.5960.887*1.221*1.094
      Levopimaric_acid (μg·g−1)21.6313.21614.8681.341*2.703*1.212
      ** denotes the significance level of p < 0.01, while * indicates p < 0.05.

      The variation in growth traits primarily stems from the interaction between family and block effects (Fig. 1). Notably, the largest effect influenced by environmental factors is the Crown variance at 66.28%, indicating a predominant influence of the environment on Slash pine growth traits. Conversely, variation in wood properties is primarily attributed to family effects, with a proportion of 78.08% for AMP, suggesting a lesser impact of environmental factors. For resin traits, environmental effects emerge as the primary source of variation, while the OY trait shows nearly equal contributions from genetic and environmental factors at 47.40% and 48.33%, respectively.

      Figure 1. 

      The proportion of variance components of different effects for each trait.

    • Based on the variance components, genetic parameters for various traits were estimated (Table 2). Among the three types of traits, Crown, AMP, and Beta-pinene exhibited relatively greater phenotypic and genetic variability, while Ht, Pd, and Levopimaric acid showed lower variability. Ht, Pd, and Beta-pinene demonstrated relatively high family heritabilities of 0.794, 0.855, and 0.695, respectively, indicating significant genetic control. Notably, all three wood properties exhibited relatively high genetic control. The individual heritability of all traits aligned with the family heritability. Particularly, Ht, Pd, and Beta-pinene achieved the highest genetic gain, reaching 11.4%, 11.1%, and 14.8%, respectively, at a 5% selection rate. These findings underscore the exceptional breeding potential of Ht, Pd, and Beta-pinene among all traits in this study.

      Table 2.  Genetic parameter estimation results for various traits.

      Traits $ {I}_{p} $ $ {I}_{g} $ $ {h}_{i}^{2} $ $ {H}_{i}^{2} $ $ \Delta G $
      r = 0.05
      i = 2.06
      r = 0.10
      i = 1.75
      r = 0.20
      i = 1.40
      DBH 0.205 0.041 0.172 (0.017) 0.593 (0.296) 0.107 0.091 0.073
      Ht 0.155 0.092 0.355 (0.021) 0.794 (0.151) 0.114 0.097 0.077
      Huc 0.240 0.091 0.145 (0.017) 0.543 (0.033) 0.052 0.045 0.036
      Crown 0.268 0.115 0.183 (0.017) 0.610 (0.284) 0.026 0.022 0.017
      ARW 0.194 0.113 0.338 (0.024) 0.781 (0.181) 0.035 0.030 0.024
      Pd 0.099 0.072 0.524 (0.024) 0.885 (0.100) 0.111 0.094 0.075
      MOE 0.121 0.063 0.273 (0.019) 0.870 (0.219) 0.015 0.012 0.010
      AMP 0.297 0.170 0.281 (0.020) 0.874 (0.214) 0.005 0.005 0.004
      OY 0.379 0.150 0.157 (0.023) 0.566 (0.042) 0.017 0.015 0.012
      Alpha_pinene 0.290 0.143 0.244 (0.020) 0.694 (0.234) 0.146 0.124 0.099
      Beta_pinene 0.421 0.202 0.245 (0.020) 0.695 (0.251) 0.148 0.126 0.101
      Abietic_acid 0.366 0.061 0.173 (0.013) 0.285 (0.049) 0.022 0.019 0.015
      Levopimaric_acid 0.147 0.061 0.172 (0.010) 0.592 (0.031) 0.084 0.071 0.057
      r represents the selection rate, and i represents the selection intensity.
    • Further phenotypic and genetic correlation analyses were conducted among the 13 traits (Fig. 2). Positive correlations were observed within growth traits and wood properties, while some negative correlations were noted within resin traits. Notably, Beta-pinene showed the highest phenotypic and genetic correlations with OY (rp = 0.62; rg = 0.67), while Alpha-pinene exhibited negative phenotypic (rp = −0.03) and genetic (rg = −0.02) correlations with Beta-pinene and Abietic acid. Among wood properties and resin traits, the lowest phenotypic correlation was between MOE and abietic acid (rp = −0.28). Among growth and resin traits, the highest phenotypic correlation was between DBH and Beta-pinene (rp = 0.79), with a higher genetic correlation (rg = 0.45). In addition, negative phenotypic (rp = −0.32) and genetic (rg = −0.09) correlations were observed with DBH and MOE.

      Figure 2. 

      The analysis of phenotypic and genetic correlation within and among various traits. The upper triangular matrix displays phenotypic correlation, while the lower triangular matrix represents the genetic correlation results. AP, BP, AA, and LA represent Alpha-pinene, Beta-pinene, Abietic acid, and Levopimaric acid, respectively.

    • The BVs for growth traits are outlined in Supplemental Table S2. Family 8-126 showed the highest BV (0.719) for DBH, whereas family 11-26 had the lowest (BV = −0.821). Elite families for DBH, selected at a 20% rate, included 8-126, 8-131, 4-49, 7-258, 0-53, 2-296, and 0-1339. Similarly, elite families for Ht comprised 0-1027, 0-53, 8-126, 0-1339, 7-258, 10-105, and 11-6. For Huc, selected elite families were 3-1, 2-90, 0-636, 8-126, 7-258, 10-73, and 5-12. For ARW, the top 20% selected families were 2-296, 8-13, 11-6, 8-49, 3-1, 0-636, and 10-73. Supplemental Table S3 presents the BVs for wood properties. For MOE, selected elite families were 2-325, 0-464, 0-53, 0-373, 5-12, 4-9, and 11-26. Selected elite families for AMP included 2-325, 10-105, 8-126, 0-465, 0-1339, 4-49, and 0-1077. Supplemental Table S4 provides BVs for resin traits. The top 20% selected families for OY were 0-53, 7-77, 3-1, 0-1339, 8-131, 2-325, and 3-1. Alpha-pinene had the highest BV (1.314) in family 7-77 and the lowest BV (−1.226) in family 2-101. Selected families for Alpha-pinene were 7-77, 8-131, 7-258, 0-373, 0-636, 0-510, and 0-1027. Families selected for Beta-pinene matched those for Alpha-pinene.

    • The results of PCA for growth traits are summarized in Table 3. Following the principle of achieving a cumulative contribution rate of over 85%, the study selected the top three PCs, which accounted for 91.5% and represented the main features of five growth traits. The eigenvalues for these three PCs were 1.689, 1.088, and 0.944, with contributions of 66.4%, 17.3%, and 7.7%, respectively. The weighted sum of individual PC scores, calculated based on the proportion of their eigenvalues to the total selected eigenvalues was utilized to determine the comprehensive scores for growth traits. The comprehensive evaluation model was defined as: FTotal = 0.454F1 + 0.292F2 + 0.254F3. Utilizing this model, comprehensive scores for each family were computed (Fig. 3), ranging from −4.032 to 3.431. The top seven families, identified as elite for growth traits based on a 20% selection rate, were 3-1, 8-126, 2-90, 0-636, 0-1027, 10-73, and 0-53.

      Table 3.  PCA of growth traits of P. elliottii half-sib families.

      Summary Traits PCs
      PC1 PC2 PC3
      Feature vector DBH −0.457 0.654 0.000
      Ht 0.740 0.322 −0.300
      Huc 0.425 0.695 0.177
      Crown −0.769 0.246 0.249
      ARW 0.400 −0.116 0.872
      Eigenvalue 1.689 1.088 0.944
      Contribution rate (%) 0.664 0.173 0.077
      Cummulative contribution rate (%) 0.664 0.837 0.915

      Figure 3. 

      PC factor scores of growth traits in P. elliottii 33 half-sib families. The color of the squares transitions from blue to red, indicating increasing PC scores.

      For wood properties, PCA results are presented in Table 4. The cumulative contribution rate of the top three principal components reached 100%, representing all features of the three wood properties. Eigenvalues for these PCs were 1.182, 0.984, and 0.833, contributing 39.4%, 32.8%, and 27.8%, respectively. The comprehensive evaluation model for wood properties was defined as: FTotal = 0.394F1 + 0.328F2 + 0.278F3. Comprehensive scores ranged from −4.167 to 2.730 (Fig. 4), with elite families identified as 2-325, 0-1077, 0-373, 0-465, 10-105, 5-12, and 0-510.

      Table 4.  PCA of wood properties of P. elliottii half-sib families.

      Summary Traits PCs
      PC1 PC2 PC3
      Feature vector Pd 0.750 0.000 0.659
      MOE 0.673 0.464 −0.576
      AMP −0.409 0.875 0.260
      Eigenvalue 1.182 0.984 0.833
      Contribution rate (%) 0.394 0.328 0.278
      Cummulative contribution rate (%) 0.394 0.722 1.000

      Figure 4. 

      PC factor scores of wood quality traits in P. elliottii 33 half-sib families. The color of the squares transitions from blue to red, indicating increasing PC scores.

      The PCA results for resin traits are presented in Table 5, where the cumulative contribution rate of the top three PCs reached 91.5%. Eigenvalues for these PCs were 2.290, 1.291, and 0.993, contributing 45.8%, 25.8%, and 19.9%, respectively. The comprehensive evaluation model for resin traits was defined as: FTotal = 0.501F1 + 0.282F2 + 0.217F3. Comprehensive score for each family ranged from −3.243 to 3.816 (Fig. 5), with elite families identified as 0-1077, 8-131, 7-258, 7-77, 0-1027, 0-510, and 3-1, based on a 20% selection rate.

      Table 5.  PCA of resin traits of P. elliottii half-sib families.

      Summary Traits PCs
      PC1 PC2 PC3
      Feature vector OY −0.334 0.717 0.494
      Alpha_pinene 0.934 0.330 0.000
      Beta_pinene 0.933 0.331 0.000
      Abietic_acid −0.660 0.572 0.000
      Levopimaric_acid 0.000 −0.481 0.856
      Eigenvalue 2.290 1.291 0.993
      Contribution rate (%) 0.458 0.258 0.199
      Cummulative contribution (%) 0.458 0.716 0.915

      Figure 5. 

      PC factor scores of resin traits in P. elliottii 33 half-sib families. The color of the squares transitions from blue to red, indicating increasing PC scores.

    • Genetic variations are pivotal, directly influencing genetic gains and guiding breeding strategies in forests[12]. The study revealed significant differences among different family lines in traits. Growth traits were notably more influenced by environmental factors, whereas timber traits were primarily shaped by familial effects and resin traits were impacted by both family and environment in comparable proportions. These findings align with previous research on Slash pine[18] and Japanese black pine[52,53]. Crown, AMP, and Beta-pinene display higher variability in this study when contrasted to Korean pine, whereas DBH and Ht show comparatively lower variations[54].

      Additionally, the heritability of growth traits in Slash pine surpassed previous findings[18], while that of wood properties was slightly higher than those related to radiata pine[55], suggesting a stronger genetic influence on wood properties compared to growth traits. Family heritability estimates for monoterpenes were akin to the study of Lei et al.[56] but slightly lower than the result of Li et al.[57], likely influenced by location-specific factors. Moreover, Ht, Pd, and Beta-pinene exhibit higher genetic gains at varying selection rates and intensities, slightly exceeding those reported in previous studies for growth and wood properties[18,58], but demonstrating lower gains for resin traits[12].

    • In practical breeding, considering multiple traits simultaneously is essential for selecting improved genes, yet past studies indicate a trade-off between growth and wood quality[59], underscoring the need for a comprehensive understanding of trait interrelationships to devise rational breeding strategies. This study reveals either non-significant or negative correlations between growth and wood properties, consistent with findings in other tree species like Norway spruce[60] and poplar[61], posing a significant challenge to genetic improvement in forestry[62]. Additionally, selecting trees at different growth stages may yield varying results[63]. Correlative research on Japanese black pine identified significant genetic correlations between growth, morphological, and OY, such as a correlation of 0.73 between DBH and OY[53]. Similarly, this study shows comparable results, with a correlation of 0.79 between DBH and Beta-pinene content, and 0.77 between Ht and resin yield, possibly due to vigorous growth promoting resin duct formation[64,65]. Furthermore, weak positive or negative correlation between wood properties and resin traits, uncommon in other tree species, were observed. Overall, breeders must meticulously consider trait correlations when selecting for improved genes, and the observed correlations in Slash pine offer valuable insights for effective breeding strategies.

    • For a considerable duration, combined selection of multiple traits has been central in both animal and plant breeding[66,67], with models addressing the impact of LD between traits. PCA, a multivariate selection technique, facilitates the exploration of relationships between explanatory variables and correlated traits[68]. In this study, a combined approach using single-trait and PCA-based selection identified elite families for growth, wood properties, and resin traits. Achievable average genetic gains for the three trait categories were 7.4%, 7.2%, and 8.8%, respectively, consistent with a prior study on Slash pine[18]. Genetic gains for growth traits exceeded the 10% for Scots pine, while those for wood properties were comparable to its WD related traits[69]. However, genetic gains for resin traits were lower than those reported for Slash pine[45]. Quantitative traits are subject to varying degrees of gene-environment interactions, with diverse outcomes influenced by factors such as experimental design, the number of tested families, and different growth periods.

    • This study conducted a comprehensive assessment of genetic variation among 33 half-sib families of Slash pine, aiming to identify elite families suitable for industrial use through a multi-trait combined selection approach. Significant differences were noted among families across 13 traits, with growth traits primarily influenced by block × family interaction, wood properties mainly affected by family effects, and resin traits showing variation attributed to both family effects and block × family interaction. Strong genetic control was evident for several traits, notably Beta-pinene, which exhibited the highest variations, and genetic gains, indicating significant breeding potential. Negative correlations were observed between growth and wood properties, while positive correlations were found between growth and resin traits. The multi-trait combined selection successfully identified elite families for growth, wood properties, and resin traits. This study provides important references for the long-term breeding strategies of Slash pine, offering rich genetic resources for genomic breeding and molecular breeding.

    • The authors confirm contribution to the paper as follows: data curation: Ding X, Zhang Y, Sun J; writing - original draft: Ding X; resources: Luan Q, Jiang J; formal analysis: Ding X; investigation: Tan Z, Huang Q, Diao S, Wu Y; writing - reviewing and editing: Luan Q, Jiang J. All authors reviewed the results and approved the final version of the manuscript.

    • The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

      • We wish to extend our gratitude to Haoran Wu, Haonan Ji, and Xiaoshu Zhang for their invaluable assistance throughout the duration of this study. Additionally, we express our sincere appreciation to the dedicated staff at Changle Forest Farms for their diligent efforts in establishing and managing the experimental site. This research was funded by the Science and Technology Innovation 2030 Program, Ministry of Science and Technology, People's Republic of China (2023ZD0405805) and the National Key Research and Development Program of Ministry of Science and Technology of the People's Republic of China (2022YFD2200023-2).

      • The authors declare that they have no conflict of interest.

      • # Authors contributed equally: Xianyin Ding, Yini Zhang, Jiaming Sun

      • Copyright: © 2024 by the author(s). Published by Maximum Academic Press, Fayetteville, GA. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (5)  Table (5) References (69)
  • About this article
    Cite this article
    Ding X, Zhang Y, Sun J, Tan Z, Huang Q, et al. 2024. Genetic selection for growth, wood quality and resin traits of potential Slash pine for multiple industrial uses. Forestry Research 4: e023 doi: 10.48130/forres-0024-0020
    Ding X, Zhang Y, Sun J, Tan Z, Huang Q, et al. 2024. Genetic selection for growth, wood quality and resin traits of potential Slash pine for multiple industrial uses. Forestry Research 4: e023 doi: 10.48130/forres-0024-0020

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return