Comparison of ddRADseq and EUChip60K SNP genotyping systems for population genetics and genomic selection in Eucalyptus dunnii (Maiden)

Front Genet. 2024 Mar 26:15:1361418. doi: 10.3389/fgene.2024.1361418. eCollection 2024.

Abstract

Eucalyptus dunnii is one of the most important Eucalyptus species for short-fiber pulp production in regions where other species of the genus are affected by poor soil and climatic conditions. In this context, E. dunnii holds promise as a resource to address and adapt to the challenges of climate change. Despite its rapid growth and favorable wood properties for solid wood products, the advancement of its improvement remains in its early stages. In this work, we evaluated the performance of two single nucleotide polymorphism, (SNP), genotyping methods for population genetics analysis and Genomic Selection in E. dunnii. Double digest restriction-site associated DNA sequencing (ddRADseq) was compared with the EUChip60K array in 308 individuals from a provenance-progeny trial. The compared SNP set included 8,011 and 19,008 informative SNPs distributed along the 11 chromosomes, respectively. Although the two datasets differed in the percentage of missing data, genome coverage, minor allele frequency and estimated genetic diversity parameters, they revealed a similar genetic structure, showing two subpopulations with little differentiation between them, and low linkage disequilibrium. GS analyses were performed for eleven traits using Genomic Best Linear Unbiased Prediction (GBLUP) and a conventional pedigree-based model (ABLUP). Regardless of the SNP dataset, the predictive ability (PA) of GBLUP was better than that of ABLUP for six traits (Cellulose content, Total and Ethanolic extractives, Total and Klason lignin content and Syringyl and Guaiacyl lignin monomer ratio). When contrasting the SNP datasets used to estimate PAs, the GBLUP-EUChip60K model gave higher and significant PA values for six traits, meanwhile, the values estimated using ddRADseq gave higher values for three other traits. The PAs correlated positively with narrow sense heritabilities, with the highest correlations shown by the ABLUP and GBLUP-EUChip60K. The two genotyping methods, ddRADseq and EUChip60K, are generally comparable for population genetics and genomic prediction, demonstrating the utility of the former when subjected to rigorous SNP filtering. The results of this study provide a basis for future whole-genome studies using ddRADseq in non-model forest species for which SNP arrays have not yet been developed.

Keywords: ABLUP; GBLUP; SNP array; double digest restriction-site associated DNA sequencing; genomic prediction; genotyping by sequencing.

Grants and funding

The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This research was supported by the BIOTECH II platform (grant number 373-780), which funded EUChip60K genotyping and NIR phenotyping; Programa de Sustentabilidad y Competitividad Forestal (PROMEF), Unidad para el Cambio Rural (UCAR), which funded staff salaries; Programa Nacional Forestal de INTA (PNFOR-1104064), which funded trial maintenance and phenotyping of growth traits; Proyecto de Investigación Científica y Tecnológica del Ministerio de Ciencia, Tecnología e Innovación (PICT 2014-0795) and Proyecto INTA 2019-PE-E6-I114-001, which funded the genotyping by ddRADseq; Programa Nacional de Biotecnología de INTA (PNBIO 1131042), which funded DNA extraction and set up ddRADseq protocol. NA was supported by doctoral fellowships from Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET). This work used computational resources from BIOCAD, IABiMo-UEDD INTA-CONICET, Programa de Sustentabilidad y Competitividad Forestal/BID 2853/Oc-Ar, Consorcio Argentino de Tecnología Genómica, MinCyT PPL 2011 004.