Development of 51 K liquid-phased probe array for Loblolly and Slash pines and its application to GWAS of Slash pine breeding population

Loblolly pine (Pinus taeda) and Slash pine (Pinus elliottii) are important economically and ecologically tree species worldwide. To implement genomic based breeding strategy, genotyping large number of samples is needed. However, huge and complex pines genome make it difficult to obtain abundant and cost-effectively markers. In this study, we designed a SNP array containing 51 K liquid-phased probes based on the published SNPs and probes of the two species, novel Slash pine-specific long-read transcript, and SNPs obtained from a double digest Restriction-Site Associated DNA sequencing (ddRADseq). We performed targeted capture sequencing using the developed SNP array on three pine species of Loblolly, Slash, and Caribbean pine (Pinus caribaea). The probe capture had an efficiency of 68.26% on average. Principal Component Analysis (PCA), phylogenetic analysis, and genetic structure analysis of 28,670 SNPs genotyped for the three pine species demonstrated the transferability of the SNP array for the three pines. Furthermore, we genotyped 221 range-wide Slash pine genotypes, analyzed their genetic diversity and performed a genome-wide association study (GWAS). We identified 95 SNPs significantly associated with growth and wood quality traits of Slash pine. We have identified potential roles of the MYB2 and MYB46 transcription factors in regulating cell cycle and differentiation in Slash pine, while also suggesting that AP2/ERF and Glycosyl hydrolase may play significant roles in wood formation processes. This new SNP array can be widely applied in future molecular breeding research of conifers.