Genomic signatures of ecological divergence between savanna and forest populations of a Neotropical tree.

BACKGROUND AND AIMS: In eastern Neotropical South America, the Cerrado, a large savanna vegetation, and the Atlantic Forest harbour high biodiversity levels, and their habitats are rather different from each other. The biomes have intrinsic evolutionary relationships, with high lineage exchange that can be attributed, in part, to a large contact zone between them. The genomic study of ecotypes, i.e. populations adapted to divergent habitats, can be a model to study the genomic signatures of ecological divergence. Here, we investigated two ecotypes of the tree Plathymenia reticulata, one from the Cerrado and the other from the Atlantic Forest, which have a hybrid zone in the ecotonal zone of Atlantic Forest-Cerrado. METHODS: The ecotypes were sampled in the two biomes and their ecotone. The evolutionary history of the divergence of the species was analysed with double-digest restriction site-associated DNA sequencing. The genetic structure and the genotypic composition of the hybrid zone were determined. Genotype-association analyses were performed, and the loci under putative selection and their functions were investigated. KEY RESULTS: High divergence between the two ecotypes was found, and only early-generation hybrids were found in the hybrid zone, suggesting a partial reproductive barrier. Ancient introgression between the Cerrado and Atlantic Forest was not detected. The soil and climate were associated with genetic divergence in Plathymenia ecotypes and outlier loci were found to be associated with the stress response, with stomatal and root development and with reproduction. CONCLUSIONS: The high genomic, ecological and morphophysiological divergence between ecotypes, coupled with partial reproductive isolation, indicate that the ecotypes represent two species and should be managed as different evolutionary lineages. We advise that the forest species should be re-evaluated and restated as vulnerable. Our results provide insights into the genomic mechanisms underlying the diversification of species across savanna and forest habitats and the evolutionary forces acting in the species diversification in the Neotropics.

Historical biogeography and local adaptation explain population genetic structure in a widespread terrestrial orchid.

BACKGROUND AND AIMS: Historical changes in environmental conditions and colonization-extinction dynamics have a direct impact on the genetic structure of plant populations. However, understanding how past environmental conditions influenced the evolution of species with high gene flow is challenging when signals for genetic isolation and adaptation are swamped by gene flow. We investigated the spatial distribution and genetic structure of the widespread terrestrial orchid Epipactis helleborine to identify glacial refugia, characterize postglacial population dynamics and assess its adaptive potential. METHODS: Ecological niche modelling was used to locate possible glacial refugia and postglacial recolonization opportunities of E. helleborine. A large single-nucleotide polymorphism (SNP) dataset obtained through genotyping by sequencing was used to define population genetic diversity and structure and to identify sources of postglacial gene flow. Outlier analyses were used to elucidate how adaptation to the local environment contributed to population divergence. KEY RESULTS: The distribution of climatically suitable areas was restricted during the Last Glacial Maximum to the Mediterranean, south-western Europe and small areas in the Alps and Carpathians. Within-population genetic diversity was high in E. helleborine (mean expected heterozygosity, 0.373 ± 0.006; observed heterozygosity, 0.571 ± 0.012; allelic richness, 1.387 ± 0.007). Italy and central Europe are likely to have acted as important genetic sources during postglacial recolonization. Adaptive SNPs were associated with temperature, elevation and precipitation. CONCLUSIONS: Forests in the Mediterranean and Carpathians are likely to have acted as glacial refugia for Epipactis helleborine. Postglacial migration northwards and to higher elevations resulted in the dispersal and diversification of E. helleborine in central Europe and Italy, and to geographical isolation and divergent adaptation in Greek and Italian populations. Distinguishing adaptive from neutral genetic diversity allowed us to conclude that E. helleborine has a high adaptive potential to climate change and demonstrates that signals of adaptation and historical isolation can be identified even in species with high gene flow.

Population genomics of a forest fungus reveals high gene flow and climate adaptation signatures.

Genome sequencing of spatially distributed individuals sheds light on how evolution structures genetic variation. Populations of Phellopilus nigrolimitatus, a red-listed wood-inhabiting fungus associated with old-growth coniferous forests, have decreased in size over the last century due to a loss of suitable habitats. We assessed the population genetic structure and investigated local adaptation in P. nigrolimitatus, by establishing a reference genome and genotyping 327 individuals sampled from 24 locations in Northern Europe by RAD sequencing. We revealed a shallow population genetic structure, indicating large historical population sizes and high levels of gene flow. Despite this weak substructuring, two genetic groups were recognized; a western group distributed mostly in Norway and an eastern group covering most of Finland, Poland and Russia. This substructuring may reflect coimmigration with the main host, Norway spruce (Picea abies), into Northern Europe after the last ice age. We found evidence of low levels of genetic diversity in southwestern Finland, which has a long history of intensive forestry and urbanization. Numerous loci were significantly associated with one or more environmental factors, indicating adaptation to specific environments. These loci clustered into two groups with different associations with temperature and precipitation. Overall, our findings indicate that the current population genetic structure of P. nigrolimitatus results from a combination of gene flow, genetic drift and selection. The acquisition of similar knowledge especially over broad geographic scales, linking signatures of adaptive genetic variation to evolutionary processes and environmental variation, for other fungal species will undoubtedly be useful for assessment of the combined effects of habitat fragmentation and climate change on fungi strongly bound to old-growth forests.

Finding the adaptive needles in a population-structured haystack: A case study in a New Zealand mollusc.

Genetic adaptation to future environmental conditions is crucial to help species persist as the climate changes. Genome scans are powerful tools to understand adaptive landscapes, enabling us to correlate genetic diversity with environmental gradients while disentangling neutral from adaptive variation. However, low gene flow can lead to both local adaptation and highly structured populations, and is a major confounding factor for genome scans, resulting in an inflated number of candidate loci. Here, we compared candidate locus detection in a marine mollusc (Onithochiton neglectus), taking advantage of a natural geographical contrast in the levels of genetic structure between its populations. O. neglectus is endemic to New Zealand and distributed throughout an environmental gradient from the subtropical north to the subantarctic south. Due to a brooding developmental mode, populations tend to be locally isolated. However, adult hitchhiking on rafting kelp increases connectivity among southern populations. We applied two genome scans for outliers (Bayescan and PCAdapt) and two genotype-environment association (GEA) tests (BayeScEnv and RDA). To limit issues with false positives, we combined results using the geometric mean of q-values and performed association tests with random environmental variables. This novel approach is a compromise between stringent and relaxed approaches widely used before, and allowed us to classify candidate loci as low confidence or high confidence. Genome scans for outliers detected a large number of significant outliers in strong and moderately structured populations. No high-confidence GEA loci were detected in the context of strong population structure. However, 86 high-confidence loci were associated predominantly with latitudinally varying abiotic factors in the less structured southern populations. This suggests that the degree of connectivity driven by kelp rafting over the southern scale may be insufficient to counteract local adaptation in this species. Our study supports the expectation that genome scans may be prone to errors in highly structured populations. Nonetheless, it also empirically demonstrates that careful statistical controls enable the identification of candidate loci that invite more detailed investigations. Ultimately, genome scans are valuable tools to help guide further research aiming to determine the potential of non-model species to adapt to future environments.

Genotype-Environment mismatch of kelp forests under climate change.

Climate change is increasingly impacting ecosystems globally. Understanding adaptive genetic diversity and whether it will keep pace with projected climatic change is necessary to assess species' vulnerability and design efficient mitigation strategies such as assisted adaptation. Kelp forests are the foundations of temperate reefs globally but are declining in many regions due to climate stress. A lack of knowledge of kelp's adaptive genetic diversity hinders assessment of vulnerability under extant and future climates. Using 4245 single nucleotide polymorphisms (SNPs), we characterized patterns of neutral and putative adaptive genetic diversity for the dominant kelp in the southern hemisphere (Ecklonia radiata) from ~1000 km of coastline off Western Australia. Strong population structure and isolation-by-distance was underpinned by significant signatures of selection related to temperature and light. Gradient forest analysis of temperature-linked SNPs under selection revealed a strong association with mean annual temperature range, suggesting adaptation to local thermal environments. Critically, modelling revealed that predicted climate-mediated temperature changes will probably result in high genomic vulnerability via a mismatch between current and future predicted genotype-environment relationships such that kelp forests off Western Australia will need to significantly adapt to keep pace with projected climate change. Proactive management techniques such as assisted adaptation to boost resilience may be required to secure the future of these kelp forests and the immense ecological and economic values they support.

Genetic signatures of divergent selection in European beech (Fagus sylvatica L.) are associated with the variation in temperature and precipitation across its distribution range.

High genetic variation and extensive gene flow may help forest trees with adapting to ongoing climate change, yet the genetic bases underlying their adaptive potential remain largely unknown. We investigated range-wide patterns of potentially adaptive genetic variation in 64 populations of European beech (Fagus sylvatica L.) using 270 SNPs from 139 candidate genes involved either in phenology or in stress responses. We inferred neutral genetic structure and processes (drift and gene flow) and performed differentiation outlier analyses and gene-environment association (GEA) analyses to detect signatures of divergent selection. Beech range-wide genetic structure was consistent with the species' previously identified postglacial expansion scenario and recolonization routes. Populations showed high diversity and low differentiation along the major expansion routes. A total of 52 loci were found to be putatively under selection and 15 of them turned up in multiple GEA analyses. Temperature and precipitation related variables were equally represented in significant genotype-climate associations. Signatures of divergent selection were detected in the same proportion for stress response and phenology-related genes. The range-wide adaptive genetic structure of beech appears highly integrated, suggesting a balanced contribution of phenology and stress-related genes to local adaptation, and of temperature and precipitation regimes to genetic clines. Our results imply a best-case scenario for the maintenance of high genetic diversity during range shifts in beech (and putatively other forest trees) with a combination of gene flow maintaining within-population neutral diversity and selection maintaining between-population adaptive differentiation.

Environment and pathogens shape local and regional adaptations to climate change in the chocolate tree, Theobroma cacao L.

Predicting the potential fate of a species in the face of climate change requires knowing the distribution of molecular adaptations across the geographic range of the species. In this work, we analysed 79 genomes of Theobroma cacao, an Amazonian tree known for the fruit from which chocolate is produced, to evaluate how local and regional molecular signatures of adaptation are distributed across the natural range of the species. We implemented novel techniques that incorporate summary statistics from multiple selection scans to infer selective sweeps. The majority of the molecular adaptations in the genome are not shared among populations. We show that ~71.5% of genes under selection also show significant associations with changes in environmental variables. Our results support the interpretation that these genes contribute to local adaptation of the populations in response to abiotic factors. We also found strong patterns of molecular adaptation in a diverse array of disease resistance genes (6.5% of selective sweeps), suggesting that differential adaptation to pathogens also contributes significantly to local adaptations. Our results are consistent with the interpretation that local selective pressures are more important than regional selective pressures in explaining adaptation across the range of a species.

Adaptive introgression as a driver of local adaptation to climate in European white oaks.

Latitudinal and elevational gradients provide valuable experimental settings for studies of the potential impact of global warming on forest tree species. The availability of long-term phenological surveys in common garden experiments for traits associated with climate, such as bud flushing for sessile oaks (Quercus petraea), provide an ideal opportunity to investigate this impact. We sequenced 18 sessile oak populations and used available sequencing data for three other closely related European white oak species (Quercus pyrenaica, Quercus pubescens, and Quercus robur) to explore the evolutionary processes responsible for shaping the genetic variation across latitudinal and elevational gradients in extant sessile oaks. We used phenotypic surveys in common garden experiments and climatic data for the population of origin to perform genome-wide scans for population differentiation and genotype-environment and genotype-phenotype associations. The inferred historical relationships between Q. petraea populations suggest that interspecific gene flow occurred between Q. robur and Q. petraea populations from cooler or wetter areas. A genome-wide scan of differentiation between Q. petraea populations identified single nucleotide polymorphisms (SNPs) displaying strong interspecific relative divergence between these two species. These SNPs followed genetic clines along climatic or phenotypic gradients, providing further support for the likely contribution of introgression to the adaptive divergence of Q. petraea populations. Overall, the results indicate that outliers and associated SNPs are Q. robur ancestry-informative. We discuss the results of this study in the framework of the postglacial colonization scenario, in which introgression and diversifying selection have been proposed as essential drivers of Q. petraea microevolution.

Rapid repeatable phenotypic and genomic adaptation following multiple introductions.

Uncovering the genomic basis of repeated adaption can provide important insights into the constraints and biases that limit the diversity of genetic responses. Demographic processes such as admixture or bottlenecks affect genetic variation underlying traits experiencing selection. The impact of these processes on the genetic basis of adaptation remains, however, largely unexamined empirically. We here test repeatability in phenotypes and genotypes along parallel climatic clines within the native North American and introduced European and Australian Ambrosia artemisiifolia ranges. To do this, we combined multiple lines of evidence from phenotype-environment associations, FST -like outlier tests, genotype-environment associations and genotype-phenotype associations. We used 853 individuals grown in common garden from 84 sampling locations, targeting 19 phenotypes, >83 k SNPs and 22 environmental variables. We found that 17%-26% of loci with adaptive signatures were repeated among ranges, despite alternative demographic histories shaping genetic variation and genetic associations. Our results suggest major adaptive changes can occur on short timescales, with seemingly minimum impacts due to demographic changes linked to introduction. These patterns reveal some predictability of evolutionary change during range expansion, key in a world facing ongoing climate change, and rapid invasive spread.

Landscape genomic signatures indicate reduced gene flow and forest-associated adaptive divergence in an endangered neotropical turtle.

Human-induced transformations of ecosystems usually result in fragmented populations subject to increased extinction risk. Fragmentation is also often associated with novel environmental heterogeneity, which in combination with restricted gene flow may increase the opportunity for local adaptation. To manage at-risk populations in these landscapes, it is important to understand how gene flow is changing, and how populations respond to habitat loss. We conducted a landscape genomics analysis using Restriction-site Associated DNA sequencing to investigate the evolutionary response of the critically endangered Dahl's Toad-headed turtle (Mesoclemmys dahli) to severe habitat modification. The species has lost almost all of its natural habitat in the southwestern part of its range and about 70% in the northeast. Based on least cost path analysis across different resistance surfaces for 3,211 SNPs, we found that the landscape matrix is restricting gene flow, causing the fragmentation of the species into at least six populations. Genome scans and allele-environment association analyses indicate that the population fragments in the deforested grasslands of the southwest are adaptively different from those in the more forested northeast. Populations in areas with no forest had low levels of adaptive genetic diversity and the fixation of ancestrally-polymorphic SNPs, consistent with directional selection in this novel environment. Our results suggest that this forest-stream specialist is adapting to pond-grassland conditions, but it is also suffering from negative consequences of habitat loss, including genetic erosion, isolation, small effective population sizes, and inbreeding. We recommend gene flow restoration via genetic rescue to counteract these threats, and provide guidance for this strategy.

Distinct genecological patterns in seedlings of Norway spruce and silver fir from a mountainous landscape.

Understanding the genecology of forest trees is critical for gene conservation, for predicting the effects of climate change and climate change adaptation, and for successful reforestation. Although common genecological patterns have emerged, species-specific details are also important. Which species are most vulnerable to climate change? Which are the most important adaptive traits and environmental drivers of natural selection? Even though species have been classified as adaptive specialists vs. adaptive generalists, large-scale studies comparing different species in the same experiment are rare. We studied the genecology of Norway spruce (Picea abies) and silver fir (Abies alba), two co-occurring but ecologically distinct European conifers in Central Europe. For each species, we collected seed from more than 90 populations across Switzerland, established a seedling common-garden test, and developed genecological models that associate population variation in seedling growth and phenology to climate, soil properties, and site water balance. Population differentiation and associations between seedling traits and environmental variables were much stronger for Norway spruce than for silver fir, and stronger for seedling height growth than for bud phenology. In Norway spruce, height growth and second flushing were strongly associated with temperature and elevation, with seedlings from the lowlands being taller and more prone to second flush than seedlings from the Alps. In silver fir, height growth was more weakly associated with temperature and elevation, but also associated with water availability. Soil characteristics explained little population variation in both species. We conclude that Norway spruce has become an adaptive specialist because trade-offs between rapid juvenile growth and frost avoidance have subjected it to strong diversifying natural selection based on temperature. In contrast, because silver fir has a more conservative growth habit, it has evolved to become an adaptive generalist. This study demonstrates that co-occurring tree species can develop very different adaptive strategies under identical environmental conditions, and suggests that Norway spruce might be more vulnerable to future maladaptation due to rapid climate change than silver fir.

Finding needles in a genomic haystack: targeted capture identifies clear signatures of selection in a nonmodel plant species.

Teasing apart neutral and adaptive genomic processes and identifying loci that are targets of selection can be difficult, particularly for nonmodel species that lack a reference genome. However, identifying such loci and the factors driving selection have the potential to greatly assist conservation and restoration practices, especially for the management of species in the face of contemporary and future climate change. Here, we focus on assessing adaptive genomic variation within a nonmodel plant species, the narrow-leaf hopbush (Dodonaea viscosa ssp. angustissima), commonly used for restoration in Australia. We used a hybrid-capture target enrichment approach to selectively sequence 970 genes across 17 populations along a latitudinal gradient from 30°S to 36°S. We analysed 8462 single-nucleotide polymorphisms (SNPs) for FST outliers as well as associations with environmental variables. Using three different methods, we found 55 SNPs with significant correlations to temperature and water availability, and 38 SNPs to elevation. Genes containing SNPs identified as under environmental selection were diverse, including aquaporin and abscisic acid genes, as well as genes with ontologies relating to responses to environmental stressors such as water deprivation and salt stress. Redundancy analysis demonstrated that only a small proportion of the total genetic variance was explained by environmental variables. We demonstrate that selection has led to clines in allele frequencies in a number of functional genes, including those linked to leaf shape and stomatal variation, which have been previously observed to vary along the sampled environmental cline. Using our approach, gene regions subject to environmental selection can be readily identified for nonmodel organisms.

Intraspecific variation in stomatal traits, leaf traits and physiology reflects adaptation along aridity gradients in a South African shrub.

BACKGROUND AND AIMS: Trait-environment relationships are commonly interpreted as evidence for local adaptation in plants. However, even when selection analyses support this interpretation, the mechanisms underlying differential benefits are often unknown. This study addresses this gap in knowledge using the broadly distributed South African shrub Protea repens. Specifically, the study examines whether broad-scale patterns of trait variation are consistent with spatial differences in selection and ecophysiology in the wild. METHODS: In a common garden study of plants sourced from 19 populations, associations were measured between five morphological traits and three axes describing source climates. Trait-trait and trait-environment associations were analysed in a multi-response model. Within two focal populations in the wild, selection and path analyses were used to test associations between traits, fecundity and physiological performance. KEY RESULTS: Across 19 populations in a common garden, stomatal density increased with the source population's mean annual temperature and decreased with its average amount of rainfall in midsummer. Concordantly, selection analysis in two natural populations revealed positive selection on stomatal density at the hotter, drier site, while failing to detect selection at the cooler, moister site. Dry-site plants with high stomatal density also had higher stomatal conductances, cooler leaf temperatures and higher light-saturated photosynthetic rates than those with low stomatal density, but no such relationships were present among wet-site plants. Leaf area, stomatal pore index and specific leaf area in the garden also co-varied with climate, but within-population differences were not associated with fitness in either wild population. CONCLUSIONS: The parallel patterns of broad-scale variation, differences in selection and differences in trait-ecophysiology relationships suggest a mechanism for adaptive differentiation in stomatal density. Densely packed stomata may improve performance by increasing transpiration and cooling, but predominately in drier, hotter climates. This study uniquely shows context-dependent benefits of stomatal density--a trait rarely linked to local adaptation in plants.

Climate change vulnerability and conservation strategies for tertiary relict tree species: Insights from landscape genomics of Taxus cuspidata.

The unprecedented habitat fragmentation or loss has threatened the existence of many species. Therefore, it is essential to understand whether and how these species can pace with the environmental changes. Recent advantages in landscape genomics enabled us to identify molecular signatures of adaptation and predict how populations will respond to changing environments, providing new insights into the conservation of species. Here, we investigated the pattern of neutral and putative adaptive genetic variation and its response to changing environments in a tertiary relict tree species, Taxus cuspidata Sieb. et Zucc, which is distributed in northeast China and adjacent regions. We investigated the pattern of genetic diversity and differentiation using restriction site-associated DNA sequencing (RAD-seq) and seven nuclear microsatellites (nSSRs) datasets. We further explored the endangered mechanism, predicted its vulnerability in the future, and provided guidelines for the conservation and management of this species. RAD-seq identified 16,087 single nucleotide polymorphisms (SNPs) in natural populations. Both the SNPs and nSSRs datasets showed high levels of genetic diversity and low genetic differentiation in T. cuspidata. Outlier detection by F ST outlier analysis and genotype-environment associations (GEAs) revealed 598 outlier SNPs as putative adaptive SNPs. Linear redundancy analysis (RDA) and nonlinear gradient forest (GF) showed that the contribution of climate to genetic variation was greater than that of geography, and precipitation played an important role in putative adaptive genetic variation. Furthermore, the genetic offset and risk of non-adaptedness (RONA) suggested that the species at the northeast edge may be more vulnerable in the future. These results suggest that although the species has maintained high current genetic diversity in the face of recent habitat loss and fragmentation, future climate change is likely to threaten the survival of the species. Temperature (Bio03) and precipitation (Prec05) variables can be potentially used as predictors of response of T. cuspidata under future climate. Together, this study provides a theoretical framework for conservation and management strategies for wildlife species in the context of future climate change.

The impact of geography and climate on the population structure and local adaptation in a wild bee.

Deciphering processes that contribute to genetic differentiation and divergent selection of natural populations is useful for evaluating the adaptive potential and resilience of organisms faced with various anthropogenic stressors. Insect pollinator species, including wild bees, provide critical ecosystem services but are highly susceptible to biodiversity declines. Here, we use population genomics to infer the genetic structure and test for evidence of local adaptation in an economically important native pollinator, the small carpenter bee (Ceratina calcarata). Using genome-wide SNP data (n = 8302), collected from specimens across the species' entire distribution, we evaluated population differentiation and genetic diversity and identified putative signatures of selection in the context of geographic and environmental variation. Results of the analyses of principal component and Bayesian clustering were concordant with the presence of two to three genetic clusters, associated with landscape features and inferred phylogeography of the species. All populations examined in our study demonstrated a heterozygote deficit, along with significant levels of inbreeding. We identified 250 robust outlier SNPs, corresponding to 85 annotated genes with known functional relevance to thermoregulation, photoperiod, and responses to various abiotic and biotic stressors. Taken together, these data provide evidence for local adaptation in a wild bee and highlight genetic responses of native pollinators to landscape and climate features.

Genetic and environmental drivers of migratory behavior in western burrowing owls and implications for conservation and management.

Migration is driven by a combination of environmental and genetic factors, but many questions remain about those drivers. Potential interactions between genetic and environmental variants associated with different migratory phenotypes are rarely the focus of study. We pair low coverage whole genome resequencing with a de novo genome assembly to examine population structure, inbreeding, and the environmental factors associated with genetic differentiation between migratory and resident breeding phenotypes in a species of conservation concern, the western burrowing owl (Athene cunicularia hypugaea). Our analyses reveal a dichotomy in gene flow depending on whether the population is resident or migratory, with the former being genetically structured and the latter exhibiting no signs of structure. Among resident populations, we observed significantly higher genetic differentiation, significant isolation-by-distance, and significantly elevated inbreeding. Among migratory breeding groups, on the other hand, we observed lower genetic differentiation, no isolation-by-distance, and substantially lower inbreeding. Using genotype-environment association analysis, we find significant evidence for relationships between migratory phenotypes (i.e., migrant versus resident) and environmental variation associated with cold temperatures during the winter and barren, open habitats. In the regions of the genome most differentiated between migrants and residents, we find significant enrichment for genes associated with the metabolism of fats. This may be linked to the increased pressure on migrants to process and store fats more efficiently in preparation for and during migration. Our results provide a significant contribution toward understanding the evolution of migratory behavior and vital insight into ongoing conservation and management efforts for the western burrowing owl.

Gene×environment associations in orofacial clefting.

This chapter reviews the evidence of gene×environment interactions (G×E) in the etiology of orofacial cleft birth defects (OFCs), specifically cleft lip (CL), cleft palate (CP), and cleft lip with or without cleft palate (CL/P). We summarize the current state of our understanding of the genetic architecture of nonsyndromic OFCs and the evidence that maternal exposures during pregnancy influence risk of OFCs. Further, we present possible candidate gene pathways for these exposures including metabolism of folates, metabolism of retinoids, retinoic acid receptor signaling, aryl hydrocarbon receptor signaling, glucocorticoid receptor signaling, and biotransformation and transport. We review genes in these pathways with prior evidence of association with OFCs, genes with evidence from prior candidate gene G×E studies, and genes identified from genome-wide searches specifically for identifying G×E. Finally, we suggest future directions for G×E research in OFCs.

Re-thinking the environment in landscape genomics.

Detecting the extrinsic selective pressures shaping genomic variation is critical for a better understanding of adaptation and for forecasting evolutionary responses of natural populations to changing environmental conditions. With increasing availability of geo-referenced environmental data, landscape genomics provides unprecedented insights into how genomic variation and underlying gene functions affect traits potentially under selection. Yet, the robustness of genotype-environment associations used in landscape genomics remains tempered due to various limitations, including the characteristics of environmental data used, sampling designs employed, and statistical frameworks applied. Here, we argue that using complementary or new environmental data sources and well-informed sampling designs may help improve the detection of selective pressures underlying patterns of local adaptation in various organisms and environments.

Population genomic analysis reveals geographic structure and climatic diversification for Macrophomina phaseolina isolated from soybean and dry bean across the United States, Puerto Rico, and Colombia.

Macrophomina phaseolina causes charcoal rot, which can significantly reduce yield and seed quality of soybean and dry bean resulting from primarily environmental stressors. Although charcoal rot has been recognized as a warm climate-driven disease of increasing concern under global climate change, knowledge regarding population genetics and climatic variables contributing to the genetic diversity of M. phaseolina is limited. This study conducted genome sequencing for 95 M. phaseolina isolates from soybean and dry bean across the continental United States, Puerto Rico, and Colombia. Inference on the population structure using 76,981 single nucleotide polymorphisms (SNPs) revealed that the isolates exhibited a discrete genetic clustering at the continental level and a continuous genetic differentiation regionally. A majority of isolates from the United States (96%) grouped in a clade with a predominantly clonal genetic structure, while 88% of Puerto Rican and Colombian isolates from dry bean were assigned to a separate clade with higher genetic diversity. A redundancy analysis (RDA) was used to estimate the contributions of climate and spatial structure to genomic variation (11,421 unlinked SNPs). Climate significantly contributed to genomic variation at a continental level with temperature seasonality explaining the most variation while precipitation of warmest quarter explaining the most when spatial structure was accounted for. The loci significantly associated with multivariate climate were found closely to the genes related to fungal stress responses, including transmembrane transport, glycoside hydrolase activity and a heat-shock protein, which may mediate climatic adaptation for M. phaseolina. On the contrary, limited genome-wide differentiation among populations by hosts was observed. These findings highlight the importance of population genetics and identify candidate genes of M. phaseolina that can be used to elucidate the molecular mechanisms that underly climatic adaptation to the changing climate.

Weak genetic differentiation but strong climate-induced selective pressure toward the rear edge of mountain pine in north-eastern Spain.

Local differentiation at distribution limits may influence species' adaptive capacity to environmental changes. However, drivers, such gene flow and local selection, are still poorly understood. We focus on the role played by range limits in mountain forests to test the hypothesis that relict tree populations are subjected to genetic differentiation and local adaptation. Two alpine treelines of mountain pine (Pinus uncinata Ram. ex DC) were investigated in the Spanish Pyrenees. Further, an isolated relict population forming the species' southernmost distribution limit in north-eastern Spain was also investigated. Using genotyping by sequencing, a genetic matrix conformed by single nucleotide polymorphisms (SNPs) was obtained. This matrix was used to perform genotype-environment and genotype-phenotype associations, as well as to model risk of non-adaptedness. Increasing climate seasonality appears as an essential element in the interpretation of SNPs subjected to selective pressures. Genetic differentiations were overall weak. The differences in leaf mass area and radial growth rate, as well as the identification of several SNPs subjected to selective pressures, exceeded neutral predictions of differentiation among populations. Despite genetic drift might prevail in the isolated population, the Fst values (0.060 and 0.066) showed a moderate genetic drift and Nm values (3.939 and 3.555) indicate the presence of gene flow between the relict population and both treelines. Nonetheless, the SNPs subjected to selection pressures provide evidences of possible selection in treeline ecotones. Persistence in range boundaries seems to involve several selective pressures in species' traits, which were significantly related to enhanced drought seasonality at the limit of P. uncinata distribution range. We conclude that gene flow is unlikely to constrain adaptation in the P. uncinata rear edge, although this species shows vulnerability to future climate change scenarios involving warmer and drier conditions.