Early evolution of the land plant circadian clock.

While angiosperm clocks can be described as an intricate network of interlocked transcriptional feedback loops, clocks of green algae have been modelled as a loop of only two genes. To investigate the transition from a simple clock in algae to a complex one in angiosperms, we performed an inventory of circadian clock genes in bryophytes and charophytes. Additionally, we performed functional characterization of putative core clock genes in the liverwort Marchantia polymorpha and the hornwort Anthoceros agrestis. Phylogenetic construction was combined with studies of spatiotemporal expression patterns and analysis of M. polymorpha clock gene mutants. Homologues to core clock genes identified in Arabidopsis were found not only in bryophytes but also in charophytes, albeit in fewer copies. Circadian rhythms were detected for most identified genes in M. polymorpha and A. agrestis, and mutant analysis supports a role for putative clock genes in M. polymorpha. Our data are in line with a recent hypothesis that adaptation to terrestrial life occurred earlier than previously expected in the evolutionary history of charophyte algae. Both gene duplication and acquisition of new genes was important in the evolution of the plant circadian clock, but gene loss has also contributed to shaping the clock of bryophytes.

No time for spruce: rapid dampening of circadian rhythms in Picea abies (L. Karst).

The identification and cloning of full-length homologs of circadian clock genes from Picea abies represent a first step to study the function and evolution of the circadian clock in gymnosperms. Phylogenetic analyses suggest that the sequences of key circadian clock genes are conserved between angiosperms and gymnosperms. though fewer homologous copies were found for most gene families in P. abies. We detected diurnal cycling of circadian clock genes in P. abies using quantitative real-time PCR; however, cycling appeared to be rapidly dampened under free-running conditions. Given the unexpected absence of transcriptional cycling during constant conditions, we employed a complementary method to assay circadian rhythmic outputs and measured delayed fluorescence in seedlings of Norway spruce. Neither of the two approaches to study circadian rhythms in Norway spruce could detect robust ∼24 h cycling behavior under constant conditions. These data suggest gene conservation but fundamental differences in clock function between gymnosperms and other plant taxa.

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 significant fraction of 21-nucleotide small RNA originates from phased degradation of resistance genes in several perennial species.

Small RNAs (sRNAs), including microRNA (miRNA) and short-interfering RNA (siRNA), are important in the regulation of diverse biological processes. Comparative studies of sRNAs from plants have mainly focused on miRNA, even though they constitute a mere fraction of the total sRNA diversity. In this study, we report results from an in-depth analysis of the sRNA population from the conifer spruce (Picea abies) and compared the results with those of a range of plant species. The vast majority of sRNA sequences in spruce can be assigned to 21-nucleotide-long siRNA sequences, of which a large fraction originate from the degradation of transcribed sequences related to nucleotide-binding site-leucine-rich repeat-type resistance genes. Over 90% of all genes predicted to contain either a Toll/interleukin-1 receptor or nucleotide-binding site domain showed evidence of siRNA degradation. The data further suggest that this phased degradation of resistance-related genes is initiated from miRNA-guided cleavage, often by an abundant 22-nucleotide miRNA. Comparative analysis over a range of plant species revealed a huge variation in the abundance of this phenomenon. The process seemed to be virtually absent in several species, including Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa), and nonvascular plants, while particularly high frequencies were observed in spruce, grape (Vitis vinifera), and poplar (Populus trichocarpa). This divergent pattern might reflect a mechanism to limit runaway transcription of these genes in species with rapidly expanding nucleotide-binding site-leucine-rich repeat gene families. Alternatively, it might reflect variation in a counter-counter defense mechanism between plant species.

FLOWERING LOCUS T/TERMINAL FLOWER1-like genes affect growth rhythm and bud set in Norway spruce.

The timing of bud set, as one determinant of the annual growth rhythm, is critical for local adaptation of the conifer Norway spruce (Picea abies). Previous gene expression and population genetic studies have suggested a role for P. abies FLOWERING LOCUS T/TERMINAL FLOWER1-Like2 (PaFTL2) in the control of growth cessation and bud set in Norway spruce as well as in local adaptation resulting in clinal variation for timing of bud set. Using transgenic plants with PaFTL2 driven by an inducible promoter, we found that PaFTL2 indeed induces bud set and most probably also growth cessation. PaFTL2 shows high expression around the procambium and vascular tissue and in the crown region in buds of both seedlings and older trees. Furthermore, PaFTL2 expression is induced in vegetative shoots and all bud types in late summer, when growth cessation occurs. This supports the notion that PaFTL2 is involved in growth cessation. A close paralog to PaFTL2, PaFTL1, is strongly expressed in meristems during the summer, possibly to repress meristem activity and the formation of needle primordia during this period. The temporal and spatial expression of PaFTL1 and PaFTL2 largely complement each other, which suggests that they act in concert to control perennial growth in Norway spruce.

Sequencing of the needle transcriptome from Norway spruce (Picea abies Karst L.) reveals lower substitution rates, but similar selective constraints in gymnosperms and angiosperms.

BACKGROUND: A detailed knowledge about spatial and temporal gene expression is important for understanding both the function of genes and their evolution. For the vast majority of species, transcriptomes are still largely uncharacterized and even in those where substantial information is available it is often in the form of partially sequenced transcriptomes. With the development of next generation sequencing, a single experiment can now simultaneously identify the transcribed part of a species genome and estimate levels of gene expression. RESULTS: mRNA from actively growing needles of Norway spruce (Picea abies) was sequenced using next generation sequencing technology. In total, close to 70 million fragments with a length of 76 bp were sequenced resulting in 5 Gbp of raw data. A de novo assembly of these reads, together with publicly available expressed sequence tag (EST) data from Norway spruce, was used to create a reference transcriptome. Of the 38,419 PUTs (putative unique transcripts) longer than 150 bp in this reference assembly, 83.5% show similarity to ESTs from other spruce species and of the remaining PUTs, 3,704 show similarity to protein sequences from other plant species, leaving 4,167 PUTs with limited similarity to currently available plant proteins. By predicting coding frames and comparing not only the Norway spruce PUTs, but also PUTs from the close relatives Picea glauca and Picea sitchensis to both Pinus taeda and Taxus mairei, we obtained estimates of synonymous and non-synonymous divergence among conifer species. In addition, we detected close to 15,000 SNPs of high quality and estimated gene expression differences between samples collected under dark and light conditions. CONCLUSIONS: Our study yielded a large number of single nucleotide polymorphisms as well as estimates of gene expression on transcriptome scale. In agreement with a recent study we find that the synonymous substitution rate per year (0.6 × 10-09 and 1.1 × 10-09) is an order of magnitude smaller than values reported for angiosperm herbs. However, if one takes generation time into account, most of this difference disappears. The estimates of the dN/dS ratio (non-synonymous over synonymous divergence) reported here are in general much lower than 1 and only a few genes showed a ratio larger than 1.

Evolution of the PEBP gene family in plants: functional diversification in seed plant evolution.

The phosphatidyl ethanolamine-binding protein (PEBP) gene family is present in all eukaryote kingdoms, with three subfamilies identified in angiosperms (FLOWERING LOCUS T [FT], MOTHER OF FT AND TFL1 [MFT], and TERMINAL FLOWER1 [TFL1] like). In angiosperms, PEBP genes have been shown to function both as promoters and suppressors of flowering and to control plant architecture. In this study, we focus on previously uncharacterized PEBP genes from gymnosperms. Extensive database searches suggest that gymnosperms possess only two types of PEBP genes, MFT-like and a group that occupies an intermediate phylogenetic position between the FT-like and TFL1-like (FT/TFL1-like). Overexpression of Picea abies PEBP genes in Arabidopsis (Arabidopsis thaliana) suggests that the FT/TFL1-like genes (PaFTL1 and PaFTL2) code for proteins with a TFL1-like function. However, PaFTL1 and PaFTL2 also show highly divergent expression patterns. While the expression of PaFTL2 is correlated with annual growth rhythm and mainly confined to needles and vegetative and reproductive buds, the expression of PaFTL1 is largely restricted to microsporophylls of male cones. The P. abies MFT-like genes (PaMFT1 and PaMFT2) show a predominant expression during embryo development, a pattern that is also found for many MFT-like genes from angiosperms. P. abies PEBP gene expression is primarily detected in tissues undergoing physiological changes related to growth arrest and dormancy. A first duplication event resulting in two families of plant PEBP genes (MFT-like and FT/TFL1-like) seems to coincide with the evolution of seed plants, in which independent control of bud and seed dormancy was required, and the second duplication resulting in the FT-like and TFL1-like clades probably coincided with the evolution of angiosperms.

Does the core circadian clock in the moss Physcomitrella patens (Bryophyta) comprise a single loop?

BACKGROUND: The endogenous circadian clock allows the organism to synchronize processes both to daily and seasonal changes. In plants, many metabolic processes such as photosynthesis, as well as photoperiodic responses, are under the control of a circadian clock. Comparative studies with the moss Physcomitrella patens provide the opportunity to study many aspects of land plant evolution. Here we present a comparative overview of clock-associated components and the circadian network in the moss P. patens. RESULTS: The moss P. patens has a set of conserved circadian core components that share genetic relationship and gene expression patterns with clock genes of vascular plants. These genes include Myb-like transcription factors PpCCA1a and PpCCA1b, pseudo-response regulators PpPRR1-4, and regulatory elements PpELF3, PpLUX and possibly PpELF4. However, the moss lacks homologs of AtTOC1, AtGI and the AtZTL-family of genes, which can be found in all vascular plants studied here. These three genes constitute essential components of two of the three integrated feed-back loops in the current model of the Arabidopsis circadian clock mechanism. Consequently, our results suggest instead a single loop circadian clock in the moss. Possibly as a result of this, temperature compensation of core clock gene expression appears to be decreased in P. patens. CONCLUSIONS: This study is the first comparative overview of the circadian clock mechanism in a basal land plant, the moss P. patens. Our results indicate that the moss clock mechanism may represent an ancestral state in contrast to the more complex and partly duplicated structure of subsequent land plants. These findings may provide insights into the understanding of the evolution of circadian network topology.

Population genomics reveals lack of greater white-fronted introgression into the Swedish lesser white-fronted goose.

Interspecific introgression is considered a potential threat to endangered taxa. One example where this has had a major impact on conservation policy is the lesser white-fronted goose (LWfG). After a dramatic decline in Sweden, captive breeding birds were released between 1981-1999 with the aim to reinforce the population. However, the detection of greater white-fronted goose (GWfG) mitochondrial DNA in the LWfG breeding stock led to the release program being dismantled, even though the presence of GWfG introgression in the actual wild Swedish LWfG population was never documented. To examine this, we sequenced the complete genomes of 21 LWfG birds from the Swedish, Russian and Norwegian populations, and compared these with genomes from other goose species, including the GWfG. We found no evidence of interspecific introgression into the wild Swedish LWfG population in either nuclear genomic or mitochondrial data. Moreover, Swedish LWfG birds are genetically distinct from the Russian and Norwegian populations and display comparatively low genomic diversity and high levels of inbreeding. Our findings highlight the utility of genomic approaches in providing scientific evidence that can help improve conservation management as well as policies for breeding and reinforcement programmes.

Multilocus patterns of nucleotide diversity, linkage disequilibrium and demographic history of Norway spruce [Picea abies (L.) Karst].

DNA polymorphism at 22 loci was studied in an average of 47 Norway spruce [Picea abies (L.) Karst.] haplotypes sampled in seven populations representative of the natural range. The overall nucleotide variation was limited, being lower than that observed in most plant species so far studied. Linkage disequilibrium was also restricted and did not extend beyond a few hundred base pairs. All populations, with the exception of the Romanian population, could be divided into two main domains, a Baltico-Nordic and an Alpine one. Mean Tajima's D and Fay and Wu's H across loci were both negative, indicating the presence of an excess of both rare and high-frequency-derived variants compared to the expected frequency spectrum in a standard neutral model. Multilocus neutrality tests based on D and H led to the rejection of the standard neutral model and exponential growth in the whole population as well as in the two main domains. On the other hand, in all three cases the data are compatible with a severe bottleneck occurring some hundreds of thousands of years ago. Hence, demographic departures from equilibrium expectations and population structure will have to be accounted for when detecting selection at candidate genes and in association mapping studies, respectively.

Association mapping in Salix viminalis L. (Salicaceae) - identification of candidate genes associated with growth and phenology.

Willow species (Salix) are important as short-rotation biomass crops for bioenergy, which creates a demand for faster genetic improvement and breeding through deployment of molecular marker-assisted selection (MAS). To find markers associated with important adaptive traits, such as growth and phenology, for use in MAS, we genetically dissected the trait variation of a Salix viminalis (L.) population of 323 accessions. The accessions were sampled throughout northern Europe and were established at two field sites in Pustnäs, Sweden, and at Woburn, UK, offering the opportunity to assess the impact of genotype-by-environment interactions (G × E) on trait-marker associations. Field measurements were recorded for growth and phenology traits. The accessions were genotyped using 1536 SNP markers developed from phenology candidate genes and from genes previously observed to be differentially expressed in contrasting environments. Association mapping between 1233 of these SNPs and the measured traits was performed taking into account population structure and threshold selection bias. At a false discovery rate (FDR) of 0.2, 29 SNPs were associated with bud burst, leaf senescence, number of shoots or shoot diameter. The percentage of accession variation (Radj2) explained by these associations ranged from 0.3% to 4.4%, suggesting that the studied traits are controlled by many loci of limited individual impact. Despite this, a SNP in the EARLY FLOWERING 3 gene was repeatedly associated (FDR < 0.2) with bud burst. The rare homozygous genotype exhibited 0.4-1.0 lower bud burst scores than the other genotype classes on a five-grade scale. Consequently, this marker could be promising for use in MAS and the gene deserves further study. Otherwise, associations were less consistent across sites, likely due to their small Radj2 estimates and to considerable G × E interactions indicated by multivariate association analyses and modest trait accession correlations across sites (0.32-0.61).

Coexistence of the social types: genetic population structure in the ant Formica exsecta.

The ant Formica exsecta has two types of colonies that exist in sympatry but usually as separate subpopulations: colonies with simple social organization and single queens (M type) or colonial networks with multiple queens (P type). We used both nuclear (DNA microsatellites) and mitochondrial markers to study the transition between the social types, and the contribution of males and females in gene flow within and between the types. Our results showed that the social types had different spatial genetic structures. The M subpopulations formed a fairly uniform population, whereas the P subpopulations were, on average, more differentiated from each other than from the nearby M subpopulations and could have been locally established from the M-type colonies, followed by philopatric behavior and restricted emigration of females. Thus, the relationship between the two social types resembles that of source (M type) and sink (P type) populations. The comparison of mitochondrial (phiST) and nuclear (FST) differentiation indicates that the dispersal rate of males is four to five times larger than that of females both among the P-type subpopulations and between the social types. Our results suggest that evolution toward complex social organization can have an important effect on genetic population structure through changes in dispersal behavior associated with different sociogenetic organizations.

Patterns of nucleotide diversity at photoperiod related genes in Norway spruce [Picea abies (L.) Karst].

The ability of plants to track seasonal changes is largely dependent on genes assigned to the photoperiod pathway, and variation in those genes is thereby important for adaptation to local day length conditions. Extensive physiological data in several temperate conifer species suggest that populations are adapted to local light conditions, but data on the genes underlying this adaptation are more limited. Here we present nucleotide diversity data from 19 genes putatively involved in photoperiodic response in Norway spruce (Picea abies). Based on similarity to model plants the genes were grouped into three categories according to their presumed position in the photoperiod pathway: photoreceptors, circadian clock genes, and downstream targets. An HKA (Hudson, Kreitman and Aquade) test showed a significant excess of diversity at photoreceptor genes, but no departure from neutrality at circadian genes and downstream targets. Departures from neutrality were also tested with Tajima's D and Fay and Wu's H statistics under three demographic scenarios: the standard neutral model, a population expansion model, and a more complex population split model. Only one gene, the circadian clock gene PaPRR3 with a highly positive Tajima's D value, deviates significantly from all tested demographic scenarios. As the PaPRR3 gene harbours multiple non-synonymous variants it appears as an excellent candidate gene for control of photoperiod response in Norway spruce.

Conserved function of core clock proteins in the gymnosperm Norway spruce (Picea abies L. Karst).

From studies of the circadian clock in the plant model species Arabidopsis (Arabidopsis thaliana), a number of important properties and components have emerged. These include the genes CIRCADIAN CLOCK ASSOCIATED 1 (CCA1), GIGANTEA (GI), ZEITLUPE (ZTL) and TIMING OF CAB EXPRESSION 1 (TOC1 also known as PSEUDO-RESPONSE REGULATOR 1 (PRR1)) that via gene expression feedback loops participate in the circadian clock. Here, we present results from ectopic expression of four Norway spruce (Picea abies) putative homologs (PaCCA1, PaGI, PaZTL and PaPRR1) in Arabidopsis, their flowering time, circadian period length, red light response phenotypes and their effect on endogenous clock genes were assessed. For PaCCA1-ox and PaZTL-ox the results were consistent with Arabidopsis lines overexpressing the corresponding Arabidopsis genes. For PaGI consistent results were obtained when expressed in the gi2 mutant, while PaGI and PaPRR1 expressed in wild type did not display the expected phenotypes. These results suggest that protein function of PaCCA1, PaGI and PaZTL are at least partly conserved compared to Arabidopsis homologs, however further studies are needed to reveal the protein function of PaPRR1. Our data suggest that components of the three-loop network typical of the circadian clock in angiosperms were present before the split of gymnosperms and angiosperms.

Distribution of long-range linkage disequilibrium and Tajima's D values in Scandinavian populations of Norway Spruce (Picea abies).

The site frequency spectrum of mutations (SFS) and linkage disequilibrium (LD) are the two major sources of information in population genetics studies. In this study we focus on the levels of LD and the SFS and on the effect of sample size on summary statistics in 10 Scandinavian populations of Norway spruce. We found that previous estimates of a low level of LD were highly influenced by both sampling strategy and the fact that data from multiple loci were analyzed jointly. Estimates of LD were in fact heterogeneous across loci and increased within individual populations compared with the estimate from the total data. The variation in levels of LD among populations most likely reflects different demographic histories, although we were unable to detect population structure by using standard approaches. As in previous studies, we also found that the SFS-based test Tajima's D was highly sensitive to sample size, revealing that care should be taken to draw strong conclusions from this test when sample size is small. In conclusion, the results from this study are in line with recent studies in other conifers that have revealed a more complex and variable pattern of LD than earlier studies suggested and with studies in trees and humans that suggest that Tajima's D is sensitive to sample size. This has large consequences for the design of future association and population genetic studies in Norway spruce.

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.

Non-radioactive in situ hybridization protocol applicable for Norway spruce and a range of plant species.

The high-throughput expression analysis technologies available today give scientists an overflow of expression profiles but their resolution in terms of tissue specific expression is limited because of problems in dissecting individual tissues. Expression data needs to be confirmed and complemented with expression patterns using e.g. in situ hybridization, a technique used to localize cell specific mRNA expression. The in situ hybridization method is laborious, time-consuming and often requires extensive optimization depending on species and tissue. In situ experiments are relatively more difficult to perform in woody species such as the conifer Norway spruce (Picea abies). Here we present a modified DIG in situ hybridization protocol, which is fast and applicable on a wide range of plant species including P. abies. With just a few adjustments, including altered RNase treatment and proteinase K concentration, we could use the protocol to study tissue specific expression of homologous genes in male reproductive organs of one gymnosperm and two angiosperm species; P. abies, Arabidopsis thaliana and Brassica napus. The protocol worked equally well for the species and genes studied. AtAP3 and BnAP3 were observed in second and third whorl floral organs in A. thaliana and B. napus and DAL13 in microsporophylls of male cones from P. abies. For P. abies the proteinase K concentration, used to permeablize the tissues, had to be increased to 3 g/ml instead of 1 g/ml, possibly due to more compact tissues and higher levels of phenolics and polysaccharides. For all species the RNase treatment was removed due to reduced signal strength without a corresponding increase in specificity. By comparing tissue specific expression patterns of homologous genes from both flowering plants and a coniferous tree we demonstrate that the DIG in situ protocol presented here, with only minute adjustments, can be applied to a wide range of plant species. Hence, the protocol avoids both extensive species specific optimization and the laborious use of radioactively labeled probes in favor of DIG labeled probes. We have chosen to illustrate the technically demanding steps of the protocol in our film.

A Norway spruce FLOWERING LOCUS T homolog is implicated in control of growth rhythm in conifers.

Growth in perennial plants possesses an annual cycle of active growth and dormancy that is controlled by environmental factors, mainly photoperiod and temperature. In conifers and other nonangiosperm species, the molecular mechanisms behind these responses are currently unknown. In Norway spruce (Picea abies L. Karst.) seedlings, growth cessation and bud set are induced by short days and plants from southern latitudes require at least 7 to 10 h of darkness, whereas plants from northern latitudes need only 2 to 3 h of darkness. Bud burst, on the other hand, is almost exclusively controlled by temperature. To test the possible role of Norway spruce FLOWERING LOCUS T (FT)-like genes in growth rhythm, we have studied expression patterns of four Norway spruce FT family genes in two populations with a divergent bud set response under various photoperiodic conditions. Our data show a significant and tight correlation between growth rhythm (both bud set and bud burst), and expression pattern of one of the four Norway spruce phosphatidylethanolamine-binding protein gene family members (PaFT4) over a variety of experimental conditions. This study strongly suggests that one Norway spruce homolog to the FT gene, which controls flowering in angiosperms, is also a key integrator of photoperiodic and thermal signals in the control of growth rhythms in gymnosperms. The data also indicate that the divergent adaptive bud set responses of northern and southern Norway spruce populations, both to photoperiod and light quality, are mediated through PaFT4. These results provide a major advance in our understanding of the molecular control of a major adaptive trait in conifers and a tool for further molecular studies of adaptive variation in plants.