Association genetics of wood physical traits in the conifer white spruce and relationships with gene expression

Genotyping SNPs in lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2) exhibits association with wood density in teak (Tectona grandis L.f.)

BACKGROUND

Teak (Tectona grandis L.f.), an important source of tropical timber with immense economic value, is a highly outcrossing forest tree species. 150 unrelated accessions of teak (Tectona grandis L.f.) plus trees assembled as clones at National Teak Germplasm Bank, Chandrapur, Maharashtra, India was investigated for association mapping of candidate lignin biosynthesis gene (CAD1) and transcription factors (MYB1 and MYB2).

METHODS AND RESULTS

The CAD1, MYB1 and MYB2 were amplified using specifically designed primers. The amplified sequences were then sequenced and genotyped for 112 SNPs/11 indels. We evaluated the association between SNPs and wood density in teak accessions using GLM and MLM statistical models, with Bonferroni correction applied. The teak accessions recorded an average wood density of 416.69 kg.m-3 (CV 4.97%) and comprised of three loosely structured admixed sub-populations (K = 3), containing 72.05% genetic variation within sub-populations with low intragenic LD (0-21% SNP pairs) at P < 0.05 and high LD decay (33-934 bp) at R2 = 0.1. GLM and MLM models discounting systematic biases (Q and K matrices) to avoid false discovery revealed five loci at rare variants (MAF 0.003) and three loci at common variants (MAF 0.05) to be significantly (P < 0.05) associated with the wood density. However, the stringent Bonferroni correction (4.06-7.04 × 10-4) yielded only a single associated locus (B1485C/A) from exon of MYB1 transcription factor, contributing to about 10.35% phenotypic variation in wood density trait.

CONCLUSION

Scored SNP locus (B1485C/A) can be developed as a molecular probe for selection of improved planting stock with proven wood density trait for a large-scale teak plantation.

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.

Genetic architecture behind developmental and seasonal control of tree growth and wood properties in Norway spruce

Genetic control of tree growth and wood formation varies depending on the age of the tree and the time of the year. Single-locus, multi-locus, and multi-trait genome-wide association studies (GWAS) were conducted on 34 growth and wood property traits in 1,303 Norway spruce individuals using exome capture to cover ~130K single-nucleotide polymorphisms (SNPs). GWAS identified associations to the different wood traits in a total of 85 gene models, and several of these were validated in a progenitor population. A multi-locus GWAS model identified more SNPs associated with the studied traits than single-locus or multivariate models. Changes in tree age and annual season influenced the genetic architecture of growth and wood properties in unique ways, manifested by non-overlapping SNP loci. In addition to completely novel candidate genes, SNPs were located in genes previously associated with wood formation, such as cellulose synthases and a NAC transcription factor, but that have not been earlier linked to seasonal or age-dependent regulation of wood properties. Interestingly, SNPs associated with the width of the year rings were identified in homologs of Arabidopsis thaliana BARELY ANY MERISTEM 1 and rice BIG GRAIN 1, which have been previously shown to control cell division and biomass production. The results provide tools for future Norway spruce breeding and functional studies.

Transcriptome and association mapping revealed functional genes respond to drought stress in Populus

Drought frequency and severity are exacerbated by global climate change, which could compromise forest ecosystems. However, there have been minimal efforts to systematically investigate the genetic basis of the response to drought stress in perennial trees. Here, we implemented a systems genetics approach that combines co-expression analysis, association genetics, and expression quantitative trait nucleotide (eQTN) mapping to construct an allelic genetic regulatory network comprising four key regulators (PtoeIF-2B, PtoABF3, PtoPSB33, and PtoLHCA4) under drought stress conditions. Furthermore, Hap_01PtoeIF-2B, a superior haplotype associated with the net photosynthesis, was revealed through allelic frequency and haplotype analysis. In total, 75 candidate genes related to drought stress were identified through transcriptome analyses of five Populus cultivars (P. tremula × P. alba, P. nigra, P. simonii, P. trichocarpa, and P. tomentosa). Through association mapping, we detected 92 unique SNPs from 38 genes and 104 epistatic gene pairs that were associated with six drought-related traits by association mapping. eQTN mapping unravels drought stress-related gene loci that were significantly associated with the expression levels of candidate genes for drought stress. In summary, we have developed an integrated strategy for dissecting a complex genetic network, which facilitates an integrated population genomics approach that can assess the effects of environmental threats.

Identifying and testing marker‐trait associations for growth and phenology in three pine species: implications for genomic prediction

In tree species, genomic prediction offers the potential to forecast mature trait values in early growth stages, if robust marker–trait associations can be identified. Here we apply a novel multispecies approach using genotypes from a new genotyping array, based on 20,795 single nucleotide polymorphisms (SNPs) from three closely related pine species (Pinus sylvestris, Pinus uncinata and Pinus mugo), to test for associations with growth and phenology data from a common garden study. Predictive models constructed using significantly associated SNPs were then tested and applied to an independent multisite field trial of P. sylvestris and the capability to predict trait values was evaluated. One hundred and eighteen SNPs showed significant associations with the traits in the pine species. Common SNPs (MAF > 0.05) associated with bud set were only found in genes putatively involved in growth and development, whereas those associated with growth and budburst were also located in genes putatively involved in response to environment and, to a lesser extent, reproduction. At one of the two independent sites, the model we developed produced highly significant correlations between predicted values and observed height data (YA, height 2020: r = 0.376, p < 0.001). Predicted values estimated with our budburst model were weakly but positively correlated with duration of budburst at one of the sites (GS, 2015: r = 0.204, p = 0.034; 2018: r = 0.205, p = 0.034–0.037) and negatively associated with budburst timing at the other (YA: r = −0.202, p = 0.046). Genomic prediction resulted in the selection of sets of trees whose mean height was taller than the average for each site. Our results provide tentative support for the capability of prediction models to forecast trait values in trees, while highlighting the need for caution in applying them to trees grown in different environments.

Genetic basis of growth, spring phenology, and susceptibility to biotic stressors in maritime pine

Forest ecosystems are increasingly challenged by extreme events, for example, drought, storms, pest attacks, and fungal pathogen outbreaks, causing severe ecological and economic losses. Understanding the genetic basis of adaptive traits in tree species is of key importance to preserve forest ecosystems, as genetic variation in a trait (i.e., heritability) determines its potential for human-mediated or evolutionary change. Maritime pine (Pinus pinaster Aiton), a conifer widely distributed in southwestern Europe and northwestern Africa, grows under contrasted environmental conditions promoting local adaptation. Genetic variation at adaptive phenotypes, including height, spring phenology, and susceptibility to two fungal pathogens (Diplodia sapinea and Armillaria ostoyae) and an insect pest (Thaumetopoea pityocampa), was assessed in a range-wide clonal common garden of maritime pine. Broad-sense heritability was significant for height (0.219), spring phenology (0.165-0.310), and pathogen susceptibility (necrosis length caused by D. sapinea, 0.152; and by A. ostoyae, 0.021, measured on inoculated, excised branches under controlled conditions), but not for pine processionary moth incidence in the common garden. The correlations of trait variation among populations revealed contrasting trends for pathogen susceptibility to D. sapinea and A. ostoyae with respect to height. Taller trees showed longer necrosis length caused by D. sapinea while shorter trees were more affected by A. ostoyae. Moreover, maritime pine populations from areas with high summer temperatures and frequent droughts were less susceptible to D. sapinea but more susceptible to A. ostoyae. Finally, an association study using 4227 genome-wide SNPs revealed several loci significantly associated with each trait (range of 3-26), including a possibly disease-induced translation initiation factor, eIF-5, associated with needle discoloration caused by D. sapinea. This study provides important insights to develop genetic conservation and breeding strategies integrating species responses to biotic stressors.

Population structure and the influence of microenvironment and genetic similarity on individual growth at Alaskan white spruce treelines

Knowledge on the adaptation of trees to rapid environmental changes is essential to preserve forests and their ecosystem services under climate change. Treeline populations are particularly suitable for studying adaptation processes in trees, as environmental stress together with reduced gene flow can enhance local adaptation. We investigated white spruce (Picea glauca) populations in Alaska on one moisture-limited and two cold-limited treeline sites with a paired plot design of one forest and one treeline population each, resulting in six plots. Additionally, one forest plot in the middle of the distribution range complements the study design. We combined spatial, climatic and dendrochronological data with neutral genetic marker of 2203 trees to investigate population genetic structure and drivers of tree growth. We used several individual-based approaches including random slope mixed-effects models to test the influence of genetic similarity and microenvironment on growth performance. A high degree of genetic diversity was found within each of the seven plots associated with high rates of gene flow. We discovered a low genetic differentiation between the three sites which was better explained by geographic distances than by environmental differences, indicating genetic drift as the main driver of population differentiation. Our findings indicated that microenvironmental features had an overall larger influence on growth performances than genetic similarity among individuals. The effects of climate on growth differed between sites but were smaller than the effect of tree size. Overall, our results suggest that the high genetic diversity of white spruce may result in a wider range of phenotypes which enhances the efficiency of selection when the species is facing rapid climatic changes. In addition, the large intra-individual variability in growth responses may indicate the high phenotypic plasticity of white spruce which can buffer short-term environmental changes and, thus, allow enduring the present changing climate conditions.

Selective Sweeps and Polygenic Adaptation Drive Local Adaptation along Moisture and Temperature Gradients in Natural Populations of Coast Redwood and Giant Sequoia

Dissecting the genomic basis of local adaptation is a major goal in evolutionary biology and conservation science. Rapid changes in the climate pose significant challenges to the survival of natural populations, and the genomic basis of long-generation plant species is still poorly understood. Here, we investigated genome-wide climate adaptation in giant sequoia and coast redwood, two iconic and ecologically important tree species. We used a combination of univariate and multivariate genotype-environment association methods and a selective sweep analysis using non-overlapping sliding windows. We identified genomic regions of potential adaptive importance, showing strong associations to moisture variables and mean annual temperature. Our results found a complex architecture of climate adaptation in the species, with genomic regions showing signatures of selective sweeps, polygenic adaptation, or a combination of both, suggesting recent or ongoing climate adaptation along moisture and temperature gradients in giant sequoia and coast redwood. The results of this study provide a first step toward identifying genomic regions of adaptive significance in the species and will provide information to guide management and conservation strategies that seek to maximize adaptive potential in the face of climate change.

Dissection of Allelic Variation Underlying Floral and Fruit Traits in Flare Tree Peony (Paeonia rockii) Using Association Mapping

Allelic variation in floral quantitative traits, including the elements of flowers and fruits, is caused by extremely complex regulatory processes. In the genetic improvement of flare tree peony (Paeonia rockii), a unique ornamental and edible oil woody species in the genus Paeonia, a better understanding of the genetic composition of these complex traits related to flowers and fruits is needed. Therefore, we investigated the genetic diversity and population structure of 160 P. rockii accessions and conducted single-marker association analysis for 19 quantitative flower and fruit traits using 81 EST-SSR markers. The results showed that the population had a high phenotypic diversity (coefficients of variation, 11.87-110.64%) and a high level of genetic diversity (mean number of alleles, N _A = 6.09). These accessions were divided into three subgroups by STRUCTURE analysis and a neighbor-joining tree. Furthermore, we also found a low level of linkage disequilibrium between these EST-SSRs and, by single-marker association analysis, identified 134 significant associations, including four flower traits with 11 EST-SSRs and 10 fruit traits with 32 EST-SSRs. Finally, based on the sequence alignment of the associated markers, P280, PS2, PS12, PS27, PS118, PS131, and PS145 may be considered potential loci to increase the yield of flare tree peony. These results laid the foundation for further analysis of the genetic structure of some key traits in P. rockii and had an obvious potential application value in marker-assisted selection breeding.

Connecting tree-ring phenotypes, genetic associations and transcriptomics to decipher the genomic architecture of drought adaptation in a widespread conifer

As boreal forests face significant threats from climate change, understanding evolutionary trajectories of coniferous species has become fundamental to adapting management and conservation to a drying climate. We examined the genomic architecture underlying adaptive variation related to drought tolerance in 43 populations of a widespread boreal conifer, white spruce (Piceaglauca [Moench] Voss), by combining genotype–environment associations, genotype–phenotype associations, and transcriptomics. Adaptive genetic variation was identified by correlating allele frequencies for 6,153 single nucleotide polymorphisms from 2,606 candidate genes with temperature, precipitation and aridity gradients, and testing for significant associations between genotypes and 11 dendrometric and drought‐related traits (i.e., anatomical, growth response and climate‐sensitivity traits) using a polygenic model. We identified a set of 285 genes significantly associated with a climatic factor or a phenotypic trait, including 110 that were differentially expressed in response to drought under greenhouse‐controlled conditions. The interlinked phenotype–genotype–environment network revealed eight high‐confidence genes involved in white spruce adaptation to drought, of which four were drought‐responsive in the expression analysis. Our findings represent a significant step toward the characterization of the genomic basis of drought tolerance and adaptation to climate in conifers, which is essential to enable the establishment of resilient forests in view of new climate conditions.

see also the Perspective by Lars Opgenoorth and Christian Rellstab

Development of a highly efficient 50K single nucleotide polymorphism genotyping array for the large and complex genome of Norway spruce (Picea abies L. Karst) by whole genome resequencing and its transferability to other spruce species

Norway spruce (Picea abies L. Karst) is one of the most important forest tree species with significant economic and ecological impact in Europe. For decades, genomic and genetic studies on Norway spruce have been challenging due to the large and repetitive genome (19.6 Gb with more than 70% being repetitive). To accelerate genomic studies, including population genetics, genome‐wide association studies (GWAS) and genomic selection (GS), in Norway spruce and related species, we here report on the design and performance of a 50K single nucleotide polymorphism (SNP) genotyping array for Norway spruce. The array is developed based on whole genome resequencing (WGS), making it the first WGS‐based SNP array in any conifer species so far. After identifying SNPs using genome resequencing data from 29 trees collected in northern Europe, we adopted a two‐step approach to design the array. First, we built a 450K screening array and used this to genotype a population of 480 trees sampled from both natural and breeding populations across the Norway spruce distribution range. These samples were then used to select high‐confidence probes that were put on the final 50K array. The SNPs selected are distributed over 45,552 scaffolds from the P. abies version 1.0 genome assembly and target 19,954 unique gene models with an even coverage of the 12 linkage groups in Norway spruce. We show that the array has a 99.5% probe specificity, >98% Mendelian allelic inheritance concordance, an average sample call rate of 96.30% and an SNP call rate of 98.90% in family trios and haploid tissues. We also observed that 23,797 probes (50%) could be identified with high confidence in three other spruce species (white spruce [Picea glauca], black spruce [P. mariana] and Sitka spruce [P. sitchensis]). The high‐quality genotyping array will be a valuable resource for genetic and genomic studies in Norway spruce as well as in other conifer species of the same genus.

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.

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.

Development of EST-SSR markers and association mapping with floral traits in Syringa oblata

BACKGROUND

Lilac (Syringa oblata) is an important woody plant with high ornamental value. However, very limited genetic marker resources are currently available, and little is known about the genetic architecture of important ornamental traits for S. oblata, which is hindering its genetic studies. Therefore, it is of great significance to develop effective molecular markers and understand the genetic architecture of complex floral traits for the genetic research of S. oblata.

RESULTS

In this study, a total of 10,988 SSRs were obtained from 9864 unigene sequences with an average of one SSR per 8.13 kb, of which di-nucleotide repeats were the dominant type (32.86%, 3611). A set of 2042 primer pairs were validated, out of which 932 (45.7%) exhibited successful amplifications, and 248 (12.1%) were polymorphic in eight S. oblata individuals. In addition, 30 polymorphic EST-SSR markers were further used to assess the genetic diversity and the population structure of 192 cultivated S. oblata individuals. Two hundred thirty-four alleles were detected, and the PIC values ranged from 0.23 to 0.88 with an average of 0.51, indicating a high level of genetic diversity within this cultivated population. The analysis of population structure showed two major subgroups in the association population. Finally, 20 significant associations were identified involving 17 markers with nine floral traits using the mixed linear model. Moreover, marker SO104, SO695 and SO790 had significant relationship with more than one trait.

CONCLUSION

The results showed newly developed markers were valuable resource and provided powerful tools for genetic breeding of lilac. Beyond that, our study could serve an efficient foundation for further facilitate genetic improvement of floral traits for lilac.

Genetic interactions among Pto-miR319 family members and their targets influence growth and wood properties in Populus tomentosa

MicroRNAs (miRNAs) play crucial roles in all aspects of plant growth and development, but the genetic interactions of miRNAs and their target genes in woody plants are largely unknown. Here, we integrated association genetics and expression profiling to decipher the allelic variations and interactions of the Pto-MIR319 family of miRNAs and 12 putative Pto-miR319 target genes related to wood formation in 435 unrelated individuals of Populus tomentosa Carrière (Chinese white poplar). Expression pattern analysis showed that among all pairings between expressions of pre-miRNA of Pto-MIR319 members and targets, 70.0% showed negative correlation of expression levels (r = - 0.944 to 0.674, P < 0.01) in eight tissues and organs of poplar, suggesting that Pto-miR319 may participate in the regulatory network of wood formation. Single SNP-based association studies identified 137 significant associations (P < 0.01, Q < 0.1), representing 126 unique SNPs from Pto-MIR319 members and their targets, with 10 tree growth traits, revealing that these genetic factors have common roles related to wood formation. Epistasis analysis uncovered 105 significant SNP-SNP associations (P < 0.01) influencing the 10 traits, demonstrating the close genetic interactions between Pto-MIR319 family members and the 12 Pto-miR319 target genes. Notably, one common SNP, in the precursor region of Pto-MIR319e, affected the stability of Pto-MIR319e's secondary structure by altering the stem-loop structure and minimum free energy, contributing to variations in the expression of Pto-MIR319e and Pto-miR319e target genes. This study enriches the understanding of the functions of miR319 family miRNAs in poplar and exemplifies a feasible approach to exploring the genetic effects underlying miRNA-mRNA interactions related to complex traits in trees.

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.

The Cytoskeleton and Its Role in Determining Cellulose Microfibril Angle in Secondary Cell Walls of Woody Tree Species

Recent advances in our understanding of the molecular control of secondary cell wall (SCW) formation have shed light on molecular mechanisms that underpin domestication traits related to wood formation. One such trait is the cellulose microfibril angle (MFA), an important wood quality determinant that varies along tree developmental phases and in response to gravitational stimulus. The cytoskeleton, mainly composed of microtubules and actin filaments, collectively contribute to plant growth and development by participating in several cellular processes, including cellulose deposition. Studies in Arabidopsis have significantly aided our understanding of the roles of microtubules in xylem cell development during which correct SCW deposition and patterning are essential to provide structural support and allow for water transport. In contrast, studies relating to SCW formation in xylary elements performed in woody trees remain elusive. In combination, the data reviewed here suggest that the cytoskeleton plays important roles in determining the exact sites of cellulose deposition, overall SCW patterning and more specifically, the alignment and orientation of cellulose microfibrils. By relating the reviewed evidence to the process of wood formation, we present a model of microtubule participation in determining MFA in woody trees forming reaction wood (RW).

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.

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.

The Resistance of Scots Pine (Pinus sylvestris L.) Half-sib Families to Heterobasidion annosum

This study observed the genotypic variation among Scots pine (Pinus sylvestris L.) half-sib families’ susceptibility to Heterobasidion annosum. Scots pine susceptibility was tested in 12 half-sib families by inoculating them with four different H. annosum strains. At two, six, and ten months after the inoculations, the susceptibility indicators (incidence rate, pathogen spread, and mortality rate) were compared and the total phenolic compounds (TPC) in the inoculated and control groups determined using the Folin–Ciocalteu method. Among half-sib families, significant differences were found for seedling mortality (range: 1.3%–21.2%); however, the differences in incidence rate (range: 54%–77%) and pathogen spread (range: 24–53 mm) were not significant. The incidence rate among half-sib families correlated positively and significantly (r = 0.72, p < 0.05) with the mortality rate, while the pathogen spread correlated negatively with mortality, although the correlation was not significant (r = −0.29, p > 0.05). The TPC comparison with susceptibility indicators showed that the half-sib families with lower susceptibility were characterized by the ability to increase TPC after inoculation compared to the control group. This tendency was most apparent in stems and roots six and two months after inoculation, respectively. Correlation analyses revealed that higher TPC in stems six months after inoculation determined a lower incidence rate (r = −0.32, p < 0.05), while higher Change in concentration of total phenolic compounds (TPCΔ) indicated a lower pathogen spread (r = −0.60, p < 0.05). The lowest incidence of the pathogen was determined in half-sib families with the highest TPCΔ in the roots two months after inoculation. The lower susceptibility of Scots pine half-sib families is based on a combination of enhanced constitutive and inducible phenolic defense mechanisms. The data may facilitate the selection of Scots pine half-sib families with low susceptibility for breeding programs and forest management strategies.

Association mapping of ectomycorrhizal traits in loblolly pine (Pinus taeda L.)

To understand how diverse mutualisms coevolve and how species adapt to complex environments, a description of the underlying genetic basis of the traits involved must be provided. For example, in diverse coevolving mutualisms, such as the interaction of host plants with a suite of symbiotic mycorrhizal fungi, a key question is whether host plants can coevolve independently with multiple species of symbionts, which depends on whether those interactions are governed independently by separate genes or pleiotropically by shared genes. To provide insight into this question, we employed an association mapping approach in a clonally replicated field experiment of loblolly pine (Pinus taeda L.) to identify genetic components of host traits governing ectomycorrhizal (EM) symbioses (mycorrhizal traits). The relative abundances of different EM fungi as well as the total number of root tips per cm root colonized by EM fungi were analyzed as separate mycorrhizal traits of loblolly pine. Single-nucleotide polymorphisms (SNPs) within candidate genes of loblolly pine were associated with loblolly pine mycorrhizal traits, mapped to the loblolly pine genome, and their putative protein function obtained when available. The results support the hypothesis that ectomycorrhiza formation is governed by host genes of large effect that apparently have independent influences on host interactions with different symbiont species.

Genetic architecture underlying the lignin biosynthesis pathway involves noncoding RNAs and transcription factors for growth and wood properties in Populus

Lignin provides structural support in perennial woody plants and is a complex phenolic polymer derived from phenylpropanoid pathway. Lignin biosynthesis is regulated by coordinated networks involving transcription factors (TFs), microRNAs (miRNAs) and long noncoding RNAs (lncRNAs). However, the genetic networks underlying the lignin biosynthesis pathway for tree growth and wood properties remain unknown. Here, we used association genetics (additive, dominant and epistasis) and expression quantitative trait nucleotide (eQTN) mapping to decipher the genetic networks for tree growth and wood properties in 435 unrelated individuals of Populus tomentosa. We detected 124 significant associations (P ≤ 6.89E-05) for 10 growth and wood property traits using 30 265 single nucleotide polymorphisms from 203 lignin biosynthetic genes, 81 TF genes, 36 miRNA genes and 71 lncRNA loci, implying their common roles in wood formation. Epistasis analysis uncovered 745 significant pairwise interactions, which helped to construct proposed genetic networks of lignin biosynthesis pathway and found that these regulators might affect phenotypes by linking two lignin biosynthetic genes. eQTNs were used to interpret how causal genes contributed to phenotypes. Lastly, we investigated the possible functions of the genes encoding 4-coumarate: CoA ligase and cinnamate-4-hydroxylase in wood traits using epistasis, eQTN mapping and enzymatic activity assays. Our study provides new insights into the lignin biosynthesis pathway in poplar and enables the novel genetic factors as biomarkers for facilitating genetic improvement of trees.

Genome-Wide Association Studies to Improve Wood Properties: Challenges and Prospects

Wood formation is an excellent model system for quantitative trait analysis due to the strong associations between the transcriptional and metabolic traits that contribute to this complex process. Investigating the genetic architecture and regulatory mechanisms underlying wood formation will enhance our understanding of the quantitative genetics and genomics of complex phenotypic variation. Genome-wide association studies (GWASs) represent an ideal statistical strategy for dissecting the genetic basis of complex quantitative traits. However, elucidating the molecular mechanisms underlying many favorable loci that contribute to wood formation and optimizing GWAS design remain challenging in this omics era. In this review, we summarize the recent progress in GWAS-based functional genomics of wood property traits in major timber species such as Eucalyptus, Populus, and various coniferous species. We discuss several appropriate experimental designs for extensive GWAS in a given undomesticated tree population, such as omics-wide association studies and high-throughput phenotyping technologies. We also explain why more attention should be paid to rare allelic and major structural variation. Finally, we explore the potential use of GWAS for the molecular breeding of trees. Such studies will help provide an integrated understanding of complex quantitative traits and should enable the molecular design of new cultivars.

Secondary Wall Regulating NACs Differentially Bind at the Promoter at a CELLULOSE SYNTHASE A4 Cis-eQTL

Arabidopsis thaliana CELLULOSE SYNTHASE A4/7/8 (CESA4/7/8) are three non-redundant subunits of the secondary cell wall cellulose synthase complex. Transcript abundance of these genes can vary among genotypes and expression quantitative trait loci (eQTL) were identified in a recombinant population of the accessions Bay-0 and Shahdara. Genetic mapping and analysis of the transcript levels of CESAs between two distinct near isogenic lines (NILs) confirmed a change in CESA4 expression that segregates within that interval. We sequenced the promoters and identified 16 polymorphisms differentiating CESA4^Sha and CESA4^Bay . In order to determine which of these SNPs could be responsible for this eQTL, we screened for transcription factor protein affinity with promoter fragments of CESA4^Bay, CESA4^Sha , and the reference genome CESA4^Col . The wall thickening activator proteins NAC SECONDARY WALL THICKENING PROMOTING FACTOR2 (NST2) and NST3 exhibited a decrease in binding with the CESA4^Sha promoter with a tracheary element-regulating cis-element (TERE) polymorphism. While NILs harboring the TERE polymorphisms exhibited significantly different CESA4 expression, cellulose crystallinity and cell wall thickness were indistinguishable. These results suggest that the TERE polymorphism resulted in differential transcription factor binding and CESA4 expression; yet A. thaliana is able to tolerate this transcriptional variability without compromising the structural elements of the plant, providing insight into the elasticity of gene regulation as it pertains to cell wall biosynthesis and regulation. We also explored available DNA affinity purification sequencing data to resolve a core binding site, C(G/T)TNNNNNNNA(A/C)G, for secondary wall NACs referred to as the VNS element.

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.

Genetic variation in transcription factors and photosynthesis light-reaction genes regulates photosynthetic traits

Transcription factors (TFs) play crucial roles in regulating the production of the components required for photosynthesis; elucidating the mechanisms by which underlying genetic variation in TFs affects complex photosynthesis-related traits may improve our understanding of photosynthesis and identify ways to improve photosynthetic efficiency. Promoter analysis of 96 nuclear-encoded Populus tomentosa Carr. genes within this pathway revealed 47 motifs responsive to light, stress, hormones and organ-specific regulation, as well as 86 TFs that might bind these motifs. Using phenotype-genotype associations, we identified 244 single-nucleotide polymorphisms (SNPs) within 105 genes associated with 12 photosynthesis-related traits. Most (30.33%) of these SNPs were located in intronic regions and these SNPs explained 18.66% of the mean phenotypic variation in the photosynthesis-related traits. Additionally, expression quantitative trait loci (eQTL) mapping identified 216 eQTLs associated with 110 eGenes (genes regulated by eQTLs), explaining 14.12% of the variability of gene expression. The lead SNPs of 12.04% of the eQTLs also contributed to phenotypic variation. Among these, a SNP in zf-Dof 5.6 (G120_9287) affected photosynthesis by modulating the expression of a sub-regulatory network of eight other TFs, which in turn regulate 55 photosynthesis-related genes. Furthermore, epistasis analysis identified a large interacting network representing 732 SNP-SNP pairs, of which 354 were photosynthesis gene-TF pairs, emphasizing the important roles of TFs in affecting photosynthesis-related traits. We combined eQTL and epistasis analysis and found 32 TFs harboring eQTLs being epistatic to their targets (identified by eQTL analysis), of which 15 TFs were also associated with photosynthesis traits. We therefore constructed a schematic model of TFs involved in regulating the photosynthetic light reaction pathway. Taken together, our results provide insight into the genetic regulation of photosynthesis, and may drive progress in the marker-assisted selection of desirable P. tomentosa genotypes with more efficient photosynthesis.

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.

Single-nucleotide polymorphisms(SNPs) in a sucrose synthase gene are associated with wood properties in Catalpa fargesii bur

BACKGROUND

Association study is a powerful means for identifying molecular markers, such as single-nucleotide polymorphisms (SNPs) associated with important traits in forest trees. Catalpa fargesii Bur is a valuable commercial tree in China and identifying SNPs that associate with wood property would make a foundation of the marker-assisted breeding in the future. However, related work has not been reported yet.

RESULTS

We cloned a 2887 bp long sucrose synthase (SUS) gene from the genome of C. fargesii, which is a key enzyme in sucrose metabolism and also associated to wood formation in trees, coding 806 amino acids that expressed mainly in young branches, xylem, and leaves according to real-time quantitative PCR. Then we identified allelic variations of CfSUS associated with nine wood quality associated traits in Catalpa fargesii Bur. Totally, 135 SNPs were identified through cloning and sequencing the CfSUS locus from a mapping population (including 93 unrelated individuals) and 47 of which were genotyped as common SNPs (minor allele frequency > 5%) in the association population that comprised of 125 unrelated individuals collected from main distribution area. Nucleotide diversity and linkage disequilibrium (LD) analysis showed CfSUS has a relative low SNP diversity (πT = 0.0034) and low LD (r2 dropped below 0.1 within 1600 bp). Using the association analysis, we found 11 common SNPs and 14 haplotypes were significantly associated with the traits (false discovery rate, Q<0.1), explaining 3.21-12.41% of the phenotypic variance. These results provide molecular markers above associated with wood basic density, pore rate, and six other traits of wood, which have potential applications in breeding of Catalpa fargesii Bur.

CONCLUSION

We first cloned a SUS gene in C. fargesii, then identified several SNPs and haplotypes that associated with wood properties within this gene, suggesting CfSUS participates in the wood formation of C. fargesii. Moreover, molecular markers we identified in this study may be applied into marker-assisted breeding of C. fargesii in the future.

Association genetics of acetophenone defence against spruce budworm in mature white spruce

BACKGROUND

Outbreaks of spruce budworm (SBW, Choristoneura fumiferana Clem.) cause major recurrent damage in boreal conifers such as white spruce (Picea glauca [Moench] Voss) and large losses of forest biomass in North America. Although defensive phenolic compounds have recently been linked to chemical resistance against SBW, their genetic basis remains poorly understood in forest trees, especially in conifers. Here, we used diverse association genetics approaches to discover genes and their variants that may control the accumulation of acetophenones, and dissect the genetic architecture of these defence compounds against SBW in white spruce mature trees.

RESULTS

Out of 4747 single nucleotide polymorphisms (SNPs) from 2312 genes genotyped in a population of 211 unrelated individuals, genetic association analyses identified 35 SNPs in 33 different genes that were significantly associated with the defence traits by using single-locus, multi-locus and multi-trait approaches. The multi-locus approach was particularly effective at detecting SNP-trait associations that explained a large fraction of the phenotypic variance (from 20 to 43%). Significant genes were regulatory including the NAC transcription factor, or they were involved in carbohydrate metabolism, falling into the binding, catalytic or transporter activity functional classes. Most of them were highly expressed in foliage. Weak positive phenotypic correlations were observed between defence and growth traits, indicating little or no evidence of defence-growth trade-offs.

CONCLUSIONS

This study provides new insights on the genetic architecture of tree defence traits, contributing to our understanding of the physiology of resistance mechanisms to biotic factors and providing a basis for the genetic improvement of the constitutive defence of white spruce against SBW.

Genetic variants in microRNA biogenesis genes as novel indicators for secondary growth in Populus

MicroRNAs (miRNAs) function as key regulators of complex traits, but how genetic alterations in miRNA biogenesis genes (miRBGs) affect quantitative variation has not been elucidated. We conducted transcript analyses and association genetics to investigate how miRBGs, miRNA genes (MIRNAs) and their respective targets contribute to secondary growth in a natural population of 435 Populus tomentosa individuals. This analysis identified 29 843 common single-nucleotide polymorphisms (SNPs; frequency > 0.10) within 682 genes (80 miRBGs, 152 MIRNAs, and 457 miRNA targets). Single-SNP association analysis found SNPs in 234 candidate genes exhibited significant additive/dominant effects on phenotypes. Among these, specific candidates that associated with the same traits produced 791 miRBG-MIRNA-target combinations, suggesting possible genetic miRBG-MIRNA and MIRNA-target interactions, providing an important clue for the regulatory mechanisms of miRBGs. Multi-SNP association found 4672 epistatic pairs involving 578 genes that showed significant associations with traits and identified 106 miRBG-MIRNA-target combinations. Two multi-hierarchical networks were constructed based on correlations of miRBG-miRNA and miRNA-target expression to further probe the mechanisms of trait diversity underlying changes in miRBGs. Our study opens avenues for the investigation of miRNA function in perennial plants and underscored miRBGs as potentially modulating quantitative variation in traits.

Identification of novel putative causative genes and genetic marker for male sterility in Japanese cedar (Cryptomeria japonica D.Don)

Background

Japanese cedar (Cryptomeria japonica) is an important tree for Japanese forestry. Male-sterile marker development in Japanese cedar would facilitate selection of male-sterile plus trees, addressing the widespread social problem of pollinosis and facilitating the identification of heterozygotes, which are useful for breeding.

Results

This study used next-generation sequencing for single-nucleotide polymorphism discovery in libraries constructed from several organs, including male-sterile and male-fertile strobili. The single-nucleotide polymorphisms obtained were used to construct a high-density linkage map, which enabled identification of a locus on linkage group 9 strongly correlated with male-sterile trait. Expressed sequence tags corresponding to 11 marker loci from 5 isotigs were associated with this locus within 33.4-34.5 cM. These marker loci explained 100% of the phenotypic variation. Several homologs of these sequences are associated with male sterility in rice or Arabidopsis, including a pre-mRNA splicing factor, a DEAD-box protein, a glycosyl hydrolase, and a galactosyltransferase. These proteins are thus candidates for the causal male-sterile gene at the ms-1 locus. After we used a SNaPshot assay to develop markers for marker-assisted selection (MAS), we tested F₂ progeny between male-sterile and wild-type plus trees to validate the markers and extrapolated the testing to a larger plus-tree population. We found that two developed from one of the candidates for the causal gene were suitable for MAS.

Conclusions

More than half of the ESTs and SNPs we collected were new, enlarging the genomic basis for genetic research on Japanese cedar. We developed two SNP markers aimed at MAS that distinguished individuals carrying the male-sterile trait with 100% accuracy, as well as individuals heterozygous at the male-sterile locus, even outside the mapping population. These markers should enable practical MAS for conifer breeding.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4581-5) contains supplementary material, which is available to authorized users.

Transcriptional Roadmap to Seasonal Variation in Wood Formation of Norway Spruce

Seasonal cues influence several aspects of the secondary growth of tree stems, including cambial activity, wood chemistry, and transition to latewood formation. We investigated seasonal changes in cambial activity, secondary cell wall formation, and tracheid cell death in woody tissues of Norway spruce (Picea abies) throughout one seasonal cycle. RNA sequencing was performed simultaneously in both the xylem and cambium/phloem tissues of the stem. Principal component analysis revealed gradual shifts in the transcriptomes that followed a chronological order throughout the season. A notable remodeling of the transcriptome was observed in the winter, with many genes having maximal expression during the coldest months of the year. A highly coexpressed set of monolignol biosynthesis genes showed high expression during the period of secondary cell wall formation as well as a second peak in midwinter. This midwinter peak in expression did not trigger lignin deposition, as determined by pyrolysis-gas chromatography/mass spectrometry. Coexpression consensus network analyses suggested the involvement of transcription factors belonging to the ASYMMETRIC LEAVES2/LATERAL ORGAN BOUNDARIES and MYELOBLASTOSIS-HOMEOBOX families in the seasonal control of secondary cell wall formation of tracheids. Interestingly, the lifetime of the latewood tracheids stretched beyond the winter dormancy period, correlating with a lack of cell death-related gene expression. Our transcriptomic analyses combined with phylogenetic and microscopic analyses also identified the cellulose and lignin biosynthetic genes and putative regulators for latewood formation and tracheid cell death in Norway spruce, providing a toolbox for further physiological and functional assays of these important phase transitions.

Drought Sensitivity of Norway Spruce at the Species' Warmest Fringe: Quantitative and Molecular Analysis Reveals High Genetic Variation Among and Within Provenances

Norway spruce (Picea abies) is by far the most important timber species in Europe, but its outstanding role in future forests is jeopardized by its high sensitivity to drought. We analyzed drought response of Norway spruce at the warmest fringe of its natural range. Based on a 35-year old provenance experiment we tested for genetic variation among and within seed provenances across consecutively occurring strong drought events using dendroclimatic time series. Moreover, we tested for associations between ≈1,700 variable SNPs and traits related to drought response, wood characteristics and climate-growth relationships. We found significant adaptive genetic variation among provenances originating from the species’ Alpine, Central and Southeastern European range. Genetic variation between individuals varied significantly among provenances explaining up to 44% of the phenotypic variation in drought response. Varying phenotypic correlations between drought response and wood traits confirmed differences in selection intensity among seed provenances. Significant associations were found between 29 SNPs and traits related to drought, climate-growth relationships and wood properties which explained between 11 and 43% of trait variation, though 12 of them were due to single individuals having extreme phenotypes of the respective trait. The majority of these SNPs are located within exons of genes and the most important ones are preferentially expressed in cambium and xylem expansion layers. Phenotype-genotype associations were stronger if only provenances with significant quantitative genetic variation in drought response were considered. The present study confirms the high adaptive variation of Norway spruce in Central and Southeastern Europe and demonstrates how quantitative genetic, dendroclimatic and genomic data can be linked to understand the genetic basis of adaptation to climate extremes in trees.

Genomic prediction accuracies in space and time for height and wood density of Douglas-fir using exome capture as the genotyping platform

Background

Genomic selection (GS) can offer unprecedented gains, in terms of cost efficiency and generation turnover, to forest tree selective breeding; especially for late expressing and low heritability traits. Here, we used: 1) exome capture as a genotyping platform for 1372 Douglas-fir trees representing 37 full-sib families growing on three sites in British Columbia, Canada and 2) height growth and wood density (EBVs), and deregressed estimated breeding values (DEBVs) as phenotypes. Representing models with (EBVs) and without (DEBVs) pedigree structure. Ridge regression best linear unbiased predictor (RR-BLUP) and generalized ridge regression (GRR) were used to assess their predictive accuracies over space (within site, cross-sites, multi-site, and multi-site to single site) and time (age-age/ trait-trait).

Results

The RR-BLUP and GRR models produced similar predictive accuracies across the studied traits. Within-site GS prediction accuracies with models trained on EBVs were high (RR-BLUP: 0.79–0.91 and GRR: 0.80–0.91), and were generally similar to the multi-site (RR-BLUP: 0.83–0.91, GRR: 0.83–0.91) and multi-site to single-site predictive accuracies (RR-BLUP: 0.79–0.92, GRR: 0.79–0.92). Cross-site predictions were surprisingly high, with predictive accuracies within a similar range (RR-BLUP: 0.79–0.92, GRR: 0.78–0.91). Height at 12 years was deemed the earliest acceptable age at which accurate predictions can be made concerning future height (age-age) and wood density (trait-trait). Using DEBVs reduced the accuracies of all cross-validation procedures dramatically, indicating that the models were tracking pedigree (family means), rather than marker-QTL LD.

Conclusions

While GS models’ prediction accuracies were high, the main driving force was the pedigree tracking rather than LD. It is likely that many more markers are needed to increase the chance of capturing the LD between causal genes and markers.

Quantitative trait locus mapping of Populus bark features and stem diameter

BACKGROUND

Bark plays important roles in photosynthate transport and storage, along with physical and chemical protection. Bark texture varies extensively among species, from smooth to fissured to deeply furrowed, but its genetic control is unknown. This study sought to determine the main genomic regions associated with natural variation in bark features and stem diameter. Quantitative trait loci (QTL) were mapped using an interspecific pseudo-backcross pedigree (Populus trichocarpa x P. deltoides and P. deltoides) for bark texture, bark thickness and diameter collected across three years, two sites and three biological replicates per site.

RESULTS

QTL specific to bark texture were highly reproducible in shared intervals across sites, years and replicates. Significant positive correlations and co-localization between trait QTL suggest pleiotropic regulators or closely linked genes. A list of candidate genes with related putative function, location close to QTL maxima and with the highest expression level in the phloem, xylem and cambium was identified.

CONCLUSION

Candidate genes for bark texture included an ortholog of Arabidopsis ANAC104 (PopNAC128), which plays a role in lignified fiber cell and ray development, as well as Pinin and Fasciclin (PopFLA) genes with a role in cell adhesion, cell shape and migration. The results presented in this study provide a basis for future genomic characterization of genes found within the QTL for bark texture, bark thickness and diameter in order to better understand stem and bark development in Populus and other woody perennial plants. The QTL mapping approach identified a list of prime candidate genes for further validation using functional genomics or forward genetics approaches.

Genetic Characterization of Chinese fir from Six Provinces in Southern China and Construction of a Core Collection

Large ex situ germplasm collections of plants generally contain significant diversity. A set of 700 well-conserved Chinese fir (Cunninghamia lanceolata (Lamb.) Hook) clones from six provinces in southern China in the ex situ gene bank of Longshan State Forest, was analyzed using 21 simple sequence repeat markers, with the aim of assessing the genetic diversity of these germplasm resources. Genetic analysis revealed extensive genetic variation among the accessions, with an average of 8.31 alleles per locus and a mean Shannon index of 1.331. Excluding loci with null alleles, we obtained a low level of genetic differentiation among provinces, consistent with the interpopulation genetic variation (1%). Three clusters were identified by STRUCTURE, which did not match the individuals’ geographical provenances. Ten traits related to growth and wood properties were quantified in these individuals, and there was substantial variation in all traits across individuals, these provide a potential source of variation for genetic improvement of the Chinese fir. Screening large collections for multiple-trait selective breeding programs is laborious and expensive; a core collection of 300 accessions, representative of the germplasm, was established, based on genotypic and phenotypic data. The identified small, but diverse, collections will be useful for further genome-wide association studies.

Association genetics of phenolic needle compounds in Norway spruce with variable susceptibility to needle bladder rust

KEY MESSAGE

Accumulation of phenolic needle metabolites in Norway spruce is regulated by many genes with small and additive effects and is correlated with the susceptibility against fungal attack. Norway spruce accumulates high foliar concentrations of secondary phenolic metabolites, with important functions for pathogen defence responses. However, the molecular genetic basis underlying the quantitative variation of phenolic compounds and their role in enhanced resistance of spruce to infection by needle bladder rust are unknown. To address these questions, a set of 1035 genome-wide single nucleotide polymorphisms (SNPs) was associated to the quantitative variation of four simple phenylpropanoids, eight stilbenes, nine flavonoids, six related arithmetic parameters and the susceptibility to infection by Chrysomyxa rhododendri in an unstructured natural population of Norway spruce. Thirty-one significant genetic associations for the flavonoids gallocatechin, kaempferol 3-glucoside and quercetin 3-glucoside and the stilbenes resveratrol, piceatannol, astringin and isorhapontin were discovered, explaining 22-59% of phenotypic variation, and indicating a regulation of phenolic accumulation by many genes with small and additive effects. The phenolics profile differed between trees with high and low susceptibility to the fungus, underlining the importance of phenolic compounds in the defence mechanisms of Norway spruce to C. rhododendri. Results highlight the utility of association studies in non-model tree species and may enable marker-assisted selection of Norway spruce adapted to severe pathogen attack.

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

Background

Genomic selection (GS) uses information from genomic signatures consisting of thousands of genetic markers to predict complex traits. As such, GS represents a promising approach to accelerate tree breeding, which is especially relevant for the genetic improvement of boreal conifers characterized by long breeding cycles. In the present study, we tested GS in an advanced-breeding population of the boreal black spruce (Picea mariana [Mill.] BSP) for growth and wood quality traits, and concurrently examined factors affecting GS model accuracy.

Results

The study relied on 734 25-year-old trees belonging to 34 full-sib families derived from 27 parents and that were established on two contrasting sites. Genomic profiles were obtained from 4993 Single Nucleotide Polymorphisms (SNPs) representative of as many gene loci distributed among the 12 linkage groups common to spruce. GS models were obtained for four growth and wood traits. Validation using independent sets of trees showed that GS model accuracy was high, related to trait heritability and equivalent to that of conventional pedigree-based models. In forward selection, gains per unit of time were three times higher with the GS approach than with conventional selection. In addition, models were also accurate across sites, indicating little genotype-by-environment interaction in the area investigated. Using information from half-sibs instead of full-sibs led to a significant reduction in model accuracy, indicating that the inclusion of relatedness in the model contributed to its higher accuracies. About 500 to 1000 markers were sufficient to obtain GS model accuracy almost equivalent to that obtained with all markers, whether they were well spread across the genome or from a single linkage group, further confirming the implication of relatedness and potential long-range linkage disequilibrium (LD) in the high accuracy estimates obtained. Only slightly higher model accuracy was obtained when using marker subsets that were identified to carry large effects, indicating a minor role for short-range LD in this population.

Conclusions

This study supports the integration of GS models in advanced-generation tree breeding programs, given that high genomic prediction accuracy was obtained with a relatively small number of markers due to high relatedness and family structure in the population. In boreal spruce breeding programs and similar ones with long breeding cycles, much larger gain per unit of time can be obtained from genomic selection at an early age than by the conventional approach. GS thus appears highly profitable, especially in the context of forward selection in species which are amenable to mass vegetative propagation of selected stock, such as spruces.

Allelic Interactions among Pto-MIR475b and Its Four Target Genes Potentially Affect Growth and Wood Properties in Populus

MicroRNAs (miRNAs) play crucial roles in plant growth and development, but few studies have illuminated the allelic interactions among miRNAs and their targets in perennial plants. Here, we combined analysis of expression patterns and single-nucleotide polymorphism (SNP)-based association studies to explore the interactions between Pto-MIR475b and its four target genes (Pto-PPR1, Pto-PPR2, Pto-PPR3, and Pto-PPR4) in 435 unrelated individuals of Populus tomentosa. Expression patterns showed a significant negative correlation (r = -0.447 to -0.411, P < 0.01) between Pto-MIR475b and its four targets in eight tissues of P. tomentosa, suggesting that Pto-miR475b may negatively regulate the four targets. Single SNP-based association studies identified 93 significant associations (P < 0.01, Q < 0.1) representing associations of 80 unique SNPs in Pto-MIR475b and its four targets with nine traits, revealing their potential roles in tree growth and wood formation. Moreover, one common SNP in the precursor region significantly altered the secondary structure of the pre-Pto-miR475b and changed the expression level of Pto-MIR475b. Analysis of epistatic interactions identified 115 significant SNP-SNP associations (P < 0.01) representing 45 unique SNPs from Pto-MIR475b and its four targets for 10 traits, revealing that genetic interactions between Pto-MIR475b and its targets influence quantitative traits of perennial plants. Our study provided a feasible strategy to study population genetics in forest trees and enhanced our understanding of miRNAs by dissecting the allelic interactions between this miRNA and its targets in P. tomentosa.

Linkage and Association Mapping for Two Major Traits Used in the Maritime Pine Breeding Program: Height Growth and Stem Straightness

Background

Increasing our understanding of the genetic architecture of complex traits, through analyses of genotype-phenotype associations and of the genes/polymorphisms accounting for trait variation, is crucial, to improve the integration of molecular markers into forest tree breeding. In this study, two full-sib families and one breeding population of maritime pine were used to identify quantitative trait loci (QTLs) for height growth and stem straightness, through linkage analysis (LA) and linkage disequilibrium (LD) mapping approaches.

Results

The populations used for LA consisted of two unrelated three-generation full-sib families (n = 197 and n = 477). These populations were assessed for height growth or stem straightness and genotyped for 248 and 217 markers, respectively. The population used for LD mapping consisted of 661 founders of the first and second generations of the breeding program. This population was phenotyped for the same traits and genotyped for 2,498 single-nucleotide polymorphism (SNP) markers corresponding to 1,652 gene loci. The gene-based reference genetic map of maritime pine was used to localize and compare the QTLs detected by the two approaches, for both traits. LA identified three QTLs for stem straightness and two QTLs for height growth. The LD study yielded seven significant associations (P ≤ 0.001): four for stem straightness and three for height growth. No colocalisation was found between QTLs identified by LA and SNPs detected by LD mapping for the same trait.

Conclusions

This study provides the first comparison of LA and LD mapping approaches in maritime pine, highlighting the complementary nature of these two approaches for deciphering the genetic architecture of two mandatory traits of the breeding program.

Association mapping for floral traits in cultivated Paeonia rockii based on SSR markers

Tree peony (Paeonia Sect. Moutan) is an economically important ornamental plant, but little is known about the genetic architecture of important ornamental traits. To effectively improve ornamental value, we require a better understanding of genetic architecture in the complex traits of the tree peony. Association mapping is a powerful tool for detection of variation associated with traits. Thus, we examined the genetic diversity and the population structure of 462 unrelated cultivated P. rockii individuals, then performed association mapping to identify simple sequence repeat (SSR) markers associated with 12 floral traits. We observed a moderate level of genetic diversity (PIC = 0.459) and low linkage disequilibrium (LD) between markers, demonstrating that the potential value of an LD approach in elucidating the molecular basis of the quantitative variation in this species. An analysis of population structure revealed three subgroups in the association population. Subsequent single-marker association analysis identified 46 significant associations, involving the 11 traits with 37 SSRs. These loci explained a small proportion of the phenotypic variance, ranging from 2.68 to 23.97% (mean 5.50%). We also validated 15 of the 46 associations in a linkage mapping population of 159 individuals. Finally, five associations were further confirmed in the linkage mapping population, involving the four traits with four SSRs. These results can serve as a foundation for further analyses of the genetic architecture of floral traits, and the SSRs associated in this work have potential applications in marker-assisted breeding in tree peony.

Allelic variation in PtoPsbW associated with photosynthesis, growth, and wood properties in Populus tomentosa

Photosynthesis is one of the most important reactions on earth. PsbW, a nuclear-encoded subunit of photosystem II (PSII), stabilizes PSII structure and plays an important role in photosynthesis. Here, we used candidate gene-based linkage disequilibrium (LD) mapping to detect significant associations between allelic variations of PtoPsbW and traits related to photosynthesis, growth, and wood properties in Populus tomentosa. PtoPsbW showed the highest expression in leaves and it increased during the development of these leaves, suggesting that PtoPsbW may play an important role in plant growth and development. Analysis of nucleotide diversity and LD revealed that PtoPsbW has low single-nucleotide polymorphism (SNP) diversity (π _tot = 0.0048 and θ _w = 0.0050) and relatively low average value of LD (0.1500), indicating that PtoPsbW is conserved due to its indispensable function. Using single-SNP associations in an association population of 435 individuals, we identified five significant associations at the threshold of P ≤ 0.05, explaining 3.28-15.98 % of the phenotypic variation. Haplotype-based association analyses indicated that 13 haplotypes (P ≤ 0.05) from six blocks were associated with photosynthesis, growth, and wood properties. Our work shows that identifying allelic variation and LD can help to decipher the genetic basis of photosynthesis and could potentially be applied for molecular marker-assisted selection in Populus.

Genetic variations and miRNA-target interactions contribute to natural phenotypic variations in Populus

Variation in regulatory factors, including microRNAs (miRNAs), contributes to variation in quantitative and complex traits. However, in plants, variants in miRNAs and their target genes that contribute to natural phenotypic variation, and the underlying regulatory networks, remain poorly characterized. We investigated the associations and interactions of single-nucleotide polymorphisms (SNPs) in miRNAs and their target genes with phenotypes in 435 individuals from a natural population of Populus. We used RNA-seq to identify 217 miRNAs differentially expressed in a tension wood system, and identified 1196 candidate target genes; degradome sequencing confirmed 60 of the target sites. In addition, 72 miRNA-target pairs showed significant co-expression. Gene ontology (GO) term analysis showed that most of the genes in the co-regulated pairs participate in biological regulation. Genome resequencing found 5383 common SNPs (frequency ≥ 0.05) in 139 miRNAs and 31 037 SNPs in 819 target genes. Single-SNP association analyses identified 232 significant associations between wood traits (P ≤ 0.05) and SNPs in 102 miRNAs and 1387 associations with 478 target genes. Among these, 102 miRNA-target pairs associated with the same traits. Multi-SNP associations found 102 epistatic pairs associated with traits. Furthermore, a reconstructed regulatory network contained 12 significantly co-expressed pairs, including eight miRNAs and nine targets associated with traits. Lastly, both expression and genetic association showed that miR156i, miR156j, miR396a and miR6445b were involved in the formation of tension wood. This study shows that variants in miRNAs and target genes contribute to natural phenotypic variation and annotated roles and interactions of miRNAs and their target genes by genetic association analysis.

Morpho-Physiological Variation of White Spruce Seedlings from Various Seed Sources and Implications for Deployment under Climate Change

Because of changes in climatic conditions, tree seeds originating from breeding programs may no longer be suited to sites where they are currently sent. As a consequence, new seed zones may have to be delineated. Assisted migration consists of transferring seed sources that match the future climatic conditions to which they are currently adapted. It represents a strategy that could be used to mitigate the potential negative consequences of climate change on forest productivity. Decisions with regard to the choice of the most appropriate seed sources have to rely on appropriate knowledge of morpho-physiological responses of trees. To meet this goal, white spruce (Picea glauca [Moench] Voss) seedlings from eight seed orchards were evaluated during two years in a forest nursery, and at the end of the first growing season on three plantation sites located in different bioclimatic domains in Quebec. The morpho-physiological responses obtained at the end of the second growing season (2+0) in the nursery made it possible to cluster the orchards into three distinct groups. Modeling growth curves of these different groups showed that the height growth of seedlings from the second-generation and southern first-generation seed orchards was significantly higher than that of those from other orchards, by at least 6%. A multiple regression model with three climatic variables (average growing season temperature, average July temperature, length of the growing season) showed that the final height of seedlings (2+0) from the first-generation seed orchards was significantly related to the local climatic conditions at the orchard sites of origin where parental trees from surrounding natural populations were sampled to provide grafts for orchard establishment. Seedling height growth was significantly affected by both seed source origins and planting sites, but the relative ranking of the different seed sources was maintained regardless of reforestation site. This knowledge could be used, in conjunction with transfer models, to refine operational seed transfer rules and select the most suitable sites in an assisted migration strategy.

Dissection of allelic interactions among Pto-miR257 and its targets and their effects on growth and wood properties in Populus

MicroRNAs (miRNAs) have important roles in the regulation of genes; however, for trees few studies have explored the potential impact of the interactions between miRNAs and their target genes. Here, we performed transcript profiling and association genetics (single-SNP, haplotype-based and multi-SNP associations) to study the genetic regulatory relationship of Pto-miR257 and its 12 target genes in 435 individuals of a natural population of Populus tomentosa. Expression profiling of Pto-miR257 and its targets showed a negative relationship between their expression levels. Of the 61 single-nucleotide polymorphisms (SNPs) detected in Pto-miR257, 6 in the pre-mature region strongly affected its secondary stability and 1 in the mature region could alter its target spectrum. Among the 1029 SNPs in the targets, 3 were located in target sites that could change the binding affinity of Pto-miR257. Single-SNP association analysis revealed that SNPs in Pto-miR257 and target genes associated with both growth and wood property traits, in agreement with haplotype-based identifications. Multi-SNP association found that 10 targets shared at least one common trait with Pto-miR257, with phenotypic variance from 0.5 to 8.5%, suggesting a possible internal genetic interaction between them. Epistasis analysis showed significant epistatic interactions among Pto-miR257 and its targets. Therefore, our study demonstrated Pto-miR257 and its 12 targets had roles in wood formation and revealed the genetic interaction network between the miRNA and its targets under additive, dominant and epistatic models. Thus, association genetics can be used to decipher the interactions between miRNAs and their target genes and to help understand the genetic architecture of complex traits.

Genetic architecture of wood properties based on association analysis and co-expression networks in white spruce

Association studies are widely utilized to analyze complex traits but their ability to disclose genetic architectures is often limited by statistical constraints, and functional insights are usually minimal in nonmodel organisms like forest trees. We developed an approach to integrate association mapping results with co-expression networks. We tested single nucleotide polymorphisms (SNPs) in 2652 candidate genes for statistical associations with wood density, stiffness, microfibril angle and ring width in a population of 1694 white spruce trees (Picea glauca). Associations mapping identified 229-292 genes per wood trait using a statistical significance level of P < 0.05 to maximize discovery. Over-representation of genes associated for nearly all traits was found in a xylem preferential co-expression group developed in independent experiments. A xylem co-expression network was reconstructed with 180 wood associated genes and several known MYB and NAC regulators were identified as network hubs. The network revealed a link between the gene PgNAC8, wood stiffness and microfibril angle, as well as considerable within-season variation for both genetic control of wood traits and gene expression. Trait associations were distributed throughout the network suggesting complex interactions and pleiotropic effects. Our findings indicate that integration of association mapping and co-expression networks enhances our understanding of complex wood traits.

Implementation of the Realized Genomic Relationship Matrix to Open-Pollinated White Spruce Family Testing for Disentangling Additive from Nonadditive Genetic Effects

The open-pollinated (OP) family testing combines the simplest known progeny evaluation and quantitative genetics analyses as candidates’ offspring are assumed to represent independent half-sib families. The accuracy of genetic parameter estimates is often questioned as the assumption of “half-sibling” in OP families may often be violated. We compared the pedigree- vs. marker-based genetic models by analysing 22-yr height and 30-yr wood density for 214 white spruce [Picea glauca (Moench) Voss] OP families represented by 1694 individuals growing on one site in Quebec, Canada. Assuming half-sibling, the pedigree-based model was limited to estimating the additive genetic variances which, in turn, were grossly overestimated as they were confounded by very minor dominance and major additive-by-additive epistatic genetic variances. In contrast, the implemented genomic pairwise realized relationship models allowed the disentanglement of additive from all nonadditive factors through genetic variance decomposition. The marker-based models produced more realistic narrow-sense heritability estimates and, for the first time, allowed estimating the dominance and epistatic genetic variances from OP testing. In addition, the genomic models showed better prediction accuracies compared to pedigree models and were able to predict individual breeding values for new individuals from untested families, which was not possible using the pedigree-based model. Clearly, the use of marker-based relationship approach is effective in estimating the quantitative genetic parameters of complex traits even under simple and shallow pedigree structure.

Genetic Adaptation to Climate in White Spruce Involves Small to Moderate Allele Frequency Shifts in Functionally Diverse Genes

Understanding the genetic basis of adaptation to climate is of paramount importance for preserving and managing genetic diversity in plants in a context of climate change. Yet, this objective has been addressed mainly in short-lived model species. Thus, expanding knowledge to nonmodel species with contrasting life histories, such as forest trees, appears necessary. To uncover the genetic basis of adaptation to climate in the widely distributed boreal conifer white spruce (Picea glauca), an environmental association study was conducted using 11,085 single nucleotide polymorphisms representing 7,819 genes, that is, approximately a quarter of the transcriptome.

Linear and quadratic regressions controlling for isolation-by-distance, and the Random Forest algorithm, identified several dozen genes putatively under selection, among which 43 showed strongest signals along temperature and precipitation gradients. Most of them were related to temperature. Small to moderate shifts in allele frequencies were observed. Genes involved encompassed a wide variety of functions and processes, some of them being likely important for plant survival under biotic and abiotic environmental stresses according to expression data. Literature mining and sequence comparison also highlighted conserved sequences and functions with angiosperm homologs.

Our results are consistent with theoretical predictions that local adaptation involves genes with small frequency shifts when selection is recent and gene flow among populations is high. Accordingly, genetic adaptation to climate in P. glauca appears to be complex, involving many independent and interacting gene functions, biochemical pathways, and processes. From an applied perspective, these results shall lead to specific functional/association studies in conifers and to the development of markers useful for the conservation of genetic resources.

Association mapping in Populus reveals the interaction between Pto-miR530a and its target Pto-KNAT1

We used transcript profiling and multi-SNP association to investigate the genetic regulatory relationship between miRNA Pto-miR530a and its target Pto-KNAT1, identifying additive, dominant, and epistatic effects. MicroRNAs (miRNAs) play crucial roles in the post-transcriptional regulation of plant growth and development; indeed, many studies have described the importance of miRNA-target interactions in herbaceous species. However, elucidation of the miRNA-target interactions in trees may require novel strategies. In the present study, we describe a strategy combining expression profiling by reverse transcription quantitative PCR (RT-qPCR) and association mapping with multiple single nucleotide polymorphisms (SNPs) to evaluate the interaction between Pto-miR530a and its target Pto-KNAT1 in Populus tomentosa. RT-qPCR analysis showed a negative correlation (r = -0.62, P < 0.05) between expression levels of Pto-miR530a and Pto-KNAT1 in eight tissues. We used a Bayesian hierarchical model to identify allelic variants of Pto-miR530a and Pto-KNAT1 that associated with eight traits related to growth and wood properties, in a population of 460 unrelated individuals of P. tomentosa. This analysis identified 27 associations, with the proportions of phenotypic variance (R (2)) contributed by each SNP ranging of 0.82-15.81 %, the additive effects of each SNP ranging of 0.16-18.09, and the dominant effects ranging from -14.09 to 19.00. Epistatic interaction models showed a strong interaction among SNPs in the miRNA target with R (2) of 0.1-3.56 %, and information gain of significant SNP pairs of -3.09 to 0.93 %, representing the regulatory interactions between the miRNA and the mRNA. Thus, we used a new strategy that combines association genetics and expression profiling based on SNPs to study the regulatory relationship between this miRNA and its target mRNA, thereby providing novel advances in our understanding of the genetic architecture of important traits.

Association genetics in Populus reveals the interactions between Pt-miR397a and its target genes

Recent studies have revealed associations between single nucleotide polymorphisms (SNPs) in microRNA (miRNA) genes and diseases. However, association studies to decipher the interactions between miRNAs and their target genes remain to be conducted. Here, we investigated the association of growth and wood traits with SNPs in Pt-miR397a and its targets, in 261 individuals from a natural population of Populus tomentosa. Of the 57 SNPs identified in Pt-miR397a, three strongly affect its secondary stability, and SNPs in target sites in Pt-LAC20 and Pt-HSP40 changed the binding affinity of Pt-miR397a. Single-SNP association analysis revealed that SNPs in Pt-miR397a significantly associated with α-cellulose content and stem volume, and SNPs in target genes also associated with growth and wood-property traits. Multi-SNP association analysis with additive and dominant models found that SNPs in six potential target genes associated with at least one trait in common with Pt-miR397a, revealing a possible genetic interaction between Pt-miR397a and its targets. Furthermore, epistasis analysis revealed epistatic interactions between SNPs in Pt-miR397a and its target genes. Thus, our study indicated that SNPs in Pt-miR397a and six target genes affect wood formation and that association studies can reveal the interactions between miRNAs and their target genes.