High density linkage maps, genetic architecture, and genomic prediction of growth and wood properties in Pinus radiata.

BACKGROUND: The growing availability of genomic resources in radiata pine paves the way for significant advances in fundamental and applied genomic research. We constructed robust high-density linkage maps based on exome-capture genotyping in two F1 populations, and used these populations to perform quantitative trait locus (QTL) scans, genomic prediction and quantitative analyses of genetic architecture for key traits targeted by tree improvement programmes. RESULTS: Our mapping approach used probabilistic error correction of the marker data, followed by an iterative approach based on stringent parameters. This approach proved highly effective in producing high-density maps with robust marker orders and realistic map lengths (1285-4674 markers per map, with sizes ranging from c. 1643-2292 cM, and mean marker intervals of 0.7-2.1 cM). Colinearity was high between parental linkage maps, although there was evidence for a large chromosomal rearrangement (affecting ~ 90 cM) in one of the parental maps. In total, 28 QTL were detected for growth (stem diameter) and wood properties (wood density and fibre properties measured by Silviscan) in the QTL discovery population, with 1-3 QTL of small to moderate effect size detected per trait in each parental map. Four of these QTL were validated in a second, unrelated F1 population. Results from genomic prediction and analyses of genetic architecture were consistent with those from QTL scans, with wood properties generally having moderate to high genomic heritabilities and predictive abilities, as well as somewhat less complex genetic architectures, compared to growth traits. CONCLUSIONS: Despite the economic importance of radiata pine as a plantation forest tree, robust high-density linkage maps constructed from reproducible, sequence-anchored markers have not been published to date. The maps produced in this study will be a valuable resource for several applications, including the selection of marker panels for genomic prediction and anchoring a recently completed de novo whole genome assembly. We also provide the first map-based evidence for a large genomic rearrangement in radiata pine. Finally, results from our QTL scans, genomic prediction, and genetic architecture analyses are informative about the genomic basis of variation in important phenotypic traits.

Effect of trait's expression level on single-step genomic evaluation of resistance to Dothistroma needle blight.

BACKGROUND: Many conifer breeding programs are paying increasing attention to breeding for resistance to needle disease due to the increasing importance of climate change. Phenotyping of traits related to resistance has many biological and temporal constraints that can often confound the ability to achieve reliable phenotypes and consequently, reliable genetic progress. The development of next generation sequencing platforms has also enabled implementation of genomic approaches in species lacking robust reference genomes. Genomic selection is, therefore, a promising strategy to overcome the constraints of needle disease phenotyping. RESULTS: We found high accuracy in the prediction of genomic breeding values in the disease-related traits that were well characterized, reaching 0.975 for genotyped individuals and 0.587 for non-genotyped individuals. This compared well with pedigree-based accuracies of up to 0.746. Surprisingly, poorly phenotyped disease traits also showed very high accuracy in terms of correlation of predicted genomic breeding values with pedigree-based counterparts. However, this was likely caused by the fact that both were clustered around the population mean, while deviations from the population mean caused by genetic effects did not appear to be well described. Caution should therefore be taken with the interpretation of results in poorly phenotyped traits. CONCLUSIONS: Implementation of genomic selection in this test population of Pinus radiata resulted in a relatively high prediction accuracy of needle loss due to Dothistroma septosporum compared with a pedigree-based approach. Using genomics to avoid biological/temporal constraints where phenotyping is reliable appears promising. Unsurprisingly, reliable phenotyping, resulting in good heritability estimates, is a fundamental requirement for the development of a reliable prediction model. Furthermore, our results are also specific to the single pathogen mating-type that is present in New Zealand, and may change with future incursion of other pathogen varieties. There is no doubt, however, that once a robust genomic prediction model is built, it will be invaluable to not only select for host tolerance, but for other economically important traits simultaneously. This tool will thus future-proof our forests by mitigating the risk of disease outbreaks induced by future changes in climate.

Genotype-by-environment interaction in coast redwood outside natural distribution - search for environmental cues.

BACKGROUND: Effective matching of genotypes and environments is required for the species to reach optimal productivity and act effectively for carbon sequestration. A common garden experiment across five different environments was undertaken to assess genotype x environment interaction (GxE) of coast redwood in order to understand the performance of genotypes across environments. RESULTS: The quantitative genetic analysis discovered no GxE between investigated environments for diameter at breast height (DBH). However, no genetic component was detected at one environment possibly due to stressful conditions. The implementation of universal response function allowed for the identification of important environmental factors affecting species productivity. Additionally, this approach enabled us to predict the performance of species across the New Zealand environmental conditions. CONCLUSIONS: In combination with quantitative genetic analysis which identified genetically superior material, the URF model can directly identify the optimal geographical regions to maximize productivity. However, the finding of ideally uncorrelated climatic variables for species with narrow ecological amplitude is rather challenging, which complicates construction of informative URF model. This, along with a small number of tested environments, tended to overfit a prediction model which resulted in extreme predictions in untested environments.

Genomic selection for non-key traits in radiata pine when the documented pedigree is corrected using DNA marker information.

BACKGROUND: Non-key traits (NKTs) in radiata pine (Pinus radiata D. Don) refer to traits other than growth, wood density and stiffness, but still of interest to breeders. Branch-cluster frequency, stem straightness, external resin bleeding and internal checking are examples of such traits and are targeted for improvement in radiata pine research programmes. Genomic selection can be conducted before the performance of selection candidates is available so that generation intervals can be reduced. Radiata pine is a species with a long generation interval, which if reduced could significantly increase genetic gain per unit of time. The aim of this study was to evaluate the accuracy and predictive ability of genomic selection and its efficiency over traditional forward selection in radiata pine for the following NKTs: branch-cluster frequency, stem straightness, internal checking, and external resin bleeding. RESULTS: Nine hundred and eighty-eight individuals were genotyped using exome capture genotyping by sequencing (GBS) and 67,168 single nucleotide polymorphisms (SNPs) used to develop genomic estimated breeding values (GEBVs) with genomic best linear unbiased prediction (GBLUP). The documented pedigree was corrected using a subset of 704 SNPs. The percentage of trio parentage confirmed was about 49% and about 50% of parents were re-assigned. The accuracy of GEBVs was 0.55-0.75 when using the documented pedigree and 0.61-0.80 when using the SNP-corrected pedigree. A higher percentage of additive genetic variance was explained and a higher predictive ability was observed when using the SNP-corrected pedigree than using the documented pedigree. With the documented pedigree, genomic selection was similar to traditional forward selection when assuming a generation interval of 17 years, but worse than traditional forward selection when assuming a generation interval of 14 years. After the pedigree was corrected, genomic selection led to 37-115% and 13-77% additional genetic gain over traditional forward selection when generation intervals of 17 years and 14 years were assumed, respectively. CONCLUSION: It was concluded that genomic selection with a pedigree corrected by SNP information was an efficient way of improving non-key traits in radiata pine breeding.

Modelling of population structure through contemporary groups in genetic evaluation.

BACKGROUND: Forest trees can occupy extensive geography and environmentally highly variable areas which result in high genetic variability in the direction of pressure from natural selection. At the same time, the majority of conifer species are wind-pollinated from both short and long distances, resulting in wide-spread gene flow, which can lead to maladaptation to local conditions. Quantitative analyses of provenance/progeny tests correct for genetic differences between populations to ensure unbiased genetic parameters are obtained. Commonly, the provenance effect is fitted as a fixed term or can be implemented as a contemporary group in the pedigree. RESULTS: The use of a provenance effect, either as a fixed term or as the same contemporary groups in both maternal and paternal sides of the pedigree, resulted in fairly similar precision of genetic parameters in our case. However, when we developed a phantom contemporary group for the paternal side of the pedigree that considered a different genetic quality of pollen compared with the maternal contribution from trees in the local environment, the model fit and accuracy of breeding values increased. CONCLUSION: Consideration of the mating dynamics and the vector of gene flow are important factors in modelling contemporary genetic groups, particularly when implementing pedigrees within a mixed model framework to obtain unbiased estimates of genetic parameters. This approach is especially important in traits involved in local adaptation.

Correction to: Efficiency of genomic prediction across two Eucalyptus nitens seed orchards with different selection histories.

Published online: 6 July 2018This article was originally published under standard licence, but has now been made available under a [CC BY 4.0] license. The PDF and HTML versions of the paper have been modified accordingly.

Efficiency of genomic prediction across two Eucalyptus nitens seed orchards with different selection histories.

Genomic selection is expected to enhance the genetic improvement of forest tree species by providing more accurate estimates of breeding values through marker-based relationship matrices compared with pedigree-based methodologies. When adequately robust genomic prediction models are available, an additional increase in genetic gains can be made possible with the shortening of the breeding cycle through elimination of the progeny testing phase and early selection of parental candidates. The potential of genomic selection was investigated in an advanced Eucalyptus nitens breeding population focused on improvement for solid wood production. A high-density SNP chip (EUChip60K) was used to genotype 691 individuals in the breeding population, which represented two seed orchards with different selection histories. Phenotypic records for growth and form traits at age six, and for wood quality traits at age seven were available to build genomic prediction models using GBLUP, which were compared to the traditional pedigree-based alternative using BLUP. GBLUP demonstrated that breeding value accuracy would be improved and substantial increases in genetic gains towards solid wood production would be achieved. Cross-validation within and across two different seed orchards indicated that genomic predictions would likely benefit in terms of higher predictive accuracy from increasing the size of the training data sets through higher relatedness and better utilization of LD.

Effect of Hidden Relatedness on Single-Step Genetic Evaluation in an Advanced Open-Pollinated Breeding Program.

Open-pollinated (OP) mating is frequently used in forest tree breeding due to the relative temporal and financial efficiency of the approach. The trade-off is the lower precision of the estimated genetic parameters. Pedigree/sib-ship reconstruction has been proven as a tool to correct and complete pedigree information and to improve the precision of genetic parameter estimates. Our study analyzed an advanced generation Eucalyptus population from an OP breeding program using single-step genetic evaluation. The relationship matrix inferred from sib-ship reconstruction was used to rescale the marker-based relationship matrix (G matrix). This was compared with a second scenario that used rescaling based on the documented pedigree. The proposed single-step model performed better with respect to both model fit and the theoretical accuracy of breeding values. We found that the prediction accuracy was superior when using the pedigree information only when compared with using a combination of the pedigree and genomic information. This pattern appeared to be mainly a result of accumulated unrecognized relatedness over several breeding cycles, resulting in breeding values being shrunk toward the population mean. Using biased, pedigree-based breeding values as the base with which to correlate predicted GEBVs, resulted in the underestimation of prediction accuracies. Using breeding values estimated on the basis of sib-ship reconstruction resulted in increased prediction accuracies of the genotyped individuals. Therefore, selection of the correct base for estimation of prediction accuracy is critical. The beneficial impact of sib-ship reconstruction using G matrix rescaling was profound, especially in traits with inbreeding depression, such as stem diameter.

Chasing genetic correlation breakers to stimulate population resilience to climate change.

Global climate change introduces new combinations of environmental conditions, which is expected to increase stress on plants. This could affect many traits in multiple ways that are as yet unknown but will likely require the modification of existing genetic relationships among functional traits potentially involved in local adaptation. Theoretical evolutionary studies have determined that it is an advantage to have an excess of recombination events under heterogeneous environmental conditions. Our study, conducted on a population of radiata pine (Pinus radiata D. Don), was able to identify individuals that show high genetic recombination at genomic regions, which potentially include pleiotropic or collocating QTLs responsible for the studied traits, reaching a prediction accuracy of 0.80 in random cross-validation and 0.72 when whole family was removed from the training population and predicted. To identify these highly recombined individuals, a training population was constructed from correlation breakers, created through tandem selection of parents in the previous generation and their consequent mating. Although the correlation breakers showed lower observed heterogeneity possibly due to direct selection in both studied traits, the genomic regions with statistically significant differences in the linkage disequilibrium pattern showed higher level of heretozygosity, which has the effect of decomposing unfavourable genetic correlation. We propose undertaking selection of correlation breakers under current environmental conditions and using genomic predictions to increase the frequency of these 'recombined' individuals in future plantations, ensuring the resilience of planted forests to changing climates. The increased frequency of such individuals will decrease the strength of the population-level genetic correlations among traits, increasing the opportunity for new trait combinations to be developed in the future.

Genetic Variation in Drought-Tolerance Traits and Their Relationships to Growth in Pinus radiata D. Don Under Water Stress.

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 (δ13C), 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 δ13C 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 δ13C and Y, the genetic correlations were significant and stronger for GBLUP. The positive genetic correlations between δ13C and tree biomass traits under water stress indicated that intraspecific variation in δ13C was likely driven by differences in the genotype's photosynthetic capacity. The results show that foliar δ13C 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.

Indication of Quantitative Multiple Disease Resistance to Foliar Pathogens in Pinus radiata D.Don in New Zealand.

Increasing resistance against foliar diseases is an important goal in the Pinus radiata D.Don breeding program in New Zealand, and screening for resistance has been in place for some time, since the late 1960s. The current study presents results of four progeny trials within the breeding program to investigate whether multiple disease resistance could be detected against three different needle diseases in P. radiata: Dothistroma needle blight (DNB) caused by Dothistroma septosporum, Cyclaneusma needle cast (CNC) caused by Cyclaneusma minus, and red needle cast (RNC) caused by Phytophthora pluvialis. Four progeny trials in the North Island of New Zealand were available to estimate heritabilities and between-trait genetic correlations. Two of the trials were assessed for DNB, involving 63 full-sib families. A third trial was assessed for CNC, involving 172 half-sib families, and a fourth trial was assessed for RNC, involving 170 half-sib families. Disease resistances had moderate estimates of heritability (0.28-0.48) in all trials. We investigated the potential for multiple disease resistance to the three foliar diseases by estimating genetic correlations between disease resistances using a spatial linear mixed model. The correlation between DNB and CNC resistance was favorable and strong (0.81), indicating that genotypes that are highly resistant to DNB also have a high resistance to CNC. These results suggest that selection based on resistance to DNB could allow for simultaneous indirect selection for resistance to CNC, usually only expressed at a later age. This would allow selections to be made earlier due to the earlier expression of DNB than CNC and reduce the number of expensive disease assessments being undertaken. Conversely, genetic correlation estimates for RNC with DNB and CNC were close to zero, and very imprecise. As such, later-age assessments for this disease would still be required.

Marker Selection in Multivariate Genomic Prediction Improves Accuracy of Low Heritability Traits.

Multivariate analysis using mixed models allows for the exploration of genetic correlations between traits. Additionally, the transition to a genomic based approach is simplified by substituting classic pedigrees with a marker-based relationship matrix. It also enables the investigation of correlated responses to selection, trait integration and modularity in different kinds of populations. This study investigated a strategy for the construction of a marker-based relationship matrix that prioritized markers using Partial Least Squares. The efficiency of this strategy was found to depend on the correlation structure between investigated traits. In terms of accuracy, we found no benefit of this strategy compared with the all-marker-based multivariate model for the primary trait of diameter at breast height (DBH) in a radiata pine (Pinus radiata) population, possibly due to the presence of strong and well-estimated correlation with other highly heritable traits. Conversely, we did see benefit in a shining gum (Eucalyptus nitens) population, where the primary trait had low or only moderate genetic correlation with other low/moderately heritable traits. Marker selection in multivariate analysis can therefore be an efficient strategy to improve prediction accuracy for low heritability traits due to improved precision in poorly estimated low/moderate genetic correlations. Additionally, our study identified the genetic diversity as a factor contributing to the efficiency of marker selection in multivariate approaches due to higher precision of genetic correlation estimates.

A high-density exome capture genotype-by-sequencing panel for forestry breeding in Pinus radiata.

Development of genome-wide resources for application in genomic selection or genome-wide association studies, in the absence of full reference genomes, present a challenge to the forestry industry, where longer breeding cycles could benefit from the accelerated selection possible through marker-based breeding value predictions. In particular, large conifer megagenomes require a strategy to reduce complexity, whilst ensuring genome-wide coverage is achieved. Using a transcriptome-based reference template, we have successfully developed a high density exome capture genotype-by-sequencing panel for radiata pine (Pinus radiata D.Don), capable of capturing in excess of 80,000 single nucleotide polymorphism (SNP) markers with a minor allele frequency above 0.03 in the population tested. This represents approximately 29,000 gene models from a core set of 48,914 probes. A set of 704 SNP markers capable of pedigree reconstruction and differentiating individual genotypes were tested within two full-sib mapping populations. While as few as 70 markers could reconstruct parentage in almost all cases, the impact of missing genotypes was noticeable in several offspring. Therefore, 60 sets of 110 randomly selected SNP markers were compared for both parentage reconstruction and clone differentiation. The performance in parentage reconstruction showed little variation over 60 iterations. However, there was notable variation in discriminatory power between closely related individuals, indicating a higher density SNP marker panel may be required to elucidate hidden relationships in complex pedigrees.

Approaches to variant discovery for conifer transcriptome sequencing.

There is a wide diversity of bioinformatic tools available for the assembly of next generation sequence and subsequence variant calling to identify genetic markers at scale. Integration of genomics tools such as genomic selection, association studies, pedigree analysis and analysis of genetic diversity, into operational breeding is a goal for New Zealand's most widely planted exotic tree species, Pinus radiata. In the absence of full reference genomes for large megagenomes such as in conifers, RNA sequencing in a range of genotypes and tissue types, offers a rich source of genetic markers for downstream application. We compared nine different assembler and variant calling software combinations in a single transcriptomic library and found that Single Nucleotide Polymorphism (SNPs) discovery could vary by as much as an order of magnitude (8,061 SNPs up to 86,815 SNPs). The assembler with the best realignment of the packages trialled, Trinity, in combination with several variant callers was then applied to a much larger multi-genotype, multi-tissue transcriptome and identified 683,135 in silico SNPs across a predicted 449,951 exons when mapped to the Pinus taeda ver 1.01e reference.

Phenotyping Whole Forests Will Help to Track Genetic Performance.

Phenotyping is the accurate and precise physical description of organisms. Accurate and quantitative phenotyping underpins the delivery of benefits from genetic improvement programs in agriculture. In forest trees, phenotyping at an equivalent precision has been impossible because trees and forests are large, long-lived, and highly variable. These facts have restricted the delivery of genetic gains in forestry compared to other agricultural sectors. We describe a landscape-scale phenotyping platform that integrates remote sensing, spatial information systems, and genomics to facilitate the delivery of greater gains enabling forestry to catch up with other sectors. Combining remote sensing at a range of spatial and temporal scales with genomics will ultimately impact on tree breeding globally.

Exploration of genetic architecture through sib-ship reconstruction in advanced breeding population of Eucalyptus nitens.

Accurate inference of relatedness between individuals in breeding population contributes to the precision of genetic parameter estimates, effectiveness of inbreeding management and the amount of genetic progress delivered from breeding programs. Pedigree reconstruction has been proven to be an efficient tool to correct pedigree errors and recover hidden relatedness in open pollinated progeny tests but the method can be limited by the lack of parental genotypes and the high proportion of alien pollen from outside the breeding population. Our study investigates the efficiency of sib-ship reconstruction in an advanced breeding population of Eucalyptus nitens with only partially tracked pedigree. The sib-ship reconstruction allowed the identification of selfs (4% of the sample) and the exploration of their potential effect on inbreeding depression in the traits studied. We detected signs of inbreeding depression in diameter at breast height and growth strain while no indications were observed in wood density, wood stiffness and tangential air-dry shrinkage. After the application of a corrected sib-ship relationship matrix, additive genetic variance and heritability were observed to increase where signs of inbreeding depression were initially detected. Conversely, the same genetic parameters for traits that appeared to be free of inbreeding depression decreased in size. It therefore appeared that greater genetic variance may be due, at least in part, to contributions from inbreeding in these studied populations rather than a removal of inbreeding as is traditionally thought.

Parentage Reconstruction in Eucalyptus nitens Using SNPs and Microsatellite Markers: A Comparative Analysis of Marker Data Power and Robustness.

Pedigree reconstruction using molecular markers enables efficient management of inbreeding in open-pollinated breeding strategies, replacing expensive and time-consuming controlled pollination. This is particularly useful in preferentially outcrossed, insect pollinated Eucalypts known to suffer considerable inbreeding depression from related matings. A single nucleotide polymorphism (SNP) marker panel consisting of 106 markers was selected for pedigree reconstruction from the recently developed high-density Eucalyptus Infinium SNP chip (EuCHIP60K). The performance of this SNP panel for pedigree reconstruction in open-pollinated progenies of two Eucalyptus nitens seed orchards was compared with that of two microsatellite panels with 13 and 16 markers respectively. The SNP marker panel out-performed one of the microsatellite panels in the resolution power to reconstruct pedigrees and out-performed both panels with respect to data quality. Parentage of all but one offspring in each clonal seed orchard was correctly matched to the expected seed parent using the SNP marker panel, whereas parentage assignment to less than a third of the expected seed parents were supported using the 13-microsatellite panel. The 16-microsatellite panel supported all but one of the recorded seed parents, one better than the SNP panel, although there was still a considerable level of missing and inconsistent data. SNP marker data was considerably superior to microsatellite data in accuracy, reproducibility and robustness. Although microsatellites and SNPs data provide equivalent resolution for pedigree reconstruction, microsatellite analysis requires more time and experience to deal with the uncertainties of allele calling and faces challenges for data transferability across labs and over time. While microsatellite analysis will continue to be useful for some breeding tasks due to the high information content, existing infrastructure and low operating costs, the multi-species SNP resource available with the EuCHIP60k, opens a whole new array of opportunities for high-throughput, genome-wide or targeted genotyping in species of Eucalyptus.