Adaptive evolution of the Populus tremula photoperiod pathway

Genomic Adaptive Evolution of Sand Rice (Agriophyllum squarrosum) and Its Implications for Desert Ecosystem Restoration

Natural selection is a significant driver of population divergence and speciation of plants. Due to local adaptation to geographic regions with ecological gradients, plant populations harbored a wide range of adaptive genetic variation to enable them to survive the heterogeneous habitats. This is all the more necessary for desert plants, as they must tolerant more striking gradients of abiotic stresses. However, the genomic mechanism by which desert plants adapt to ecological heterogeneity remains unclear, which could help to guide the sustainability of desert ecosystems. Here, using restriction-site-associated DNA sequencing in 38 natural populations, we investigated the genomic divergence and environmental adaptation of sand rice, Agriophyllum squarrosum, an annual pioneer species that covers sand dunes in northern China. Population genetic structure analyses showed that sand rice could be divided into three geographically distinct lineages, namely, Northwest, Central, and East. Phylogeographic analyses revealed that the plant might originate locally in Bergen County and further differentiated into the East lineage and then the Central lineage. Ecological niche modeling found that different lineages occupied distinct ecological niches, suggesting that the ecological gradient would have triggered genomic differentiation among sand rice lineages. Ecological association study supported that the three SNPs under divergent selection were closely correlated with precipitation gradients, indicating that precipitation might be the most important stress trigger for lineage diversity in sand rice. These adaptive SNPs could be used to genotype suitable germplasms for the ecological restoration of specific desertified lands. Further analyses found that genetic structure could significantly overestimate the signals for balancing selection. Within the Central lineage, we still found that 175 SNPs could be subject to balancing selection, which could be the means by which sand rice maintains genetic diversity and adapts to multiple stresses across heterogeneous deserts and sandy lands. From a genomic point of view, this study highlighted the local and global adaptation patterns of a desert plant to extreme and heterogeneous habitats. Our data provide molecular guidance for the restoration of desertified lands in the arid and semi-arid regions of China and could facilitate the marker assistant breeding of this potential crop to mitigate climate change.

Pathway position constrains the evolution of an ecologically important pathway in aspens (Populus tremula L.)

Many ecological interactions of aspens and their relatives (Populus spp.) are affected by products of the phenylpropanoid pathway synthesizing condensed tannins (CTs), whose production involves trade-offs with other ecologically important compounds and with growth. Genes of this pathway are candidates for investigating the role of selection on ecologically important, polygenic traits. We analysed sequences from 25 genes representing 10 steps of the CT synthesis pathway, which produces CTs used in defence and lignins used for growth, in 12 individuals of European aspen (Populus tremula). We compared these to homologs from P. trichocarpa, to a control set of 77 P. tremula genes, to genome-wide resequencing data and to RNA-seq expression levels, in order to identify signatures of selection distinct from those of demography. In Populus, pathway position exerts a strong influence on the evolution of these genes. Nonsynonymous diversity, divergence and allele frequency shifts (Tajima's D) were much lower than for synonymous measures. Expression levels were higher, and the direction of selection more negative, for upstream genes than for those downstream. Selective constraints act with increasing intensity on upstream genes, despite the presence of multiple paralogs in most gene families. Pleiotropy, expression level, flux control and codon bias appear to interact in determining levels and patterns of variation in genes of this pathway, whose products mediate a wide array of ecological interactions for this widely distributed species.

Selection for population-specific adaptation shaped patterns of variation in the photoperiod pathway genes in Arabidopsis lyrata during post-glacial colonization

Spatially varying selection can lead to population-specific adaptation, which is often recognized at the phenotypic level; however, the genetic evidence is weaker in many groups of organisms. In plants, environmental shifts that occur due to colonization of a novel environment may require adaptive changes in the timing of growth and flowering, which are often governed by location-specific environmental cues such as day length. We studied locally varying selection in 19 flowering time loci in nine populations of the perennial herb Arabidopsis lyrata, which has a wide but patchy distribution in temperate and boreal regions of the northern hemisphere. The populations differ in their recent population demographic and colonization histories and current environmental conditions, especially in the growing season length. We searched for population-specific molecular signatures of directional selection by comparing a set of candidate flowering time loci with a genomic reference set within each population using multiple approaches and contrasted the patterns of different populations. The candidate loci possessed approximately 20% of the diversity of the reference loci. On average the flowering time loci had more rare alleles (a smaller Tajima's D) and an excess of highly differentiated sites relative to the reference, suggesting positive selection. The strongest signal of selection was detected in photoperiodic pathway loci in the colonizing populations of Northwestern Europe, whereas no evidence of positive selection was detected in the Central European populations. These findings emphasized the population-specific nature of selection and suggested that photoperiodic adaptation was important during postglacial colonization of the species.

Using archaeogenomic and computational approaches to unravel the history of local adaptation in crops

Our understanding of the evolution of domestication has changed radically in the past 10 years, from a relatively simplistic rapid origin scenario to a protracted complex process in which plants adapted to the human environment. The adaptation of plants continued as the human environment changed with the expansion of agriculture from its centres of origin. Using archaeogenomics and computational models, we can observe genome evolution directly and understand how plants adapted to the human environment and the regional conditions to which agriculture expanded. We have applied various archaeogenomics approaches as exemplars to study local adaptation of barley to drought resistance at Qasr Ibrim, Egypt. We show the utility of DNA capture, ancient RNA, methylation patterns and DNA from charred remains of archaeobotanical samples from low latitudes where preservation conditions restrict ancient DNA research to within a Holocene timescale. The genomic level of analyses that is now possible, and the complexity of the evolutionary process of local adaptation means that plant studies are set to move to the genome level, and account for the interaction of genes under selection in systems-level approaches. This way we can understand how plants adapted during the expansion of agriculture across many latitudes with rapidity.

Archaeogenomic insights into the adaptation of plants to the human environment: pushing plant-hominin co-evolution back to the Pliocene

The colonization of the human environment by plants, and the consequent evolution of domesticated forms is increasingly being viewed as a co-evolutionary plant-human process that occurred over a long time period, with evidence for the co-evolutionary relationship between plants and humans reaching ever deeper into the hominin past. This developing view is characterized by a change in emphasis on the drivers of evolution in the case of plants. Rather than individual species being passive recipients of artificial selection pressures and ultimately becoming domesticates, entire plant communities adapted to the human environment. This evolutionary scenario leads to systems level genetic expectations from models that can be explored through ancient DNA and Next Generation Sequencing approaches. Emerging evidence suggests that domesticated genomes fit well with these expectations, with periods of stable complex evolution characterized by large amounts of change associated with relatively small selective value, punctuated by periods in which changes in one-half of the plant-hominin relationship cause rapid, low-complexity adaptation in the other. A corollary of a single plant-hominin co-evolutionary process is that clues about the initiation of the domestication process may well lie deep within the hominin lineage.

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra- and interspecific genetic variation at 25 climate-related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two-thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. FST outlier analysis identified six SNPs in five climate-related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation-with-migration model indicated that gene flow between the species has been more restricted for climate-related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species-specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.

Exome resequencing reveals signatures of demographic and adaptive processes across the genome and range of black cottonwood (Populus trichocarpa)

Extant variation in temperate and boreal plant species has been influenced by both demographic histories associated with Pleistocene glacial cycles and adaptation to local climate. We used sequence capture to investigate the role of these neutral and adaptive processes in shaping diversity in black cottonwood (Populus trichocarpa). Nucleotide diversity and Tajima's D were lowest at replacement sites and highest at intergenic sites, while LD showed the opposite pattern. With samples grouped into three populations arrayed latitudinally, effective population size was highest in the north, followed by south and centre, and LD was highest in the south followed by the north and centre, suggesting a possible northern glacial refuge. FST outlier analysis revealed that promoter, 5'-UTR and intronic sites were enriched for outliers compared with coding regions, while no outliers were found among intergenic sites. Codon usage bias was evident, and genes with synonymous outliers had 30% higher average expression compared with genes containing replacement outliers. These results suggest divergent selection related to regulation of gene expression is important to local adaptation in P. trichocarpa. Finally, within-population selective sweeps were much more pronounced in the central population than in putative northern and southern refugia, which may reflect the different demographic histories of the populations and concomitant effects on signatures of genetic hitchhiking from standing variation.

Clinal variation at phenology-related genes in spruce: parallel evolution in FTL2 and Gigantea?

Parallel clines in different species, or in different geographical regions of the same species, are an important source of information on the genetic basis of local adaptation. We recently detected latitudinal clines in SNPs frequencies and gene expression of candidate genes for growth cessation in Scandinavian populations of Norway spruce (Picea abies). Here we test whether the same clines are also present in Siberian spruce (P. obovata), a close relative of Norway spruce with a different Quaternary history. We sequenced nine candidate genes and 27 control loci and genotyped 14 SSR loci in six populations of P. obovata located along the Yenisei river from latitude 56°N to latitude 67°N. In contrast to Scandinavian Norway spruce that both departs from the standard neutral model (SNM) and shows a clear population structure, Siberian spruce populations along the Yenisei do not depart from the SNM and are genetically unstructured. Nonetheless, as in Norway spruce, growth cessation is significantly clinal. Polymorphisms in photoperiodic (FTL2) and circadian clock (Gigantea, GI, PRR3) genes also show significant clinal variation and/or evidence of local selection. In GI, one of the variants is the same as in Norway spruce. Finally, a strong cline in gene expression is observed for FTL2, but not for GI. These results, together with recent physiological studies, confirm the key role played by FTL2 and circadian clock genes in the control of growth cessation in spruce species and suggest the presence of parallel adaptation in these two species.

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

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

The perennial clock is an essential timer for seasonal growth events and cold hardiness

Over the last several decades, changes in global temperatures have led to changes in local environments affecting the growth conditions for many species. This is a trend that makes it even more important to understand how plants respond to local variations and seasonal changes in climate. To detect daily and seasonal changes as well as acute stress factors such as cold and drought, plants rely on a circadian clock. This chapter introduces the current knowledge and literature about the setup and function of the circadian clock in various tree and perennial species, with a focus on the Populus genus.

A recent local sweep at the PHYA locus in the Northern European Spiterstulen population of Arabidopsis lyrata

Northern and central European Arabidopsis lyrata ssp. petraea populations are locally adapted to prevailing climatic conditions through differences in timing of life history events. The timing of flowering and, in perennials, the timing of growth cessation influence fitness. Phytochrome A may have an important role in regulating these life history traits as it perceives changes in daylength. We asked whether PHYA has contributed to local adaptation to the northern conditions in A. l. petraea. To search for signals of directional selection at the PHYA locus, we resequenced PHYA and 9 short fragments around PHYA from a 57-kb region from a German (Plech) and a Norwegian (Spiterstulen) population and compared patterns of differentiation and diversity to a set of 19 reference loci around the genome. First, we found that the populations were highly differentiated: there were three nonsynonymous fixed differences at the PHYA locus, which was in stark contrast with the total four fixed differences in the 19 reference loci. Compatible with a sweep hypothesis, variation was almost completely removed from the 9.4-kb region around PHYA in the northern Spiterstulen population. The overall level of linkage disequilibrium (LD) was higher in Spiterstulen, but there was no LD across the PHYA locus in the population, which is also a known consequence of a selective sweep. The sweep has likely occurred after the last glacial maximum, which suggests that it has contributed to adaptation to the northern conditions.

Revisiting the sequencing of the first tree genome: Populus trichocarpa

Ten years ago, it was announced that the Joint Genome Institute with funds provided by the Department of Energy, Office of Science, Biological and Environmental Research would sequence the black cottonwood (Populus trichocarpa Torr. & Gray) genome. This landmark decision was the culmination of work by the forest science community to develop Populus as a model system. Since its public release in late 2006, the availability of the Populus genome has spawned research in plant biology, morphology, genetics and ecology. Here we address how the tree physiologist has used this resource. More specifically, we revisit our earlier contention that the rewards of sequencing the Populus genome would depend on how quickly scientists working with woody perennials could adopt molecular approaches to investigate the mechanistic underpinnings of basic physiological processes. Several examples illustrate the integration of functional and comparative genomics into the forest sciences, especially in areas that target improved understanding of the developmental differences between woody perennials and herbaceous annuals (e.g., phase transitions). Sequencing the Populus genome and the availability of genetic and genomic resources has also been instrumental in identifying candidate genes that underlie physiological and morphological traits of interest. Genome-enabled research has advanced our understanding of how phenotype and genotype are related and provided insights into the genetic mechanisms whereby woody perennials adapt to environmental stress. In the future, we anticipate that low-cost, high-throughput sequencing will continue to facilitate research in tree physiology and enhance our understanding at scales of individual organisms and populations. A challenge remains, however, as to how genomic resources, including the Populus genome, can be used to understand ecosystem function. Although examples are limited, progress in this area is encouraging and will undoubtedly improve as future research targets the many unique aspects of Populus as a keystone species in terrestrial ecosystems.

The dynamic nature of bud dormancy in trees: environmental control and molecular mechanisms

In tree species native to temperate and boreal regions, the activity-dormancy cycle is an important adaptive trait both for survival and growth. We discuss recent research on mechanisms controlling the overlapping developmental processes that define the activity-dormancy cycle, including cessation of apical growth, bud development, induction, maintenance and release of dormancy, and bud burst. The cycle involves an extensive reconfiguration of metabolism. Environmental control of the activity-dormancy cycle is based on perception of photoperiodic and temperature signals, reflecting adaptation to prevailing climatic conditions. Several molecular actors for control of growth cessation have been identified, with the CO/FT regulatory network and circadian clock having important coordinating roles in control of growth and dormancy. Other candidate regulators of bud set, dormancy and bud burst have been identified, such as dormancy-associated MADS-box factors, but their exact roles remain to be discovered. Epigenetic mechanisms also appear to factor in control of the activity-dormancy cycle. Despite evidence for gibberellins as negative regulators in growth cessation, and ABA and ethylene in bud formation, understanding of the roles that plant growth regulators play in controlling the activity-dormancy cycle is still very fragmentary. Finally, some of the challenges for further research in bud dormancy are discussed.

Extreme conservation and non-neutral evolution of the cpmA Circadian locus in a globally distributed Chroococcidiopsis sp. from naturally stressful habitats

Cyanobacteria are among the most ancient organisms known to have circadian rhythms. The cpmA gene is involved in controlling the circadian output signal. We studied polymorphism and divergence of this gene in six populations of a stress-tolerant cyanobacterium, Chroococcidiopsis sp., sampled in extreme habitats across the globe. Despite high haplotype diversity (0.774), nucleotide diversity of cpmA is very low (π = 0.0034): the gene appears to be even more conserved than housekeeping genes. Even though the populations were sampled thousands kilometers apart, they manifested virtually no genetic differentiation at this locus (F(ST) = 0.0228). Using various tests for neutrality, we determined that evolution of cpmA significantly departures from the neutral model and is governed by episodic positive selection.

Disentangling the roles of history and local selection in shaping clinal variation of allele frequencies and gene expression in Norway spruce (Picea abies)

Understanding the genetic basis of local adaptation is challenging due to the subtle balance among conflicting evolutionary forces that are involved in its establishment and maintenance. One system with which to tease apart these difficulties is clines in adaptive characters. Here we analyzed genetic and phenotypic variation in bud set, a highly heritable and adaptive trait, among 18 populations of Norway spruce (Picea abies), arrayed along a latitudinal gradient ranging from 47°N to 68°N. We confirmed that variation in bud set is strongly clinal, using a subset of five populations. Genotypes for 137 single-nucleotide polymorphisms (SNPs) chosen from 18 candidate genes putatively affecting bud set and 308 control SNPs chosen from 264 random genes were analyzed for patterns of genetic structure and correlation to environment. Population genetic structure was low (F(ST) = 0.05), but latitudinal patterns were apparent among Scandinavian populations. Hence, part of the observed clinal variation should be attributable to population demography. Conditional on patterns of genetic structure, there was enrichment of SNPs within candidate genes for correlations with latitude. Twenty-nine SNPs were also outliers with respect to F(ST). The enrichment for clinal variation at SNPs within candidate genes (i.e., SNPs in PaGI, PaPhyP, PaPhyN, PaPRR7, and PaFTL2) indicated that local selection in the 18 populations, and/or selection in the ancestral populations from which they were recently derived, shaped the observed cline. Validation of these genes using expression studies also revealed that PaFTL2 expression is significantly associated with latitude, thereby confirming the central role played by this gene in the control of phenology in plants.

Evolutionary history of pearl millet (Pennisetum glaucum [L.] R. Br.) and selection on flowering genes since its domestication

The plant domestication process is associated with considerable modifications of plant phenotype. The identification of the genetic basis of this adaptation is of great interest for evolutionary biology. One of the methods used to identify such genes is the detection of signatures of selection. However, domestication is generally associated with major demographic effects. It is therefore crucial to disentangle the effects of demography and selection on diversity. In this study, we investigated selection in a flowering time pathway during domestication of pearl millet. We first used a random set of 20 genes to model pearl millet domestication using approximate Bayesian computation. This analysis showed that a model with exponential growth and wild-cultivated gene flow was well supported by our data set. Under this model, the domestication date of pearl millet is estimated at around 4,800 years ago. We assessed selection in 15 pearl millet DNA sequences homologous to flowering time genes and showed that these genes underwent selection more frequently than expected. We highlighted significant signatures of selection in six pearl millet flowering time genes associated with domestication or improvement of pearl millet. Moreover, higher deviations from neutrality were found for circadian clock-associated genes. Our study provides new insights into the domestication process of pearl millet and shows that a category of genes of the flowering pathway were preferentially selected during pearl millet domestication.

Local selection across a latitudinal gradient shapes nucleotide diversity in balsam poplar, Populus balsamifera L

Molecular studies of adaptive evolution often focus on detecting selective sweeps driven by positive selection on a species-wide scale; however, much adaptation is local, particularly of ecologically important traits. Here, we look for evidence of range-wide and local adaptation at candidate genes for adaptive phenology in balsam poplar, Populus balsamifera, a widespread forest tree whose range extends across environmental gradients of photoperiod and growing season length. We examined nucleotide diversity of 27 poplar homologs of the flowering-time network-a group of genes that control plant developmental phenology through interactions with environmental cues such as photoperiod and temperature. Only one gene, ZTL2, showed evidence of reduced diversity and an excess of fixed replacement sites, consistent with a species-wide selective sweep. Two other genes, LFY and FRI, harbored high levels of nucleotide diversity and exhibited elevated differentiation between northern and southern accessions, suggesting local adaptation along a latitudinal gradient. Interestingly, FRI has also been identified as a target of local selection between northern and southern accessions of Arabidopsis thaliana, indicating that this gene may be commonly involved in ecological adaptation in distantly related species. Our findings suggest an important role for local selection shaping molecular diversity and reveal limitations of inferring molecular adaptation from analyses designed only to detect species-wide selective sweeps.