Factors affecting the accuracy of genomic selection for growth and wood quality traits in an advanced-breeding population of black spruce (Picea mariana)

Integrating evolutionary genomics of forest trees to inform future tree breeding amid rapid climate change

Global climate change is leading to rapid and drastic shifts in environmental conditions, posing threats to biodiversity and nearly all life forms worldwide. Forest trees serve as foundational components of terrestrial ecosystems and play a crucial and leading role in combating and mitigating the adverse effects of extreme climate events, despite their own vulnerability to these threats. Therefore, understanding and monitoring how natural forests respond to rapid climate change is a key priority for biodiversity conservation. Recent progress in evolutionary genomics, driven primarily by cutting-edge multi-omics technologies, offers powerful new tools to address several key issues. These include precise delineation of species and evolutionary units, inference of past evolutionary histories and demographic fluctuations, identification of environmentally adaptive variants, and measurement of genetic load levels. As the urgency to deal with more extreme environmental stresses grows, understanding the genomics of evolutionary history, local adaptation, future responses to climate change, and conservation and restoration of natural forest trees will be critical for research at the nexus of global change, population genomics, and conservation biology. In this review, we explore the application of evolutionary genomics to assess the effects of global climate change using multi-omics approaches and discuss the outlook for breeding of climate-adapted trees.

GWAS elucidated grain yield genetics in Indian spring wheat under diverse water conditions

KEY MESSAGE

Underpinned natural variations and key genes associated with yield under different water regimes, and identified genomic signatures of genetic gain in the Indian wheat breeding program. A novel KASP marker for TKW under water stress was developed and validated. A comprehensive genome-wide association study was conducted on 300 spring wheat genotypes to elucidate the natural variations associated with grain yield and its eleven contributing traits under fully irrigated, restricted water, and simulated no water conditions. Utilizing the 35K Wheat Breeders' Array, we identified 1155 quantitative trait nucleotides (QTNs), with 207 QTNs exhibiting stability across diverse conditions. These QTNs were further delimited into 539 genomic regions using a genome-wide LD value of 3.0 Mbp, revealing pleiotropic control across traits and conditions. Sub-genome A was significantly associated with traits under irrigated conditions, while sub-genome B showed more QTNs under water stressed conditions. Favourable alleles with significantly associated QTNs were delineated, with a notable pyramiding effect for enhancing trait performance. Additionally, allele of only 921 QTNs significantly affected the population mean. Allele profiling highlighted C-306 as a most potential source of drought tolerance. Moreover, 762 genes overlapping significant QTNs were identified, narrowing down to 27 putative candidate genes overlapping 29 novel and functional SNPs expressing (≥ 0.5 tpm) relevance across various growth conditions. A new KASP assay was developed, targeting a gene TraesCS2A03G1123700 regulating thousand kernel weight under severe drought condition. Genomic selection models (GBLUP, BayesB, MxE, and R-Norm) demonstrated an average prediction accuracy of 0.06-0.58 across environments, indicating potential for trait selection. Retrospective analysis of the Indian wheat breeding program supported a genetic gain in GY at the rate of ca. 0.56% per breeding cycle, since 1960, supporting the identification of genomic signatures driving trait selection and genetic gain. These findings offer insight into improving the rate of genetic gain in wheat breeding programs globally.

Genomic selection for growth characteristics in Korean red pine (Pinus densiflora Seibold & Zucc.)

Traditionally, selective breeding has been used to improve tree growth. However, traditional selection methods are time-consuming and limit annual genetic gain. Genomic selection (GS) offers an alternative to progeny testing by estimating the genotype-based breeding values of individuals based on genomic information using molecular markers. In the present study, we introduced GS to an open-pollinated breeding population of Korean red pine (Pinus densiflora), which is in high demand in South Korea, to shorten the breeding cycle. We compared the prediction accuracies of GS for growth characteristics (diameter at breast height [DBH], height, straightness, and volume) in Korean red pines under various conditions (marker set, model, and training set) and evaluated the selection efficiency of GS compared to traditional selection methods. Training the GS model to include individuals from various environments using genomic best linear unbiased prediction (GBLUP) and markers with a minor allele frequency larger than 0.05 was effective. The optimized model had an accuracy of 0.164-0.498 and a predictive ability of 0.018-0.441. The predictive ability of GBLUP against that of additive best linear unbiased prediction (ABLUP) was 0.86-5.10, and against the square root of heritability was 0.19-0.76, indicating that GS for Korean red pine was as efficient as in previous studies on forest trees. Moreover, the response to GS was higher than that to traditional selection regarding the annual genetic gain. Therefore, we conclude that the trained GS model is more effective than the traditional breeding methods for Korean red pines. We anticipate that the next generation of trees selected by GS will lay the foundation for the accelerated breeding of Korean red pine.

Assembly and annotation of the black spruce genome provide insights on spruce phylogeny and evolution of stress response

Black spruce (Picea mariana [Mill.] B.S.P.) is a dominant conifer species in the North American boreal forest that plays important ecological and economic roles. Here, we present the first genome assembly of P. mariana with a reconstructed genome size of 18.3 Gbp and NG50 scaffold length of 36.0 kbp. A total of 66,332 protein-coding sequences were predicted in silico and annotated based on sequence homology. We analyzed the evolutionary relationships between P. mariana and 5 other spruces for which complete nuclear and organelle genome sequences were available. The phylogenetic tree estimated from mitochondrial genome sequences agrees with biogeography; specifically, P. mariana was strongly supported as a sister lineage to P. glauca and 3 other taxa found in western North America, followed by the European Picea abies. We obtained mixed topologies with weaker statistical support in phylogenetic trees estimated from nuclear and chloroplast genome sequences, indicative of ancient reticulate evolution affecting these 2 genomes. Clustering of protein-coding sequences from the 6 Picea taxa and 2 Pinus species resulted in 34,776 orthogroups, 560 of which appeared to be specific to P. mariana. Analysis of these specific orthogroups and dN/dS analysis of positive selection signatures for 497 single-copy orthogroups identified gene functions mostly related to plant development and stress response. The P. mariana genome assembly and annotation provides a valuable resource for forest genetics research and applications in this broadly distributed species, especially in relation to climate adaptation.

Increasing genomic prediction accuracy for unphenotyped full-sib families by modeling additive and dominance effects with large datasets in white spruce

Introduction

Genomic selection is becoming a standard technique in plant breeding and is now being introduced into forest tree breeding. Despite promising results to predict the genetic merit of superior material based on their additive breeding values, many studies and operational programs still neglect non-additive effects and their potential for enhancing genetic gains.

Methods

Using two large comprehensive datasets totaling 4,066 trees from 146 full-sib families of white spruce (Picea glauca (Moench) Voss), we evaluated the effect of the inclusion of dominance on the precision of genetic parameter estimates and on the accuracy of conventional pedigree-based (ABLUP-AD) and genomic-based (GBLUP-AD) models.

Results

While wood quality traits were mostly additively inherited, considerable non-additive effects and lower heritabilities were detected for growth traits. For growth, GBLUP-AD better partitioned the additive and dominance effects into roughly equal variances, while ABLUP-AD strongly overestimated dominance. The predictive abilities of breeding and total genetic value estimates were similar between ABLUP-AD and GBLUP-AD when predicting individuals from the same families as those included in the training dataset. However, GBLUP-AD outperformed ABLUP-AD when predicting for new unphenotyped families that were not represented in the training dataset, with, on average, 22% and 53% higher predictive ability of breeding and genetic values, respectively. Resampling simulations showed that GBLUP-AD required smaller sample sizes than ABLUP-AD to produce precise estimates of genetic variances and accurate predictions of genetic values. Still, regardless of the method used, large training datasets were needed to estimate additive and non-additive genetic variances precisely.

Discussion

This study highlights the different quantitative genetic architectures between growth and wood traits. Furthermore, the usefulness of genomic additive-dominance models for predicting new families should allow practicing mating allocation to maximize the total genetic values for the propagation of elite material.

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.

Genomic Prediction of Complex Traits in Perennial Plants: A Case for Forest Trees

This chapter provides an overview of the genomic selection progress in long-lived forest tree species. Factors affecting the prediction accuracy in genomic prediction are assessed with examples from empirical studies. Infrastructure and resources required for the implementation of genomic selection are evaluated. Some general guidelines are provided for the successful application of genomic selection in forest tree breeding programs.

A bulked segregant analysis tool for out-crossing species (BSATOS) and QTL-based genomics-assisted prediction of complex traits in apple

INTRODUCTION

Genomic heterozygosity, self-incompatibility, and rich-in somatic mutations hinder the molecular breeding efficiency of outcrossing plants.

OBJECTIVES

We attempted to develop an efficient integrated strategy to identify quantitative trait loci (QTLs) and trait-associated genes, to develop gene markers, and to construct genomics-assisted prediction (GAP) modes.

METHODS

A novel protocol, bulked segregant analysis tool for out-crossing species (BSATOS), is presented here, which is characterized by taking full advantage of all segregation patterns (including AB × AB markers) and haplotype information. To verify the effectiveness of the protocol in dealing with the complex traits of outbreeding species, three apple cross populations with 9,654 individuals were adopted.

RESULTS

By using BSATOS, 90, 60, and 77 significant QTLs were identified successfully and candidate genes were predicted for apple fruit weight (FW), fruit ripening date (FRD), and fruit soluble solid content (SSC), respectively. The gene-based markers were developed and genotyped for 1,396 individuals in a training population, including 145 Malus accessions and 1,251 F1 plants of the three full-sib families. GAP models were trained using marker genotype effect estimates of the training population. The prediction accuracy was 0.7658, 0.6455, and 0.3758 for FW, FRD, and SSC, respectively.

CONCLUSION

The BSATOS and GAP models provided a convenient and efficient methodology for candidate gene mining and molecular breeding in out-crossing plant species. The BSATOS pipeline can be freely downloaded from: https://github.com/maypoleflyn/BSATOS.

Integrating genomic information and productivity and climate-adaptability traits into a regional white spruce breeding program

Tree improvement programs often focus on improving productivity-related traits; however, under present climate change scenarios, climate change-related (adaptive) traits should also be incorporated into such programs. Therefore, quantifying the genetic variation and correlations among productivity and adaptability traits, and the importance of genotype by environment interactions, including defense compounds involved in biotic and abiotic resistance, is essential for selecting parents for the production of resilient and sustainable forests. Here, we estimated quantitative genetic parameters for 15 growth, wood quality, drought resilience, and monoterpene traits for Picea glauca (Moench) Voss (white spruce). We sampled 1,540 trees from three open-pollinated progeny trials, genotyped with 467,224 SNP markers using genotyping-by-sequencing (GBS). We used the pedigree and SNP information to calculate, respectively, the average numerator and genomic relationship matrices, and univariate and multivariate individual-tree models to obtain estimates of (co)variance components. With few site-specific exceptions, all traits examined were under genetic control. Overall, higher heritability estimates were derived from the genomic- than their counterpart pedigree-based relationship matrix. Selection for height, generally, improved diameter and water use efficiency, but decreased wood density, microfibril angle, and drought resistance. Genome-based correlations between traits reaffirmed the pedigree-based correlations for most trait pairs. High and positive genetic correlations between sites were observed (average 0.68), except for those pairs involving the highest elevation, warmer, and moister site, specifically for growth and microfibril angle. These results illustrate the advantage of using genomic information jointly with productivity and adaptability traits, and defense compounds to enhance tree breeding selection for changing climate.

Metadata analysis indicates biased estimation of genetic parameters and gains using conventional pedigree information instead of genomic-based approaches in tree breeding

Forest tree improvement helps provide adapted planting stock to ensure growth productivity, fibre quality and carbon sequestration through reforestation and afforestation activities. However, there is increasing doubt that conventional pedigree provides the most accurate estimates for selection and prediction of performance of improved planting stock. When the additive genetic relationships among relatives is estimated using pedigree information, it is not possible to take account of Mendelian sampling due to the random segregation of parental alleles. The use of DNA markers distributed genome-wide (multi-locus genotypes) makes it possible to estimate the realized additive genomic relationships, which takes account of the Mendelian sampling and possible pedigree errors. We reviewed a series of papers on conifer and broadleaf tree species in which both pedigree-based and marker-based estimates of genetic parameters have been reported. Using metadata analyses, we show that for heritability and genetic gains, the estimates obtained using only the pedigree information are generally biased upward compared to those obtained using DNA markers distributed genome-wide, and that genotype-by-environment (GxE) interaction can be underestimated for low to moderate heritability traits. As high-throughput genotyping becomes economically affordable, we recommend expanding the use of genomic selection to obtain more accurate estimates of genetic parameters and gains.

Breeding for adaptation to climate change: genomic selection for drought response in a white spruce multi-site polycross test

With climate change, increasingly intense and frequent drought episodes will be affecting water availability for boreal tree species, prompting tree breeders and forest managers to consider adaptation to drought stress as a priority in their reforestation efforts. We used a 19-year-old polycross progeny test of the model conifer white spruce (Picea glauca) replicated on two sites affected by distinct drought episodes at different ages to estimate the genetic control and the potential for improvement of drought response in addition to conventional cumulative growth and wood quality traits. Drought response components were measured from dendrochronological signatures matching drought episodes in wood ring increment cores. We found that trees with more vigorous growth during their lifespan resisted better during the current year of a drought episode when the drought had more severe effects. Phenotypic data were also analyzed using genomic prediction (GBLUP) relying on the genomic relationship matrix of multi-locus gene SNP marker information, and conventional analysis (ABLUP) based on validated pedigree information. The accuracy of predicted breeding values for drought response components was marginally lower than that for conventional traits and comparable between GBLUP and ABLUP. Genetic correlations were generally low and nonsignificant between drought response components and conventional traits, except for resistance which was positively correlated to tree height. Heritability estimates for the components of drought response were slightly lower than for conventional traits, but similar single-trait genetic gains could be obtained. Multi-trait genomic selection simulations indicated that it was possible to improve simultaneously for all traits on both sites while sacrificing little on gain in tree height. In a context of rapid climate change, our results suggest that with careful phenotypic assessment, drought response may be considered in multi-trait improvement of white spruce, with accelerated screening of large numbers of candidates and selection at young age with genomic selection.

Improving lodgepole pine genomic evaluation using spatial correlation structure and SNP selection with single-step GBLUP

Modeling environmental spatial heterogeneity can improve the efficiency of forest tree genomic evaluation. Furthermore, genotyping costs can be lowered by reducing the number of markers needed. We investigated the impact on variance components, breeding value accuracy, and bias of two phenotypic data adjustments (experimental design and autoregressive spatial models), and a relationship matrix calculated from a subset of markers selected for their ability to infer ancestry. Using a multiple-trait multiple-site single-step Genomic Best Linear Unbiased Prediction (ssGBLUP) approach, four scenarios (2 phenotype adjustments × 2 marker sets) were applied to diameter at breast height (DBH), height (HT), and resistance to western gall rust (WGR) in four open-pollinated progeny trials of lodgepole pine, with 1490 (out of 11,188) trees genotyped with 25,099 SNPs. As a control, we fitted the conventional ABLUP model using pedigree information. The highest heritability estimates were achieved for the ABLUP followed closely by the ssGBLUP with the full marker set and using the spatial phenotype adjustments. The highest predictive ability was obtained by using a reduced marker subset (8000 SNPs) when either the spatial (DBH: 0.429, and WGR: 0.513) or design (HT: 0.467) phenotype corrections were used. No significant difference was detected in prediction bias among the six fitted models, and all values were close to 1 (0.918-1.014). Results demonstrated that selecting informative markers, such as those capturing ancestry, can improve the predictive ability. The use of spatial correlation structure increased traits' heritability and reduced prediction bias, while increases in predictive ability were trait-dependent.

Genomics-Enabled Management of Genetic Resources in Radiata Pine

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.

Prediction ability of genome-wide markers in Pinus taeda L. within and between population is affected by relatedness to the training population and trait genetic architecture

Genomic prediction has the potential to significantly increase the rate of genetic gain in tree breeding programs. In this study, a clonally replicated population (n = 2063) was used to train a genomic prediction model. The model was validated both within the training population and in a separate population (n = 451). The prediction abilities from random (20% vs 80%) cross validation within the training population were 0.56 and 0.78 for height and stem form, respectively. Removal of all full-sib relatives within the training population resulted in ∼50% reduction in their genomic prediction ability for both traits. The average prediction ability for all 451 individual trees was 0.29 for height and 0.57 for stem form. The degree of genetic linkage (full-sib family, half sib family, unrelated) between the training and validation sets had a strong impact on prediction ability for stem form but not for height. A dominant dwarfing allele, the first to be reported in a conifer species, was discovered via genome-wide association studies on linkage Group 5 that conferred a 0.33-m mean height reduction. However, the QTL was family specific. The rapid decay of linkage disequilibrium, large genome size, and inconsistencies in marker-QTL linkage phase suggest that large, diverse training populations are needed for genomic selection in Pinus taeda L.

A Semi-Automated SNP-Based Approach for Contaminant Identification in Biparental Polyploid Populations of Tropical Forage Grasses

Artificial hybridization plays a fundamental role in plant breeding programs since it generates new genotypic combinations that can result in desirable phenotypes. Depending on the species and mode of reproduction, controlled crosses may be challenging, and contaminating individuals can be introduced accidentally. In this context, the identification of such contaminants is important to avoid compromising further selection cycles, as well as genetic and genomic studies. The main objective of this work was to propose an automated multivariate methodology for the detection and classification of putative contaminants, including apomictic clones (ACs), self-fertilized individuals, half-siblings (HSs), and full contaminants (FCs), in biparental polyploid progenies of tropical forage grasses. We established a pipeline to identify contaminants in genotyping-by-sequencing (GBS) data encoded as allele dosages of single nucleotide polymorphism (SNP) markers by integrating principal component analysis (PCA), genotypic analysis (GA) measures based on Mendelian segregation, and clustering analysis (CA). The combination of these methods allowed for the correct identification of all contaminants in all simulated progenies and the detection of putative contaminants in three real progenies of tropical forage grasses, providing an easy and promising methodology for the identification of contaminants in biparental progenies of tetraploid and hexaploid species. The proposed pipeline was made available through the polyCID Shiny app and can be easily coupled with traditional genetic approaches, such as linkage map construction, thereby increasing the efficiency of breeding programs.

Achievements and Challenges of Genomics-Assisted Breeding in Forest Trees: From Marker-Assisted Selection to Genome Editing

Forest tree breeding efforts have focused mainly on improving traits of economic importance, selecting trees suited to new environments or generating trees that are more resilient to biotic and abiotic stressors. This review describes various methods of forest tree selection assisted by genomics and the main technological challenges and achievements in research at the genomic level. Due to the long rotation time of a forest plantation and the resulting long generation times necessary to complete a breeding cycle, the use of advanced techniques with traditional breeding have been necessary, allowing the use of more precise methods for determining the genetic architecture of traits of interest, such as genome-wide association studies (GWASs) and genomic selection (GS). In this sense, main factors that determine the accuracy of genomic prediction models are also addressed. In turn, the introduction of genome editing opens the door to new possibilities in forest trees and especially clustered regularly interspaced short palindromic repeats and CRISPR-associated protein 9 (CRISPR/Cas9). It is a highly efficient and effective genome editing technique that has been used to effectively implement targetable changes at specific places in the genome of a forest tree. In this sense, forest trees still lack a transformation method and an inefficient number of genotypes for CRISPR/Cas9. This challenge could be addressed with the use of the newly developing technique GRF-GIF with speed breeding.

Genome editing with CRISPR/Cas9 in Pinus radiata (D. Don)

BACKGROUND

To meet increasing demand for forest-based products and protect natural forests from further deforestation requires increased productivity from planted forests. Genetic improvement of conifers by traditional breeding is time consuming due to the long juvenile phase and genome complexity. Genetic modification (GM) offers the opportunity to make transformational changes in shorter time frames but is challenged by current genetically modified organism (GMO) regulations. Genome editing, which can be used to generate site-specific mutations, offers the opportunity to rapidly implement targeted improvements and is globally regulated in a less restrictive way than GM technologies.

RESULTS

We have demonstrated CRISPR/Cas9 genome editing in P. radiata targeting a single-copy cell wall gene GUX1 in somatic embryogenic tissue and produced plantlets from the edited tissue. We generated biallelic INDELs with an efficiency of 15 % using a single gRNA. 12 % of the transgenic embryogenic tissue was edited when two gRNAs were used and deletions of up to 1.3 kb were identified. However, the regenerated plants did not contain large deletions but had single nucleotide insertions at one of the target sites. We assessed the use of CRISPR/Cas9 ribonucleoproteins (RNPs) for their ability to accomplish DNA-free genome editing in P. radiata. We chose a hybrid approach, with RNPs co-delivered with a plasmid-based selectable marker. A two-gRNA strategy was used which produced an editing efficiency of 33 %, and generated INDELs, including large deletions. Using the RNP approach, deletions found in embryogenic tissue were also present in the plantlets. But, all plants produced using the RNP strategy were monoallelic.

CONCLUSIONS

We have demonstrated the generation of biallelic and monoallelic INDELs in the coniferous tree P. radiata with the CRISPR/Cas9 system using plasmid expressed Cas9 gRNA and RNPs respectively. This opens the opportunity to apply genome editing in conifers to rapidly modify key traits of interest.

Genomic Predictions With Nonadditive Effects Improved Estimates of Additive Effects and Predictions of Total Genetic Values in Pinus sylvestris

Genomic selection study (GS) focusing on nonadditive genetic effects of dominance and the first order of epistatic effects, in a full-sib family population of 695 Scots pine (Pinus sylvestris L.) trees, was undertaken for growth and wood quality traits, using 6,344 single nucleotide polymorphism markers (SNPs) generated by genotyping-by-sequencing (GBS). Genomic marker-based relationship matrices offer more effective modeling of nonadditive genetic effects than pedigree-based models, thus increasing the knowledge on the relevance of dominance and epistatic variation in forest tree breeding. Genomic marker-based models were compared with pedigree-based models showing a considerable dominance and epistatic variation for growth traits. Nonadditive genetic variation of epistatic nature (additive × additive) was detected for growth traits, wood density (DEN), and modulus of elasticity (MOEd) representing between 2.27 and 34.5% of the total phenotypic variance. Including dominance variance in pedigree-based Best Linear Unbiased Prediction (PBLUP) and epistatic variance in genomic-based Best Linear Unbiased Prediction (GBLUP) resulted in decreased narrow-sense heritability and increased broad-sense heritability for growth traits, DEN and MOEd. Higher genetic gains were reached with early GS based on total genetic values, than with conventional pedigree selection for a selection intensity of 1%. This study indicates that nonadditive genetic variance may have a significant role in the variation of selection traits of Scots pine, thus clonal deployment could be an attractive alternative for the species. Additionally, confidence in the role of nonadditive genetic effects in this breeding program should be pursued in the future, using GS.

Limited haplotype diversity underlies polygenic trait architecture across 70 years of wheat breeding

Background

Selection has dramatically shaped genetic and phenotypic variation in bread wheat. We can assess the genomic basis of historical phenotypic changes, and the potential for future improvement, using experimental populations that attempt to undo selection through the randomizing effects of recombination.

Results

We bred the NIAB Diverse MAGIC multi-parent population comprising over 500 recombinant inbred lines, descended from sixteen historical UK bread wheat varieties released between 1935 and 2004. We sequence the founders’ genes and promoters by capture, and the MAGIC population by low-coverage whole-genome sequencing. We impute 1.1 M high-quality SNPs that are over 99% concordant with array genotypes. Imputation accuracy only marginally improves when including the founders’ genomes as a haplotype reference panel. Despite capturing 73% of global wheat genetic polymorphism, 83% of genes cluster into no more than three haplotypes. We phenotype 47 agronomic traits over 2 years and map 136 genome-wide significant associations, concentrated at 42 genetic loci with large and often pleiotropic effects. Around half of these overlap known quantitative trait loci. Most traits exhibit extensive polygenicity, as revealed by multi-locus shrinkage modelling.

Conclusions

Our results are consistent with a gene pool of low haplotypic diversity, containing few novel loci of large effect. Most past, and projected future, phenotypic changes arising from existing variation involve fine-scale shuffling of a few haplotypes to recombine dozens of polygenic alleles of small effect. Moreover, extensive pleiotropy means selection on one trait will have unintended consequences, exemplified by the negative trade-off between yield and protein content, unless selection and recombination can break unfavorable trait-trait associations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-021-02354-7.

Factors Influencing Genomic Prediction Accuracies of Tropical Maize Resistance to Fall Armyworm and Weevils

Genomic selection (GS) can accelerate variety improvement when training set (TS) size and its relationship with the breeding set (BS) are optimized for prediction accuracies (PAs) of genomic prediction (GP) models. Sixteen GP algorithms were run on phenotypic best linear unbiased predictors (BLUPs) and estimators (BLUEs) of resistance to both fall armyworm (FAW) and maize weevil (MW) in a tropical maize panel. For MW resistance, 37% of the panel was the TS, and the BS was the remainder, whilst for FAW, random-based training sets (RBTS) and pedigree-based training sets (PBTSs) were designed. PAs achieved with BLUPs varied from 0.66 to 0.82 for MW-resistance traits, and for FAW resistance, 0.694 to 0.714 for RBTS of 37%, and 0.843 to 0.844 for RBTS of 85%, and these were at least two-fold those from BLUEs. For PBTS, FAW resistance PAs were generally higher than those for RBTS, except for one dataset. GP models generally showed similar PAs across individual traits whilst the TS designation was determinant, since a positive correlation (R = 0.92***) between TS size and PAs was observed for RBTS, and for the PBTS, it was negative (R = 0.44**). This study pioneered the use of GS for maize resistance to insect pests in sub-Saharan Africa.

Adoption and Optimization of Genomic Selection To Sustain Breeding for Apricot Fruit Quality

Genomic selection (GS) is a breeding approach which exploits genome-wide information and whose unprecedented success has shaped several animal and plant breeding schemes through delivering their genetic progress. This is the first study assessing the potential of GS in apricot (Prunus armeniaca) to enhance postharvest fruit quality attributes. Genomic predictions were based on a F1 pseudo-testcross population, comprising 153 individuals with contrasting fruit quality traits. They were phenotyped for physical and biochemical fruit metrics in contrasting climatic conditions over two years. Prediction accuracy (PA) varied from 0.31 for glucose content with the Bayesian LASSO (BL) to 0.78 for ethylene production with RR-BLUP, which yielded the most accurate predictions in comparison to Bayesian models and only 10% out of 61,030 SNPs were sufficient to reach accurate predictions. Useful insights were provided on the genetic architecture of apricot fruit quality whose integration in prediction models improved their performance, notably for traits governed by major QTL. Furthermore, multivariate modeling yielded promising outcomes in terms of PA within training partitions partially phenotyped for target traits. This provides a useful framework for the implementation of indirect selection based on easy-to-measure traits. Thus, we highlighted the main levers to take into account for the implementation of GS for fruit quality in apricot, but also to improve the genetic gain in perennial species.

Genomic selection for resistance to spruce budworm in white spruce and relationships with growth and wood quality traits

With climate change, the pressure on tree breeding to provide varieties with improved resilience to biotic and abiotic stress is increasing. As such, pest resistance is of high priority but has been neglected in most tree breeding programs, given the complexity of phenotyping for these traits and delays to assess mature trees. In addition, the existing genetic variation of resistance and its relationship with productivity should be better understood for their consideration in multitrait breeding. In this study, we evaluated the prospects for genetic improvement of the levels of acetophenone aglycones (AAs) in white spruce needles, which have been shown to be tightly linked to resistance to spruce budworm. Furthermore, we estimated the accuracy of genomic selection (GS) for these traits, allowing selection at a very early stage to accelerate breeding. A total of 1,516 progeny trees established on five sites and belonging to 136 full-sib families from a mature breeding population in New Brunswick were measured for height growth and genotyped for 4,148 high-quality SNPs belonging to as many genes along the white spruce genome. In addition, 598 trees were assessed for levels of AAs piceol and pungenol in needles, and 578 for wood stiffness. GS models were developed with the phenotyped trees and then applied to predict the trait values of unphenotyped trees. AAs were under moderate-to-high genetic control (h 2: 0.43-0.57) with null or marginally negative genetic correlations with other traits. The prediction accuracy of GS models (GBLUP) for AAs was high (PAAC: 0.63-0.67) and comparable or slightly higher than pedigree-based (ABLUP) or BayesCπ models. We show that AA traits can be improved and that GS speeds up the selection of improved trees for insect resistance and for growth and wood quality traits. Various selection strategies were tested to optimize multitrait gains.

Evaluation of the efficiency of genomic versus pedigree predictions for growth and wood quality traits in Scots pine

BACKGROUND

Genomic selection (GS) or genomic prediction is a promising approach for tree breeding to obtain higher genetic gains by shortening time of progeny testing in breeding programs. As proof-of-concept for Scots pine (Pinus sylvestris L.), a genomic prediction study was conducted with 694 individuals representing 183 full-sib families that were genotyped with genotyping-by-sequencing (GBS) and phenotyped for growth and wood quality traits. 8719 SNPs were used to compare different genomic with pedigree prediction models. Additionally, four prediction efficiency methods were used to evaluate the impact of genomic breeding value estimations by assigning diverse ratios of training and validation sets, as well as several subsets of SNP markers.

RESULTS

Genomic Best Linear Unbiased Prediction (GBLUP) and Bayesian Ridge Regression (BRR) combined with expectation maximization (EM) imputation algorithm showed slightly higher prediction efficiencies than Pedigree Best Linear Unbiased Prediction (PBLUP) and Bayesian LASSO, with some exceptions. A subset of approximately 6000 SNP markers, was enough to provide similar prediction efficiencies as the full set of 8719 markers. Additionally, prediction efficiencies of genomic models were enough to achieve a higher selection response, that varied between 50-143% higher than the traditional pedigree-based selection.

CONCLUSIONS

Although prediction efficiencies were similar for genomic and pedigree models, the relative selection response was doubled for genomic models by assuming that earlier selections can be done at the seedling stage, reducing the progeny testing time, thus shortening the breeding cycle length roughly by 50%.

Genomic Selection for Forest Tree Improvement: Methods, Achievements and Perspectives

The breeding of forest trees is only a few decades old, and is a much more complicated, longer, and expensive endeavor than the breeding of agricultural crops. One breeding cycle for forest trees can take 20–30 years. Recent advances in genomics and molecular biology have revolutionized traditional plant breeding based on visual phenotype assessment: the development of different types of molecular markers has made genotype selection possible. Marker-assisted breeding can significantly accelerate the breeding process, but this method has not been shown to be effective for selection of complex traits on forest trees. This new method of genomic selection is based on the analysis of all effects of quantitative trait loci (QTLs) using a large number of molecular markers distributed throughout the genome, which makes it possible to assess the genomic estimated breeding value (GEBV) of an individual. This approach is expected to be much more efficient for forest tree improvement than traditional breeding. Here, we review the current state of the art in the application of genomic selection in forest tree breeding and discuss different methods of genotyping and phenotyping. We also compare the accuracies of genomic prediction models and highlight the importance of a prior cost-benefit analysis before implementing genomic selection. Perspectives for the further development of this approach in forest breeding are also discussed: expanding the range of species and the list of valuable traits, the application of high-throughput phenotyping methods, and the possibility of using epigenetic variance to improve of forest trees.

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.

Favorable Conditions for Genomic Evaluation to Outperform Classical Pedigree Evaluation Highlighted by a Proof-of-Concept Study in Poplar

Forest trees like poplar are particular in many ways compared to other domesticated species. They have long juvenile phases, ongoing crop-wild gene flow, extensive outcrossing, and slow growth. All these particularities tend to make the conduction of breeding programs and evaluation stages costly both in time and resources. Perennials like trees are therefore good candidates for the implementation of genomic selection (GS) which is a good way to accelerate the breeding process, by unchaining selection from phenotypic evaluation without affecting precision. In this study, we tried to compare GS to pedigree-based traditional evaluation, and evaluated under which conditions genomic evaluation outperforms classical pedigree evaluation. Several conditions were evaluated as the constitution of the training population by cross-validation, the implementation of multi-trait, single trait, additive and non-additive models with different estimation methods (G-BLUP or weighted G-BLUP). Finally, the impact of the marker densification was tested through four marker density sets. The population under study corresponds to a pedigree of 24 parents and 1,011 offspring, structured into 35 full-sib families. Four evaluation batches were planted in the same location and seven traits were evaluated on 1 and 2 years old trees. The quality of prediction was reported by the accuracy, the Spearman rank correlation and prediction bias and tested with a cross-validation and an independent individual test set. Our results show that genomic evaluation performance could be comparable to the already well-optimized pedigree-based evaluation under certain conditions. Genomic evaluation appeared to be advantageous when using an independent test set and a set of less precise phenotypes. Genome-based methods showed advantages over pedigree counterparts when ranking candidates at the within-family levels, for most of the families. Our study also showed that looking at ranking criteria as Spearman rank correlation can reveal benefits to genomic selection hidden by biased predictions.

Whole-exome sequencing reveals a long-term decline in effective population size of red spruce (Picea rubens)

Understanding the factors influencing the current distribution of genetic diversity across a species range is one of the main questions of evolutionary biology, especially given the increasing threat to biodiversity posed by climate change. Historical demographic processes such as population expansion or bottlenecks and decline are known to exert a predominant influence on past and current levels of genetic diversity, and revealing this demo-genetic history can have immediate conservation implications. We used a whole-exome capture sequencing approach to analyze polymorphism across the gene space of red spruce (Picea rubens Sarg.), an endemic and emblematic tree species of eastern North America high-elevation forests that are facing the combined threat of global warming and increasing human activities. We sampled a total of 340 individuals, including populations from the current core of the range in northeastern USA and southeastern Canada and from the southern portions of its range along the Appalachian Mountains, where populations occur as highly fragmented mountaintop "sky islands." Exome capture baits were designed from the closely relative white spruce (P. glauca Voss) transcriptome, and sequencing successfully captured most regions on or near our target genes, resulting in the generation of a new and expansive genomic resource for studying standing genetic variation in red spruce applicable to its conservation. Our results, based on over 2 million exome-derived variants, indicate that red spruce is structured into three distinct ancestry groups that occupy different geographic regions of its highly fragmented range. Moreover, these groups show small Ne , with a temporal history of sustained population decline that has been ongoing for thousands (or even hundreds of thousands) of years. These results demonstrate the broad potential of genomic studies for revealing details of the demographic history that can inform management and conservation efforts of nonmodel species with active restoration programs, such as red spruce.

Genome-Wide Association Study towards Genomic Predictive Power for High Production and Quality of Milk in American Alpine Goats

This paper reports an exploratory study based on quantitative genomic analysis in dairy traits of American Alpine goats. The dairy traits are quality-determining components in goat milk, cheese, ice cream, etc. Alpine goat phenotypes for quality components have been routinely recorded for many years and deposited in the Council on Dairy Cattle Breeding (CDCB) repository. The data collected were used to conduct an exploratory genome-wide association study (GWAS) from 72 female Alpine goats originating from locations throughout the U.S. Genotypes were identified with the Illumina Goat 50K single-nucleotide polymorphisms (SNP) BeadChip. The analysis used a polygenic model where the dropping criterion was a call rate ≥ 0.95. The initial dataset was composed of ~60,000 rows of SNPs and 21 columns of phenotypic traits and composed of 53,384 scaffolds containing other informative data points used for genomic predictive power. Phenotypic association with the 50K BeadChip revealed 26,074 reads of candidate genes. These candidate genes segregated as separate novel SNPs and were identified as statistically significant regions for genome and chromosome level trait associations. Candidate genes associated differently for each of the following phenotypic traits: test day milk yield (13,469 candidate genes), test day protein yield (25,690 candidate genes), test day fat yield (25,690 candidate genes), percentage protein (25,690 candidate genes), percentage fat (25,690 candidate genes), and percentage lactose content (25,690 candidate genes). The outcome of this study supports elucidation of novel genes that are important for livestock species in association to key phenotypic traits. Validation towards the development of marker-based selection that provides precision breeding methods will thereby increase the breeding value.

Linkage disequilibrium vs. pedigree: Genomic selection prediction accuracy in conifer species

BACKGROUND

The presupposition of genomic selection (GS) is that predictive accuracies should be based on population-wide linkage disequilibrium (LD). However, in species with large, highly complex genomes the limitation of marker density may preclude the ability to resolve LD accurately enough for GS. Here we investigate such an effect in two conifer species with ~ 20 Gbp genomes, Douglas-fir (Pseudotsuga menziesii Mirb. (Franco)) and Interior spruce (Picea glauca (Moench) Voss x Picea engelmannii Parry ex Engelm.). Random sampling of markers was performed to obtain SNP sets with totals in the range of 200-50,000, this was replicated 10 times. Ridge Regression Best Linear Unbiased Predictor (RR-BLUP) was deployed as the GS method to test these SNP sets, and 10-fold cross-validation was performed on 1,321 Douglas-fir trees, representing 37 full-sib F1 families and on 1,126 Interior spruce trees, representing 25 open-pollinated (half-sib) families. Both trials are located on 3 sites in British Columbia, Canada.

RESULTS

As marker number increased, so did GS predictive accuracy for both conifer species. However, a plateau in the gain of accuracy became apparent around 10,000-15,000 markers for both Douglas-fir and Interior spruce. Despite random marker selection, little variation in predictive accuracy was observed across replications. On average, Douglas-fir prediction accuracies were higher than those of Interior spruce, reflecting the difference between full- and half-sib families for Douglas-fir and Interior spruce populations, respectively, as well as their respective effective population size.

CONCLUSIONS

Although possibly advantageous within an advanced breeding population, reducing marker density cannot be recommended for carrying out GS in conifers. Significant LD between markers and putative causal variants was not detected using 50,000 SNPS, and GS was enabled only through the tracking of relatedness in the populations studied. Dramatically increasing marker density would enable said markers to better track LD with causal variants in these large, genetically diverse genomes; as well as providing a model that could be used across populations, breeding programs, and 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.

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.

Effect of number of annual rings and tree ages on genomic predictive ability for solid wood properties of Norway spruce

Background

Genomic selection (GS) or genomic prediction is considered as a promising approach to accelerate tree breeding and increase genetic gain by shortening breeding cycle, but the efforts to develop routines for operational breeding are so far limited. We investigated the predictive ability (PA) of GS based on 484 progeny trees from 62 half-sib families in Norway spruce (Picea abies (L.) Karst.) for wood density, modulus of elasticity (MOE) and microfibril angle (MFA) measured with SilviScan, as well as for measurements on standing trees by Pilodyn and Hitman instruments.

Results

GS predictive abilities were comparable with those based on pedigree-based prediction. Marker-based PAs were generally 25–30% higher for traits density, MFA and MOE measured with SilviScan than for their respective standing tree-based method which measured with Pilodyn and Hitman. Prediction accuracy (PC) of the standing tree-based methods were similar or even higher than increment core-based method. 78–95% of the maximal PAs of density, MFA and MOE obtained from coring to the pith at high age were reached by using data possible to obtain by drilling 3–5 rings towards the pith at tree age 10–12.

Conclusions

This study indicates standing tree-based measurements is a cost-effective alternative method for GS. PA of GS methods were comparable with those pedigree-based prediction. The highest PAs were reached with at least 80–90% of the dataset used as training set. Selection for trait density could be conducted at an earlier age than for MFA and MOE. Operational breeding can also be optimized by training the model at an earlier age or using 3 to 5 outermost rings at tree age 10 to 12 years, thereby shortening the cycle and reducing the impact on the tree.

Potential of Genome-Wide Association Studies and Genomic Selection to Improve Productivity and Quality of Commercial Timber Species in Tropical Rainforest, a Case Study of Shorea platyclados

Shorea platyclados (Dark Red Meranti) is a commercially important timber tree species in Southeast Asia. However, its stocks have dramatically declined due, inter alia, to excessive logging, insufficient natural regeneration and a slow recovery rate. Thus, there is a need to promote enrichment planting and develop effective technique to support its rehabilitation and improve timber production through implementation of Genome-Wide Association Studies (GWAS) and Genomic Selection (GS). To assist such efforts, plant materials were collected from a half-sib progeny population in Sari Bumi Kusuma forest concession, Kalimantan, Indonesia. Using 5900 markers in sequences obtained from 356 individuals, we detected high linkage disequilibrium (LD) extending up to >145 kb, suggesting that associations between phenotypic traits and markers in LD can be more easily and feasibly detected with GWAS than with analysis of quantitative trait loci (QTLs). However, the detection power of GWAS seems low, since few single nucleotide polymorphisms linked to any focal traits were detected with a stringent false discovery rate, indicating that the species’ phenotypic traits are mostly under polygenic quantitative control. Furthermore, Machine Learning provided higher prediction accuracies than Bayesian methods. We also found that stem diameter, branch diameter ratio and wood density were more predictable than height, clear bole, branch angle and wood stiffness traits. Our study suggests that GS has potential for improving the productivity and quality of S. platyclados, and our genomic heritability estimates may improve the selection of traits to target in future breeding of this species.

Genomic prediction for hastening and improving efficiency of forward selection in conifer polycross mating designs: an example from white spruce

Genomic selection (GS) has a large potential for improving the prediction accuracy of breeding values and significantly reducing the length of breeding cycles. In this context, the choice of mating designs becomes critical to improve the efficiency of breeding operations and to obtain the largest genetic gains per time unit. Polycross mating designs have been traditionally used in tree and plant breeding to perform backward selection of the female parents. The possibility to use genetic markers for paternity identification and for building genomic prediction models should allow for a broader use of polycross tests in forward selection schemes. We compared the accuracies of genomic predictions of offspring's breeding values from a polycross and a full-sib (partial diallel) mating design with similar genetic background in white spruce (Picea glauca). Trees were phenotyped for growth and wood quality traits, and genotyped for 4092 SNPs representing as many gene loci distributed across the 12 spruce chromosomes. For the polycross progeny test, heritability estimates were smaller, but more precise using the genomic BLUP (GBLUP) model as compared with pedigree-based models accounting for the maternal pedigree or for the reconstructed full pedigree. Cross-validations showed that GBLUP predictions were 22-52% more accurate than predictions based on the maternal pedigree, and 5-7% more accurate than predictions using the reconstructed full pedigree. The accuracies of GBLUP predictions were high and in the same range for most traits between the polycross (0.61-0.70) and full-sib progeny tests (0.61-0.74). However, higher genetic gains per time unit were expected from the polycross mating design given the shorter time needed to conduct crosses. Considering the operational advantages of the polycross design in terms of easier handling of crosses and lower associated costs for test establishment, we believe that this mating scheme offers great opportunities for the development and operational application of forward GS.

Genomic Predictions Using Low-Density SNP Markers, Pedigree and GWAS Information: A Case Study with the Non-Model Species Eucalyptus cladocalyx

High-throughput genotyping techniques have enabled large-scale genomic analysis to precisely predict complex traits in many plant species. However, not all species can be well represented in commercial SNP (single nucleotide polymorphism) arrays. In this study, a high-density SNP array (60 K) developed for commercial Eucalyptus was used to genotype a breeding population of Eucalyptus cladocalyx, yielding only ~3.9 K informative SNPs. Traditional Bayesian genomic models were investigated to predict flowering, stem quality and growth traits by considering the following effects: (i) polygenic background and all informative markers (GS model) and (ii) polygenic background, QTL-genotype effects (determined by GWAS) and SNP markers that were not associated with any trait (GSq model). The estimates of pedigree-based heritability and genomic heritability varied from 0.08 to 0.34 and 0.002 to 0.5, respectively, whereas the predictive ability varied from 0.19 (GS) and 0.45 (GSq). The GSq approach outperformed GS models in terms of predictive ability when the proportion of the variance explained by the significant marker-trait associations was higher than those explained by the polygenic background and non-significant markers. This approach can be particularly useful for plant/tree species poorly represented in the high-density SNP arrays, developed for economically important species, or when high-density marker panels are not available.

275 years of forestry meets genomics in Pinus sylvestris

Pinus sylvestris has a long history of basic and applied research that is relevant for both forestry and evolutionary studies. Its patterns of adaptive variation and role in forest economic and ecological systems have been studied extensively for nearly 275 years, detailed demography for a 100 years and mating system more than 50 years. However, its reference genome sequence is not yet available and genomic studies have been lagging compared to, for example, Pinus taeda and Picea abies, two other economically important conifers. Despite the lack of reference genome, many modern genomic methods are applicable for a more detailed look at its biological characteristics. For example, RNA-seq has revealed a complex transcriptional landscape and targeted DNA sequencing displays an excess of rare variants and geographically homogenously distributed molecular genetic diversity. Current DNA and RNA resources can be used as a reference for gene expression studies, SNP discovery, and further targeted sequencing. In the future, specific consequences of the large genome size, such as functional effects of regulatory open chromatin regions and transposable elements, should be investigated more carefully. For forest breeding and long-term management purposes, genomic data can help in assessing the genetic basis of inbreeding depression and the application of genomic tools for genomic prediction and relatedness estimates. Given the challenges of breeding (long generation time, no easy vegetative propagation) and the economic importance, application of genomic tools has a potential to have a considerable impact. Here, we explore how genomic characteristics of P. sylvestris, such as rare alleles and the low extent of linkage disequilibrium, impact the applicability and power of the tools.

Multi-trait genomic selection for weevil resistance, growth, and wood quality in Norway spruce

Plantation-grown trees have to cope with an increasing pressure of pest and disease in the context of climate change, and breeding approaches using genomics may offer efficient and flexible tools to face this pressure. In the present study, we targeted genetic improvement of resistance of an introduced conifer species in Canada, Norway spruce (Picea abies (L.) Karst.), to the native white pine weevil (Pissodes strobi Peck). We developed single- and multi-trait genomic selection (GS) models and selection indices considering the relationships between weevil resistance, intrinsic wood quality, and growth traits. Weevil resistance, acoustic velocity as a proxy for mechanical wood stiffness, and average wood density showed moderate-to-high heritability and low genotype-by-environment interactions. Weevil resistance was genetically positively correlated with tree height, height-to-diameter at breast height (DBH) ratio, and acoustic velocity. The accuracy of the different GS models tested (GBLUP, threshold GBLUP, Bayesian ridge regression, BayesCπ) was high and did not differ among each other. Multi-trait models performed similarly as single-trait models when all trees were phenotyped. However, when weevil attack data were not available for all trees, weevil resistance was more accurately predicted by integrating genetically correlated growth traits into multi-trait GS models. A GS index that corresponded to the breeders' priorities achieved near maximum gains for weevil resistance, acoustic velocity, and height growth, but a small decrease for DBH. The results of this study indicate that it is possible to breed for high-quality, weevil-resistant Norway spruce reforestation stock with high accuracy achieved from single-trait or multi-trait GS.

Enhancing genomic selection by fitting large-effect SNPs as fixed effects and a genotype-by-environment effect using a maize BC1F3:4 population

The popularity of genomic selection (GS) has increased owing to its prospects in commercial breeding. It is necessary to enhance GS to increase its efficiency. In this study, a maize BC₁F_3:4 population, consisting of 481 families, was evaluated for days to anthesis in four environments, and genotyped with DNA chips including 55,000 single nucleotide polymorphisms (SNPs). This population was used to investigate whether GS could be enhanced by borrowing information from the genetic basis and genotype-by-environment (G × E) interaction. The results showed that: 1) fitting the top four large-effect SNPs as fixed effects could increase prediction accuracy, including three minor-effect SNPs explaining less than 10% phenotypic variance; 2) the increase of prediction accuracy when fitting large-effect SNPs as fixed effects was related to the decrease of genetic variance; 3) generally, the GS model fitting large-effect SNPs as fixed effects and G × E component enhanced GS. Therefore, we propose fitting large-effect markers as fixed effects and G × E effect for crop breeding projects in order to obtain accurately predicted phenotypic data and conduct efficient selection of desired plants.

SNP and Haplotype-Based Genomic Selection of Quantitative Traits in Eucalyptus globulus

Eucalyptus globulus (Labill.) is one of the most important cultivated eucalypts in temperate and subtropical regions and has been successfully subjected to intensive breeding. In this study, Bayesian genomic models that include the effects of haplotype and single nucleotide polymorphisms (SNP) were assessed to predict quantitative traits related to wood quality and tree growth in a 6-year-old breeding population. To this end, the following markers were considered: (a) ~14 K SNP markers (SNP), (b) ~3 K haplotypes (HAP), and (c) haplotypes and SNPs that were not assigned to a haplotype (HAP-SNP). Predictive ability values (PA) were dependent on the genomic prediction models and markers. On average, Bayesian ridge regression (BRR) and Bayes C had the highest PA for the majority of traits. Notably, genomic models that included the haplotype effect (either HAP or HAP-SNP) significantly increased the PA of low-heritability traits. For instance, BRR based on HAP had the highest PA (0.58) for stem straightness. Consistently, the heritability estimates from genomic models were higher than the pedigree-based estimates for these traits. The results provide additional perspectives for the implementation of genomic selection in Eucalyptus breeding programs, which could be especially beneficial for improving traits with low heritability.

Genome-Wide SNP Data Revealed the Extent of Linkage Disequilibrium, Persistence of Phase and Effective Population Size in Purebred and Crossbred Buffalo Populations

Linkage disequilibrium (LD) is a useful parameter for guiding the accuracy and power of both genome-wide association studies (GWAS) and genomic selection (GS) among different livestock species. The present study evaluated the extent of LD, persistence of phase and effective population size (Ne) for the purebred (Mediterranean buffalo; n = 411) and crossbred [Mediterranean × Jianghan × Nili-Ravi buffalo, n = 9; Murrah × Nili-Ravi × local (Xilin or Fuzhong) buffalo, n = 36] buffalo populations using the 90K Buffalo SNP genotyping array. The results showed that the average square of correlation coefficient (r ²) between adjacent SNP was 0.13 ± 0.19 across all autosomes for purebred and 0.09 ± 0.13 for crossbred, and the most rapid decline in LD was observed over the first 200 kb. Estimated r ² ≥ 0.2 extended up to ~50 kb in crossbred and 170 kb in purebred populations, while average r ² values ≥0.3 were respectively observed in the ~10 and 60 kb in the crossbred and purebred populations. The largest phase correlation (R _{P, C} = 0.47) was observed at the distance of 100 kb, suggesting that this phase was not actively preserved between the two populations. Estimated Ne for the purebred and crossbred population at the current generation was 387 and 113 individuals, respectively. These findings may provide useful information to guide the GS and GWAS in buffaloes.

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

BACKGROUND

Genomic selection (GS) can increase genetic gain by reducing the length of breeding cycle in forest trees. Here we genotyped 1370 control-pollinated progeny trees from 128 full-sib families in Norway spruce (Picea abies (L.) Karst.), using exome capture as genotyping platform. We used 116,765 high-quality SNPs to develop genomic prediction models for tree height and wood quality traits. We assessed the impact of different genomic prediction methods, genotype-by-environment interaction (G × E), genetic composition, size of the training and validation set, relatedness, and number of SNPs on accuracy and predictive ability (PA) of GS.

RESULTS

Using G matrix slightly altered heritability estimates relative to pedigree-based method. GS accuracies were about 11-14% lower than those based on pedigree-based selection. The efficiency of GS per year varied from 1.71 to 1.78, compared to that of the pedigree-based model if breeding cycle length was halved using GS. Height GS accuracy decreased to more than 30% while using one site as training for GS prediction and using this model to predict the second site, indicating that G × E for tree height should be accommodated in model fitting. Using a half-sib family structure instead of full-sib structure led to a significant reduction in GS accuracy and PA. The full-sib family structure needed only 750 markers to reach similar accuracy and PA, as compared to 100,000 markers required for the half-sib family, indicating that maintaining the high relatedness in the model improves accuracy and PA. Using 4000-8000 markers in full-sib family structure was sufficient to obtain GS model accuracy and PA for tree height and wood quality traits, almost equivalent to that obtained with all markers.

CONCLUSIONS

The study indicates that GS would be efficient in reducing generation time of breeding cycle in conifer tree breeding program that requires long-term progeny testing. The sufficient number of trees within-family (16 for growth and 12 for wood quality traits) and number of SNPs (8000) are required for GS with full-sib family relationship. GS methods had little impact on GS efficiency for growth and wood quality traits. GS model should incorporate G × E effect when a strong G × E is detected.

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

BACKGROUND

Genomic selection (GS) can increase genetic gain by reducing the length of breeding cycle in forest trees. Here we genotyped 1370 control-pollinated progeny trees from 128 full-sib families in Norway spruce (Picea abies (L.) Karst.), using exome capture as genotyping platform. We used 116,765 high-quality SNPs to develop genomic prediction models for tree height and wood quality traits. We assessed the impact of different genomic prediction methods, genotype-by-environment interaction (G × E), genetic composition, size of the training and validation set, relatedness, and number of SNPs on accuracy and predictive ability (PA) of GS.

RESULTS

Using G matrix slightly altered heritability estimates relative to pedigree-based method. GS accuracies were about 11-14% lower than those based on pedigree-based selection. The efficiency of GS per year varied from 1.71 to 1.78, compared to that of the pedigree-based model if breeding cycle length was halved using GS. Height GS accuracy decreased to more than 30% while using one site as training for GS prediction and using this model to predict the second site, indicating that G × E for tree height should be accommodated in model fitting. Using a half-sib family structure instead of full-sib structure led to a significant reduction in GS accuracy and PA. The full-sib family structure needed only 750 markers to reach similar accuracy and PA, as compared to 100,000 markers required for the half-sib family, indicating that maintaining the high relatedness in the model improves accuracy and PA. Using 4000-8000 markers in full-sib family structure was sufficient to obtain GS model accuracy and PA for tree height and wood quality traits, almost equivalent to that obtained with all markers.

CONCLUSIONS

The study indicates that GS would be efficient in reducing generation time of breeding cycle in conifer tree breeding program that requires long-term progeny testing. The sufficient number of trees within-family (16 for growth and 12 for wood quality traits) and number of SNPs (8000) are required for GS with full-sib family relationship. GS methods had little impact on GS efficiency for growth and wood quality traits. GS model should incorporate G × E effect when a strong G × E is detected.

A catalog of annotated high-confidence SNPs from exome capture and sequencing reveals highly polymorphic genes in Norway spruce (Picea abies)

BACKGROUND

Norway spruce [Picea abies (L.) Karst.] is ecologically and economically one of the most important conifer worldwide. Our main goal was to develop a large catalog of annotated high confidence gene SNPs that should sustain the development of genomic tools for the conservation of natural and domesticated genetic diversity resources, and hasten tree breeding efforts in this species.

RESULTS

Targeted sequencing was achieved by capturing P. abies exome with probes previously designed from the sequenced transcriptome of white spruce (Picea glauca (Moench) Voss). Capture efficiency was high (74.5%) given a high level of exome conservation between the two species. Using stringent criteria, we delimited a set of 61,771 high-confidence SNPs across 13,543 genes. To validate SNPs, a high-throughput genotyping array was developed for a subset of 5571 predicted SNPs representing as many different gene loci, and was used to genotype over 1000 trees. The estimated true positive rate of the resource was 84.2%, which was comparable with the genotyping success rate obtained for P. abies control SNPs recycled from previous genotyping efforts. We also analyzed SNP abundance across various gene functional categories. Several GO terms and gene families involved in stress response were found over-represented in highly polymorphic genes.

CONCLUSION

The annotated high-confidence SNP catalog developed herein represents a valuable genomic resource, being representative of over 13 K genes distributed across the P. abies genome. This resource should serve a variety of population genomics and breeding applications in Norway spruce.

Expected benefit of genomic selection over forward selection in conifer breeding and deployment

Genomic selection is a proven technology in animal and plant breeding to accelerate genetic gain, but as yet is to be fully realised in forest tree breeding. This paper examines, through stochastic simulation, the potential benefits of genomic selection (GS) over forward selection (FS) in a typical conifer breeding program. Methods of speeding the deployment of selected material were also considered, including top-grafting onto mature seed orchard ortets, using additional replicates of clones in archives for crossing, and embryogenesis and clonal propagation. Genetic gain per generation was found to increase considerably when the size of the training population was larger (800 c.f. 3000 clones), or when the heritability was higher (0.2 c.f. 0.5). The largest genetic gain, of 24% was achieved where large training populations (3000 clones) and high heritability traits (0.5) were combined. The accuracy of genomic breeding values (GEBVs) increased with the increase in the number of clones in the training population, the heritability of the trait and the density of the SNP markers. Calculated accuracies of simulated GEBVs and genetic gain per unit of time suggested that 2000 clones in the training population is the minimum size for effective genomic selection for conifers. Compared with forward selection, genomic selection with 2000 clones in the training population, and a 60K SNP panel, an increase of 1.58 mm per year in diameter-at-breast-height (DBH) and 2.44 kg/m³ per year for wood density can be expected. After one generation (9-years), this would be equivalent to 14.23 mm and 21.97 kg/m³ for DBH and wood density respectively. Deploying clones of the selected individuals always resulted in higher additional genetic gain than deploying progeny/seedlings. Deploying genetic material selected from genomic selection with top-grafting for early coning appeared to be the best option. Application of genomic selection to conifer breeding programs, combined with deployment tools such as top-grafting and embryogenesis are powerful tools to speed the delivery of genetic gain to the forest.

Quantitative Genetics and Genomics Converge to Accelerate Forest Tree Breeding

Forest tree breeding has been successful at delivering genetically improved material for multiple traits based on recurrent cycles of selection, mating, and testing. However, long breeding cycles, late flowering, variable juvenile-mature correlations, emerging pests and diseases, climate, and market changes, all pose formidable challenges. Genetic dissection approaches such as quantitative trait mapping and association genetics have been fruitless to effectively drive operational marker-assisted selection (MAS) in forest trees, largely because of the complex multifactorial inheritance of most, if not all traits of interest. The convergence of high-throughput genomics and quantitative genetics has established two new paradigms that are changing contemporary tree breeding dogmas. Genomic selection (GS) uses large number of genome-wide markers to predict complex phenotypes. It has the potential to accelerate breeding cycles, increase selection intensity and improve the accuracy of breeding values. Realized genomic relationships matrices, on the other hand, provide innovations in genetic parameters' estimation and breeding approaches by tracking the variation arising from random Mendelian segregation in pedigrees. In light of a recent flow of promising experimental results, here we briefly review the main concepts, analytical tools and remaining challenges that currently underlie the application of genomics data to tree breeding. With easy and cost-effective genotyping, we are now at the brink of extensive adoption of GS in tree breeding. Areas for future GS research include optimizing strategies for updating prediction models, adding validated functional genomics data to improve prediction accuracy, and integrating genomic and multi-environment data for forecasting the performance of genetic material in untested sites or under changing climate scenarios. The buildup of phenotypic and genome-wide data across large-scale breeding populations and advances in computational prediction of discrete genomic features should also provide opportunities to enhance the application of genomics to tree breeding.

Prospects and Challenges of Applied Genomic Selection-A New Paradigm in Breeding for Grain Yield in Bread Wheat

Genomic selection (GS) has been promising for increasing genetic gains in several species. Therefore, we evaluated the potential integration of GS for grain yield (GY) in bread wheat ( L.) in CIMMYT's elite yield trial nurseries. We observed that the genomic prediction accuracies within nurseries (0.44 and 0.35) were substantially higher than across-nursery accuracies (0.15 and 0.05) for GY evaluated in the bed and flat planting systems, respectively. The accuracies from using only a subset of 251 genotyping-by-sequencing markers were comparable to the accuracies using all 2038 markers. We also used the item-based collaborative filtering approach for incorporating other related traits in predicting GY and observed that it outperformed genomic predictions across nurseries, but was less predictive when trait correlations with GY were low. Furthermore, we compared GS and phenotypic selections (PS) and observed that at a selection intensity of 0.5, GS could select a maximum of 70.9 and 61.5% of the top lines and discard 71.5 and 60.5% of the poor lines selected or discarded by PS within and across nurseries, respectively. Comparisons of GS and pedigree-based predictions revealed that the advantage of GS over the pedigree was moderate in populations without full-sibs. However, GS was less advantageous for within-family selections in elite families with few full-sibs and minimal Mendelian sampling variance. Overall, our results demonstrate the importance of applying GS for GY at the appropriate stage of the breeding cycle, and we speculate that gains can be maximized if it is implemented in early-generation within-family selections.

Potential of Genome-Wide Studies in Unrelated Plus Trees of a Coniferous Species, Cryptomeria japonica (Japanese Cedar)

A genome-wide association study (GWAS) was conducted on more than 30,000 single nucleotide polymorphisms (SNPs) in unrelated first-generation plus tree genotypes from three populations of Japanese cedar Cryptomeria japonica D. Don with genomic prediction for traits of growth, wood properties and male fecundity. Among the assessed populations, genetic characteristics including the extent of linkage disequilibrium (LD) and genetic structure differed and these differences are considered to be due to differences in genetic background. Through population-independent GWAS, several significant SNPs found close to the regions associated with each of these traits and shared in common across the populations were identified. The accuracies of genomic predictions were dependent on the traits and populations and reflected the genetic architecture of traits and genetic characteristics. Prediction accuracies using SNPs selected based on GWAS results were similar to those using all SNPs for several combinations of traits and populations. We discussed the application of genome-wide studies for C. japonica improvement.

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.

Genomic relationships reveal significant dominance effects for growth in hybrid Eucalyptus

Non-additive genetic effects can be effectively exploited in control-pollinated families with the availability of genome-wide markers. We used 41,304 SNP markers and compared pedigree vs. marker-based genetic models by analysing height, diameter, basic density and pulp yield for Eucalyptus urophylla × E.grandis control-pollinated families represented by 949 informative individuals. We evaluated models accounting for additive, dominance, and first-order epistatic interactions (additive by additive, dominance by dominance, and additive by dominance). We showed that the models can capture a large proportion of the genetic variance from dominance and epistasis for growth traits as those components are typically not independent. We also showed that we could partition genetic variances more precisely when using relationship matrices derived from markers compared to using only pedigree information. In addition, phenotypic prediction accuracies were only slightly increased by including dominance effects for growth traits since estimates of non-additive variances yielded rather high standard errors. This novel result improves our current understanding of the architecture of quantitative traits and recommends accounting for dominance variance when developing genomic selection strategies in hybrid Eucalyptus.