A Genome Scan for Genes Underlying Microgeographic-Scale Local Adaptation in a Wild Arabidopsis Species

The genetic architecture of repeated local adaptation to climate in distantly related plants

Closely related species often use the same genes to adapt to similar environments. However, we know little about why such genes possess increased adaptive potential and whether this is conserved across deeper evolutionary lineages. Adaptation to climate presents a natural laboratory to test these ideas, as even distantly related species must contend with similar stresses. Here, we re-analyse genomic data from thousands of individuals from 25 plant species as diverged as lodgepole pine and Arabidopsis (~300 Myr). We test for genetic repeatability based on within-species associations between allele frequencies in genes and variation in 21 climate variables. Our results demonstrate significant statistical evidence for genetic repeatability across deep time that is not expected under randomness, identifying a suite of 108 gene families (orthogroups) and gene functions that repeatedly drive local adaptation to climate. This set includes many orthogroups with well-known functions in abiotic stress response. Using gene co-expression networks to quantify pleiotropy, we find that orthogroups with stronger evidence for repeatability exhibit greater network centrality and broader expression across tissues (higher pleiotropy), contrary to the 'cost of complexity' theory. These gene families may be important in helping wild and crop species cope with future climate change, representing important candidates for future study.

Genomics of ecological adaptation in Canary Island Descurainia (Brassicaceae) and comparisons with other Brassicaceae

Oceanic archipelagos provide striking examples of lineages that have radiated over pronounced ecological gradients. Accompanying this diversification, lineages have evolved adaptations allowing survival in extreme environments. Here, we investigate the genomic basis of ecological adaptation in Canary Island Descurainia (Brassicaceae), an island relative of Arabidopsis. The seven endemic species have diversified in situ along an elevational and ecological gradient, from low-elevation scrub to high-elevation sub-alpine desert. We first generated a reference genome for Descurainia millefolia, phylogenetic analysis of which placed it as sister to D. sophioides. Ninety-six gene families were found to be specific to D. millefolia and a further 1087 and 1469 gene families have expanded or contracted in size, respectively, along the D. millefolia branch. We then employed genome re-sequencing to sample 14 genomes across the seven species of Canary Island Descurainia and an outgroup. Phylogenomic analyses were consistent with previous reconstructions of Canary Island Descurainia in resolving low- and high-elevation clades. Using the branch-site dN/dS method, we detected positive selection for 275 genes on the branch separating the low- and high-elevation species and these positively selected genes (PSGs) were significantly enriched for functions related to reproduction and stress tolerance. Comparing PSGs to those in analyses of adaptation to elevation and/or latitude in other Brassicaceae, we found little evidence of widespread convergence and gene reuse, except for two examples, one of which was a significant overlap between Descurainia and Draba nivalis, a species restricted to high latitudes. The study of Canary Island Descurainia suggests that the transition to high-elevation environments such as that found in the high mountains of the Canary Islands involves selection on genes related to reproduction and stress tolerance but that repeated evolution across different lineages that have evolved into similar habitats is limited, indicating substantially different molecular trajectories to adaptation.

Integrative analysis of physiology and genomics provides insights into freeze tolerance adaptations of Acacia koa along an elevational cline

Natural selection for plant species in heterogeneous environments creates genetic variation for traits such as cold tolerance. While physiological or molecular analyses have been used to evaluate stress tolerance adaptations, combining these approaches may provide deeper insight. Acacia koa (koa) occurs from sea level to 2300 m in Hawai'i, USA. At high elevations, natural koa populations have declined due to deforestation, and freeze tolerance is a limiting factor for tree regeneration. We used physiology and molecular analyses to evaluate cold tolerance of koa populations from low (300-750 m), middle (750-1500 m), and high elevations (1500-2100 m). Half of the seedlings were cold acclimated by exposure to progressively lowered air temperatures for eight weeks (from 25.6/22.2°C to 8/4°C, day/night). Using the whole plant physiology-freezing test and koa C-repeat Binding Factor CBF genes, our results indicated that koa can be cold-acclimated when exposed to low, non-freezing temperatures. Seedlings from high elevations had consistently higher expression of Koa CBF genes associated with cold tolerance, helping to explain variation in cold-hardy phenotypes. Evaluation of the genetic background of 22 koa families across the elevations with low coverage RNA sequencing indicated that high elevation koa had relatively low values of heterozygosity, suggesting that adaptation is more likely to arise in the middle and low elevation sources. This physiology and molecular data for cold tolerance of koa across the elevation gradient of the Hawaiian Islands provides insights into natural selection processes and may help to support guidelines for conservation and seed transfer in forest restoration efforts.

Does selection occur at the intermediate zone of two insufficiently isolated populations? A whole-genome analysis along an altitudinal gradient

Adaptive divergence occurs even between insufficiently isolated populations when there is a great difference in environments between their habitats. Individuals present in an intermediate zone of the two divergent populations are expected to have an admixed genetic structure due to gene flow. A selective pressure that acts on the genetically admixed individuals may limit the gene flow and maintain the adaptive divergence. Here, we addressed a question whether selection occurs in the genetically admixed individuals between two divergent populations. Arabidopsis halleri is a perennial montane plant, which has clear phenotypic dimorphisms between highland and lowland habitats in Mt. Ibuki, central Japan. We obtained the whole-genome sequences of Arabidopsis halleri plants along an altitudinal gradient of 359-1,317 m with a high spatial resolution (mean altitudinal interval of 20 m). We found a zone where the highland and lowland genes were mixing (intermediate subpopulation). In the intermediate subpopulation, we identified 5 and 13 genome regions, which included 3 and 8 genes, that had a high frequency of alleles that are accumulated in highland and lowland subpopulations, respectively. In addition, we also found that the frequency of highland alleles of these selected genome regions was smaller in the lowland subpopulation compared with that of the non-selected regions. These results suggest that the selection in the intermediate subpopulation might limit the gene flow and contribute to the adaptive divergence between altitudes. We also identified 7 genome regions that had low heterozygote frequencies in the intermediate subpopulation. We conclude that different types of selection in addition to gene flow occur at the intermediate altitude and shape the genetic structure across altitudes.

Whole genome sequencing reveals footprints of adaptive genetic variation in populations of Eruca sativa

Populations of Eruca sativa (Brassicaceae) derived from arid and Mediterranean habitats exhibit ecotypic differentiation. Here, pooled DNA sequencing was used to assess adaptive genome differentiation in the two ecotypes. Differentiated SNP loci were scanned with the empirical FST outlier method and by correlating allele frequencies with environmental parameters. Genetic diversity values were relatively higher in the pooled arid genome, whereas the pooled Mediterranean genome exhibited stronger directional selection, indicating the impact of climatic conditions on genetic diversity. GO enrichment analysis categorized the annotated differentiated loci according to biological processes, revealing a large set of candidate genes related to abiotic and biotic stress responses. Allelic variation was detected in regulatory elements and coding regions (synonymous and non-synonymous mutations) of genes belonging to different transcription factors and phytohormone signaling, suggesting adaptation to both abiotic and biotic conditions. Furthermore, SNP mutations were also found in genic regions belonging to the synthesis of secondary metabolites, including aliphatic glucosinolates and their hydrolyzed bioactive compounds, among others. The results of this eco-genomic study demonstrate the role of divergent abiotic and biotic selection factors in evolutionary processes leading to adaptive ecotypic differentiation.

Multi-Omics Approaches Revealed the Associations of Host Metabolism and Gut Microbiome With Phylogeny and Environmental Adaptation in Mountain Dragons

The molecular basis enabling the adaptation of animals to spatially heterogeneous environments is a critical clue for understanding the variation, formation, and maintenance of biodiversity in the context of global climate change. Mountain dragons (Agamidae: Diploderma) thrive in the Hengduan Mountain Region, a biodiversity hotspot and a typical spatially heterogeneous environment. Here, we compare the liver and muscle metabolome and gut microbiome of 11 geographical populations from three Diploderma species (D. iadinum, D. yulongsense, and D. vela) after 7 days acclimation in the same laboratory conditions. Amino acid metabolism, particularly the products of the glutathione cycle, accounted for major interspecies variations, implying its significance in genetic differentiation among mountain dragons. Notably, the cold-dwelling D. vela and D. yulongense populations tended to have higher glycerophosphate, glycerol-3-phosphocholine, and kinetin levels in their liver, higher carnosine levels in their muscle, and higher Lachnospiraceae levels in their gut. Phylogeny, net primary productivity (NPP), and the temperature had the highest explanation rate to the variations in muscle metabolome, liver metabolome, and gut microbiome, respectively, suggesting heterogeneity of biological systems in response to climatic variations. Therefore, we suggested that the organ heterogeneity in environmental responsiveness might be substantial for mountain dragons to thrive in complicated environments.

A cross-scale approach to unravel the molecular basis of plant phenology in temperate and tropical climates

Plants have evolved to time their leafing, flowering and fruiting in appropriate seasons for growth, reproduction and resting. As a consequence of their adaptation to geographically different environments, there is a rich diversity in plant phenology from temperate and tropical climates. Recent progress in genetic and molecular studies will provide numerous opportunities to study the genetic basis of phenological traits and the history of adaptation of phenological traits to seasonal and aseasonal environments. Integrating molecular data with long-term phenology and climate data into predictive models will be a powerful tool to forecast future phenological changes in the face of global environmental change. Here, we review the cross-scale approach from genes to plant communities from three aspects: the latitudinal gradient of plant phenology at the community level, the environmental and genetic factors underlying the diversity of plant phenology, and an integrated approach to forecast future plant phenology based on genetically informed knowledge. Synthesizing the latest knowledge about plant phenology from molecular, ecological and mathematical perspectives will help us understand how natural selection can lead to the further evolution of the gene regulatory mechanisms in phenological traits in future forest ecosystems.

Commensal microbiota and host metabolic divergence are associated with the adaptation of Diploderma vela to spatially heterogeneous environments

Heterogeneous environment adaptation is critical to understand the species evolution and response to climate change. However, how narrow-range species adapt to micro-geographic heterogeneity has been overlooked, and there is a lack of insights from metabolism and commensal microbiota. Here, we studied the environmental adaptation for 3 geographic populations (>40 km apart) of Diploderma vela, a lizard endemic to dry-hot valleys of the Hengduan Mountain Region. The climatic boundary caused a cooler, droughtier, and barren environment for northernmost population (RM) than the middle (QZK) and southernmost populations (FS). Correspondingly, significant divergences in liver and muscle metabolism and commensal microbiota were detected between RM and QZK or FS individuals, but not between QZK and FS individuals. Phospholipid composition, coenzyme level (i.e. pyridoxal and NAD⁺ ), and cholesterol metabolism (e.g. androgen and estriol synthesis) constituted the major metabolic difference between RM and QZK/FS groups. FS and QZK individuals kept abundant Proteobacteria and antifungal strains, while RM individuals maintained more Firmicutes and Bacteroidota. Strong associations existed between varied host metabolite and gut microbes. How were these interpopulation variations associated to the environment adaptation were discussed. These results provided some novel insights into the environmental adaptation and implicated the consequence of climate change on narrow-range species.

Genomic basis of parallel adaptation varies with divergence in Arabidopsis and its relatives

Parallel adaptation provides valuable insight into the predictability of evolutionary change through replicated natural experiments. A steadily increasing number of studies have demonstrated genomic parallelism, yet the magnitude of this parallelism varies depending on whether populations, species, or genera are compared. This led us to hypothesize that the magnitude of genomic parallelism scales with genetic divergence between lineages, but whether this is the case and the underlying evolutionary processes remain unknown. Here, we resequenced seven parallel lineages of two Arabidopsis species, which repeatedly adapted to challenging alpine environments. By combining genome-wide divergence scans with model-based approaches, we detected a suite of 151 genes that show parallel signatures of positive selection associated with alpine colonization, involved in response to cold, high radiation, short season, herbivores, and pathogens. We complemented these parallel candidates with published gene lists from five additional alpine Brassicaceae and tested our hypothesis on a broad scale spanning ∼0.02 to 18 My of divergence. Indeed, we found quantitatively variable genomic parallelism whose extent significantly decreased with increasing divergence between the compared lineages. We further modeled parallel evolution over the Arabidopsis candidate genes and showed that a decreasing probability of repeated selection on the same standing or introgressed alleles drives the observed pattern of divergence-dependent parallelism. We therefore conclude that genetic divergence between populations, species, and genera, affecting the pool of shared variants, is an important factor in the predictability of genome evolution.

An overview of current population genomics methods for the analysis of whole-genome resequencing data in eukaryotes

Characterizing the population history of a species and identifying loci underlying local adaptation is crucial in functional ecology, evolutionary biology, conservation and agronomy. The constant improvement of high-throughput sequencing techniques has facilitated the production of whole genome data in a wide range of species. Population genomics now provides tools to better integrate selection into a historical framework, and take into account selection when reconstructing demographic history. However, this improvement has come with a profusion of analytical tools that can confuse and discourage users. Such confusion limits the amount of information effectively retrieved from complex genomic data sets, and impairs the diffusion of the most recent analytical tools into fields such as conservation biology. It may also lead to redundancy among methods. To address these isssues, we propose an overview of more than 100 state-of-the-art methods that can deal with whole genome data. We summarize the strategies they use to infer demographic history and selection, and discuss some of their limitations. A website listing these methods is available at www.methodspopgen.com.

Population genomic structure of Eurasian and African foxtail millet landrace accessions inferred from genotyping-by-sequencing

Foxtail millet [Setaria italica (L.) P. Beauv.] is the second most important millet species globally and is adapted to cultivation in diverse environments. Like its wild progenitor, green foxtail [S. viridis (L.) P. Beauv.], it is a model species for C4 photosynthetic pathways and stress tolerance genes in related bioenergy crops. We addressed questions regarding the evolution and spread of foxtail millet through a population genomic study of landraces from across its cultivated range in Europe, Asia, and Africa. We sought to determine population genomic structure and the relationship of domesticated lineages relative to green foxtail. Further, we aimed to identify genes involved in environmental stress tolerance that have undergone differential selection between geographical and genetic groups. Foxtail millet landrace accessions (n = 328) and green foxtail accessions (n = 12) were sequenced by genotyping-by-sequencing (GBS). After filtering, 5,677 single nucleotide polymorphisms (SNPs) were retained for the combined foxtail millet-green foxtail dataset and 5,020 for the foxtail millet dataset. We extended geographic coverage of green foxtail by including previously published GBS sequence tags, yielding a 4,515-SNP dataset for phylogenetic reconstruction. All foxtail millet samples were monophyletic relative to green foxtail, suggesting a single origin of foxtail millet, although no group of foxtail millet was clearly the most ancestral. Four genetic clusters were found within foxtail millet, each with a distinctive geographical distribution. These results, together with archaeobotanical evidence, suggest plausible routes of spread of foxtail millet. Selection scans identified nine candidate genes potentially involved in environmental adaptations, particularly to novel climates encountered, as domesticated foxtail millet spread to new altitudes and latitudes.

Altitudinal differentiation in the leaf wax-mediated flowering bud protection against frost in a perennial Arabidopsis

An altitudinal gradient of leaf water repellency is often observed between and within species. In a previous study of Arabidopsis halleri, cauline leaves (stem leaves that wrap flowering buds) showed higher water repellency in exposed semi-alpine plants than in understory low-elevation plants. Here, we examined altitudinal variations in the cuticular wax content of the leaf surface and experimentally evaluated the role of high water repellency of cauline leaves. Leaf cuticular wax was analysed using comprehensive two-dimensional gas chromatography (GC)-mass spectrometry and a GC-flame ionisation detector. Young flowering buds wrapped by cauline leaves were exposed to freezing temperatures with or without water, and frost damage to the flowering buds was compared between plants from semi-alpine and low-elevation habitats. Higher amounts of C29, C31, and C33 alkanes were observed in the cauline leaves of semi-alpine plants than in those of low-elevation plants. In the freezing experiment, water application increased damage to the flowering buds of low-elevation plants, and the extent of damage to the flowering buds was lower in semi-alpine plants than in low-elevation plants when water was applied to the plant surface. Genetic variations in the amounts of alkanes on the leaf surface depending on the altitude occurred specifically in cauline leaves. Our results indicate that the water repellency of cauline leaves presumably minimises frost damage to flowering buds at high altitudes.

Parallel Alpine Differentiation in Arabidopsis arenosa

Parallel evolution provides powerful natural experiments for studying repeatability of evolution and genomic basis of adaptation. Well-documented examples from plants are, however, still rare, as are inquiries of mechanisms driving convergence in some traits while divergence in others. Arabidopsis arenosa, a predominantly foothill species with scattered morphologically distinct alpine occurrences is a promising candidate. Yet, the hypothesis of parallelism remained untested. We sampled foothill and alpine populations in all regions known to harbor the alpine ecotype and used SNP genotyping to test for repeated alpine colonization. Then, we combined field surveys and a common garden experiment to quantify phenotypic parallelism. Genetic clustering by region but not elevation and coalescent simulations demonstrated parallel origin of alpine ecotype in four mountain regions. Alpine populations exhibited parallelism in height and floral traits which persisted after two generations in cultivation. In contrast, leaf traits were distinctive only in certain region(s), reflecting a mixture of plasticity and genetically determined non-parallelism. We demonstrate varying degrees and causes of parallelism and non-parallelism across populations and traits within a plant species. Parallel divergence along a sharp elevation gradient makes A. arenosa a promising candidate for studying genomic basis of adaptation.

Genomic signatures of convergent adaptation to Alpine environments in three Brassicaceae species

It has long been discussed to what extent related species develop similar genetic mechanisms to adapt to similar environments. Most studies documenting such convergence have either used different lineages within species or surveyed only a limited portion of the genome. Here, we investigated whether similar or different sets of orthologous genes were involved in genetic adaptation of natural populations of three related plant species to similar environmental gradients in the Alps. We used whole-genome pooled population sequencing to study genome-wide SNP variation in 18 natural populations of three Brassicaceae (Arabis alpina, Arabidopsis halleri, and Cardamine resedifolia) from the Swiss Alps. We first de novo assembled draft reference genomes for all three species. We then ran population and landscape genomic analyses with ~3 million SNPs per species to look for shared genomic signatures of selection and adaptation in response to similar environmental gradients acting on these species. Genes with a signature of convergent adaptation were found at significantly higher numbers than expected by chance. The most closely related species pair showed the highest relative over-representation of shared adaptation signatures. Moreover, the identified genes of convergent adaptation were enriched for nonsynonymous mutations, suggesting functional relevance of these genes, even though many of the identified candidate genes have hitherto unknown or poorly described functions based on comparison with Arabidopsis thaliana. We conclude that adaptation to heterogeneous Alpine environments in related species is partly driven by convergent evolution, but that most of the genomic signatures of adaptation remain species-specific.

Genetic consequences of being a dwarf: do evolutionary changes in life-history traits influence gene flow patterns in populations of the world's smallest goldenrod?

BACKGROUND AND AIMS

Contrasting life-history traits can evolve through generations of dwarf plant ecotypes, yet such phenotypic changes often involve decreased plant size and reproductive allocation, which can configure seed dispersal patterns and, subsequently, population demography. Therefore, evolutionary transitions to dwarfism can represent good study systems to test the roles of life-history traits in population demography by comparing genetic structure between related but phenotypically divergent ecotypes.

METHODS

In this study, we examined an ecotypic taxon pair of the world's smallest goldenrod (stem height 2.6 cm) in alpine habitats and its closely related lowland taxon (30-40 cm) found on Yakushima Island, Japan. Genetic variation in chloroplast DNA sequences, nuclear microsatellites and genome-wide single-nucleotide polymorphisms were used to investigate 197 samples from 16 populations, to infer the population genetic demography and compare local genetic structure of the ecotypes.

KEY RESULTS

We found a pronounced level of genetic differentiation among alpine dwarf populations, which were much less geographically isolated than their lowland counterparts. In particular, several neighbouring dwarf populations (located ~500 m apart) harboured completely different sets of chloroplast haplotypes and nuclear genetic clusters. Demographic modelling revealed that the dwarf populations have not exchanged genes at significant levels after population divergence.

CONCLUSIONS

These lines of evidence suggest that substantial effects of genetic drift have operated on these dwarf populations. The low-growing stature and reduced fecundity (only 3.1 heads per plant) of the dwarf plants may have reduced gene flow and rare long-distance seed dispersal among habitat patches, although the effects of life-history traits require further evaluation using ecological approaches.

Arabidopsis halleri: a perennial model system for studying population differentiation and local adaptation

Local adaptation is assumed to occur when populations differ in a phenotypic trait or a set of traits, and such variation has a genetic basis. Here, we introduce Arabidopsis halleri and its life history as a perennial model system to study population differentiation and local adaptation. Studies on altitudinal adaptation have been conducted in two regions: Mt. Ibuki in Japan and the European Alps. Several studies have demonstrated altitudinal adaptation in ultraviolet-B (UV-B) tolerance, leaf water repellency against spring frost and anti-herbivore defences. Studies on population differentiation in A. halleri have also focused on metal hyperaccumulation and tolerance to heavy metal contamination. In these study systems, genome scans to identify candidate genes under selection have been applied. Lastly, we briefly discuss how RNA-Seq can broaden phenotypic space and serve as a link to underlying mechanisms. In conclusion, A. halleri provides us with opportunities to study population differentiation and local adaptation, and relate these to the genetic systems underlying target functional traits.

Differential retention of transposable element-derived sequences in outcrossing Arabidopsis genomes

Background

Transposable elements (TEs) are genomic parasites with major impacts on host genome architecture and host adaptation. A proper evaluation of their evolutionary significance has been hampered by the paucity of short scale phylogenetic comparisons between closely related species. Here, we characterized the dynamics of TE accumulation at the micro-evolutionary scale by comparing two closely related plant species, Arabidopsis lyrata and A. halleri.

Results

Joint genome annotation in these two outcrossing species confirmed that both contain two distinct populations of TEs with either 'recent' or 'old' insertion histories. Identification of rare segregating insertions suggests that diverse TE families contribute to the ongoing dynamics of TE accumulation in the two species. Orthologous TE fragments (i.e. those that have been maintained in both species), tend to be located closer to genes than those that are retained in one species only. Compared to non-orthologous TE insertions, those that are orthologous tend to produce fewer short interfering RNAs, are less heavily methylated when found within or adjacent to genes and these tend to have lower expression levels. These findings suggest that long-term retention of TE insertions reflects their frequent acquisition of adaptive roles and/or the deleterious effects of removing nearly neutral TE insertions when they are close to genes.

Conclusion

Our results indicate a rapid evolutionary dynamics of the TE landscape in these two outcrossing species, with an important input of a diverse set of new insertions with variable propensity to resist deletion.

Genomic Patterns of Local Adaptation under Gene Flow in Arabidopsis lyrata

Short-scale local adaptation is a complex process involving selection, migration, and drift. The expected effects on the genome are well grounded in theory but examining these on an empirical level has proven difficult, as it requires information about local selection, demographic history, and recombination rate variation. Here, we use locally adapted and phenotypically differentiated Arabidopsis lyrata populations from two altitudinal gradients in Norway to test these expectations at the whole-genome level. Demography modeling indicates that populations within the gradients diverged <2 kya and that the sites are connected by gene flow. The gene flow estimates are, however, highly asymmetric with migration from high to low altitudes being several times more frequent than vice versa. To detect signatures of selection for local adaptation, we estimate patterns of lineage-specific differentiation among these populations. Theory predicts that gene flow leads to concentration of adaptive loci in areas of low recombination; a pattern we observe in both lowland-alpine comparisons. Although most selected loci display patterns of conditional neutrality, we found indications of genetic trade-offs, with one locus particularly showing high differentiation and signs of selection in both populations. Our results further suggest that resistance to solar radiation is an important adaptation to alpine environments, while vegetative growth and bacterial defense are indicated as selected traits in the lowland habitats. These results provide insights into genetic architectures and evolutionary processes driving local adaptation under gene flow. We also contribute to understanding of traits and biological processes underlying alpine adaptation in northern latitudes.

Plasticity of functional traits and optimality of biomass allocation in elevational ecotypes of Arabidopsis halleri grown at different soil nutrient availabilities

In mountainous areas, plant distribution is constrained by various environmental stresses. Plasticity and constancy in plant functional traits may relate to optimal strategies at respective habitats and to ecotypic differentiation along elevation. Although plant biomass allocation has been extensively studied in relation to adaptation to soil nutrient availability along elevation, its optimality is still poorly understood. We examined soil nutrient availability in the field and conducted growth analysis for two elevational ecotypes of Arabidopsis halleri grown under different nutrient availabilities. We determined plasticity in morphological and physiological traits and evaluated optimal biomass allocation using an optimality model. Our field investigation indicated that soil nitrogen (N) availability increased rather than decreased with increasing elevation. Our growth analysis revealed that lowland ecotype was more plastic in morphological variables and N concentrations, whereas the highland ecotype was more plastic in other physiological variables such as the net assimilation rate (NAR). The leaf mass ratio (LMR) in the lowland ecotype was moderately plastic at the whole range of N availabilities, whereas LMR in the highland ecotype was very plastic at higher N availabilities only. The optimality model indicated that the LMR of the lowland ecotype was nearly optimal throughout the range of studied N availabilities, whereas that of the highland ecotype was suboptimal at low N availability. These results suggest that highland ecotype is adapted only to high N availability, whereas the lowland ecotype is adapted to a relatively wide range of N availabilities as a result of natural selection in their respective habitats. We conclude that an adaptive differentiation has occurred between the two ecotypes and plasticity in the biomass allocation is directly related to its optimization in changing environments.

Trait differentiation and adaptation of plants along elevation gradients

Studies of genetic adaptation in plant populations along elevation gradients in mountains have a long history, but there has until now been neither a synthesis of how frequently plant populations exhibit adaptation to elevation nor an evaluation of how consistent underlying trait differences across species are. We reviewed studies of adaptation along elevation gradients (i) from a meta-analysis of phenotypic differentiation of three traits (height, biomass and phenology) from plants growing in 70 common garden experiments; (ii) by testing elevation adaptation using three fitness proxies (survival, reproductive output and biomass) from 14 reciprocal transplant experiments; (iii) by qualitatively assessing information at the molecular level, from 10 genomewide surveys and candidate gene approaches. We found that plants originating from high elevations were generally shorter and produced less biomass, but phenology did not vary consistently. We found significant evidence for elevation adaptation in terms of survival and biomass, but not for reproductive output. Variation in phenotypic and fitness responses to elevation across species was not related to life history traits or to environmental conditions. Molecular studies, which have focussed mainly on loci related to plant physiology and phenology, also provide evidence for adaptation along elevation gradients. Together, these studies indicate that genetically based trait differentiation and adaptation to elevation are widespread in plants. We conclude that a better understanding of the mechanisms underlying adaptation, not only to elevation but also to environmental change, will require more studies combining the ecological and molecular approaches.

Phylogeographic analysis of the East Asian goldenrod (Solidago virgaurea complex, Asteraceae) reveals hidden ecological diversification with recurrent formation of ecotypes

BACKGROUND AND AIMS

The processes and mechanisms underlying lineage diversification are major topics in evolutionary biology. Eurasian goldenrod species of the Solidago virgaurea complex show remarkable morphological and ecological diversity in the Japanese Archipelago, with ecotypic taxa well adapted to specific environments (climate, edaphic conditions and disturbance regimes). The species complex is a suitable model to investigate the evolutionary processes of actively speciating plant groups, due to its ability to evolve in relation to environmental adaptation and its historical population dynamics.

METHODS

Two chloroplast markers, 18 nuclear microsatellite markers and ddRAD-sequencing were used to infer population genetic demography of S. virgaurea complex with its related species/genera.

KEY RESULTS

Our analysis showed that populations in Japan form an evolutionary unit, which was genetically diverged from adjacent continental populations. The phylogenetic structure within the archipelago strongly corresponds to the geography, but interestingly there is no concordance between genetic structure and ecotypic boundaries; neighbouring populations of distinct ecotypes share a genetic background.

CONCLUSIONS

We propose that the traits specific to the ecotypic entities are maintained by natural selection or are very recently generated and have little effect on the genomes, making genome-wide genetic markers unsuitable for detecting ecotypic differentiation. Furthermore, some sporadically distributed taxa (found as rheophytes and alpine plants) were repeatedly generated from a more widespread taxon in geographically distant areas by means of selection. Overall, this study showed that the goldenrod complex has a high ability to evolve, enabling rapid ecological diversification over a recent timeframe.

A Genomic Map of Climate Adaptation in Arabidopsis thaliana at a Micro-Geographic Scale

Understanding the genetic bases underlying climate adaptation is a key element to predict the potential of species to face climate warming. Although substantial climate variation is observed at a micro-geographic scale, most genomic maps of climate adaptation have been established at broader geographical scales. Here, by using a Pool-Seq approach combined with a Bayesian hierarchical model that control for confounding by population structure, we performed a genome-environment association (GEA) analysis to investigate the genetic basis of adaptation to six climate variables in 168 natural populations of Arabidopsis thaliana distributed in south-west of France. Climate variation among the 168 populations represented up to 24% of climate variation among 521 European locations where A. thaliana inhabits. We identified neat and strong peaks of association, with most of the associated SNPs being significantly enriched in likely functional variants and/or in the extreme tail of genetic differentiation among populations. Furthermore, genes involved in transcriptional mechanisms appear predominant in plant functions associated with local climate adaptation. Globally, our results suggest that climate adaptation is an important driver of genomic variation in A. thaliana at a small spatial scale and mainly involves genome-wide changes in fundamental mechanisms of gene regulation. The identification of climate-adaptive genetic loci at a micro-geographic scale also highlights the importance to include within-species genetic diversity in ecological niche models for projecting potential species distributional shifts over short geographic distances.

Simultaneous evaluation of the effects of geographic, environmental and temporal isolation in ecotypic populations of Solidago virgaurea

Early stages of ecological speciation can create populations with an ecology and reproduction timing distinct from those of related populations. Landscape genetic models incorporating environmental heterogeneity and population-specific reproductive traits enable the processes of population genetic differentiation to be inferred. We investigated genome-wide genetic variation in ecotypic populations of Solidago virgaurea sensu lato, a herbaceous plant inhabiting a wide range of habitats (woodlands, serpentine barrens and alpine grasslands) and displaying remarkable variation in flowering time. Simultaneous evaluation of environmental factors revealed an overwhelming effect of soil type differences on neutral genetic differentiation, compared with elevational differences. This result probably reflects the abrupt environmental changes generated by geological boundaries, whereas mountain slopes exhibit clinal changes, facilitating gene exchange between neighbouring populations. Temporal isolation was positively associated with genetic differentiation, with some early-flowering serpentine populations having allele frequencies distinct from adjacent nonserpentine populations. Overall, this study highlights the importance of ecological processes and of evolution of flowering time to promote genetic differentiation of S. virgaurea populations in a complex landscape.

RADseq provides evidence for parallel ecotypic divergence in the autotetraploid Cochlearia officinalis in Northern Norway

Speciation encompasses a continuum over time from freely interbreeding populations to reproductively isolated species. Along this process, ecotypes - the result of local adaptation - may be on the road to new species. We investigated whether three autotetraploid Cochlearia officinalis ecotypes, adapted to different habitats (beach, estuary, spring), are genetically differentiated and result from parallel ecotypic divergence in two distinct geographical regions. We obtained genetic data from thousands of single nucleotide polymorphisms (SNPs) from restriction-site associated DNA sequencing (RADseq) and from six microsatellite markers for 12 populations to assess genetic divergence at ecotypic, geographic and population level. The genetic patterns support differentiation among ecotypes as suggested by morphology and ecology. The data fit a scenario where the ancestral beach ecotype has recurrently and polytopically given rise to the estuary and spring ecotypes. Several ecologically-relevant loci with consistent non-random segregating patterns are identified across the recurrent origins, in particular around genes related to salt stress. Despite being ecologically distinct, the Cochlearia ecotypes still represent an early stage in the process of speciation, as reproductive isolation has not (yet) developed. A sequenced annotated genome is needed to specifically target candidate genes underlying local adaptation.

Fine-scale frequency differentiation along a herbivory gradient in the trichome dimorphism of a wild Arabidopsis

Geographic variation is commonly observed in plant resistance traits, where plant species might experience different selection pressure across a heterogeneous landscape. Arabidopsis halleri subsp. gemmifera is dimorphic for trichome production, generating two morphs, trichome-producing (hairy) and trichomeless (glabrous) plants. Trichomes of A. halleri are known to confer resistance against the white butterfly, cabbage sawfly, and brassica leaf beetle, but not against flea beetles. We combined leaf damage, microclimate, and microsatellite loci data of 26 A. halleri populations in central Japan, to explore factors responsible for fine-scale geographic variation in the morph frequency. We found that hairy plants were less damaged than glabrous plants within populations, but the among-site variation was the most significant source of variation in the individual-level damage. Fixation index (Gst″) of a putative trichome locus exhibited a significant divergence along population-level damage with an exception of an outlier population, inferring the local adaptation to herbivory. Notably, this outlier was a population wherein our previous study reported a balancing role of the brassica leaf beetle Phaedon brassicae on the morph frequency. This differentiation of the trichome locus was unrelated to neutral genetic differentiation (evaluated by Gst″ of microsatellite loci) and meteorological factors (including temperature and solar radiation). The present findings, combined with those of our previous work, provide suggestive evidence that herbivore-driven divergence and occasional outbreak of a specific herbivore have jointly contributed to the ecogeographic pattern in the frequency of two morphs.

The genomic basis of adaptation in plants

Plants are powerful models for the study of adaptive evolution. Since they are rooted in place, they must directly face environmental insults, making adaptation to local conditions vital. In addition to adaptation to natural conditions, some plant species have held a central role in human subsistence over the past several thousand years. In these species, humans exerted strong selective pressures on traits of agricultural importance. Recently, an increasing number of studies have aimed to identify the genomic basis of adaptation. These studies have provided insights into the mechanisms through which the raw materials of adaptation were introduced as well as the modes of adaptation in wild and domesticated species.

Genomic studies of adaptive evolution in outcrossing Arabidopsis species

Large-scale population genomic approaches have very recently been fruitfully applied to the Arabidopsis relatives Arabidopsis halleri, A. lyrata and especially A. arenosa. In contrast to A. thaliana, these species are obligately outcrossing and thus the footprints of natural selection are more straightforward to detect. Furthermore, both theoretical and empirical studies indicate that outcrossers are better able to evolve in response to selection pressure. As a result, recent work in these species serves as a paradigm of population genomic studies of adaptation both to environmental as well as intracellular challenges.

Genomic Scans across Three Eucalypts Suggest that Adaptation to Aridity is a Genome-Wide Phenomenon

Widespread species spanning strong environmental (e.g., climatic) gradients frequently display morphological and physiological adaptations to local conditions. Some adaptations are common to different species that occupy similar environments. However, the genomic architecture underlying such convergent traits may not be the same between species. Using genomic data from previous studies of three widespread eucalypt species that grow along rainfall gradients in southern Australia, our probabilistic approach provides evidence that adaptation to aridity is a genome-wide phenomenon, likely to involve multiple and diverse genes, gene families and regulatory regions that affect a multitude of complex genetic and biochemical processes.

Field Guide to Plant Model Systems

For the past several decades, advances in plant development, physiology, cell biology, and genetics have relied heavily on the model (or reference) plant Arabidopsis thaliana. Arabidopsis resembles other plants, including crop plants, in many but by no means all respects. Study of Arabidopsis alone provides little information on the evolutionary history of plants, evolutionary differences between species, plants that survive in different environments, or plants that access nutrients and photosynthesize differently. Empowered by the availability of large-scale sequencing and new technologies for investigating gene function, many new plant models are being proposed and studied.

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site-specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole-genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree-based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.

The genetics of extreme microgeographic adaptation: an integrated approach identifies a major gene underlying leaf trichome divergence in Yellowstone Mimulus guttatus

Microgeographic adaptation provides a particularly interesting context for understanding the genetic basis of phenotypic divergence and may also present unique empirical challenges. In particular, plant adaptation to extreme soil mosaics may generate barriers to gene flow or shifts in mating system that confound simple genomic scans for adaptive loci. Here, we combine three approaches - quantitative trait locus (QTL) mapping of candidate intervals in controlled crosses, population resequencing (PoolSeq) and analyses of wild recombinant individuals - to investigate one trait associated with Mimulus guttatus (yellow monkeyflower) adaptation to geothermal soils in Yellowstone National Park. We mapped a major QTL causing dense leaf trichomes in thermally adapted plants to a <50-kb region of linkage Group 14 (Tr14) previously implicated in trichome divergence between independent M. guttatus populations. A PoolSeq scan of Tr14 region revealed a cluster of six genes, coincident with the inferred QTL peak, with high allele frequency differences sufficient to explain observed phenotypic differentiation. One of these, the R2R3 MYB transcription factor Migut.N02661, is a plausible functional candidate and was also strongly associated (r²  = 0.27) with trichome phenotype in analyses of wild-collected admixed individuals. Although functional analyses will be necessary to definitively link molecular variants in Tr14 with trichome divergence, our analyses are a major step in that direction. They point to a simple, and parallel, genetic basis for one axis of Mimulus guttatus adaptation to an extreme habitat, suggest a broadly conserved genetic basis for trichome variation across flowering plants and pave the way for further investigations of this challenging case of microgeographic incipient speciation.

Genomic and phenotypic differentiation of Arabidopsis thaliana along altitudinal gradients in the North Italian Alps

Altitudinal gradients in mountain regions are short-range clines of different environmental parameters such as temperature or radiation. We investigated genomic and phenotypic signatures of adaptation to such gradients in five Arabidopsis thaliana populations from the North Italian Alps that originated from 580 to 2350 m altitude by resequencing pools of 19-29 individuals from each population. The sample includes two pairs of low- and high-altitude populations from two different valleys. High-altitude populations showed a lower nucleotide diversity and negative Tajima's D values and were more closely related to each other than to low-altitude populations from the same valley. Despite their close geographic proximity, demographic analysis revealed that low- and high-altitude populations split between 260 000 and 15 000 years before present. Single nucleotide polymorphisms whose allele frequencies were highly differentiated between low- and high-altitude populations identified genomic regions of up to 50 kb length where patterns of genetic diversity are consistent with signatures of local selective sweeps. These regions harbour multiple genes involved in stress response. Variation among populations in two putative adaptive phenotypic traits, frost tolerance and response to light/UV stress was not correlated with altitude. Taken together, the spatial distribution of genetic diversity reflects a potentially adaptive differentiation between low- and high-altitude populations, whereas the phenotypic differentiation in the two traits investigated does not. It may resemble an interaction between adaptation to the local microhabitat and demographic history influenced by historical glaciation cycles, recent seed dispersal and genetic drift in local populations.

Genome sequencing of Metrosideros polymorpha (Myrtaceae), a dominant species in various habitats in the Hawaiian Islands with remarkable phenotypic variations

Whole genome sequences, which can be provided even for non-model organisms owing to high-throughput sequencers, are valuable in enhancing the understanding of adaptive evolution. Metrosideros polymorpha, a tree species endemic to the Hawaiian Islands, occupies a wide range of ecological habitats and shows remarkable polymorphism in phenotypes among/within populations. The biological functions of genetic variations observed within this species could provide significant insights into the adaptive radiation found in a single species. Here de novo assembled genome sequences of M. polymorpha are presented to reveal basic genomic parameters about this species and to develop our knowledge of ecological divergences. The assembly yielded 304-Mbp genome sequences, half of which were covered by 19 scaffolds with >5 Mbp, and contained 30 K protein-coding genes. Demographic history inferred from the genome-wide heterozygosity indicated that this species experienced a dramatic rise and fall in the effective population size, possibly owing to past geographic or climatic changes in the Hawaiian Islands. This M. polymorpha genome assembly represents a high-quality genome resource useful for future functional analyses of both intra- and interspecies genetic variations or comparative genomics.

High light intensity plays a major role in emergence of population level variation in Arabidopsis thaliana along an altitudinal gradient

Environmental conditions play an important role in the emergence of genetic variations in natural populations. We identified genome-wide patterns of nucleotide variations in the coding regions of natural Arabidopsis thaliana populations. These populations originated from 700 m to 3400 m a.m.s.l. in the Western Himalaya. Using a pooled RNA-Seq approach, we identified the local and global level population-specific SNPs. The biological functions of the SNP-containing genes were primarily related to the high light intensity prevalent at high-altitude regions. The novel SNPs identified in these genes might have arisen de novo in these populations. In another approach, the F_STs of SNP-containing genes were correlated with the corresponding climatic factors. ‘Radiation in the growing season’ was the only environmental factor found to be strongly correlated with the gene-level F_STs. In both the approaches, the high light intensity was identified as the primary abiotic stress associated with the variations in these populations. The differential gene expression analysis between field and controlled condition grown plants also showed high light intensity as the primary abiotic stress, particularly for the high altitude populations. Our results provide a genome-wide perspective of nucleotide variations in populations along altitudinal gradient and their putative role in emergence of these variations.