Growing up aspen: ontogeny and trade-offs shape growth, defence and reproduction in a foundation species.

BACKGROUND AND AIMS: Intraspecific variation in foundation species of forest ecosystems can shape community and ecosystem properties, particularly when that variation has a genetic basis. Traits mediating interactions with other species are predicted by simple allocation models to follow ontogenetic patterns that are rarely studied in trees. The aim of this research was to identify the roles of genotype, ontogeny and genotypic trade-offs shaping growth, defence and reproduction in aspen. METHODS: We established a common garden replicating >500 aspen genets in Wisconsin, USA. Trees were measured through the juvenile period into the onset of reproduction, for growth, defence chemistry (phenolic glycosides and condensed tannins), nitrogen, extrafloral nectaries, leaf morphology (specific leaf area), flower production and foliar herbivory and disease. We also assayed the TOZ19 sex marker and heterozygosity at ten microsatellite loci. KEY RESULTS: We found high levels of genotypic variation for all traits, and high heritabilities for both the traits and their ontogenetic trajectories. Ontogeny strongly shaped intraspecific variation, and trade-offs among growth, defence and reproduction supported some predictions while contradicting others. Both direct resistance (chemical defence) and indirect defence (extrafloral nectaries) declined during the juvenile stage, prior to the onset of reproduction. Reproduction was higher in trees that were larger, male and had higher individual heterozygosity. Growth was diminished by genotypic allocation to both direct and indirect defence as well as to reproduction, but we found no evidence of trade-offs between defence and reproduction. CONCLUSIONS: Key traits affecting the ecological communities of aspen have high levels of genotypic variation and heritability, strong patterns of ontogeny and clear trade-offs among growth, defence and reproduction. The architecture of aspen's community genetics - its ontogeny, trade-offs and especially its great variability - is shaped by both its broad range and the diverse community of associates, and in turn further fosters that diversity.

Pathway position constrains the evolution of an ecologically important pathway in aspens (Populus tremula L.).

Many ecological interactions of aspens and their relatives (Populus spp.) are affected by products of the phenylpropanoid pathway synthesizing condensed tannins (CTs), whose production involves trade-offs with other ecologically important compounds and with growth. Genes of this pathway are candidates for investigating the role of selection on ecologically important, polygenic traits. We analysed sequences from 25 genes representing 10 steps of the CT synthesis pathway, which produces CTs used in defence and lignins used for growth, in 12 individuals of European aspen (Populus tremula). We compared these to homologs from P. trichocarpa, to a control set of 77 P. tremula genes, to genome-wide resequencing data and to RNA-seq expression levels, in order to identify signatures of selection distinct from those of demography. In Populus, pathway position exerts a strong influence on the evolution of these genes. Nonsynonymous diversity, divergence and allele frequency shifts (Tajima's D) were much lower than for synonymous measures. Expression levels were higher, and the direction of selection more negative, for upstream genes than for those downstream. Selective constraints act with increasing intensity on upstream genes, despite the presence of multiple paralogs in most gene families. Pleiotropy, expression level, flux control and codon bias appear to interact in determining levels and patterns of variation in genes of this pathway, whose products mediate a wide array of ecological interactions for this widely distributed species.

Heterozygosity, gender, and the growth-defense trade-off in quaking aspen.

Although plant growth is generally recognized to be influenced by allocation to defense, genetic background (e.g., inbreeding), and gender, rarely have those factors been addressed collectively. In quaking aspen (Populus tremuloides Michx.), phenolic glycosides (PGs) and condensed tannins (CTs) constitute up to 30 % of leaf dry weight. To quantify the allocation cost of this chemical defense, we measured growth, defense chemistry, and individual heterozygosity (H obs at 16 microsatellite loci) for male and female trees in both controlled and natural environments. The controlled environment consisted of 12 juvenile genets grown for 3 years in a common garden, with replication. The natural environment consisted of 51 mature genets in wild populations, from which we sampled multiple ramets (trees) per genet. Concentrations of PGs and CTs were negatively correlated. PGs were uncorrelated with growth, but CT production represented a major cost. Across the range of CT levels found in wild-grown trees, growth rates varied by 2.6-fold, such that a 10 % increase in CT concentration occurred with a 38.5 % decrease in growth. H obs had a marked effect on aspen growth: for wild trees, a 10 % increase in H obs corresponded to a 12.5 % increase in growth. In wild trees, this CT effect was significant only in females, in which reproduction seems to exacerbate the cost of defense, while the H obs effect was significant only in males. Despite the lower growth rate of low-H obs trees, their higher CT levels may improve survival, which could account for the deficit of heterozygotes repeatedly found in natural aspen populations.

Genomic and transcriptomic analyses reveal polygenic architecture for ecologically important traits in aspen (Populus tremuloides Michx.).

Intraspecific genetic variation in foundation species such as aspen (Populus tremuloides Michx.) shapes their impact on forest structure and function. Identifying genes underlying ecologically important traits is key to understanding that impact. Previous studies, using single-locus genome-wide association (GWA) analyses to identify candidate genes, have identified fewer genes than anticipated for highly heritable quantitative traits. Mounting evidence suggests that polygenic control of quantitative traits is largely responsible for this "missing heritability" phenomenon. Our research characterized the genetic architecture of 30 ecologically important traits using a common garden of aspen through genomic and transcriptomic analyses. A multilocus association model revealed that most traits displayed a highly polygenic architecture, with most variation explained by loci with small effects (likely below the detection levels of single-locus GWA methods). Consistent with a polygenic architecture, our single-locus GWA analyses found only 38 significant SNPs in 22 genes across 15 traits. Next, we used differential expression analysis on a subset of aspen genets with divergent concentrations of salicinoid phenolic glycosides (key defense traits). This complementary method to traditional GWA discovered 1243 differentially expressed genes for a polygenic trait. Soft clustering analysis revealed three gene clusters (241 candidate genes) involved in secondary metabolite biosynthesis and regulation. Our work reveals that ecologically important traits governing higher-order community- and ecosystem-level attributes of a foundation forest tree species have complex underlying genetic structures and will require methods beyond traditional GWA analyses to unravel.

Aspen phenylpropanoid genes' expression levels correlate with genets' tannin richness and vary both in responses to soil nitrogen and associations with phenolic profiles.

Condensed tannin (CT) contents of European aspen (Populus tremula L.) vary among genotypes, and increases in nitrogen (N) availability generally reduce plants' tannin production in favor of growth, through poorly understood mechanisms. We hypothesized that intrinsic tannin production rates may co-vary with gene expression responses to soil N and resource allocation within the phenylpropanoid pathway (PPP). Thus, we examined correlations between soil N levels and both expression patterns of eight PPP genes (measured by quantitative-reverse transcription PCR) and foliar phenolic compounds (measured by liquid chromatography-mass spectrometry) in young aspen genets with intrinsically extreme CT levels. Monitored phenolics included salicinoids, lignins, flavones, flavonols, CT precursors and CTs. The PPP genes were consistently expressed more strongly in high-CT trees. Low N supplements reduced expression of genes throughout the PPP in all genets, while high N doses restored expression of genes at the beginning and end of the pathway. These PPP changes were not reflected in pools of tannin precursors, but varying correlations between gene expression and foliar phenolic pools were detected in young and mature leaves, suggesting that processes linking gene expression and the resulting phenolics vary spatially and temporally. Precursor fluxes suggested that CT-related metabolic rate or sink controls are linked to intrinsic carbon allocation strategies associated with N responses. Overall, we found more negative correlations (indicative of allocation trade-offs) between PPP gene expression and phenolic products following N additions in low-CT plants than in high-CT plants. The tannin-related expression dynamics suggest that, in addition to defense, relative tannin levels may also be indicative of intraspecific variations in the way aspen genets respond to soil fertility.

Allelic and population variation of microsatellite loci in aspen (Populus tremuloides).

To develop a robust basis for inferences about population genetics and evolution, this work assayed 192 aspens (Populus tremuloides) from 11 sites in Wisconsin, USA, for allelic and population variation at 16 microsatellite loci distributed across the Populus genome. Frequency distributions of fluorochrome-labeled alleles resolved by capillary electrophoresis were analyzed for relationships to repeat size and number. Population-level statistics were compared with those of other studies, especially in Populus. All loci were polymorphic, varying widely in the number of alleles per locus (mean = 8.25, range 2-20). Expected and observed heterozygosities were high (0.45 and 0.41, respectively), with little differentiation among populations (F(ST) = 0.006-0.045) and a moderate level of inbreeding (F(IS) = 0.09), intermediate among levels reported in studies based on isozymes. Contrary to several other reports, allele frequencies clustered tightly around the modal frequency, and the genetic diversity (measured as alleles per locus or as expected heterozygosity) was not related to either the repeat unit size or to the number of repeats.