Yousry A. El-Kassaby

Single-copy, species-transferable microsatellite markers developed from loblolly pine ESTs

Microsatellites, or simple sequence repeats (SSRs), are usually regarded as the “markers of choice” in population genetics research because they exhibit high variability. The development cost of these markers is usually high. In addition, microsatellite primers developed for one species often do not cross-amplify in related species, requiring separate development for each species. However, microsatellites found in expressed sequence tags (ESTs) might better cross-amplify as they reside in or near conserved coding DNA. In this study, we identified 14 Pinus taeda (loblolly pine) EST-SSRs from public EST databases and tested for their cross-species transferability to P. contorta ssp. latifolia, P. ponderosa, and P. sylvestris. As part of our development of a P. contorta microsatellite set, we also compared their transferability to that of 99 traditional microsatellite markers developed in P. taeda and tested on P. contorta ssp. latifolia. Compared to traditional microsatellites, EST-SSRs had higher transfer rates across pine species; however, the level of polymorphism of microsatellites derived from ESTs was lower. Sequence analyses revealed that the frequencies of insertions/deletions and base substitutions were lower in EST-SSRs than in other types of microsatellites, confirming that EST-SSRs are more conserved than traditional SSRs. Our results also provide a battery of 23 polymorphic, robust microsatellite primer pairs for lodgepole pine.

The nature of inbreeding in a seed orchard of Douglas fir as shown by an efficient multilocus model

SummaryThe amounts of self-fertilization versus consanguineous matings (as measured by effective selfing) was estimated in a seed orchard of Douglas-fir, using progeny array data at six allozyme loci. The orchard is family structured, consisting several grafts (clones) and/or open-pollinated (o-p) progeny from each of several ‘plus-trees’. Population-wide selfing rates were found to be 7% for the o-p trees and 2% for the cloned trees. Estimates of mating system parameters for individual trees showed this difference for average outcrossing rate t (1) still largely remained when outcrossingpollen gene frequency p was not allowed to vary among trees and (2) disappeared when p was allowed to vary among trees. Under this joint t and p estimation, o-p trees showed both significant variation of t (based upon a one-way ANOVA grouped by common plus-tree) and significant regressions of p on ovule genotype (indicative of consanguineous matings); cloned trees showed neither. This higher rate of consanguineous mating for o-p trees might be explained by the larger and more variable size of o-p families in the orchard. Estimates of outcrossing rate t and outcrossingpollen gene frequency p were based upon a multilocus model which makes full use of the information in the data. The increased information it gives over ‘observed outcross’ models is equivalent to adding 30–50% more loci, and it gives enough degrees of freedom to jointly estimate t and p for individual trees (individual progeny arrays) under certain conditions. In addition, inclusion of megagametophyte data nearly doubles the information about the mating system of individual trees.

Genome resequencing reveals multiscale geographic structure and extensive linkage disequilibrium in the forest tree Populus trichocarpa

• 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.

Geographical and environmental gradients shape phenotypic trait variation and genetic structure in Populus trichocarpa

• Populus 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. H(2) ranged from 0.3 to 0.9 in phenology, 0.3 to 0.8 in biomass and 0.1 to 0.8 in ecophysiology traits. Most traits correlated strongly with latitude, maximum daylength and temperature of tree origin, but not necessarily with elevation, precipitation or heat : moisture indices. Trait H(2) 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.

Inheritance of chloroplast and mitochondrial DNA in Picea and composition of hybrids from introgression zones

SummaryThe cloning of white spruce (Picea glauca) mitochondrial DNA homologous to the cytochrome oxidase II and ATPaseα genes of maize is described. These probes were used to define restriction fragment length polymorphisms which distinguish the white, Engelmann (P. engelmannii) and Sitka spruce (P. sitchensis) populations that occur in British Columbia. Analysis of progeny from crosses between the species revealed that mitochondrial DNA was maternally inherited in all cases (32 progeny from five independent crosses). The inheritance of chloroplast DNA was determined using a probe described previously; in this case, all progeny exhibited paternal inheritance (27 progeny from four crosses). Mitochondrial and chloroplast probes were used to test trees from zones of introgression between coastal (Sitka) and interior spruces (white and Engelmann). In most cases mitochondria and chloroplasts within individuals were contributed by different species. The data shows that there is a significant Sitka spruce component in trees east of the coastal watershed in British Columbia.

Genome‐wide association implicates numerous genes underlying ecological trait variation in natural populations of Populus trichocarpa

In 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×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.

Genome-wide association mapping for wood characteristics in Populus identifies an array of candidate single nucleotide polymorphisms.

Establishing 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 genome-wide associations. The analysis revealed 141 significant SNPs (α ≤ 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 traits 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.

A 34K SNP genotyping array for Populus trichocarpa: design, application to the study of natural populations and transferability to other Populus species.

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.

Allozyme inheritance, heterozygosity and outcrossing rate among Pinus monticola near Ladysmith, British Columbia

Analysis of megagametophytic and embryonic allozyme variants in nine enzymatic systems encoded by 14 loci was conducted on 30 western white pine (Pinus monticola Dougl.) trees from a natural stand on Vancouver Island, B.C. The segregation of allozymes in megagametophytes of heterozygous trees indicated distinct, simple, Mendelian inheritance. The stand heterozygosity parameters (proportion of polymorphic loci (0·64), average number of alleles per locus (1·79), and average expected heterozygosity (0·18)) were identical to those reported from a broadly-based survey of the species (Steinhoff et al., 1983). Comparison of progeny (viable embryo) and parental genotypic distributions indicated a major genetic shift between the two life-cycles phases. Among the progeny, heterozygotes occurred less frequently than expected under panmixia, while the reverse was observed for the parental population. Single-locus estimates of outcrossing rate varied between 0·751 and 1·043 and were significantly heterogeneous. Comparison between the single-locus estimates and the multilocus estimate (0·977) of outcrossing rates indicated that most of the inbreeding detected was due to consanguineous matings, rather than selfing. This form of inbreeding was consistent with the clustering of genotypes in the stand. Single-tree estimates of outcrossing rate varied considerably among trees and ranged between 0·682 and 1·207, indicating that some trees might possess high self compatibility and tolerate high selfing, or that violation of mixed-mating assumptions occurred. The implications of the apparent inbreeding on the derivation of genetic estimates from open-pollinated seeds and on seed crops for production plantations are discussed.

Populus trichocarpa cell wall chemistry and ultrastructure trait variation, genetic control and genetic correlations

The 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 Hardy-Weinberg 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.