Radiata pine (Pinus radiata D.Don) is one of the world’s most domesticated pines and a key economic species in New Zealand. Thus, the development of genomic resources for radiata pine has been a high priority for both research and commercial breeding. Leveraging off a previously developed exome capture panel, we tested the performance of 438,744 single nucleotide polymorphisms (SNPs) on a screening array (NZPRAD01) and then selected 36,285 SNPs for a final genotyping array (NZPRAD02). These SNPs aligned to 15,372 scaffolds from the Pinus taeda L. v. 1.01e assembly, and 20,039 contigs from the radiata pine transcriptome assembly. The genotyping array was tested on more than 8000 samples, including material from archival progenitors, current breeding trials, nursery material, clonal lines, and material from Australia. Our analyses indicate that the array is performing well, with sample call rates greater than 98% and a sample reproducibility of 99.9%. Genotyping in two linkage mapping families indicated that the SNPs are well distributed across the 12 linkage groups. Using genotypic data from this array, we were also able to differentiate representatives of the five recognized provenances of radiata pine, Año Nuevo, Monterey, Cambria, Cedros and Guadalupe. Furthermore, principal component analysis of genotyped trees revealed clear patterns of population structure, with the primary axis of variation driven by provenance ancestry and the secondary axis reflecting breeding activities. This represents the first commercial use of genomics in a radiata pine breeding program.

The selection of drought-tolerant genotypes is globally recognized as an effective strategy to maintain the growth and survival of commercial tree species exposed to future drought periods. New genomic selection tools that reduce the time of progeny trials are required to substitute traditional tree breeding programs. We investigated the genetic variation of water stress tolerance in New Zealand-grown Pinus radiata D. Don using 622 commercially-used genotypes from 63 families. We used quantitative pedigree-based (Genomic Best Linear Unbiased Prediction or ABLUP) and genomic-based (Genomic Best Linear Unbiased Prediction or GBLUP) approaches to examine the heritability estimates associated with water stress tolerance in P. radiata. Tree seedling growth traits, foliar carbon isotope composition (δ¹³C), and dark-adapted chlorophyll fluorescence (Y) were monitored before, during and after 10 months of water stress. Height growth showed a constant and moderate heritability level, while the heritability estimate for diameter growth and δ¹³C decreased with water stress. In contrast, chlorophyll fluorescence exhibited low heritability after 5 and 10 months of water stress. The GBLUP approach provided less breeding value accuracy than ABLUP, however, the relative selection efficiency of GBLUP was greater compared with ABLUP selection techniques. Although there was no significant relationship directly between δ¹³C and Y, the genetic correlations were significant and stronger for GBLUP. The positive genetic correlations between δ¹³C and tree biomass traits under water stress indicated that intraspecific variation in δ¹³C was likely driven by differences in the genotype’s photosynthetic capacity. The results show that foliar δ¹³C can predict P. radiata genotype tolerance to water stress using ABLUP and GBLUP approaches and that such approaches can provide a faster screening and selection of drought-tolerant genotypes for forestry breeding programs.

Traditional tree improvement is cumbersome and costly. Our main objective was to assess the extent to which genomic data can currently accelerate and improve decision making in this field. We used diameter at breast height (DBH) and wood density (WD) data for 4430 tree genotypes and single-nucleotide polymorphism (SNP) data for 2446 tree genotypes. Pedigree reconstruction was performed using a combination of maximum likelihood parentage assignment and matching based on identity-by-state (IBS) similarity. In addition, we used best linear unbiased prediction (BLUP) methods to predict phenotypes using SNP markers (GBLUP), recorded pedigree information (ABLUP), and single-step “blended” BLUP (HBLUP) combining SNP and pedigree information. We substantially improved the accuracy of pedigree records, resolving the inconsistent parental information of 506 tree genotypes. This led to substantially increased predictive ability (i.e., by up to 87%) in HBLUP analyses compared to a baseline from ABLUP. Genomic prediction was possible across populations and within previously untested families with moderately large training populations (N = 800–1200 tree genotypes) and using as few as 2000–5000 SNP markers. HBLUP was generally more effective than traditional ABLUP approaches, particularly after dealing appropriately with pedigree uncertainties. Our study provides evidence that genome-wide marker data can significantly enhance tree improvement. The operational implementation of genomic selection has started in radiata pine breeding in New Zealand, but further reductions in DNA extraction and genotyping costs may be required to realise the full potential of this approach.