SummaryPopulus trichocarpa is widespread across western North America spanning extensive variation in photoperiod, growing season and climate. We investigated trait variation in P. trichocarpa using over 2000 trees from a common garden at Vancouver, Canada, representing replicate plantings of 461 genotypes originating from 136 provenance localities.We measured 40 traits encompassing phenological events, biomass accumulation, growth rates, and leaf, isotope and gas exchange-based ecophysiology traits. With replicated plantings and 29 354 single nucleotide polymorphisms (SNPs) from 3518 genes, we estimated both broad-sense trait heritability (H 2 ) and overall population genetic structure from principal component analysis.Populus trichocarpa had high phenotypic variation and moderate/high H 2 for many traits. values reflected trait correlation strength with geoclimate variables. The population genetic structure had one significant principal component (PC1) which correlated with daylength and showed enrichment for genes relating to circadian rhythm and photoperiod. Robust relationships between traits, population structure and geoclimate in P. trichocarpa reflect patterns which suggest that range-wide geographical and environment gradients have shaped its genotypic and phenotypic variability.

Summary• Plant population genomics informs evolutionary biology, breeding, conservation and bioenergy feedstock development. For example, the detection of reliable phenotype-genotype associations and molecular signatures of selection requires a detailed knowledge about genome-wide patterns of allele frequency variation, linkage disequilibrium and recombination.• We resequenced 16 genomes of the model tree Populus trichocarpa and genotyped 120 trees from 10 subpopulations using 29 213 single-nucleotide polymorphisms.• Significant geographic differentiation was present at multiple spatial scales, and range-wide latitudinal allele frequency gradients were strikingly common across the genome. The decay of linkage disequilibrium with physical distance was slower than expected from previous studies in Populus, with r 2 dropping below 0.2 within 3-6 kb. Consistent with this, estimates of recent effective population size from linkage disequilibrium (N e 4000-6000) were remarkably low relative to the large census sizes of P. trichocarpa stands. Fine-scale rates of recombination varied widely across the genome, but were largely predictable on the basis of DNA sequence and methylation features.• Our results suggest that genetic drift has played a significant role in the recent evolutionary history of P. trichocarpa. Most importantly, the extensive linkage disequilibrium detected suggests that genome-wide association studies and genomic selection in undomesticated populations may be more feasible in Populus than previously assumed.

SummaryIn order to uncover the genetic basis of phenotypic trait variation, we used 448 unrelated wild accessions of black cottonwood (Populus trichocarpa) from much of its range in western North America. Extensive data from large-scale trait phenotyping (with spatial and temporal replications within a common garden) and genotyping (with a 34 K Populus single nucleotide polymorphism (SNP) array) of all accessions were used for gene discovery in a genome-wide association study (GWAS).We performed GWAS with 40 biomass, ecophysiology and phenology traits and 29 355 filtered SNPs representing 3518 genes. The association analyses were carried out using a Unified Mixed Model accounting for population structure effects among accessions.We uncovered 410 significant SNPs using a Bonferroni-corrected threshold (P < 1.7 9 10 À6). Markers were found across 19 chromosomes, explained 1-13% of trait variation, and implicated 275 unique genes in trait associations. Phenology had the largest number of associated genes (240 genes), followed by biomass (53 genes) and ecophysiology traits (25 genes).The GWAS results propose numerous loci for further investigation. Many traits had significant associations with multiple genes, underscoring their genetic complexity. Genes were also identified with multiple trait associations within and/or across trait categories. In some cases, traits were genetically correlated while in others they were not.

SummaryEstablishing links between phenotypes and molecular variants is of central importance to accelerate genetic improvement of economically important plant species. Our work represents the first genome-wide association study to the inherently complex and currently poorly understood genetic architecture of industrially relevant wood traits.Here, we employed an Illumina Infinium 34K single nucleotide polymorphism (SNP) genotyping array that generated 29 233 high-quality SNPs in c. 3500 broad-based candidate genes within a population of 334 unrelated Populus trichocarpa individuals to establish genomewide associations.The analysis revealed 141 significant SNPs (a ≤ 0.05) associated with 16 wood chemistry/ ultrastructure traits, individually explaining 3-7% of the phenotypic variance. A large set of associations (41% of all hits) occurred in candidate genes preselected for their suggested a priori involvement with secondary growth. For example, an allelic variant in the FRA8 ortholog explained 21% of the total genetic variance in fiber length, when the trait's heritability estimate was considered. The remaining associations identified SNPs in genes not previously implicated in wood or secondary wall formation.Our findings provide unique insights into wood trait architecture and support efforts for population improvement based on desirable allelic variants.

Stomata are essential for diffusive entry of gases to support photosynthesis, but may also expose internal leaf tissues to pathogens. To uncover trade-offs in range-wide adaptation relating to stomata, we investigated the underlying genetics of stomatal traits and linked variability in these traits with geoclimate, ecophysiology, condensed foliar tannins and pathogen susceptibility in black cottonwood (Populus trichocarpa). Upper (adaxial) and lower (abaxial) leaf stomatal traits were measured from 454 accessions collected throughout much of the species range. We calculated broad-sense heritability (H(2) ) of stomatal traits and, using SNP data from a 34K Populus SNP array, performed a genome-wide association studies (GWAS) to uncover genes underlying stomatal trait variation. H(2) values for stomatal traits were moderate (average H(2) = 0.33). GWAS identified genes associated primarily with adaxial stomata, including polarity genes (PHABULOSA), stomatal development genes (BRASSINOSTEROID-INSENSITIVE 2) and disease/wound-response genes (GLUTAMATE-CYSTEINE LIGASE). Stomatal traits correlated with latitude, gas exchange, condensed tannins and leaf rust (Melampsora) infection. Latitudinal trends of greater adaxial stomata numbers and guard cell pore size corresponded with higher stomatal conductance (gs ) and photosynthesis (Amax ), faster shoot elongation, lower foliar tannins and greater Melampsora susceptibility. This suggests an evolutionary trade-off related to differing selection pressures across the species range. In northern environments, more adaxial stomata and larger pore sizes reflect selection for rapid carbon gain and growth. By contrast, southern genotypes have fewer adaxial stomata, smaller pore sizes and higher levels of condensed tannins, possibly linked to greater pressure from natural leaf pathogens, which are less significant in northern ecosystems.

Genetic mapping of quantitative traits requires genotypic data for large numbers of markers in many individuals. For such studies, the use of large single nucleotide polymorphism (SNP) genotyping arrays still offers the most cost-effective solution. Herein we report on the design and performance of a SNP genotyping array for Populus trichocarpa (black cottonwood). This genotyping array was designed with SNPs pre-ascertained in 34 wild accessions covering most of the species latitudinal range. We adopted a candidate gene approach to the array design that resulted in the selection of 34 131 SNPs, the majority of which are located in, or within 2 kb of, 3543 candidate genes. A subset of the SNPs on the array (539) was selected based on patterns of variation among the SNP discovery accessions. We show that more than 95% of the loci produce high quality genotypes and that the genotyping error rate for these is likely below 2%. We demonstrate that even among small numbers of samples (n = 10) from local populations over 84% of loci are polymorphic. We also tested the applicability of the array to other species in the genus and found that the number of polymorphic loci decreases rapidly with genetic distance, with the largest numbers detected in other species in section Tacamahaca. Finally, we provide evidence for the utility of the array to address evolutionary questions such as intraspecific studies of genetic differentiation, species assignment and the detection of natural hybrids.

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.

SummaryThe increasing ecological and economical importance of Populus species and hybrids has stimulated research into the investigation of the natural variation of the species and the estimation of the extent of genetic control over its wood quality traits for traditional forestry activities as well as the emerging bioenergy sector. A realized kinship matrix based on informative, high-density, biallelic single nucleotide polymorphism (SNP) genetic markers was constructed to estimate trait variance components, heritabilities, and genetic and phenotypic correlations.Seventeen traits related to wood chemistry and ultrastructure were examined in 334 9-yr-old Populus trichocarpa grown in a common-garden plot representing populations spanning the latitudinal range 44°to 58.6°. In these individuals, 9342 SNPs that conformed to HardyWeinberg expectations were employed to assess the genomic pair-wise kinship to estimate narrow-sense heritabilities and genetic correlations among traits.The range-wide phenotypic variation in all traits was substantial and several trait heritabilities were > 0.6. In total, 61 significant genetic and phenotypic correlations and a network of highly interrelated traits were identified.The high trait variation, the evidence for moderate to high heritabilities and the identification of advantageous trait combinations of industrially important characteristics should aid in providing the foundation for the enhancement of poplar tree breeding strategies for modern industrial use.