RESUMEN
Forests face many threats. While traditional breeding may be too slow to deliver well-adapted trees, genomic selection (GS) can accelerate the process. We describe a comprehensive study of GS from proof of concept to operational application in western redcedar (WRC, Thuja plicata). Using genomic data, we developed models on a training population (TrP) of trees to predict breeding values (BVs) in a target seedling population (TaP) for growth, heartwood chemistry, and foliar chemistry traits. We used cross-validation to assess prediction accuracy (PACC) in the TrP; we also validated models for early-expressed foliar traits in the TaP. Prediction accuracy was high across generations, environments, and ages. PACC was not reduced to zero among unrelated individuals in TrP and was only slightly reduced in the TaP, confirming strong linkage disequilibrium and the ability of the model to generate accurate predictions across breeding generations. Genomic BV predictions were correlated with those from pedigree but displayed a wider range of within-family variation due to the ability of GS to capture the Mendelian sampling term. Using predicted TaP BVs in multi-trait selection, we functionally implemented and integrated GS into an operational tree-breeding program.
RESUMEN
Western redcedar (WRC; Thuja plicata) is a conifer of the Pacific Northwest of North America prized for its durable and rot-resistant wood. WRC has naturally low outcrossing rates and readily self-fertilizes in nature. Challenges faced in WRC breeding and propagation involve selecting trees for accelerated growth while also ensuring enhanced heartwood rot resistance and resistance to ungulate browsing, as well as mitigating potential effects of inbreeding depression. Terpenes, a large and diverse class of specialized metabolites, confer both rot and browse resistance in the wood and foliage of WRC, respectively. Using a Bayesian modelling approach, we isolated single nucleotide polymorphism (SNP) markers estimated to be associated with three different foliar terpene traits and four different heartwood terpene traits, as well as two growth traits. We found that all traits were complex, being associated with between 1700 and 3600 SNPs linked with putatively causal loci, with significant polygenic components. Growth traits tended to have a larger polygenic component while terpene traits had larger major gene components; SNPs with small or polygenic effect were spread across the genome, while larger-effect SNPs tended to be localized to specific linkage groups. To determine whether there was inbreeding depression for terpene chemistry or growth traits, we used mixed linear models for a genomic selection training population to estimate the effect of the inbreeding coefficient F on foliar terpenes, heartwood terpenes and several growth and dendrochronological traits. We did not find significant inbreeding depression for any assessed trait. We further assessed inbreeding depression across four generations of complete selfing and found that not only was inbreeding depression not significant but that selection for height growth was the only significant predictor for growth during selfing, suggesting that inbreeding depression due to selfing during operational breeding can be mitigated by increased selection intensity.
RESUMEN
We assembled the 9.8-Gbp genome of western redcedar (WRC; Thuja plicata), an ecologically and economically important conifer species of the Cupressaceae. The genome assembly, derived from a uniquely inbred tree produced through five generations of self-fertilization (selfing), was determined to be 86% complete by BUSCO analysis, one of the most complete genome assemblies for a conifer. Population genomic analysis revealed WRC to be one of the most genetically depauperate wild plant species, with an effective population size of approximately 300 and no significant genetic differentiation across its geographic range. Nucleotide diversity, π, is low for a continuous tree species, with many loci showing zero diversity, and the ratio of π at zero- to fourfold degenerate sites is relatively high (approximately 0.33), suggestive of weak purifying selection. Using an array of genetic lines derived from up to five generations of selfing, we explored the relationship between genetic diversity and mating system. Although overall heterozygosity was found to decline faster than expected during selfing, heterozygosity persisted at many loci, and nearly 100 loci were found to deviate from expectations of genetic drift, suggestive of associative overdominance. Nonreference alleles at such loci often harbor deleterious mutations and are rare in natural populations, implying that balanced polymorphisms are maintained by linkage to dominant beneficial alleles. This may account for how WRC remains responsive to natural and artificial selection, despite low genetic diversity.
Asunto(s)
Tracheophyta , Tracheophyta/genética , Autofecundación/genética , Alelos , Heterocigoto , Polimorfismo Genético , Variación Genética , Selección GenéticaRESUMEN
Western redcedar (WRC) is an ecologically and economically important forest tree species characterized by low genetic diversity with high self-compatibility and high heartwood durability. Using sequence capture genotyping of target genic and non-genic regions, we genotyped 44 parent trees and 1520 offspring trees representing 26 polycross (PX) families collected from three progeny test sites using 45,378 SNPs. Trees were phenotyped for eight traits related to growth, heartwood and foliar chemistry associated with wood durability and deer browse resistance. We used the genomic realized relationship matrix for paternity assignment, maternal pedigree correction, and to estimate genetic parameters. We compared genomics-based (GBLUP) and two pedigree-based (ABLUP: polycross and reconstructed full-sib [FS] pedigrees) models. Models were extended to estimate dominance genetic effects. Pedigree reconstruction revealed significant unequal male contribution and separated the 26 PX families into 438 FS families. Traditional maternal PX pedigree analysis resulted in up to 51% overestimation in genetic gain and 44% in diversity. Genomic analysis resulted in up to 22% improvement in offspring breeding value (BV) theoretical accuracy, 35% increase in expected genetic gain for forward selection, and doubled selection intensity for backward selection. Overall, all traits showed low to moderate heritability (0.09-0.28), moderate genotype by environment interaction (type-B genetic correlation: 0.51-0.80), low to high expected genetic gain (6.01%-55%), and no significant negative genetic correlation reflecting no large trade-offs for multi-trait selection. Only three traits showed a significant dominance effect. GBLUP resulted in smaller but more accurate heritability estimates for five traits, but larger estimates for the wood traits. Comparison between all, genic-coding, genic-non-coding and intergenic SNPs showed little difference in genetic estimates. In summary, we show that GBLUP overcomes the PX limitations, successfully captures expected historical and hidden relatedness as well as linkage disequilibrium (LD), and results in increased breeding efficiency in WRC.