Implications of accounting for marker-based population structure in the quantitative genetic evaluation of genetic parameters related to growth and wood properties in Norway spruce.

BACKGROUND: Forest geneticists typically use provenances to account for population differences in their improvement schemes; however, the historical records of the imported materials might not be very precise or well-aligned with the genetic clusters derived from advanced molecular techniques. The main objective of this study was to assess the impact of marker-based population structure on genetic parameter estimates related to growth and wood properties and their trade-offs in Norway spruce, by either incorporating it as a fixed effect (model-A) or excluding it entirely from the analysis (model-B). RESULTS: Our results indicate that models incorporating population structure significantly reduce estimates of additive genetic variance, resulting in substantial reduction of narrow-sense heritability. However, these models considerably improve prediction accuracies. This was particularly significant for growth and solid-wood properties, which showed to have the highest population genetic differentiation (QST) among the studied traits. Additionally, although the pattern of correlations remained similar across the models, their magnitude was slightly lower for models that included population structure as a fixed effect. This suggests that selection, consistently performed within populations, might be less affected by unfavourable genetic correlations compared to mass selection conducted without pedigree restrictions. CONCLUSION: We conclude that the results of models properly accounting for population structure are more accurate and less biased compared to those neglecting this effect. This might have practical implications for breeders and forest managers where, decisions based on imprecise selections can pose a high risk to economic efficiency.

Divergent selection predating the Last Glacial Maximum mainly acted on macro-phenotypes in Norway spruce.

The current distribution and population structure of many species were, to a large extent, shaped by cycles of isolation in glacial refugia and subsequent population expansions. Isolation in and postglacial expansion through heterogeneous environments led to either neutral or adaptive divergence. Norway spruce is no exception, and its current distribution is the consequence of a constant interplay between evolutionary and demographic processes. We investigated population differentiation and adaptation of Norway spruce for juvenile growth, diameter of the stem, wood density, and tracheid traits at breast height. Data from 4461 phenotyped and genotyped Norway spruce from 396 half-sib families in two progeny tests were used to test for divergent selection in the framework of Q ST vs. F ST. We show that the macroscopic resultant trait (stem diameter), unlike its microscopic components (tracheid dimensions) and juvenile growth, was under divergent selection that predated the Last Glacial Maximum. Altogether, the current variation in these phenotypic traits in Norway spruce is better explained by local adaptation to ancestral environments than to current ones, where populations were partly preadapted, mainly through growth-related traits.

Genetic control of tracheid properties in Norway spruce wood.

Through the use of genome-wide association studies (GWAS) mapping it is possible to establish the genetic basis of phenotypic trait variation. Our GWAS study presents the first such effort in Norway spruce (Picea abies (L). Karst.) for the traits related to wood tracheid characteristics. The study employed an exome capture genotyping approach that generated 178 101 Single Nucleotide Polymorphisms (SNPs) from 40 018 probes within a population of 517 Norway spruce mother trees. We applied a least absolute shrinkage and selection operator (LASSO) based association mapping method using a functional multi-locus mapping approach, with a stability selection probability method as the hypothesis testing approach to determine significant Quantitative Trait Loci (QTLs). The analysis has provided 30 significant associations, the majority of which show specific expression in wood-forming tissues or high ubiquitous expression, potentially controlling tracheids dimensions, their cell wall thickness and microfibril angle. Among the most promising candidates based on our results and prior information for other species are: Picea abies BIG GRAIN 2 (PabBG2) with a predicted function in auxin transport and sensitivity, and MA_373300g0010 encoding a protein similar to wall-associated receptor kinases, which were both associated with cell wall thickness. The results demonstrate feasibility of GWAS to identify novel candidate genes controlling industrially-relevant tracheid traits in Norway spruce.

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.

Genome-wide association study identified novel candidate loci affecting wood formation in Norway spruce.

Norway spruce is a boreal forest tree species of significant ecological and economic importance. Hence there is a strong imperative to dissect the genetics underlying important wood quality traits in the species. We performed a functional genome-wide association study (GWAS) of 17 wood traits in Norway spruce using 178 101 single nucleotide polymorphisms (SNPs) generated from exome genotyping of 517 mother trees. The wood traits were defined using functional modelling of wood properties across annual growth rings. We applied a Least Absolute Shrinkage and Selection Operator (LASSO-based) association mapping method using a functional multilocus mapping approach that utilizes latent traits, with a stability selection probability method as the hypothesis testing approach to determine a significant quantitative trait locus. The analysis provided 52 significant SNPs from 39 candidate genes, including genes previously implicated in wood formation and tree growth in spruce and other species. Our study represents a multilocus GWAS for complex wood traits in Norway spruce. The results advance our understanding of the genetics influencing wood traits and identifies candidate genes for future functional studies.

Downregulating aspen xylan biosynthetic GT43 genes in developing wood stimulates growth via reprograming of the transcriptome.

Xylan is one of the main compounds determining wood properties in hardwood species. The xylan backbone is thought to be synthesized by a synthase complex comprising two members of the GT43 family. We downregulated all GT43 genes in hybrid aspen (Populus tremula × tremuloides) to understand their involvement in xylan biosynthesis. All three clades of the GT43 family were targeted for downregulation using RNA interference individually or in different combinations, either constitutively or specifically in developing wood. Simultaneous downregulation in developing wood of the B (IRX9) and C (IRX14) clades resulted in reduced xylan Xyl content relative to reducing end sequence, supporting their role in xylan backbone biosynthesis. This was accompanied by a higher lignocellulose saccharification efficiency. Unexpectedly, GT43 suppression in developing wood led to an overall growth stimulation, xylem cell wall thinning and a shift in cellulose orientation. Transcriptome profiling of these transgenic lines indicated that cell cycling was stimulated and secondary wall biosynthesis was repressed. We suggest that the reduced xylan elongation is sensed by the cell wall integrity surveying mechanism in developing wood. Our results show that wood-specific suppression of xylan-biosynthetic GT43 genes activates signaling responses, leading to increased growth and improved lignocellulose saccharification.

A collection of genetically engineered Populus trees reveals wood biomass traits that predict glucose yield from enzymatic hydrolysis.

Wood represents a promising source of sugars to produce bio-based renewables, including biofuels. However, breaking down lignocellulose requires costly pretreatments because lignocellulose is recalcitrant to enzymatic saccharification. Increasing saccharification potential would greatly contribute to make wood a competitive alternative to petroleum, but this requires improving wood properties. To identify wood biomass traits associated with saccharification, we analyzed a total of 65 traits related to wood chemistry, anatomy and structure, biomass production and saccharification in 40 genetically engineered Populus tree lines. These lines exhibited broad variation in quantitative traits, allowing for multivariate analyses and mathematical modeling. Modeling revealed that seven wood biomass traits associated in a predictive manner with saccharification of glucose after pretreatment. Four of these seven traits were also negatively associated with biomass production, suggesting a trade-off between saccharification potential and total biomass, which has previously been observed to offset the overall sugar yield from whole trees. We therefore estimated the "total-wood glucose yield" (TWG) from whole trees and found 22 biomass traits predictive of TWG after pretreatment. Both saccharification and TWG were associated with low abundant, often overlooked matrix polysaccharides such as arabinose and rhamnose which possibly represent new markers for improved Populus feedstocks.

Wood properties of Populus and Betula in long-term exposure to elevated CO2 and O3.

We studied the interactive effects of elevated concentrations of CO2 and O3 on radial growth and wood properties of four trembling aspen (Populus tremuloides Michx.) clones and paper birch (Betula papyrifera Marsh.) saplings. The material for the study was collected from the Aspen FACE (free-air CO2 enrichment) experiment in Rhinelander (WI, USA). Trees had been exposed to four treatments [control, elevated CO2 (560 ppm), elevated O3 (1.5 times ambient) and combined CO2 + O3 ] during growing seasons 1998-2008. Most treatment responses were observed in the early phase of experiment. Our results show that the CO2- and O3-exposed aspen trees displayed a differential balance between efficiency and safety of water transport. Under elevated CO2, radial growth was enhanced and the trees had fewer but hydraulically more efficient larger diameter vessels. In contrast, elevated O3 decreased radial growth and the diameters of vessels and fibres. Clone-specific decrease in wood density and cell wall thickness was observed under elevated CO2 . In birch, the treatments had no major impacts on wood anatomy or wood density. Our study indicates that short-term impact studies conducted with young seedlings may not give a realistic view of long-term ecosystem responses.