Quantitative trait loci and candidate gene mapping of bud set and bud flush in populus

Induced and natural variation affect traits independently in hybrid Populus

The genetic control of many plant traits can be highly complex. Both allelic variation (sequence change) and dosage variation (copy number change) contribute to a plant's phenotype. While numerous studies have investigated the effect of allelic or dosage variation, very few have documented both within the same system, leaving their relative contribution to phenotypic effects unclear. The Populus genome is highly polymorphic, and poplars are fairly tolerant of gene dosage variation. Here, using a previously established Populus hybrid F1 population, we assessed and compared the effect of natural allelic variation and induced dosage variation on biomass, phenology, and leaf morphology traits. We identified QTLs for many of these traits, but our results indicate limited overlap between the QTLs associated with natural allelic variation and induced dosage variation. Additionally, the integration of data from both allelic and dosage variation identifies a larger set of QTLs that together explain a larger percentage of the phenotypic variance. Finally, our results suggest that the effect of the large indels might mask that of allelic QTLs. Our study helps clarify the relationship between allelic and dosage variation and their effects on quantitative traits.

ELONGATED HYPOCOTYL 5a modulates FLOWERING LOCUS T2 and gibberellin levels to control dormancy and bud break in poplar

Photoperiod is a crucial environmental cue for phenological responses, including growth cessation and winter dormancy in perennial woody plants. Two regulatory modules within the photoperiod pathway explain bud dormancy induction in poplar (Populus spp.): the circadian oscillator LATE ELONGATED HYPOCOTYL 2 (LHY2) and GIGANTEA-like genes (GIs) both regulate the key target for winter dormancy induction FLOWERING LOCUS T2 (FT2). However, modification of LHY2 and GIs cannot completely prevent growth cessation and bud set under short-day (SD) conditions, indicating that additional regulatory modules are likely involved. We identified PtoHY5a, an orthologs of the photomorphogenesis regulatory factor ELONGATED HYPOCOTYL 5 (HY5) in poplar (Populus tomentosa), that directly activates PtoFT2 expression and represses the circadian oscillation of LHY2, indirectly activating PtoFT2 expression. Thus, PtoHY5a suppresses SD-induced growth cessation and bud set. Accordingly, PtoHY5a knockout facilitates dormancy induction. PtoHY5a also inhibits bud-break in poplar by controlling gibberellic acid (GA) levels in apical buds. Additionally, PtoHY5a regulates the photoperiodic control of seasonal growth downstream of phytochrome PHYB2. Thus, PtoHY5a modulates seasonal growth in poplar by regulating the PtoPHYB2-PtoHY5a-PtoFT2 module to determine the onset of winter dormancy, and by fine-tuning GA levels to control bud-break.

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.

Genetic association of toxin production in the dinoflagellate Alexandrium minutum

Dinoflagellates of the genus Alexandrium are responsible for harmful algal blooms and produce paralytic shellfish toxins (PSTs). Their very large and complex genomes make it challenging to identify the genes responsible for toxin synthesis. A family-based genomic association study was developed to determine the inheritance of toxin production in Alexandrium minutum and identify genomic regions linked to this production. We show that the ability to produce toxins is inheritable in a Mendelian way, while the heritability of the toxin profile is more complex. We developed the first dinoflagellate genetic linkage map. Using this map, several major results were obtained: 1. A genomic region related to the ability to produce toxins was identified. 2. This region does not contain any polymorphic sxt genes, known to be involved in toxin production in cyanobacteria. 3. The sxt genes, known to be present in a single cluster in cyanobacteria, are scattered on different linkage groups in A. minutum. 4. The expression of two sxt genes not assigned to any linkage group, sxtI and sxtG, may be regulated by the genomic region related to the ability to produce toxins. Our results provide new insights into the organization of toxicity-related genes in A. minutum, suggesting a dissociated genetic mechanism for the production of the different analogues and the ability to produce toxins. However, most of the newly identified genes remain unannotated. This study therefore proposes new candidate genes to be further explored to understand how dinoflagellates synthesize their toxins.

Evidence of climate-driven selection on tree traits and trait plasticity across the climatic range of a riparian foundation species

Selection on quantitative traits by heterogeneous climatic conditions can lead to substantial trait variation across a species range. In the context of rapidly changing environments, however, it is equally important to understand selection on trait plasticity. To evaluate the role of selection in driving divergences in traits and their associated plasticities within a widespread species, we compared molecular and quantitative trait variation in Populus fremontii (Fremont cottonwood), a foundation riparian distributed throughout Arizona. Using SNP data and genotypes from 16 populations reciprocally planted in three common gardens, we first performed Q_ST -F_ST analyses to detect selection on traits and trait plasticity. We then explored the environmental drivers of selection using trait-climate and plasticity-climate regressions. Three major findings emerged: 1) There was significant genetic variation in traits expressed in each of the common gardens and in the phenotypic plasticity of traits across gardens, both of which were heritable. 2) Based on Q_ST -F_ST comparisons, there was evidence of selection in all traits measured; however, this result varied from no effect in one garden to highly significant in another, indicating that detection of past selection is environmentally dependent. We also found strong evidence of divergent selection on plasticity across environments for two traits. 3) Traits and/or their plasticity were often correlated with population source climate (R² up to 0.77 and 0.66, respectively). These results suggest that steep climate gradients across the Southwest have played a major role in shaping the evolution of divergent phenotypic responses in populations and genotypes now experiencing climate change.

Photoperiod Control of Plant Growth: Flowering Time Genes Beyond Flowering

Fluctuations in environmental conditions greatly influence life on earth. Plants, as sessile organisms, have developed molecular mechanisms to adapt their development to changes in daylength, or photoperiod. One of the first plant features that comes to mind as affected by the duration of the day is flowering time; we all bring up a clear image of spring blossom. However, for many plants flowering happens at other times of the year, and many other developmental aspects are also affected by changes in daylength, which range from hypocotyl elongation in Arabidopsis thaliana to tuberization in potato or autumn growth cessation in trees. Strikingly, many of the processes affected by photoperiod employ similar gene networks to respond to changes in the length of light/dark cycles. In this review, we have focused on developmental processes affected by photoperiod that share similar genes and gene regulatory networks.

Phytochrome B and PHYTOCHROME INTERACTING FACTOR8 modulate seasonal growth in trees

The seasonally synchronized annual growth cycle that is regulated mainly by photoperiod and temperature cues is a crucial adaptive strategy for perennial plants in boreal and temperate ecosystems. Phytochrome B (phyB), as a light and thermal sensor, has been extensively studied in Arabidopsis. However, the specific mechanisms for how the phytochrome photoreceptors control the phenology in tree species remain poorly understood. We characterized the functions of PHYB genes and their downstream PHYTOCHROME INTERACTING FACTOR (PIF) targets in the regulation of shade avoidance and seasonal growth in hybrid aspen trees. We show that while phyB1 and phyB2, as phyB in other plants, act as suppressors of shoot elongation during vegetative growth, they act as promoters of tree seasonal growth. Furthermore, while the Populus homologs of both PIF4 and PIF8 are involved in the shade avoidance syndrome (SAS), only PIF8 plays a major role as a suppressor of seasonal growth. Our data suggest that the PHYB-PIF8 regulon controls seasonal growth through the regulation of FT and CENL1 expression while a genome-wide transcriptome analysis suggests how, in Populus trees, phyB coordinately regulates SAS responses and seasonal growth cessation.

Phenotypic and Genomic Local Adaptation across Latitude and Altitude in Populus trichocarpa

Local adaptation to climate allows plants to cope with temporally and spatially heterogeneous environments, and parallel phenotypic clines provide a natural experiment to uncover the genomic architecture of adaptation. Though extensive effort has been made to investigate the genomic basis of local adaptation to climate across the latitudinal range of tree species, less is known for altitudinal clines. We used exome capture to genotype 451 Populus trichocarpa genotypes across altitudinal and latitudinal gradients spanning the natural species range, and phenotyped these trees for a variety of adaptive traits in two common gardens. We observed clinal variation in phenotypic traits across the two transects, which indicates climate-driven selection, and coupled gene-based genotype–phenotype and genotype–environment association scans to identify imprints of climatic adaptation on the genome. Although many of the phenotype- and climate-associated genes were unique to one transect, we found evidence of parallelism between latitude and altitude, as well as significant convergence when we compared our outlier genes with those putatively involved in climatic adaptation in two gymnosperm species. These results suggest that not only genomic constraint during adaptation to similar environmental gradients in poplar but also different environmental contexts, spatial scale, and perhaps redundant function among potentially adaptive genes and polymorphisms lead to divergent adaptive architectures.

A high-density exome capture genotype-by-sequencing panel for forestry breeding in Pinus radiata

Development of genome-wide resources for application in genomic selection or genome-wide association studies, in the absence of full reference genomes, present a challenge to the forestry industry, where longer breeding cycles could benefit from the accelerated selection possible through marker-based breeding value predictions. In particular, large conifer megagenomes require a strategy to reduce complexity, whilst ensuring genome-wide coverage is achieved. Using a transcriptome-based reference template, we have successfully developed a high density exome capture genotype-by-sequencing panel for radiata pine (Pinus radiata D.Don), capable of capturing in excess of 80,000 single nucleotide polymorphism (SNP) markers with a minor allele frequency above 0.03 in the population tested. This represents approximately 29,000 gene models from a core set of 48,914 probes. A set of 704 SNP markers capable of pedigree reconstruction and differentiating individual genotypes were tested within two full-sib mapping populations. While as few as 70 markers could reconstruct parentage in almost all cases, the impact of missing genotypes was noticeable in several offspring. Therefore, 60 sets of 110 randomly selected SNP markers were compared for both parentage reconstruction and clone differentiation. The performance in parentage reconstruction showed little variation over 60 iterations. However, there was notable variation in discriminatory power between closely related individuals, indicating a higher density SNP marker panel may be required to elucidate hidden relationships in complex pedigrees.

Genome-Wide Association Study to Identify Favorable SNP Allelic Variations and Candidate Genes That Control the Timing of Spring Bud Flush of Tea (Camellia sinensis) Using SLAF-seq

The timing of spring bud flush (TBF) is of economic importance for tea plant (Camellia sinensis) breeding. We employed a genome-wide association study (GWAS) to identify favorable single nucleotide polymorphism (SNP) allelic variations as well as candidate genes that control TBF of C. sinensis using specific-locus-amplified fragment sequencing (SLAF-seq) in a diversity panel comprising 151 tea plant germplasm resources. GWAS analysis revealed 26 SNPs associated with TBF in three years, and we eventually identified a final significant SNP for TBF. To identify candidate genes possibly related to TBF, we screened seven candidate genes within 100 kb regions surrounding the trait-related SNP loci. Furthermore, the favorable allelic variation, the "TT" genotype in the SNP loci, was discovered, and a derived cleaved amplified polymorphism (dCAPS) marker was designed that cosegregated with TBF, which could be used for marker-assisted selection (MAS) breeding in C. sinensis. The results obtained from this study can provide a theoretical and applied basis for the MAS of early breeding in tea plants in the future.

A perspective on ecologically relevant plant-UV research and its practical application

Plants perceive ultraviolet-B (UV-B) radiation through the UV-B photoreceptor UV RESISTANCE LOCUS 8 (UVR8), and initiate regulatory responses via associated signalling networks, gene expression and metabolic pathways. Various regulatory adaptations to UV-B radiation enable plants to harvest information about fluctuations in UV-B irradiance and spectral composition in natural environments, and to defend themselves against UV-B exposure. Given that UVR8 is present across plant organs and tissues, knowledge of the systemic signalling involved in its activation and function throughout the plant is important for understanding the context of specific responses. Fine-scale understanding of both UV-B irradiance and perception within tissues and cells requires improved application of knowledge about UV-attenuation in leaves and canopies, warranting greater consideration when designing experiments. In this context, reciprocal crosstalk among photoreceptor-induced pathways also needs to be considered, as this appears to produce particularly complex patterns of physiological and morphological response. Through crosstalk, plant responses to UV-B radiation go beyond simply UV-protection or amelioration of damage, but may give cross-protection over a suite of environmental stressors. Overall, there is emerging knowledge showing how information captured by UVR8 is used to regulate molecular and physiological processes, although understanding of upscaling to higher levels of organisation, i.e. organisms, canopies and communities remains poor. Achieving this will require further studies using model plant species beyond Arabidopsis, and that represent a broad range of functional types. More attention should also be given to plants in natural environments in all their complexity, as such studies are needed to acquire an improved understanding of the impact of climate change in the context of plant-UV responses. Furthermore, broadening the scope of experiments into the regulation of plant-UV responses will facilitate the application of UV radiation in commercial plant production. By considering the progress made in plant-UV research, this perspective highlights prescient topics in plant-UV photobiology where future research efforts can profitably be focussed. This perspective also emphasises burgeoning interdisciplinary links that will assist in understanding of UV-B effects across organisational scales and gaps in knowledge that need to be filled so as to achieve an integrated vision of plant responses to UV-radiation.

A comprehensive genomic scan reveals gene dosage balance impacts on quantitative traits in Populus trees

Gene dosage variation and the associated changes in gene expression influence a wide variety of traits, ranging from cancer in humans to yield in plants. It is also expected to affect important traits of ecological and agronomic importance in forest trees, but this variation has not been systematically characterized or exploited. Here we performed a comprehensive scan of the Populus genome for dosage-sensitive loci affecting quantitative trait variation for spring and fall phenology and biomass production. The study population was a large collection of clonally propagated F1 hybrid lines of Populus that saturate the genome 10-fold with deletions and insertions (indels) of known sizes and positions. As a group, the phenotypic means of the indel lines consistently differed from control nonindel lines, with an overall negative effect of both insertions and deletions on all biomass-related traits but more diverse effects and an overall wider phenotypic distribution of the indel lines for the phenology-related traits. We also investigated the correlation between gene dosage at specific chromosomal locations and phenotype, to identify dosage quantitative trait loci (dQTL). Such dQTL were detected for most phenotypes examined, but stronger effect dQTL were identified for the phenology-related traits than for the biomass traits. Our genome-wide screen for dosage sensitivity in a higher eukaryote demonstrates the importance of global genomic balance and the impact of dosage on life history traits.

The influence of spectral composition on spring and autumn phenology in trees

Several recent reviews highlight the molecular mechanisms that underpin phenological responses to temperature and photoperiod; however, these have mostly overlooked the influence of solar radiation and its spectral composition on these processes. For instance, solar radiation in the blue and ultraviolet (UV) regions of the spectrum, as well as the red/far-red (R:FR) ratio, can influence spring and autumn phenology. Solar radiation reaching the Earth changes diurnally and seasonally; however, rising global temperatures, latitudinal range shifts and light pollution are likely to produce novel combinations of phenological cues for tree species. Here, we review the literature on phenological responses to spectral composition. Our objective was to explore the natural variation in spectral composition using radiative transfer models and to reveal any species-specific or ecotype-specific responses relating to latitudinal origin. These responses are likely to be most pronounced at high latitudes where spectral composition varies most throughout the year. For instance, trees from high latitudes tend to be more sensitive to changes in R:FR than those from low latitudes. The effects of blue light and UV radiation on phenology have not been studied as much as those of R:FR, but the limited results available suggest both could be candidate cues affecting autumn leaf colouration and senescence. Failure of more-southern species and ecotypes to adapt and use spectral cues during northwards range shifts could result in mistimed phenology, potentially resulting in frost damage, reduced fitness and limited range expansion. Future areas for research should look to establish how consistently different functional types of tree respond to spectral cues and identify photoreceptor-mediated mechanisms that allow plants to combine information from multiple light cues to coordinate the timing of phenological events. It should then be feasible to consider the synchronous or sequential action of light cues within a hierarchy of environmental factors regulating phenology.

Genotypic variation in phenological plasticity: Reciprocal common gardens reveal adaptive responses to warmer springs but not to fall frost

Species faced with rapidly shifting environments must be able to move, adapt, or acclimate in order to survive. One mechanism to meet this challenge is phenotypic plasticity: altering phenotype in response to environmental change. Here, we investigated the magnitude, direction, and consequences of changes in two key phenology traits (fall bud set and spring bud flush) in a widespread riparian tree species, Populus fremontii. Using replicated genotypes from 16 populations from throughout the species' thermal range, and reciprocal common gardens at hot, warm, and cool sites, we identified four major findings: (a) There are significant genetic (G), environmental (E), and GxE components of variation for both traits across three common gardens; (b) The magnitude of phenotypic plasticity is correlated with provenance climate, where trees from hotter, southern populations exhibited up to four times greater plasticity compared to the northern, frost-adapted populations; (c) Phenological mismatches are correlated with higher mortality as the transfer distances between provenance and garden increase; and (d) The relationship between plasticity and survival depends not only on the magnitude and direction of environmental transfer, but also on the type of environmental stress (i.e., heat or freezing), and how particular traits have evolved in response to that stress. Trees transferred to warmer climates generally showed small to moderate shifts in an adaptive direction, a hopeful result for climate change. Trees experiencing cooler climates exhibited large, non-adaptive changes, suggesting smaller transfer distances for assisted migration. This study is especially important as it deconstructs trait responses to environmental cues that are rapidly changing (e.g., temperature and spring onset) and those that are fixed (photoperiod), and that vary across the species' range. Understanding the magnitude and adaptive nature of phenotypic plasticity of multiple traits responding to multiple environmental cues is key to guiding restoration management decisions as climate continues to change.

Construction of a high-density genetic map and its application for leaf shape QTL mapping in poplar

High-quality and dense genetic maps were constructed, and leaf shape variation was dissected by QTL mapping in poplar. Species in the genus Populus, also known as poplars, are important woody species and considered model plants for perennial trees. High-density genetic maps are valuable genomic resources for population genetics. Here, we generated a high-quality and dense genetic map for an F1 poplar population using high-throughput NGS-based genotyping. A total of 92,097 high-quality SNP markers were developed by stringent filtering and identification. In total, 889 and 1650 SNPs formed the female and male genetic maps, respectively. To test the application of the genetic maps, QTL mapping of leaf shape was conducted for this F1 population. A total of nine parameters were scored for leaf shape variation in three different environments. Combining genetic maps and measurements of the nine leaf shape parameters, we mapped a total of 42 significant QTLs. The highest LOD score of all QTLs was 9.2, and that QTL explained the most (15.13%) trait variation. A total of nine QTLs could be detected in at least two environments, and they were located in two genomic regions. Within these two QTL regions, some candidate genes for regulating leaf shape were predicted through functional annotation. The successful mapping of leaf shape QTLs demonstrated the utility of our genetic maps. According to the performance of this study, we were able to provide high-quality and dense genetic maps and dissect the leaf shape variation in poplar.

Natural variation of DNA methylation and gene expression may determine local adaptations of Scots pine populations

Long-lived conifers are vulnerable to climate change because classical evolutionary processes are slow in developing adaptive responses. Therefore, the capacity of a genotype to adopt different phenotypes is important. Gene expression is the primary mechanism that converts genome-encoded information into phenotypes, and DNA methylation is employed in the epigenetic regulation of gene expression. We investigated variations in global DNA methylation and gene expression between three Scots pine (Pinus sylvestris L.) populations located in northern and southern Finland using mature seeds. Gene expression levels were studied in six DNA methyltransferase (DNMT) genes, which were characterized in this study, and in 19 circadian clock genes regulating adaptive traits. In embryos, expression diversity was found for three DNMT genes, which maintain DNA methylation. The expression of two DNMT genes was strongly correlated with climate variables, which suggests a role for DNA methylation in local adaptation. For adaptation-related genes, expression levels showed between-population variation in 11 genes in megagametophytes and in eight genes in embryos, and many of these genes were linked to climate factors. Altogether, our results suggest that differential DNA methylation and gene expression contribute to local adaptation in Scots pine populations and may enhance the fitness of trees under rapidly changing climatic conditions.

Mechanism of Overwintering in Trees

Boreal trees possess very high freezing resistance, which is induced by short-day length and low temperatures, in order to survive severe subzero temperatures in winter. During autumn, cooperation of photoreceptors and circadian clock system perceiving photoperiod shortening results in growth cessation, dormancy development, and first induction of freezing resistance. The freezing resistance is further enhanced by subsequent low temperature during seasonal cold acclimation with concomitant changes in various morphological and physiological features including accumulation of sugars and late embryogenesis abundant proteins. The mechanism of adaptation to freezing temperatures differs depending on the type of tissue in boreal trees. For example, bark, cambium, and leaf cells tolerate freezing-induced dehydration by extracellular freezing, whereas xylem parenchyma cells avoid intracellular freezing by deep supercooling. In addition, dormant buds in some trees respond by extraorgan freezing. Boreal trees have evolved overwintering mechanisms such as dormancy and high freezing resistance in order to survive freezing temperatures in winter.

Ecological genomics of variation in bud-break phenology and mechanisms of response to climate warming in Populus trichocarpa

Spring bud-break phenology is a critical adaptive feature common to temperate perennial woody plants. Understanding the molecular underpinnings of variation in bud-break is important for elucidating adaptive evolution and predicting outcomes relating to climate change. Field and controlled growth chamber tests were used to assess population-wide patterns in bud-break from wild-sourced black cottonwood (Populus trichocarpa) genotypes. We conducted a genome-wide association study (GWAS) with single nucleotide polymorphisms (SNPs) derived from whole genome sequencing to test for loci underlying variation in bud-break. Bud-break had a quadratic relationship with latitude, where southern- and northern-most provenances generally broke bud earlier than those from central parts of the species' range. Reduced winter chilling increased population-wide variation in bud-break, whereas greater chilling decreased variation. GWAS uncovered 16 loci associated with bud-break. Phenotypic changes connected with allelic variation were replicated in an independent set of P. trichocarpa trees. Despite phenotypic similarities, genetic profiles between southern- and northern-most genotypes were dissimilar based on our GWAS-identified SNPs. We propose that the GWAS-identified loci underpin the geographical pattern in P. trichocarpa and that variation in bud-break reflects different selection for winter chilling and heat sum accumulation, both of which can be affected by climate warming.

Photoperiodic Regulation of Shoot Apical Growth in Poplar

Woody perennials adapt their genetic traits to local climate conditions. Day length plays an essential role in the seasonal growth of poplar trees. When photoperiod falls below a given critical day length, poplars undergo growth cessation and bud set. A leaf-localized mechanism of photoperiod measurement triggers the transcriptional modulation of a long distance signaling molecule, FLOWERING LOCUS T (FT). This molecule targets meristem function giving rise to these seasonal responses. Studies over the past decade have identified conserved orthologous genes involved in photoperiodic flowering in Arabidopsis that regulate poplar vegetative growth. However, phenological and molecular examination of key photoperiod signaling molecules reveals functional differences between these two plant model systems suggesting alternative components and/or regulatory mechanisms operating during poplar vegetative growth. Here, we review current knowledge and provide new data regarding the molecular components of the photoperiod measuring mechanism that regulates annual growth in poplar focusing on main achievements and new perspectives.

Development of F1 hybrid population and the high-density linkage map for European aspen (Populus tremula L.) using RADseq technology

BACKGROUND

Restriction-site associated DNA sequencing (RADseq) technology was recently employed to identify a large number of single nucleotide polymorphisms (SNP) for linkage mapping of a North American and Eastern Asian Populus species. However, there is also the need for high-density genetic linkage maps for the European aspen (P. tremula) as a tool for further mapping of quantitative trait loci (QTLs) and marker-assisted selection of the Populus species native to Europe.

RESULTS

We established a hybrid F1 population from the cross of two aspen parental genotypes diverged in their phenological and morphological traits. We performed RADseq of 122 F1 progenies and two parents yielding 15,732 high-quality SNPs that were successfully identified using the reference genome of P. trichocarpa. 2055 SNPs were employed for the construction of maternal and paternal linkage maps. The maternal linkage map was assembled with 1000 SNPs, containing 19 linkage groups and spanning 3054.9 cM of the genome, with an average distance of 3.05 cM between adjacent markers. The paternal map consisted of 1055 SNPs and the same number of linkage groups with a total length of 3090.56 cM and average interval distance of 2.93 cM. The linkage maps were employed for QTL mapping of one-year-old seedlings height variation. The most significant QTL (LOD = 5.73) was localized to LG5 (96.94 cM) of the male linkage map, explaining 18% of the phenotypic variation.

CONCLUSIONS

The set of 15,732 SNPs polymorphic in aspen and high-density genetic linkage maps constructed for the P. tremula intra-specific cross will provide a valuable source for QTL mapping and identification of candidate genes facilitating marker-assisted selection in European aspen.

The shifting phenological landscape: Within- and between-species variation in leaf emergence in a mixed-deciduous woodland

Many organisms rely on synchronizing the timing of their life-history events with those of other trophic levels-known as phenological matching-for survival or successful reproduction. In temperate deciduous forests, the extent of matching with the budburst date of key tree species is of particular relevance for many herbivorous insects and, in turn, insectivorous birds. In order to understand the ecological and evolutionary forces operating in these systems, we require knowledge of the factors influencing leaf emergence of tree communities. However, little is known about how phenology at the level of individual trees varies across landscapes, or how consistent this spatial variation is between different tree species. Here, we use field observations, collected over 2 years, to characterize within- and between-species differences in spring phenology for 825 trees of six species (Quercus robur, Fraxinus excelsior, Fagus sylvatica, Betula pendula, Corylus avellana, and Acer pseudoplatanus) in a 385-ha woodland. We explore environmental predictors of individual variation in budburst date and bud development rate and establish how these phenological traits vary over space. Trees of all species showed markedly consistent individual differences in their budburst timing. Bud development rate also varied considerably between individuals and was repeatable in oak, beech, and sycamore. We identified multiple predictors of budburst date including altitude, local temperature, and soil type, but none were universal across species. Furthermore, we found no evidence for interspecific covariance of phenology over space within the woodland. These analyses suggest that phenological landscapes are highly complex, varying over small spatial scales both within and between species. Such spatial variation in vegetation phenology is likely to influence patterns of selection on phenology within populations of consumers. Knowledge of the factors shaping the phenological environments experienced by animals is therefore likely to be key in understanding how these evolutionary processes operate.

Bud phenology and growth are subject to divergent selection across a latitudinal gradient in Populus angustifolia and impact adaptation across the distributional range and associated arthropods

Temperate forest tree species that span large geographical areas and climatic gradients often have high levels of genetic variation. Such species are ideal for testing how neutral demographic factors and climate‐driven selection structure genetic variation within species, and how this genetic variation can affect ecological communities. Here, we quantified genetic variation in vegetative phenology and growth traits in narrowleaf cottonwood, Populus angustifolia, using three common gardens planted with genotypes originating from source populations spanning the species' range along the Rocky Mountains of North America (ca. 1700 km). We present three main findings. First, we found strong evidence of divergent selection (Q_ST > F_ST) on fall phenology (bud set) with adaptive consequences for frost avoidance. We also found evidence for selection on bud flush duration, tree height, and basal diameter, resulting in population differentiation. Second, we found strong associations with climate variables that were strongly correlated with latitude of origin. More strongly differentiated traits also showed stronger climate correlations, which emphasizes the role that climate has played in divergent selection throughout the range. We found population × garden interaction effects; for some traits, this accounted for more of the variance than either factor alone. Tree height was influenced by the difference in climate of the source and garden locations and declined with increasing transfer distance. Third, growth traits were correlated with dependent arthropod community diversity metrics. Synthesis. Overall, we conclude that climate has influenced genetic variation and structure in phenology and growth traits and leads to local adaptation in P. angustifolia, which can then impact dependent arthropod species. Importantly, relocation of genotypes far northward or southward often resulted in poor growth, likely due to a phenological mismatch with photoperiod, the proximate cue for fall growth cessation. Genotypes moved too far southward suffer from early growth cessation, whereas those moved too far northward are prone to fall frost and winter dieback. In the face of current and forecasted climate change, habitat restoration, forestry, and tree breeding efforts should utilize these findings to better match latitudinal and climatic source environments with management locations for optimal future outcomes.

Construction of High-Density Linkage Maps of Populus deltoides × P. simonii Using Restriction-Site Associated DNA Sequencing

Although numerous linkage maps have been constructed in the genus Populus, they are typically sparse and thus have limited applications due to low throughput of traditional molecular markers. Restriction-site associated DNA sequencing (RADSeq) technology allows us to identify a large number of single nucleotide polymorphisms (SNP) across genomes of many individuals in a fast and cost-effective way, and makes it possible to construct high-density genetic linkage maps. We performed RADSeq for 299 progeny and their two parents in an F1 hybrid population generated by crossing the female Populus deltoides 'I-69' and male Populus simonii 'L3'. A total of 2,545 high quality SNP markers were obtained and two parent-specific linkage maps were constructed. The female genetic map contained 1601 SNPs and 20 linkage groups, spanning 4,249.12 cM of the genome with an average distance of 2.69 cM between adjacent markers, while the male map consisted of 940 SNPs and also 20 linkage groups with a total length of 3,816.24 cM and an average marker interval distance of 4.15 cM. Finally, our analysis revealed that synteny and collinearity are highly conserved between the parental linkage maps and the reference genome of P. trichocarpa. We demonstrated that RAD sequencing is a powerful technique capable of rapidly generating a large number of SNPs for constructing genetic maps in outbred forest trees. The high-quality linkage maps constructed here provided reliable genetic resources to facilitate locating quantitative trait loci (QTLs) that control growth and wood quality traits in the hybrid population.

Molecular regulation of phenology in trees-because the seasons they are a-changin'

The perennial trees, in contrast to the much more studied annual plants, have to adapt their vegetative growth and development to the sometimes extremely contrasting environmental conditions that occur over the different seasons. Recently, studies of the molecular framework underlying this adaptation in Populus trees is reinforcing the notion that the genetic pathways controlling growth and dormancy cycles have a remarkable conservation with the pathways controlling the regulation of flowering time in annual plants. Insight into these mechanisms will be important for our understanding of how trees will respond to various future global climate scenarios.

EARLY BUD-BREAK1 (EBB1) defines a conserved mechanism for control of bud-break in woody perennials

Bud-break is an environmentally and economically important trait in trees, shrubs and vines from temperate latitudes. Poor synchronization of bud-break timing with local climates can lead to frost injuries, susceptibility to pests and pathogens and poor crop yields in fruit trees and vines. The rapid climate changes outpace the adaptive capacities of plants to respond through natural selection. This is particularly true for trees which have long generation cycle and thus the adaptive changes are significantly delayed. Therefore, to devise appropriate breeding and conservation strategies, it is imperative to understand the molecular underpinnings that govern dormancy mechanisms. We have recently identified and characterized the poplar EARLY BUD-BREAK 1 (EBB1) gene. EBB1 is a positive regulator of bud-break and encodes a transcription factor from the AP2/ERF family. Here, using comparative and functional genomics approaches we show that EBB1 function in regulation of bud-break is likely conserved across wide range of woody perennial species with importance to forestry and agriculture.

Biomass traits and candidate genes for bioenergy revealed through association genetics in coppiced European Populus nigra (L.)

BACKGROUND

Second generation (2G) bioenergy from lignocellulosic feedstocks has the potential to develop as a sustainable source of renewable energy; however, significant hurdles still remain for large-scale commercialisation. Populus is considered as a promising 2G feedstock and understanding the genetic basis of biomass yield and feedstock quality are a research priority in this model tree species.

RESULTS

We report the first coppiced biomass study for 714 members of a wide population of European black poplar (Populus nigra L.), a native European tree, selected from 20 river populations ranging in latitude and longitude between 40.5 and 52.1°N and 1.0 and 16.4°E, respectively. When grown at a single site in southern UK, significant Site of Origin (SO) effects were seen for 14 of the 15 directly measured or derived traits including biomass yield, leaf area and stomatal index. There was significant correlation (p < 0.001) between biomass yield traits over 3 years of harvest which identified leaf size and cell production as strong predictors of biomass yield. A 12 K Illumina genotyping array (constructed from 10,331 SNPs in 14 QTL regions and 4648 genes) highlighted significant population genetic structure with pairwise FST showing strong differentiation (p < 0.001) between the Spanish and Italian subpopulations. Robust associations reaching genome-wide significance are reported for main stem height and cell number per leaf; two traits tightly linked to biomass yield. These genotyping and phenotypic data were also used to show the presence of significant isolation by distance (IBD) and isolation by adaption (IBA) within this population.

CONCLUSIONS

The three associations identified reaching genome-wide significance at p < 0.05 include a transcription factor; a putative stress response gene and a gene of unknown function. None of them have been previously linked to bioenergy yield; were shown to be differentially expressed in a panel of three selected genotypes from the collection and represent exciting, novel candidates for further study in a bioenergy tree native to Europe and Euro-Asia. A further 26 markers (22 genes) were found to reach putative significance and are also of interest for biomass yield, leaf area, epidermal cell expansion and stomatal patterning. This research on European P. nigra provides an important foundation for the development of commercial native trees for bioenergy and for advanced, molecular breeding in these species.

From genotypes to phenotypes: expression levels of genes encompassing adaptive SNPs in black spruce

Measuring transcript levels for adaptive genes revealed polymorphisms having cis -effect upon gene expression levels related to phenotype variation in a black spruce natural population. Trees growing in temperate and boreal regions must acclimate to changes in climatic factors such as low winter temperatures to survive to seasonal variations. Common garden studies have shown that genetic variation in quantitative traits helps species to survive and adapt to environmental changes and local conditions. Twenty-four genes carrying SNPs were previously associated with genetic adaptation in black spruce (Picea mariana [Mill.] BSP). The objectives of this study were to investigate the potential role of these genes in regulation of winter acclimation and adaptation by studying their patterns of expression as a function of the physiological stage during the annual growth cycle, tissue type, and their SNP genotypic class. Considerable variability in gene expression was observed between different vegetative tissues or organs, and between physiological stages. The genes were expressed predominantly in tissues that could be linked more directly to winter acclimation and adaptation. The expression levels of several of the genes were significantly related to variation in tree height growth or budset timing and expression level variation related to SNP genotypic classes was observed in four of the genes. An interaction between genotypic classes and physiological stages was also observed for some genes, indicating genotypes with different reaction norms in terms of gene expression.

Extensive Transcriptome Changes During Natural Onset and Release of Vegetative Bud Dormancy in Populus

To survive winter, many perennial plants become endodormant, a state of suspended growth maintained even in favorable growing environments. To understand vegetative bud endodormancy, we collected paradormant, endodormant, and ecodormant axillary buds from Populus trees growing under natural conditions. Of 44,441 Populus gene models analyzed using NimbleGen microarrays, we found that 1,362 (3.1%) were differentially expressed among the three dormancy states, and 429 (1.0%) were differentially expressed during only one of the two dormancy transitions (FDR p-value < 0.05). Of all differentially expressed genes, 69% were down-regulated from paradormancy to endodormancy, which was expected given the lower metabolic activity associated with endodormancy. Dormancy transitions were accompanied by changes in genes associated with DNA methylation (via RNA-directed DNA methylation) and histone modifications (via Polycomb Repressive Complex 2), confirming and extending knowledge of chromatin modifications as major features of dormancy transitions. Among the chromatin-associated genes, two genes similar to SPT (SUPPRESSOR OF TY) were strongly up-regulated during endodormancy. Transcription factor genes and gene sets that were atypically up-regulated during endodormancy include a gene that seems to encode a trihelix transcription factor and genes associated with proteins involved in responses to ethylene, cold, and other abiotic stresses. These latter transcription factors include ETHYLENE INSENSITIVE 3 (EIN3), ETHYLENE-RESPONSIVE ELEMENT BINDING PROTEIN (EBP), ETHYLENE RESPONSE FACTOR (ERF), ZINC FINGER PROTEIN 10 (ZAT10), ZAT12, and WRKY DNA-binding domain proteins. Analyses of phytohormone-associated genes suggest important changes in responses to ethylene, auxin, and brassinosteroids occur during endodormancy. We found weaker evidence for changes in genes associated with salicylic acid and jasmonic acid, and little evidence for important changes in genes associated with gibberellins, abscisic acid, and cytokinin. We identified 315 upstream sequence motifs associated with eight patterns of gene expression, including novel motifs and motifs associated with the circadian clock and responses to photoperiod, cold, dehydration, and ABA. Analogies between flowering and endodormancy suggest important roles for genes similar to SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL), DORMANCY ASSOCIATED MADS-BOX (DAM), and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS 1 (SOC1).

EARLY BUD-BREAK 1 (EBB1) is a regulator of release from seasonal dormancy in poplar trees

Trees from temperate latitudes transition between growth and dormancy to survive dehydration and freezing stress during winter months. We used activation tagging to isolate a dominant mutation affecting release from dormancy and identified the corresponding gene EARLY BUD-BREAK 1 (EBB1). We demonstrate through positioning of the tag, expression analysis, and retransformation experiments that EBB1 encodes a putative APETALA2/Ethylene responsive factor transcription factor. Transgenic up-regulation of the gene caused early bud-flush, whereas down-regulation delayed bud-break. Native EBB1 expression was highest in actively growing apices, undetectable during the dormancy period, but rapidly increased before bud-break. The EBB1 transcript was localized in the L1/L2 layers of the shoot meristem and leaf primordia. EBB1-overexpressing transgenic plants displayed enlarged shoot meristems, open and poorly differentiated buds, and a higher rate of cell division in the apex. Transcriptome analyses of the EBB1 transgenics identified 971 differentially expressed genes whose expression correlated with the EBB1 expression changes in the transgenic plants. Promoter analysis among the differentially expressed genes for the presence of a canonical EBB1-binding site identified 65 putative target genes, indicative of a broad regulatory context of EBB1 function. Our results suggest that EBB1 has a major and integrative role in reactivation of meristem activity after winter dormancy.

The perennial clock is an essential timer for seasonal growth events and cold hardiness

Over the last several decades, changes in global temperatures have led to changes in local environments affecting the growth conditions for many species. This is a trend that makes it even more important to understand how plants respond to local variations and seasonal changes in climate. To detect daily and seasonal changes as well as acute stress factors such as cold and drought, plants rely on a circadian clock. This chapter introduces the current knowledge and literature about the setup and function of the circadian clock in various tree and perennial species, with a focus on the Populus genus.

The embryonic shoot: a lifeline through winter

The tiny vascular axis of the embryo emerges post-embryonically as an elaborate and critical infrastructure, pervading the entire plant system. Its expansive nature is especially impressive in trees, where growth and development continue for extended periods. While the shoot apical meristem (SAM) orchestrates primary morphogenesis, the vascular system is mapped out in its wake in the provascular cylinder, situated just below the emerging leaf primordia and surrounding the rib meristem. Formation of leaf primordia and provascular tissues is incompatible with the harsh conditions of winter. Deciduous trees of boreal and temperate climates therefore enter a survival mode at the end of the season. However, to be competitive, they need to maximize their growth period while avoiding cellular frost damage. Trees achieve this by monitoring photoperiod, and by timely implementation of a survival strategy that schedules downstream events, including growth cessation, terminal bud formation, dormancy assumption, acquisition of freezing tolerance, and shedding of leaves. Of central importance are buds, which contain an embryonic shoot that allows shoot development and elongation in spring. The genetic and molecular processes that drive the cycle in synchrony with the seasons are largely elusive. Here, we review what is known about the signals and signal conduits that are involved, the processes that are initiated, and the developmental transitions that ensue in a terminal bud. We propose that addressing dormancy as a property of the SAM and the bud as a unique shoot type will facilitate our understanding of winter dormancy.

The sex-linked region in Populus tremuloides Turesson 141 corresponds to a pericentromeric region of about two million base pairs on P. trichocarpa chromosome 19

In the dioecious genus Populus, sex determination has been located to chromosome 19. However, despite a high degree of genome collinearity, various Populus species seem to differ with regard to the location of the sex-determining region on the respective chromosome and the apparent heterogametic sex. In this study, the boundaries of the recombination-suppressed, sex-linked region of the male P. tremuloides clone Turesson 141 were localised by genetic mapping using new SNP and InDel markers. The respective region seems to be located in a pericentromeric position. The corresponding P. trichocarpa genome region spans about two million bp and comprises 65 gene loci, which were bioinformatically evaluated for their potential as candidate genes for sex determination. Three putative transcription factor genes and four genes that are potentially involved in flower development processes, e.g. meristem transition from the vegetative to the reproductive phase, were identified. Populus tremuloides sequence data of the sex-linked region is required for a final search for candidate genes.

Quantitative trait loci affecting reproductive phenology in peach

BACKGROUND

The reproductive phenology of perennial plants in temperate climates is largely conditioned by the duration of bud dormancy, and fruit developmental processes. Bud dormancy release and bud break depends on the perception of cumulative chilling and heat during the bud development. The objective of this work was to identify new quantitative trait loci (QTLs) associated to temperature requirements for bud dormancy release and flowering and to fruit harvest date, in a segregating population of peach.

RESULTS

We have identified QTLs for nine traits related to bud dormancy, flowering and fruit harvest in an intraspecific hybrid population of peach in two locations differing in chilling time accumulation. QTLs were located in a genetic linkage map of peach based on single nucleotide polymorphism (SNP) markers for eight linkage groups (LGs) of the peach genome sequence. QTLs for chilling requirements for dormancy release and blooming clustered in seven different genomic regions that partially coincided with loci identified in previous works. The most significant QTL for chilling requirements mapped to LG1, close to the evergrowing locus. QTLs for heat requirement related traits were distributed in nine genomic regions, four of them co-localizing with QTLs for chilling requirement trait. Two major loci in LG4 and LG6 determined fruit harvest time.

CONCLUSIONS

We identified QTLs associated to nine traits related to the reproductive phenology in peach. A search of candidate genes for these QTLs rendered different genes related to flowering regulation, chromatin modification and hormone signalling. A better understanding of the genetic factors affecting crop phenology might help scientists and breeders to predict changes in genotype performance in a context of global climate change.

Genetic architecture of spring and autumn phenology in Salix

BACKGROUND

In woody plants from temperate regions, adaptation to the local climate results in annual cycles of growth and dormancy, and optimal regulation of these cycles are critical for growth, long-term survival, and competitive success. In this study we have investigated the genetic background to growth phenology in a Salix pedigree by assessing genetic and phenotypic variation in growth cessation, leaf senescence and bud burst in different years and environments. A previously constructed linkage map using the same pedigree and anchored to the annotated genome of P. trichocarpa was improved in target regions and used for QTL analysis of the traits. The major aims in this study were to map QTLs for phenology traits in Salix, and to identify candidate genes in QTL hot spots through comparative mapping with the closely related Populus trichocarpa.

RESULTS

All traits varied significantly among genotypes and the broad-sense heritabilities ranged between 0.5 and 0.9, with the highest for leaf senescence. In total across experiment and years, 80 QTLs were detected. For individual traits, the QTLs explained together from 21.5 to 56.5% of the variation. Generally each individual QTL explained a low amount of the variation but three QTLs explained above 15% of the variation with one QTL for leaf senescence explaining 34% of the variation. The majority of the QTLs were recurrently identified across traits, years and environments. Two hotspots were identified on linkage group (LG) II and X where narrow QTLs for all traits co-localized.

CONCLUSIONS

This study provides the most detailed analysis of QTL detection for phenology in Salix conducted so far. Several hotspot regions were found where QTLs for different traits and QTLs for the same trait but identified during different years co-localised. Many QTLs co-localised with QTLs found in poplar for similar traits that could indicate common pathways for these traits in Salicaceae. This study is an important first step in identifying QTLs and candidate genes for phenology traits in Salix.

A molecular framework for seasonal growth-dormancy regulation in perennial plants

The timing of the onset and release of dormancy impacts the survival, productivity and spatial distribution of temperate horticultural and forestry perennials and is mediated by at least three main regulatory programs involving signal perception and processing by phytochromes (PHYs) and PHY-interacting transcription factors (PIFs). PIF4 functions as a key regulator of plant growth in response to both external and internal signals. In poplar, the expression of PIF4 and PIF3-LIKE1 is upregulated in response to short days, while PHYA and PHYB are not regulated at the transcriptional level. Integration of light and environmental signals is achieved by gating the expression and transcriptional activity of PIF4. During this annual cycle, auxin promotes the degradation of Aux/IAA transcriptional repressors through the SKP-Cullin-F-boxTIR1 complex, relieving the repression of auxin-responsive genes by allowing auxin response factors (ARFs) to activate the transcription of auxin-responsive genes involved in growth responses. Analyses of transcriptome changes during dormancy transitions have identified MADS-box transcription factors associated with endodormancy induction. Previous studies show that poplar dormancy-associated MADS-box (DAM) genes PtMADS7 and PtMADS21 are differentially regulated during the growth-dormancy cycle. Endodormancy may be regulated by internal factors, which are specifically localized in buds. PtMADS7/PtMADS21 may function as an internal regulator in poplar. The control of flowering time shares certain regulatory hierarchies with control of the dormancy/growth cycle. However, the particularities of different stages of the dormancy/growth cycle warrant comprehensive approaches to identify the causative genes for the entire cycle. A growing body of knowledge also indicates epigenetic regulation plays a role in these processes in perennial horticultural and forestry plants. The increased knowledge contributes to better understanding of the dormancy process and consequently to precise manipulation of dormancy-related horticultural traits, such as flowering time.

Genome-wide expression profiles of seasonal bud dormancy at four critical stages in Japanese apricot

Dormancy is one of the most important adaptive mechanisms developed by perennial plants. To reveal the comprehensive mechanism of seasonal bud dormancy at four critical stages in Japanese apricot (Prunus persica), we applied Illumina sequencing to study differentially expressed genes (DEGs) at the transcriptional level. As a result, 19,759, 16,375, 19,749 and 20,800 tag-mapped genes were sequenced from libraries of paradormancy (R1), endodormancy (R2), ecodormancy (R3) and dormancy release (R4) stages based on the P. persica genome. Moreover, 6,199, 5,539, and 5,317 genes were differentially expressed in R1 versus R2, R2 versus R3, and R3 versus R4, respectively. Gene Ontology analysis of dormancy-related genes showed that these were mainly related to the cytoplasm, cytoplasmic part metabolism, intracellular metabolism and membrane-bound organelle metabolism. Pathway-enrichment annotation revealed that highly ranked genes were involved in ribosome pathways and protein processing in the endoplasmic reticulum. The results demonstrated that hormone response genes such as auxin, abscisic acid, ethylene and jasmonic acid, as well as zinc finger family protein genes are possibly involved in seasonal bud dormancy in Japanese apricot. The expression patterns of DEGs were verified using real-time quantitative RT-PCR. These results contribute to further understanding of the mechanism of bud dormancy in Japanese apricot.

The genomic architecture and association genetics of adaptive characters using a candidate SNP approach in boreal black spruce

BACKGROUND

The genomic architecture of adaptive traits remains poorly understood in non-model plants. Various approaches can be used to bridge this gap, including the mapping of quantitative trait loci (QTL) in pedigrees, and genetic association studies in non-structured populations. Here we present results on the genomic architecture of adaptive traits in black spruce, which is a widely distributed conifer of the North American boreal forest. As an alternative to the usual candidate gene approach, a candidate SNP approach was developed for association testing.

RESULTS

A genetic map containing 231 gene loci was used to identify QTL that were related to budset timing and to tree height assessed over multiple years and sites. Twenty-two unique genomic regions were identified, including 20 that were related to budset timing and 6 that were related to tree height. From results of outlier detection and bulk segregant analysis for adaptive traits using DNA pool sequencing of 434 genes, 52 candidate SNPs were identified and subsequently tested in genetic association studies for budset timing and tree height assessed over multiple years and sites. A total of 34 (65%) SNPs were significantly associated with budset timing, or tree height, or both. Although the percentages of explained variance (PVE) by individual SNPs were small, several significant SNPs were shared between sites and among years.

CONCLUSIONS

The sharing of genomic regions and significant SNPs between budset timing and tree height indicates pleiotropic effects. Significant QTLs and SNPs differed quite greatly among years, suggesting that different sets of genes for the same characters are involved at different stages in the tree's life history. The functional diversity of genes carrying significant SNPs and low observed PVE further indicated that a large number of polymorphisms are involved in adaptive genetic variation. Accordingly, for undomesticated species such as black spruce with natural populations of large effective size and low linkage disequilibrium, efficient marker systems that are predictive of adaptation should require the survey of large numbers of SNPs. Candidate SNP approaches like the one developed in the present study could contribute to reducing these numbers.

Initial characterization of shade avoidance response suggests functional diversity between Populus phytochrome B genes

Shade avoidance signaling involves perception of incident red/far-red (R/FR) light by phytochromes (PHYs) and modulation of downstream transcriptional networks. Although these responses are well studied in Arabidopsis, little is known about the role of PHYs and the transcriptional responses to shade in the woody perennial Populus. Tissue expression and subcellular localization of Populus PHYs was studied by quantitative real-time PCR (qRT-PCR) and protoplast transient assay. Transgenic lines with altered PHYB1 and/or PHYB2 expression were used in phenotypic assays and transcript profiling with qRT-PCR. RNA-Seq was used to identify transcriptional responses to enriched FR light. All three PHYs were differentially expressed among tissue types and PHYBs were targeted to the nucleus under white light. Populus PHYB1 rescued Arabidopsis phyB mutant phenotypes. Phenotypes of Populus transgenic lines and the expression of candidate shade response genes suggested that PHYB1 and PHYB2 have distinct yet overlapping functions. RNA-Seq analysis indicated that genes associated with cell wall modification and brassinosteroid signaling were induced under enriched FR light in Populus. This study is an initial attempt at deciphering the role of Populus PHYs by evaluating transcriptional reprogramming to enriched FR and demonstrates functional diversity and overlap of the Populus PHYB1 and PHYB2 in regulating shade responses.

The dynamic nature of bud dormancy in trees: environmental control and molecular mechanisms

In tree species native to temperate and boreal regions, the activity-dormancy cycle is an important adaptive trait both for survival and growth. We discuss recent research on mechanisms controlling the overlapping developmental processes that define the activity-dormancy cycle, including cessation of apical growth, bud development, induction, maintenance and release of dormancy, and bud burst. The cycle involves an extensive reconfiguration of metabolism. Environmental control of the activity-dormancy cycle is based on perception of photoperiodic and temperature signals, reflecting adaptation to prevailing climatic conditions. Several molecular actors for control of growth cessation have been identified, with the CO/FT regulatory network and circadian clock having important coordinating roles in control of growth and dormancy. Other candidate regulators of bud set, dormancy and bud burst have been identified, such as dormancy-associated MADS-box factors, but their exact roles remain to be discovered. Epigenetic mechanisms also appear to factor in control of the activity-dormancy cycle. Despite evidence for gibberellins as negative regulators in growth cessation, and ABA and ethylene in bud formation, understanding of the roles that plant growth regulators play in controlling the activity-dormancy cycle is still very fragmentary. Finally, some of the challenges for further research in bud dormancy are discussed.

Genetic variation in functional traits influences arthropod community composition in aspen (Populus tremula L.)

We conducted a study of natural variation in functional leaf traits and herbivory in 116 clones of European aspen, Populus tremula L., the Swedish Aspen (SwAsp) collection, originating from ten degrees of latitude across Sweden and grown in a common garden. In surveys of phytophagous arthropods over two years, we found the aspen canopy supports nearly 100 morphospecies. We identified significant broad-sense heritability of plant functional traits, basic plant defence chemistry, and arthropod community traits. The majority of arthropods were specialists, those coevolved with P. tremula to tolerate and even utilize leaf defence compounds. Arthropod abundance and richness were more closely related to plant growth rates than general chemical defences and relationships were identified between the arthropod community and stem growth, leaf and petiole morphology, anthocyanins, and condensed tannins. Heritable genetic variation in plant traits in young aspen was found to structure arthropod community; however no single trait drives the preferences of arthropod folivores among young aspen genotypes. The influence of natural variation in plant traits on the arthropod community indicates the importance of maintaining genetic variation in wild trees as keystone species for biodiversity. It further suggests that aspen can be a resource for the study of mechanisms of natural resistance to herbivores.

Epicormic buds in trees: a review of bud establishment, development and dormancy release

The formation of epicormic sprouts on the boles of trees is a phenomenon that has, until recently, been poorly understood. Renewed interest in the topic in the last two decades has led to significant advances in our knowledge of the subject, especially in regard to bud anatomy, morphology and ontogeny. There exists, however, no comprehensive synthesis of results from different disciplines across genera and geographical areas; this review seeks to fill that void and provide a comprehensive framework capable of guiding future research. A tree's potential for producing epicormic branches is dependent on the number of buds that are produced on a growing shoot, the development of those buds and associated meristems over time and the factors that promote sprout formation or bud death. Based on the descriptions of a limited number of researched species, we were able to describe four different developmental strategies for epicormics based on characteristics of meristem development. Control over epicormic bud dormancy is complex, but it is clear that the traditional view of auxin-mediated dormancy release is incomplete. Genetic control over epicormic development is yet to be empirically proven. Future research should focus on clarifying these physiological and genetic controls of epicormic bud development as well as developing more robust methods for tracking epicormics in ecological and silvicultural studies.

Phenotypic plasticity, QTL mapping and genomic characterization of bud set in black poplar

BACKGROUND

The genetic control of important adaptive traits, such as bud set, is still poorly understood in most forest trees species. Poplar is an ideal model tree to study bud set because of its indeterminate shoot growth. Thus, a full-sib family derived from an intraspecific cross of P. nigra with 162 clonally replicated progeny was used to assess the phenotypic plasticity and genetic variation of bud set in two sites of contrasting environmental conditions.

RESULTS

Six crucial phenological stages of bud set were scored. Night length appeared to be the most important signal triggering the onset of growth cessation. Nevertheless, the effect of other environmental factors, such as temperature, increased during the process. Moreover, a considerable role of genotype × environment (G × E) interaction was found in all phenological stages with the lowest temperature appearing to influence the sensitivity of the most plastic genotypes.Descriptors of growth cessation and bud onset explained the largest part of phenotypic variation of the entire process. Quantitative trait loci (QTL) for these traits were detected. For the four selected traits (the onset of growth cessation (date2.5), the transition from shoot to bud (date1.5), the duration of bud formation (subproc1) and bud maturation (subproc2)) eight and sixteen QTL were mapped on the maternal and paternal map, respectively. The identified QTL, each one characterized by small or modest effect, highlighted the complex nature of traits involved in bud set process. Comparison between map location of QTL and P. trichocarpa genome sequence allowed the identification of 13 gene models, 67 bud set-related expressional and six functional candidate genes (CGs). These CGs are functionally related to relevant biological processes, environmental sensing, signaling, and cell growth and development. Some strong QTL had no obvious CGs, and hold great promise to identify unknown genes that affect bud set.

CONCLUSIONS

This study provides a better understanding of the physiological and genetic dissection of bud set in poplar. The putative QTL identified will be tested for associations in P. nigra natural populations. The identified QTL and CGs will also serve as useful targets for poplar breeding.

Local selection across a latitudinal gradient shapes nucleotide diversity in balsam poplar, Populus balsamifera L

Molecular studies of adaptive evolution often focus on detecting selective sweeps driven by positive selection on a species-wide scale; however, much adaptation is local, particularly of ecologically important traits. Here, we look for evidence of range-wide and local adaptation at candidate genes for adaptive phenology in balsam poplar, Populus balsamifera, a widespread forest tree whose range extends across environmental gradients of photoperiod and growing season length. We examined nucleotide diversity of 27 poplar homologs of the flowering-time network-a group of genes that control plant developmental phenology through interactions with environmental cues such as photoperiod and temperature. Only one gene, ZTL2, showed evidence of reduced diversity and an excess of fixed replacement sites, consistent with a species-wide selective sweep. Two other genes, LFY and FRI, harbored high levels of nucleotide diversity and exhibited elevated differentiation between northern and southern accessions, suggesting local adaptation along a latitudinal gradient. Interestingly, FRI has also been identified as a target of local selection between northern and southern accessions of Arabidopsis thaliana, indicating that this gene may be commonly involved in ecological adaptation in distantly related species. Our findings suggest an important role for local selection shaping molecular diversity and reveal limitations of inferring molecular adaptation from analyses designed only to detect species-wide selective sweeps.