Genomic selection for growth and wood properties in multi-generation hybrid populations of Populus deltoides

Abstract

The approximately 20-year breeding cycle has severely restricted the progress of genetic improvement in poplar. Genomic selection (GS) breeding has been demonstrated as an effective approach to accelerate this process. However, its application in forest tree species remains at an early stage. To advance the genetic improvement of target traits in Populus deltoides, the primary species of poplar plantations in China, we systematically implemented GS breeding using 765 hybrid progenies from 32 multi-generational full-sib families. Firstly, we assembled a high-quality genome of one core parent P. deltoides ‘Danhong’, with a genome size of 419.4 Mb and scaffold N50 of 22.0 Mb, which is also the first telomere-to-telomere (T2T) level genome of P. deltoides. Through comparative genomic analysis, we identified 1395 specific structural variants closely associated with growth and development. Subsequently, through genome-wide association studies (GWAS), we identified 135 quantitative trait nucleotides (QTNs) associated with growth and wood quality traits. By systematically evaluating reference genomes, statistical models, and various marker selection strategies, we developed optimal genomic prediction (GP) models for six traits, with the highest prediction accuracy (PA) reaching 0.730 for DBH. Compared with using all markers, the PA was improved by an average of 136.34%. Furthermore, by integrating GP, GWAS, and RNA-seq results, we identified core breeding parents and elite clones for P. deltoides genetic improvement and discovered important candidate genes. Our results provide a promising strategy for accelerating breeding cycles and genetic improvement, offering valuable breeding and genetic resources for forest tree improvement.

Acknowledgements

This work was supported by the Major Project of Agricultural Biological Breeding (2022ZD0401501), the National Key Research and Development Program of China (2021YFD2200201) and the National Natural Science Foundation (32071797).

Author contributions

X.Z., L.Z. and J.H. designed the project and edited the manuscript. X.Z., M.Z., Y.B., and H.W. collected the experimental materials and collect data. X.Z. and L.Z. performed the genome analyses. X.Z. performed the resequencing and genomic selection analyses. X.Z. and J.T. contributed to the interpretation of results. X.Z. and J.H. wrote and revised the manuscript. All authors contributed to the article and approved the submitted version.

Data availability

DHY genome sequencing data have been deposited under National Genomics Data Center (https://ngdc.cncb.ac.cn/?lang=en) under BioProject PRJCA033198. The whole-genome resequencing data have been deposited under BioProject PRJCA033212. The transcriptome sequencing data have been deposited under BioProject PRJCA033213. All other data related to this study are provided in the supplementary files of the article.

Conflict of interest statement

The authors declare no conflict of interest.

Supplementary data

Supplementary data is available at Horticulture Research online.

References

[1] Chen Y, Wei S, Guanzheng Q. et al. The key and core technologies for accelerating the tree breeding process. J Nanjing For Univ (Nat Sci Ed). 2022; 46:1-9

[2] Ma T, Hou J. Trait assessment of 1122 Populus deltoides clones: unveiling correlations among growth, wood properties, and dis-ease susceptibility. Forests. 2024; 15:1250

[3] Zhang C, Li S, Zhao Z. et al. A new poplar variety Populus deltoides CL. ’Danhong’. Sci Silvae Sin. 2008; 44:1

[4] Dong M, He J. Linking the past to the future: a reality check on cross-border timber trade from Myanmar (Burma) to China. Forest Policy Econ. 2018; 87:11-9

[5] Solovieff N, Cotsapas C, Lee PH. et al. Pleiotropy in complex traits: challenges and strategies. Nat Rev Genet. 2013; 14:483-95

[6] Alemu A, Åstrand J, Montesinos-López OA. et al. Genomic selec-tion in plant breeding: key factors shaping two decades of progress. Mol Plant. 2024; 17:552-78

[7] Wientjes YCJ, Bijma P, Calus MPL. Optimizing genomic reference populations to improve crossbred performance. Genet Sel Evol. 2020; 52:65

[8] XuY, LiuX, FuJ. et al. Enhancing genetic gain through genomic selection: from livestock to plants. Plant Commun. 2020; 1:100005

[9] Xu Y, Zhang X, Li H. et al. Smart breeding driven by big data, arti-ficial intelligence, and integrated genomic-enviromic prediction. Mol Plant. 2022; 15:1664-95

[10] Ahmar S, Ballesta P, Ali M. et al. Achievements and challenges of genomics-assisted breeding in forest trees: from marker-assisted selection to genome editing. Int J Mol Sci. 2021; 22:10583

[11] Crossa J, Pérez-Rodríguez P, Cuevas J. et al. Genomic selection in plant breeding: methods, models, and perspectives. Trends Plant Sci. 2017; 22:961-75

[12] Wang K, Abid MA, Rasheed A. et al. DNNGP, a deep neural network-based method for genomic prediction using multi-omics data in plants. Mol Plant. 2023; 16:279-93

[13] Juliana P, Poland J, Huerta-Espino J. et al. Improving grain yield, stress resilience and quality of bread wheat using large-scale genomics. Nat Genet. 2019; 51:1530-9

[14] Rabier CE, Barre P, Asp T. et al. On the accuracy of genomic selection. PLoS One. 2016; 11:e0156086

[15] Wang J, Zhang Z. GAPIT version 3: boosting power and accu-racy for genomic association and prediction. Genomics Proteomics Bioinformatics. 2021; 19:629-40

[16] Akutsu H, Na’iem M, Widiyatno et al. Comparing modeling meth-ods of genomic prediction for growth traits of a tropical timber species, Shorea macrophylla. Front Plant Sci 2023; 14: 1241908.

[17] Bernhardsson C, Zan Y, Chen Z. et al. Development of a highly efficient 50K single nucleotide polymorphism genotyping array for the large and complex genome of Norway spruce (Picea abies L. Karst) by whole genome resequencing and its transferability to other spruce species. Mol Ecol Resour. 2021; 21:880-96

[18] Mphahlele MM, Isik F, Mostert-O’Neill MM. et al. Expected ben-efits of genomic selection for growth and wood quality traits in Eucalyptus grandis. Tree Genet Genomes. 2020; 16:49

[19] Gamal El-Dien O, Shalev TJ, Yuen MMS. et al. Genomic selection in western redcedar: from proof of concept to operational appli-cation. New Phytol. 2024; 244:588-602

[20] Sun M, Zhang M, Kumar S. et al. Genomic selection of eight fruit traits in pear. Hortic Plant J. 2024; 10:318-26

[21] Xia H, Hao Z, Shen Y. et al. Genome-wide association study of multiyear dynamic growth traits in hybrid Liriodendron identifies robust genetic loci associated with growth trajectories. Plant J. 2023; 115:1544-63

[22] Alves FC, Balmant KM, Resende MFR Jr. et al. Accelerating forest tree breeding by integrating genomic selection and greenhouse phenotyping. Plant Genome. 2020; 13:e20048

[23] Tuskan GA, DiFazio S, Jansson S. et al. The genome of black cot-tonwood, populus trichocarpa (Torr. & Gray). Science. 2006; 313: 1596-604

[24] Bai S, Wu H, Zhang J. et al. Genome assembly of salicaceae populus deltoides (eastern cottonwood) I-69 based on nanopore sequencing and Hi-C technologies. J Hered. 2021; 112:303-10

[25] Shi T-L, Jia K-H, Bao Y-T. et al. High-quality genome assembly enables prediction of allele-specific gene expression in hybrid poplar. Plant Physiol. 2024; 195:652-70

[26] Liu W, Liu C, Chen S. et al. A nearly gapless, highly contigu-ous reference genome for a doubled haploid line of Populus ussuriensis, enabling advanced genomic studies. For Res. 2024; 4:e019

[27] Li H, Durbin R. Genome assembly in the telomere-to-telomere era. Nat Rev Genet. 2024; 25:658-70

[28] Li F, Xu S, Xiao Z. et al. Gap-free genome assembly and comparative analysis reveal the evolution and anthocyanin accumulation mechanism of Rhodomyrtus tomentosa. Hortic Res. 2023;10:uhad005

[29] Han Y, Wu J, Zhu Q-H. et al. Assembling a reference-quality genome and resequencing diverse accessions of Beckmannia syzigachne provide insights into population structure and gene family evolution. Plant Commun. 2024; 13:101174

[30] Wang Y, Zhao L, Wang D. et al. Four near-complete genome assemblies reveal the landscape and evolution of centromeres in Salicaceae. Genome Biol. 2025; 26:111

[31] Shi T, Zhang X, Hou Y. et al. The super-pangenome of Populus unveils genomic facets for its adaptation and diversification in widespread forest trees. Mol Plant. 2024; 17:725-46

[32] Letunic I, Bork P. Interactive Tree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49:W293-6

[33] Chen Z, Baison J, Pan J. et al. Accuracy of genomic selection for growth and wood quality traits in two control-pollinated progeny trials using exome capture as the genotyping platform in Norway spruce. BMC Genomics. 2018; 19:946

[34] Hickey JM, Chiurugwi T, Mackay I. et al. Genomic prediction unifies animal and plant breeding programs to form platforms for biological discovery. Nat Genet. 2017; 49:1297-303

[35] Resende MDV, Resende MFR Jr, Sansaloni CP. et al. Genomic selection for growth and wood quality in eucalyptus: capturing the missing heritability and accelerating breeding for complex traits in forest trees. New Phytol. 2012; 194:116-28

[36] Grattapaglia D. Twelve years into genomic selection in forest trees: climbing the slope of enlightenment of marker assisted tree breeding. Forests. 2022; 13:1554

[37] Liu X, Wang H, Wang H. et al. Factors affecting genomic selection revealed by empirical evidence in maize. Crop J. 2018; 6:341-52

[38] Grattapaglia D. tatus and perspectives of genomic selection in forest tree breeding. In: Varshney RK, Roorkiwal M, Sorrells ME, Springer International Publishing,eds. Genomic Selection for Crop Improvement: New Molecular Breeding Strategies for Crop Improvement. 2017,199-249

[39] Sallam AH, Endelman JB, Jannink J-L. et al. Assessing genomic selection prediction accuracy in a dynamic barley breeding population. Plant Genome. 2015; 8:1-15

[40] Clark SA, Hickey JM, Daetwyler HD. et al. The importance of information on relatives for the prediction of genomic breeding values and the implications for the makeup of reference data sets in livestock breeding schemes. Genet Sel Evol. 2012; 44:4

[41] Weber SE, Frisch M, Snowdon RJ. et al. Haplotype blocks for genomic prediction: a comparative evaluation in multiple crop datasets. Front Plant Sci. 2023; 14:1217589

[42] Xu H, Wang Z, Wang F. et al. Genome-wide association study and genomic selection of spike-related traits in bread wheat. Theor Appl Genet. 2024; 137:131

[43] Minamikawa MF, Nonaka K, Kaminuma E. et al. Genome-wide association study and genomic prediction in citrus: potential of genomics-assisted breeding for fruit quality traits. Sci Rep. 2017; 7:4721

[44] Ali M, Zhang Y, Rasheed A. et al. Genomic prediction for grain yield and yield-related traits in Chinese winter wheat. Int J Mol Sci. 2020; 21:1342

[45] Wang Z, Hu H, Sun T. et al. Genomic selection improves inner shell purpleness in triangle sail mussel Hyriopsis cumingii (Lea, 1852). Aquaculture. 2023; 575:739815

[46] Estopa RA, Paludeto JGZ, Müller BSF. et al. Genomic prediction of growth and wood quality traits in Eucalyptus benthamii using different genomic models and variable SNP genotyping density. New For. 2023; 54:343-62

[47] Li P, Lu W, Xiao L. et al. Progress and prospect of genome wide association study (GWAS) in forest trees. Sci Sin Vitae. 2020; 50: 144-53

[48] Zhou X, Xiang X, Cao D. et al. Selective sweep and GWAS provide insights into adaptive variation of Populus cathayana leaves. For Res. 2024; 4:e012

[49] Osorio-Guarín JA, Garzón-Martínez GA, Delgadillo-Duran P. et al. Genome-wide association study (GWAS) for morpho-logical and yield-related traits in an oil palm hybrid (Elaeis oleifera x Elaeis guineensis) population. BMC Plant Biol. 2019; 19: 533

[50] O’Connor K, Hayes B, Hardner C. et al. Genome-wide association studies for yield component traits in a macadamia breeding population. BMC Genomics. 2020; 21:199

[51] Zhang F, Liu X, Xia H. et al. Identification of genetic loci for growth and stem form traits in hybrid Liriodendron via a genome-wide association study. For Res. 2025; 5:e001

[52] Li P, Xiao L, Du Q. et al. Genomic insights into selection for heterozygous alleles and woody traits in Populus tomentosa. Plant Biotechnol J. 2023; 21:2002-18

[53] Cheng H, Concepcion GT, Feng X. et al. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods. 2021; 18:170-5

[54] Heng L. Minimap2: pairwise alignment for nucleotide sequences. Bioinformatics. 2018; 18:3094-3100

[55] Zhang X, Zhang S, Zhao Q. et al. Assembly of allele-aware, chromosomal-scale autopolyploid genomes based on Hi-C data. Nat Plants. 2019; 5:833-45

[56] Durand NC, Robinson JT, Shamim MS. et al. Juicebox provides a visualization system for Hi-C contact maps with unlimited zoom. Cell Syst. 2016; 3:99-101

[57] Li H, Durbin R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics. 2010; 26:589-95

[58] Rhie A, Walenz BP, Koren S. et al. Merqury: reference-free quality, completeness, and phasing assessment for genome assemblies. Genome Biol. 2020; 21:245

[59] Lin Y, Ye C, Li X. et al. quarTeT: a telomere-to-telomere toolkit for gap-free genome assembly and centromeric repeat identifi-cation. Hortic Res. 2023;10:uhad127

[60] Gary B. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999; 27:573-80

[61] Danecek P, McCarthy SA. BCFtools/csq: haplotype-aware variant consequences. Bioinformatics. 2017; 33:2037-9

[62] Yu G, Wang LG, Han Y. et al. clusterProfiler: an R package for comparing biological themes among gene clusters. Omics. 2012; 16:284-7

[63] Raphael B, Zhi D, Tang H. et al. A novel method for multiple alignment of sequences with repeated and shuffled elements. Genome Res. 2004; 14:2336-46

[64] Tarailo-Graovac M, Chen N.Using RepeatMasker to identify repetitive elements in genomic sequences. Curr Protoc Bioinfor-matics. 2009; Chapter 4:4.10.11-14.10.14

[65] Keller O, Kollmar M, Stanke M. et al. A novel hybrid gene predic-tion method employing protein multiple sequence alignments. Bioinformatics. 2011; 27:757-63

[66] Leskovec J, Sosicˇ R. SNAP: a general purpose network analysis and graph mining library. ACM Trans Intell Syst Technol. 2016; 8: 1-20

[67] Altschul SF, Gish W, Miller W. et al. Basic local alignment search tool. J Mol Biol. 1990; 215:403-10

[68] Birney E, Clamp M, Durbin R. GeneWise and genomewise. Genome Res. 2004; 14:988-95

[69] John B, Alexandre L, Mark B. GeneMarkS: a self-training method for prediction of gene starts in microbial genomes. Implications for finding sequence motifs in regulatory regions. Nucleic Acids Res. 2001; 29:2607-18

[70] Haas BJ, Salzberg SL, Zhu W. et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 2008;9:R7

[71] Haas BJ, Delcher AL, Mount SM. et al. Improving the Arabidopsis genome annotation using maximal transcript alignment assem-blies. Nucleic Acids Res. 2003; 31:5654-66

[72] Nawrocki EP, Eddy SR. Infernal 1.1: 100-fold faster RNA homol-ogy searches. Bioinformatics. 2013; 29:2933-5

[73] Chan PP, Lowe TM. tRNAscan-SE: searching for tRNA genes in genomic sequences. Methods Mol Biol. 2019; 1962:1-14

[74] Wang Y, Tang H, Debarry JD. et al. MCScanX: a toolkit for detec-tion and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 2012; 40:e49

[75] Riva G, Mauri M. Mummer: how robotics can reboot social inter-action and customer engagement in shops and malls. Cyberpsy-chol Behav Soc Netw. 2021; 24:210-1

[76] Goel M, Sun H, Jiao WB. et al. SyRI: finding genomic rear-rangements and local sequence differences from whole-genome assemblies. Genome Biol. 2019; 20:277

[77] Goel M, Schneeberger K. Plotsr: visualizing structural similar-ities and rearrangements between multiple genomes. Bioinfor-matics. 2022; 38:2922-6

[78] Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010; 26:841-2

[79] McKenna A, Hanna M, Banks E. et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010; 20:1297-303

[80] Danecek P, Auton A, Abecasis G. et al. The variant call format and VCFtools. Bioinformatics. 2011; 27:2156-8

[81] Browning BL, Zhou Y, Browning SR. A one-penny imputed genome from next-generation reference panels. Am J Hum Genet. 2018; 103:338-48

[82] Alexander DH, Novembre J, Lange K. Fast model-based estima-tion of ancestry in unrelated individuals. Genome Res. 2009; 19: 1655-64

[83] Retief JD. Phylogenetic analysis using PHYLIP. Methods Mol Biol. 2000; 132:243-58

[84] Zhang C, Dong SS, Xu JY. et al. PopLDdecay: a fast and effective tool for linkage disequilibrium decay analysis based on variant call format files. Bioinformatics. 2019; 35:1786-8

[85] Li M, Zhang YW, Xiang Y. et al. IIIVmrMLM: the R and C++ tools associated with 3VmrMLM, a comprehensive GWAS method for dissecting quantitative traits. Mol Plant. 2022; 15: 1251-3

[86] Purcell S, Neale B, Todd-Brown K. et al. PLINK: a tool set for whole-genome association and population-based linkage anal-yses. Am J Hum Genet. 2007; 81:559-75

[87] Garrison E, Kronenberg ZN, Dawson ET. et al. A spectrum of free software tools for processing the VCF variant call format: vcflib, bio-VCF, cyvcf2, hts-nim and slivar. PLoS Comput Biol. 2022; 18:e1009123

[88] Liu S, Wang Z, Zhu R. et al. Three differential expression analysis methods for RNA sequencing: limma, EdgeR, DESeq2. JVisExp. 2021; 18:e62528