Nucleotide polymorphism and phenotypic associations within and around the phytochrome B2 Locus in European aspen (Populus tremula, Salicaceae)

We investigated the utility of association mapping to dissect the genetic basis of naturally occurring variation in bud phenology in European aspen (Populus tremula). With this aim, we surveyed nucleotide polymorphism in 13 fragments spanning an 80-kb region surrounding the phytochrome B2 (phyB2) locus. Although polymorphism varies substantially across the phyB2 region, we detected no signs for deviations from neutral expectations. We also identified a total of 41 single nucleotide polymorphisms (SNPs) that were subsequently scored in a mapping population consisting of 120 trees. We identified two nonsynonymous SNPs in the phytochrome B2 gene that were independently associated with variation in the timing of bud set and that explained between 1.5 and 5% of the observed phenotypic variation in bud set. Earlier studies have shown that the frequencies of both these SNPs vary clinally with latitude. Linkage disequilibrium across the region was low, suggesting that the SNPs we identified are strong candidates for being causally linked to variation in bud set in our mapping populations. One of the SNPs (T608N) is located in the "hinge region," close to the chromophore binding site of the phyB2 protein. The other SNP (L1078P) is located in a region supposed to mediate downstream signaling from the phyB2 locus. The lack of population structure, combined with low levels of linkage disequilibrium, suggests that association mapping is a fruitful method for dissecting naturally occurring variation in Populus tremula.     

Adaptive population differentiation in phenology across a latitudinal gradient in European aspen (Populus tremula, L.): a comparison of neutral markers, candidate genes and phenotypic traits

A correct timing of growth cessation and dormancy induction represents a critical ecological and evolutionary trade-off between survival and growth in most forest trees (Rehfeldt et al. 1999; Horvath et al. 2003; Howe et al. 2003). We have studied the deciduous tree European Aspen (Populus tremula) across a latitudinal gradient and compared genetic differentiation in phenology traits with molecular markers. Trees from 12 different areas covering 10 latitudinal degrees were cloned and planted in two common gardens. Several phenology traits showed strong genetic differentiation and clinal variation across the latitudinal gradient, with Q(ST) values generally exceeding 0.5. This is in stark contrast to genetic differentiation at several classes of genetic markers (18 neutral SSRs, 7 SSRs located close to phenology candidate genes and 50 SNPs from five phenology candidate genes) that all showed F(ST) values around 0.015. We thus find strong evidence for adaptive divergence in phenology traits across the latitudinal gradient. However, the strong population structure seen at the quantitative traits is not reflected in underlying candidate genes. This result fit theoretical expectations that suggest that genetic differentiation at candidate loci is better described by F(ST) at neutral loci rather than by Q(ST) at the quantitative traits themselves.     

Disentangling the roles of history and local selection in shaping clinal variation of allele frequencies and gene expression in Norway spruce (Picea abies)

Understanding the genetic basis of local adaptation is challenging due to the subtle balance among conflicting evolutionary forces that are involved in its establishment and maintenance. One system with which to tease apart these difficulties is clines in adaptive characters. Here we analyzed genetic and phenotypic variation in bud set, a highly heritable and adaptive trait, among 18 populations of Norway spruce (Picea abies), arrayed along a latitudinal gradient ranging from 47°N to 68°N. We confirmed that variation in bud set is strongly clinal, using a subset of five populations. Genotypes for 137 single-nucleotide polymorphisms (SNPs) chosen from 18 candidate genes putatively affecting bud set and 308 control SNPs chosen from 264 random genes were analyzed for patterns of genetic structure and correlation to environment. Population genetic structure was low (F(ST) = 0.05), but latitudinal patterns were apparent among Scandinavian populations. Hence, part of the observed clinal variation should be attributable to population demography. Conditional on patterns of genetic structure, there was enrichment of SNPs within candidate genes for correlations with latitude. Twenty-nine SNPs were also outliers with respect to F(ST). The enrichment for clinal variation at SNPs within candidate genes (i.e., SNPs in PaGI, PaPhyP, PaPhyN, PaPRR7, and PaFTL2) indicated that local selection in the 18 populations, and/or selection in the ancestral populations from which they were recently derived, shaped the observed cline. Validation of these genes using expression studies also revealed that PaFTL2 expression is significantly associated with latitude, thereby confirming the central role played by this gene in the control of phenology in plants.     

Bud set in poplar--genetic dissection of a complex trait in natural and hybrid populations

? The seasonal timing of growth events is crucial to tree distribution and conservation. The seasonal growth cycle is strongly adapted to the local climate that is changing because of global warming. We studied bud set as one cornerstone of the seasonal growth cycle in an integrative approach. ? Bud set was dissected at the phenotypic level into several components, and phenotypic components with most genetic variation were identified. While phenotypic variation resided in the timing of growth cessation, and even so more in the duration from growth cessation to bud set, the timing of growth cessation had a stronger genetic component in both natural and hybrid populations. ? Quantitative trait loci (QTL) were identified for the most discriminative phenotypic bud-set components across four poplar pedigrees. The QTL from different pedigrees were recurrently detected in six regions of the poplar genome. ? These regions of 1.83-4.25 Mbp in size, containing between 202 and 394 genes, form the basis for further molecular-genetic dissection of bud set.     

The Control of Autumn Senescence in European Aspen

The initiation, progression, and natural variation of autumn senescence in European aspen (Populus tremula) was investigated by monitoring chlorophyll degradation in (1) trees growing in natural stands and (2) cloned trees growing in a greenhouse under various light regimes. The main trigger for the initiation of autumn senescence in aspen is the shortening photoperiod, but there was a large degree of variation in the onset of senescence, both within local populations and among trees originating from different populations, where it correlated with the latitude of their respective origins. The variation for onset of senescence with a population was much larger than the variation of bud set. Once started, autumn senescence was accelerated by low temperature and longer nights, and clones that started to senescence late had a faster senescence. Bud set and autumn senescence appeared to be under the control of two independent critical photoperiods, but senescence could not be initiated until a certain time after bud set, suggesting that bud set and growth arrest are important for the trees to acquire competence to respond to the photoperiodic trigger to undergo autumn senescence. A timetable of events related to bud set and autumn senescence is presented.     

The adaptive potential of Populus balsamifera L. to phenology requirements in a warmer global climate

The manner in which organisms adapt to climate change informs a broader understanding of the evolution of biodiversity as well as conservation and mitigation plans. We apply common garden and association mapping approaches to quantify genetic variance and identify loci affecting bud flush and bud set, traits that define a tree's season for height growth, in the boreal forest tree Populus balsamifera L. (balsam poplar). Using data from 478 genotypes grown in each of two common gardens, one near the southern edge and another near the northern edge of P. balsamifera's range, we found that broad‐sense heritability for bud flush and bud set was generally high (H² > 0.5 in most cases), suggesting that abundant genetic variation exists for phenological response to changes in the length of the growing season. To identify the molecular genetic basis of this variation, we genotyped trees for 346 candidate single nucleotide polymorphisms (SNPs) from 27 candidate genes for the CO/FT pathway in poplar. Mixed‐model analyses of variance identified SNPs in 10 genes to be associated with variation in either bud flush or bud set. Multiple SNPs within FRIGIDA were associated with bud flush, whereas multiple SNPs in LEAFY and GIGANTEA 5 were associated with bud set. Although there was strong population structure in stem phenology, the geographic distribution of multilocus association SNP genotypes was widespread except at the most northern populations, indicating that geographic regions may harbour sufficient diversity in functional genes to facilitate adaption to future climatic conditions in many sites.     

Imprints of Natural Selection Along Environmental Gradients in Phenology-Related Genes of Quercus petraea

We explored single nucleotide polymorphism (SNP) variation in candidate genes for bud burst from Quercus petraea populations sampled along gradients of latitude and altitude in Western Europe. SNP diversity was monitored for 106 candidate genes, in 758 individuals from 32 natural populations. We investigated whether SNP variation reflected the clinal pattern of bud burst observed in common garden experiments. We used different methods to detect imprints of natural selection (F_ST outlier, clinal variation at allelic frequencies, association tests) and compared the results obtained for the two gradients. F_ST outlier SNPs were found in 15 genes, 5 of which were common to both gradients. The type of selection differed between the two gradients (directional or balancing) for 3 of these 5. Clinal variations were observed for six SNPs, and one cline was conserved across both gradients. Association tests between the phenotypic or breeding values of trees and SNP genotypes identified 14 significant associations, involving 12 genes. The results of outlier detection on the basis of population differentiation or clinal variation were not very consistent with the results of association tests. The discrepancies between these approaches may reflect the different hierarchical levels of selection considered (inter- and intrapopulation selection). Finally, we obtained evidence for convergent selection (similar for gradients) and clinal variation for a few genes, suggesting that comparisons between parallel gradients could be used to screen for major candidate genes responding to natural selection in trees.     

Natural phenological variation in aspen (Populus tremula): the SwAsp collection

The genus Populus is currently the main model system for genetic, genomic, and physiological research in trees. Phenotypic variation in aspen (Populus tremula) populations growing in different environments across Sweden is expected to reflect genetic variation that is important for local adaptation. To analyze such natural phenotypic and genetic variation, the Swedish Aspen (SwAsp) Collection was established. Trees were taken from 12 different populations across Sweden, from 56° to 66° latitude north and planted in two common gardens in Ekebo (55.9°N) and Sävar (63.4°N). Data related to phenological and growth traits were collected during the second year of growth. Some traits like the date of bud set and leaf area duration showed strong clinal variation patterns with latitude in both field trials, but the date of bud flush did not change along a latitudinal cline. The phenological traits showed moderate within-populations heritabilities, although growth traits showed weaker clinal patterns and lower heritabilities than the phenological traits. This research forms the starting point for the development of the SwAsp collection, a resource facilitating analysis of the natural genetic variation in aspen, the elucidation of the structure and dynamics of aspen populations, and the future identification of the genes controlling adaptive traits using association mapping of selected candidate genes.     

Low night temperature and inhibition of gibberellin biosynthesis override phytochrome action and induce bud set and cold acclimation, but not dormancy in PHYA overexpressors and wild-type of hybrid aspen

Juvenile trees of temperate and boreal regions cease growth and set buds in autumn in response to short day-lengths (SD) detected by phytochrome. Growth cessation and bud set are prerequisites for the development of winter dormancy and full cold hardiness. In this study we show that the SD-requirement for bud set and cold hardening can be overcome in hybrid aspen (Populus tremula L. × tremuloides Michx.) by low night temperature and inhibition of gibberellin (GA) biosynthesis. Bud set and increased cold hardiness were observed under normally non-inductive long day-length (LD) in wild-type plants, when exposed to low night temperature and paclobutrazol. In addition, the effect of PHYA overexpression could be overcome in transgenic plants, producing bud set and cold acclimation by treatment with: SD, low night temperature and paclobutrazol. After cold acclimation, the degree of bud dormancy was lower for wild-type plants prior treated with LD and transgenic plants (overexpressing PHYA), than SD-treated, wild-type plants. Thus, low night temperature in combination with reduced GA content induced bud set and promoted cold hardiness under normally non-inductive photoperiods in hybrid aspen, but was unable to affect development of dormancy. This might suggest separate signalling pathways from phytochrome regulating the induction of cold/cold hardiness and bud dormancy in hybrid aspen or alternatively, there might be one pathway that fails to complete its action in the transgenic and paclobutrazol treated plants.     

Admixture facilitates adaptation from standing variation in the European aspen (Populus tremula L.), a widespread forest tree

Adaptation to new environments can start from new mutations or from standing variation already present in natural populations. Whether admixture constrains or facilitates adaptation from standing variation is largely unknown, especially in ecological keystone or foundation species. We examined patterns of neutral and adaptive population divergence in Populus tremula L., a widespread forest tree, using mapped molecular genetic markers. We detected the genetic signature of postglacial admixture between a Western and an Eastern lineage of P. tremula in Scandinavia, an area suspected to represent a zone of postglacial contact for many species of animals and plants. Stringent divergence‐based neutrality tests provided clear indications for locally varying selection at the European scale. Six of 12 polymorphisms under selection were located less than 1 kb away from the nearest gene predicted by the Populus trichocarpa genome sequence. Few of these loci exhibited a signature of ‘selective sweeps’ in diversity‐based tests, which is to be expected if adaptation occurs primarily from standing variation. In Scandinavia, admixture explained genomic patterns of ancestry and the nature of clinal variation and strength of selection for bud set, a phenological trait of great adaptive significance in temperate trees, measured in a common garden trial. Our data provide a hitherto missing direct link between past range shifts because of climatic oscillations, and levels of standing variation currently available for selection and adaptation in a terrestrial foundation species.     

Timing of photoperiodic competency causes phenological mismatch in balsam poplar (Populus balsamifera L.)

Plant phenology is expected to be sensitive to climate warming. In boreal trees, spring flush is primarily temperature driven, whereas height growth cessation and autumn leaf senescence are predominantly controlled by photoperiod. Cuttings of 525 genotypes from the full range of balsam poplar were planted into two common gardens (Vancouver and Indian Head, Canada) at similar latitudes, but with differing winter temperatures and growing seasons. There was clinal variation in spring and, particularly, summer and fall phenology. Bud flush and, despite milder climate, bud set and leaf drop were earlier at Vancouver than at Indian Head by 44, 28 and 7 d, respectively. Although newly flushed growth is insensitive to photoperiod, many genotypes at both sites became competent before the summer solstice. At Vancouver, high‐latitude genotypes set dormant terminal buds in mid‐spring. Most other genotypes grew until midsummer or set bud temporarily and then experienced a second flush. In both gardens and in a growth chamber experiment, earlier bud set was associated with reduced height growth and higher root/shoot ratios. Shoots attained competency ～5 weeks after flushing, which would normally prevent dormancy induction before the solstice, but may be insufficient if spring advances by more than a few weeks.     

Photoperiodic responses of a northern and southern ecotype of black cottonwood

Photoperiod is an important signal controlling the onset of dormancy in perennial plants. Short days typically induce growth cessation, the initiation of cold acclimation, the formation of a terminal bud. bud dormancy and other adaptive responses. Photoperiodic ecotypes have evolved in many species with large latitudinal distributions. The photoperiodic responses of two northern (53°35′ and 53°50′N) and two southern (34°10′ and 40°32′N) genotypes of black cottonwood (Populus trichocarpa Torr. & Gray) were characterized by growing trees under a range of photoperiods in the greenhouse and growth chamber. Short days induced bud set in both ecotypes. resulting in trees with fewer leaves and less height growth than trees grown under long days. Short days also enhanced anthocyanin accumulation in the northern ecotype and decreased branching of the southernmost genotype. Two aspects of the photoperiodic response were evaluated for each trail: critical photoperiod. which was defined as the longest photoperiod that elicited a short-day response, and photoperiodic sensitivity, which was defined as the change in response per unit change in photoperiod. For each of the traits analyzed, the northern ecotype had a longer critical photoperiod and greater photoperiodic sensitivity than did the southern ecotype. The short critical photoperiod and reduced photoperiodic sensitivity of the southern ecotype resulted in a significant delay in bud set compared to that of the northern ecotype, even under a 9-h photoperiod. Typically, photoperiodic ecotypes have been characterized as having different critical photoperiods. Ecotypic differences in photoperiodic sensitivity, however, indicate that differences in the photoperiodic response curves cannot be completely described by the critical photoperiod alone. These results also suggest that the critical photoperiod. photoperiodic sensitivity and speed of bud set have a common physiological basis. Bud set occurred earlier hi the northern ecotype primarily because bud scale leaves were initiated earlier. For one of the northern genotypes, leaf primordia that were initialed under long days subsequently differentiated into bud scale leaves after the trees were transferred to a 9-h photoperiod. This demonstrates that primordia initiated under long days are not necessarily committed to becoming foliage leaves. The response to photoperiod did not differ appreciably between the greenhouse and growth chamber conditions that were tested.     

Phytochrome control of short-day-induced bud set in black cottonwood

In trees and other woody perennial plants, short days (SDs) typically induce growth cessation, the initiation of cold acclimation, the formation of a terminal bud and bud dormancy. Phytochrome control of SD-induced bud set was investigated in two northern clones of black cottonwood (Populus trichocarpa Torr. & Gray) by using night breaks with red light (R) and far-red light (FR). For both clones (BC-1 and BC-2), SD-induced bud set was prevented when R night breaks as short as 2 min were given in the middle of the night. When night breaks with 2 min of R were immediately followed by 2 min of FR, substantial reversibility of bud set was observed for BC-1 but not for BC-2. By comparing the effects of the R night breaks on bud set and the length of specific internodes, we determined that the R night breaks influenced internode elongation in two opposing ways. First, the addition of a R night break to the SD treatment prevented the cessation of internode elongation that is associated with bud set. Those internodes that would not have elongated under SDs (and would have been found within the terminal bud) elongated in the R treatment. Second, the R night breaks decreased internode length relative to the long-day (LD) control. In contrast to the clonal differences in reversibility that we observed for bud set, the decrease in internode length (i.e. the second effect of R) was R/FR reversible in both clones. Based on these results, we conclude that internode elongation is influenced by two distinct types of phytochrome-mediated response. The first response is a typical response to photpperiod, whereas the second response is a typical “end-of-day” response to light quality. Our results demonstrate that SD-induced bud set in black cottonwood is controlled by phytochrome but that clonal differences have an important influence on the R/FR reversibility of this response. The availability of an experimental system in which SD-induced bud set is R/FR reversible will be valuable for studying the physiological genetics of photoperiodism in trees.     

Clinal Variation at Phenology-Related Genes in Spruce: Parallel Evolution in FTL2 and Gigantea?

Parallel clines in different species, or in different geographical regions of the same species, are an important source of information on the genetic basis of local adaptation. We recently detected latitudinal clines in SNPs frequencies and gene expression of candidate genes for growth cessation in Scandinavian populations of Norway spruce (Picea abies). Here we test whether the same clines are also present in Siberian spruce (P. obovata), a close relative of Norway spruce with a different Quaternary history. We sequenced nine candidate genes and 27 control loci and genotyped 14 SSR loci in six populations of P. obovata located along the Yenisei river from latitude 56°N to latitude 67°N. In contrast to Scandinavian Norway spruce that both departs from the standard neutral model (SNM) and shows a clear population structure, Siberian spruce populations along the Yenisei do not depart from the SNM and are genetically unstructured. Nonetheless, as in Norway spruce, growth cessation is significantly clinal. Polymorphisms in photoperiodic (FTL2) and circadian clock (Gigantea, GI, PRR3) genes also show significant clinal variation and/or evidence of local selection. In GI, one of the variants is the same as in Norway spruce. Finally, a strong cline in gene expression is observed for FTL2, but not for GI. These results, together with recent physiological studies, confirm the key role played by FTL2 and circadian clock genes in the control of growth cessation in spruce species and suggest the presence of parallel adaptation in these two species.     

Geographic structure in metabolome and herbivore community co‐occurs with genetic structure in plant defence genes

Plant–herbivore interactions vary across the landscape and have been hypothesised to promote local adaption in plants to the prevailing herbivore regime. Herbivores that feed on European aspen (Populus tremula) change across regional scales and selection on host defence genes may thus change at comparable scales. We have previously observed strong population differentiation in a set of inducible defence genes in Swedish P. tremula. Here, we study the geographic patterns of abundance and diversity of herbivorous insects, the untargeted metabolome of the foliage and genetic variation in a set of wound‐induced genes and show that the geographic structure co‐occurs in all three data sets. In response to this structure, we observe local maladaptation of herbivores, with fewer herbivores on local trees than on trees originated from more distant localities. Finally, we also identify 28 significant associations between single nucleotide polymorphisms SNPs from defence genes and a number of the herbivore traits and metabolic profiles.     

Effects of red, far-red and blue light in maintaining growth in latitudinal populations of Norway spruce (Picea abies)

Seedlings of trees with a free growth pattern cease growth when night-lengths become shorter than a critical value, and this critical night-length (CNL) decreases with increasing latitude of origin. In northern populations, the light quality also appears to play an important role and a clinal variation in requirement for far-red (FR) light has been documented. In this study we dissected the light quality requirements for maintaining growth in different latitudinal populations of Norway spruce (Picea abies (L.) H. Karst.) using light emitting diodes for red (R), FR and blue (B) light, as 12 h day extension to provide 24 h photoperiod. At equal spectral photon flux, FR light was more effective than R light in maintaining growth, and the requirement of both R and FR increased with northern latitude of origin. One-to-one mixtures of R and FR light were more effective in maintaining growth than either FR or R light alone, indicating a possible interaction between R and FR light maintaining growth. Using the blue light as day extension could not prevent growth cessation in any of the populations, but delayed the bud set slightly in all populations. Our results suggest that phytochrome(s) are the primary photoreceptors in high irradiance responses maintaining growth in Norway spruce seedlings.     

Secondary contact and local adaptation contribute to genome‐wide patterns of clinal variation in Drosophila melanogaster

Populations arrayed along broad latitudinal gradients often show patterns of clinal variation in phenotype and genotype. Such population differentiation can be generated and maintained by both historical demographic events and local adaptation. These evolutionary forces are not mutually exclusive and can in some cases produce nearly identical patterns of genetic differentiation among populations. Here, we investigate the evolutionary forces that generated and maintain clinal variation genome‐wide among populations of Drosophila melanogaster sampled in North America and Australia. We contrast patterns of clinal variation in these continents with patterns of differentiation among ancestral European and African populations. Using established and novel methods we derive here, we show that recently derived North America and Australia populations were likely founded by both European and African lineages and that this hybridization event likely contributed to genome‐wide patterns of parallel clinal variation between continents. The pervasive effects of admixture mean that differentiation at only several hundred loci can be attributed to the operation of spatially varying selection using an F_ST outlier approach. Our results provide novel insight into the well‐studied system of clinal differentiation in D. melanogaster and provide a context for future studies seeking to identify loci contributing to local adaptation in a wide variety of organisms, including other invasive species as well as temperate endemics.