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

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 (MFTlike 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.

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.L OCAL adaptation is a key process in the evolution of species. Understanding how local adaptation is established and maintained, however, is especially difficult as its establishment is contingent upon historical conditions and its maintenance depends on the balance among conflicting evolutionary forces (e.g., Yeaman and Otto 2011). It is a particularly challenging task in forest trees, because they have long generation times and therefore cannot be easily manipulated experimentally. For instance, transfer experiments are theoretically possible but practically difficult to implement. On the other hand, the analysis of the strong latitudinal clines displayed by forest trees for potentially adaptive traits such as bud set (Dormling 1973;Savolainen et al. 2007;Aitken et al. 2008) can provide crucial information on the forces involved in local adaptation and, in particular, on the relative parts played by demography and selection in the establishment of the cline. Furthermore, phenology in general, and flowering time and bud set in particular, have been extensively studied and strong candidate genes are available, many of which belong to the photoperiodic pathway including the circadian clock (Gyllenstrand et al. 2007;Albani and Coupland 2010;Bergelson and Roux 2010;Fornara et al. 2...

Nucleotide variation at 12-16 nuclear loci was studied in three spruce species from the Qinghai-Tibetan Plateau (QTP), Picea likiangensis, P. wilsonii, and P. purpurea, and one species from the Tian Shan mountain range, P. schrenkiana. Silent nucleotide diversity was limited in P. schrenkiana and high in the three species from the QTP, with values higher than in boreal spruce species, despite their much more restricted distributions compared with that of the boreal species. In contrast to European boreal species that have experienced severe bottlenecks in the past, coalescent-based analysis suggests that DNA polymorphism in the species from the QTP and adjacent areas is compatible with the standard neutral model (P. likiangensis, P. wilsonii, and P. schrenkiana) or with population growth (P. purpurea). In order to test if P. purpurea is a diploid hybrid of P. likiangensis and P. wilsonii, we used a combination of approaches, including model-based inference of population structure, isolation-with-migration models, and recent theoretical results on the effect of introgression on the geographic distribution of diversity. In contrast to the three other species, each of which was predominantly assigned to a single cluster in the Structure analysis, P. purpurea individuals were scattered over the three main clusters and not, as we had expected, confined to the P. likiangensis and P. wilsonii clusters. Furthermore, the contribution of P. schrenkiana was by far the largest one. In agreement with this, the divergence between P. purpurea and P. schrenkiana was lower than the divergence of either P. likiangensis or P. wilsonii from P. schrenkiana. These results, together with previous ones showing that P. purpurea and P. wilsonii share the same haplotypes at both chloroplast and mitochondrial markers, suggest that P. purpurea has a complex origin, possibly involving additional species.

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.Trees and other perennial plants must adapt their growth rhythm to seasonal changes in the environment. To a large extent, this adaptation is genetically controlled (Howe et al., 2003). A clear example is the strong clinal pattern of growth cessation and bud set in seedlings of the conifer Norway spruce (Picea abies L. Karst.; Ekberg et al., 1976).Like most conifers, Norway spruce has a long juvenile phase of about 20 years before the first cones are formed. In first-year seedlings, the annual cycle (Fig. 1) can be summarized as (1) shoot extension stops and terminal buds are set in late summer in response to a shortening photoperiod, after which the cambium ceases growth, needle primordia are initiated within the buds, and frost tolerance begins to increase; (2) rest dormancy (endodormancy) develops in the meristems during autumn after bud set and, with exposure to chilling temperatures (2°C-10°C), changes into quiescence dormancy (ectodormancy) by midwinter, when frost tolerance is maximal; and (3) opening of the bud scales (bud burst) occurs in spring after a temperature sum (TS) has been attained. The extension growth of first-year seedlings consists of the expansion of stem units formed in the current season. This free growth is in the following years (Fig.

Angiosperm genes sharing a conserved phosphatidylethanolamine-binding (PEPB) domain have been shown to be involved in the control of shoot meristem identity and flowering time. The family is divided into three subfamilies, FT-like, TFL1-like and MFT-like. This study is focused on the evolution of the MFT-like clade, suggested to be ancestral to the two other clades. We report that the bryophyte Physcomitrella patens and the lycopod Selaginella moellendorfii contain four and two MFT-like genes respectively. Neither species have any FT or TFL1-like genes. Furthermore, we have identified a new subclade of MFT-like genes in Angiosperms. Quantitative expression analysis of MFT-like genes in Physcomitrella patens reveals that the expression patterns are circadian and reaches maximum in gametangia and sporophytes. Our data suggest that the occurrence FT and TFL1-like genes, is associated with the evolution of seed plants. Expression data for Physcomitrella MFT-like genes implicates an involvement in the development of reproductive tissues in the moss.

Boreal species were repeatedly exposed to ice ages and went through cycles of contraction and expansion while sister species alternated periods of contact and isolation. The resulting genetic structure is consequently complex, and demographic inferences are intrinsically challenging. The range of Norway spruce (Picea abies) and Siberian spruce (Picea obovata) covers most of northern Eurasia; yet their geographical limits and histories remain poorly understood. To delineate the hybrid zone between the two species and reconstruct their joint demographic history, we analysed variation at nuclear SSR and mitochondrial DNA in 102 and 88 populations, respectively. The dynamics of the hybrid zone was analysed with approximate Bayesian computation (ABC) followed by posterior predictive structure plot reconstruction and the presence of barriers across the range tested with estimated effective migration surfaces. To estimate the divergence time between the two species, nuclear sequences from two well-separated populations of each species were analysed with ABC. Two main barriers divide the range of the two species: one corresponds to the hybrid zone between them, and the other separates the southern and northern domains of Norway spruce. The hybrid zone is centred on the Urals, but the genetic impact of Siberian spruce extends further west. The joint distribution of mitochondrial and nuclear variation indicates an introgression of mitochondrial DNA from Norway spruce into Siberian spruce. Overall, our data reveal a demographic history where the two species interacted frequently and where migrants originating from the Urals and the West Siberian Plain recolonized northern Russia and Scandinavia using scattered refugial populations of Norway spruce as stepping stones towards the west.

A consensus linkage map of Picea abies, an economically important conifer, was constructed based on the segregation of 686 SNP markers in a F1 progeny population consisting of 247 individuals. The total length of 1889.2 cM covered 96.5% of the estimated genome length and comprised 12 large linkage groups, corresponding to the number of haploid P. abies chromosomes. The sizes of the groups (from 5.9 to 9.9% of the total map length) correlated well with previous estimates of chromosome sizes (from 5.8 to 10.8% of total genome size). Any locus in the genome has a 97% probability to be within 10 cM from a mapped marker, which makes the map suited for QTL mapping. Infecting the progeny trees with the root rot pathogen Heterobasidion parviporum allowed for mapping of four different resistance traits: lesion length at the inoculation site, fungal spread within the sapwood, exclusion of the pathogen from the host after initial infection, and ability to prevent the infection from establishing at all. These four traits were associated with two, four, four and three QTL regions respectively of which none overlapped between the traits. Each QTL explained between 4.6 and 10.1% of the respective traits phenotypic variation. Although the QTL regions contain many more genes than the ones represented by the SNP markers, at least four markers within the confidence intervals originated from genes with known function in conifer defence; a leucoanthocyanidine reductase, which has previously been shown to upregulate during H. parviporum infection, and three intermediates of the lignification process; a hydroxycinnamoyl CoA shikimate/quinate hydroxycinnamoyltransferase, a 4-coumarate CoA ligase, and a R2R3-MYB transcription factor.