Understanding the genetic basis of the timing of bud set, an important trait in conifers, is relevant for adaptation and forestry practice. In common garden experiments, both Scots pine (Pinus sylvestris) and Norway spruce (Picea abies) show a latitudinal cline in the trait. We compared the regulation of their bud set biology by examining the expression of PsFTL2, a Pinus sylvestris homolog to PaFTL2, a FLOWERING LOCUS T/TERMINAL FLOWER 1 (FT/TFL1)-like gene, the expression levels of which have been found previously to be associated with the timing of bud set in Norway spruce. In a common garden study, we analyzed the relationship of bud phenology under natural and artificial photoperiods and the expression of PsFTL2 in a set of Scots pine populations from different latitudes. The expression of PsFTL2 increased in the needles preceding bud set and decreased during bud burst. In the northernmost population, even short night periods were efficient to trigger this expression, which also increased earlier under all photoperiodic regimes compared with the southern populations. Despite the different biology, with few limitations, the two conifers that diverged 140 million yr ago probably share an association of FTL2 with bud set, pointing to a common mechanism for the timing of growth cessation in conifers.
To investigate the potential of association genetics for willow breeding, Salix viminalis germplasm was assembled from UK and Swedish collections (comprising accessions from several European countries) and new samples collected from nature. A subset of the germplasm was planted at two sites (UK and Sweden), genotyped using 38 SSR markers and assessed for phenological and biomass traits. Population structure, genetic differentiation (F-ST) and quantitative trait differentiation (Q(ST)) were investigated. The extent and patterns of trait adaptation were assessed by comparing F-ST and Q(ST) parameters. Of the 505 genotyped diploid accessions, 27 % were not unique. Genetic diversity was high: 471 alleles was amplified; the mean number of alleles per locus was 13.46, mean observed heterozygosity was 0.55 and mean expected heterozygosity was 0.62. Bayesian clustering identified four subpopulations which generally corresponded to Western Russia, Western Europe, Eastern Europe and Sweden. All pairwise F-ST values were highly significant (p<0.001) with the greatest genetic differentiation detected between the Western Russian and the Western European subpopulations (F-ST = 0.12), and the smallest between the Swedish and Eastern European populations (F-ST = 0.04). The Swedish population also had the highest number of identical accessions, supporting the view that S. viminalis was introduced into this country and has been heavily influenced by humans. Q(ST) values were high for growth cessation and leaf senescence, and to some extent stem diameter, but low for bud burst time and shoot number. Overall negative clines between longitudinal coordinates and leaf senescence, bud burst and stem diameter were also found.
We investigated species divergence, present and past gene flow, levels of nucleotide polymorphism, and linkage disequilibrium in two willows from the plant genus Salix. Salix belongs together with Populus to the Salicaceae family; however, most population genetic studies of Salicaceae have been performed in Populus, the model genus in forest biology. Here we present a study on two closely related willow species Salix viminalis and S. schwerinii, in which we have resequenced 33 and 32 nuclear gene segments representing parts of 18 nuclear loci in 24 individuals for each species. We used coalescent simulations and estimated the split time to around 600,000 years ago and found that there is currently limited gene flow between the species. Mean intronic nucleotide diversity across gene segments was slightly higher in S. schwerinii (πi = 0.00849) than in S. viminalis (πi = 0.00655). Compared with other angiosperm trees, the two willows harbor intermediate levels of silent polymorphisms. The decay of linkage disequilibrium was slower in S. viminalis compared with S. schwerinii, and we speculate that this is due to different demographic histories as S. viminalis has been partly domesticated in Europe.
BACKGROUND: Salix (willow) and Populus (poplar) are members of the Salicaceae family and they share many ecological as well as genetic and genomic characteristics. The interest of using willow for biomass production is growing, which has resulted in increased pressure on breeding of high yielding and resistant clones adapted to different environments. The main purpose of this work was to develop dense genetic linkage maps for mapping of traits related to yield and resistance in willow. We used the Populus trichocarpa genome to extract evenly spaced markers and mapped the orthologous loci in the willow genome. The marker positions in the two genomes were used to study genome evolution since the divergence of the two lineages some 45 mya. RESULTS: We constructed two linkage maps covering the 19 linkage groups in willow. The most detailed consensus map, S1, contains 495 markers with a total genetic distance of 2477 cM and an average distance of 5.0 cM between the markers. The S3 consensus map contains 221 markers and has a total genetic distance of 1793 cM and an average distance of 8.1 cM between the markers. We found high degree of synteny and gene order conservation between willow and poplar. There is however evidence for two major interchromosomal rearrangements involving poplar LG I and XVI and willow LG Ib, suggesting a fission or a fusion in one of the lineages, as well as five intrachromosomal inversions. The number of silent substitutions were three times lower (median: 0.12) between orthologs than between paralogs (median: 0.37 - 0.41). CONCLUSIONS: The relatively slow rates of genomic change between willow and poplar mean that the genomic resources in poplar will be most useful in genomic research in willow, such as identifying genes underlying QTLs of important traits. Our data suggest that the whole-genome duplication occurred long before the divergence of the two genera, events which have until now been regarded as contemporary. Estimated silent substitution rates were 1.28 x 10-9 and 1.68 x 10-9 per site and year, which are close to rates found in other perennials but much lower than rates in annuals.
The evolution of land flora transformed the terrestrial environment. Land plants evolved from an ancestral charophycean alga from which they inherited developmental, biochemical, and cell biological attributes. Additional biochemical and physiological adaptations to land, and a life cycle with an alternation between multicellular haploid and diploid generations that facilitated efficient dispersal of desiccation tolerant spores, evolved in the ancestral land plant. We analyzed the genome of the liverwort Marchantia polymorpha, a member of a basal land plant lineage. Relative to charophycean algae, land plant genomes are characterized by genes encoding novel biochemical pathways, new phytohormone signaling pathways (notably auxin), expanded repertoires of signaling pathways, and increased diversity in some transcription factor families. Compared with other sequenced land plants, M. polymorpha exhibits low genetic redundancy in most regulatory pathways, with this portion of its genome resembling that predicted for the ancestral land plant. PAPERCLIP.
To gain new insights into the mechanism underlying cytoplasmic male sterility (CMS), we compared the nuclear gene expression profiles of flowers of a Brassica napus CMS line with that of the fertile B. napus maintainer line using Arabidopsis thaliana flower-specific cDNA microarrays. The CMS line used has a B. napus nuclear genome, but has a rearranged mitochondrial (mt) genome consisting of both B. napus and A. thaliana DNA. Gene expression profiling revealed that a large number of genes differed in expression between the two lines. For example, nuclear genes coding for proteins that are involved in protein import into organelles, genes expressed in stamens and pollen, as well as genes implicated in either cell-wall remodeling or architecture, were repressed in the CMS line compared with B. napus. These results show that the mt genome of the CMS line strongly influences nuclear gene expression, and thus reveal the importance of retrograde signalling between the mitochondria and the nucleus. Furthermore, flowers of the CMS line are characterized by a replacement of stamens with carpelloid organs, and thus partially resemble the APETALA3 (AP3) and PISTILLATA (PI) mutants. In accordance with this phenotype, AP3 expression was downregulated in the stamens, shortly before these organs developed carpelloid characteristics, even though it was initiated correctly. Repression of PI succeeded that of AP3 and might be a consequence of a loss of AP3 activity. These results suggest that AP3 expression in stamens depends on proper mt function and a correct nuclear-mt interaction, and that mt alterations cause the male sterility phenotype of the CMS line.
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.
Background: A detailed knowledge about which genes are expressed in which tissues and at which developmental stage 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 give both a snap-shot of the transcribed part of a species genome and simultaneously 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 were, together with publicly available expressed sequence tag (EST) data from Norway spruce, used to create a reference transcriptome. Of the 38,419 PUTs (putative unique transcripts) longer than 150 bp in this reference assembly, 59% 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 difference between samples collected during 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, but if one takes generation time in to account, most of this difference disappear. The estimates of the non-synonymous over the synonymous divergence (dN/dS ratio) reported here is in general much lower than 1 and only a few genes showed a ratio larger than 1.
Gene duplication followed by functional divergence in the event of polyploidization is a major contributor to evolutionary novelties. The Brassica genus evolved from a common ancestor after whole-genome triplication. Here, we studied the evolutionary and functional features of Brassica spp. homologs to Tic40 (for translocon at the inner membrane of chloroplasts with 40 kDa). Four Tic40 loci were identified in allotetraploid Brassica napus and two loci in each of three basic diploid Brassica spp. Although these Tic40 homologs share high sequence identities and similar expression patterns, they exhibit altered functional features. Complementation assays conducted on Arabidopsis thaliana tic40 and the B. napus male-sterile line 7365A suggested that all Brassica spp. Tic40 homologs retain an ancestral function similar to that of AtTic40, whereas BolC9.Tic40 in Brassica oleracea and its ortholog in B. napus, BnaC9.Tic40, in addition, evolved a novel function that can rescue the fertility of 7365A. A homologous chromosomal rearrangement placed bnac9.tic40 originating from the A genome (BraA10.Tic40) as an allele of BnaC9.Tic40 in the C genome, resulting in phenotypic variation for male sterility in the B. napus near-isogenic two-type line 7365AB. Assessment of the complementation activity of chimeric B. napus Tic40 domain-swapping constructs in 7365A suggested that amino acid replacements in the carboxyl terminus of BnaC9.Tic40 cause this functional divergence. The distribution of these amino acid replacements in 59 diverse Brassica spp. accessions demonstrated that the neofunctionalization of Tic40 is restricted to B. oleracea and its derivatives and thus occurred after the divergence of the Brassica spp. A, B, and C genomes.
The plant hormone auxin (indole-3-acetic acid [IAA]) has previously been suggested to regulate diverse forms of dormancy in both seed plants and liverworts. Here, we use loss-and gain-of-function alleles for auxin synthesis-and signaling-related genes, as well as pharmacological approaches, to study how auxin regulates development and dormancy in the gametophyte generation of the liverwort Marchantia polymorpha. We found that M. polymorpha possess the smallest known toolkit for the indole-3-pyruvic acid (IPyA) pathway in any land plant and that this auxin synthesis pathway mainly is active in meristematic regions of the thallus. Previously a Trp-independent auxin synthesis pathway has been suggested to produce a majority of IAA in bryophytes. Our results indicate that the Trp-dependent IPyA pathway produces IAA that is essential for proper development of the gametophyte thallus of M. polymorpha. Furthermore, we show that dormancy of gemmae is positively regulated by auxin synthesized by the IPyA pathway in the apex of the thallus. Our results indicate that auxin synthesis, transport, and signaling, in addition to its role in growth and development, have a critical role in regulation of gemmae dormancy in M. polymorpha.
Dormancy is a key process allowing land plants to adapt to changing conditions in the terrestrial habitat, allowing the cessation of growth in response to environmental or physiological cues, entrance into a temporary quiescent state, and subsequent reactivation of growth in more favorable environmental conditions [1-3]. Dormancy may be induced seasonally, sporadically (e.g., in response to drought), or developmentally (e.g., seeds and apical dominance). Asexual propagules, known as gemmae, derived via clonal reproduction in bryophytes, are often dormant until displaced from the parent plant. In the liverwort Marchantia polymorpha, gemmae are produced within specialized receptacles, gemma cups, located on the dorsal side of the vegetative thallus [4]. Mature gemmae are detached from the parent plant but may remain in the cup, with gemma growth suppressed as long as the gemmae remain in the gemma cup and the parental plant is alive [5]. Following dispersal of gemmae from gemma cups by rain, the gemmae germinate in the presence of light and moisture, producing clonal offspring [6]. In land plants, the plant hormone abscisic acid (ABA) regulates many aspects of dormancy and water balance [7]. Here, we demonstrate that ABA plays a central role in the control of gemma dormancy as transgenic M. polymorpha gemmae with reduced sensitivity to ABA fail to establish and/or maintain dormancy. Thus, the common ancestor of land plants used the ABA signaling module to regulate germination of progeny in response to environmental cues, with both gemmae and seeds being derived structures co-opting an ancestral response system.
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.
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.
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 Pustnas, Sweden, and at Woburn, UK, offering the opportunity to assess the impact of genotype-by-environment interactions (GxE) 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 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 estimates and to considerable GxE interactions indicated by multivariate association analyses and modest trait accession correlations across sites (0.32-0.61).
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.
Growth and defence are plant traits that determine attack success by sessile herbivorous insects. Models of plant resistance suggest induced resistance to be an especially important mechanism in growth-dominated plant species. Little is known, however, about the effects of absence of induced responses as an alternative resistance mechanism. Large genetically based intraspecific variation in resistance (neonate larval mortality) occurs in Salix viminalis against the gall midge Dasineura marginemtorquens. In certain resistant genotypes, larval mortality is associated with formation of necrotic lesions and biochemical changes commonly found in hypersensitive response (HR). More detailed studies, however, have shown that resistance can occur without HR symptoms. By means of a quantitative trait locus (QTL) analysis, this study tested whether (1) resistance is caused by the HR or (2) by the failure of the midge to initiate a gall (non-responsiveness). One QTL explained 68 and 18 % of the variation in resistance in the experimental population during 2 years of investigation. Co-location of QTLs for resistance and HR was found in only one of the years, i.e. when HR was suppressed. The data suggest that the mechanism of resistance involves non-responsiveness, implying that the insect starves to death due to lack of resources. Such symptomless defence may be more important than previously thought because it is difficult to observe in nature.
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.
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 partlyconserved compared to Arabidopsis homologs, however further studies are needed to reveal the protein function of PaPRR1. Our data suggest that components of thethree-loop network typical of the circadian clock in angiosperms were present beforethe split of gymnosperms and angiosperms.
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.
FLOWERING LOCUS C (FLC) is one of the main genes influencing the vernalization requirement and natural flowering time variation in the annual Arabidopsis thaliana. Here we studied the effects of vernalization on flowering and its genetic basis in the perennial Arabidopsis lyrata. Two tandemly duplicated FLC genes (FLC1 and FLC2) were compared with respect to expression and DNA sequence. The effect of vernalization on flowering and on the expression of FLC1 was studied in three European populations. The genetic basis of the FLC1 expression difference between two of the populations was further studied by expression quantitative trait locus (eQTL) mapping and sequence analysis. FLC1 was shown to have a likely role in the vernalization requirement for flowering in A. lyrata. Vernalization decreased its expression and the northern study populations showed higher FLC1 expression than the southern one. eQTL mapping between two of the populations revealed one eQTL affecting FLC1 expression in the genomic region containing the FLC genes. Most FLC1 sequence differences between the study populations were found in the promoter region and in the first intron. Variation in the FLC1 sequence may cause differences in FLC1 expression between late- and early-flowering A. lyrata populations.
Population structure is a potential problem when testing for adaptive phenotypic differences among populations. The observed phenotypic differences among populations can simply be due to genetic drift, and if the genetic distance between them is not considered, the differentiation may be falsely interpreted as adaptive. Conversely, adaptive and demographic processes might have been tightly associated and correcting for the population structure may lead to false negatives. Here, we evaluated this problem in the cosmopolitan weed Capsella bursa-pastoris. We used RNA-Seq to analyse gene expression differences among 24 accessions, which belonged to a much larger group that had been previously characterized for flowering time and circadian rhythm and were genotyped using genotyping-by-sequencing (GBS) technique. We found that clustering of accessions for gene expression retrieved the same three clusters that were obtained with GBS data previously, namely Europe, the Middle East and Asia. Moreover, the three groups were also differentiated for both flowering time and circadian rhythm variation. Correction for population genetic structure when analysing differential gene expression analysis removed all differences among the three groups. This may suggest that most differences are neutral and simply reflect population history. However, geographical variation in flowering time and circadian rhythm indicated that the distribution of adaptive traits might be confounded by population structure. To bypass this confounding effect, we compared gene expression differentiation between flowering ecotypes within the genetic groups. Among the differentially expressed genes, FLOWERING LOCUS C was the strongest candidate for local adaptation in regulation of flowering time.
Allopolyploidy is generally perceived as a major source of evolutionary novelties and as an instantaneous way to create isolation barriers. However, we do not have a clear understanding of how two subgenomes evolve and interact once they have fused in an allopolyploid species nor how isolated they are from their relatives. Here, we address these questions by analyzing genomic and transcriptomic data of allotetraploid Capsella bursa-pastoris in three differentiated populations, Asia, Europe, and the Middle East. We phased the two subgenomes, one descended from the outcrossing and highly diverse Capsella grandiflora (Cbp(Cg)) and the other one from the selfing and genetically depauperate Capsella orientalis (Cbp(Co)). For each subgenome, we assessed its relationship with the diploid relatives, temporal changes of effective population size (N-e), signatures of positive and negative selection, and gene expression patterns. In all three regions, N-e of the two subgenomes decreased gradually over time and the Cbp(Co) subgenome accumulated more deleterious changes than Cbp(Cg). There were signs of widespread admixture between C. bursa-pastoris and its diploid relatives. The two subgenomes were impacted differentially depending on geographic region suggesting either strong interploidy gene flow or multiple origins of C. bursa-pastoris. Selective sweeps were more common on the Cbp(Cg) subgenome in Europe and the Middle East, and on the Cbp(Co) subgenome in Asia. In contrast, differences in expression were limited with the Cbp(Cg) subgenome slightly more expressed than Cbp(Co) in Europe and the Middle-East. In summary, after more than 100,000 generations of co-existence, the two subgenomes of C. bursa-pastoris still retained a strong signature of parental legacy but their evolutionary trajectory strongly varied across geographic regions. Author summary Allopolyploid species have two or more sets of chromosomes that originate from hybridization of different species. It remains largely unknown how the two genomes evolve in the same organism and how strongly their evolutionary trajectory depends on the initial differences between the two parental species and the specific demographic history of the newly formed allopolyploid species. To address these questions, we analyzed the genomic and gene expression variation of the shepherd's purse, a recent allopolyploid species, in three regions of its natural range. After approximate to 100,000 generations of co-existence within the same species, the two subgenomes had still retained part of the initial difference between the two parental species in the number of deleterious mutations reflecting a history of mating system differences. This difference, as well as differences in patterns of positive selection and levels of gene expression, also strongly depended on the specific histories of the three regions considered. Most strikingly, and unexpectedly, the allopolyploid species showed signs of hybridization with different diploid relatives or multiple origins in different parts of its range. Regardless if it was hybridization or multiple origins, this profoundly altered the relationship between the two subgenomes in different regions. Hence, our study illustrates how both the genomic structure and ecological arena interact to determine the evolutionary trajectories of allopolyploid species.
Allopolyploidy is generally perceived as a major source of evolutionary novelties and as an instantaneous way to create isolation barriers. However, we do not have a clear understanding of how two subgenomes evolve and interact once they have fused in an allopolyploid species and how isolated they are from their relatives. Here, we address these questions by analyzing genomic and transcriptomic data of allotetraploid Capsella bursa-pastoris in three differentiated populations, Asia, Europe and the Middle East. We phased the two subgenomes, one descended from the outcrossing and highly diverse Capsella grandiflora (Cg) and the other one from the selfing and genetically depauperate Capsella orientalis (Co). For each subgenome, we assessed its relationship with the diploid relatives, temporal change of effective population size Ne, signatures of positive and negative selection, and gene expression patterns. Introgression between C. bursa-pastoris and its diploid relatives was widespread and the two subgenomes were impacted differentially depending on geographic region. In all three regions, Ne of the two subgenomes decreased gradually and the Co subgenome accumulated more deleterious changes than Cg. Selective sweeps were more common on the Cg subgenome in Europe and the Middle East, and on the Co subgenome in Asia. In contrast, differences in expression were limited with the Cg subgenome slightly more expressed than Co in Europe and the Middle-East. In summary, after more than 100,000 generations of co-existence, the two subgenomes of C. bursa-pastoris still retained a strong signature of parental legacy and were differentially affected by introgression and selection.
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.
Small RNA (sRNA) including miRNA and siRNA are important in the regulation of diverse biological processes. Comparative studies of sRNA from plants have mainly focused on miRNA even though this group in most cases is just a mere fraction of the total sRNA diversity. In the present study we report results from an in-depth analysis of the sRNA population from the conifer Picea abies and compared the results to those of a range of species from the plant kingdom. The vast majority of sRNA in P. abies was 21 nucleotides long siRNA, of which a large fraction originate from degradation of transcribed sequences related to TIR-NBS-LRR (Toll/Interleukin-1-Nucleotide Binding Site-Leucine Rich Repeats) type resistance genes. Over 90% of all genes predicted to contain either a TIR or an NBS domain showed evidence of siRNA degradation. Data further suggests that those phased siRNA are initiated from miRNA guided cleavage, often by an abundant 22 nt miRNA. Comparative analysis over a range of divergent plant species revealed a large variation between species in the abundance of this phenomenon. The process seemed to be virtually absent in several species, including Arabidopsis thaliana, Oryza sativa and non-vascular plants, while particularly high frequencies were observed in Vitis vinifera and Populus thrichocarpa. This divergent pattern between species might reflect a mechanism to limit runaway transcription of these genes in species with rapidly expanding NBS-LRR gene families. Alternatively it might reflect variation in a counter-counter defence mechanism partly affected by differences in life history traits, e.g. perennial versus annual life cycles. A major difference between annuals and perennials, affecting the evolution of resistance, is that perennials will with almost certainty encounter many different pathogens before reproduction, and that the long generation times confers problems in matching the evolutionary rates of the pathogens.
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
Previous studies in the liverwort Marchantia polymorpha have shown that the putative evening complex (EC) genes LUX ARRHYTHMO (LUX) and ELF4-LIKE (EFL) have a function in the liverwort circadian clock. Here, we studied the growth phenotypes of MpLUX and MpEFL loss-of-function mutants, to establish if PHYTOCHROME-INTERACTING FACTOR (PIF) and auxin act downstream of the M. polymorpha EC in a growth-related pathway similar to the one described for the flowering plant Arabidopsis. We examined growth rates and cell properties of loss-of-function mutants, analyzed protein-protein interactions and performed gene expression studies using reporter genes. Obtained data indicate that an EC can form in M. polymorpha and that this EC regulates growth of the thallus. Altered auxin levels in Mplux mutants could explain some of the phenotypes related to an increased thallus surface area. However, because MpPIF is not regulated by the EC, and because Mppif mutants do not show reduced growth, the growth phenotype of EC-mutants is likely not mediated via MpPIF. In Arabidopsis, the circadian clock regulates elongation growth via PIF and auxin, but this is likely not an evolutionarily conserved growth mechanism in land plants. Previous inventories of orthologs to Arabidopsis clock genes in various plant lineages showed that there is high levels of structural differences between clocks of different plant lineages. Here, we conclude that there is also variation in the output pathways used by the different plant clocks to control growth and development.
Previous studies of plant circadian clock evolution have often relied on clock models and genes defined in Arabidopsis. These studies identified homologues with seemingly conserved function, as well as frequent gene loss. In the present study, we aimed to identify candidate clock genes in the liverwort Marchantia polymorpha using a more unbiased approach.
To identify genes with circadian rhythm we sequenced the transcriptomes of gemmalings in a time series in constant light conditions. Subsequently, we performed functional studies using loss-of-function mutants and gene expression reporters.
Among the genes displaying circadian rhythm, a homologue to the transcriptional co-repressor Arabidopsis DE-ETIOLATED1 showed high amplitude and morning phase. Because AtDET1 is arrhythmic and associated with the morning gene function of AtCCA1/LHY, that lack a homologue in liverworts, we functionally studied DET1 in M. polymorpha.
We found that the circadian rhythm of MpDET1 expression is disrupted in loss-of-function mutants of core clock genes and putative evening-complex genes. MpDET1 knock-down in turn results in altered circadian rhythm of nyctinastic thallus movement and clock gene expression. We could not detect any effect of MpDET1 knock-down on circadian response to light, suggesting that MpDET1 has a yet unknown function in the M. polymorpha circadian clock.
The circadian clock coordinates an organism’s growth, development and physiology with environmental factors. One illuminating example is the rhythmic growth of hypocotyls and cotyledons in Arabidopsis thaliana. Such daily oscillations in leaf position are often referred to as sleep movements or nyctinasty. Here, we report that plantlets of the liverwort Marchantia polymorpha show analogous rhythmic movements of thallus lobes, and that the circadian clock controls this rhythm, with auxin a likely output pathway affecting these movements. The mechanisms of this circadian clock are partly conserved as compared to angiosperms, with homologs to the core clock genes PRR, RVE and TOC1 forming a core transcriptional feedback loop also in M. polymorpha.
A Brassica nigra genetic linkage map was developed from a highly polymorphic cross analyzed with a set of low copy number Brassica RFLP probes. The Brassica genome is extensively duplicated with eight distinct sets of chromosomal segments, each present in
Plants commonly referred to as "bryophytes" belong to three major lineages of non-vascular plants: the liverworts, the hornworts and the mosses. They are unique among land plants in having a dominant haploid generation and a short-lived diploid sporophytic generation. The dynamics of selection acting on a haploid genome differs from those acting on a diploid genome: new mutations are directly exposed to selection. The general aim of this paper is to investigate the diversification rate of bryophytes - measured as silent site substitution rate representing neutral evolution (mutation rate) and the nonsynonymous to synonymous substitution rate ratio (dN/dS) representing selective evolution - and compare it with earlier studies on vascular plants. Results show that the silent site substitution rate is lower for liverworts as compared to angiosperms, but not as low as for gymnosperms. The selection pressure, measured as dN/dS, is not remarkably lower for bryophytes as compared to other diploid dominant plants as would be expected by the masking hypothesis, indicating that other factors are more important than ploidy.
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.
The potential role of introgression in evolution has gained increased interest in recent years. Although some fascinating examples have been reported, more information is needed to generalize the importance of hybridization and introgression for adaptive divergence. As limited data exist on haploid dominant species, we analyzed genomes of three subspecies of the liverwortMarchantia polymorpha. We used available genomic data for subsp.ruderalisand carried out whole-genome (PacBio) sequencing for one individual each of subsp.montivagansand subsp.polymorphaas well as Illumina resequencing of additional genomes for all three subspecies. The three subspecies were compared againstM. paleaceaas outgroup. Our analyses revealed separation of the three taxa, but all three possible topologies were richly represented across the genomes, and the underlying divergence order less obvious. This uncertainty could be the result of the divergence of the three subspecies close in time, or that introgression has been frequent since divergence. In particular, we found that pseudo-chromosome 2 in subsp.montivaganswas much more diverged than other parts of the genomes. This could either be explained by specific capture of chromosome 2 from an unknown related species through hybridization or by conservation of chromosome 2 despite intermittent or ongoing introgression affecting more permeable parts of the genomes. A higher degree of chromosomal rearrangements on pseudo-chromosome 2 support the second hypothesis. Species tree analyses recovered an overall topology where subsp.montivagansdiverged first and subsp.ruderalisand subsp.polymorphaappeared as sister lineages. Each subspecies was associated with its own chloroplast and mitochondrial haplotype group. Our data suggest introgression but refute a previous hypothesis that subsp.ruderalisis a new stabilized hybrid between the other two subspecies.
Why do some genomes stay small and simple, while others become huge, and why are some genomes more stable? In contrast to angiosperms and gymnosperms, liverworts are characterized by small genomes with low variation in size and conserved chromosome numbers. We quantified genome evolution among five Marchantiophyta (liverworts), measuring gene characteristics, transposable element (TE) landscape, collinearity, and sex chromosome evolution that might explain the small size and limited variability of liverwort genomes. No genome duplications were identified among examined liverworts and levels of duplicated genes are low. Among the liverwort species, Lunularia cruciata stands out with a genome size almost twice that of the other liverwort species investigated here, and most of this increased size is due to bursts of Ty3/Gypsy retrotransposons. Intrachromosomal rearrangements between examined liverworts are abundant but occur at a slower rate compared with angiosperms. Most genes on L. cruciata scaffolds have their orthologs on homologous Marchantia polymorpha chromosomes, indicating a low degree of rearrangements between chromosomes. Still, translocation of a fragment of the female U chromosome to an autosome was predicted from our data, which might explain the uniquely small U chromosome in L. cruciata. Low levels of gene duplication, TE activity, and chromosomal rearrangements might contribute to the apparent slow rate of morphological evolution in liverworts.
Plants are the foundation of terrestrial ecosystems, and their colonization of land was probably facilitated by mutualistic associations with arbuscular mycorrhizal fungi. Following this founding event, plant diversification has led to the emergence of a tremendous diversity of mutualistic symbioses with microorganisms, ranging from extracellular associations to the most intimate intracellular associations, where fungal or bacterial symbionts are hosted inside plant cells. Here, through analysis of 271 transcriptomes and 116 plant genomes spanning the entire land-plant diversity, we demonstrate that a common symbiosis signalling pathway co-evolved with intracellular endosymbioses, from the ancestral arbuscular mycorrhiza to the more recent ericoid and orchid mycorrhizae in angiosperms and ericoid-like associations of bryophytes. By contrast, species forming exclusively extracellular symbioses, such as ectomycorrhizae, and those forming associations with cyanobacteria, have lost this signalling pathway. This work unifies intracellular symbioses, revealing conservation in their evolution across 450 million yr of plant diversification. An extensive phylogenomics study based on hundreds of genomes and transcriptomes provides a new interpretation of the evolution of different types of symbiotic associations in land plants, and reveals a conserved ancestral symbiosis pathway.
Although land plant germ cells have received much attention, knowledge about their specification is still limited. We thus identified transcripts enriched in egg cells of the bryophyte model species Physcomitrium patens, compared the results with angiosperm egg cells, and selected important candidate genes for functional analysis. We used laser-assisted microdissection to perform a cell-type-specific transcriptome analysis on egg cells for comparison with available expression profiles of vegetative tissues and male reproductive organs. We made reporter lines and knockout mutants of the two BONOBO (PbBNB) genes and studied their role in reproduction. We observed an overlap in gene activity between bryophyte and angiosperm egg cells, but also clear differences. Strikingly, several processes that are male-germline specific in Arabidopsis are active in the P. patens egg cell. Among those were the moss PbBNB genes, which control proliferation and identity of both female and male germlines. Pathways shared between male and female germlines were most likely present in the common ancestors of land plants, besides sex-specifying factors. A set of genes may also be involved in the switches between the diploid and haploid moss generations. Nonangiosperm gene networks also contribute to the specification of the P. patens egg cell.