The genomic basis of adaptation in plants

Genome-wide patterns of local adaptation associated with transposable elements in Tetrastigma hemsleyanum (Vitaceae)

The mobility of transposable elements (TEs) partly drives genome evolution, potentially leading to either adaptive or deleterious effects. However, it remains far from clear whether and how TEs contribute to adaptation to changing environments, especially in plants. We analyzed whole-genome sequencing data from 29 ecologically diverse Tetrastigma hemsleyanum populations to infer the species' demographic history and its impact on TE polymorphisms. Integrated selective sweep and genome-environment association (GEA) approaches were employed to examine the contribution of TEs to environmental adaptation. The ancestor of T. hemsleyanum diverged during the late Miocene/Pliocene, forming two lineages that further split into four sublineages. These (sub)lineages underwent periodic population declines and recoveries during the late-Pleistocene climatic oscillations, with most polymorphic TEs transposing during the last glacial period. A small fraction of these TEs (0.033-0.40%) showed signatures of positive selection, while a broader subset (0.081-0.76%) correlated significantly with climatic variables. Notably, these selected or climate-linked TE polymorphisms were preferentially retained in gene-poor regions and frequently linked to genes involved in organ development and stress/defense response. Our findings demonstrate that TEs played a key regulatory and adaptive role in T. hemsleyanum's response to environmental change, underscoring their importance in better understanding the genomic mechanisms underlying adaptation.

Comparative transcriptomics reveals potential regulators of climate adaptation in Santalum album L. (Indian Sandalwood)

Santalum album L. (Indian Sandalwood), a valued tree species known for its fragrant heartwood and essential oil is facing increasing threat due to severe anthropogenic pressures compounded by climate change which has resulted in depletion of its adaptive gene pool. The present study investigates the transcriptome-level responses of nine sandalwood genotypes sourced from diverse climatic zones to identify adaptive genes in the species. Comparative transcriptomics predicted 727, 1141 and 479 differentially expressed transcripts (DETs) across wet vs. dry; monsoon vs. dry and wet vs. monsoon conditions, respectively, and majority of DETs were up-regulated in samples sourced from high rainfall areas. Transcripts including heat shock proteins, Zinc finger binding protein, ribosomal proteins, transcription factors and protein kinase were identified as probable regulators of climate adaptation in S. album. The expression changes of eight selected transcripts were further validated by real-time quantitative PCR. Protein-protein interaction analysis revealed key hub transcripts involved in climate response, while alternative splicing events in transcripts such as SURP and G-patch domain-containing protein 1-like protein, G-type lectin S-receptor-like serine/threonine protein kinase B120, Tetraspanin-3 and ARM repeat superfamily protein indicated the probable role of alternate splicing in increasing the transcript diversity during adaptation. This study presents the first insight into the molecular mechanisms of climate adaptation in the species and can form the basis for specific interventions such as selective breeding, genetic manipulation, and habitat management for conservation and long-term survival of sandalwood.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-025-04218-4.

The Conditioned Environmental Center-Periphery Hypothesis of Biogeography: Statistical Evidence From Tree Species

It has been discussed for decades whether species occur most frequently at their geographic center, and more recently at their environmental niches' center. The aim here is to analyze for each environmental gradient separately the ecological niche of 12 Mexican tree species and 16 abiotic environmental gradients, in the form of statistical probability density functions. Is a symmetrically positioned center always possible by searching for additional data? For each species-variable combination, the occurrences along an environmental gradient were grouped in histograms. Logistic regression was used to fit a polynomial equation, whose degree depended on the number of significantly different bins. A highest-probability interval on the gradient was determined, where 25% of the individuals were found with the highest probability. The relative distance from the center (midpoint) of the variable's range was calculated, and the feasibility of expanding the encountered interval on the environmental gradient for symmetry was analyzed. For 183 species-variable combinations, in only 22 cases (12.0%) did the highest-probability intervals include the midpoint of the environmental gradient. Furthermore, for 55% of the species-variable combinations, the truncation of the environmental gradients for species makes it impossible to expand the measured range with additional data for the shorter tail. For example, precipitation cannot be negative. This truncation frequently causes asymmetry around the highest-probability intervals. In those cases, the classical environmental center-periphery hypothesis turns out to be wrong, whereas in the remaining cases it could apply. This has implications for biogeographical assumptions, such as where to identify the best areas for conservation or how to predict the effects of climate change. We propose a new conditioned environmental center-periphery hypothesis: "On an environmental gradient, a given species is able to cover a certain range. For environmental gradients, where natural truncation of the environmental gradient is not limiting, the highest probability of occurrence is found away from the range's endpoints, tending towards its midpoint."

Repeated global adaptation across plant species

Global adaptation occurs when all populations of a species undergo selection toward a common optimum. This can occur by a hard selective sweep with the emergence of a new globally advantageous allele that spreads throughout a species' natural range until reaching fixation. This evolutionary process leaves a temporary trace in the region affected, which is detectable using population genomic methods. While selective sweeps have been identified in many species, there have been few comparative and systematic studies of the genes involved in global adaptation. Building upon recent findings showing repeated genetic basis of local adaptation across independent populations and species, we asked whether certain genes play a more significant role in driving global adaptation across plant species. To address this question, we scanned the genomes of 17 plant species to identify signals of repeated global selective sweeps. Despite the substantial evolutionary distance between the species analyzed, we identified several gene families with strong evidence of repeated positive selection. These gene families tend to be enriched for reduced pleiotropy, consistent with predictions from Fisher's evolutionary model and the cost of complexity hypothesis. We also found that genes with repeated sweeps exhibit elevated levels of gene duplication. Our findings contrast with recent observations of increased pleiotropy in genes driving local adaptation, consistent with predictions based on the theory of migration-selection balance.

Multi-omics analysis reveals genetic architecture and local adaptation of coumarins metabolites in Populus

BACKGROUND

Accumulation of coumarins plays key roles in response to immune and abiotic stress in plants, but the genetic adaptation basis of controlling coumarins in perennial woody plants remain unclear.

RESULTS

We detected 792 SNPs within 334 genes that were significantly associated with the phenotypic variations of 15 single-metabolic traits and multiple comprehensive index, such as principal components (PCs) of coumarins metabolites. Expression quantitative trait locus mapping uncovered that 337 eQTLs associated with the expression levels of 132 associated genes. Selective sweep revealed 55 candidate genes have potential selective signature among three geographical populations, highlighting that the coumarins biosynthesis have been encountered forceful local adaptation. Furthermore, we constructed a genetic network of seven candidate genes that coordinately regulate coumarins biosynthesis, revealing the multiple regulatory patterns affecting coumarins accumulation in Populus tomentosa. Validation of candidate gene variations in a drought-tolerated population and DUF538 heterologous transformation experiments verified the function of candidate genes and their roles in adapting to the different geographical conditions in poplar.

CONCLUSIONS

Our study uncovered the genetic regulation of the coumarins metabolic biosynthesis of Populus, and offered potential clues for drought-tolerance evaluation and regional improvement in woody plants.

Genomic Insights into High-Altitude Adaptation: A Comparative Analysis of Roscoea alpina and R. purpurea in the Himalayas

Environmental stress at high altitudes drives the development of distinct adaptive mechanisms in plants. However, studies exploring the genetic adaptive mechanisms of high-altitude plant species are scarce. In the present study, we explored the high-altitude adaptive mechanisms of plants in the Himalayas through whole-genome resequencing. We studied two widespread members of the Himalayan endemic alpine genus Roscoea (Zingiberaceae): R. alpina (a selfing species) and R. purpurea (an outcrossing species). These species are distributed widely in the Himalayas with distinct non-overlapping altitude distributions; R. alpina is distributed at higher elevations, and R. purpurea occurs at lower elevations. Compared to R. purpurea, R. alpina exhibited higher levels of linkage disequilibrium, Tajima's D, and inbreeding coefficient, as well as lower recombination rates and genetic diversity. Approximately 96.3% of the genes in the reference genome underwent significant genetic divergence (FST ≥ 0.25). We reported 58 completely divergent genes (FST = 1), of which only 17 genes were annotated with specific functions. The functions of these genes were primarily related to adapting to the specific characteristics of high-altitude environments. Our findings provide novel insights into how evolutionary innovations promote the adaptation of mountain alpine species to high altitudes and harsh habitats.

Genomic insights into selection for heterozygous alleles and woody traits in Populus tomentosa

Heterozygous alleles are widespread in outcrossing and clonally propagated woody plants. The variation in heterozygosity that underlies population adaptive evolution and phenotypic variation, however, remains largely unknown. Here, we describe a de novo chromosome-level genome assembly of Populus tomentosa, an economic and ecologically important native tree in northern China. By resequencing 302 natural accessions, we determined that the South subpopulation (Pop_S) encompasses the ancestral strains of P. tomentosa, while the Northwest subpopulation (Pop_NW) and Northeast subpopulation (Pop_NE) experienced different selection pressures during population evolution, resulting in significant population differentiation and a decrease in the extent of heterozygosity. Analysis of heterozygous selective sweep regions (HSSR) suggested that selection for lower heterozygosity contributed to the local adaptation of P. tomentosa by dwindling gene expression and genetic load in the Pop_NW and Pop_NE subpopulations. Genome-wide association studies (GWAS) revealed that 88 single nucleotide polymorphisms (SNPs) within 63 genes are associated with nine wood composition traits. Among them, the selection for the homozygous AA allele in PtoARF8 is associated with reductions in cellulose and hemicellulose contents by attenuating PtoARF8 expression, and the increase in lignin content is attributable to the selection for decreases in exon heterozygosity in PtoLOX3 during adaptive evolution of natural populations. This study provides novel insights into allelic variations in heterozygosity associated with adaptive evolution of P. tomentosa in response to the local environment and identifies a series of key genes for wood component traits, thereby facilitating genomic-based breeding of important traits in perennial woody plants.

Towards the Investigation of the Adaptive Divergence in a Species of Exceptional Ecological Plasticity: Chromosome-Scale Genome Assembly of Chouardia litardierei (Hyacinthaceae)

One of the central goals of evolutionary biology is to understand the genomic basis of adaptive divergence. Different aspects of evolutionary processes should be studied through genome-wide approaches, therefore maximizing the investigated genomic space. However, in-depth genome-scale analyses often are restricted to a model or economically important species and their closely related wild congeners with available reference genomes. Here, we present the high-quality chromosome-level genome assembly of Chouardia litardierei, a plant species with exceptional ecological plasticity. By combining PacBio and Hi-C sequencing technologies, we generated a 3.7 Gbp genome with a scaffold N50 size of 210 Mbp. Over 80% of the genome comprised repetitive elements, among which the LTR retrotransposons prevailed. Approximately 86% of the 27,257 predicted genes were functionally annotated using public databases. For the comparative analysis of different ecotypes' genomes, the whole-genome sequencing of two individuals, each from a distinct ecotype, was performed. The detected above-average SNP density within coding regions suggests increased adaptive divergence-related mutation rates, therefore confirming the assumed divergence processes within the group. The constructed genome presents an invaluable resource for future research activities oriented toward the investigation of the genetics underlying the adaptive divergence that is likely unfolding among the studied species' ecotypes.

Fusarium equiseti-inoculation altered rhizosphere soil microbial community, potentially driving perennial ryegrass growth and salt tolerance

Fusarium equiseti is an effective plant growth-promoting fungi that induce systemic disease resistance in plants. However, the role of F. equiseti in regulating salt stress response and the underlying mechanisms remain largely unknown. Here, we investigated the effect of F. equiseti Z7 strain on the growth and salt stress response in perennial ryegrass. Additionally, the role of Z7 in regulating the abundance, composition, and structure of native microbial communities in the rhizosphere soil was determined. We observed that Z7 could produce indole-3-acetic acid (IAA) and siderophores. Hence, Z7 inoculation further enhanced plant growth and salt tolerance in perennial ryegrass. Inoculating Z7 increased K⁺ and decreased Na⁺ in plant tissues. Z7 inoculation also enhanced soil quality by reducing soluble salt and increasing available phosphorus. Moreover, inoculating Z7 altered the compositions of bacterial and fungal communities in the rhizosphere soil. For instance, beneficial bacterial genera, such as Flavobacterium, Enterobacter, Agrobacterium, and Burkholderiales were dominantly enriched in Z7-inoculated soil. Interestingly, the relative abundance of these genera showed significantly positive correlations with the fresh weight of perennial ryegrass. Our results demonstrate that Z7 could remarkably promote plant growth and salt tolerance by regulating ion homeostasis in plant tissues and microbial communities in the rhizosphere soil. This study provides a scientific foundation for applying microbes to improve plant growth under extreme salt stress conditions.

Genomic vulnerability to climate change in Quercus acutissima, a dominant tree species in East Asian deciduous forests

Understanding the evolutionary processes that shape the landscape of genetic variation and influence the response of species to future climate change is critical for biodiversity conservation. Here, we sampled 27 populations across the distribution range of a dominant forest tree, Quercus acutissima, in East Asia, and applied genome-wide analyses to track the evolutionary history and predict the fate of populations under future climate. We found two genetic groups (East and West) in Q. acutissima that diverged during Pliocene. We also found a heterogeneous landscape of genomic variation in this species, which may have been shaped by population demography and linked selections. Using genotype-environment association analyses, we identified climate-associated SNPs in a diverse set of genes and functional categories, indicating a model of polygenic adaptation in Q. acutissima. We further estimated three genetic offset metrics to quantify genomic vulnerability of this species to climate change due to the complex interplay between local adaptation and migration. We found that marginal populations are under higher risk of local extinction because of future climate change, and may not be able to track suitable habitats to maintain the gene-environment relationships observed under the current climate. We also detected higher reverse genetic offsets in northern China, indicating that genetic variation currently present in the whole range of Q. acutissima may not adapt to future climate conditions in this area. Overall, this study illustrates how evolutionary processes have shaped the landscape of genomic variation, and provides a comprehensive genome-wide view of climate maladaptation in Q. acutissima.

Genetic diversity goals and targets have improved, but remain insufficient for clear implementation of the post-2020 global biodiversity framework

Genetic diversity among and within populations of all species is necessary for people and nature to survive and thrive in a changing world. Over the past three years, commitments for conserving genetic diversity have become more ambitious and specific under the Convention on Biological Diversity's (CBD) draft post-2020 global biodiversity framework (GBF). This Perspective article comments on how goals and targets of the GBF have evolved, the improvements that are still needed, lessons learned from this process, and connections between goals and targets and the actions and reporting that will be needed to maintain, protect, manage and monitor genetic diversity. It is possible and necessary that the GBF strives to maintain genetic diversity within and among populations of all species, to restore genetic connectivity, and to develop national genetic conservation strategies, and to report on these using proposed, feasible indicators.

Gene family evolution and natural selection signatures in Datura spp. (Solanaceae)

Elucidating the diversification process of congeneric species makes it necessary to identify the factors promoting species variation and diversification. Comparative gene family analysis allows us to elucidate the evolutionary history of species by identifying common genetic/genomic mechanisms underlying species responses to biotic and abiotic environments at the genomic level. In this study, we analyzed the high-quality transcriptomes of four Datura species, D. inoxia, D. pruinosa, D. stramonium, and D. wrightii. We performed a thorough comparative gene family analysis to infer the role of selection in molecular variation, changes in protein physicochemical properties, and gain/loss of genes during their diversification processes. The results revealed common and species-specific signals of positive selection, physicochemical divergence and/or expansion of metabolic genes (e.g., transferases and oxidoreductases) associated with terpene and tropane metabolism and some resistance genes (R genes). The gene family analysis presented here is a valuable tool for understanding the genome evolution of economically and ecologically significant taxa such as the Solanaceae family.

Natural variation in root suberization is associated with local environment in Arabidopsis thaliana

Genetic signature of climate adaptation has been widely recognized across the genome of many organisms; however, the eco-physiological basis for linking genomic polymorphisms with local adaptations remains largely unexplored. Using a panel of 218 world-wide Arabidopsis accessions, we characterized the natural variation in root suberization by quantifying 16 suberin monomers. We explored the associations between suberization traits and 126 climate variables. We conducted genome-wide association analysis and integrated previous genotype-environment association (GEA) to identify the genetic bases underlying suberization variation and their involvements in climate adaptation. Root suberin content displays extensive variation across Arabidopsis populations and significantly correlates with local moisture gradients and soil characteristics. Specifically, enhanced suberization is associated with drier environments, higher soil cation-exchange capacity, and lower soil pH; higher proportional levels of very-long-chain suberin is negatively correlated with moisture availability, lower soil gravel content, and higher soil silt fraction. We identified 94 putative causal loci and experimentally proved that GPAT6 is involved in C16 suberin biosynthesis. Highly significant associations between the putative genes and environmental variables were observed. Roots appear highly responsive to environmental heterogeneity via regulation of suberization, especially the suberin composition. The patterns of suberization-environment correlation and the suberin-related GEA fit the expectations of local adaptation for the polygenic suberization trait.

Evolution in agricultural systems: Moving toward the understanding of complexity

Agricultural fields are typically simplified ecosystems compared to natural sites, a characteristic that has long-attracted researchers in Ecology and Evolution. In recent years, there has been a rising interest in understanding how agricultural systems are shaped by evolution in the context of changing agricultural practices by integrating biological information of crop systems. This editorial introduces the special issue "Evolution in agricultural systems," incorporating the articles published within this issue into three general areas of research: phenotypic and genetic responses to the environment, biotic interactions and the role of microbes. Together, this body of work unveils unforeseen complexity at all levels, from microbes to trophic chains. Understanding such complexity is critical not only to better understand natural systems, but also if we wish to improve the sustainability of the food system.

Hard versus soft selective sweeps during domestication and improvement in soybean

Genome scans for selection can provide an efficient way to dissect the genetic basis of domestication traits and understand mechanisms of adaptation during crop evolution. Selection involving soft sweeps (simultaneous selection for multiple alleles) is probably common in plant genomes but is under-studied, and few if any studies have systematically scanned for soft sweeps in the context of crop domestication. Using genome resequencing data from 302 wild and domesticated soybean accessions, we conducted selection scans using five widely employed statistics to identify selection candidates under classical (hard) and soft sweeps. Across the genome, inferred hard sweeps are predominant in domesticated soybean landraces and improved varieties, whereas soft sweeps are more prevalent in a representative subpopulation of the wild ancestor. Six domestication-related genes, representing both hard and soft sweeps and different stages of domestication, were used as positive controls to assess the detectability of domestication-associated sweeps. Performance of various test statistics suggests that differentiation-based (FST ) methods are robust for detecting complete hard sweeps, and that LD-based strategies perform well for identifying recent/ongoing sweeps; however, none of the test statistics detected a known soft sweep we previously documented at the domestication gene Dt1. Genome scans yielded a set of 66 candidate loci that were identified by both differentiation-based and LD-based (iHH) methods; notably, this shared set overlaps with many previously identified QTLs for soybean domestication/improvement traits. Collectively, our results will help to advance genetic characterizations of soybean domestication traits and shed light on selection modes involved in adaptation in domesticated plant species.

Evolutionary Ecology of Plant-Arthropod Interactions in Light of the "Omics" Sciences: A Broad Guide

Aboveground plant-arthropod interactions are typically complex, involving herbivores, predators, pollinators, and various other guilds that can strongly affect plant fitness, directly or indirectly, and individually, synergistically, or antagonistically. However, little is known about how ongoing natural selection by these interacting guilds shapes the evolution of plants, i.e., how they affect the differential survival and reproduction of genotypes due to differences in phenotypes in an environment. Recent technological advances, including next-generation sequencing, metabolomics, and gene-editing technologies along with traditional experimental approaches (e.g., quantitative genetics experiments), have enabled far more comprehensive exploration of the genes and traits involved in complex ecological interactions. Connecting different levels of biological organization (genes to communities) will enhance the understanding of evolutionary interactions in complex communities, but this requires a multidisciplinary approach. Here, we review traditional and modern methods and concepts, then highlight future avenues for studying the evolution of plant-arthropod interactions (e.g., plant-herbivore-pollinator interactions). Besides promoting a fundamental understanding of plant-associated arthropod communities' genetic background and evolution, such knowledge can also help address many current global environmental challenges.

Globally Relaxed Selection and Local Adaptation in Boechera stricta

The strength of selection varies among populations and across the genome, but the determinants of efficacy of selection remain unclear. In this study, we used whole-genome sequencing data from 467 Boechera stricta accessions to quantify the strength of selection and characterize the pattern of local adaptation. We found low genetic diversity on 0-fold degenerate sites and conserved non-coding sites, indicating functional constraints on these regions. The estimated distribution of fitness effects and the proportion of fixed substitutions suggest relaxed negative and positive selection in B. stricta. Among the four population groups, the NOR and WES groups have smaller effective population size (Ne), higher proportions of effectively neutral sites, and lower rates of adaptive evolution compared with UTA and COL groups, reflecting the effect of Ne on the efficacy of natural selection. We also found weaker selection on GC-biased sites compared with GC-conservative (unbiased) sites, suggested that GC-biased gene conversion has affected the strength of selection in B. stricta. We found mixed evidence for the role of the recombination rate on the efficacy of selection. The positive and negative selection was stronger in high-recombination regions compared with low-recombination regions in COL but not in other groups. By scanning the genome, we found different subsets of selected genes suggesting differential adaptation among B. stricta groups. These results show that differences in effective population size, nucleotide composition, and recombination rate are important determinants of the efficacy of selection. This study enriches our understanding of the roles of natural selection and local adaptation in shaping genomic variation.

Climate Adaptation, Drought Susceptibility, and Genomic-Informed Predictions of Future Climate Refugia for the Australian Forest Tree Eucalyptus globulus

Understanding the capacity of forest tree species to adapt to climate change is of increasing importance for managing forest genetic resources. Through a genomics approach, we modelled spatial variation in climate adaptation within the Australian temperate forest tree Eucalyptus globulus, identified putative climate drivers of this genomic variation, and predicted locations of future climate refugia and populations at-risk of future maladaptation. Using 812,158 SNPs across 130 individuals from 30 populations (i.e., localities) spanning the species’ natural range, a gradientForest algorithm found 1177 SNPs associated with locality variation in home-site climate (climate-SNPs), putatively linking them to climate adaptation. Very few climate-SNPs were associated with population-level variation in drought susceptibility, signalling the multi-faceted nature and complexity of climate adaptation. Redundancy analysis (RDA) showed 24% of the climate-SNP variation could be explained by annual precipitation, isothermality, and maximum temperature of the warmest month. Spatial predictions of the RDA climate vectors associated with climate-SNPs allowed mapping of genomically informed climate selective surfaces across the species’ range under contemporary and projected future climates. These surfaces suggest over 50% of the current distribution of E. globulus will be outside the modelled adaptive range by 2070 and at risk of climate maladaptation. Such surfaces present a new integrated approach for natural resource managers to capture adaptive genetic variation and plan translocations in the face of climate change.

Linked selection shapes the landscape of genomic variation in three oak species

Natural selection shapes genome-wide patterns of diversity within species and divergence between species. However, quantifying the efficacy of selection and elucidating the relative importance of different types of selection in shaping genomic variation remain challenging. We sequenced whole genomes of 101 individuals of three closely related oak species to track the divergence history, and to dissect the impacts of selective sweeps and background selection on patterns of genomic variation. We estimated that the three species diverged around the late Neogene and experienced a bottleneck during the Pleistocene. We detected genomic regions with elevated relative differentiation ('F_ST -islands'). Population genetic inferences from the site frequency spectrum and ancestral recombination graph indicated that F_ST -islands were formed by selective sweeps. We also found extensive positive selection; the fixation of adaptive mutations and reduction neutral diversity around substitutions generated a signature of selective sweeps. Prevalent negative selection and background selection have reduced genetic diversity in both genic and intergenic regions, and contributed substantially to the baseline variation in genetic diversity. Our results demonstrate the importance of linked selection in shaping genomic variation, and illustrate how the extent and strength of different selection models vary across the genome.

Salt-induced recruitment of specific root-associated bacterial consortium capable of enhancing plant adaptability to salt stress

Salinity is a major abiotic stress threatening crop production. Root-derived bacteria (RDB) are hypothesized to play a role in enhancing plant adaptability to various stresses. However, it is still unclear whether and how plants build up specific RDB when challenged by salinity. In this study, we measured the composition and variation in the rhizosphere and endophyte bacteria of salt-sensitive (SSs) and salt-resistant (SRs) plants under soil conditions with/without salinity. The salt-induced RDB (both rhizobiomes and endophytes) were isolated to examine their effects on the physiological responses of SSs and SRs to salinity challenge. Moreover, we examined whether functional redundancy exists among salt-induced RDB in enhancing plant adaptability to salt stress. We observed that although SSs and SRs recruited distinct RDB and relevant functions when challenged by salinity, salt-induced recruitment of specific RDB led to a consistent growth promotion in plants regardless of their salinity tolerance capacities. Plants employed a species-specific strategy to recruit beneficial soil bacteria in the rhizosphere rather than in the endosphere. Furthermore, we demonstrated that the consortium, but not individual members of the salt-induced RDB, provided enduring resistance against salt stress. This study confirms the critical role of salt-induced RDB in enhancing plant adaptability to salt stress.

Evidence of local adaptation despite strong drift in a Neotropical patchily distributed bromeliad

Both genetic drift and divergent selection are predicted to be drivers of population differentiation across patchy habitats, but the extent to which these forces act on natural populations to shape traits is strongly affected by species' ecological features. In this study, we infer the genomic structure of Pitcairnia lanuginosa, a widespread herbaceous perennial plant with a patchy distribution. We sampled populations in the Brazilian Cerrado and the Central Andean Yungas and discovered and genotyped SNP markers using double-digest restriction-site associated DNA sequencing. In addition, we analyzed ecophysiological traits obtained from a common garden experiment and compared patterns of phenotypic and genetic divergence (P_ST-F_ST comparisons) in a subset of populations from the Cerrado. Our results from molecular analyses pointed to extremely low genetic diversity and a remarkable population differentiation, supporting a major role of genetic drift. Approximately 0.3% of genotyped SNPs were flagged as differentiation outliers by at least two distinct methods, and Bayesian generalized linear mixed models revealed a signature of isolation by environment in addition to isolation by distance for high-differentiation outlier SNPs among the Cerrado populations. P_ST-F_ST comparisons suggested divergent selection on two ecophysiological traits linked to drought tolerance. We showed that these traits vary among populations, although without any particular macro-spatial pattern, suggesting local adaptation to differences in micro-habitats. Our study shows that selection might be a relevant force, particularly for traits involved in drought stress, even for populations experiencing strong drift, which improves our knowledge on eco-evolutionary processes acting on non-continuously distributed species.

Population transcriptomic sequencing reveals allopatric divergence and local adaptation in Pseudotaxus chienii (Taxaceae)

BACKGROUND

Elucidating the effects of geography and selection on genetic variation is critical for understanding the relative importance of adaptation in driving differentiation and identifying the environmental factors underlying its occurrence. Adaptive genetic variation is common in tree species, especially widely distributed long-lived species. Pseudotaxus chienii can occupy diverse habitats with environmental heterogeneity and thus provides an ideal material for investigating the process of population adaptive evolution. Here, we characterize genetic and expression variation patterns and investigate adaptive genetic variation in P. chienii populations.

RESULTS

We generated population transcriptome data and identified 13,545 single nucleotide polymorphisms (SNPs) in 5037 unigenes across 108 individuals from 10 populations. We observed lower nucleotide diversity (π = 0.000701) among the 10 populations than observed in other gymnosperms. Significant negative correlations between expression diversity and nucleotide diversity in eight populations suggest that when the species adapts to the surrounding environment, gene expression and nucleotide diversity have a reciprocal relationship. Genetic structure analyses indicated that each distribution region contains a distinct genetic group, with high genetic differentiation among them due to geographical isolation and local adaptation. We used FST outlier, redundancy analysis, and latent factor mixed model methods to detect molecular signatures of local adaptation. We identified 244 associations between 164 outlier SNPs and 17 environmental variables. The mean temperature of the coldest quarter, soil Fe and Cu contents, precipitation of the driest month, and altitude were identified as the most important determinants of adaptive genetic variation. Most candidate unigenes with outlier signatures were related to abiotic and biotic stress responses, and the monoterpenoid biosynthesis and ubiquitin-mediated proteolysis KEGG pathways were significantly enriched in certain populations and deserve further attention in other long-lived trees.

CONCLUSIONS

Despite the strong population structure in P. chienii, genomic data revealed signatures of divergent selection associated with environmental variables. Our research provides SNPs, candidate unigenes, and biological pathways related to environmental variables to facilitate elucidation of the genetic variation in P. chienii in relation to environmental adaptation. Our study provides a promising tool for population genomic analyses and insights into the molecular basis of local adaptation.

Pests, diseases, and aridity have shaped the genome of Corymbia citriodora

Corymbia citriodora is a member of the predominantly Southern Hemisphere Myrtaceae family, which includes the eucalypts (Eucalyptus, Corymbia and Angophora; ~800 species). Corymbia is grown for timber, pulp and paper, and essential oils in Australia, South Africa, Asia, and Brazil, maintaining a high-growth rate under marginal conditions due to drought, poor-quality soil, and biotic stresses. To dissect the genetic basis of these desirable traits, we sequenced and assembled the 408 Mb genome of Corymbia citriodora, anchored into eleven chromosomes. Comparative analysis with Eucalyptus grandis reveals high synteny, although the two diverged approximately 60 million years ago and have different genome sizes (408 vs 641 Mb), with few large intra-chromosomal rearrangements. C. citriodora shares an ancient whole-genome duplication event with E. grandis but has undergone tandem gene family expansions related to terpene biosynthesis, innate pathogen resistance, and leaf wax formation, enabling their successful adaptation to biotic/abiotic stresses and arid conditions of the Australian continent.

High-depth resequencing of 312 accessions reveals the local adaptation of foxtail millet

Based on the high-density variation map, we identified genome-level evidence for local adaptation and demonstrated that Siprr37 with transposon insertion contributes to the fitness of foxtail millet in the northeastern ecoregion. Adaptation is a robust way through which plants are able to overcome environmental constraints. The mechanisms of adaptation in heterogeneous natural environments are largely unknown. Deciphering the genomic basis of local adaptation will contribute to further improvement in domesticated plants. To this end, we describe a high-depth (19.4 ×) haplotype map of 3.02 million single nucleotide polymorphisms in foxtail millet (Setaria italica) from whole-genome resequencing of 312 accessions. In the genome-wide scan, we identified a set of improvement signals (including the homologous gene of OsIPA1, a key gene controlling ideal plant architecture) related to the geographical adaptation to four ecoregions in China. In particular, based on the genome-wide association analysis results, we identified the contribution of a pseudo-response regulator gene, SiPRR37, to heading date adaptation in foxtail millet. We observed the expression changes of SiPRR37 resulted from a key Tc1-Mariner transposon insertion in the first intron. Positive selection analyses revealed that SiPRR37 mainly contributed to the adaptation of northeastern ecoregions. Taken together, foxtail millet adapted to the northeastern region by regulating the function of SiPRR37, which sheds lights on genome-level evidence for adaptive geographical divergence. Besides, our data provide a nearly complete catalog of genomic variation aiding the identification of functionally important variants.

Locally adapted gut microbiomes mediate host stress tolerance

While evidence for the role of the microbiome in shaping host stress tolerance is becoming well-established, to what extent this depends on the interaction between the host and its local microbiome is less clear. Therefore, we investigated whether locally adapted gut microbiomes affect host stress tolerance. In the water flea Daphnia magna, we studied if the host performs better when receiving a microbiome from their source region than from another region when facing a stressful condition, more in particular exposure to the toxic cyanobacteria Microcystis aeruginosa. Therefore, a reciprocal transplant experiment was performed in which recipient, germ-free D. magna, isolated from different ponds, received a donor microbiome from sympatric or allopatric D. magna that were pre-exposed to toxic cyanobacteria or not. We tested for effects on host life history traits and gut microbiome composition. Our data indicate that Daphnia interact with particular microbial strains mediating local adaptation in host stress tolerance. Most recipient D. magna individuals performed better when inoculated with sympatric than with allopatric microbiomes. This effect was most pronounced when the donors were pre-exposed to the toxic cyanobacteria, but this effect was also pond and genotype dependent. We discuss how this host fitness benefit is associated with microbiome diversity patterns.

Review: Plant eco-evolutionary responses to climate change: Emerging directions

Contemporary climate change is exposing plant populations to novel combinations of temperatures, drought stress, [CO₂] and other abiotic and biotic conditions. These changes are rapidly disrupting the evolutionary dynamics of plants. Despite the multifactorial nature of climate change, most studies typically manipulate only one climatic factor. In this opinion piece, we explore how climate change factors interact with each other and with biotic pressures to alter evolutionary processes. We evaluate the ramifications of climate change across life history stages,and examine how mating system variation influences population persistence under rapid environmental change. Furthermore, we discuss how spatial and temporal mismatches between plants and their mutualists and antagonists could affect adaptive responses to climate change. For example, plant-virus interactions vary from highly pathogenic to mildly facilitative, and are partly mediated by temperature, moisture availability and [CO₂]. Will host plants exposed to novel, stressful abiotic conditions be more susceptible to viral pathogens? Finally, we propose novel experimental approaches that could illuminate how plants will cope with unprecedented global change, such as resurrection studies combined with experimental evolution, genomics or epigenetics.

Temperature, rainfall and wind variables underlie environmental adaptation in natural populations of Drosophila melanogaster

While several studies in a diverse set of species have shed light on the genes underlying adaptation, our knowledge on the selective pressures that explain the observed patterns lags behind. Drosophila melanogaster is a valuable organism to study environmental adaptation because this species originated in Southern Africa and has recently expanded worldwide, and also because it has a functionally well-annotated genome. In this study, we aimed to decipher which environmental variables are relevant for adaptation of D. melanogaster natural populations in Europe and North America. We analysed 36 whole-genome pool-seq samples of D. melanogaster natural populations collected in 20 European and 11 North American locations. We used the BayPass software to identify single nucleotide polymorphisms (SNPs) and transposable elements (TEs) showing signature of adaptive differentiation across populations, as well as significant associations with 59 environmental variables related to temperature, rainfall, evaporation, solar radiation, wind, daylight hours, and soil type. We found that in addition to temperature and rainfall, wind related variables are also relevant for D. melanogaster environmental adaptation. Interestingly, 23%-51% of the genes that showed significant associations with environmental variables were not found overly differentiated across populations. In addition to SNPs, we also identified 10 reference transposable element insertions associated with environmental variables. Our results showed that genome-environment association analysis can identify adaptive genetic variants that are undetected by population differentiation analysis while also allowing the identification of candidate environmental drivers of adaptation.

Genomic Analyses of Phenotypic Differences Between Native and Invasive Populations of Diffuse Knapweed (Centaurea diffusa)

Invasive species represent excellent opportunities to study the evolutionary potential of traits important to success in novel environments. Although some ecologically important traits have been identified in invasive species, little is typically known about the genetic mechanisms that underlie invasion success in non-model species. Here, we use a genome-wide association (GWAS) approach to identify the genetic basis of trait variation in the non-model, invasive, diffuse knapweed [Centaurea diffusa Lam. (Asteraceae)]. To assist with this analysis, we have assembled the first draft genome reference and fully annotated plastome assembly for this species, and one of the first from this large, weedy, genus, which is of major ecological and economic importance. We collected phenotype data from 372 individuals from four native and four invasive populations of C. diffusa grown in a common environment. Using these individuals, we produced reduced-representation genotype-by-sequencing (GBS) libraries and identified 7,058 SNPs. We identify two SNPs associated with leaf width in these populations, a trait which significantly varies between native and invasive populations. In this rosette forming species, increased leaf width is a major component of increased biomass, a common trait in invasive plants correlated with increased fitness. Finally, we use annotations from Arabidopsis thaliana to identify 98 candidate genes that are near the associated SNPs and highlight several good candidates for leaf width variation.

Microgeographical adaptation corresponds to elevational distributions of congeneric montane grasshoppers

Local adaptation can occur at small spatial scales relative to the dispersal capacity of species. Alpine ecosystems have sharp environmental clines that offer an opportunity to investigate the effects of fine-scale shifts in species' niche breadth on adaptive genetic processes. Here we examine two grasshopper species endemic to the Australian Alps (Kosciuscola spp.) that differ in elevational niche breadth: one broader, K. usitatus (1400-2200 m), and one narrower, K. tristis (1600-2000 m). We examine signatures of selection with respect to environmental and morphological variables in two mountain regions using F_ST outlier tests and environmental association analyses (EAAs) applied to single nucleotide polymorphism (SNP) data (K. usitatus: 9017 SNPs, n = 130; K. tristis: 7363 SNPs, n = 135). Stronger genetic structure was found in the more narrowly distributed K. tristis, which showed almost twice the number of SNPs under putative selection (10.8%) compared with K. usitatus (5.3%). When examining SNPs in common across species (n = 3058), 260 SNPs (8.5%) were outliers shared across species, and these were mostly associated with elevation, a proxy for temperature, suggesting parallel adaptive processes in response to climatic drivers. Additive polygenic scores (an estimate of the cumulative signal of selection across all candidate loci) were nonlinearly and positively correlated with elevation in both species. However, a steeper correlation in K. tristis indicated a stronger signal of spatially varying selection towards higher elevations. Our study illustrates that the niche breadth of co-occurring and related species distributed along the same environmental cline is associated with differences in patterns of microgeographical adaptation.

Genomic and chemical evidence for local adaptation in resistance to different herbivores in Datura stramonium

Because most species are collections of genetically variable populations distributed to habitats differing in their abiotic/biotic environmental factors and community composition, the pattern and strength of natural selection imposed by species on each other's traits are also expected to be highly spatially variable. Here, we used genomic and quantitative genetic approaches to understand how spatially variable selection operates on the genetic basis of plant defenses to herbivores. To this end, an F₂ progeny was generated by crossing Datura stramonium (Solanaceae) parents from two populations differing in their level of chemical defense. This F₂ progeny was reciprocally transplanted into the parental plants' habitats and by measuring the identity by descent (IBD) relationship of each F₂ plant to each parent, we were able to elucidate how spatially variable selection imposed by herbivores operated on the genetic background (IBD) of resistance to herbivory, promoting local adaptation. The results highlight that plants possessing the highest total alkaloid concentrations (sum of all alkaloid classes) were not the most well-defended or fit. Instead, specific alkaloids and their linked loci/alleles were favored by selection imposed by different herbivores. This has led to population differentiation in plant defenses and thus, to local adaptation driven by plant-herbivore interactions.

Genomic insight into the developmental history of southern highbush blueberry populations

Interspecific hybridization is a common breeding approach for introducing novel traits and genetic diversity to breeding populations. Southern highbush blueberry (SHB) is a blueberry cultivar group that has been intensively bred over the last 60 years. Specifically, it was developed by multiple interspecific crosses between northern highbush blueberry [NHB, Vaccinium corymbosum L. (2n = 4x = 48)] and low-chill Vaccinium species to expand the geographic limits of highbush blueberry production. In this study, we genotyped polyploid blueberries, including 105 SHB, 17 NHB, and 10 rabbiteye blueberry (RE) (Vaccinium virgatum Aiton), from the accessions planted at Poplarville, Mississippi, and accessions distributed in Japan, based on the double-digest restriction site-associated DNA sequencing. The genome-wide SNP data clearly indicated that RE cultivars were genetically distinct from SHB and NHB cultivars, whereas NHB and SHB were genetically indistinguishable. The population structure results appeared to reflect the differences in the allele selection strategies that breeders used for developing germplasm adapted to local climates. The genotype data implied that there are no or very few genomic segments that were commonly introgressed from low-chill Vaccinium species to the SHB genome. Principal component analysis-based outlier detection analysis found a few loci associated with a variable that could partially differentiate NHB and SHB. These SNP loci were detected in Mb-scale haplotype blocks and may be close to the functional genes related to SHB development. Collectively, the data generated in this study suggest a polygenic adaptation of SHB to the southern climate, and may be relevant for future population-scale genome-wide analyses of blueberry.

Ten Years of the Maize Nested Association Mapping Population: Impact, Limitations, and Future Directions

It has been just over a decade since the release of the maize (Zea mays) Nested Association Mapping (NAM) population. The NAM population has been and continues to be an invaluable resource for the maize genetics community and has yielded insights into the genetic architecture of complex traits. The parental lines have become some of the most well-characterized maize germplasm, and their de novo assemblies were recently made publicly available. As we enter an exciting new stage in maize genomics, this retrospective will summarize the design and intentions behind the NAM population; its application, the discoveries it has enabled, and its influence in other systems; and use the past decade of hindsight to consider whether and how it will remain useful in a new age of genomics.

Genetic Variation Related to High Elevation Adaptation Revealed by Common Garden Experiments in Pinus yunnanensis

Local adaptation, adaptation to specialized niches and environmental clines have been extensively reported for forest trees. Investigation of the adaptive genetic variation is crucial for forest resource management and breeding, especially in the context of global climate change. Here, we utilized a Pinus yunnanensis common garden experiments established at high and low elevation sites to assess the differences in growth and survival among populations and between the two common garden sites. The studied traits showed significant variation between the two test sites and among populations, suggesting adaptive divergence. To detect genetic variation related to environment, we captured 103,608 high quality SNPs based on RNA sequencing, and used them to assess the genetic diversity and population structure. We identified 321 outlier SNPs from 131 genes showing significant divergence in allelic frequency between survival populations of two sites. Functional categories associated with adaptation to high elevation were found to be related to flavonoid biosynthesis, response to UV, DNA repair, response to reactive oxygen species, and membrane lipid metabolic process. Further investigation of the outlier genes showed overrepresentation of the flavonoid biosynthesis pathway, suggesting that this pathway may play a key role in P. yunnanensis adaptation to high elevation environments. The outlier genes identified, and their variants, provide a basic reference for advanced investigations.

Evolutionary genomics can improve prediction of species' responses to climate change

Global climate change (GCC) increasingly threatens biodiversity through the loss of species, and the transformation of entire ecosystems. Many species are challenged by the pace of GCC because they might not be able to respond fast enough to changing biotic and abiotic conditions. Species can respond either by shifting their range, or by persisting in their local habitat. If populations persist, they can tolerate climatic changes through phenotypic plasticity, or genetically adapt to changing conditions depending on their genetic variability and census population size to allow for de novo mutations. Otherwise, populations will experience demographic collapses and species may go extinct. Current approaches to predicting species responses to GCC begin to combine ecological and evolutionary information for species distribution modelling. Including an evolutionary dimension will substantially improve species distribution projections which have not accounted for key processes such as dispersal, adaptive genetic change, demography, or species interactions. However, eco-evolutionary models require new data and methods for the estimation of a species' adaptive potential, which have so far only been available for a small number of model species. To represent global biodiversity, we need to devise large-scale data collection strategies to define the ecology and evolutionary potential of a broad range of species, especially of keystone species of ecosystems. We also need standardized and replicable modelling approaches that integrate these new data to account for eco-evolutionary processes when predicting the impact of GCC on species' survival. Here, we discuss different genomic approaches that can be used to investigate and predict species responses to GCC. This can serve as guidance for researchers looking for the appropriate experimental setup for their particular system. We furthermore highlight future directions for moving forward in the field and allocating available resources more effectively, to implement mitigation measures before species go extinct and ecosystems lose important functions.

Postglacial ecotype formation under outcrossing and self-fertilization in Arabidopsis lyrata

The formation of ecotypes has been invoked as an important driver of postglacial biodiversity, because many species colonized heterogeneous habitats and experienced divergent selection. Ecotype formation has been predominantly studied in outcrossing taxa, while far less attention has been paid to the implications of mating system shifts. Here, we addressed whether substrate-related ecotypes exist in selfing and outcrossing populations of Arabidopsis lyrata subsp. lyrata and whether the genomic footprint differs between mating systems. The North American subspecies colonized both rocky and sandy habitats during postglacial range expansion and shifted the mating system from predominantly outcrossing to predominantly selfing in a number of regions. We performed an association study on pooled whole-genome sequence data of 20 selfing or outcrossing populations, which suggested genes involved in adaptation to substrate. Motivated by enriched gene ontology terms, we compared root growth between plants from the two substrates in a common environment and found that plants originating from sand grew roots faster and produced more side roots, independent of mating system. Furthermore, single nucleotide polymorphisms associated with substrate-related ecotypes were more clustered among selfing populations. Our study provides evidence for substrate-related ecotypes in A. lyrata and divergence in the genomic footprint between mating systems. The latter is the likely result of selfing populations having experienced divergent selection on larger genomic regions due to higher genome-wide linkage disequilibrium.

Using natural variation to understand the evolutionary pressures on plant photosynthesis

Photosynthesis is the gateway of the Sun's energy into the biosphere and the source of the ozone layer; thus it is both provider and protector of life as we know it. Despite its pivotal role we know surprisingly little about the genetic basis of variation in photosynthesis and the selective pressures giving rise to or maintaining this variation. In this review, I will briefly summarise our current knowledge of intraspecific and interspecific variation in photosynthesis to understand the main selective constraints on photosynthesis and what this means for the future of nature and agriculture in a changing world.

WhoGEM: an admixture-based prediction machine accurately predicts quantitative functional traits in plants

The explosive growth of genomic data provides an opportunity to make increased use of sequence variations for phenotype prediction. We have developed a prediction machine for quantitative phenotypes (WhoGEM) that overcomes some of the bottlenecks limiting the current methods. We demonstrated its performance by predicting quantitative disease resistance and quantitative functional traits in the wild model plant species, Medicago truncatula, using geographical locations as covariates for admixture analysis. The method's prediction reliability equals or outperforms all existing algorithms for quantitative phenotype prediction. WhoGEM analysis produces evidence that variation in genome admixture proportions explains most of the phenotypic variation for quantitative phenotypes.

Resolving Form–Structure–Function Relationships in Plants with MRI for Biomimetic Transfer

In many biomimetic approaches, a deep understanding of the form–structure–function relationships in living and functionally intact organisms, which act as biological role models, is essential. This knowledge is a prerequisite for the identification of parameters that are relevant for the desired technical transfer of working principles. Hence, non-invasive and non-destructive techniques for static (3D) and dynamic (4D) high-resolution plant imaging and analysis on multiple hierarchical levels become increasingly important. In this study we demonstrate that magnetic resonance imaging (MRI) can be used to resolve the plants inner tissue structuring and functioning on the example of four plant concept generators with sizes larger than 5 mm used in current biomimetic research projects: Dragon tree (Dracaena reflexa var. angustifolia), Venus flytrap (Dionaea muscipula), Sugar pine (Pinus lambertiana) and Chinese witch hazel (Hamamelis mollis). Two different MRI sequences were applied for high-resolution 3D imaging of the differing material composition (amount, distribution, and density of various tissues) and condition (hydrated, desiccated, and mechanically stressed) of the four model organisms. Main aim is to better understand their biomechanics, development, and kinematics. The results are used as inspiration for developing novel design and fabrication concepts for bio-inspired technical fiber-reinforced branchings and smart biomimetic actuators.

Stress response, behavior, and development are shaped by transposable element-induced mutations in Drosophila

Most of the current knowledge on the genetic basis of adaptive evolution is based on the analysis of single nucleotide polymorphisms (SNPs). Despite increasing evidence for their causal role, the contribution of structural variants to adaptive evolution remains largely unexplored. In this work, we analyzed the population frequencies of 1,615 Transposable Element (TE) insertions annotated in the reference genome of Drosophila melanogaster, in 91 samples from 60 worldwide natural populations. We identified a set of 300 polymorphic TEs that are present at high population frequencies, and located in genomic regions with high recombination rate, where the efficiency of natural selection is high. The age and the length of these 300 TEs are consistent with relatively young and long insertions reaching high frequencies due to the action of positive selection. Besides, we identified a set of 21 fixed TEs also likely to be adaptive. Indeed, we, and others, found evidence of selection for 84 of these reference TE insertions. The analysis of the genes located nearby these 84 candidate adaptive insertions suggested that the functional response to selection is related with the GO categories of response to stimulus, behavior, and development. We further showed that a subset of the candidate adaptive TEs affects expression of nearby genes, and five of them have already been linked to an ecologically relevant phenotypic effect. Our results provide a more complete understanding of the genetic variation and the fitness-related traits relevant for adaptive evolution. Similar studies should help uncover the importance of TE-induced adaptive mutations in other species as well.

Target enrichment sequencing of 307 germplasm accessions identified ancestry of ancient and modern hybrids and signatures of adaptation and selection in sugarcane (Saccharum spp.), a 'sweet' crop with 'bitter' genomes

Sugarcane (Saccharum spp.) is a highly energy-efficient crop primarily for sugar and bio-ethanol production. Sugarcane genetics and cultivar improvement have been extremely challenging largely due to its complex genomes with high polyploidy levels. In this study, we deeply sequenced the coding regions of 307 sugarcane germplasm accessions. Nearly five million sequence variations were catalogued. The average of 98× sequence depth enabled different allele dosages of sequence variation to be differentiated in this polyploid collection. With selected high-quality genome-wide SNPs, we performed population genomic studies and environmental association analysis. Results illustrated that the ancient sugarcane hybrids, S. barberi and S. sinense, and modern sugarcane hybrids are significantly different in terms of genomic compositions, hybridization processes and their potential ancestry contributors. Linkage disequilibrium (LD) analysis showed a large extent of LD in sugarcane, with 962.4 Kbp, 2739.2 Kbp and 3573.6 Kbp for S. spontaneum, S. officinarum and modern S. hybrids respectively. Candidate selective sweep regions and genes were identified during domestication and historical selection processes of sugarcane in addition to genes associated with environmental variables at the original locations of the collection. This research provided an extensive amount of genomic resources for sugarcane community and the in-depth population genomic analyses shed light on the breeding and evolution history of sugarcane, a highly polyploid species.

The genomic basis of adaptation to calcareous and siliceous soils in Arabidopsis lyrata

Edaphic conditions are important determinants of plant fitness. While much has been learnt in recent years about plant adaptation to heavy metal contaminated soils, the genomic basis underlying adaptation to calcareous and siliceous substrates remains largely unknown. We performed a reciprocal germination experiment and whole-genome resequencing in natural calcareous and siliceous populations of diploid Arabidopsis lyrata to test for edaphic adaptation and detect signatures of selection at loci associated with soil-mediated divergence. In parallel, genome scans on respective diploid ecotypes from the Arabidopsis arenosa species complex were undertaken, to search for shared patterns of adaptive genetic divergence. Soil ecotypes of A. lyrata display significant genotype-by-treatment responses for seed germination. Sequence (SNPs) and copy-number variants (CNVs) point towards loci involved in ion transport as the main targets of adaptive genetic divergence. Two genes exhibiting high differentiation among soil types in A. lyrata further share trans-specific single nucleotide polymorphisms with A. arenosa. This work applies experimental and genomic approaches to study edaphic adaptation in A. lyrata and suggests that physiological response to elemental toxicity and deficiency underlies the evolution of calcareous and siliceous ecotypes. The discovery of shared adaptive variation between sister species indicates that ancient polymorphisms contribute to adaptive evolution.

Drought adaptation in Arabidopsis thaliana by extensive genetic loss-of-function

Interdisciplinary syntheses are needed to scale up discovery of the environmental drivers and molecular basis of adaptation in nature. Here we integrated novel approaches using whole genome sequences, satellite remote sensing, and transgenic experiments to study natural loss-of-function alleles associated with drought histories in wild Arabidopsis thaliana. The genes we identified exhibit population genetic signatures of parallel molecular evolution, selection for loss-of-function, and shared associations with flowering time phenotypes in directions consistent with longstanding adaptive hypotheses seven times more often than expected by chance. We then confirmed predicted phenotypes experimentally in transgenic knockout lines. These findings reveal the importance of drought timing to explain the evolution of alternative drought tolerance strategies and further challenge popular assumptions about the adaptive value of genetic loss-of-function in nature. These results also motivate improved species-wide sequencing efforts to better identify loss-of-function variants and inspire new opportunities for engineering climate resilience in crops.

Water shortages caused by droughts lead to crop losses that affect billions of people around the world each year. By discovering how wild plants adapt to drought, it may be possible to identify traits and genes that help to improve the growth of crop plants when water is scarce. It has been suggested that plants have adapted to droughts by flowering at times of the year when droughts are less likely to occur. For example, if droughts are more likely to happen in spring, the plants may delay flowering until the summer.

Arabidopsis thaliana is a small plant that is found across Eurasia, Africa and North America, including in areas that are prone to drought at different times of the year. Individual plants of the same species may carry different versions of the same gene (known as alleles). Some of these alleles may not work properly and are referred to as loss-of-function alleles. Monroe et al. investigated whether A. thaliana plants carry any loss-of-function alleles that are associated with droughts happening in the spring or summer, and whether they are linked to when those plants will flower.

Monroe et al. analyzed satellite images collected over the last 30 years to measure when droughts have occurred. Next, they searched genome sequences of Arabidopsis thaliana for alleles that might help the plants to adapt to droughts in the spring or summer. Combining the two approaches revealed that loss-of-function alleles associated with spring droughts were strongly predicted to be associated with the plants flowering later in the year. Similarly, loss-of-function alleles associated with summer droughts were predicted to be associated with the plants flowering earlier in the year.

These findings support the idea that plants can adapt to drought by changing when they produce flowers, and suggest that loss-of-function alleles play a major role in this process. New techniques for editing genes mean it is easier than ever to generate new loss-of-function alleles in specific genes. Therefore, the results presented by Monroe et al. may help researchers to develop new varieties of crop plants that are better adapted to droughts.

Genome-wide patterns of local adaptation in Western European Drosophila melanogaster natural populations

Signatures of spatially varying selection have been investigated both at the genomic and transcriptomic level in several organisms. In Drosophila melanogaster, the majority of these studies have analyzed North American and Australian populations, leading to the identification of several loci and traits under selection. However, several studies based mainly in North American populations showed evidence of admixture that likely contributed to the observed population differentiation patterns. Thus, disentangling demography from selection might be challenging when analyzing these populations. European populations could help identify loci under spatially varying selection provided that no recent admixture from African populations would have occurred. In this work, we individually sequence the genome of 42 European strains collected in populations from contrasting environments: Stockholm (Sweden) and Castellana Grotte (Southern Italy). We found low levels of population structure and no evidence of recent African admixture in these two populations. We thus look for patterns of spatially varying selection affecting individual genes and gene sets. Besides single nucleotide polymorphisms, we also investigated the role of transposable elements in local adaptation. We concluded that European populations are a good dataset to identify candidate loci under spatially varying selection. The analysis of the two populations sequenced in this work in the context of all the available D. melanogaster data allowed us to pinpoint genes and biological processes likely to be relevant for local adaptation. Identifying and analyzing populations with low levels of population structure and admixture should help to disentangle selective from non-selective forces underlying patterns of population differentiation in other species as well.

Genome-wide scans of selection highlight the impact of biotic and abiotic constraints in natural populations of the model grass Brachypodium distachyon

Grasses are essential plants for ecosystem functioning. Quantifying the selective pressures that act on natural variation in grass species is therefore essential regarding biodiversity maintenance. In this study, we investigate the selection pressures that act on two distinct populations of the grass model Brachypodium distachyon without prior knowledge about the traits under selection. We took advantage of whole-genome sequencing data produced for 44 natural accessions of B. distachyon and used complementary genome-wide selection scans (GWSS) methods to detect genomic regions under balancing and positive selection. We show that selection is shaping genetic diversity at multiple temporal and spatial scales in this species, and affects different genomic regions across the two populations. Gene ontology annotation of candidate genes reveals that pathogens may constitute important factors of positive and balancing selection in B. distachyon. We eventually cross-validated our results with quantitative trait locus data available for leaf-rust resistance in this species and demonstrate that, when paired with classical trait mapping, GWSS can help pinpointing candidate genes for further molecular validation. Thanks to a near base-perfect reference genome and the large collection of freely available natural accessions collected across its natural range, B. distachyon appears as a prime system for studies in ecology, population genomics and evolutionary biology.

Measuring the dynamic photosynthome

Background

Photosynthesis underpins plant productivity and yet is notoriously sensitive to small changes in environmental conditions, meaning that quantitation in nature across different time scales is not straightforward. The 'dynamic' changes in photosynthesis (i.e. the kinetics of the various reactions of photosynthesis in response to environmental shifts) are now known to be important in driving crop yield.

Scope

It is known that photosynthesis does not respond in a timely manner, and even a small temporal 'mismatch' between a change in the environment and the appropriate response of photosynthesis toward optimality can result in a fall in productivity. Yet the most commonly measured parameters are still made at steady state or a temporary steady state (including those for crop breeding purposes), meaning that new photosynthetic traits remain undiscovered.

Conclusions

There is a great need to understand photosynthesis dynamics from a mechanistic and biological viewpoint especially when applied to the field of 'phenomics' which typically uses large genetically diverse populations of plants. Despite huge advances in measurement technology in recent years, it is still unclear whether we possess the capability of capturing and describing the physiologically relevant dynamic features of field photosynthesis in sufficient detail. Such traits are highly complex, hence we dub this the 'photosynthome'. This review sets out the state of play and describes some approaches that could be made to address this challenge with reference to the relevant biological processes involved.