Genomics of local adaptation with gene flow

Barriers and corridors: Assessment of gene flow and movement among red panda populations in eastern Himalayas

Landscape features can impede dispersal, gene flow, and population demography, resulting in the formation of several meta-populations within a continuous landscape. Understanding a species' ability to overcome these barriers is critical for predicting genetic connectivity and population persistence, and implementing effective conservation strategies. In the present study, we conducted a fine-scale spatial genetic analysis to understand the contemporary gene flow within red panda populations in the Eastern Himalayas. Employing geometric aspects of reserve design, we delineated the critical core habitats for red pandas, which comprise 14.5 % of the landscape (12,189.75 Km2), with only a mere 443 Km2 falling within the protected areas. We identified corridors among the core habitats, which may be vital for the species' long-term genetic viability. Furthermore, we identified substantial landscape barriers, including Sela Pass in the western region, Siang river in the central region, and the Dibang river, Lohit river, along with Dihang, Dipher, and Kumjawng passes in the eastern region, which hinder gene flow. We suggest managing red panda populations through the creation of Community Conservation Reserves in the identified core habitats, following landscape-level management planning based on the core principles of geometric reserve design. This includes a specific emphasis on identified core habitats of red panda (CH-RP 5 and CH-RP 8) to facilitate corridors and implement meta-population dynamics. We propose the development of a comprehensive, long-term conservation and management plan for red pandas in the transboundary landscape, covering China, Nepal, and Bhutan.

Does selection favour the maintenance of porous species boundaries?

The endpoint of speciation has been viewed as complete isolation and the absence of gene flow between species. If the influx of genes from another species is maladaptive because species have different adaptations and genetic backgrounds, selection should favour the closing of species boundaries and zero gene flow, a process known as reinforcement. Recently, numerous cases of gene flow between species have been identified, many of which involved adaptive introgression of beneficial alleles. These cases could reflect transient states on the way to closed species boundaries or the result of declining strength or efficacy of selection for reinforcement as the level of gene flow approaches zero. An alternative hypothesis, however, is that selection favours porous species boundaries that allow beneficial alleles to cross, especially in changing environments. This perspective evaluates the conditions that would be needed for selection to favour porous species boundaries and the evidence for them. A contrast is made between hybridization in sexual eukaryotes and gene transfer via homologous recombination in bacteria. Current evidence is inconclusive on whether non-zero gene flow is favoured by selection. Studies are needed that quantify selection gradients on rates of gene flow and test for evolution towards intermediate values, especially experiments that manipulate conditions and track evolution for multiple generations. Estimation of gene flow networks for more clades and regional assemblages using emerging genome data will also allow the evolutionary determinants of interspecific gene flow to be better understood.

A cryptic syngameon within Betula shrubs revealed: Implications for conservation in changing subarctic environments

Arctic and subarctic ecosystems are rapidly transforming due to global warming, emphasizing the need to understand the genetic diversity and adaptive strategies of northern plant species for effective conservation. This study focuses on Betula glandulosa, a native North American tundra shrub known as dwarf birch, which demonstrates an apparent capacity to adapt to changing climate conditions. To address the taxonomic challenges associated with shrub birches and logistical complexities of sampling in the northernmost areas where species' ranges overlap, we adopted a multicriteria approach. Incorporating molecular data, ploidy level assessment and leaf morphology, we aimed to distinguish B. glandulosa individuals from other shrub birch species sampled. Our results revealed three distinct species and their hybrids within the 537 collected samples, suggesting the existence of a shrub birch syngameon, a reproductive network of interconnected species. Additionally, we identified two discrete genetic clusters within the core species, B. glandulosa, that likely correspond to two different glacial lineages. A comparison between the nuclear and chloroplast SNP data emphasizes a long history of gene exchange between different birch species and genetic clusters. Furthermore, our results highlight the significance of incorporating interfertile congeneric species in conservation strategies and underscores the need for a holistic approach to conservation in the context of climate change, considering the complex dynamics of species interactions. While further research will be needed to describe this shrub birches syngameon and its constituents, this study is a first step in recognizing its existence and disseminating awareness among ecologists and conservation practitioners. This biological phenomenon, which offers evolutionary flexibility and resilience beyond what its constituent species can achieve individually, may have significant ecological implications.

Patterns of genetic variation and local adaptation of a native herbivore to a lethal invasive plant

Understanding the evolutionary processes that influence fitness is critical to predicting species' responses to selection. Interactions among evolutionary processes including gene flow, drift and the strength of selection can lead to either local adaptation or maladaptation, especially in heterogenous landscapes. Populations experiencing novel environments or resources are ideal for understanding the mechanisms underlying adaptation or maladaptation, specifically in locally co-evolved interactions. We used the interaction between a native herbivore that oviposits on a patchily distributed introduced plant that in turn causes significant mortality to the larvae to test for signatures of local adaptation in areas where the two co-occurred. We used whole-genome sequencing to explore population structure, patterns of gene flow and signatures of local adaptation. We found signatures of local adaptation in response to the introduced plant in the absence of strong population structure with no genetic differentiation and low genetic variation. Additionally, we found localized allele frequency differences within a single population between habitats with and without the lethal plant, highlighting the effects of strong selection. Finally, we identified that selection was acting on larval ability to feed on the plant rather than on females' ability to avoid oviposition, thus uncovering the specific ontogenetic target of selection. Our work highlights the potential for adaptation to occur in a fine-grained landscape in the presence of gene flow and low genetic variation.

Multigenerational Fitness Effects of Natural Immigration Indicate Strong Heterosis and Epistatic Breakdown in a Wild Bird Population

AbstractThe fitness of immigrants and their descendants produced within recipient populations fundamentally underpins the genetic and population dynamic consequences of immigration. Immigrants can in principle induce contrasting genetic effects on fitness across generations, reflecting multifaceted additive, dominance, and epistatic effects. Yet full multigenerational and sex-specific fitness effects of regular immigration have not been quantified within naturally structured systems, precluding inference on underlying genetic architectures and population outcomes. We used four decades of song sparrow (Melospiza melodia) life history and pedigree data to quantify fitness of natural immigrants, natives, and their F1, F2, and backcross descendants and test for evidence of nonadditive genetic effects. Values of key fitness components (including adult lifetime reproductive success and zygote survival) of F1 offspring of immigrant-native matings substantially exceeded their parent mean, indicating strong heterosis. Meanwhile, F2 offspring of F1-F1 matings had notably low values, indicating surprisingly strong epistatic breakdown. Furthermore, magnitudes of effects varied among fitness components and differed between female and male descendants. These results demonstrate that strong nonadditive genetic effects on fitness can arise within weakly structured and fragmented populations experiencing frequent natural immigration. Such effects will substantially affect the net degree of effective gene flow and resulting local genetic introgression and adaptation.

Genomic architecture controls multivariate adaptation to climate change

As climate change advances, environmental gradients may decouple, generating novel multivariate environments that stress wild populations. A commonly invoked mechanism of evolutionary rescue is adaptive gene flow tracking climate shifts, but gene flow from populations inhabiting similar conditions on one environmental axis could cause maladaptive introgression when populations are adapted to different environmental variables that do not shift together. Genomic architecture can play an important role in determining the effectiveness and relative magnitudes of adaptive gene flow and in situ adaptation. This may have direct consequences for how species respond to climate change but is often overlooked. Here, we simulated microevolutionary responses to environmental change under scenarios defined by variation in the polygenicity, linkage, and genetic redundancy of two independent traits, one of which is adapted to a gradient that shifts under climate change. We used these simulations to examine how genomic architecture influences evolutionary outcomes under climate change. We found that climate-tracking (up-gradient) gene flow, though present in all scenarios, was strongly constrained under scenarios of lower linkage and higher polygenicity and redundancy, suggesting in situ adaptation as the predominant mechanism of evolutionary rescue under these conditions. We also found that high polygenicity caused increased maladaptation and demographic decline, a concerning result given that many climate-adapted traits may be polygenic. Finally, in scenarios with high redundancy, we observed increased adaptive capacity. This finding adds to the growing recognition of the importance of redundancy in mediating in situ adaptive capacity and suggests opportunities for better understanding the climatic vulnerability of real populations.

Secondary Contact, Introgressive Hybridization, and Genome Stabilization in Sticklebacks

Advances in genomic studies have revealed that hybridization in nature is pervasive and raised questions about the dynamics of different genetic and evolutionary factors following the initial hybridization event. While recent research has proposed that the genomic outcomes of hybridization might be predictable to some extent, many uncertainties remain. With comprehensive whole-genome sequence data, we investigated the genetic introgression between 2 divergent lineages of 9-spined sticklebacks (Pungitius pungitius) in the Baltic Sea. We found that the intensity and direction of selection on the introgressed variation has varied across different genomic elements: while functionally important regions displayed reduced rates of introgression, promoter regions showed enrichment. Despite the general trend of negative selection, we identified specific genomic regions that were enriched for introgressed variants, and within these regions, we detected footprints of selection, indicating adaptive introgression. Geographically, we found the selection against the functional changes to be strongest in the vicinity of the secondary contact zone and weaken as a function of distance from the initial contact. Altogether, the results suggest that the stabilization of introgressed variation in the genomes is a complex, multistage process involving both negative and positive selection. In spite of the predominance of negative selection against introgressed variants, we also found evidence for adaptive introgression variants likely associated with adaptation to Baltic Sea environmental conditions.

Unraveling Adaptive Evolutionary Divergence at Microgeographic Scales

AbstractStriking examples of local adaptation at fine geographic scales are increasingly being documented in natural populations. However, the relative contributions made by natural selection, phenotype-dependent dispersal (when individuals disperse with respect to a habitat preference), and mate preference in generating and maintaining microgeographic adaptation and divergence are not well studied. Here, we develop quantitative genetics models and individual-based simulations (IBSs) to uncover the evolutionary forces that possibly drive microgeographic divergence. We also perform Bayesian estimation of the parameters in our IBS using empirical data on habitat-specific variation in bill morphology in the island scrub-jay (Aphelocoma insularis) to apply our models to a natural system. We find that natural selection and phenotype-dependent dispersal can generate the patterns of divergence we observe in the island scrub-jay. However, mate preference for a mate with similar bill morphology, even though observed in the species, does not play a significant role in driving divergence. Our modeling approach provides insights into phenotypic evolution occurring over small spatial scales relative to dispersal ranges, suggesting that adaptive divergence at microgeographic scales may be common across a wider range of taxa than previously thought. Our quantitative genetic models help to inform future theoretical and empirical work to determine how selection, habitat preference, and mate preference contribute to local adaptation and microgeographic divergence.

Adaptation to distinct habitats is maintained by contrasting selection at different life stages in sunflower ecotypes

Conspecific populations living in adjacent but contrasting microenvironments represent excellent systems for studying natural selection. These systems are valuable because gene flow is expected to force genetic homogeneity except at loci experiencing divergent selection. A history of reciprocal transplant and common garden studies in such systems, and a growing number of genomic studies, have contributed to understanding how selection operates in natural populations. While selection can vary across different fitness components and life stages, few studies have investigated how this ultimately affects allele frequencies and the maintenance of divergence between populations. Here, we study two sunflower ecotypes in distinct, adjacent habitats by combining demographic models with genome-wide sequence data to estimate fitness and allele frequency change at multiple life stages. This framework allows us to estimate that only local ecotypes are likely to experience positive population growth (λ > 1) and that the maintenance of divergent adaptation appears to be mediated via habitat- and life stage-specific selection. We identify genetic variation, significantly driven by loci in chromosomal inversions, associated with different life history strategies in neighbouring ecotypes that optimize different fitness components and may contribute to the maintenance of distinct ecotypes.

Genetically Encoded Lizard Color Divergence for Camouflage and Thermoregulation

Local adaptation is critical in speciation and evolution, yet comprehensive studies on proximate and ultimate causes of local adaptation are generally scarce. Here, we integrated field ecological experiments, genome sequencing, and genetic verification to demonstrate both driving forces and molecular mechanisms governing local adaptation of body coloration in a lizard from the Qinghai-Tibet Plateau. We found dark lizards from the cold meadow population had lower spectrum reflectance but higher melanin contents than light counterparts from the warm dune population. Additionally, the colorations of both dark and light lizards facilitated the camouflage and thermoregulation in their respective microhabitat simultaneously. More importantly, by genome resequencing analysis, we detected a novel mutation in Tyrp1 that underpinned this color adaptation. The allele frequencies at the site of SNP 459# in the gene of Tyrp1 are 22.22% G/C and 77.78% C/C in dark lizards and 100% G/G in light lizards. Model-predicted structure and catalytic activity showed that this mutation increased structure flexibility and catalytic activity in enzyme TYRP1, and thereby facilitated the generation of eumelanin in dark lizards. The function of the mutation in Tyrp1 was further verified by more melanin contents and darker coloration detected in the zebrafish injected with the genotype of Tyrp1 from dark lizards. Therefore, our study demonstrates that a novel mutation of a major melanin-generating gene underpins skin color variation co-selected by camouflage and thermoregulation in a lizard. The resulting strong selection may reinforce adaptive genetic divergence and enable the persistence of adjacent populations with distinct body coloration.

Understanding the joint evolution of dispersal and host specialisation using phytophagous arthropods as a model group

Theory generally predicts that host specialisation and dispersal should evolve jointly. Indeed, many models predict that specialists should be poor dispersers to avoid landing on unsuitable hosts while generalists will have high dispersal abilities. Phytophagous arthropods are an excellent group to test this prediction, given extensive variation in their host range and dispersal abilities. Here, we explore the degree to which the empirical literature on this group is in accordance with theoretical predictions. We first briefly outline the theoretical reasons to expect such a correlation. We then report empirical studies that measured both dispersal and the degree of specialisation in phytophagous arthropods. We find a correlation between dispersal and levels of specialisation in some studies, but with wide variation in this result. We then review theoretical attributes of species and environment that may blur this correlation, namely environmental grain, temporal heterogeneity, habitat selection, genetic architecture, and coevolution between plants and herbivores. We argue that theoretical models fail to account for important aspects, such as phenotypic plasticity and the impact of selective forces stemming from other biotic interactions, on both dispersal and specialisation. Next, we review empirical caveats in the study of this interplay. We find that studies use different measures of both dispersal and specialisation, hampering comparisons. Moreover, several studies do not provide independent measures of these two traits. Finally, variation in these traits may occur at scales that are not being considered. We conclude that this correlation is likely not to be expected from large-scale comparative analyses as it is highly context dependent and should not be considered in isolation from the factors that modulate it, such as environmental scale and heterogeneity, intrinsic traits or biotic interactions. A stronger crosstalk between theoretical and empirical studies is needed to understand better the prevalence and basis of the correlation between dispersal and specialisation.

Divergence and gene flow history at two large chromosomal inversions underlying ecotype differentiation in the long-snouted seahorse

Chromosomal inversions can play an important role in divergence and reproductive isolation by building and maintaining distinct allelic combinations between evolutionary lineages. Alternatively, they can take the form of balanced polymorphisms that segregate within populations until one arrangement becomes fixed. Many questions remain about how inversion polymorphisms arise, how they are maintained over the long term, and ultimately, whether and how they contribute to speciation. The long-snouted seahorse (Hippocampus guttulatus) is genetically subdivided into geographic lineages and marine-lagoon ecotypes, with shared structural variation underlying lineage and ecotype divergence. Here, we aim to characterize structural variants and to reconstruct their history and suspected role in ecotype formation. We generated a near chromosome-level genome assembly and described genome-wide patterns of diversity and divergence through the analysis of 112 whole-genome sequences from Atlantic, Mediterranean, and Black Sea populations. By also analysing linked-read sequencing data, we found evidence for two chromosomal inversions that were several megabases in length and showed contrasting allele frequency patterns between lineages and ecotypes across the species range. We reveal that these inversions represent ancient intraspecific polymorphisms, one likely being maintained by divergent selection and the other by pseudo-overdominance. A possible selective coupling between the two inversions was further supported by the absence of specific haplotype combinations and a putative functional interaction between the two inversions in reproduction. Lastly, we detected gene flux eroding divergence between inverted alleles at varying levels for the two inversions, with a likely impact on their dynamics and contribution to divergence and speciation.

Local adaptation with gene flow in a highly dispersive shark

Adaptive divergence in response to environmental clines are expected to be common in species occupying heterogeneous environments. Despite numerous advances in techniques appropriate for non-model species, gene-environment association studies in elasmobranchs are still scarce. The bronze whaler or copper shark (Carcharhinus brachyurus) is a large coastal shark with a wide distribution and one of the most exploited elasmobranchs in southern Africa. Here, we assessed the distribution of neutral and adaptive genomic diversity in C. brachyurus across a highly heterogeneous environment in southern Africa based on genome-wide SNPs obtained through a restriction site-associated DNA method (3RAD). A combination of differentiation-based genome-scan (outflank) and genotype-environment analyses (redundancy analysis, latent factor mixed models) identified a total of 234 differentiation-based outlier and candidate SNPs associated with bioclimatic variables. Analysis of 26,299 putatively neutral SNPs revealed moderate and evenly distributed levels of genomic diversity across sites from the east coast of South Africa to Angola. Multivariate and clustering analyses demonstrated a high degree of gene flow with no significant population structuring among or within ocean basins. In contrast, the putatively adaptive SNPs demonstrated the presence of two clusters and deep divergence between Angola and all other individuals from Namibia and South Africa. These results provide evidence for adaptive divergence in response to a heterogeneous seascape in a large, mobile shark despite high levels of gene flow. These results are expected to inform management strategies and policy at the national and regional level for conservation of C. brachyurus populations.

Multiple genetic sources facilitate the northward range expansion of an intertidal oyster along China's coast

Coastal artificial structures on the former mudflats provide available habitats for the rocky intertidal species which can establish new populations in these emerging habitats over their former distribution range limits. As a former southern species, the oyster Crassostrea sikamea has become a pioneer and rapidly invaded the artificial shorelines in northern China. We used a seascape genomics approach to investigate the population structure and genetic sources of C. sikamea on the coastal artificial structures, which is crucial for understanding the genetic mechanisms driving species distribution range expansion and invasion pathway of intertidal species. Five C. sikamea populations, including two artificial substrate populations (WGZ and ZAP), one oyster reef population (LS), and two natural rocky shore populations (ZS and XM), were measured using single nucleotide polymorphism (SNPs) obtained from double digest restriction-site associated DNA sequencing (ddRAD-Seq). Redundancy analyses (RDA) were implemented for investigating the relationship between local temperature variables and the temperature adaptability of C. sikamea. Genetic diversity, direction and strength of gene flow, and population structure all revealed that the LS and ZS populations were the genetic sources for the oyster populations on the emerging northern coastal artificial structures. Results of RDA showed that there were different adaptive potentials for northern and southern populations to local temperature variables and the oyster reef population which frequently suffers from heat stress owned high heat adaptability. The ZS population as a genetic source nearby the Yangtze River estuary provided mass larvae for the northern populations, and the other genetic source, the heat-tolerant LS population, in the oyster reef played an important role in the post-settlement success by providing preadapted genotypes. These results highlight the importance of multiple sources with divergent adaptative capabilities for biological invasion, and also emphasize the importance of the oyster reef in coastal biodiversity and conservation.

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

Genomic signals of local adaptation across climatically heterogenous habitats in an invasive tropical fruit fly (Bactrocera tryoni)

Local adaptation plays a key role in the successful establishment of pest populations in new environments by enabling them to tolerate novel biotic and abiotic conditions experienced outside their native range. However, the genomic underpinnings of such adaptive responses remain unclear, especially for agriculturally important pests. We investigated population genomic signatures in the tropical/subtropical Queensland fruit fly, Bactrocera tryoni, which has an expanded range encompassing temperate and arid zones in Australia, and tropical zones in the Pacific Islands. Using reduced representation sequencing data from 28 populations, we detected allele frequency shifts associated with the native/invasive status of populations and identified environmental factors that have likely driven population differentiation. We also determined that precipitation, temperature, and geographic variables explain allelic shifts across the distribution range of B. tryoni. We found spatial heterogeneity in signatures of local adaptation across various climatic conditions in invaded areas. Specifically, disjunct invasive populations in the tropical Pacific Islands and arid zones of Australia were characterised by multiple significantly differentiated single nucleotide polymorphisms (SNPs), some of which were associated with genes with well-understood function in environmental stress (e.g., heat and desiccation) response. However, invasive populations in southeast Australian temperate zones showed higher gene flow with the native range and lacked a strong local adaptive signal. These results suggest that population connectivity with the native range has differentially affected local adaptive patterns in different invasive populations. Overall, our findings provide insights into the evolutionary underpinnings of invasion success of an important horticultural pest in climatically distinct environments.

Coral Reef Population Genomics in an Age of Global Change

Coral reefs are both exceptionally biodiverse and threatened by climate change and other human activities. Here, we review population genomic processes in coral reef taxa and their importance for understanding responses to global change. Many taxa on coral reefs are characterized by weak genetic drift, extensive gene flow, and strong selection from complex biotic and abiotic environments, which together present a fascinating test of microevolutionary theory. Selection, gene flow, and hybridization have played and will continue to play an important role in the adaptation or extinction of coral reef taxa in the face of rapid environmental change, but research remains exceptionally limited compared to the urgent needs. Critical areas for future investigation include understanding evolutionary potential and the mechanisms of local adaptation, developing historical baselines, and building greater research capacity in the countries where most reef diversity is concentrated.

Whole-Genome Sequencing Reveals That Regulatory and Low Pleiotropy Variants Underlie Local Adaptation to Environmental Variability in Purple Sea Urchins

AbstractOrganisms experience environments that vary across both space and time. Such environmental heterogeneity shapes standing genetic variation and may influence species' capacity to adapt to rapid environmental change. However, we know little about the kind of genetic variation that is involved in local adaptation to environmental variability. To address this gap, we sequenced the whole genomes of 140 purple sea urchins (Strongylocentrotus purpuratus) from seven populations that vary in their degree of pH variability. Despite no evidence of global population structure, we found a suite of single-nucleotide polymorphisms (SNPs) tightly correlated with local pH variability (outlier SNPs), which were overrepresented in regions putatively involved in gene regulation (long noncoding RNA and enhancers), supporting the idea that variation in regulatory regions is important for local adaptation to variability. In addition, outliers in genes were found to be (i) enriched for biomineralization and ion homeostasis functions related to low pH response, (ii) less central to the protein-protein interaction network, and (iii) underrepresented among genes highly expressed during early development. Taken together, these results suggest that loci that underlie local adaptation to pH variability in purple sea urchins fall in regions with potentially low pleiotropic effects (based on analyses involving regulatory regions, network centrality, and expression time) involved in low pH response (based on functional enrichment).

Phylogeography of the red alga Gracilariopsis tenuifrons (Gracilariales) along the Brazilian coast

Changes in the sea level during the Holocene are regarded as one of the most prevalent drivers of the diversity and distribution of macroalgae in Brazil, influenced by the emergence of the Vitória-Trindade seamount chain (VTC). Gracilariopsis tenuifrons has a wide geographic distribution along the Brazilian coast, from Maranhão state (2°48'64.3" S) to Santa Catarina state (27.5°73'83" S). The knowledge of historical processes affecting diversity may allow the development of conservation strategies in environments against anthropogenic influence. Therefore, knowledge about phylogeography and populational genetic diversity in G. tenuifrons is necessary. Six populations were sampled along the northeastern tropical (Maranhão-MA, Rio Grande do Norte-RN, Alagoas-AL, and Bahia-BA States) and southeastern subtropical (São Paulo "Ubatuba"-SP1 and São Paulo "Itanhaém"-SP2 States) regions along the Brazilian coast. The genetic diversity and structure of G. tenuifrons were inferred using mitochondrial (COI-5P and cox2-3 concatenated) DNA markers. Gracilariopsis tenuifrons populations showed an evident separation between the northeast (from 2°48'64.3" S to 14°18'23" S; 17 haplotypes) and the southeast (from 23°50'14.9" S to 24°20'04.7" S; 10 haplotypes) regions by two mutational steps between them. The main biogeographical barrier to gene flow is located nearby the VTC. The southeast region (São Paulo State) is separated by two subphylogroups (SP1, three haplotypes and SP2, six haplotypes), and Santos Bay (estuary) has been considered a biogeographical barrier between them. The presence of genetic structure and putative barriers to gene flow are in concordance with previous studies reporting biogeographic breaks in the southwest Atlantic Ocean, including the genetic isolation between northeast and southeast regions for red and brown algae in the vicinity of the VTC.

Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production

Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.

Conserved islands of divergence associated with adaptive variation in sockeye salmon are maintained by multiple mechanisms

Local adaptation is facilitated by loci clustered in relatively few regions of the genome, termed genomic islands of divergence. The mechanisms that create and maintain these islands and how they contribute to adaptive divergence is an active research topic. Here, we use sockeye salmon as a model to investigate both the mechanisms responsible for creating islands of divergence and the patterns of differentiation at these islands. Previous research suggested that multiple islands contributed to adaptive radiation of sockeye salmon. However, the low-density genomic methods used by these studies made it difficult to fully elucidate the mechanisms responsible for islands and connect genotypes to adaptive variation. We used whole genome resequencing to genotype millions of loci to investigate patterns of genetic variation at islands and the mechanisms that potentially created them. We discovered 64 islands, including 16 clustered in four genomic regions shared between two isolated populations. Characterisation of these four regions suggested that three were likely created by structural variation, while one was created by processes not involving structural variation. All four regions were small (< 600 kb), suggesting low recombination regions do not have to span megabases to be important for adaptive divergence. Differentiation at islands was not consistently associated with established population attributes. In sum, the landscape of adaptive divergence and the mechanisms that create it are complex; this complexity likely helps to facilitate fine-scale local adaptation unique to each population.

Genetic Incompatibilities and Evolutionary Rescue by Wild Relatives Shaped Grain Amaranth Domestication

Crop domestication and the subsequent expansion of crops have long been thought of as a linear process from a wild ancestor to a domesticate. However, evidence of gene flow from locally adapted wild relatives that provided adaptive alleles into crops has been identified in multiple species. Yet, little is known about the evolutionary consequences of gene flow during domestication and the interaction of gene flow and genetic load in crop populations. We study the pseudo-cereal grain amaranth that has been domesticated three times in different geographic regions of the Americas. We quantify the amount and distribution of gene flow and genetic load along the genome of the three grain amaranth species and their two wild relatives. Our results show ample gene flow between crop species and between crops and their wild relatives. Gene flow from wild relatives decreased genetic load in the three crop species. This suggests that wild relatives could provide evolutionary rescue by replacing deleterious alleles in crops. We assess experimental hybrids between the three crop species and found genetic incompatibilities between one Central American grain amaranth and the other two crop species. These incompatibilities might have created recent reproductive barriers and maintained species integrity today. Together, our results show that gene flow played an important role in the domestication and expansion of grain amaranth, despite genetic species barriers. The domestication of plants was likely not linear and created a genomic mosaic by multiple contributors with varying fitness effects for today's crops.

Out to sea: ocean currents and patterns of asymmetric gene flow in an intertidal fish species

Passive dispersal via wind or ocean currents can drive asymmetric gene flow, which influences patterns of genetic variation and the capacity of populations to evolve in response to environmental change. The mangrove rivulus fish (Kryptolebias marmoratus), hereafter "rivulus," is an intertidal fish species restricted to the highly fragmented New World mangrove forests of Central America, the Caribbean, the Bahamas, and Florida. Mangrove patches are biological islands with dramatic differences in both abiotic and biotic conditions compared to adjacent habitat. Over 1,000 individual rivulus across 17 populations throughout its range were genotyped at 32 highly polymorphic microsatellites. Range-wide population genetic structure was evaluated with five complementary approaches that found eight distinct population clusters. However, an analysis of molecular variance indicated significant population genetic structure among regions, populations within regions, sampling locations within populations, and individuals within sampling locations, indicating that rivulus has both broad- and fine-scale genetic differentiation. Integrating range-wide genetic data with biophysical modeling based on 10 years of ocean current data showed that ocean currents and the distance between populations over water drive gene flow patterns on broad scales. Directional migration estimates suggested some significant asymmetries in gene flow that also were mediated by ocean currents and distance. Specifically, populations in the center of the range (Florida Keys) were identified as sinks that received migrants (and alleles) from other populations but failed to export individuals. These populations thus harbor genetic variation, perhaps even from extirpated populations across the range, but ocean currents and complex arrangements of landmasses might prevent the distribution of that genetic variation elsewhere. Hence, the inherent asymmetry of ocean currents shown to impact both genetic differentiation and directional migration rates may be responsible for the complex distribution of genetic variation across the range and observed patterns of metapopulation structure.

The molecular background of the aspartate aminotransferase polymorphism in Littorina snails maintained by strong selection on small spatial scales

Allozymes present several classical examples of divergent selection, including the variation in the cytosolic aspartate aminotransferase (AAT) in the intertidal snails Littorina saxatilis. AAT is a part of the asparate-malate shuttle, in the interidal molluscs involved in the anaerobic respiration during desiccation. Previous allozyme studies reported the sharp gradient in the frequencies of the AAT¹⁰⁰and the AAT¹²⁰ alleles between the low and high shores in the Northern Europe and the differences in their enzymatic activity, supporting the role of AAT in adaptation to desiccation. However, the populations in the Iberian Peninsula showed the opposite allele cline. Using the mRNA sequencing and the genome pool-seq analyses we characterize DNA sequences of the different AAT alleles, report the amino acid replacements behind the allozyme variation and show that same allozyme alleles in Northern and Southern populations have different protein sequences. Gene phylogeny reveals that the AAT¹⁰⁰ and the northern AAT¹²⁰ alleles represent the old polymorphism, shared among the closely related species of Littorina, while the southern AAT¹²⁰ allele is more recently derived from AAT¹⁰⁰. Further, we show that the Aat gene is expressed at constitutive level in different genotypes and conditions, supporting the role of structural variation in regulation of enzyme activity. Finally, we report the location and the structure of the gene in the L. saxatilis genome and the presence of two additional non-functional gene copies. Altogether, we provide a missing link between the classical allozyme studies and the genome scans and bring together the results produced over decades of the genetic research.

Growth, drought response, and climate-associated genomic structure in whitebark pine in the Sierra Nevada of California

Whitebark pine (Pinus albicaulis Engelm.) has experienced rapid population declines and is listed as threatened under the Endangered Species Act in the United States. Whitebark pine in the Sierra Nevada of California represents the southernmost end of the species' distribution and, like other portions of its range, faces threats from an introduced pathogen, native bark beetles, and a rapidly warming climate. Beyond these chronic stressors, there is also concern about how this species will respond to acute stressors, such as drought. We present patterns of stem growth from 766 large (average diameter at breast height >25 cm), disease-free whitebark pine across the Sierra Nevada before and during a recent period of drought. We contextualize growth patterns using population genomic diversity and structure from a subset of 327 trees. Sampled whitebark pine generally had positive to neutral stem growth trends from 1970 to 2011, which was positively correlated with minimum temperature and precipitation. Indices of stem growth during drought years (2012 to 2015) relative to a predrought interval were mostly positive to neutral at our sampled sites. Individual tree growth response phenotypes appeared to be linked to genotypic variation in climate-associated loci, suggesting that some genotypes can take better advantage of local climatic conditions than others. We speculate that reduced snowpack during the 2012 to 2015 drought years may have lengthened the growing season while retaining sufficient moisture to maintain growth at most study sites. Growth responses may differ under future warming, however, particularly if drought severity increases and modifies interactions with pests and pathogens.

A comprehensive analysis of the genetic diversity and environmental adaptability in worldwide Merino and Merino-derived sheep breeds

BACKGROUND

To enhance and extend the knowledge about the global historical and phylogenetic relationships between Merino and Merino-derived breeds, 19 populations were genotyped with the OvineSNP50 BeadChip specifically for this study, while an additional 23 populations from the publicly available genotypes were retrieved. Three complementary statistical tests, Rsb (extended haplotype homozygosity between-populations), XP-EHH (cross-population extended haplotype homozygosity), and runs of homozygosity (ROH) islands were applied to identify genomic variants with potential impact on the adaptability of Merino genetic type in two contrasting climate zones.

RESULTS

The results indicate that a large part of the Merino's genetic relatedness and admixture patterns are explained by their genetic background and/or geographic origin, followed by local admixture. Multi-dimensional scaling, Neighbor-Net, Admixture, and TREEMIX analyses consistently provided evidence of the role of Australian, Rambouillet and German strains in the extensive gene introgression into the other Merino and Merino-derived breeds. The close relationship between Iberian Merinos and other South-western European breeds is consistent with the Iberian origin of the Merino genetic type, with traces from previous contributions of other Mediterranean stocks. Using Rsb and XP-EHH approaches, signatures of selection were detected spanning four genomic regions located on Ovis aries chromosomes (OAR) 1, 6 and 16, whereas two genomic regions on OAR6, that partially overlapped with the previous ones, were highlighted by ROH islands. Overall, the three approaches identified 106 candidate genes putatively under selection. Among them, genes related to immune response were identified via the gene interaction network. In addition, several candidate genes were found, such as LEKR1, LCORL, GHR, RBPJ, BMPR1B, PPARGC1A, and PRKAA1, related to morphological, growth and reproductive traits, adaptive thermogenesis, and hypoxia responses.

CONCLUSIONS

To the best of our knowledge, this is the first comprehensive dataset that includes most of the Merino and Merino-derived sheep breeds raised in different regions of the world. The results provide an in-depth picture of the genetic makeup of the current Merino and Merino-derived breeds, highlighting the possible selection pressures associated with the combined effect of anthropic and environmental factors. The study underlines the importance of Merino genetic types as invaluable resources of possible adaptive diversity in the context of the occurring climate changes.

Microgeographic variation in early fitness traits of Pinus sylvestris from contrasting soils

PREMISE

The possibility of fine-scale intraspecific adaptive divergence under gene flow is established by theoretical models and has been confirmed empirically in tree populations distributed along steep altitudinal clines or across extreme edaphic discontinuities. However, the possibility of microgeographic adaptive divergence due to less severe but more frequent kinds of soil variation is unclear.

METHODS

In this study, we looked for evidence of local adaptation to calcareous versus siliceous soil types in two nearby Mediterranean Pinus sylvestris populations connected via pollen flow. Using a greenhouse experiment, we tested for variation in early (up to three years of age) seedling performance among open-pollinated maternal families originating from each edaphic provenance when experimentally grown on both types of natural local substrate.

RESULTS

Although seedlings were clearly affected by the edaphic environment, exhibiting lower and slower emergence as well as higher mortality on the calcareous than in the siliceous substrate, neither the performance on each substrate nor the plasticity among substrates varied significantly with seedling edaphic provenance.

CONCLUSIONS

We found no evidence of local adaptation to a non-extreme edaphic discontinuity over a small spatial scale, at least during early stages of seedling establishment. Future studies on microgeographic soil-driven adaptation should consider long-term experiments to minimize maternal effects and allow a potentially delayed expression of edaphic adaptive divergence.

The impact of habitat loss on molecular signatures of coevolution between an iconic butterfly (Alcon blue) and its host plant (Marsh gentian)

Habitat loss is threatening natural communities worldwide. Small and isolated populations suffer from inbreeding and genetic drift, which jeopardize their long-term survival and adaptive capacities. However, the consequences of habitat loss for reciprocal coevolutionary interactions remain poorly studied. In this study, we investigated the effects of decreasing habitat patch size and connectivity associated with habitat loss on molecular signatures of coevolution in the Alcon blue butterfly (Phengaris alcon) and its most limited host, the marsh gentian (Gentiana pneumonanthe). Because reciprocal coevolution is characterized by negative frequency-dependent selection as a particular type of balancing selection, we investigated how signatures of balancing selection vary along a gradient of patch size and connectivity, using single nucleotide polymorphisms (SNPs). We found that signatures of coevolution were unaffected by patch characteristics in the host plants. On the other hand, more pronounced signatures of coevolution were observed in both spatially isolated and in large Alcon populations, together with pronounced spatial variation in SNPs that are putatively involved in coevolution. These findings suggest that habitat loss can facilitate coevolution in large butterfly populations through limiting swamping of locally beneficial alleles by maladaptive ones. We also found that allelic richness (Ar) of the coevolutionary SNPs is decoupled from neutral Ar in the butterfly, indicating that habitat loss has different effects on coevolutionary as compared with neutral processes. We conclude that this specialized coevolutionary system requires particular conservation interventions aiming at generating a spatial mosaic of both connected and of isolated habitat to maintain coevolutionary dynamics.

Facilitated Adaptation as A Conservation Tool in the Present Climate Change Context: A Methodological Guide

Climate change poses a novel threat to biodiversity that urgently requires the development of adequate conservation strategies. Living organisms respond to environmental change by migrating to locations where their ecological niche is preserved or by adapting to the new environment. While the first response has been used to develop, discuss and implement the strategy of assisted migration, facilitated adaptation is only beginning to be considered as a potential approach. Here, we present a review of the conceptual framework for facilitated adaptation, integrating advances and methodologies from different disciplines. Briefly, facilitated adaptation involves a population reinforcement that introduces beneficial alleles to enable the evolutionary adaptation of a focal population to pressing environmental conditions. To this purpose, we propose two methodological approaches. The first one (called pre-existing adaptation approach) is based on using pre-adapted genotypes existing in the focal population, in other populations, or even in closely related species. The second approach (called de novo adaptation approach) aims to generate new pre-adapted genotypes from the diversity present in the species through artificial selection. For each approach, we present a stage-by-stage procedure, with some techniques that can be used for its implementation. The associated risks and difficulties of each approach are also discussed.

Contrasting patterns of population structure of Bulwer's petrel (Bulweria bulwerii) between oceans revealed by statistical phylogeography

The patterns of population divergence of mid-latitude marine birds are impacted by only a few biogeographic barriers to dispersal and the effect of intrinsic factors, such as fidelity to natal colonies or wintering grounds, may become more conspicuous. Here we describe, for the first time, the phylogeographic patterns and historical demography of Bulwer's petrel Bulweria bulwerii and provide new insights regarding the drivers of species diversification in the marine environment. We sampled Bulwer's petrels from the main breeding colonies and used a statistical phylogeography approach based on surveying nuclear and mitochondrial loci (~ 9100 bp) to study its mechanisms of global diversification. We uncovered three highly differentiated groups including the Western Pacific, the Central Pacific and the Atlantic. The older divergence occurred within the Pacific Ocean, ca. 850,000 ya, and since then the W Pacific group has been evolving in isolation. Conversely, divergence between the Central Pacific and Atlantic populations occurred within the last 200,000 years. While the Isthmus of Panama is important in restricting gene flow between oceans in Bulwer's petrels, the deepest phylogeographic break is within the Pacific Ocean, where oceanographic barriers are key in driving and maintaining the remarkable structure found in this highly mobile seabird. This is in contrast with the Atlantic, where no structure was detected. Further data will provide insights regarding the extent of lineage divergence of Bulwer's petrels in the Western Pacific.

Dealing with disjunct populations of vascular plants: implications for assessing the effect of climate change

Species distribution models are the most widely used tool to predict species distributions for species conservation and assessment of climate change impact. However, they usually do not consider intraspecific ecological variation exhibited by many species. Overlooking the potential differentiation among groups of populations may lead to misplacing any conservation actions. This issue may be particularly relevant in species in which few populations with potential local adaptation occur, as in species with disjunct populations. Here, we used ecological niche modeling to analyze how the projections of current and future climatically suitable areas of 12 plant species can be affected using the whole taxa occurrences compared to occurrences from geographically disjunct populations. Niche analyses suggest that usually the disjunct group of populations selects the climatic conditions as similar as possible to the other according to climate availability. Integrating intraspecific variability only slightly increases models' ability to predict species occurrences. However, it results in different predictions of the magnitude of range change. In some species, integrating or not integrating intraspecific variability may lead to opposite trend in projected range change. Our results suggest that integrating intraspecific variability does not strongly improve overall models' accuracy, but it can result in considerably different conclusions about future range change. Consequently, accounting for intraspecific differentiation may enable the detection of potential local adaptations to new climate and so to design targeted conservation strategies.

Quantitative genetic-by-soil microbiome interactions in a perennial grass affect functional traits

Plants interact with diverse microbiomes that can impact plant growth and performance. Recent studies highlight the potential beneficial aspects of plant microbiomes, including the possibility that microbes facilitate the process of local adaptation in their host plants. Microbially mediated local adaptation in plants occurs when local host genotypes have higher fitness than foreign genotypes because of their affiliation with locally beneficial microbes. Here, plant adaptation results from genetic interactions of the host with locally beneficial microbes (e.g. host genotype-by-microbiome interactions). We used a recombinant inbred line (RIL) mapping population derived from upland and lowland ecotypes of the diploid C4 perennial bunch grass Panicum hallii to explore quantitative genetic responses to soil microbiomes focusing on functional root and shoot traits involved in ecotypic divergence. We show that the growth and development of ecotypes and their trait divergence depends on soil microbiomes. Moreover, we find that the genetic architecture is modified by soil microbiomes, revealing important plant genotype-by-microbiome interactions for quantitative traits. We detected a number of quantitative trait loci (QTL) that interact with the soil microbiome. Our results highlight the importance of microbial interactions in ecotypic divergence and trait genetic architecture in C4 perennial grasses.

Climate adaptation and vulnerability of foundation species in a global change hotspot

Climate change is altering species ranges, and relative abundances within ranges, as populations become differentially adapted and vulnerable to the climates they face. Understanding present species ranges, whether species harbour and exchange adaptive variants, and how variants are distributed across landscapes undergoing rapid change, is therefore crucial to predicting responses to future climates and informing conservation strategies. Such insights are nonetheless lacking for most species of conservation concern. We assess genomic patterns of neutral variation, climate adaptation and climate vulnerability (offsets in predicted distributions of putatively adaptive variants across present and future landscapes) for sister foundation species, the marine tubeworms Galeolaria caespitosa and Galeolaria gemineoa, in a sentinel region for climate change impacts. We find that species are genetically isolated despite uncovering sympatry in their ranges, show parallel and nonparallel signals of thermal adaptation on spatial scales smaller than gene flow across their ranges, and are predicted to face different risks of maladaptation under future temperatures across their ranges. Our findings have implications for understanding local adaptation in the face of gene flow, and generate spatially explicit predictions for climatic disruption of adaptation and species distributions in coastal ecosystems that could guide experimental validation and conservation planning.

The Biological Hierarchy, Time, and Temporal 'Omics in Evolutionary Biology: A Perspective

Sequencing data-genomics, transcriptomics, epigenomics, proteomics, and metabolomics-have revolutionized biological research, enabling a more detailed study of processes, ranging from subcellular to evolutionary, that drive biological organization. These processes, collectively, are responsible for generating patterns of phenotypic variation and can operate over dramatically different timescales (milliseconds to billions of years). While researchers often study phenotypic variation at specific levels of biological organization to isolate processes operating at that particular scale, the varying types of sequence data, or 'omics, can also provide complementary inferences to link molecular and phenotypic variation to produce an integrated view of evolutionary biology, ranging from molecular pathways to speciation. We briefly describe how 'omics has been used across biological levels and then demonstrate the utility of integrating different types of sequencing data across multiple biological levels within the same study to better understand biological phenomena. However, single-time-point studies cannot evaluate the temporal dynamics of these biological processes. Therefore, we put forward temporal 'omics as a framework that can better enable researchers to study the temporal dynamics of target processes. Temporal 'omics is not infallible, as the temporal sampling regime directly impacts inferential ability. Thus, we also discuss the role the temporal sampling regime plays in deriving inferences about the environmental conditions driving biological processes and provide examples that demonstrate the impact of the sampling regime on biological inference. Finally, we forecast the future of temporal 'omics by highlighting current methodological advancements that will enable temporal 'omics to be extended across species and timescales. We extend this discussion to using temporal multi-omics to integrate across the biological hierarchy to evaluate and link the temporal dynamics of processes that generate phenotypic variation.

Effects of Vertically Heterogeneous Soil Salinity on Genetic Polymorphism and Productivity of the Widespread Halophyte Bassia prostrata

Salinity is one of the environmental factors that affects both productivity and genetic diversity in plant species. Within the soil profile, salinity is a dynamic indicator and significantly changes with depth. The present study examined the effects of the vertical heterogeneity of soil salinity chemistry on the plant height, fresh and dry biomass accumulation, water content, level of genetic polymorphism, and observed and expected heterozygosity in seven populations of halophyte Bassia prostrata in natural habitats. Soil salinity ranged from slight (S_salts = 0.11-0.25%) to extreme (S_salts = 1.35-2.57%). The main contributors to salinity were Na⁺, Ca²⁺, and Mg²⁺. Multivariate analysis revealed that biomass accumulation is positively affected by moderate/high salinity in 20-60 cm soil layers, which may be associated with the salt required for the optimal growth of the halophyte B. prostrata. The formation of seed genetic diversity is negatively affected by slight/moderate salinity in the 0-40 cm layers. An increase in divalent ion content can reduce genetic diversity and increase the local adaptation of B. prostrata to magnesium-calcium sulfate salinity. The effect of the in-depth distribution of soil salinity on productivity and genetic diversity may be related to seasonal variables during biomass accumulation (summer) and seed formation (autumn).

The genomic basis of reproductive and migratory behaviour in a polymorphic salmonid

Recent ecotypic differentiation provides unique opportunities to investigate the genomic basis and architecture of local adaptation, while offering insights into how species form and persist. Sockeye salmon (Oncorhynchus nerka) exhibit migratory and resident ("kokanee") ecotypes, which are further distinguished into shore-spawning and stream-spawning reproductive ecotypes. Here, we analysed 36 sockeye (stream-spawning) and kokanee (stream- and shore-spawning) genomes from a system where they co-occur and have recent common ancestry (Okanagan Lake/River in British Columbia, Canada) to investigate the genomic basis of reproductive and migratory behaviour. Examination of the genomic landscape of differentiation, differences in allele frequencies and genotype-phenotype associations revealed three main blocks of sequence differentiation on chromosomes 7, 12 and 20, associated with migratory behaviour, spawning location and spawning timing. Structural variants identified in these same areas suggest they could contribute to ecotypic differentiation directly as causal variants or via maintenance of their genomic architecture through recombination suppression mechanisms. Genes in these regions were related to spatial memory and swimming endurance (SYNGAP, TPM3), as well as eye and brain development (including SIX6), potentially associated with differences in migratory behaviour and visual habitats across spawning locations, respectively. Additional genes (GREB1L, ROCK1) identified here have been associated with timing of migration in other salmonids and could explain variation in timing of O. nerka spawning. Together, these results based on the joint analysis of sequence and structural variation represent a significant advance in our understanding of the genomic landscape of ecotypic differentiation at different stages in the speciation continuum.

The role of genomic signatures of directional selection and demographic history in the population structure of a marine teleost with high gene flow

Recent studies have uncovered patterns of genomic divergence in marine teleosts where panmixia due to high gene flow has been the general paradigm. These signatures of divergent selection are often impacted by structural variants, acting as "supergenes" facilitating local adaptation. The highly dispersing European plaice (Pleuronectes platessa)-in which putative structural variants (i.e., inversions) have been identified-has successfully colonized the brackish water ecosystem of the Baltic Sea. Thus, the species represents an ideal opportunity to investigate how the interplay of gene flow, structural variants, natural selection, past demographic history, and gene flow impacts on population (sub)structuring in marine systems. Here, we report on the generation of an annotated draft plaice genome assembly in combination with population sequencing data-following the salinity gradient from the Baltic Sea into the North Sea together with samples from Icelandic waters-to illuminate genome-wide patterns of divergence. Neutral markers pointed at large-scale panmixia across the European continental shelf associated with high gene flow and a common postglacial colonization history of shelf populations. However, based on genome-wide outlier loci, we uncovered signatures of population substructuring among the European continental shelf populations, i.e., suggesting signs of ongoing selection. Genome-wide selection analyses (xp-EHH) and the identification of genes within genomic regions of recent selective sweeps-overlapping with the outlier loci-suggest that these represent the signs of divergent selection. Our findings provide support for genomic divergence driven by local adaptation in the face of high gene flow and elucidate the relative importance of demographic history versus adaptive divergence in shaping the contemporary population genetic structure of a marine teleost. The role of the putative inversion(s) in the substructuring-and potentially ongoing adaptation-was seemingly not substantial.

Fragmentation by major dams and implications for the future viability of platypus populations

The evolutionarily unique platypus (Ornithorhynchus anatinus) has experienced major declines and extinctions from a range of historical and recent interacting human-mediated threats. Although spending most of their time in the water, platypuses can move over land. Nevertheless, uncertainties remain whether dams are barriers to movement, thus limiting gene flow and dispersal, essential to evolution and ecology. Here we examined disruption of gene flow between platypus groups below and above five major dams, matched to four adjacent rivers without major dams. Genetic differentiation (F_ST) across dams was 4- to 20-fold higher than along similar stretches of adjacent undammed rivers; F_ST across dams was similar to differentiation between adjacent river systems. This indicates that major dams represent major barriers for platypus movements. Furthermore, F_ST between groups was correlated with the year in which the dam was built, increasing by 0.011 every generation, reflecting the effects of these barriers on platypus genetics. This study provides evidence of gene flow restriction, which jeopardises the long-term viability of platypus populations when groups are fragmented by major dams. Mitigation strategies, such as building of by-pass structures and translocation between upstream and downstream of the dam, should be considered in conservation and management planning.

Risks, benefits, and knowledge gaps of non-native tree species in Europe

Changing ecosystem conditions and diverse socio-economical events have contributed to an ingrained presence of non-native tree species (NNTs) in the natural and cultural European landscapes. Recent research endeavors have focused on different aspects of NNTs such as legislation, benefits, and risks for forestry, emphasizing that large knowledge gaps remain. As an attempt to fulfill part of these gaps, within the PEN-CAFoRR COST Action (CA19128) network, we established an open-access questionnaire that allows both academic experts and practitioners to provide information regarding NNTs from 20 European countries. Then, we integrated the data originating from the questionnaire, related to the country-based assessment of both peer-reviewed and grey literature, with information from available datasets (EUFORGEN and EU-Forest), which gave the main structure to the study and led to a mixed approach review. Finally, our study provided important insights into the current state of knowledge regarding NNTs. In particular, we highlighted NNTs that have shown to be less commonly addressed in research, raising caution about those characterized by an invasive behavior and used for specific purposes (e.g., wood production, soil recultivation, afforestation, and reforestation). NNTs were especially explored in the context of resilient and adaptive forest management. Moreover, we emphasized the assisted and natural northward migration of NNTs as another underscored pressing issue, which needs to be addressed by joint efforts, especially in the context of the hybridization potential. This study represents an additional effort toward the knowledge enhancement of the NNTs situation in Europe, aiming for a continuously active common source deriving from interprofessional collaboration.

Genetics reveals shifts in reproductive behaviour of the invasive bird parasite Philornis downsi collected from Darwin’s finch nests

Due to novel or dynamic fluctuations in environmental conditions and resources, host and parasite relationships can be subject to diverse selection pressures that may lead to significant changes during and after invasion of a parasite. Genomic analyses are useful for elucidating evolutionary processes in invasive parasites following their arrival to a new area and host. Philornis downsi (Diptera: Muscidae), the avian vampire fly, was introduced to the Galápagos Islands circa 1964 and has since spread across the archipelago, feeding on the blood of developing nestlings of endemic land birds. Since its discovery, there have been significant changes to the dynamics of P. downsi and its novel hosts, such as shifting mortality rates and changing oviposition behaviour, however no temporal genetic studies have been conducted. We collected P. downsi from nests and traps from a single island population over a 14-year period, and genotyped flies at 469 single nucleotide polymorphisms (SNPs) using restriction-site associated DNA sequencing (RADSeq). Despite significant genetic differentiation (FST) between years, there was no evidence for genetic clustering within or across four sampling years between 2006 and 2020, suggesting a lack of population isolation. Sibship reconstructions from P. downsi collected from 10 Darwin’s finch nests sampled in 2020 showed evidence for shifts in reproductive behaviour compared to a similar genetic analysis conducted in 2004–2006. Compared with this previous study, females mated with fewer males, individual females oviposited fewer offspring per nest, but more unique females oviposited per nest. These findings are important to consider within reproductive control techniques, and have fitness implications for both parasite evolution and host fitness.

Massive genome investigations reveal insights of prevalent introgression for environmental adaptation and triterpene biosynthesis in Ganoderma

Genome introgression is one of the driving forces that can increase species and genetic diversity and facilitate the adaptive evolution of organisms and biodiversity conservation. However, the genomic introgression and its contribution to biodiversity of macrofungi are still unclear. The genus Ganoderma is a typical macrofungal group that plays crucial roles in forest ecosystem as saprophytic organisms and plant pathogens, and is also involved in human health as medicinal mushrooms. Most public Ganoderma genomes are fragmented, and reference genomes and whole-genome information of diverse germplasm resources for many Ganoderma species are lacking, thus hindering functional and evolutionary genomic investigations among Ganoderma species. In this study, we provide high-quality genomes of 10 Ganoderma species and whole-genome variants data of 224 individuals from various ecoregions, enabling us to infer the phylogeny of Ganoderma species and their historical population dynamics. Based on whole-genome variants, widespread and genome-wide introgression among Ganoderma species is revealed. Genes with significant introgression signals were related to stress response, digestive absorption, and secondary metabolite synthesis, factors that may contribute to environmental adaptation and important biocomponent metabolism. CYP512U6, an essential functional gene in the CYP450 family related to Ganoderma triterpene synthesis, was detected with significant introgression and selection signals combined with Ganoderma metabolomic analysis, indicating that both ancient gene exchange and recent domestication have contributed to the categories and content of secondary metabolites of Ganoderma. The reference genomes, whole-genome variants, and metabolite profiles could serve as abundant and valuable genetic resources for evolution, ecology, and conservation investigations of Ganoderma species and other macrofungi.

Genetic architecture of repeated phenotypic divergence in Littorina saxatilis ecotype evolution

Chromosomal inversions have been shown to play a major role in a local adaptation by suppressing recombination between alternative arrangements and maintaining beneficial allele combinations. However, so far, their importance relative to the remaining genome remains largely unknown. Understanding the genetic architecture of adaptation requires better estimates of how loci of different effect sizes contribute to phenotypic variation. Here, we used three Swedish islands where the marine snail Littorina saxatilis has repeatedly evolved into two distinct ecotypes along a habitat transition. We estimated the contribution of inversion polymorphisms to phenotypic divergence while controlling for polygenic effects in the remaining genome using a quantitative genetics framework. We confirmed the importance of inversions but showed that contributions of loci outside inversions are of similar magnitude, with variable proportions dependent on the trait and the population. Some inversions showed consistent effects across all sites, whereas others exhibited site-specific effects, indicating that the genomic basis for replicated phenotypic divergence is only partly shared. The contributions of sexual dimorphism as well as environmental factors to phenotypic variation were significant but minor compared to inversions and polygenic background. Overall, this integrated approach provides insight into the multiple mechanisms contributing to parallel phenotypic divergence.

Selection in the city: Rapid and fine-scale evolution of urban eastern water dragons

Oceanic archipelagos have long been treated as a Petri dish for studies of evolutionary and ecological processes. Like archipelagos, cities exhibit similar patterns and processes, such as the rapid phenotypic divergence of a species between urban and nonurban environments. However, on a local scale, cities can be highly heterogenous, where geographically close populations can experience dramatically different environmental conditions. Nevertheless, we are yet to understand the evolutionary and ecological implications for populations spread across a heterogenous cityscape. To address this, we compared neutral genetic divergence to quantitative trait divergence within three native riparian and four city park populations of an iconic urban adapter, the eastern water dragon. We demonstrated that selection is likely acting to drive divergence of snout-vent length and jaw width across native riparian populations that are geographically isolated and across city park populations that are geographically close yet isolated by urbanization. City park populations as close as 0.9 km exhibited signs of selection-driven divergence to the same extent as native riparian populations isolated by up to 114.5 km. These findings suggest that local adaptation may be occurring over exceptionally small geographic and temporal scales within a single metropolis, demonstrating that city parks can act as archipelagos for the study of rapid evolution.

Genomic basis of insularity and ecological divergence in barn owls (Tyto alba) of the Canary Islands

Islands, and the particular organisms that populate them, have long fascinated biologists. Due to their isolation, islands offer unique opportunities to study the effect of neutral and adaptive mechanisms in determining genomic and phenotypical divergence. In the Canary Islands, an archipelago rich in endemics, the barn owl (Tyto alba), present in all the islands, is thought to have diverged into a subspecies (T. a. gracilirostris) on the eastern ones, Fuerteventura and Lanzarote. Taking advantage of 40 whole-genomes and modern population genomics tools, we provide the first look at the origin and genetic makeup of barn owls of this archipelago. We show that the Canaries hold diverse, long-standing and monophyletic populations with a neat distinction of gene pools from the different islands. Using a new method, less sensitive to structure than classical F_ST, to detect regions involved in local adaptation to insular environments, we identified a haplotype-like region likely under selection in all Canaries individuals and genes in this region suggest morphological adaptations to insularity. In the eastern islands, where the subspecies is present, genomic traces of selection pinpoint signs of adapted body proportions and blood pressure, consistent with the smaller size of this population living in a hot arid climate. In turn, genomic regions under selection in the western barn owls from Tenerife showed an enrichment in genes linked to hypoxia, a potential response to inhabiting a small island with a marked altitudinal gradient. Our results illustrate the interplay of neutral and adaptive forces in shaping divergence and early onset speciation.

Population Epigenetics: The Extent of DNA Methylation Variation in Wild Animal Populations

Population epigenetics explores the extent of epigenetic variation and its dynamics in natural populations encountering changing environmental conditions. In contrast to population genetics, the basic concepts of this field are still in their early stages, especially in animal populations. Epigenetic variation may play a crucial role in phenotypic plasticity and local adaptation as it can be affected by the environment, it is likely to have higher spontaneous mutation rate than nucleotide sequences do, and it may be inherited via non-mendelian processes. In this review, we aim to bring together natural animal population epigenetic studies to generate new insights into ecological epigenetics and its evolutionary implications. We first provide an overview of the extent of DNA methylation variation and its autonomy from genetic variation in wild animal population. Second, we discuss DNA methylation dynamics which create observed epigenetic population structures by including basic population genetics processes. Then, we highlight the relevance of DNA methylation variation as an evolutionary mechanism in the extended evolutionary synthesis. Finally, we suggest new research directions by highlighting gaps in the knowledge of the population epigenetics field. As for our results, DNA methylation diversity was found to reveal parameters that can be used to characterize natural animal populations. Some concepts of population genetics dynamics can be applied to explain the observed epigenetic structure in natural animal populations. The set of recent advancements in ecological epigenetics, especially in transgenerational epigenetic inheritance in wild animal population, might reshape the way ecologists generate predictive models of the capacity of organisms to adapt to changing environments.

Quantitative trait loci mapping reveals an oligogenic architecture of a rapidly adapting trait during the European invasion of common ragweed

Biological invasions offer a unique opportunity to investigate evolution over contemporary timescales. Rapid adaptation to local climates during range expansion can be a major determinant of invasion success, yet fundamental questions remain about its genetic basis. This study sought to investigate the genetic basis of climate adaptation in invasive common ragweed (Ambrosia artemisiifolia). Flowering time adaptation is key to this annual species' invasion success, so much so that it has evolved repeated latitudinal clines in size and phenology across its native and introduced ranges despite high gene flow among populations. Here, we produced a high-density linkage map (4493 SNPs) and paired this with phenotypic data from an F2 mapping population (n = 336) to identify one major and two minor quantitative trait loci (QTL) underlying flowering time and height differentiation in this species. Within each QTL interval, several candidate flowering time genes were also identified. Notably, the major flowering time QTL detected in this study was found to overlap with a previously identified haploblock (putative inversion). Multiple genetic maps of this region identified evidence of suppressed recombination in specific genotypes, consistent with inversions. These discoveries support the expectation that a concentrated genetic architecture with fewer, larger, and more tightly linked alleles should underlie rapid local adaptation during invasion, particularly when divergently adapting populations experience high levels of gene flow.

Assessment of Compound-Specific Fatty Acid δ13C and δ2H Values to Track Fish Mobility in a Small Sub-alpine Catchment

Methods for identifying origin, movement, and foraging areas of animals are essential for understanding ecosystem connectivity, nutrient flows, and other ecological processes. Telemetric methods can provide detailed spatial coverage but are limited to a minimum body size of specimen for tagging. In recent years, stable isotopes have been increasingly used to track animal migration by linking landscape isotope patterns into movement (isoscapes). However, compared to telemetric methods, the spatial resolution of bulk stable isotopes is low. Here, we examined a novel approach by evaluating the use of compound-specific hydrogen and carbon stable isotopes of fatty acids (δ²H_FA and δ¹³C_FA) from fish liver, muscle, brain, and eye tissues for identifying site specificity in a 254 km² sub-alpine river catchment. We analyzed 208 fish (European bullhead, rainbow trout, and brown trout) collected in 2016 and 2018 at 15 different sites. δ¹³C_FA values of these fish tissues correlated more among each other than those of δ²H_FA values. Both δ²H_FA and δ¹³C_FA values showed tissue-dependent isotopic fractionation, while fish taxa had only small effects. The highest site specificity was for δ¹³C_DHA values, while the δ²H isotopic difference between linoleic acid and alpha-linolenic acid resulted in the highest site specificity. Using linear discrimination analysis of FA isotope values, over 90% of fish could be assigned to their location of origin; however, the accuracy dropped to about 56% when isotope data from 2016 were used to predict the sites for samples collected in 2018, suggesting temporal shifts in site specificity of δ²H_FA and δ¹³C_FA. However, the predictive power of δ²H_FA and δ¹³C_FA over this time interval was still higher than site specificity of bulk tissue isotopes for a single time point. In summary, compound-specific isotope analysis of fatty acids may become a highly effective tool for assessing fine and large-scale movement and foraging areas of animals.

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

Across many species where inversions have been implicated in local adaptation, genomes often evolve to contain multiple, large inversions that arise early in divergence. Why this occurs has yet to be resolved. To address this gap, we built forward-time simulations in which inversions have flexible characteristics and can invade a metapopulation undergoing spatially divergent selection for a highly polygenic trait. In our simulations, inversions typically arose early in divergence, captured standing genetic variation upon mutation, and then accumulated many small-effect loci over time. Under special conditions, inversions could also arise late in adaptation and capture locally adapted alleles. Polygenic inversions behaved similarly to a single supergene of large effect and were detectable by genome scans. Our results show that characteristics of adaptive inversions found in empirical studies (e.g. multiple large, old inversions that are F_ST outliers, sometimes overlapping with other inversions) are consistent with a highly polygenic architecture, and inversions do not need to contain any large-effect genes to play an important role in local adaptation. By combining a population and quantitative genetic framework, our results give a deeper understanding of the specific conditions needed for inversions to be involved in adaptation when the genetic architecture is polygenic. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.