Mapping the genetic architecture of forest tree traits is important in order to understand the evolutionary forces that have shaped these traits and to facilitate the development of genomic-based breeding strategies. We examined the number, size, and distribution of allelic effects influencing eight types of traits using 30 published mapping studies (linkage and association mapping) in forest trees. The sizes of allelic effects, measured as the phenotypic variance explained, generally showed a severely right-skewed distribution. We estimated the numbers of underlying causal effects (n qtl ) for different trait categories by improving a method previously developed by Otto and Jones (Genetics 156:2093(Genetics 156: -2107(Genetics 156: , 2000. Estimates of n qtl based on association mapping studies were generally higher (median at 643) than those based on linkage mapping (median at 33). Comparisons with simulated linkage and association mapping data suggested that the lower n qtl estimates for the linkage mapping studies could partly be explained by fewer causal loci segregating within the full-sib family populations normally used, but also by the cosegregation of causal loci due to limited recombination. Disease resistance estimates based on linkage mapping studies had the lowest median of four underlying effects, while growth traits based on association mapping had about 580 effects. Theoretically, the capture of 50% of the genetic variation would thus require a population size of about 200 for disease resistance in linkage mapping, while growth traits in association mapping would require about 25,000. The adequacy and reliability of the improved method was successfully verified by applying it to the simulated data.

To investigate the potential of association genetics for willow breeding, Salix viminalis germplasm was assembled from UK and Swedish collections (comprising accessions from several European countries) and new samples collected from nature. A subset of the germplasm was planted at two sites (UK and Sweden), genotyped using 38 SSR markers and assessed for phenological and biomass traits. Population structure, genetic differentiation (F ST ) and quantitative trait differentiation (Q ST ) were investigated. The extent and patterns of trait adaptation were assessed by comparing F ST and Q ST parameters. Of the 505 genotyped diploid accessions, 27 % were not unique. Genetic diversity was high: 471 alleles was amplified; the mean number of alleles per locus was 13.46, mean observed heterozygosity was 0.55 and mean expected heterozygosity was 0.62. Bayesian clustering identified four subpopulations which generally corresponded to Western Russia, Western Europe, Eastern Europe and Sweden. All pairwise F ST values were highly significant (p<0.001) with the greatest genetic differentiation detected between the Western Russian and the Western European subpopulations (F ST = 0.12), and the smallest between the Swedish and Eastern European populations (F ST = 0.04). The Swedish population also had the highest number of identical accessions, supporting the view that S. viminalis was introduced into this country and has been heavily influenced by humans. Q ST values were high for growth cessation and leaf senescence, and to some extent stem diameter, but low for bud burst time and shoot number. Overall negative clines between longitudinal coordinates and leaf senescence, bud burst and stem diameter were also found.

Quantitative trait loci (QTL) mapping of wood properties in conifer species has focused on single time point measurements or on trait means based on heterogeneous wood samples (e.g., increment cores), thus ignoring systematic within-tree trends. In this study, functional QTL mapping was performed for a set of important wood properties in increment cores from a 17-yr-old Scots pine (Pinus sylvestris L.) full-sib family with the aim of detecting wood trait QTL for general intercepts (means) and for linear slopes by increasing cambial age. Two multi-locus functional QTL analysis approaches were proposed and their performances were compared on trait datasets comprising 2 to 9 time points, 91 to 455 individual tree measurements and genotype datasets of amplified length polymorphisms (AFLP), and single nucleotide polymorphism (SNP) markers. The first method was a multilevel LASSO analysis whereby trend parameter estimation and QTL mapping were conducted consecutively; the second method was our Bayesian linear mixed model whereby trends and underlying genetic effects were estimated simultaneously. We also compared several different hypothesis testing methods under either the LASSO or the Bayesian framework to perform QTL inference. In total, five and four significant QTL were observed for the intercepts and slopes, respectively, across wood traits such as earlywood percentage, wood density, radial fiberwidth, and spiral grain angle. Four of these QTL were represented by candidate gene SNPs, thus providing promising targets for future research in QTL mapping and molecular function. Bayesian and LASSO methods both detected similar sets of QTL given datasets that comprised large numbers of individuals.

Background Salix spp. are high-productivity crops potentially used for lignocellulosic biofuels such as bioethanol. In general, pretreatment is needed to facilitate the enzymatic depolymerization process. Biomass resistance to degradation, i.e., biomass recalcitrance , is a trait which can be assessed by measuring the sugar released after combined pretreatment and enzymatic hydrolysis. We have examined genetic parameters of enzymatic sugar release and other traits related to biorefinery use in a population of 286 natural Salix viminalis clones. Furthermore, we have evaluated phenotypic and genetic correlations between these traits and performed a genomewide association mapping analysis using a set of 19,411 markers. Results Sugar release (glucose and xylose) after pretreatment and enzymatic saccharification proved highly variable with large genetic and phenotypic variations, and chip heritability estimates ( h 2 ) of 0.23–0.29. Lignin syringyl/guaiacyl (S/G) ratio and wood density were the most heritable traits ( h 2 = 0.42 and 0.59, respectively). Sugar release traits were positively correlated, phenotypically and genetically, with biomass yield and lignin S/G ratio. Association mapping revealed seven marker–trait associations below a suggestive significance threshold, including one marker associated with glucose release. Conclusions We identified lignin S/G ratio and shoot diameter as heritable traits that could be relatively easily evaluated by breeders, making them suitable proxy traits for developing low-recalcitrance varieties. One marker below the suggestive threshold for marker associations was identified for sugar release, meriting further investigation while also highlighting the difficulties in employing genomewide association mapping for complex traits. Electronic supplementary material The online version of this article (10.1186/s13068-019-1479-7) contains supplementary material, which is available to authorized users.

Willow species (Salix) are important as short‐rotation biomass crops for bioenergy, which creates a demand for faster genetic improvement and breeding through deployment of molecular marker‐assisted selection (MAS). To find markers associated with important adaptive traits, such as growth and phenology, for use in MAS, we genetically dissected the trait variation of a Salix viminalis (L.) population of 323 accessions. The accessions were sampled throughout northern Europe and were established at two field sites in Pustnäs, Sweden, and at Woburn, UK, offering the opportunity to assess the impact of genotype‐by‐environment interactions (G × E) on trait–marker associations. Field measurements were recorded for growth and phenology traits. The accessions were genotyped using 1536 SNP markers developed from phenology candidate genes and from genes previously observed to be differentially expressed in contrasting environments. Association mapping between 1233 of these SNPs and the measured traits was performed taking into account population structure and threshold selection bias. At a false discovery rate (FDR) of 0.2, 29 SNPs were associated with bud burst, leaf senescence, number of shoots or shoot diameter. The percentage of accession variation (Rnormaladj2) explained by these associations ranged from 0.3% to 4.4%, suggesting that the studied traits are controlled by many loci of limited individual impact. Despite this, a SNP in the EARLY FLOWERING 3 gene was repeatedly associated (FDR < 0.2) with bud burst. The rare homozygous genotype exhibited 0.4–1.0 lower bud burst scores than the other genotype classes on a five‐grade scale. Consequently, this marker could be promising for use in MAS and the gene deserves further study. Otherwise, associations were less consistent across sites, likely due to their small Rnormaladj2 estimates and to considerable G × E interactions indicated by multivariate association analyses and modest trait accession correlations across sites (0.32–0.61).

The oxidation of melatonin by the mammalian myeloperoxidase (MPO) provides protection against the damaging effects of reactive oxygen species. Indole derivatives, such as melatonin and serotonin, are also substrates of the plant horseradish peroxidase (HRP), but this enzyme exhibits remarkable differences from MPO in the specificity and reaction rates for these compounds. A structural understanding of the determinants of the reactivity of these enzymes to indole derivatives would greatly aid their exploitation for biosynthetic and drug design applications. Consequently, after validation of the docking procedure, we performed computational docking of melatonin and serotonin to structural models of the ferric and compound I and II (co I and co II, respectively) states of HRP and MPO. The substrates dock at the heme edge on the distal side, but with different orientations in the two proteins. The distal cavity is larger in MPO than in HRP; however, in MPO, the substrates make closer contacts with the heme involving ring stacking, whereas in HRP, no ring stacking is observed. The observed differences in substrate binding may contribute to the higher reaction rates and lower substrate specificity of MPO relative to those of HRP. The docking results, along with the previously measured heme-protein reduction potentials, suggest that the differentially lowered reaction rates of co II of HRP and MPO with respect to those of co I could stem from as yet undetermined conformational or electrostatic differences between the co I and co II states of MPO, which are absent in HRP.

BackgroundBiased sex ratios are common among dioecious plant species despite the theoretical prediction of selective advantage of even sex ratios. Albeit the high prevalence of deviations from even sex ratios, the genetic causes to sex biases are rarely known outside of a few model species. Here we present a mechanism underlying the female biased sex ratio in the dioecious willow species Salix viminalis.ResultsWe compared the segregation pattern of genome-wide single nucleotide polymorphism markers in two contrasting bi-parental pedigree populations, the S3 with even sex ratio and the S5 with a female biased sex ratio. With the segregation analysis and comparison between the two populations, we were able to demonstrate that sex determination and sex ratio distortion are controlled by different genetic mechanisms. We furthermore located the sex ratio distorter locus to a Z/W-gametologous region on chromosome 15, which was in close linkage with the sex determination locus. Interestingly, all males in the population with biased sex ratio have in this sex ratio distorter locus the same genotype, meaning that males with the Z1/Z3-genotype were missing from the population, thereby creating the 2:1 female biased sex ratio.ConclusionsWe attribute the absence of Z1/Z3 males to an allelic incompatibility between maternally and paternally inherited alleles in this sex ratio distorter locus. Due to the tight linkage with the sex determination locus only male individuals are purged from the population at an early age, presumably before or during seed development. We showed that such allelic incompatibility could be stably maintained over evolutionary times through a system of overdominant or pseudooverdominant alleles. Thus, it is possible that the same mechanism generates the female biased sex ratio in natural willow populations.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-017-3634-5) contains supplementary material, which is available to authorized users.

The short rotation biomass crop willow ( Salix genera) has been of interest for bioenergy but recently also for biofuel production. For a faster development of new varieties molecular markers could be used as selection tool in an early stage of the breeding cycle. To identify markers associated with growth traits, genome-wide association mapping was conducted using a population of 291 Salix viminalis accessions collected across Europe and Russia and a large set of genotyping-by-sequencing markers. The accessions were vegetatively propagated and planted in replicated field experiments, one in Southern Sweden and one in Central Sweden. Phenology data, including bud burst and leaf senescence, as well as different growth traits were collected and measured repeatedly between 2010 and 2017 at both field environments. A value of the plasticity for each accession was calculated for all traits that were measured the same year in both environments as the normalized accession value in one environment subtracted by the corresponding value in the other environment. Broad-sense accession heritabilities and narrow-sense chip heritabilities ranged from 0.68 to 0.95 and 0.45 to 0.99, respectively for phenology traits and from 0.56 to 0.85 and 0.24 to 0.97 for growth traits indicating a considerable genetic component for most traits. Population structure and kinship between accessions were taken into account in the association analyses. In total, 39 marker-trait associations were found where four were specifically connected to plasticity and interestingly one particular marker was associated to several different plasticity growth traits. Otherwise association consistency was poor, possibly due to accession by environment interactions which were demonstrated by the low structure adjusted accession correlations across environments (ranging from 0.40 to 0.58). However, one marker association with biomass fresh weight was repeatedly observed in the same environment over two harvest years. For some traits where several associations were found, the markers jointly explained over 20% of the accession variation. The result from this study using a population of unrelated accessions has given useful information about marker-trait associations especially highlighting marker-plasticity associations and genotype-by-environment interactions as important factors to take account of in future strategies of Salix breeding.