In plants, like in humans, DNA is arranged into sections known as genes that are in turn organised into structures called chromosomes. Mutations that modify the activity of these genes can help plant species to adapt to a new environment or to extreme conditions such as drought. However, successful adaptation often requires changes in many different genes. If these sets of genes are located close to each other on the same chromosome, any mutations will likely be passed onto the next generation together. If the genes are located further away, or even on different chromosomes, they may instead be inherited separately so that the next generation does not benefit as much from the adaptation.

A chromosome inversion – when a segment of chromosome breaks off and reattaches the other way around – can increase the likelihood that sets of mutations on the same chromosome will be inherited together. Many previous studies have found that chromosome inversions tend to drive the ability of species to adapt to different environments by keeping together mutations that affect the same characteristics. However, it is not clear how inversions affect the repeatability of the adaptation, that is, if another group of closely related plants faced the same challenge in their environment would they evolve in the same way, or would they evolve a new response?

To address this question, Soudi, Jahani et al. used a genetics approach known as a genome wide association study to explore how three closely related species of sunflower have adapted to their respective environments. Two of the species grow in various environments across the centre and west of the USA that are often hot and dry, whereas the third species is restricted to the more humid coastal plain of Texas, USA.

The experiments found that a few key genes had changed in all three sunflower species. However, each species also had mutations in a larger set of unique genes that were not changed in the other species. Regions of chromosomes harbouring inversions in one of the species tended to have more of the key genes within them, compared to other genomic regions. This was also true for species that did not have inversions in those regions. This demonstrates that genes in regions affected by chromosome inversions can still help plants adapt to changes in the environment even in the absence of inversions.

Sunflowers are widely grown for their edible oily seeds. In the future, some of the key genes identified in this work may be useful candidates for plant breeding to improve the resilience of sunflowers to drought, high temperatures and other environmental challenges.

The genetic basis of local adaptation is sometimes highly repeatable, with examples of large effect genes driving responses in multiple species, such as FT affecting flowering time in numerous plants (Izawa, 2007; Auge et al., 2019) or Mc1r driving colour polymorphisms in vertebrates (Manceau et al., 2010; Rosenblum et al., 2014). Local adaptation can also be repeated for polygenic traits, with significant patterns of similar association found across many loci for comparisons of conifers (Yeaman et al., 2016), maize and its wild relative teosinte (Tittes et al., 2021; Wang et al., 2021), and Brassicaceae (Bohutínská et al., 2021), to name a few. While we have a growing number of examples of repeatability in the basis of adaptation, it is also interesting to know if species use different genes to adapt to the same selection pressure. Genotypic redundancy – the potential for many genotypes to yield a given phenotype – is one critical factor affecting the repeatability of adaptation, as high redundancy would be expected to result in lower repeatability (Yeaman et al., 2018). Another critical factor affecting repeatability is shared standing variation, whether present due to introgression or incomplete lineage sorting. In either case, variants shared among lineages are much more likely to contribute to a repeated response than new mutations (MacPherson and Nuismer 2018; Ralph and Coop 2015). Consistent with this, Bohutínská et al., 2021 found that repeatability was negatively related to phylogenetic distance. While we are now accumulating more studies about the repeatability of adaptation, we still have very few examples and much remains unknown about the relative importance of these factors (Yeaman, 2022).

Inversions have been implicated in local adaptation in many species (Wellenreuther and Bernatchez, 2018), likely due to their effect to suppress recombination among inverted and non-inverted haplotypes, and thereby maintain linkage disequilibrium (LD) among beneficial combinations of locally adapted alleles (Rieseberg, 2001; Noor et al. 2001; Kirkpatrick and Barton, 2006). This has been approached by models studying the establishment of inversions that capture combinations of locally adapted alleles present as standing variation (e.g. Kirkpatrick and Barton, 2006), as well as models examining the accumulation of locally adapted mutations within inversions (e.g. Schaal et al., 2022). If there is variation in the density of loci that can potentially contribute to local adaptation, inversions would be expected to preferentially establish and be retained in regions harbouring a high density of such loci (and this expectation would hold for both the capture and accumulation models). We would also expect to see stronger signatures of repeated local adaptation in such high-density regions. Despite mounting evidence of their importance in adaptation, it is unclear how inversions may covary with repeatability of adaptation among species. A fundamental parameter of importance in these models is the relationship between migration rate and strength of selection on individual alleles, which may not make persistent contributions to local adaptation without the suppressing effects of recombination if selection is too weak (Yeaman and Whitlock, 2011; Bürger and Akerman, 2011). If most alleles have small effects relative to migration rate and can only contribute to local adaptation via the benefit of the recombination-suppressing effect of an inversion, then we would expect little repeatability at the site of an inversion – other species lacking the inversion would not tend to use that same region for adaptation because selection would be too weak for alleles to persist. On the other hand, if some loci are particularly important for local adaptation and regularly yield mutations of large effect, with these patterns being conserved among species, repeatability within regions harbouring inversions may be substantial. Thus, studying whether adaptation at the same genomic region harbouring an inversion is observed in other species lacking the inversion can give insights about the underlying architecture of adaptation, and the evolution and maintenance of inversions.

Here, we explore the repeatability of local adaptation in three species of sunflowers, Helianthus annuus, Helianthus argophyllus, and Helianthus petiolaris (Figure 1), which harbour large regions of suppressed recombination (‘haploblocks’), most of which are inversions, and are often associated with adaptive traits (Todesco et al., 2020). H. annuus, the common sunflower, is the closest wild relative of cultivated sunflower, which was domesticated from it around 4000 years ago (Blackman et al., 2011). Populations of H. annuus are distributed throughout the central and western United States and generally found on mesic soils, but can grow in a variety of disturbed or extreme habitats, such as semi-desertic or frequently flooded areas, as well as salt marshes. H. petiolaris, the prairie sunflower, prefers sandier soils, and ecotypes of this species are adapted to sand sheets and sand dunes (Ostevik et al., 2016). Here we include samples from two subspecies: H. petiolaris ssp. petiolaris, which is commonly found in the southern Great Plains, and H. petiolaris ssp. fallax, which is limited to more arid regions in Colorado, Utah, New Mexico, and Arizona (Heiser et al., 1969). H. petiolaris and H. annuus have broad and overlapping distributions throughout the central and western United States and appear to have adapted to similar changes in temperature, moisture, and photoperiod regimes. There is also evidence indicating H. annuus and H. petiolaris have likely been exchanging genes during much of their history of divergence (Strasburg and Rieseberg, 2008), although partially isolated by strong pre- and post-zygotic barriers (Sambatti et al., 2012). The third species, H. argophyllus, the silverleaf sunflower, is found exclusively in the southeast coast of Texas and includes both an early flowering ecotype on the coastal barrier islands and a late flowering ecotype inland (Moyers and Rieseberg, 2016; Todesco et al., 2020). H. argophyllus is thought to have undergone cycles of sympatry and allopatry with H. annuus and H. petiolaris over time, but currently only overlaps with H. annuus, with thus overlap likely being a recent event (Heiser, 1951).

Figure 1

Download asset Open asset

Using broad sampling across the ranges of these species (Figure 1) and sequencing data first published by Todesco et al., 2020, we study the genomic basis of local adaptation by conducting genome-wide scans for associations with climatic and soil environmental variables and phenotypes measured in a common garden. An inherent problem in studying the basis of local adaptation is accounting for the covariance between genome-wide population genetic structure and environment. When the selection pressure driving local adaptation tends to parallel a main axis of demographic expansion or isolation by distance, neutral alleles will have spatial patterns resembling those of causal alleles. Methods that do not correct for population structure have large numbers of false positives because allele frequencies at neutral loci tend to correlate with the environment more than expected by chance (Lotterhos and Whitlock, 2014). By contrast, methods that use structure correction tend to suffer from false negatives because the causal loci have similar patterns of allele frequency variation as the genomic background, and so their statistical signatures are decreased (DeRaad et al., 2021; Booker et al., 2023a). When a non-corrected approach is applied to multiple species, the false-positive problem can be mitigated when testing for loci that contribute to repeated adaptation as the same gene should not tend to be a false positive in multiple species more often than expected by chance (Yeaman et al., 2016). Here we use a conventional genome-wide association study (GWAS) approach with structure correction to study the basis of phenotypes, which vary both within and among populations, and use an uncorrected approach to test association to environment, which only varies among populations (as all individuals within the same population have the same value of environmental variable). For the environmental association portions of the study, we also explore the effect of structure correction in some analyses to contrast the false-negative vs. false-positive problem inherent in each approach. We apply these approaches individually to H. annuus, H. argophyllus, and the two H. petiolaris subspecies, and make pairwise and higher-order comparisons among combinations of these four lineages to study repeatability.

Our overall aim is to characterize the extent of repeatability of local adaptation for different traits and environments and explore the importance of inversions and other factors that may drive differences in repeatability. This analysis builds upon the work of Todesco et al., 2020, which focused on the evolution of haploblocks within each species, but did not systematically study patterns of repeatability in signatures of association among species. We begin by using the strength of correlations between environment and phenotypes measured in a common garden to identify which environmental variables are likely driving local adaptation within each species, and use an index to quantify the relative similarity in these patterns among pairs of species (see Appendix 1—figure 1 for a schematic overview of methods and research questions). We then conduct genome-wide association tests to identify regions most strongly associated with phenotypes and environments within each species. To assess the similarity in these association statistics among species, we use methods similar to Yeaman et al., 2016 to identify regions of the genome that exhibit greater similarity in signatures of association among pairs of lineages than expected by chance, which are likely driven by repeated adaptation using the same genes (but not necessarily the same alleles). We then explore how the number and total size of these regions of repeated association covary with the index of similarity in environmental importance. We also explore whether these regions of repeated association tend to covary with patterns of shared standing variation and test if they are enriched within previously identified haploblocks to explore the role of inversions in adaptation. We then identify candidate genes that show particularly strong signatures of repeated association across multiple lineages. Finally, based on the number of signatures of association that are shared vs. non-shared among lineages, we estimate the number of regions genome-wide that could potentially contribute to local adaptation for each variable and phenotype, which is related to the genotypic redundancy. Taken together, our results show that local adaptation in sunflowers tends to involve both strong repeatability at a small number of genes, often associated with inversions, coupled with high redundancy and non-repeated responses across a much larger number of loci.

Despite its critical importance in shaping the architecture of adaptation, little is known about the extent of genotypic redundancy underlying different traits (Barghi et al., 2020; Láruson et al., 2020; Yeaman, 2022). Here, we have shown evidence of significant repeatability in the basis of local adaptation (Figures 4 and 5), but also an abundance of species-specific, non-repeated signatures (Figure 3; Appendix 1—figure 14). In particular, we find that regions of the genome that harbour inversions in one species also tend to be strongly enriched for signatures of association in other species lacking the inversion (Figure 5). Taken together, this suggests that local adaptation in these species is highly flexible – different species apparently use quite different sets of loci to adapt to the same environment – yet still involves some component that has minimal redundancy, with inversions playing a particularly important role. Some of the ‘usual suspects’ show up in the set of significantly repeated loci identified by PicMin: in addition to the homologs of the FT (FLOWERING LOCUS T) gene reported by Todesco et al., 2020, we also found hits for several other genes involved in circadian regulation and flowering time, including PRR3 (Para et al., 2007), TOE1 (Aukerman and Sakai, 2003), and PHYC (Takano et al., 2005; Chen et al., 2014), which is known to regulate photoperiodic responses in Arabidopsis accessions (Balasubramanian et al., 2006). Other top hits included genes involved in plant development and auxin transport (PIN3 [Keuskamp et al., 2010], ARF4 [Pekker et al., 2005], and MN [Bhatia et al., 2016]), and plant immunity (CRK13 [Acharya et al., 2007] and PRR2 [Cheval et al., 2017]).

The increased repeatability found in regions of the genome that harbour inversions in only one species is particularly interesting. Inversions are commonly associated with local adaptation (Wellenreuther and Bernatchez, 2018), likely because they reduce the rate at which recombination breaks up combinations of co-selected alleles (Kirkpatrick and Barton, 2006), which perhaps facilitates contributions by alleles that would be individually too weakly selected to overcome swamping by migration (Bürger and Akerman, 2011; Yeaman and Whitlock, 2011; Schaal et al., 2022). Evidence here suggests these regions still tend to contribute to adaptation in the species lacking the recombination-suppressing effect of an inversion, consistent with a strong effect of selection relative to migration on at least one locus in the region (rather than adaptation exclusively via many alleles of small effect; Yeaman, 2022). As an example, chromosome 15 harbours a large (72 Mbp) haploblock in H. annuus that is strongly associated with NFFD, and also shows some signatures of association in the other taxa, with particularly strong signatures of association on either side of the haploblock in H. argophyllus (Figure 3A). Interestingly, loci of repeated association identified by PicMin within this region include two genes whose homologs are known to regulate responses to cold: COLD-RESPONSIVE PROTEIN KINASE 1,CRPK1 (Liu et al., 2017) and LATE ELONGATED HYPOCOTYL, LHY (Mizoguchi et al., 2002; Dong et al., 2011). It seems likely that strong selection relative to migration is therefore acting upon several loci in this region, and in many others harbouring inversions.

The observed repeatability associated with inversions further supports the local adaptation model as an explanation for the long-term persistence of segregating inversions (at least in sunflower), rather than mechanisms based on dominance or meiotic drive (Rieseberg, 2001). If there is variation across the genome in the density of loci with the potential to be involved in local adaptation, then the establishment and maintenance of inversions would be biased towards regions harbouring a high density of such loci under this model. If the genomic basis for local adaptation is conserved amongst species, then these same regions are more likely to have high repeatability. Thus, our observation of genomic regions harbouring inversions also being enriched for WRAs is consistent with this general model for inversion evolution. Unfortunately, our observations do not provide much insight into whether inversions evolve through the capture (e.g. Kirkpatrick and Barton, 2006) or accumulation (e.g. Schaal et al., 2022) type of model as either model would be consistent with our results. Most of the sunflower inversions are >1 My old, and therefore predate any current local adaptation patterns, but likely do not predate the genes underlying local adaptation (which appear to be shared among the species we studied). As for the alleles underlying local adaptation, they may be younger than the inversions, but as our work suggests, these regions are prone to harbouring locally adaptive alleles, so it is possible that they also harboured other ancestral locally adaptive alleles. While many studies have demonstrated the importance of inversions for adaptation (Wellenreuther and Bernatchez, 2018; Hager et al., 2022), to our knowledge only two other studies have documented the involvement of the same loci making contributions in the absence of the recombination-suppressing effect of the inversion (Lee et al., 2017; Coughlan and Willis, 2019). This also highlights how comparative studies of a species lacking an inversion may help identify which genes are driving adaptation in another species with an inversion, as segregating inversions tend to have extensive LD that prevents identification of any potential targets of natural selection within them.

While our results suggest a large number of loci can potentially contribute to adaptation, implying high redundancy (Figure 3B; Appendix 1—figure 14), there are several factors that complicate inference. Separating the effects of drift and selection to detect signatures of local adaptation is notoriously difficult because population structure often covaries with features of the environment that drive adaptation (Lotterhos and Whitlock, 2014; Hoban et al., 2016; DeRaad et al., 2021). As found by other analyses (Yeaman et al., 2016; DeRaad et al., 2021), when we use structure correction in our genotype–environment association tests, we find many fewer signatures of repeated association (Appendix 1—figure 7 vs. Appendix 1—figure 11), likely due to the reduced power of BayPass to detect true positives when the environment covaries with population structure (Booker et al., 2023a). Here, we have side-stepped the issues involved in correction for population structure by instead relying on comparisons among species to identify loci with associations more extreme than expected by chance in multiple species, using the null-W and PicMin tests. This assumes that most loci in the genome are not involved in local adaptation, so that the relatively small proportion that are driving adaptation can therefore be picked out due to their tendency to fall into the tail of the distribution of association statistics in at least one other species (which would only occur at the rate of random chance under a purely neutral model). While this approach should be relatively robust for identifying loci with repeated patterns of adaptation, there is no way to formally estimate significance of associations found in only a single species, many of which may be false positives due to covariation of population structure and environment. However, even if we assume that 80% of the observed non-repeated loci are false positives, we still find that estimates of the effective number of loci contributing to adaptation (L_eff) are in the hundreds (Appendix 1—figure 15C). As we are unable to detect loci of small effect due to the limited power, but these loci still likely make important contributions to heritability, our estimates of L_eff will also be biased downwards by excluding these potentially important drivers. Finally, it is also difficult to exclude the contribution of introgression or incomplete lineage sorting to the observed signatures of repeated association. If a locus tended to be highly introgressed between two species in a restricted region of their range, it is possible it could also covary with environment and therefore result in a signature of repeated association that would be mistakenly interpreted as adaptation. While it is difficult to preclude this from our analysis, regions of the genome with signatures of repeated association do not tend to have higher levels of shared standing variation than background regions (Appendix 1—figure 6), suggesting that this is not a broad explanation for our observations. If anything, shared standing variation would be expected to increase repeatability of adaptation (MacPherson and Nuismer, 2017), adding further weight to the inference that there is high genotypic redundancy in these species.

In general, we find that temperature is the strongest driver of repeated adaptation at both the phenotypic and genomic levels. We quantified local adaptation at the phenotypic level using correlations between traits measured in common gardens and the home environment of the population they were sampled from. Across all phenotypes, these correlations tended to be strongest and most similar among species for temperature variables, particularly in comparisons among H. annuus and H. petiolaris subspecies (Appendix 1—figure 3). We see similar patterns of repeatability reflected in the genome, where temperature variables also tend to have the greatest repeatability (Figure 4). The similarity in these phenotypic and genomic signatures is consistent with an effect of strong selection as other artefactual or drift-based explanations for repeatability would not be expected to reflect patterns found at the phenotypic level. It should be noted that the reduced importance of soil variables in the SIPEC index might be partly driven by the fact that all traits were measured on above-ground features due to the difficulty of getting non-disruptive phenotypic measurements for roots.

Taken together, these results suggest that some fraction of the genome contributes to adaptation with low redundancy and high repeatability (which tends to be enriched within genomic regions where there exist inversions in at least one species), while the remainder of the adaptive response is driven by loci with high redundancy and species-specific contributions. Theoretical models of adaptive evolution necessarily involve simplified representations of genetics: population genetic approaches explore cases where strong selection deterministically drives a change in allele frequency at one or a few loci (usually without epistasis), whereas quantitative genetic approaches make the infinitesimal assumption that phenotypic change can be realized through small frequency changes at many loci (Barton et al., 2017; Barghi et al., 2020; Yeaman, 2022). Such models implicitly make quite different assumptions about redundancy: if a population genetic model of directional selection includes no epistasis, then it is implicitly assuming there is no redundancy as every locus can increase fitness and no locus can affect the potential of another locus to increase fitness. On the other hand, quantitative genetic models implicitly assume complete redundancy as the loci driving trait variation have interchangeable effects. Given that differences in redundancy result in quite different evolutionary dynamics (Hollinger et al., 2019; Yeaman, 2022), it is important to re-evaluate the behaviour of theoretical models in light of what we can learn about redundancy from empirical data. Our results suggest that adaptation is a complex process that does not map cleanly onto the assumptions of either approach. At least some component of trait variation experiences sufficiently strong selection and has low enough redundancy to drive the significant signatures of repeatability above the random null that we observe here. On the other hand, most of the strongest signatures of association are specific to a single species, implying high redundancy for most of the genome, even if there are still some detectable signatures of repeatability. Such observations may be consistent with classifying loci into ‘core’ and ‘peripheral’ sets similar to the omnigenic model, but would involve classification based on redundancy rather than gene expression, as advocated by Boyle et al. (2017). Theoretical models adopting a two-class representation of the genetic basis of trait variation may provide more realism than models assuming a simpler distribution, and perhaps yield new insights about qualitatively different dynamical behaviour. Given the difficulty of rigorously quantifying the contribution of alleles of small effect and loci contributing to local adaptation in only one species, it remains an open question what proportion of adaptive response is governed by loci with dynamical behaviour that follow the infinitesimal assumption. This could perhaps be quantified explicitly by decomposing heritabilities for traits driving local adaptation, as has been done for standing variation in humans (Visscher et al., 2017), but this would require a considerable increase in sample size as methods such as GCTA (Yang et al., 2011) do not yield accurate estimates at the sample sizes used here. Future work could perhaps resolve this question using multi-generational crossing designs to break up the LD that tends to accompany local adaptation, allowing a more accurate parsing of the effect sizes of genomic regions and the contribution by repeated vs. redundant components of adaptive trait variation.

All scripts used for analysis and figures are deposited in a Dryad Digital Repository. Raw sequence data are deposited under NCBI BioProject accessions PRJNA532579, PRJNA398560, and PRJNA564337, as described in Todesco et al., 2020.

1. 1. Sam Y
   2. Shaghayegh S
   3. Mojtaba J

   (2023) Dryad Digital Repository

   Data from: Repeatability of adaptation in sunflowers reveals that genomic regions harbouring inversions also drive adaptation in species lacking an inversion.

   https://doi.org/10.5061/dryad.wpzgmsbtd

1. 1. Todesco M

   (2020) NCBI BioProject

   ID PRJNA532579. Wild Helianthus GWAS and GEA.

   https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA532579
2. 1. Todesco M

   (2020) NCBI BioProject

   ID PRJNA398560. Wild and Weedy Helianthus annuus whole genome resequencing.

   https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA398560
3. 1. Todesco M

   (2020) NCBI BioProject

   ID PRJNA564337. Hi-C sequencing of wild sunflowers.

   https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA564337

1. 1. Acharya BR
   2. Raina S
   3. Maqbool SB
   4. Jagadeeswaran G
   5. Mosher SL
   6. Appel HM
   7. Schultz JC
   8. Klessig DF
   9. Raina R

   (2007) Overexpression of CRK13, an Arabidopsis cysteine-rich receptor-like kinase, results in enhanced resistance to Pseudomonas syringae

   The Plant Journal 50:488–499.

   https://doi.org/10.1111/j.1365-313X.2007.03064.x

   • PubMed
   • Google Scholar
2. Software

   1. Alexa A
   2. Rahnenfuhrer J

   (2023) topGO: enrichment analysis for gene Ontology, version 2.54.0

   R Package.

   https://bioconductor.org/packages/topGO
3. 1. Auge GA
   2. Penfield S
   3. Donohue K

   (2019) Pleiotropy in developmental regulation by flowering-pathway genes: is it an evolutionary constraint?

   The New Phytologist 224:55–70.

   https://doi.org/10.1111/nph.15901

   • PubMed
   • Google Scholar
4. 1. Aukerman MJ
   2. Sakai H

   (2003) Regulation of flowering time and floral organ identity by a MicroRNA and its APETALA2-like target genes

   The Plant Cell 15:2730–2741.

   https://doi.org/10.1105/tpc.016238

   • PubMed
   • Google Scholar
5. 1. Badouin H
   2. Gouzy J
   3. Grassa CJ
   4. Murat F
   5. Staton SE
   6. Cottret L
   7. Lelandais-Brière C
   8. Owens GL
   9. Carrère S
   10. Mayjonade B
   11. Legrand L
   12. Gill N
   13. Kane NC
   14. Bowers JE
   15. Hubner S
   16. Bellec A
   17. Bérard A
   18. Bergès H
   19. Blanchet N
   20. Boniface MC
   21. Brunel D
   22. Catrice O
   23. Chaidir N
   24. Claudel C
   25. Donnadieu C
   26. Faraut T
   27. Fievet G
   28. Helmstetter N
   29. King M
   30. Knapp SJ
   31. Lai Z
   32. Le Paslier MC
   33. Lippi Y
   34. Lorenzon L
   35. Mandel JR
   36. Marage G
   37. Marchand G
   38. Marquand E
   39. Bret-Mestries E
   40. Morien E
   41. Nambeesan S
   42. Nguyen T
   43. Pegot-Espagnet P
   44. Pouilly N
   45. Raftis F
   46. Sallet E
   47. Schiex T
   48. Thomas J
   49. Vandecasteele C
   50. Varès D
   51. Vear F
   52. Vautrin S
   53. Crespi M
   54. Mangin B
   55. Burke JM
   56. Salse J
   57. Muños S
   58. Vincourt P
   59. Rieseberg LH
   60. Langlade NB

   (2017) The sunflower genome provides insights into oil metabolism, flowering and Asterid evolution

   Nature 546:148–152.

   https://doi.org/10.1038/nature22380

   • PubMed
   • Google Scholar
6. 1. Balasubramanian S
   2. Sureshkumar S
   3. Agrawal M
   4. Michael TP
   5. Wessinger C
   6. Maloof JN
   7. Clark R
   8. Warthmann N
   9. Chory J
   10. Weigel D

   (2006) The PHYTOCHROME C photoreceptor gene mediates natural variation in flowering and growth responses of Arabidopsis thaliana

   Nature Genetics 38:711–715.

   https://doi.org/10.1038/ng1818

   • PubMed
   • Google Scholar
7. 1. Barghi N
   2. Hermisson J
   3. Schlötterer C

   (2020) Polygenic adaptation: a unifying framework to understand positive selection

   Nature Reviews. Genetics 21:769–781.

   https://doi.org/10.1038/s41576-020-0250-z

   • PubMed
   • Google Scholar
8. 1. Barton NH
   2. Etheridge AM
   3. Véber A

   (2017) The infinitesimal model: Definition, derivation, and implications

   Theoretical Population Biology 118:50–73.

   https://doi.org/10.1016/j.tpb.2017.06.001

   • PubMed
   • Google Scholar
9. 1. Belton J-M
   2. McCord RP
   3. Gibcus JH
   4. Naumova N
   5. Zhan Y
   6. Dekker J

   (2012) Hi-C: a comprehensive technique to capture the conformation of genomes

   Methods 58:268–276.

   https://doi.org/10.1016/j.ymeth.2012.05.001

   • PubMed
   • Google Scholar
10. 1. Bhatia N
    2. Bozorg B
    3. Larsson A
    4. Ohno C
    5. Jönsson H
    6. Heisler MG

    (2016) Auxin Acts through MONOPTEROS to regulate plant cell polarity and pattern phyllotaxis

    Current Biology 26:3202–3208.

    https://doi.org/10.1016/j.cub.2016.09.044

    • PubMed
    • Google Scholar
11. 1. Blackman BK
    2. Scascitelli M
    3. Kane NC
    4. Luton HH
    5. Rasmussen DA
    6. Bye RA
    7. Lentz DL
    8. Rieseberg LH

    (2011) Sunflower domestication alleles support single domestication center in eastern North America

    PNAS 108:14360–14365.

    https://doi.org/10.1073/pnas.1104853108

    • PubMed
    • Google Scholar
12. 1. Bohutínská M
    2. Vlček J
    3. Yair S
    4. Laenen B
    5. Konečná V
    6. Fracassetti M
    7. Slotte T
    8. Kolář F

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

    PNAS 118:e2022713118.

    https://doi.org/10.1073/pnas.2022713118

    • PubMed
    • Google Scholar
13. 1. Bolger AM
    2. Lohse M
    3. Usadel B

    (2014) Trimmomatic: a flexible trimmer for Illumina sequence data

    Bioinformatics 30:2114–2120.

    https://doi.org/10.1093/bioinformatics/btu170

    • PubMed
    • Google Scholar
14. 1. Booker TR
    2. Yeaman S
    3. Whitlock MC

    (2020) Variation in recombination rate affects detection of outliers in genome scans under neutrality

    Molecular Ecology 29:4274–4279.

    https://doi.org/10.1111/mec.15501

    • PubMed
    • Google Scholar
15. 1. Booker TR
    2. Yeaman S
    3. Whiting JR
    4. Whitlock MC

    (2023a) The WZA: A window-based method for characterizing genotype-environment associations

    Molecular Ecology Resources 1:13768.

    https://doi.org/10.1111/1755-0998.13768

    • PubMed
    • Google Scholar
16. 1. Booker TR
    2. Yeaman S
    3. Whitlock MC

    (2023b) Using genome scans to identify genes used repeatedly for adaptation

    Evolution; International Journal of Organic Evolution 77:801–811.

    https://doi.org/10.1093/evolut/qpac063

    • PubMed
    • Google Scholar
17. 1. Bürger R
    2. Akerman A

    (2011) The effects of linkage and gene flow on local adaptation: A two-locus continent-island model

    Theoretical Population Biology 80:272–288.

    https://doi.org/10.1016/j.tpb.2011.07.002

    • PubMed
    • Google Scholar
18. 1. Chen A
    2. Li C
    3. Hu W
    4. Lau MY
    5. Lin H
    6. Rockwell NC
    7. Martin SS
    8. Jernstedt JA
    9. Lagarias JC
    10. Dubcovsky J

    (2014) PHYTOCHROME C plays a major role in the acceleration of wheat flowering under long-day photoperiod

    PNAS 111:10037–10044.

    https://doi.org/10.1073/pnas.1409795111

    • PubMed
    • Google Scholar
19. 1. Cheval C
    2. Perez M
    3. Leba LJ
    4. Ranty B
    5. Perochon A
    6. Reichelt M
    7. Mithöfer A
    8. Robe E
    9. Mazars C
    10. Galaud JP
    11. Aldon D

    (2017) PRR2, a pseudo-response regulator, promotes salicylic acid and camalexin accumulation during plant immunity

    Scientific Reports 7:6979.

    https://doi.org/10.1038/s41598-017-07535-8

    • PubMed
    • Google Scholar
20. 1. Coughlan JM
    2. Willis JH

    (2019) Dissecting the role of a large chromosomal inversion in life history divergence throughout the Mimulus guttatus species complex

    Molecular Ecology 28:1343–1357.

    https://doi.org/10.1111/mec.14804

    • PubMed
    • Google Scholar
21. 1. DeRaad DA
    2. Cobos ME
    3. Alkishe A
    4. Ashraf U
    5. Ahadji-Dabla KM
    6. Nuñez-Penichet C
    7. Peterson AT

    (2021) Genome-environment association methods comparison supports omnigenic adaptation to ecological niche in malaria vector mosquitoes

    Molecular Ecology 30:6468–6485.

    https://doi.org/10.1111/mec.16094

    • PubMed
    • Google Scholar
22. 1. Dong MA
    2. Farré EM
    3. Thomashow MF

    (2011) Circadian clock-associated 1 and late elongated hypocotyl regulate expression of the c-repeat binding factor (CBF) pathway in Arabidopsis

    PNAS 108:7241–7246.

    https://doi.org/10.1073/pnas.1103741108

    • PubMed
    • Google Scholar
23. 1. Gautier M

    (2015) Genome-wide scan for adaptive divergence and association with population-specific covariates

    Genetics 201:1555–1579.

    https://doi.org/10.1534/genetics.115.181453

    • PubMed
    • Google Scholar
24. 1. Grimm DG
    2. Roqueiro D
    3. Salomé PA
    4. Kleeberger S
    5. Greshake B
    6. Zhu W
    7. Liu C
    8. Lippert C
    9. Stegle O
    10. Schölkopf B
    11. Weigel D
    12. Borgwardt KM

    (2017) easyGWAS: a cloud-based platform for comparing the results of genome-wide association studies

    The Plant Cell 29:5–19.

    https://doi.org/10.1105/tpc.16.00551

    • PubMed
    • Google Scholar
25. 1. Gunther T
    2. Coop G

    (2013) Robust identification of local adaptation from allele frequencies

    Genetics 195:205–220.

    https://doi.org/10.1534/genetics.113.152462

    • PubMed
    • Google Scholar
26. 1. Hager ER
    2. Harringmeyer OS
    3. Wooldridge TB
    4. Theingi S
    5. Gable JT
    6. McFadden S
    7. Neugeboren B
    8. Turner KM
    9. Jensen JD
    10. Hoekstra HE

    (2022) A chromosomal inversion contributes to divergence in multiple traits between deer mouse ecotypes

    Science 377:399–405.

    https://doi.org/10.1126/science.abg0718

    • PubMed
    • Google Scholar
27. 1. Heiser CB

    (1951) Hybridization in the annual sunflowers: helianthus annuus x H. argophyllus

    The American Naturalist 85:65–72.

    https://doi.org/10.1086/281651

    • Google Scholar
28. 1. Heiser C
    2. Smith D
    3. Clevenger S
    4. Martin W

    (1969)

    The North American sunflowers (Helianthus)

    Memoirs of the Torrey Botanical Club 22:1–218.

    • Google Scholar
29. 1. Hoban S
    2. Kelley JL
    3. Lotterhos KE
    4. Antolin MF
    5. Bradburd G
    6. Lowry DB
    7. Poss ML
    8. Reed LK
    9. Storfer A
    10. Whitlock MC

    (2016) Finding the genomic basis of local adaptation: pitfalls, practical solutions, and future directions

    The American Naturalist 188:379–397.

    https://doi.org/10.1086/688018

    • PubMed
    • Google Scholar
30. 1. Hollinger I
    2. Pennings PS
    3. Hermisson J

    (2019) Polygenic adaptation: From sweeps to subtle frequency shifts

    PLOS Genetics 15:e1008035.

    https://doi.org/10.1371/journal.pgen.1008035

    • PubMed
    • Google Scholar
31. 1. Izawa T

    (2007) Adaptation of flowering-time by natural and artificial selection in Arabidopsis and rice

    Journal of Experimental Botany 58:3091–3097.

    https://doi.org/10.1093/jxb/erm159

    • PubMed
    • Google Scholar
32. Book

    1. Jeffreys H

    (1961)

    Theory of Probability

    Oxford University Press.

    • Google Scholar
33. 1. Kang HM
    2. Sul JH
    3. Service SK
    4. Zaitlen NA
    5. Kong S-Y
    6. Freimer NB
    7. Sabatti C
    8. Eskin E

    (2010) Variance component model to account for sample structure in genome-wide association studies

    Nature Genetics 42:348–354.

    https://doi.org/10.1038/ng.548

    • PubMed
    • Google Scholar
34. 1. Keuskamp DH
    2. Pollmann S
    3. Voesenek L
    4. Peeters AJM
    5. Pierik R

    (2010) Auxin transport through PIN-FORMED 3 (PIN3) controls shade avoidance and fitness during competition

    PNAS 107:22740–22744.

    https://doi.org/10.1073/pnas.1013457108

    • PubMed
    • Google Scholar
35. 1. Kirkpatrick M
    2. Barton N

    (2006) Chromosome inversions, local adaptation and speciation

    Genetics 173:419–434.

    https://doi.org/10.1534/genetics.105.047985

    • PubMed
    • Google Scholar
36. 1. Láruson ÁJ
    2. Yeaman S
    3. Lotterhos KE

    (2020) The importance of genetic redundancy in evolution

    Trends in Ecology & Evolution 35:809–822.

    https://doi.org/10.1016/j.tree.2020.04.009

    • Google Scholar
37. 1. Lee CR
    2. Wang B
    3. Mojica JP
    4. Mandáková T
    5. Prasad K
    6. Goicoechea JL
    7. Perera N
    8. Hellsten U
    9. Hundley HN
    10. Johnson J
    11. Grimwood J
    12. Barry K
    13. Fairclough S
    14. Jenkins JW
    15. Yu Y
    16. Kudrna D
    17. Zhang J
    18. Talag J
    19. Golser W
    20. Ghattas K
    21. Schranz ME
    22. Wing R
    23. Lysak MA
    24. Schmutz J
    25. Rokhsar DS
    26. Mitchell-Olds T

    (2017) Young inversion with multiple linked QTLs under selection in a hybrid zone

    Nature Ecology & Evolution 1:119.

    https://doi.org/10.1038/s41559-017-0119

    • PubMed
    • Google Scholar
38. Preprint

    1. Li H

    (2013) Aligning Sequence Reads, Clone Sequences and Assembly Contigs with BWA-MEM

    arXiv.

    https://doi.org/10.48550/arXiv.1303.3997

    • Google Scholar
39. 1. Liu Z
    2. Jia Y
    3. Ding Y
    4. Shi Y
    5. Li Z
    6. Guo Y
    7. Gong Z
    8. Yang S

    (2017) Plasma Membrane CRPK1-Mediated Phosphorylation of 14-3-3 proteins induces their nuclear import to fine-tune cbf signaling during cold response

    Molecular Cell 66:117–128.

    https://doi.org/10.1016/j.molcel.2017.02.016

    • PubMed
    • Google Scholar
40. 1. Lotterhos KE
    2. Whitlock MC

    (2014) Evaluation of demographic history and neutral parameterization on the performance of FST outlier tests

    Molecular Ecology 23:2178–2192.

    https://doi.org/10.1111/mec.12725

    • PubMed
    • Google Scholar
41. 1. MacPherson A
    2. Nuismer SL

    (2017) The probability of parallel genetic evolution from standing genetic variation

    Journal of Evolutionary Biology 30:326–337.

    https://doi.org/10.1111/jeb.13006

    • PubMed
    • Google Scholar
42. 1. Manceau M
    2. Domingues VS
    3. Linnen CR
    4. Rosenblum EB
    5. Hoekstra HE

    (2010) Convergence in pigmentation at multiple levels: mutations, genes and function

    Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 365:2439–2450.

    https://doi.org/10.1098/rstb.2010.0104

    • PubMed
    • Google Scholar
43. 1. Marie-Nelly H
    2. Marbouty M
    3. Cournac A
    4. Flot JF
    5. Liti G
    6. Parodi DP
    7. Syan S
    8. Guillén N
    9. Margeot A
    10. Zimmer C
    11. Koszul R

    (2014) High-quality genome (re)assembly using chromosomal contact data

    Nature Communications 5:5695.

    https://doi.org/10.1038/ncomms6695

    • PubMed
    • Google Scholar
44. 1. Mizoguchi T
    2. Wheatley K
    3. Hanzawa Y
    4. Wright L
    5. Mizoguchi M
    6. Song HR
    7. Carré IA
    8. Coupland G

    (2002) LHY and CCA1 are partially redundant genes required to maintain circadian rhythms in Arabidopsis

    Developmental Cell 2:629–641.

    https://doi.org/10.1016/s1534-5807(02)00170-3

    • PubMed
    • Google Scholar
45. 1. Moyers BT
    2. Rieseberg LH

    (2016) Remarkable life history polymorphism may be evolving under divergent selection in the silverleaf sunflower

    Molecular Ecology 25:3817–3830.

    https://doi.org/10.1111/mec.13723

    • PubMed
    • Google Scholar
46. 1. Ostevik KL
    2. Andrew RL
    3. Otto SP
    4. Rieseberg LH

    (2016) Multiple reproductive barriers separate recently diverged sunflower ecotypes: reproductive barriers between sunflower ecotypes

    Evolution; International Journal of Organic Evolution 70:2322–2335.

    https://doi.org/10.1111/evo.13027

    • PubMed
    • Google Scholar
47. 1. Para A
    2. Farré EM
    3. Imaizumi T
    4. Pruneda-Paz JL
    5. Harmon FG
    6. Kay SA

    (2007) PRR3 Is a vascular regulator of TOC1 stability in the Arabidopsis circadian clock

    The Plant Cell 19:3462–3473.

    https://doi.org/10.1105/tpc.107.054775

    • PubMed
    • Google Scholar
48. 1. Pekker I
    2. Alvarez JP
    3. Eshed Y

    (2005) Auxin response factors mediate Arabidopsis organ asymmetry via modulation of KANADI activity

    The Plant Cell 17:2899–2910.

    https://doi.org/10.1105/tpc.105.034876

    • PubMed
    • Google Scholar
49. 1. Rieseberg LH

    (2001) Chromosomal rearrangements and speciation

    Trends in Ecology & Evolution 16:351–358.

    https://doi.org/10.1016/s0169-5347(01)02187-5

    • PubMed
    • Google Scholar
50. 1. Rosenblum EB
    2. Parent CE
    3. Brandt EE

    (2014) The molecular basis of phenotypic convergence

    Annual Review of Ecology, Evolution, and Systematics 45:203–226.

    https://doi.org/10.1146/annurev-ecolsys-120213-091851

    • Google Scholar
51. 1. Sambatti JBM
    2. Strasburg JL
    3. Ortiz-Barrientos D
    4. Baack EJ
    5. Rieseberg LH

    (2012) Reconciling extremely strong barriers with high levels of gene exchange in annual sunflowers

    Evolution; International Journal of Organic Evolution 66:1459–1473.

    https://doi.org/10.1111/j.1558-5646.2011.01537.x

    • PubMed
    • Google Scholar
52. 1. Schaal SM
    2. Haller BC
    3. Lotterhos KE

    (2022) Inversion invasions: when the genetic basis of local adaptation is concentrated within inversions in the face of gene flow

    Philosophical Transactions of the Royal Society of London. Series B, Biological Sciences 377:20210200.

    https://doi.org/10.1098/rstb.2021.0200

    • PubMed
    • Google Scholar
53. 1. Sedlazeck FJ
    2. Rescheneder P
    3. von Haeseler A

    (2013) NextGenMap: fast and accurate read mapping in highly polymorphic genomes

    Bioinformatics 29:2790–2791.

    https://doi.org/10.1093/bioinformatics/btt468

    • PubMed
    • Google Scholar
54. 1. Strasburg JL
    2. Rieseberg LH

    (2008) Molecular demographic history of the annual sunflowers Helianthus annuus and H. petiolaris--large effective population sizes and rates of long-term gene flow

    Evolution; International Journal of Organic Evolution 62:1936–1950.

    https://doi.org/10.1111/j.1558-5646.2008.00415.x

    • PubMed
    • Google Scholar
55. 1. Takano M
    2. Inagaki N
    3. Xie X
    4. Yuzurihara N
    5. Hihara F
    6. Ishizuka T
    7. Yano M
    8. Nishimura M
    9. Miyao A
    10. Hirochika H
    11. Shinomura T

    (2005) Distinct and cooperative functions of phytochromes A, B, and C in the control of deetiolation and flowering in rice

    The Plant Cell 17:3311–3325.

    https://doi.org/10.1105/tpc.105.035899

    • PubMed
    • Google Scholar
56. Preprint

    1. Tittes S
    2. Lorant A
    3. McGinty S
    4. Holland JB
    5. Sánchez-González J
    6. Seetharam A
    7. Tenaillon M
    8. Ross-Ibarra J

    (2021) Not so local: the population genetics of convergent adaptation in maize and teosinte

    bioRxiv.

    https://doi.org/10.1101/2021.09.09.459637

    • Google Scholar
57. 1. Todesco M
    2. Owens GL
    3. Bercovich N
    4. Légaré J-S
    5. Soudi S
    6. Burge DO
    7. Huang K
    8. Ostevik KL
    9. Drummond EBM
    10. Imerovski I
    11. Lande K
    12. Pascual-Robles MA
    13. Nanavati M
    14. Jahani M
    15. Cheung W
    16. Staton SE
    17. Muños S
    18. Nielsen R
    19. Donovan LA
    20. Burke JM
    21. Yeaman S
    22. Rieseberg LH

    (2020) Massive haplotypes underlie ecotypic differentiation in sunflowers

    Nature 584:602–607.

    https://doi.org/10.1038/s41586-020-2467-6

    • PubMed
    • Google Scholar
58. Book

    1. Van der Auwera GA
    2. O’Connor BD

    (2000)

    Genomics in the Cloud

    O’Reilly Media, Inc.

    • Google Scholar
59. 1. Visscher PM
    2. Wray NR
    3. Zhang Q
    4. Sklar P
    5. McCarthy MI
    6. Brown MA
    7. Yang J

    (2017) 10 Years of GWAS discovery: biology, function, and translation

    American Journal of Human Genetics 101:5–22.

    https://doi.org/10.1016/j.ajhg.2017.06.005

    • PubMed
    • Google Scholar
60. 1. Wang T
    2. Hamann A
    3. Spittlehouse DL
    4. Murdock TQ

    (2012) ClimateWNA—high-resolution spatial climate data for western north america

    Journal of Applied Meteorology and Climatology 51:16–29.

    https://doi.org/10.1175/JAMC-D-11-043.1

    • Google Scholar
61. 1. Wang L
    2. Josephs EB
    3. Lee KM
    4. Roberts LM
    5. Rellán-Álvarez R
    6. Ross-Ibarra J
    7. Hufford MB

    (2021) Molecular parallelism underlies convergent highland adaptation of maize landraces

    Molecular Biology and Evolution 38:3567–3580.

    https://doi.org/10.1093/molbev/msab119

    • PubMed
    • Google Scholar
62. 1. Wellenreuther M
    2. Bernatchez L

    (2018) Eco-evolutionary genomics of chromosomal inversions

    Trends in Ecology & Evolution 33:427–440.

    https://doi.org/10.1016/j.tree.2018.04.002

    • Google Scholar
63. Book

    1. Whitlock M
    2. Schluter D

    (2020)

    The Analysis of Biological Data

    Macmillan.

    • Google Scholar
64. 1. Yang J
    2. Lee SH
    3. Goddard ME
    4. Visscher PM

    (2011) GCTA: A tool for genome-wide complex trait analysis

    American Journal of Human Genetics 88:76–82.

    https://doi.org/10.1016/j.ajhg.2010.11.011

    • PubMed
    • Google Scholar
65. 1. Yeaman S
    2. Whitlock MC

    (2011) The genetic architecture of adaptation under migration-selection balance

    Evolution; International Journal of Organic Evolution 65:1897–1911.

    https://doi.org/10.1111/j.1558-5646.2011.01269.x

    • PubMed
    • Google Scholar
66. 1. Yeaman S
    2. Hodgins KA
    3. Lotterhos KE
    4. Suren H
    5. Nadeau S
    6. Degner JC
    7. Nurkowski KA
    8. Smets P
    9. Wang T
    10. Gray LK
    11. Liepe KJ
    12. Hamann A
    13. Holliday JA
    14. Whitlock MC
    15. Rieseberg LH
    16. Aitken SN

    (2016) Convergent local adaptation to climate in distantly related conifers

    Science 353:1431–1433.

    https://doi.org/10.1126/science.aaf7812

    • PubMed
    • Google Scholar
67. 1. Yeaman S
    2. Gerstein AC
    3. Hodgins KA
    4. Whitlock MC

    (2018) Quantifying how constraints limit the diversity of viable routes to adaptation

    PLOS Genetics 14:e1007717.

    https://doi.org/10.1371/journal.pgen.1007717

    • PubMed
    • Google Scholar
68. 1. Yeaman S

    (2022) Evolution of polygenic traits under global vs local adaptation

    Genetics 220:iyab134.

    https://doi.org/10.1093/genetics/iyab134

    • PubMed
    • Google Scholar

Author details

1. Shaghayegh Soudi

   Department of Biological Sciences, University of Calgary, Calgary, Canada

   Contribution

   Resources, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review and editing

   Contributed equally with

   Mojtaba Jahani

   Competing interests

   No competing interests declared

2. Mojtaba Jahani

   1. Department of Biological Sciences, University of Calgary, Calgary, Canada
   2. Department of Botany, University of British Columbia, Vancouver, Canada

   Contribution

   Resources, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review and editing

   Contributed equally with

   Shaghayegh Soudi

   Competing interests

   No competing interests declared

3. Marco Todesco

   1. Department of Botany, University of British Columbia, Vancouver, Canada
   2. Michael Smith Laboratories, University of British Columbia, Vancouver, Canada
   3. Irving K. Barber Faculty of Science, University of British Columbia Okanagan, Kelowna, Canada

   Contribution

   Visualization, Formal analysis, Investigation, Methodology, Writing – review and editing

   Competing interests

   No competing interests declared

4. Gregory L Owens

   Department of Biology, University of Victoria, Victoria, Canada

   Contribution

   Visualization, Formal analysis, Investigation, Methodology, Writing – review and editing

   Competing interests

   No competing interests declared

5. Natalia Bercovich

   Department of Botany, University of British Columbia, Vancouver, Canada

   Contribution

   Resources, Visualization, Methodology

   Competing interests

   No competing interests declared

6. Loren H Rieseberg

   Department of Botany, University of British Columbia, Vancouver, Canada

   Contribution

   Conceptualization, Supervision, Investigation, Project administration, Writing – review and editing

   Competing interests

   No competing interests declared

7. Sam Yeaman

   Department of Biological Sciences, University of Calgary, Calgary, Canada

   Contribution

   Conceptualization, Resources, Formal analysis, Supervision, Investigation, Visualization, Methodology, Writing – original draft, Project administration, Writing – review and editing

   For correspondence

   samuel.yeaman@ucalgary.ca

   Competing interests

   No competing interests declared

   "This ORCID iD identifies the author of this article:" 0000-0002-1706-8699

• 1,305

  views

• 112

  downloads

• 3

  citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.