Thomas Källman

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 Tajimas D and Fay and Wus 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.

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.

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.

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 (FST = 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 FST. 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.

Demographic histories of four spruce (Picea) species of the Qinghai-Tibetan Plateau and neighboring areas inferred from multiple nuclear loci

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.

Early evolution of the MFT-like gene family in plants.

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.

New insights on the speciation history and nucleotide diversity of three boreal spruce species and a Tertiary relict.

In all, 10 nuclear loci were re-sequenced in four spruce species. Three of the species are boreal species with very large natural ranges: Picea mariana and P. glauca are North American, and P. abies, is Eurasian. The fourth species, P. breweriana, is a Tertiary relict from Northern California, with a very small natural range. Although the boreal species population sizes have fluctuated through the Ice Ages, P. breweriana is believed to have had a rather stable population size through the Quaternary. Indeed, the average Tajimas D was close to zero in this species and negative in the three boreal ones. Reflecting differences in current population sizes, nucleotide diversity was an order of magnitude lower in P. breweriana than in the boreal species. This is in contrast to the similar and high levels of heterozygosity observed in previous studies at allozyme loci across species. As the species have very different histories and effective population sizes, selection at allozyme loci rather than demography appears to be a better explanation for this discrepancy. Parameters of Isolation-with-Migration (IM) models were also estimated for pairs of species. Shared polymorphisms were extensive and fixed polymorphisms few. Divergence times were much shorter than those previously reported. There was also evidence of historical gene flow between P. abies and P. glauca. The latter was more closely related to P. abies than to its sympatric relative P. mariana. This last result suggests that North American and Eurasian species might have been geographically much closer in the recent past than they are today.

Contrasting demographic history and population structure in Capsella rubella and Capsella grandiflora, two closely related species with different mating systems

Both mating system and population history can have large impacts on genetic diversity and population structure. Here, we use multilocus sequence data to investigate how these factors impact two closely related Brassicaceae species: the selfing Capsella rubella and the outcrossing C. grandiflora. To do this, we have sequenced 16 loci in approximately 70 individuals from 7 populations of each species. Patterns of population structure differ strongly between the two species. In C. grandiflora, we observe an isolation‐by‐distance pattern and identify three clearly delineated genetic groups. In C. rubella, where we estimate the selfing rate to be 0.90–0.94, the pattern is less clear with some sampling populations forming separate genetic clusters while others are highly mixed. The two species also have divergent histories. Our analysis gives support for a bottleneck approximately 73 kya (20–139 kya) in C. rubella, which most likely represents speciation from C. grandiflora. In C. grandiflora, there is moderate support for the standard neutral model in 2 of 3 genetic clusters, while the third cluster and the total data set show evidence of expansion. It is clear that mating system has an impact on these two species, for example affecting the level of genetic variation and the genetic structure. However, our results also clearly show that a combination of past and present processes, some of which are not affected by mating system, is needed to explain the differences between C. rubella and C. grandiflora.

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

BackgroundThe 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.ResultsThe 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.ConclusionsThis 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.

A Significant Fraction of 21-Nucleotide Small RNA Originates from Phased Degradation of Resistance Genes in Several Perennial Species

NBS-LRR genes, which constitute the major class of plant innate immune receptors, are massively degraded through microRNA-guided generation of secondary small interfering RNAs 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.