Climate structures genetic variation across a species' elevation range: a test of range limits hypotheses

Local adaptation to an altitudinal gradient: the interplay between mean phenotypic trait variation and phenotypic plasticity in Mimulus laciniatus

Organisms can adapt to environmental heterogeneity through two mechanisms: (1) expression of population genetic variation or (2) phenotypic plasticity. In this study we investigated whether patterns of variation in both trait means and phenotypic plasticity along elevational and latitudinal clines in a North American endemic plant, Mimulus laciniatus, were consistent with local adaptation. We grew inbred lines of M. laciniatus from across the species' range in two common gardens varying in day length to measure mean and plastic trait expression in several traits previously shown to be involved in adaptation to M. laciniatus's rocky outcrop microhabitat: flowering time, size-related traits, and leaf shape. We examined correlations between the mean phenotype and phenotypic plasticity, and tested for a relationship between trait variation and population elevation and latitude. We did not find a strong correlation between mean and plastic trait expression at the individual genotype level suggesting that they operate under independent genetic controls. We identified multiple traits that show patterns consistent with local adaptation to elevation: critical photoperiod, flowering time, flower size, mean leaf lobing, and leaf lobing plasticity. These trends occur along multiple geographically independent altitudinal clines indicating that selection is a more likely cause of this pattern than gene flow among nearby populations with similar trait values. We also found that population variation in mean leaf lobing is associated with latitude. Our results indicate that both having more highly lobed leaves and greater leaf shape plasticity may be adaptive at high elevation within M. laciniatus. Our data strongly suggest that traits known to be under divergent selection between M. laciniatus and close relative Mimulus guttatus are also under locally varying selection within M. laciniatus.

Individual-based landscape genomics for conservation: An analysis pipeline

Landscape genomics can harness environmental and genetic data to inform conservation decisions by providing essential insights into how landscapes shape biodiversity. The massive increase in genetic data afforded by the genomic era provides exceptional resolution for answering critical conservation genetics questions. The accessibility of genomic data for non-model systems has also enabled a shift away from population-based sampling to individual-based sampling, which now provides accurate and robust estimates of genetic variation that can be used to examine the spatial structure of genomic diversity, population connectivity and the nature of environmental adaptation. Nevertheless, the adoption of individual-based sampling in conservation genetics has been slowed due, in large part, to concerns over how to apply methods developed for population-based sampling to individual-based sampling schemes. Here, we discuss the benefits of individual-based sampling for conservation and describe how landscape genomic methods, paired with individual-based sampling, can answer fundamental conservation questions. We have curated key landscape genomic methods into a user-friendly, open-source workflow, which we provide as a new R package, A Landscape Genomics Analysis Toolkit in R (algatr). The algatr package includes novel added functionality for all of the included methods and extensive vignettes designed with the primary goal of making landscape genomic approaches more accessible and explicitly applicable to conservation biology.

Genetic variability and the ecology of geographic range: A test of the central-marginal hypothesis in Australian scincid lizards

For many species, both local abundance and regional occupancy are highest near the centre of their geographic distributions. One hypothesis for this pattern is that niche suitability declines with increasing distance from a species geographic centre, such that populations near range margins are characterized by reduced density and increased patchiness. In these smaller edge populations, genetic drift is more powerful, leading to the loss of genetic diversity. This simple verbal model has been formalized as the central-marginal hypothesis, which predicts that core populations should have greater genetic diversity than edge populations. Here, we tested the central-marginal hypothesis using a genomic data set of 25 species-level taxa of Australian scincid lizards in the genera Ctenotus and Lerista. A majority of taxa in our data set showed range-wide patterns of genetic variation consistent with central-marginal hypothesis, and eight of 25 taxa showed significantly greater genetic diversity in the centre of their range. We then explored biological, historical, and methodological factors that might predict which taxa support the central-marginal hypothesis. We found that taxa with the strongest evidence for range expansion were the least likely to follow predictions of the central-marginal hypothesis. The majority of these taxa had range expansions that originated at the range edge, which led to a gradient of decreasing genetic diversity from the range edge to the core, contrary to the central-marginal hypothesis.

Demographic, physiological, and genetic factors linked to the poleward range expansion of the snail Nerita yoldii along the shoreline of China

Species range shift is one of the most significant consequences of climate change in the Anthropocene. A comprehensive study, including demographic, physiological, and genetic factors linked to poleward range expansion, is crucial for understanding how the expanding population occupies the new habitat. In the present study, we investigated the demographic, physiological, and genetic features of the intertidal gastropod Nerita yoldii, which has extended its northern limit by ~200 km over the former biogeographic break of the Yangtze River Estuary during recent decades. The neutral SNPs data showed that the new marginal populations formed a distinct cluster established by a few founders. Demographic modelling analysis revealed that the new marginal populations experienced a strong genetic bottleneck followed by recent demographic expansion. Successful expansion that overcame the founder effect might be attributed to its high capacity of rapid population growth and multiple introductions. According to the non-neutral SNPs under diversifying selection, there were high levels of heterozygosity in the new marginal populations, which might be beneficial for adapting to the novel thermal conditions. The common garden experiment showed that the new marginal populations have evolved divergent transcriptomic and physiological responses to heat stress, allowing them to occupy and survive in the novel environment. Lower transcriptional plasticity was observed in the new marginal populations. These results suggest a new biogeographic pattern of N. yoldii has formed with the occurrence of demographic, physiologic, and genetic changes, and emphasize the roles of adaptation of marginal populations during range expansion.

Hybrid zone of a tree in a Cerrado/Atlantic Forest ecotone as a hotspot of genetic diversity and conservation

The Cerrado, the largest Neotropical savanna, and the Brazilian Atlantic Forest form large ecotonal areas where savanna and forest habitats occupy adjacent patches with closely related species occurring side by side, providing opportunities for hybridization. Here, we investigated the evolutionary divergence between the savanna and forest ecotypes of the widely distributed tree Plathymenia reticulata (n = 233 individuals). Genetic structure analysis of P. reticulata was congruent with the recognition of two ecotypes, whose divergence captured the largest proportion of genetic variance in the data (F CT = 0.222 and F ST = 0.307). The ecotonal areas between the Cerrado and the Atlantic Forest constitute a hybrid zone in which a diversity of hybrid classes was observed, most of them corresponding to second-generation hybrids (F2) or backcrosses. Gene flow occurred mainly toward the forest ecotype. The genetic structure was congruent with isolation by environment, and environmental correlates of divergence were identified. The observed pattern of high genetic divergence between ecotypes may reflect an incipient speciation process in P. reticulata. The low genetic diversity of the P. reticulata forest ecotype indicate that it is threatened in areas with high habitat loss on Atlantic Forest. In addition, the high divergence from the savanna ecotype suggests it should be treated as a different unit of management. The high genetic diversity found in the ecotonal hybrid zone supports the view of ecotones as important areas for the origin and conservation of biodiversity in the Neotropics.

Rules of Plant Species Ranges: Applications for Conservation Strategies

Earth is changing rapidly and so are many plant species’ ranges. Here, we synthesize eco-evolutionary patterns found in plant range studies and how knowledge of species ranges can inform our understanding of species conservation in the face of global change. We discuss whether general biogeographic “rules” are reliable and how they can be used to develop adaptive conservation strategies of native plant species across their ranges. Rules considered include (1) factors that set species range limits and promote range shifts; (2) the impact of biotic interactions on species range limits; (3) patterns of abundance and adaptive properties across species ranges; (4) patterns of gene flow and their implications for genetic rescue, and (5) the relationship between range size and conservation risk. We conclude by summarizing and evaluating potential species range rules to inform future conservation and management decisions. We also outline areas of research to better understand the adaptive capacity of plants under environmental change and the properties that govern species ranges. We advise conservationists to extend their work to specifically consider peripheral and novel populations, with a particular emphasis on small ranges. Finally, we call for a global effort to identify, synthesize, and analyze prevailing patterns or rules in ecology to help speed conservation efforts.

A strategic sampling design revealed the local genetic structure of cold-water fluvial sculpin: a focus on groundwater-dependent water temperature heterogeneity

A key piece of information for ecosystem management is the relationship between the environment and population genetic structure. However, it is difficult to clearly quantify the effects of environmental factors on genetic differentiation because of spatial autocorrelation and analytical problems. In this study, we focused on stream ecosystems and the environmental heterogeneity caused by groundwater and constructed a sampling design in which geographic distance and environmental differences are not correlated. Using multiplexed ISSR genotyping by sequencing (MIG-seq) method, a fine-scale population genetics study was conducted in fluvial sculpin Cottus nozawae, for which summer water temperature is the determinant factor in distribution and survival. There was a clear genetic structure in the watershed. Although a significant isolation-by-distance pattern was detected in the watershed, there was no association between genetic differentiation and water temperature. Instead, asymmetric gene flow from relatively low-temperature streams to high-temperature streams was detected, indicating the importance of low-temperature streams and continuous habitats. The groundwater-focused sampling strategy yielded insightful results for conservation.

Dispersal evolution in temporally variable environments: implications for plant range dynamics

Dispersal-the movement of an individual from the site of birth to a different site for reproduction-is an ecological and evolutionary driver of species ranges that shapes patterns of colonization, connectivity, gene flow, and adaptation. In plants, the traits that influence dispersal often vary within and among species, are heritable, and evolve in response to the fitness consequences of moving through heterogeneous landscapes. Spatial and temporal variation in the quality and quantity of habitat are important sources of selection on dispersal strategies across species ranges. While recent reviews have evaluated the interactions between spatial variation in habitat and dispersal dynamics, the extent to which geographic variation in temporal variability can also shape range-wide patterns in dispersal traits has not been synthesized. In this paper, we summarize key predictions from metapopulation models that evaluate how dispersal evolves in response to spatial and temporal habitat variability. Next, we compile empirical data that quantify temporal variability in plant demography and patterns of dispersal trait variation across species ranges to evaluate the hypothesis that higher temporal variability favors increased dispersal at plant range limits. We found some suggestive evidence supporting this hypothesis while more generally identifying a major gap in empirical work evaluating plant metapopulation dynamics across species ranges and geographic variation in dispersal traits. To address this gap, we propose several future research directions that would advance our understanding of the interplay between spatiotemporal variability and dispersal trait variation in shaping the dynamics of current and future species ranges.

Long-term persistence of experimental populations beyond a species' natural range

Ecological experiments usually infer long-term processes from short-term data, and the analysis of geographic range limits is a good example. Species' geographic ranges may be limited by low fitness due to niche constraints, a hypothesis most directly tested by comparing the fitness of populations transplanted within and beyond the range. Such studies often fail to find beyond-range fitness declines strong enough to conclude that geographic range limits are solely imposed by niche limits. However, almost all studies only follow transplants for a single generation, which will underestimate the importance of niche limitation because critical but infrequent range-limiting events may be missed and methodological issues may artificially boost the fitness of beyond-range transplants. Here, we present the first multi-generation beyond-range transplant experiment that involves adequate replication and proper experimental controls. In 2005, experimental populations of the coastal dune plant Camissoniopsis cheiranthifolia were planted at four sites within and one site beyond the northern limit. Fitness of initial transplants was high beyond the limit, suggesting that the range was limited by dispersal and not niche constraints. To better address the niche-limitation hypothesis, we quantified density and fitness of descendant C. cheiranthifolia populations 12-14 yr (˜10 generations) after transplant. Average annual fruit production and density of reproductive individuals were as high beyond the range as at four comparable experimental populations and eight natural populations within the range, and the beyond-range population had more than tripled in size since it was planted. This provides unprecedented support for the conclusion that northern range limit of C. cheiranthifolia results from something other than niche limitation, likely involving constraints on local dispersal.

How is adaptive potential distributed within species ranges?

Quantitative genetic variation (QGV) represents a major component of adaptive potential and, if reduced toward range-edge populations, could prevent a species' expansion or adaptive response to rapid ecological change. It has been hypothesized that QGV will be lower at the range edge due to small populations-often the result of poor habitat quality-and potentially decreased gene flow. However, whether central populations are higher in QGV is unknown. We used a meta-analytic approach to test for a general QGV-range position relationship, including geographic and climatic distance from range centers. We identified 35 studies meeting our criteria, yielding nearly 1000 estimates of QGV (including broad-sense heritability, narrow-sense heritability, and evolvability) from 34 species. The relationship between QGV and distance from the geographic range or climatic niche center depended on the focal trait and how QGV was estimated. We found some evidence that QGV declines from geographic centers but that it increases toward niche edges; niche and geographic distances were uncorrelated. Nevertheless, few studies have compared QGV in both central and marginal regions or environments within the same species. We call for more research in this area and discuss potential research avenues related to adaptive potential in the context of global change.

Adaptation across geographic ranges is consistent with strong selection in marginal climates and legacies of range expansion

Every species experiences limits to its geographic distribution. Some evolutionary models predict that populations at range edges are less well adapted to their local environments due to drift, expansion load, or swamping gene flow from the range interior. Alternatively, populations near range edges might be uniquely adapted to marginal environments. In this study, we use a database of transplant studies that quantify performance at broad geographic scales to test how local adaptation, site quality, and population quality change from spatial and climatic range centers toward edges. We find that populations from poleward edges perform relatively poorly, both on average across all sites (15% lower population quality) and when compared to other populations at home (31% relative fitness disadvantage), consistent with these populations harboring high genetic load. Populations from equatorial edges also perform poorly on average (18% lower population quality) but, in contrast, outperform foreign populations (16% relative fitness advantage), suggesting that populations from equatorial edges have strongly adapted to unique environments. Finally, we find that populations from sites that are thermally extreme relative to the species' niche demonstrate strong local adaptation, regardless of their geographic position. Our findings indicate that both nonadaptive processes and adaptive evolution contribute to variation in adaptation across species' ranges.

What Do We Really Know About Adaptation at Range Edges?

Recent theory and empirical evidence have provided new insights regarding how evolutionary forces interact to shape adaptation at stable and transient range margins. Predictions regarding trait divergence at leading edges are frequently supported. However, declines in fitness at and beyond edges show that trait divergence has sometimes been insufficient to maintain high fitness, so identifying constraints to adaptation at range edges remains a key challenge. Indirect evidence suggests that range expansion may be limited by adaptive genetic variation, but direct estimates of genetic constraints at and beyond range edges are still scarce. Sequence data suggest increased genetic load in edge populations in several systems, but its causes and fitness consequences are usually poorly understood. The balance between maladaptive and positive effects of gene flow on fitness at range edges deserves further study. It is becoming increasingly clear that characterizations about degree of adaptation based solely on geographical peripherality are unsupported.

Populations Are Differentiated in Biological Rhythms without Explicit Elevational Clines in the Plant Mimulus laciniatus

Environmental variation along an elevational gradient can yield phenotypic differentiation resulting from varying selection pressures on plant traits related to seasonal responses. Thus, genetic clines can evolve in a suite of traits, including the circadian clock, that drives daily cycling in varied traits and that shares its genetic background with adaptation to seasonality. We used populations of annual Mimulus laciniatus from different elevations in the Sierra Nevada in California to explore among-population differentiation in the circadian clock, flowering responses to photoperiod, and phenological traits (days to cotyledon emergence, days to flowering, and days to seed ripening) in controlled common-garden conditions. Further, we examined correlations of these traits with environmental variables related to temperature and precipitation. We observed that the circadian period in leaf movement was differentiated among populations sampled within about 100 km, with population means varying by 1.6 h. Significant local genetic variation occurred within 2 populations in which circadian period among families varied by up to 1.8 h. Replicated treatments with variable ecologically relevant photoperiods revealed marked population differentiation in critical day length for flowering that ranged from 11.0 to 14.1 h, corresponding to the time period between late February and mid-May in the wild. Flowering time varied among populations in a 14-h photoperiod. Regardless of this substantial population-level diversity, obvious linear clinality in trait variability across elevations could not be determined based on our genotypic sample; it is possible that more complex spatial patterns of variation arise in complex terrains such as those in the Sierra Nevada. Moreover, we did not find statistically significant bivariate correlations between population means of different traits. Our research contributes to the understanding of genetic variation in the circadian clock and in seasonal responses in natural populations, highlighting the need for more comprehensive investigations on the association between the clock and other adaptive traits in plants.

Range edges in heterogeneous landscapes: Integrating geographic scale and climate complexity into range dynamics

The impacts of climate change have re-energized interest in understanding the role of climate in setting species geographic range edges. Despite the strong focus on species' distributions in ecology and evolution, defining a species range edge is theoretically and empirically difficult. The challenge of determining a range edge and its relationship to climate is in part driven by the nested nature of geography and the multidimensionality of climate, which together generate complex patterns of both climate and biotic distributions across landscapes. Because range-limiting processes occur in both geographic and climate space, the relationship between these two spaces plays a critical role in setting range limits. With both conceptual and empirical support, we argue that three factors-climate heterogeneity, collinearity among climate variables, and spatial scale-interact to shape the spatial structure of range edges along climate gradients, and we discuss several ways that these factors influence the stability of species range edges with a changing climate. We demonstrate that geographic and climate edges are often not concordant across species ranges. Furthermore, high climate heterogeneity and low climate collinearity across landscapes increase the spectrum of possible relationships between geographic and climatic space, suggesting that geographic range edges and climatic niche limits correspond less frequently than we may expect. More empirical explorations of how the complexity of real landscapes shapes the ecological and evolutionary processes that determine species range edges will advance the development of range limit theory and its applications to biodiversity conservation in the context of changing climate.

Genetic evidence of a northward range expansion in the eastern Bering Sea stock of Pacific cod

Poleward species range shifts have been predicted to result from climate change, and many observations have confirmed such movement. Poleward shifts may represent a homogeneous shift in distribution, seasonal northward movement of specific populations, or colonization processes at the poleward edge of the distribution. The ecosystem of the Bering Sea has been changing along with the climate, moving from an arctic to a subarctic system. Several fish species have been observed farther north than previously reported and in increasing abundances. We examined one of these fish species, Pacific cod, in the northern Bering Sea (NBS) to assess whether they migrated from another stock in the eastern Bering Sea (EBS), Gulf of Alaska, or Aleutian Islands, or whether they represent a separate population. Genetic analyses using 3,599 single nucleotide polymorphism markers indicated that nonspawning cod collected in August 2017 in the NBS were similar to spawning stocks of cod in the EBS. This result suggests escalating northward movement of the large EBS stock during summer months. Whether the cod observed in the NBS migrate south during winter to spawn or remain in the NBS as a sink population is unknown.

Interannual variation in season length is linked to strong co-gradient plasticity of phenology in a montane annual plant

Species are commonly distributed along latitudinal and elevational gradients of growing season length to which they might respond via phenotypic plasticity and/or adaptive genetic differentiation. However, the relative contribution of these processes and whether plasticity, if it occurs, facilitates expansion along season-length gradients remain unclear, but are important for predicting species fates during anthropogenic change. We quantified phenological trait variation in the montane annual Rhinanthus minor for three generations at 12 sites across 900 m of elevation in the Canadian Rocky Mountains and conducted a reciprocal transplant experiment for two generations among nine sites. We compared clines and interannual variation of phenological traits between natural and transplanted individuals. Season length declined by c. 37% along our elevational gradient and, as expected, plants emerged, reached first flower and made their first seed in c. 41% fewer growing degree days under shorter growing seasons. Although reciprocal transplants revealed modest genetic differentiation across elevation, trait clines primarily were due to striking co-gradient plasticity that paralleled genetic differentiation. Co-gradient plasticity likely evolved in response to considerable interannual variation in season length across our elevational transect, and should prepare R. minor to make adaptive changes to phenology in response to ongoing climate change predicted for montane environments.

Landscape Features and Climatic Forces Shape the Genetic Structure and Evolutionary History of an Oak Species (Quercus chenii) in East China

Major topographic features facilitate intraspecific divergence through geographic isolation. This process may be enhanced by environmental isolation along climatic gradients, but also may be reduced by range shifts under rapid climatic changes. In this study, we examined how topography and climate have interacted over time and space to influence the genetic structure and evolutionary history of Quercus chenii, a deciduous oak species representative of the East China flora. Based on the nuclear microsatellite variation at 14 loci, we identified multiple genetic boundaries that were well associated with persistent landscape barriers of East China. Redundancy analysis indicated that both geography and climate explained similar amounts of intraspecific variation. Ecological differences along altitudinal gradients may have driven the divergence between highlands and lowlands. However, range expansions during the Last Interglacial as inferred from approximate Bayesian computation (ABC) may have increased the genetic diversity and eliminated the differentiation of lowland populations via admixture. Chloroplast (cp) DNA analysis of four intergenic spacers (2,866 bp in length) identified a total of 18 haplotypes, 15 of which were private to a single population, probably a result of long-term isolation among multiple montane habitats. A time-calibrated phylogeny suggested that palaeoclimatic changes of the Miocene underlay the lineage divergence of three major clades. In combination with ecological niche modeling (ENM), we concluded that mountainous areas with higher climatic stability are more likely to be glacial refugia that preserved higher phylogenetic diversity, while plains and basins may have acted as dispersal corridors for the post-glacial south-to-north migration. Our findings provide compelling evidence that both topography and climate have shaped the pattern of genetic variation of Q. chenii. Mountains as barriers facilitated differentiation through both geographic and environmental isolation, whereas lowlands as corridors increased the population connectivity especially when the species experienced range expansions.

Evidence for adaptive responses to historic drought across a native plant species range

As climatic conditions change, species will be forced to move or adapt to avoid extinction. Exacerbated by ongoing climate change, California recently experienced a severe and exceptional drought from 2011 to 2017. To investigate whether an adaptive response occurred during this event, we conducted a "resurrection" study of the cutleaf monkeyflower (Mimulus laciniatus), an annual plant, by comparing trait means and variances of ancestral seed collections ("pre-drought") with contemporary descendant collections ("drought"). Plants were grown under common conditions to test whether this geographically restricted species has the capacity to respond evolutionarily to climate stress across its range. We examined if traits shifted in response to the recent, severe drought and included populations across an elevation gradient, including populations at the low- and high-elevation edges of the species range. We found that time to seedling emergence in the drought generation was significantly earlier than in the pre-drought generation, a response consistent with drought adaptation. Additionally, trait variation in days to emergence was reduced in the drought generation, which suggests selection or bottleneck events. Days to first flower increased significantly by elevation, consistent with climate adaptation across the species range. Drought generation plants were larger and had greater reproduction, which was likely a carryover effect of earlier germination. These results demonstrate that rapid shifts in trait means and variances consistent with climate adaptation are occurring within populations, including peripheral populations at warm and cold climate limits, of a plant species with a relatively restricted range that has so far not shifted its elevation distribution during contemporary climate change. Thus, rapid evolution may mitigate, at least temporarily, range shifts under global climate change. This study highlights the need for better understanding rapid adaptation as a means for plant communities to cope with extraordinary climate events.

The Tanggula Mountains enhance population divergence in Carex moorcroftii: a dominant sedge on the Qinghai-Tibetan Plateau

High-altitude mountains are often geographic barriers to gene flow and play important roles in shaping population divergence. The central Qinghai-Tibetan Plateau (QTP) stands the location of the Tanggula Mountains (TM). We use the TM as a case, using Carex moorcroftii, a dominant species on the QTP to test the effects of geographic barriers on plant population divergence. We sampled 18 C. moorcroftii populations along a north-south transect crossing the TM to investigate the correlations of genetic variation and morphological traits with climate variables. The results showed this species holds high genetic diversity (H_e = 0.58) and the surveyed populations can be genetically clustered into two groups: populations from the north face of TM, and the other from the south. Gene flow between populations within groups is higher than those between groups. The traits, number and mass of seeds, mass of root and infructescence significantly varied among populations. Mantel-tests detected a weak but significantly positive correlation between genetic and geographic (R² = 0.107, p = 0.032) and climatic distance (R² = 0.162, p = 0.005), indicating both isolation by distance and isolation by environment. These findings together suggest high-altitude mountains of TM interrupt habitat continuity, result in distinct climatic conditions on both sides, increasing population divergence of plant species.

Comparative epigenetic and genetic spatial structure of the perennial herb Helleborus foetidus: Isolation by environment, isolation by distance, and functional trait divergence

PREMISE OF THE STUDY

Epigenetic variation can play a role in local adaptation; thus, there should be associations among epigenetic variation, environmental variation, and functional trait variation across populations. This study examines these relationships in the perennial herb Helleborus foetidus (Ranunculaceae).

METHODS

Plants from 10 subpopulations were characterized genetically (AFLP, SSR markers), epigenetically (MSAP markers), and phenotypically (20 functional traits). Habitats were characterized using six environmental variables. Isolation-by-distance (IBD) and isolation-by-environment (IBE) patterns of genetic and epigenetic divergence were assessed, as was the comparative explanatory value of geographical and environmental distance as predictors of epigenetic, genetic, and functional differentiation.

KEY RESULTS

Subpopulations were differentiated genetically, epigenetically, and phenotypically. Genetic differentiation was best explained by geographical distance, while epigenetic differentiation was best explained by environmental distance. Divergence in functional traits was correlated with environmental and epigenetic distances, but not with geographical and genetic distances.

CONCLUSIONS

Results are compatible with the hypothesis that epigenetic IBE and functional divergence reflected responses to environmental variation. Spatial analyses simultaneously considering epigenetic, genetic, phenotypic and environmental information provide a useful tool to evaluate the role of environmental features as drivers of natural epigenetic variation between populations.

Connecting genomic patterns of local adaptation and niche suitability in teosintes

The central abundance hypothesis predicts that local adaptation is a function of the distance to the centre of a species' geographic range. To test this hypothesis, we gathered genomic diversity data from 49 populations, 646 individuals and 33,464 SNPs of two wild relatives of maize, the teosintes Zea mays ssp. parviglumis and Zea. mays. ssp. mexicana. We examined the association between the distance to their climatic and geographic centroids and the enrichment of SNPs bearing signals of adaptation. We identified candidate adaptive SNPs in each population by combining neutrality tests and cline analyses. By applying linear regression models, we found that the number of candidate SNPs is positively associated with niche suitability, while genetic diversity is reduced at the limits of the geographic distribution. Our results suggest that overall, populations located at the limit of the species' niches are adapting locally. We argue that local adaptation to this limit could initiate ecological speciation processes and facilitate adaptation to global change.