Deleterious phenotypes in wild Arabidopsis arenosa populations are common and linked to runs of homozygosity.

In this study, we aimed to systematically assess the frequency at which potentially deleterious phenotypes appear in natural populations of the outcrossing model plant Arabidopsis arenosa, and to establish their underlying genetics. For this purpose, we collected seeds from wild A. arenosa populations and screened over 2,500 plants for unusual phenotypes in the greenhouse. We repeatedly found plants with obvious phenotypic defects, such as small stature and necrotic or chlorotic leaves, among first-generation progeny of wild A. arenosa plants. Such abnormal plants were present in about 10% of maternal sibships, with multiple plants with similar phenotypes in each of these sibships, pointing to a genetic basis of the observed defects. A combination of transcriptome profiling, linkage mapping and genome-wide runs of homozygosity patterns using a newly assembled reference genome indicated a range of underlying genetic architectures associated with phenotypic abnormalities. This included evidence for homozygosity of certain genomic regions, consistent with alleles that are identical by descent being responsible for these defects. Our observations suggest that deleterious alleles with different genetic architectures are segregating at appreciable frequencies in wild A. arenosa populations.

Evolution of the Sabulina verna group (Caryophyllaceae) in Europe: A deep split, followed by secondary contacts, multiple allopolyploidization and colonization of challenging substrates.

One of the major goals of contemporary evolutionary biology is to elucidate the relative roles of allopatric and ecological differentiation and polyploidy in speciation. In this study, we address the taxonomically intricate Sabulina verna group, which has a disjunct Arctic-alpine postglacial range in Europe and occupies a broad range of ecological niches, including substrates toxic to plants. Using genome-wide ddRAD sequencing combined with morphometric analyses based on extensive sampling of 111 natural populations, we aimed to disentangle internal evolutionary relationships and examine their correspondence with the pronounced edaphic and ploidy diversity within the group. We identified two spatially distinct groups of diploids: a widespread Arctic-alpine group and a spatially restricted yet diverse Balkan group. Most tetraploids exhibited a considerably admixed ancestry derived from both these groups, suggesting their allopolyploid origin. Four genetic clusters in congruence with geography and mostly supported by morphological traits were recognized in the diploid Arctic-alpine group. Tetraploids are split into two distinct and geographically vicariant groups, indicating their repeated polytopic origin. Furthermore, our results also revealed at least five-fold parallel colonization of toxic substrates (serpentine and metalliferous), altogether demonstrating a complex interaction between geography, challenging substrates and polyploidy in the evolution of the group. Finally, we propose a new taxonomic treatment of this complex.

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.

The Use of Flow Cytometry for Estimating Genome Sizes and DNA Ploidy Levels in Plants.

Flow cytometry has emerged as a uniquely flexible, accurate, and widely applicable technology for the analysis of plant cells. One of its most important applications centers on the measurement of nuclear DNA contents. This chapter describes the essential features of this measurement, outlining the overall methods and strategies, but going on to provide a wealth of technical details to ensure the most accurate and reproducible results. The chapter is aimed to be equally accessible to experienced plant cytometrists as well as those newly entering the field. Besides providing a step-by-step guide for estimating genome sizes and DNA-ploidy levels from fresh tissues, special attention is paid to the use of seeds and desiccated tissues for such purposes. Methodological aspects regarding field sampling, transport, and storage of plant material are also given in detail. Finally, troubleshooting information for the most common problems that may arise during the application of these methods is provided.

Repeatedly Northwards and Upwards: Southern African Grasslands Fuel the Colonization of the African Sky Islands in Helichrysum (Compositae).

The Afromontane and Afroalpine areas constitute some of the main biodiversity hotspots of Africa. They are particularly rich in plant endemics, but the biogeographic origins and evolutionary processes leading to this outstanding diversity are poorly understood. We performed phylogenomic and biogeographic analyses of one of the most species-rich plant genera in these mountains, Helichrysum (Compositae-Gnaphalieae). Most previous studies have focused on Afroalpine elements of Eurasian origin, and the southern African origin of Helichrysum provides an interesting counterexample. We obtained a comprehensive nuclear dataset from 304 species (≈50% of the genus) using target-enrichment with the Compositae1061 probe set. Summary-coalescent and concatenation approaches combined with paralog recovery yielded congruent, well-resolved phylogenies. Ancestral range estimations revealed that Helichrysum originated in arid southern Africa, whereas the southern African grasslands were the source of most lineages that dispersed within and outside Africa. Colonization of the tropical Afromontane and Afroalpine areas occurred repeatedly throughout the Miocene-Pliocene. This timing coincides with mountain uplift and the onset of glacial cycles, which together may have facilitated both speciation and intermountain gene flow, contributing to the evolution of the Afroalpine flora.

Population Genomic Analysis of Diploid-Autopolyploid Species.

This chapter outlines an empirical analysis of genome-wide single-nucleotide polymorphism (SNP) variation and its underlying drivers among multiple natural populations within a diploid-autopolyploid species. The aim is to reconstruct the genetic structure among natural populations of varying ploidy and infer footprints of selection in these populations, framed around specific questions that are typically encountered when analyzing a mixed-ploidy data set,e.g., addressing the relevance of natural whole-genome duplication for speciation and adaptation. We briefly review the options for the analysis of polyploid population genomic data involving variant calling, population structure, demographic history inference, and selection scanning approaches. Further, we provide suggestions for methods and associated software, possible caveats, and examples of their application to mixed-ploidy and autopolyploid data sets.

Inbreeding depression in polyploid species: a meta-analysis.

Whole-genome duplication is a common mutation in eukaryotes with far-reaching phenotypic effects, the resulting morphological and fitness consequences and how they affect the survival of polyploid lineages are intensively studied. Another important factor may also determine the probability of establishment and success of polyploid lineages: inbreeding depression. Inbreeding depression is expected to play an important role in the establishment of neopolyploid lineages, their capacity to colonize new environments, and in the simultaneous evolution of ploidy and other life-history traits such as self-fertilization. Both theoretically and empirically, there is no consensus on the consequences of polyploidy on inbreeding depression. In this meta-analysis, we investigated the effect of polyploidy on the evolution of inbreeding depression, by performing a meta-analysis within angiosperm species. The main results of our study are that the consequences of polyploidy on inbreeding depression are complex and depend on the time since polyploidization. We found that young polyploid lineages have a much lower amount of inbreeding depression than their diploid relatives and their established counterparts. Natural polyploid lineages are intermediate and have a higher amount of inbreeding depression than synthetic neopolyploids, and a smaller amount than diploids, suggesting that the negative effect of polyploidy on inbreeding depression decreases with time since polyploidization.

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.

Genomic basis and phenotypic manifestation of (non-)parallel serpentine adaptation in Arabidopsis arenosa.

Parallel evolution is common in nature and provides one of the most compelling examples of rapid environmental adaptation. In contrast to the recent burst of studies addressing genomic basis of parallel evolution, integrative studies linking genomic and phenotypic parallelism are scarce. Edaphic islands of toxic serpentine soils provide ideal systems for studying rapid parallel adaptation in plants, imposing strong, spatially replicated selection on recently diverged populations. We leveraged threefold independent serpentine adaptation of Arabidopsis arenosa and combined reciprocal transplants, ion uptake phenotyping, and available genome-wide polymorphisms to test if parallelism is manifested to a similar extent at both genomic and phenotypic levels. We found pervasive phenotypic parallelism in functional traits yet with varying magnitude of fitness differences that was congruent with neutral genetic differentiation between populations. Limited costs of serpentine adaptation suggest absence of soil-driven trade-offs. On the other hand, the genomic parallelism at the gene level was significant, although relatively minor. Therefore, the similarly modified phenotypes, for example, of ion uptake arose possibly by selection on different loci in similar functional pathways. In summary, we bring evidence for the important role of genetic redundancy in rapid adaptation involving traits with polygenic architecture.

Polyploidization as an opportunistic mutation: The role of unreduced gametes formation and genetic drift in polyploid establishment.

It is broadly assumed that polyploidy success reflects an increase in fitness associated with whole-genome duplication (WGD), due to higher tolerance to stressful conditions. Nevertheless, WGD also arises with several costs in neo-polyploid lineages, like genomic instability, or cellular mis-management. In addition to these costs, neo-polyploid individuals also face frequency dependent selection because of frequent low-fitness triploids formed by cross-ploidy pollinations when tetraploids are primarily rare in the population. Interestingly, the idea that polyploidy can be fixed by genetic drift as a neutral or deleterious mutation is currently underexplored in the literature. To test how and when polyploidy can fix in a population by chance, we built a theoretical model in which autopolyploidization occurs through the production of unreduced gametes, a trait modelled as a quantitative trait that is allowed to vary through time. We found that when tetraploid individuals are less or as fit as their diploid progenitors, fixation of polyploidy is only possible when genetic drift is stronger than natural selection. The necessity of drift for tetraploid fixation holds even when polyploidy confers a selective advantage, except for scenarios where tetraploids are much fitter than diploids. Finally, we found that self-fertilization is less beneficial for tetraploid establishment than previously thought, notably when polyploids harbour an initial decrease in fitness. Our results bring a novel, non-exclusive explanation for the unequal temporal and spatial distribution of polyploid species.

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.

Plant material selection, collection, preservation, and storage for nuclear DNA content estimation.

In theory, any plant tissue providing intact nuclei in sufficient quantity is suitable for nuclear DNA content estimation using flow cytometry (FCM). While this certainly opens a wide variety of possible applications of FCM, especially when compared to classical karyological techniques restricted to tissues with active cell division, tissue selection and quality may directly affect the precision (and sometimes even reliability) of FCM measurements. It is usually convenient to first consider the goals of the study to either aim for the highest possible accuracy of estimates (e.g., for inferring genome size, detecting homoploid intraspecific genome size variation, aneuploidy, among others), or to decide that histograms of reasonable resolution provide sufficient information (e.g., ploidy level screening within a single model species). Here, a set of best practices guidelines for selecting the optimal plant tissue for FCM analysis, sampling of material, and material preservation and storage are provided. In addition, factors potentially compromising the quality of FCM estimates of nuclear DNA content and data interpretation are discussed.

Genomic Signatures of Sexual Selection on Pollen-Expressed Genes in Arabis alpina.

Fertilization in angiosperms involves the germination of pollen on the stigma, followed by the extrusion of a pollen tube that elongates through the style and delivers two sperm cells to the embryo sac. Sexual selection could occur throughout this process when male gametophytes compete for fertilization. The strength of sexual selection during pollen competition should be affected by the number of genotypes deposited on the stigma. As increased self-fertilization reduces the number of mating partners, and the genetic diversity and heterozygosity of populations, it should thereby reduce the intensity of sexual selection during pollen competition. Despite the prevalence of mating system shifts, few studies have directly compared the molecular signatures of sexual selection during pollen competition in populations with different mating systems. Here we analyzed whole-genome sequences from natural populations of Arabis alpina, a species showing mating system variation across its distribution, to test whether shifts from cross- to self-fertilization result in molecular signatures consistent with sexual selection on genes involved in pollen competition. We found evidence for efficient purifying selection on genes expressed in vegetative pollen, and overall weaker selection on sperm-expressed genes. This pattern was robust when controlling for gene expression level and specificity. In agreement with the expectation that sexual selection intensifies under cross-fertilization, we found that the efficacy of purifying selection on male gametophyte-expressed genes was significantly stronger in genetically more diverse and outbred populations. Our results show that intra-sexual competition shapes the evolution of pollen-expressed genes, and that its strength fades with increasing self-fertilization rates.

Short- and long-term consequences of genome doubling: a meta-analysis.

PREMISE: Whole-genome duplication (WGD) is ubiquitous in plants. Recent reviews and meta-analyses, aiming to understand how such phenotypic transition could facilitate neopolyploid establishment, demonstrated multifarious immediate effects of WGD on fitness and reproductive traits. Yet, little is known about how short-term modifications evolve through time. Such a comparison among new and established polyploid lineages is crucial to understand which effects of WGD promote or impede polyploid survival. METHODS: We performed a meta-analysis to determine how WGD affects morphological, cellular, and fitness traits in autotetraploid individuals compared to their diploid progenitors. We studied how established tetraploids differed from diploids compared to neotetraploids, to further learn about the fate of WGD-associated phenotypic effects during polyploid establishment. RESULTS: The short-term effects of WGD were an increase in size of morphological traits and cells, accompanied by a decrease in fitness and the number of cells. After establishment, the morphological effect persisted, but cellular and fitness components reverted back to the values observed in the diploid ancestors. CONCLUSIONS: Our results suggest that the larger morphology of autotetraploids is not a constraint to establishment. However, other observable effects of genome doubling disappeared with time, suggesting that solving cellular and fitness constraints are critical aspects for polyploid establishment.

Parallel adaptation in autopolyploid Arabidopsis arenosa is dominated by repeated recruitment of shared alleles.

Relative contributions of pre-existing vs de novo genomic variation to adaptation are poorly understood, especially in polyploid organisms. We assess this in high resolution using autotetraploid Arabidopsis arenosa, which repeatedly adapted to toxic serpentine soils that exhibit skewed elemental profiles. Leveraging a fivefold replicated serpentine invasion, we assess selection on SNPs and structural variants (TEs) in 78 resequenced individuals and discover significant parallelism in candidate genes involved in ion homeostasis. We further model parallel selection and infer repeated sweeps on a shared pool of variants in nearly all these loci, supporting theoretical expectations. A single striking exception is represented by TWO PORE CHANNEL 1, which exhibits convergent evolution from independent de novo mutations at an identical, otherwise conserved site at the calcium channel selectivity gate. Taken together, this suggests that polyploid populations can rapidly adapt to environmental extremes, calling on both pre-existing variation and novel polymorphisms.

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.

Novelty and Convergence in Adaptation to Whole Genome Duplication.

Whole genome duplication (WGD) can promote adaptation but is disruptive to conserved processes, especially meiosis. Studies in Arabidopsis arenosa revealed a coordinated evolutionary response to WGD involving interacting proteins controlling meiotic crossovers, which are minimized in an autotetraploid (within-species polyploid) to avoid missegregation. Here, we test whether this surprising flexibility of a conserved essential process, meiosis, is recapitulated in an independent WGD system, Cardamine amara, 17 My diverged from A. arenosa. We assess meiotic stability and perform population-based scans for positive selection, contrasting the genomic response to WGD in C. amara with that of A. arenosa. We found in C. amara the strongest selection signals at genes with predicted functions thought important to adaptation to WGD: meiosis, chromosome remodeling, cell cycle, and ion transport. However, genomic responses to WGD in the two species differ: minimal ortholog-level convergence emerged, with none of the meiosis genes found in A. arenosa exhibiting strong signal in C. amara. This is consistent with our observations of lower meiotic stability and occasional clonal spreading in diploid C. amara, suggesting that nascent C. amara autotetraploid lineages were preadapted by their diploid lifestyle to survive while enduring reduced meiotic fidelity. However, in contrast to a lack of ortholog convergence, we see process-level and network convergence in DNA management, chromosome organization, stress signaling, and ion homeostasis processes. This gives the first insight into the salient adaptations required to meet the challenges of a WGD state and shows that autopolyploids can utilize multiple evolutionary trajectories to adapt to WGD.

Disentangling the components of triploid block and its fitness consequences in natural diploid-tetraploid contact zones of Arabidopsis arenosa.

Hybrid seed inviability (HSI) is an important mechanism of reproductive isolation and speciation. HSI varies in strength among populations of diploid species but it remains to be tested whether similar processes affect natural variation in HSI within ploidy-variable species (triploid block). Here we used extensive endosperm, seed and F1 -hybrid phenotyping to explore HSI variation within a diploid-autotetraploid species. By leveraging 12 population pairs from three ploidy contact zones, we tested for the effect of interploidy crossing direction (parent of origin), ploidy divergence and spatial arrangement in shaping reproductive barriers in a naturally relevant context. We detected strong parent-of-origin effects on endosperm development, F1 germination and survival, which was also reflected in the rates of triploid formation in the field. Endosperm cellularization failure was least severe and F1 -hybrid performance was slightly better in the primary contact zone, with genetically closest diploid and tetraploid lineages. We demonstrated overall strong parent-of-origin effects on HSI in a ploidy variable species, which translate to fitness effects and contribute to interploidy reproductive isolation in a natural context. Subtle intraspecific variation in these traits suggests the fitness consequences of HSI are predominantly a constitutive property of the species regardless of the evolutionary background of its populations.

Tracing evolutionary history and admixture in mixed-ploidy systems.

Polyploidy, resulting from whole genome duplication, is a widespread phenomenon throughout Eukaryotic kingdoms. It is estimated that 15% of speciation events in plants is due to polyploidization and 16% of plant species encompass ploidy variation. In spite of the evolutionary and economic significance of polyploidy, there is a limited set of tools that would allow routine population genetic and genomic analysis of polyploid systems, in particular for the inference of population diversity and differentiation from large genome-wide data sets (Dufresne et al., 2014; Meirmans, 2020). Such a shortage is striking especially when compared to the rapid development of such tools in diploid systems over the last decade. Consequently, population genomic research in polyploids is still lagging behind diploids, especially in autopolyploids, for example, polyploids with multiple (>2) similar genome copies. In this issue of Molecular Ecology Resources, Shastry et al. (2021) develops the first software specifically designed to infer population structure and ancestry in large genome-wide mixed-ploidy data sets, opening novel opportunities for polyploid population genomic analysis.

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.