Parallelism in gene expression between foothill and alpine ecotypes in Arabidopsis arenosa

Repeated colonisation of alpine habitats by Arabidopsis arenosa involved parallel adjustments of leaf cuticle traits

Cuticle function can be pivotal to plant success in different environments. Yet, the occurrence of intraspecific adjustments in cuticle traits resulting from acclimation or adaptation to different habitats remains poorly understood. Here, we used genetically well-characterised populations of Arabidopsis arenosa to investigate whether cuticle traits were adjusted as part of the parallel evolution from a foothill to an alpine ecotype. Six alpine and six foothill populations, representing at least three independent evolutionary origins of an alpine ecotype, were used in reciprocal transplantation experiments, to investigate cuticle traits at the eco-physiological, biochemical and structural levels. The genetic basis behind these traits was assessed by combining selection scans and differential gene expression analysis. Overall, alpine populations showed reduced cuticular transpiration in conjunction with consistently altered cuticular wax composition, with higher accumulation of two fatty alcohols and two iso-alkanes. Genomic analysis unravelled nine genes associated with cuticular wax metabolism showing allelic differentiation in alpine compared to lowland populations. In silico gene expression analysis revealed differences between ecotypes for several genes related to cuticle metabolism. Repeated ecotypic differentiation in cuticle traits together with the genetic architecture of the alpine ecotype points at an adaptive value of cuticle adjustments for the colonisation of alpine habitats.

Role of phenotypic and transcriptomic plasticity in alpine adaptation of Arabidopsis arenosa

Plasticity is an important component of the response of organism to environmental changes, but whether plasticity facilitates adaptation is still largely debated. Using transcriptomic and phenotypic data, we explored the evolution of ancestral plasticity during alpine colonization in Arabidopsis arenosa. We leveraged naturally replicated adaptation in four distinct mountain regions in Central Europe. We sampled seeds from ancestral foothill and independently formed alpine populations in each region and raised them in growth chambers under conditions approximating their natural environments. We gathered RNA-seq and genetic data of 48 and 63 plants and scored vegetative and flowering traits in 203 and 272 plants respectively. Then, we compared gene expression and trait values over two treatments differing in temperature and irradiance and elevations of origin and quantified the extent of ancestral and derived plasticity. At the transcriptomic level, initial plastic changes tended to be more reinforced than reversed in adapted alpine populations. Genes showing reinforcement were involved in the stress response, developmental processes and morphogenesis and those undergoing reversion were related to the stress response (light and biotic stress). At the phenotypic level, initial plastic changes in all but one trait were also reinforced supporting a facilitating role of phenotypic plasticity during colonization of an alpine environment. Our results contrasted with previous studies that showed generally higher reversion than reinforcement and supported the idea that ancestral plasticity tends to be reinforced in the context of alpine adaptation. However, plasticity may also be the source of potential maladaptation, especially at the transcriptomic level.

Understanding natural selection and similarity: Convergent, parallel and repeated evolution

Parallel and convergent evolution offer some of the most compelling evidence for the significance of natural selection in evolution, as the emergence of similar adaptive solutions is unlikely to occur by random chance alone. However, these terms are often employed inconsistently, leading to misinterpretation and confusion, and recently proposed definitions have unintentionally diminished the emphasis on the evolution of similar adaptive solutions. Here, I examine various conceptual frameworks and definitions related to parallel and convergent evolution and propose a consolidated framework that enhances our comprehension of these evolutionary patterns. The primary aim of this framework is to harmonize the concepts of parallel and convergent evolution together with natural selection and the idea of similarity. Both concepts involve the evolution of similar adaptive solutions as a result of environmental challenges. The distinction lies in ancestral phenotypes. Parallel evolution takes place when the ancestral phenotypes (before selection) of the lineages are similar. Convergent evolution happens when the lineages have distinct ancestral phenotypes (before selection). Because an ancestral-based distinction will inevitably lead to cases where uncertainty in the distinction may arise, the framework includes a general term, repeated evolution, which can be used as a term applying to the evolution of similar phenotypes and genotypes as well as similar responses to environmental pressures. Based on the argument that genetic similarity may frequently arise without selection, the framework posits that the similarity of genetic sequences is not of great interest unless linked to the actions of natural selection or to the origins (mutation, standing genetic variation, gene flow) and locations of the similar sequences.

Genetic assimilation of ancestral plasticity during parallel adaptation to zinc contamination in Silene uniflora

Phenotypic plasticity in ancestral populations is hypothesized to facilitate adaptation, but evidence is piecemeal and often contradictory. Further, whether ancestral plasticity increases the probability of parallel adaptive changes has not been explored. The most general finding is that ancestral responses to a new environment are reversed following adaptation (known as reversion). We investigated the contribution of ancestral plasticity to adaptive evolution of gene expression in two independently evolved lineages of zinc-tolerant Silene uniflora. We found that the general pattern of reversion is driven by the absence of a widespread stress response in zinc-adapted plants compared with zinc-sensitive plants. We show that ancestral plasticity that moves expression closer to the optimum value in the new environment influences the evolution of gene expression among genes that are likely to be involved in adaptation and increases the chance that genes are recruited repeatedly during adaptation. However, despite convergence in gene expression levels between independently adapted lineages, ancestral plasticity does not influence how similar expression values of adaptive genes become. Surprisingly, we also observed that ancestral plasticity that increases fitness often becomes genetically determined and fixed, that is, genetically assimilated. These results emphasize the important role of ancestral plasticity in parallel adaptation.

Environmental response in gene expression and DNA methylation reveals factors influencing the adaptive potential of Arabidopsis lyrata

Understanding what factors influence plastic and genetic variation is valuable for predicting how organisms respond to changes in the selective environment. Here, using gene expression and DNA methylation as molecular phenotypes, we study environmentally induced variation among Arabidopsis lyrata plants grown at lowland and alpine field sites. Our results show that gene expression is highly plastic, as many more genes are differentially expressed between the field sites than between populations. These environmentally responsive genes evolve under strong selective constraint - the strength of purifying selection on the coding sequence is high, while the rate of adaptive evolution is low. We find, however, that positive selection on cis-regulatory variants has likely contributed to the maintenance of genetically variable environmental responses, but such variants segregate only between distantly related populations. In contrast to gene expression, DNA methylation at genic regions is largely insensitive to the environment, and plastic methylation changes are not associated with differential gene expression. Besides genes, we detect environmental effects at transposable elements (TEs): TEs at the high-altitude field site have higher expression and methylation levels, suggestive of a broad-scale TE activation. Compared to the lowland population, plants native to the alpine environment harbor an excess of recent TE insertions, and we observe that specific TE families are enriched within environmentally responsive genes. Our findings provide insight into selective forces shaping plastic and genetic variation. We also highlight how plastic responses at TEs can rapidly create novel heritable variation in stressful conditions.

Parallel Differentiation and Plastic Adjustment of Leaf Anatomy in Alpine Arabidopsis arenosa Ecotypes

Functional and structural adjustments of plants in response to environmental factors, including those occurring in alpine habitats, can result in transient acclimation, plastic phenotypic adjustments and/or heritable adaptation. To unravel repeatedly selected traits with potential adaptive advantage, we studied parallel (ecotypic) and non-parallel (regional) differentiation in leaf traits in alpine and foothill ecotypes of Arabidopsis arenosa. Leaves of plants from eight alpine and eight foothill populations, representing three independent alpine colonization events in different mountain ranges, were investigated by microscopy techniques after reciprocal transplantation. Most traits clearly differed between the foothill and the alpine ecotype, with plastic adjustments to the local environment. In alpine populations, leaves were thicker, with altered proportions of palisade and spongy parenchyma, and had fewer trichomes, and chloroplasts contained large starch grains with less stacked grana thylakoids compared to foothill populations. Geographical origin had no impact on most traits except for trichome and stomatal density on abaxial leaf surfaces. The strong parallel, heritable ecotypic differentiation in various leaf traits and the absence of regional effects suggests that most of the observed leaf traits are adaptive. These trait shifts may reflect general trends in the adaptation of leaf anatomy associated with the colonization of alpine habitats.

Repeated colonization of alpine habitats by Arabidopsis arenosa viewed through freezing resistance and ice management strategies

Success or failure of plants to cope with freezing temperatures can critically influence plant distribution and adaptation to new habitats. Especially in alpine environments, frost is a likely major selective force driving adaptation. In Arabidopsis arenosa (L.) Lawalrée, alpine populations have evolved independently in different mountain ranges, enabling studying mechanisms of acclimation and adaptation to alpine environments. We tested for heritable, parallel differentiation in freezing resistance, cold acclimation potential and ice management strategies using eight alpine and eight foothill populations. Plants from three European mountain ranges (Niedere Tauern, Făgăraș and Tatra Mountains) were grown from seeds of tetraploid populations in four common gardens, together with diploid populations from the Tatra Mountains. Freezing resistance was assessed using controlled freezing treatments and measuring effective quantum yield of photosystem II, and ice management strategies by infrared video thermography and cryomicroscopy. The alpine ecotype had a higher cold acclimation potential than the foothill ecotype, whereby this differentiation was more pronounced in tetraploid than diploid populations. However, no ecotypic differentiation was found in one region (Făgăraș), where the ancient lineage had a different evolutionary history. Upon freezing, an ice lens within a lacuna between the palisade and spongy parenchyma tissues was formed by separation of leaf tissues, a mechanism not previously reported for herbaceous species. The dynamic adjustment of freezing resistance to temperature conditions may be particularly important in alpine environments characterized by large temperature fluctuations. Furthermore, the formation of an extracellular ice lens may be a useful strategy to avoid tissue damage during freezing.

Ploidy and local environment drive intraspecific variation in endoreduplication in Arabidopsis arenosa

PREMISE

Endoreduplication, nonheritable duplication of a nuclear genome, is widespread in plants and plays a role in developmental processes related to cell differentiation. However, neither ecological nor cytological factors influencing intraspecific variation in endoreduplication are fully understood.

METHODS

We cultivated plants covering the range-wide natural diversity of diploid and tetraploid populations of Arabidopsis arenosa in common conditions to investigate the effect of original ploidy level on endoreduplication. We also raised plants from several foothill and alpine populations from different lineages and of both ploidies to test for the effect of elevation. We determined the endoreduplication level in leaves of young plants by flow cytometry. Using RNA-seq data available for our populations, we analyzed gene expression analysis in individuals that differed in endoreduplication level.

RESULTS

We found intraspecific variation in endoreduplication that was mainly driven by the original ploidy level of populations, with significantly higher endoreduplication in diploids. An effect of elevation was also found within each ploidy, yet its direction exhibited rather regional-specific patterns. Transcriptomic analysis comparing individuals with high vs. low endopolyploidy revealed a majority of differentially expressed genes related to the stress and hormone response and to modifications especially in the cell wall and in chloroplasts.

CONCLUSIONS

Our results support the general assumption of higher potential of low-ploidy organisms to undergo endoreduplication and suggest that endoreduplication is further integrated within the stress response pathways for a fine-tune adjustment of the endoreduplication process to their local environment.

The role of alternative splicing in adaptation and evolution

Regulation of gene expression plays a central role in adaptive divergence and evolution. Although the role of gene regulation in microevolutionary processes is gaining wide acceptance, most studies have only investigated the evolution of transcript levels, ignoring the potentially significant role of transcript structures. We argue that variation in alternative splicing plays an important and widely unexplored role in adaptation (e.g., by increasing transcriptome and/or proteome diversity, or buffering potentially deleterious genetic variation). New studies increasingly highlight the potential for independent evolution in alternative splicing and transcript level, providing alternative paths for selection to act upon. We propose that alternative splicing and transcript levels can provide contrasting, nonredundant mechanisms of equal importance for adaptive diversification of gene function and regulation.

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.

Transcriptional activity of transposable elements along an elevational gradient in Arabidopsis arenosa

Background

Plant genomes can respond rapidly to environmental changes and transposable elements (TEs) arise as important drivers contributing to genome dynamics. Although some elements were reported to be induced by various abiotic or biotic factors, there is a lack of general understanding on how environment influences the activity and diversity of TEs. Here, we combined common garden experiment with short-read sequencing to investigate genomic abundance and expression of 2245 consensus TE sequences (containing retrotransposons and DNA transposons) in an alpine environment in Arabidopsis arenosa. To disentangle general trends from local differentiation, we leveraged four foothill-alpine population pairs from different mountain regions. Seeds of each of the eight populations were raised under four treatments that differed in temperature and irradiance, two factors varying with elevation. RNA-seq analysis was performed on leaves of young plants to test for the effect of elevation and subsequently of temperature and irradiance on expression of TE sequences.

Results

Genomic abundance of the 2245 consensus TE sequences varied greatly between the mountain regions in line with neutral divergence among the regions, representing distinct genetic lineages of A. arenosa. Accounting for intraspecific variation in abundance, we found consistent transcriptomic response for some TE sequences across the different pairs of foothill-alpine populations suggesting parallelism in TE expression. In particular expression of retrotransposon LTR Copia (e.g. Ivana and Ale clades) and LTR Gypsy (e.g. Athila and CRM clades) but also non-LTR LINE or DNA transposon TIR MuDR consistently varied with elevation of origin. TE sequences responding specifically to temperature and irradiance belonged to the same classes as well as additional TE clades containing potentially stress-responsive elements (e.g. LTR Copia Sire and Tar, LTR Gypsy Reina).

Conclusions

Our study demonstrated that the A. arenosa genome harbours a considerable diversity of TE sequences whose abundance and expression response varies across its native range. Some TE clades may contain transcriptionally active elements responding to a natural environmental gradient. This may further contribute to genetic variation between populations and may ultimately provide new regulatory mechanisms to face environmental challenges.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-021-00236-0.