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

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

Relatively weak inbreeding depression in selfing but also in outcrossing populations of North American Arabidopsis lyrata

Hermaphroditic plants can potentially self‐fertilize, but most possess adaptations that promote outcrossing. However, evolutionary transitions to higher selfing rates are frequent. Selfing comes with a transmission advantage over outcrossing, but self‐progeny may suffer from inbreeding depression, which forms the main barrier to the evolution of higher selfing rates. Here, we assessed inbreeding depression in the North American herb Arabidopsis lyrata, which is normally self‐incompatible, with a low frequency of self‐compatible plants. However, a few populations have become fixed for self‐compatibility and have high selfing rates. Under greenhouse conditions, we estimated mean inbreeding depression per seed (based on cumulative vegetative performance calculated as the product of germination, survival and aboveground biomass) to be 0.34 for six outcrossing populations, and 0.26 for five selfing populations. Exposing plants to drought and inducing defences with jasmonic acid did not magnify these estimates. For outcrossing populations, however, inbreeding depression per seed may underestimate true levels of inbreeding depression, because self‐incompatible plants showed strong reductions in seed set after (enforced) selfing. Inbreeding‐depression estimates incorporating seed set averaged 0.63 for outcrossing populations (compared to 0.30 for selfing populations). However, this is likely an overestimate because exposing plants to 5% CO₂ to circumvent self‐incompatibility to produce selfed seed might leave residual effects of self‐incompatibility that contribute to reduced seed set. Nevertheless, our estimates of inbreeding depression were clearly lower than previous estimates based on the same performance traits in outcrossing European populations of A. lyrata, which may help explain why selfing could evolve in North American A. lyrata.     

Parallel adaptation to climate above the 35th parallel

Independent or parallel evolution of similar traits is key to understanding the genetics and limitations of adaptation. Adaptation from the same genetic changes in different populations defines parallel evolution. Such genetic changes can derive from standing ancestral variation or de novo mutations and excludes instances of adaptive introgression. In this issue of Molecular Ecology, Walden et al.(2020) investigate the scale of parallel climate adaptation from standing genetic variation between two North American Arabidopsis lyrata lineages, each formed by a distinct evolutionary history during the last glacial cycle. By identifying adaptive variants correlated with three ecologically significant climatic gradients, they show that instead of the same genetic variants or even genes, parallel evolution is only observed at the level of biological processes. The evolution of independent adaptive variants to climate in two genetically close lineages is explained by their different post‐glacial demographic histories. Separate glacial refugia and strong population bottlenecks were probably sufficient to change the landscape of shared allele frequencies, hindering the possibility of parallel evolution.     

Longitudinal trends in climate drive flowering time clines in North American Arabidopsis thaliana

Introduced species frequently show geographic differentiation, and when differentiation mirrors the ancestral range, it is often taken as evidence of adaptive evolution. The mouse-ear cress (Arabidopsis thaliana) was introduced to North America from Eurasia 150-200 years ago, providing an opportunity to study parallel adaptation in a genetic model organism. Here, we test for clinal variation in flowering time using 199 North American (NA) accessions of A. thaliana, and evaluate the contributions of major flowering time genes FRI, FLC, and PHYC as well as potential ecological mechanisms underlying differentiation. We find evidence for substantial within population genetic variation in quantitative traits and flowering time, and putatively adaptive longitudinal differentiation, despite low levels of variation at FRI, FLC, and PHYC and genome-wide reductions in population structure relative to Eurasian (EA) samples. The observed longitudinal cline in flowering time in North America is parallel to an EA cline, robust to the effects of population structure, and associated with geographic variation in winter precipitation and temperature. We detected major effects of FRI on quantitative traits associated with reproductive fitness, although the haplotype associated with higher fitness remains rare in North America. Collectively, our results suggest the evolution of parallel flowering time clines through novel genetic mechanisms.     

A genome‐wide search for local adaptation in a terrestrial‐breeding frog reveals vulnerability to climate change

Terrestrial‐breeding amphibians are likely to be vulnerable to warming and drying climates, as their embryos require consistent moisture for successful development. Adaptation to environmental change will depend on sufficient genetic variation existing within or between connected populations. Here, we use Single Nucleotide Polymorphism (SNP) data to investigate genome‐wide patterns in genetic diversity, gene flow and local adaptation in a terrestrial‐breeding frog (Pseudophryne guentheri) subject to a rapidly drying climate and recent habitat fragmentation. The species was sampled across 12 central and range‐edge populations (192 samples), and strong genetic structure was apparent, as were high inbreeding coefficients. Populations showed differences in genetic diversity, and one population lost significant genetic diversity in a decade. More than 500 SNP loci were putatively under directional selection, and 413 of these loci were correlated with environmental variables such as temperature, rainfall, evaporation and soil moisture. One locus showed homology to a gene involved in the activation of maturation in Xenopus oocytes, which may facilitate rapid development of embryos in drier climates. The low genetic diversity, strong population structuring and presence of local adaptation revealed in this study shows why management strategies such as targeted gene flow may be necessary to assist isolated populations to adapt to future climates.     

Genome‐wide signatures of flowering adaptation to climate temperature: Regional analyses in a highly diverse native range of Arabidopsis thaliana

Current global change is fueling an interest to understand the genetic and molecular mechanisms of plant adaptation to climate. In particular, altered flowering time is a common strategy for escape from unfavourable climate temperature. In order to determine the genomic bases underlying flowering time adaptation to this climatic factor, we have systematically analysed a collection of 174 highly diverse Arabidopsis thaliana accessions from the Iberian Peninsula. Analyses of 1.88 million single nucleotide polymorphisms provide evidence for a spatially heterogeneous contribution of demographic and adaptive processes to geographic patterns of genetic variation. Mountains appear to be allele dispersal barriers, whereas the relationship between flowering time and temperature depended on the precise temperature range. Environmental genome‐wide associations supported an overall genome adaptation to temperature, with 9.4% of the genes showing significant associations. Furthermore, phenotypic genome‐wide associations provided a catalogue of candidate genes underlying flowering time variation. Finally, comparison of environmental and phenotypic genome‐wide associations identified known (Twin Sister of FT, FRIGIDA‐like 1, and Casein Kinase II Beta chain 1) and new (Epithiospecifer Modifier 1 and Voltage‐Dependent Anion Channel 5) genes as candidates for adaptation to climate temperature by altered flowering time. Thus, this regional collection provides an excellent resource to address the spatial complexity of climate adaptation in annual plants.     

Evidence of genomic adaptation to climate in Eucalyptus microcarpa: Implications for adaptive potential to projected climate change

Understanding whether populations can adapt in situ or whether interventions are required is of key importance for biodiversity management under climate change. Landscape genomics is becoming an increasingly important and powerful tool for rapid assessments of climate adaptation, especially in long‐lived species such as trees. We investigated climate adaptation in Eucalyptus microcarpa using the DArTseq genomic approach. A combination of F_ST outlier and environmental association analyses were performed using >4200 genomewide single nucleotide polymorphisms (SNPs) from 26 populations spanning climate gradients in southeastern Australia. Eighty‐one SNPs were identified as putatively adaptive, based on significance in F_ST outlier tests and significant associations with one or more climate variables related to temperature (70/81), aridity (37/81) or precipitation (35/81). Adaptive SNPs were located on all 11 chromosomes, with no particular region associated with individual climate variables. Climate adaptation appeared to be characterized by subtle shifts in allele frequencies, with no consistent fixed differences identified. Based on these associations, we predict adaptation under projected changes in climate will include a suite of shifts in allele frequencies. Whether this can occur sufficiently rapidly through natural selection within populations, or would benefit from assisted gene migration, requires further evaluation. In some populations, the absence or predicted increases to near fixation of particular adaptive alleles hint at potential limits to adaptive capacity. Together, these results reinforce the importance of standing genetic variation at the geographic level for maintaining species’ evolutionary potential.     

Signatures of local adaptation along environmental gradients in a range‐expanding damselfly (Ischnura elegans)

Insect distributions are shifting rapidly in response to climate change and are undergoing rapid evolutionary change. We investigate the molecular signatures underlying local adaptation in the range‐expanding damselfly, Ischnura elegans. Using a landscape genomic approach combined with generalized dissimilarity modelling (GDM), we detect selection signatures on loci via allelic frequency change along environmental gradients. We analyse 13,612 single nucleotide polymorphisms (SNPs), derived from restriction site‐associated DNA sequencing (RADseq), in 426 individuals from 25 sites spanning the I. elegans distribution in Sweden, including its expanding northern range edge. Environmental association analysis (EAA) and the magnitude of allele frequency change along the range expansion gradient revealed significant signatures of selection in relation to high maximum summer temperature, high mean annual precipitation and low wind speeds at the range edge. SNP annotations with significant signatures of selection revealed gene functions associated with ongoing range expansion, including heat shock proteins (HSP40 and HSP70), ion transport (V‐ATPase) and visual processes (long‐wavelength‐sensitive opsin), which have implications for thermal stress response, salinity tolerance and mate discrimination, respectively. We also identified environmental thresholds where climate‐mediated selection is likely to be strong, and indicate that I. elegans is rapidly adapting to the climatic environment during its ongoing range expansion. Our findings empirically validate an integrative approach for detecting spatially explicit signatures of local adaptation along environmental gradients.     

Imputation of 3 million SNPs in the Arabidopsis regional mapping population

Natural variation has become a prime resource to identify genetic variants that contribute to phenotypic variation. The regional mapping (RegMap) population is one of the most important populations for studying natural variation in Arabidopsis thaliana, and has been used in a large number of association studies and in studies on climatic adaptation. However, only 413 RegMap accessions have been completely sequenced, as part of the 1001 Genomes (1001G) Project, while the remaining 894 accessions have only been genotyped with the Affymetrix 250k chip. As a consequence, most association studies involving the RegMap are either restricted to the sequenced accessions, reducing power, or rely on a limited set of SNPs. Here we impute millions of SNPs to the 894 accessions that are exclusive to the RegMap, using the 1135 accessions of the 1001G Project as the reference panel. We assess imputation accuracy using a novel cross‐validation scheme, which we show provides a more reliable measure of accuracy than existing methods. After filtering out low accuracy SNPs, we obtain high‐quality genotypic information for 2029 accessions and 3 million markers. To illustrate the benefits of these imputed data, we reconducted genome‐wide association studies on five stress‐related traits and could identify novel candidate genes.     

Lineage‐specific sequence evolution and exon edge conservation partially explain the relationship between evolutionary rate and expression level in A. thaliana

Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the ‘edges’ of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage‐specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage‐specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage‐specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.     

Evolutionary shifts to self‐fertilisation restricted to geographic range margins in North American Arabidopsis lyrata

Cross‐fertilisation predominates in eukaryotes, but shifts to self‐fertilisation are common and ecologically and evolutionarily important. Reproductive assurance under outcross gamete limitation is one eco‐evolutionary process held responsible for the shift to selfing. Although small effective population size is a situation where selfing plants could theoretically benefit from reproductive assurance, empirical tests of the role of population size are rare. Here, we show that selfing evolved repeatedly at range margins, where historical demographic processes produced low effective population sizes. Outcrossing populations of North American Arabidopsis lyrata have low genetic diversity at geographic margins, with a signature of post‐glacial range expansion in the north and rear‐edge isolation in the south. Selfing populations occur at the margins of two genetic groups and never in their interior. These results corroborate small effective population size as the promoter of self‐fertilisation and have important implications for our understanding of species turnover, range limits and range dynamics.     

Identification of genetic variants associated with maize flowering time using an extremely large multi‐genetic background population

Flowering time is one of the major adaptive traits in domestication of maize and an important selection criterion in breeding. To detect more maize flowering time variants we evaluated flowering time traits using an extremely large multi‐ genetic background population that contained more than 8000 lines under multiple Sino‐United States environments. The population included two nested association mapping (NAM) panels and a natural association panel. Nearly 1 million single‐nucleotide polymorphisms (SNPs) were used in the analyses. Through the parallel linkage analysis of the two NAM panels, both common and unique flowering time regions were detected. Genome wide, a total of 90 flowering time regions were identified. One‐third of these regions were connected to traits associated with the environmental sensitivity of maize flowering time. The genome‐wide association study of the three panels identified nearly 1000 flowering time‐associated SNPs, mainly distributed around 220 candidate genes (within a distance of 1 Mb). Interestingly, two types of regions were significantly enriched for these associated SNPs – one was the candidate gene regions and the other was the approximately 5 kb regions away from the candidate genes. Moreover, the associated SNPs exhibited high accuracy for predicting flowering time.     

Association between putative functional variants in the PSMB9 gene and risk of melanoma – re‐analysis of published melanoma genome‐wide association studies

To mine possibly hidden causal single‐nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome‐wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer‐free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype–phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility.     

Lineage‐specific plasmid acquisition and the evolution of specialized pathogens in Bacillus thuringiensis and the Bacillus cereus group

Bacterial plasmids can vary from small selfish genetic elements to large autonomous replicons that constitute a significant proportion of total cellular DNA. By conferring novel function to the cell, plasmids may facilitate evolution but their mobility may be opposed by co‐evolutionary relationships with chromosomes or encouraged via the infectious sharing of genes encoding public goods. Here, we explore these hypotheses through large‐scale examination of the association between plasmids and chromosomal DNA in the phenotypically diverse Bacillus cereus group. This complex group is rich in plasmids, many of which encode essential virulence factors (Cry toxins) that are known public goods. We characterized population genomic structure, gene content and plasmid distribution to investigate the role of mobile elements in diversification. We analysed coding sequence within the core and accessory genome of 190 B. cereus group isolates, including 23 novel sequences and genes from 410 reference plasmid genomes. While cry genes were widely distributed, those with invertebrate toxicity were predominantly associated with one sequence cluster (clade 2) and phenotypically defined Bacillus thuringiensis. Cry toxin plasmids in clade 2 showed evidence of recent horizontal transfer and variable gene content, a pattern of plasmid segregation consistent with transfer during infectious cooperation. Nevertheless, comparison between clades suggests that co‐evolutionary interactions may drive association between plasmids and chromosomes and limit wider transfer of key virulence traits. Proliferation of successful plasmid and chromosome combinations is a feature of specialized pathogens with characteristic niches (Bacillus anthracis, B. thuringiensis) and has occurred multiple times in the B. cereus group.     

Genome‐wide scans detect adaptation to aridity in a widespread forest tree species

Patterns of adaptive variation within plant species are best studied through common garden experiments, but these are costly and time‐consuming, especially for trees that have long generation times. We explored whether genome‐wide scanning technology combined with outlier marker detection could be used to detect adaptation to climate and provide an alternative to common garden experiments. As a case study, we sampled nine provenances of the widespread forest tree species, Eucalyptus tricarpa, across an aridity gradient in southeastern Australia. Using a Bayesian analysis, we identified a suite of 94 putatively adaptive (outlying) sequence‐tagged markers across the genome. Population‐level allele frequencies of these outlier markers were strongly correlated with temperature and moisture availability at the site of origin, and with population differences in functional traits measured in two common gardens. Using the output from a canonical analysis of principal coordinates, we devised a metric that provides a holistic measure of genomic adaptation to aridity that could be used to guide assisted migration or genetic augmentation.     

Behavioral constraints on local adaptation and counter‐gradient variation: Implications for climate change

Resource allocation to growth, reproduction, and body maintenance varies within species along latitudinal gradients. Two hypotheses explaining this variation are local adaptation and counter‐gradient variation. The local adaptation hypothesis proposes that populations are adapted to local environmental conditions and are therefore less adapted to environmental conditions at other locations. The counter‐gradient variation hypothesis proposes that one population out performs others across an environmental gradient because its source location has greater selective pressure than other locations. Our study had two goals. First, we tested the local adaptation and counter‐gradient variation hypotheses by measuring effects of environmental temperature on phenotypic expression of reproductive traits in the burying beetle, Nicrophorus orbicollis Say, from three populations along a latitudinal gradient in a common garden experimental design. Second, we compared patterns of variation to evaluate whether traits covary or whether local adaptation of traits precludes adaptive responses by others. Across a latitudinal range, N. orbicollis exhibits variation in initiating reproduction and brood sizes. Consistent with local adaptation: (a) beetles were less likely to initiate breeding at extreme temperatures, especially when that temperature represents their source range; (b) once beetles initiate reproduction, source populations produce relatively larger broods at temperatures consistent with their local environment. Consistent with counter‐gradient variation, lower latitude populations were more successful at producing offspring at lower temperatures. We found no evidence for adaptive variation in other adult or offspring performance traits. This suite of traits does not appear to coevolve along the latitudinal gradient. Rather, response to selection to breed within a narrow temperature range may preclude selection on other traits. Our study highlights that N. orbicollis uses temperature as an environmental cue to determine whether to initiate reproduction, providing insight into how behavior is modified to avoid costly reproductive attempts. Furthermore, our results suggest a temperature constraint that shapes reproductive behavior.     

Metabolome‐genome‐wide association study dissects genetic architecture for generating natural variation in rice secondary metabolism

Plants produce structurally diverse secondary (specialized) metabolites to increase their fitness for survival under adverse environments. Several bioactive compounds for new drugs have been identified through screening of plant extracts. In this study, genome‐wide association studies (GWAS) were conducted to investigate the genetic architecture behind the natural variation of rice secondary metabolites. GWAS using the metabolome data of 175 rice accessions successfully identified 323 associations among 143 single nucleotide polymorphisms (SNPs) and 89 metabolites. The data analysis highlighted that levels of many metabolites are tightly associated with a small number of strong quantitative trait loci (QTLs). The tight association may be a mechanism generating strains with distinct metabolic composition through the crossing of two different strains. The results indicate that one plant species produces more diverse phytochemicals than previously expected, and plants still contain many useful compounds for human applications.     

Fitness declines towards range limits and local adaptation to climate affect dispersal evolution during climate‐induced range shifts

Dispersal ability will largely determine whether species track their climatic niches during climate change, a process especially important for populations at contracting (low‐latitude/low‐elevation) range limits that otherwise risk extinction. We investigate whether dispersal evolution at contracting range limits is facilitated by two processes that potentially enable edge populations to experience and adjust to the effects of climate deterioration before they cause extinction: (i) climate‐induced fitness declines towards range limits and (ii) local adaptation to a shifting climate gradient. We simulate a species distributed continuously along a temperature gradient using a spatially explicit, individual‐based model. We compare range‐wide dispersal evolution during climate stability vs. directional climate change, with uniform fitness vs. fitness that declines towards range limits (RLs), and for a single climate genotype vs. multiple genotypes locally adapted to temperature. During climate stability, dispersal decreased towards RLs when fitness was uniform, but increased when fitness declined towards RLs, due to highly dispersive genotypes maintaining sink populations at RLs, increased kin selection in smaller populations, and an emergent fitness asymmetry that favoured dispersal in low‐quality habitat. However, this initial dispersal advantage at low‐fitness RLs did not facilitate climate tracking, as it was outweighed by an increased probability of extinction. Locally adapted genotypes benefited from staying close to their climate optima; this selected against dispersal under stable climates but for increased dispersal throughout shifting ranges, compared to cases without local adaptation. Dispersal increased at expanding RLs in most scenarios, but only increased at the range centre and contracting RLs given local adaptation to climate.