Niclas Gyllenstrand

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.

Colony kin structure and male production in Dolichovespula wasps.

In annual hymenopteran societies headed by a single outbred queen, paternity (determined by queen mating frequency and sperm use) is the sole variable affecting colony kin structure and is therefore a key predictor of colony reproductive characteristics. Here we investigate paternity and male production in five species of Dolichovespula wasps. Twenty workers from each of 10 colonies of each of five species, 1000 workers in total, were analysed at three DNA microsatellite loci to estimate paternity. To examine the relationship between kin structure and reproductive behaviour, worker ovary activation was assessed by dissection and the maternal origin of adult males was assessed by DNA microsatellites. Effective paternity was low in all species (D.media 1.08, D. maculata 1.0, D. sylvestris 1.15, D. norwegica 1.08 and D. saxonica 1.35), leading to the prediction of queen–worker conflict over male production. In support of this, workers with full‐size eggs in their ovaries (four out of five species) and adult males that were workers’ sons (all five species) were found in queenright colonies. However, workers were only responsible for a minority of male production (D.media 7.4%, D. maculata 20.9%, D. sylvestris 9.8%, D. norwegica 2.6% and D. saxonica 34.6%) suggesting that the queen maintains considerable reproductive power over the workers. Kin structure and reproductive conflict in Dolichovespula contrast with their sister group Vespula. Dolichovespula is characterized by low paternity, worker reproduction, and queen–worker conflict and Vespula by high paternity, effective worker policing and absence of worker reproduction. The trend revealed by this comparison is as predicted by kin selection theory suggesting that colony kin structure has been pivotal in the evolution of the yellowjacket wasps.

Conservation genetics of the wood ant, Formica lugubris , in a fragmented landscape

Various intrinsic factors connected to the special features of sociality influence the persistence of social insect populations, including low effective population size, reduced amount of genetic variation easily leading to inbreeding depression, and spatially structured populations. In this work, we studied an isolated, small and fragmented population system of the red wood ant Formica lugubris, and evaluated the impact of social and genetic population structure on the persistence and conservation of the populations. The effective population size was large in our study population because all nests were polygynous. As a result, and despite the apparent isolation, the amount of nuclear genetic variability was similar to that in a nonisolated population system. Lack of inbreeding, as well as a high level of variability, indirectly suggests that this population does not suffer from inbreeding depression. The spatial distribution of genetic variation between local populations suggests intensive, but strongly male‐biased, nuclear gene flow. Thus, the persistence of this population system does not seem to be threatened by any immediate social or genetic factor, but colonization of new habitat patches may be difficult because of restricted female dispersal.

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.

COEXISTENCE OF THE SOCIAL TYPES: GENETIC POPULATION STRUCTURE IN THE ANT FORMICA EXSECTA

Abstract —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 multipe queens (P type). We used both nuclear (DNA microsatellites) and mitochondrial markers to study the taransition between the social types, and the contribution of males nd femeles in gene flow within and between the types. Our results showed that the social types had different spatial genetic stractures. The M subpoplations formed a fairly uniform population, whereas the P subpopulations were, on average, more differentiated from each other than from the nearby M bubpoplations and could have been locally established from the M‐type colonies, followed by philopatric behavior and restrcted emigration of females. Thus, the relationship between the two social types resembles that of souce (M type) and sink (P type) populations. The comparison of mitochondrial (Φst) 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 socigenetic orgamizations.

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.