Perspectives on the use of landscape genetics to detect genetic adaptive variation in the field

Historical biogeography of the Mugil cephalus species complex and its rapid global colonization

Our understanding of speciation processes in marine environments remains very limited and the role of different reproductive barriers are still debated. While physical barriers were considered important drivers causing reproductive isolation, recent studies highlight the importance of climatic and hydrological changes creating unsuitable habitat conditions as factors promoting population isolation. Although speciation in marine fish has been investigated from different perspectives, these studies often have a limited geographical extant. Therefore, data on speciation within widely distributed species are largely lacking. Species complexes offer valuable opportunities to study the initial stages of speciation. Herein we study speciation within the Mugil cephalus species complex (MCSC) which presents a unique opportunity due to its circumglobal distribution. We used a whole-genome shotgun analysis approach to identify SNPs among the 16 species within the MCSC. We inferred the phylogenetic relationships within the species complex followed by a time-calibration analysis. Subsequently, we estimated the ancestral ranges within the species complex to explore their biogeographical history. Herein, we present a fully resolved and well-supported phylogeny of the MCSC. Its origin is dated at around 3.79 Ma after which two main clades emerged: one comprising all West Atlantic and East Pacific species and the other all East Atlantic and Indo-Pacific species. Rapid dispersal following an initial founder colonization from the West to the East Atlantic led to the population of all major realms worldwide in less than 2 Myr. Physical and climatic barriers heavily impacted the ancestral distribution ranges within the MCSC and triggered the onset of speciation.

South African indigenous chickens' genetic diversity, and the adoption of ecological niche modelling and landscape genomics as strategic conservation techniques

Selection pressures found in the prevailing production environments have shaped the genetic structure of indigenous chickens we see today. Indigenous chickens, raised in villages, provide essential genetic resources and income for poverty alleviation by providing affordable protein. However, they are threatened by predators, emerging diseases, and market demand for ideal breeds and fast production which causes loss of their valuable traits. The lack of knowledge about genetic diversity and genetic mechanisms underlying adaptive variants may compromise the goal of conserving indigenous chicken breeds. The main insights of the study are that indigenous chickens are highly diversified, and environmental factors play a key role in enabling chicken adaptation and distribution. Genomic and spatial technologies have made it possible to explore the genetic structure and fully comprehend the mechanism underlying the local adaptation of indigenous chickens. These technologies can aid in creating programs that enhance productivity and promote climate-resilient breeds. This review explores the impact of natural selection on indigenous chicken, genetic diversity, population size, and the advancement of technologies in understanding local adaptation drivers. In conclusion, this review highlights the importance of studying the habitats and how this will guide in conserving local breeds in their intended production environment.

Genomic Analysis of Sarda Sheep Raised at Diverse Temperatures Highlights Several Genes Involved in Adaptations to the Environment and Heat Stress Response

Livestock expresses complex traits influenced by several factors. The response of animals to variations in climatic factors, such as increases in temperature, may induce heat stress conditions. In this study, animals living at different temperatures were compared using the genome-wide Wright fixation index (FST). A total of 825 genotypes of Sarda breed ewes were divided into two groups based on the flocks' average temperature over a 20-year period to compute the FST: 395 and 430 sheep were represented in colder and hotter groups, respectively. After LOWESS regression and CONTROL CHART application, 623 significant markers and 97 selection signatures were found. A total of 280 positional candidate genes were retrieved from a public database. Among these genomic regions, we found 51 annotated genes previously associated with heat stress/tolerance in ruminants (FCGR1A, MDH1, UGP2, MYO1G, and HSPB3), as well as immune response and cellular mechanisms related to how animals cope with thermal stress (RIPK1, SERPINB1, SERPINB9, and PELI1). Moreover, other genes were associated with milk fat (SCD, HERC3, SCFD2, and CHUK), body weight, body fat, and intramuscular fat composition (AGPAT2, ABCD2, MFAP32, YTHDC1, SIRT3, SCD, and RNF121), which might suggest the influence of environmental conditions on the genome of Sarda sheep.

Adaptation of High-Altitude Plants to Harsh Environments: Application of Phenotypic-Variation-Related Methods and Multi-Omics Techniques

High-altitude plants face extreme environments such as low temperature, low oxygen, low nutrient levels, and strong ultraviolet radiation, causing them to adopt complex adaptation mechanisms. Phenotypic variation is the core manifestation of ecological adaptation and evolution. Many plants have developed a series of adaptive strategies through long-term natural selection and evolution, enabling them to survive and reproduce under such harsh conditions. This article reviews the techniques and methods used in recent years to study the adaptive evolution of high-altitude plants, including transplantation techniques, genomics, transcriptomics, proteomics, and metabolomics techniques, and their applications in high-altitude plant adaptive evolution. Transplantation technology focuses on phenotypic variation, which refers to natural variations in morphological, physiological, and biochemical characteristics, exploring their key roles in nutrient utilization, photosynthesis optimization, and stress-resistance protection. Multiple omics technologies, including genomics, transcriptomics, proteomics, and metabolomics, have revealed genes, regulatory pathways, and metabolic networks associated with phenotypic variations at the genetic and molecular levels. At the same time, the limitations and deficiencies of current technologies used to study plant adaptation to high-altitude environments were discussed. In addition, we propose future improvements to existing technologies and advocate for the integration of different technologies at multiple levels to study the molecular mechanisms of plant adaptation to high-altitude environments, thus providing insights for future research in this field.

Climatic factors, but not geographic distance, promote genetic structure and differentiation of Cleistogenes squarrosa (Trin.) Keng populations

Climate can shape plant genetic diversity and genetic structure, and genetic diversity and genetic structure can reflect the adaptation of plants to climate change. We used rbcl and trnL-trnF sequences to analyze the genetic diversity and genetic structure of C. squarrosa under the influence of different environmental factors in Inner Mongolia grassland. The results showed that the genetic diversity of this species was low. (The trnL-trnF sequences have higher genetic diversity than rbcl sequences.) C. squarrosa had low genetic diversity compared to other prairie plants, but had a more pronounced genetic structure. The haplotype network diagram of the combined sequences could be divided into two categories, and the results of the NJ, MP, and ML trees also showed that the haplotypes were divided into two branches. The results of genetic structure analysis showed that that the populations located in the desert steppe fall into exactly one cluster, and the populations located in the typical steppe fall into exactly another cluster. The neutrality tests were all negative and the mismatch distribution also showed a single peak across the population, suggesting that C. squarrosa had undergone population expansion and was well adapted to the local environment. The results of the mantel test showed that climate had a greater influence on the genetic distance of C. squarrosa, with annual precipitation having a higher influence than mean annual temperature. This study provided basic genetic information on the genetic structure of C. squarrosa and contributes to the study of genetic adaptation mechanisms in grassland plants.

Conservation genomics of an endangered montane amphibian reveals low population structure, low genomic diversity and selection pressure from disease

Wildlife diseases are a major global threat to biodiversity. Boreal toads (Anaxyrus [Bufo] boreas) are a state-endangered species in the southern Rocky Mountains of Colorado and New Mexico, and a species of concern in Wyoming, largely due to lethal skin infections caused by the amphibian chytrid fungus Batrachochytrium dendrobatidis (Bd). We performed conservation and landscape genomic analyses using single nucleotide polymorphisms from double-digest, restriction site-associated DNA sequencing in combination with the development of the first boreal toad (and first North American toad) reference genome to investigate population structure, genomic diversity, landscape connectivity and adaptive divergence. Genomic diversity (π = 0.00034-0.00040) and effective population sizes (Ne = 8.9-38.4) were low, likely due to post-Pleistocene founder effects and Bd-related population crashes over the last three decades. Population structure was also low, likely due to formerly high connectivity among a higher density of geographically proximate populations. Boreal toad gene flow was facilitated by low precipitation, cold minimum temperatures, less tree canopy, low heat load and less urbanization. We found >8X more putatively adaptive loci related to Bd intensity than to all other environmental factors combined, and evidence for genes under selection related to immune response, heart development and regulation and skin function. These data suggest boreal toads in habitats with Bd have experienced stronger selection pressure from disease than from other, broad-scale environmental variations. These findings can be used by managers to conserve and recover the species through actions including reintroduction and supplementation of populations that have declined due to Bd.

Genomic insights into the conservation of wild and domestic animal diversity: A review

Due to environmental change and anthropogenic activities, global biodiversity has suffered an unprecedented loss, and the world is now heading toward the sixth mass extinction event. This urges the need to step up our efforts to promote the sustainable use of animal genetic resources and plan effective strategies for their conservation. Although habitat preservation and restoration are the primary means of conserving biodiversity, genomic technologies offer a variety of novel tools for identifying biodiversity hotspots and thus, support conservation efforts. Conservation genomics is a broad area of science that encompasses the application of genomic data from thousands or tens of thousands of genome-wide markers to address important conservation biology concerns. Genomic approaches have revolutionized the way we understand and manage animal populations, providing tools to identify and preserve unique genetic variants and alleles responsible for adaptive genetic variation, reducing the deleterious consequences of inbreeding, and increasing the adaptive potential of threatened species. The advancement of genomic technologies, particularly comparative genomic approaches, and the increased accessibility of genomic resources in the form of genome-enabled taxa for non-model organisms, provides a distinct advantage in defining conservation units over traditional genetics approaches. The objective of this review is to provide an exhaustive overview of the concept of conservation genomics, discuss the rationale behind the transition from conservation genetics to genomic approaches, and emphasize the potential applications of genomic techniques for conservation purposes. We also highlight interesting case studies in both livestock and wildlife species where genomic techniques have been used to accomplish conservation goals. Finally, we address some challenges and future perspectives in this field.

Individual-based landscape genomics for conservation: An analysis pipeline

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

Urbanization drives adaptive evolution in a Neotropical bird

Urbanization has dramatic impacts on natural habitats and such changes may potentially drive local adaptation of urban populations. Behavioral change has been specifically shown to facilitate the fast adaptation of birds to changing environments, but few studies have investigated the genetic mechanisms of this process. Such investigations could provide insights into questions about both evolutionary theory and management of urban populations. In this study, we investigated whether local adaptation has occurred in urban populations of a Neotropical bird species, Coereba flaveola, specifically addressing whether observed behavioral adaptations are correlated to genetic signatures of natural selection. To answer this question, we sampled 24 individuals in urban and rural environments, and searched for selected loci through a genome-scan approach based on RADseq genomic data, generated and assembled using a reference genome for the species. We recovered 46 loci as putative selection outliers, and 30 of them were identified as associated with biological processes possibly related to urban adaptation, such as the regulation of energetic metabolism, regulation of genetic expression, and changes in the immunological system. Moreover, genes involved in the development of the nervous system showed signatures of selection, suggesting a link between behavioral and genetic adaptations. Our findings, in conjunction with similar results in previous studies, support the idea that cities provide a similar selective pressure on urban populations and that behavioral plasticity may be enhanced through genetic changes in urban populations.

Tropical forest loss and geographic location drive the functional genomic diversity of an endangered palm tree

Human activity has diminished forests in different terrestrial ecosystems. This is well illustrated in the Brazilian Atlantic Forest, which still hosts high levels of species richness and endemism, even with only 28% of its original extent remaining. The consequences of such forest loss in remaining populations can be investigated with several approaches, including the genomic perspective, which allows a broader understanding of how human disturbance influences the genetic variability in natural populations. In this context, our study investigated the genomic responses of Euterpe edulis Martius, an endangered palm tree, in forest remnants located in landscapes presenting different forest cover amount and composed by distinct bird assemblage that disperse its seeds. We sampled 22 areas of the Brazilian Atlantic Forest in four regions using SNP markers inserted into transcribed regions of the genome of E. edulis, distinguishing neutral loci from those putatively under natural selection (outlier). We demonstrate that populations show patterns of structure and genetic variability that differ between regions, as a possible reflection of deforestation and biogeographic histories. Deforested landscapes still maintain high neutral genetic diversity due to gene flow over short distances. Overall, we not only support previous evidence with microsatellite markers, but also show that deforestation can influence the genetic variability outlier, in the scenario of selective pressures imposed by these stressful environments. Based on our findings, we suggest that, to protect genetic diversity in the long term, it is necessary to reforest and enrich deforested areas, using seeds from populations in the same management target region.

Landscape Genetics for Brazilian Equines

Optimization of DNA collection for National gene bank and conservation programs requires information on spatial and genetic distribution of animals countrywide. The relationship between genetic and geographic distances were examined in 8 Brazilian horse breeds (Baixadeiro, Crioulo, Campeiro, Lavradeiro, Marajoara, Mangalarga Marchador, Pantaneiro and Puruca) using Single Nucleotide Polymorphism markers and collection point locations. Mantel correlations, Genetic Landscape Shape Interpolation, Allelic Aggregation Index Analyses and Spatial autocorrelation tests indicated a nonrandom distribution of horses throughout the country. Minimum collection distances for the national Gene Bank should be 530km, with clear divisions seen in genetic structure of horse populations in both North/South and East/West directions. Comparing Pantaneiro and North/Northeastern breeds, physical distance is not necessarily the defining factor for genetic differentiation. This should be considered when sampling these local breeds. These data can help optimise GenBank collection routines and conservation strategies for these breeds.

Deciphering local adaptation of native Indian cattle (Bos indicus) breeds using landscape genomics and in-silico prediction of deleterious SNP effects on protein structure and function

India has 50 registered breeds of native cattle (Bos indicus) which are locally adapted to diverse environmental conditions. This study aimed to investigate the genomic basis of adaptation of native Indian cattle and to predict the impact of key SNPs on the amino acid changes that affect protein function. The Illumina 777 K BovineHD BeadChip was used to genotype 178 native cattle belonging to contrasting landscapes and agro-climatic conditions. The genotype-environment association was investigated with R. SamBada, using 5,74,382 QC passed SNPs and 11 predictor variables (10 multi-collinearity controlled environmental variables and 1 variable as "score of PCA" on ancestry coefficients of individuals). In total, 1,12,780 models were selected as significant (q < 0.05) based on G score. The pathway ontology of the annotated genes revealed many important pathways and genes having a direct and indirect role in cold and hot adaptation. Only ten SNP variants had a SIFT score of < 0.05 (deleterious), and only two of them, each lying in the genes CRYBA1 and USP18, were predicted to be deleterious with high confidence. RaptorX predicted the tertiary structures of proteins encoded by wild and mutant variants of these genes. The quality of the models was determined using Ramachandran plots and RaptorX parameters, indicating that they are accurate. RaptorX and I-Mutant 2.0 softwares revealed significant differences among wild and mutant proteins. Adaptive alleles identified in the present investigation might be responsible for the local adaptation of these cattle breeds.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03493-3.

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.

Seascape genomics of common dolphins (Delphinus delphis) reveals adaptive diversity linked to regional and local oceanography

Background

High levels of standing genomic variation in wide-ranging marine species may enhance prospects for their long-term persistence. Patterns of connectivity and adaptation in such species are often thought to be influenced by spatial factors, environmental heterogeneity, and oceanographic and geomorphological features. Population-level studies that analytically integrate genome-wide data with environmental information (i.e., seascape genomics) have the potential to inform the spatial distribution of adaptive diversity in wide-ranging marine species, such as many marine mammals. We assessed genotype-environment associations (GEAs) in 214 common dolphins (Delphinus delphis) along > 3000 km of the southern coast of Australia.

Results

We identified 747 candidate adaptive SNPs out of a filtered panel of 17,327 SNPs, and five putatively locally-adapted populations with high levels of standing genomic variation were disclosed along environmentally heterogeneous coasts. Current velocity, sea surface temperature, salinity, and primary productivity were the key environmental variables associated with genomic variation. These environmental variables are in turn related to three main oceanographic phenomena that are likely affecting the dispersal of common dolphins: (1) regional oceanographic circulation, (2) localised and seasonal upwellings, and (3) seasonal on-shelf circulation in protected coastal habitats. Signals of selection at exonic gene regions suggest that adaptive divergence is related to important metabolic traits.

Conclusion

To the best of our knowledge, this represents the first seascape genomics study for common dolphins (genus Delphinus). Information from the associations between populations and their environment can assist population management in forecasting the adaptive capacity of common dolphins to climate change and other anthropogenic impacts.

The Past, Present, and Future of Host Plant Resistance in Cotton: An Australian Perspective

Cotton is a key global fiber crop. However, yield potential is limited by the presence of endemic and introduced pests and diseases. The introduction of host plant resistance (HPR), defined as the purposeful use of resistant crop cultivars to reduce the impact of pests and diseases, has been a key breeding target for the Commonwealth Scientific and Industrial Research Organisation (CSIRO) cotton breeding program. The program has seen success in releasing cultivars resistant to Bacterial blight, Verticillium wilt, Fusarium wilt, and Cotton bunchy top. However, emerging biotic threats such as Black root rot and secondary pests, are becoming more frequent in Australian cotton production systems. The uptake of tools and breeding methods, such as genomic selection, high throughput phenomics, gene editing, and landscape genomics, paired with the continued utilization of sources of resistance from Gossypium germplasm, will be critical for the future of cotton breeding. This review celebrates the success of HPR breeding activities in the CSIRO cotton breeding program and maps a pathway for the future in developing resistant cultivars.

Population and Landscape Genetics Provide Insights Into Species Conservation of Two Evergreen Oaks in Qinghai-Tibet Plateau and Adjacent Regions

The combination of population and landscape genetics can facilitate the understanding of conservation strategy under the changing climate. Here, we focused on the two most diverse and ecologically important evergreen oaks: Quercus aquifolioides and Quercus spinosa in Qinghai-Tibetan Plateau (QTP), which is considered as world's biodiversity hotspot. We genotyped 1,657 individuals of 106 populations at 15 nuclear microsatellite loci throughout the species distribution range. Spatial patterns of genetic diversity were identified by mapping the allelic richness (AR) and locally common alleles (LCA) according to the circular neighborhood methodology. Migration routes from QTP were detected by historical gene flow estimation. The response pattern of genetic variation to environmental gradient was assessed by the genotype-environment association (GEA) analysis. The overall genetic structure showed a high level of intra-species genetic divergence of a strong west-east pattern. The West-to-East migration route indicated the complex demographic history of two oak species. We found evidence of isolation by the environment in Q. aqu-East and Q. spi-West lineage but not in Q. aqu-West and Q. spi-East lineage. Furthermore, priority for conservation should be given to populations that retain higher spatial genetic diversity or isolated at the edge of the distribution range. Our findings indicate that knowledge of spatial diversity and migration route can provide valuable information for the conservation of existing populations. This study provides an important guide for species conservation for two oak species by the integration of population and landscape genetic methods.

Traces of Genetic but Not Epigenetic Adaptation in the Invasive Goldenrod Solidago canadensis Despite the Absence of Population Structure

Biological invasions may result from multiple introductions, which might compensate for reduced gene pools caused by bottleneck events, but could also dilute adaptive processes. A previous common-garden experiment showed heritable latitudinal clines in fitness-related traits in the invasive goldenrod Solidago canadensis in Central Europe. These latitudinal clines remained stable even in plants chemically treated with zebularine to reduce epigenetic variation. However, despite the heritability of traits investigated, genetic isolation-by-distance was non-significant. Utilizing the same specimens, we applied a molecular analysis of (epi)genetic differentiation with standard and methylation-sensitive (MSAP) AFLPs. We tested whether this variation was spatially structured among populations and whether zebularine had altered epigenetic variation. Additionally, we used genome scans to mine for putative outlier loci susceptible to selection processes in the invaded range. Despite the absence of isolation-by-distance, we found spatial genetic neighborhoods among populations and two AFLP clusters differentiating northern and southern Solidago populations. Genetic and epigenetic diversity were significantly correlated, but not linked to phenotypic variation. Hence, no spatial epigenetic patterns were detected along the latitudinal gradient sampled. Applying genome-scan approaches (BAYESCAN, BAYESCENV, RDA, and LFMM), we found 51 genetic and epigenetic loci putatively responding to selection. One of these genetic loci was significantly more frequent in populations at the northern range. Also, one epigenetic locus was more frequent in populations in the southern range, but this pattern was lost under zebularine treatment. Our results point to some genetic, but not epigenetic adaptation processes along a large-scale latitudinal gradient of S. canadensis in its invasive range.

Habitat-linked genetic variation supports microgeographic adaptive divergence in an island-endemic bird species

We investigated the potential mechanisms driving habitat-linked genetic divergence within a bird species endemic to a single 250 km² island. The island scrub-jay (Aphelocoma insularis) exhibits microgeographic divergence in bill morphology across pine-oak ecotones on Santa Cruz Island, California (USA) similar to adaptive differences described in mainland congeners over much larger geographic scales. To test whether individuals exhibit genetic differentiation related to habitat type and divergence in bill length, we genotyped over 3,000 single nucleotide polymorphisms (SNPs) in 123 adult island scrub-jay males from across Santa Cruz Island using restriction site-associated DNA sequencing (RADseq). Neutral landscape genomic analyses revealed that genome-wide genetic differentiation was primarily related to geographic distance and differences in habitat composition. We also found 168 putatively adaptive loci associated with habitat type using multivariate redundancy analysis (RDA) while controlling for spatial effects. Finally, two genome-wide association analyses revealed a polygenic basis to variation in bill length with multiple loci detected in or near genes known to affect bill morphology in other birds. Our findings support the hypothesis that divergent selection at microgeographic scales can cause adaptive divergence in the presence of ongoing gene flow.

Applying genomics in assisted migration under climate change: Framework, empirical applications, and case studies

The rate of global climate change is projected to outpace the ability of many natural populations and species to adapt. Assisted migration (AM), which is defined as the managed movement of climate-adapted individuals within or outside the species ranges, is a conservation option to improve species' adaptive capacity and facilitate persistence. Although conservation biologists have long been using genetic tools to increase or maintain diversity of natural populations, genomic techniques could add extra benefit in AM that include selectively neutral and adaptive regions of the genome. In this review, we first propose a framework along with detailed procedures to aid collaboration among scientists, agencies, and local and regional managers during the decision-making process of genomics-guided AM. We then summarize the genomic approaches for applying AM, followed by a literature search of existing incorporation of genomics in AM across taxa. Our literature search initially identified 729 publications, but after filtering returned only 50 empirical studies that were either directly applied or considered genomics in AM related to climate change across taxa of plants, terrestrial animals, and aquatic animals; 42 studies were in plants. This demonstrated limited application of genomic methods in AM in organisms other than plants, so we provide further case studies as two examples to demonstrate the negative impact of climate change on non-model species and how genomics could be applied in AM. With the rapidly developing sequencing technology and accumulating genomic data, we expect to see more successful applications of genomics in AM, and more broadly, in the conservation of biodiversity.

Land use and life history constrain adaptive genetic variation and reduce the capacity for climate change adaptation in turtles

BACKGROUND

Rapid anthropogenic climate change will require species to adapt to shifting environmental conditions, with successful adaptation dependent upon current patterns of genetic variation. While landscape genomic approaches allow for exploration of local adaptation in non-model systems, most landscape genomics studies of adaptive capacity are limited to exploratory identification of potentially important functional genes, often without a priori expectations as to the gene functions that may be most important for climate change responses. In this study, we integrated targeted sequencing of genes of known function and genotyping of single-nucleotide polymorphisms to examine spatial, environmental, and species-specific patterns of potential local adaptation in two co-occuring turtle species: the Blanding's turtle (Emydoidea blandingii) and the snapping turtle (Chelydra serpentina).

RESULTS

We documented divergent patterns of spatial clustering between neutral and putatively adaptive genetic variation in both species. Environmental associations varied among gene regions and between species, with stronger environmental associations detected for genes involved in stress response and for the more specialized Blanding's turtle. Land cover appeared to be more important than climate in shaping spatial variation in functional genes, indicating that human landscape alterations may affect adaptive capacity important for climate change responses.

CONCLUSIONS

Our study provides evidence that responses to climate change will be contingent on species-specific adaptive capacity and past history of exposure to human land cover change.

Climate differently influences the genomic patterns of two sympatric marine fish species

Climate influences population genetic variation in marine species. Capturing these impacts remains challenging for marine fishes which disperse over large geographical scales spanning steep environmental gradients. It requires the extensive spatial sampling of individuals or populations, representative of seascape heterogeneity, combined with a set of highly informative molecular markers capable of revealing climatic-associated genetic variations. We explored how space, dispersal and environment shape the genomic patterns of two sympatric fish species in the Mediterranean Sea, which ranks among the oceanic basins most affected by climate change and human pressure. We hypothesized that the population structure and climate-associated genomic signatures of selection would be stronger in the less mobile species, as restricted gene flow tends to facilitate the fixation of locally adapted alleles. To test our hypothesis, we genotyped two species with contrasting dispersal abilities: the white seabream Diplodus sargus and the striped red mullet Mullus surmuletus. We collected 823 individuals and used genotyping by sequencing (GBS) to detect 8,206 single nucleotide polymorphisms (SNPs) for the seabream and 2,794 for the mullet. For each species, we identified highly differentiated genomic regions (i.e. outliers) and disentangled the relative contribution of space, dispersal and environmental variables (climate, marine primary productivity) on the outliers' genetic structure to test the prevalence of gene flow and local adaptation. We observed contrasting patterns of gene flow and adaptive genetic variation between the two species. The seabream showed a distinct Alboran sea population and panmixia across the Mediterranean Sea. The mullet revealed additional differentiation within the Mediterranean Sea that was significantly correlated to summer and winter temperatures, as well as marine primary productivity. Functional annotation of the climate-associated outlier SNPs then identified candidate genes involved in heat tolerance that could be examined to further predict species' responses to climate change. Our results illustrate the key steps of a comparative seascape genomics study aiming to unravel the evolutionary processes at play in marine species, to better anticipate their response to climate change. Defining population adaptation capacities and environmental niches can then serve to incorporate evolutionary processes into species conservation planning.

Assessing signals of selection and historical demography to develop conservation strategies in the Chilean emblematic Araucaria araucana

Loss of genetic diversity reduces the ability of species to evolve and respond to environmental change. Araucaria araucana is an emblematic conifer species from southern South America, with important ethnic value for the Mapuche people (Pehuenche); the Chilean Government has catalogued its conservation status as vulnerable. Climatic fluctuations were potentially a major impact in the genetic variation within many tree species. In this context, the restricted geographic distribution of A. araucana in Chile appears to be a consequence of the Last Glacial Maximum (LGM). During the past two centuries, strong human intervention has also affected the geographical distribution and population sizes of A. araucana. Reduction of population size may cause loss of genetic diversity, which could affect frequency of adaptive loci. The aims of this study were to know the existence of potential loci under selection and populations with genetic, demographic disequilibrium in the Chilean distribution of A. araucana. Based on 268 polymorphic AFLP loci, we have investigated potential loci under selection and genetic, demographic disequilibrium within seven Chilean populations of Araucaria araucana. Correlation of 41 outlier loci with the environmental variables of precipitation and temperature reveals signatures of selection, whereas 227 neutral loci provide estimates of demographic equilibrium and genetic population structure. Three populations are recommended as priorities for conservation.

Prospects for the natural distribution of crop wild-relatives with limited adaptability: The case of the wild pea Pisum fulvum

Plant breeders and conservationist depend on knowledge about the genetic variation of their species of interest. Pisum fulvum, a wild relative of domesticated pea, has attracted attention as a genetic resource for crop improvement, yet little information about its diversity in the wild has been published hitherto. We sampled 15 populations of P. fulvum from Israeli natural habitats and conducted genotyping by sequencing to analyse their genetic diversity and adaptive state. We also attempted to evaluate the species past demography and the prospects of its future reaction to environmental changes. The results suggest that genetic diversity of P. fulvum is low to medium and is distributed between well diverged populations. Surprisingly, with 56 % in the total population the selfing rate was found to be significantly lower than expected from a species that is commonly assumed to be a predominant selfer. We found a strong genetic bottleneck during the last glacial period and only limited patterns of isolation by distance and environment, which explained 13 %-18 % of the genetic variation. Despite the weak signatures of genome-wide IBE, 1,354 markers were significantly correlated with environmental factors, 1,233 of which were located within known genes with a nonsynonymous to synonymous ratio of 0.382. Species distribution modelling depicted an ongoing fragmentation and decreased habitable area over the next 80 years under two different socio-economic pathways. Our results suggest that complex interactions of substantial drift and selection shaped the genome of P. fulvum. Climate changeis likely to cause further erosion of genetic diversity in P. fulvum. Systematic ex-situ conservation may be advisable to safeguard genetic variability for future utilization of this species.

Adaptive divergence despite low effective population size in a peripherally isolated population of the pygmy rabbit, Brachylagus idahoensis

Local adaptation can occur when spatially separated populations are subjected to contrasting environmental conditions. Historically, understanding the genetic basis of adaptation has been difficult, but increased availability of genome-wide markers facilitates studies of local adaptation in non-model organisms of conservation concern. The pygmy rabbit (Brachylagus idahoensis) is an imperiled lagomorph that relies on sagebrush for forage and cover. This reliance has led to widespread population declines following reductions in the distribution of sagebrush, leading to geographic separation between populations. In this study, we used >20,000 single nucleotide polymorphisms, genotype-environment association methods, and demographic modeling to examine neutral genetic variation and local adaptation in the pygmy rabbit in Nevada and California. We identified 308 loci as outliers, many of which had functional annotations related to metabolism of plant secondary compounds. Likewise, patterns of spatial variation in outlier loci were correlated with landscape and climatic variables including proximity to streams, sagebrush cover, and precipitation. We found that populations in the Mono Basin of California probably diverged from other Great Basin populations during late Pleistocene climate oscillations, and that this region is adaptively differentiated from other regions in the southern Great Basin despite limited gene flow and low effective population size. Our results demonstrate that peripherally isolated populations can maintain adaptive divergence.

Multi-Approach Analysis Reveals Local Adaptation in a Widespread Forest Tree of Reunion Island

Detecting processes of local adaptation in forest trees and identifying environmental selective drivers are of primary importance for forest management and conservation. Transplant experiments, functional genomics and population genomics are complementary tools to efficiently characterize heritable phenotypic traits and to decipher the genetic bases of adaptive traits. Using an integrative approach combining phenotypic assessment in common garden, transcriptomics and landscape genomics, we investigated leaf adaptive traits in Coffea mauritiana, a forest tree endemic to Reunion Island. Eight populations of C. mauritiana originating from sites with contrasted environmental conditions were sampled in common garden to assess several leaf morphological traits, to analyze the leaf transcriptome and leaf cuticular wax composition. The relative alkane content of cuticular waxes was significantly correlated with major climatic gradients, paving the way for further transcriptome-based analyses. The expression pattern of cuticle biosynthetic genes was consistent with a modulation of alkane accumulation across the population studied, supporting the hypothesis that the composition of cuticular wax is involved in the local adaptation of C. mauritiana. Association tests in landscape genomics performed using RNA-seq-derived single-nucleotide polymorphisms revealed that genes associated with cell wall remodeling also likely play an adaptive role. By combining these different approaches, this study efficiently identified local adaptation processes in a non-model species. Our results provide the first evidence for local adaptation in trees endemic to Reunion Island and highlight the importance of cuticle composition for the adaptation of trees to the high evaporative demand in warm climates.

Anthropogenic load іs a leading factor in the morphological variability of Chondrula tridens (Gastropoda, Enidae) in the northwestern Azov Sea region

Morphometric data are widely used in biology to assess intraspecific and inter-population variability and for bioindication and environmental condition assessment. The following hypotheses have been experimentally tested in the paper: 1) the vegetation type affects the change in the shell shape of Chondrula tridens martynovi Gural-Sverlova & Gural, 2010; 2) the change in the shell shape of this species is influenced by the biotope moisture regime; 3) the shell shape changes depending on the anthropogenic load level. The material in the form of empty, fully formed Ch. tridens shells was collected in 2019 in the north-western Azov region within the basin of the Molochna River. The collection points were located in settlements and outside them and differed in vegetation, moisture regime and level of anthropogenic load. The vegetation has been expertly attributed to two alternative types: herbaceous vegetation and tree plantations. By moisture level, the locations have been assessed as xerophytic and mesoxerophytic. The anthropogenic load levels have been assessed as low, medium and high. The study revealed that the morphological characteristics of Ch. tridens demonstrate a significant component of variability, which is due to the shell size. The shell size depends on the anthropogenic impact level. Under conditions of high anthropogenic impact, the shell size increases. Mollusks from locations with low and medium anthropogenic impact levels did not differ in shell size. After extraction of the size component, morphological properties develop three main trends of variability. The mouth apparatus development of mollusks does not depend on the vegetation type, but depends on the biotope moisture level and the anthropogenic transformation level. The mollusk shell elongation was observed to have the opposite dynamics of the height parameters in relation to the width and depended on the level of anthropogenic load. Rearrangement in the mouth apparatus depended on the biotope moisture level and the anthropogenic load level. There were distinguished four clusters, the quantitative morphological features of which allowed us to identify them as morphotypes. Each location was characterized by a combination of different morphotypes, according to which the sampling points may be classified. Morphotype 1 corresponds to biotopes with low level of anthropogenic load, morphotype 4 corresponded to biotopes with high anthropogenic load. Morphotypes 2 and 3 corresponded to moderate level of anthropogenic load. Vegetation type is not an important factor in determining the morphotypic diversity of populations. Under xerophytic conditions, morphotypes 2 and 3 are more common, and under mesoxerophytic conditions, morphotypes 1 and 4 are more common. The range of molluscs in different habitats needs to be expanded in the future to clarify climatic and other patterns.

Comprehensive genotyping of a Brazilian cassava (Manihot esculenta Crantz) germplasm bank: insights into diversification and domestication

KEY MESSAGE

Brazilian cassava diversity was characterized through population genetics and clustering approaches, highlighting contrasted genetic groups and spatial genetic differentiation. Cassava (Manihot esculenta Crantz) is a major staple root crop of the tropics, originating from the Amazonian region. In this study, 3354 cassava landraces and modern breeding lines from the Embrapa Cassava Germplasm Bank (CGB) were characterized. All individuals were subjected to genotyping-by-sequencing (GBS), identifying 27,045 single-nucleotide polymorphisms (SNPs). Identity-by-state and population structure analyses revealed a unique set of 1536 individuals and 10 distinct genetic groups with heterogeneous linkage disequilibrium (LD). On this basis, a density of 1300-4700 SNP markers were selected for large-effect quantitative trait loci (QTL) detection. Identified genetic groups were further characterized for population genetics parameters including minor allele frequency (MAF), observed heterozygosity [Formula: see text], effective population size estimate [Formula: see text]) and polymorphism information content (PIC). Selection footprints and introgressions of M. glaziovii were detected. Spatial population structure analysis revealed five ancestral populations related to distinct Brazilian ecoregions. Estimation of historical relationships among identified populations suggests an early population split from Amazonian to Atlantic forest and Caatinga ecoregions and active gene flows. This study provides a thorough genetic characterization of ex situ germplasm resources from cassava's center of origin, South America, with results shedding light on Brazilian cassava characteristics and its biogeographical landscape. These findings support and facilitate the use of genetic resources in modern breeding programs including implementation of association mapping and genomic selection strategies.

Comprehensive genotyping of a Brazilian cassava (Manihot esculenta Crantz) germplasm bank: insights into diversification and domestication

Brazilian cassava diversity was characterized through population genetics and clustering approaches, highlighting contrasted genetic groups and spatial genetic differentiation. Cassava (Manihot esculenta Crantz) is a major staple root crop of the tropics, originating from the Amazonian region. In this study, 3354 cassava landraces and modern breeding lines from the Embrapa Cassava Germplasm Bank (CGB) were characterized. All individuals were subjected to genotyping-by-sequencing (GBS), identifying 27,045 single-nucleotide polymorphisms (SNPs). Identity-by-state and population structure analyses revealed a unique set of 1536 individuals and 10 distinct genetic groups with heterogeneous linkage disequilibrium (LD). On this basis, a density of 1300-4700 SNP markers were selected for large-effect quantitative trait loci (QTL) detection. Identified genetic groups were further characterized for population genetics parameters including minor allele frequency (MAF), observed heterozygosity [Formula: see text], effective population size estimate [Formula: see text]) and polymorphism information content (PIC). Selection footprints and introgressions of M. glaziovii were detected. Spatial population structure analysis revealed five ancestral populations related to distinct Brazilian ecoregions. Estimation of historical relationships among identified populations suggests an early population split from Amazonian to Atlantic forest and Caatinga ecoregions and active gene flows. This study provides a thorough genetic characterization of ex situ germplasm resources from cassava's center of origin, South America, with results shedding light on Brazilian cassava characteristics and its biogeographical landscape. These findings support and facilitate the use of genetic resources in modern breeding programs including implementation of association mapping and genomic selection strategies.

Can We Compare Effect Size of Spatial Genetic Structure Between Studies and Species Using Moran Eigenvector Maps?

As the field of landscape genetics is progressing toward comparative empirical studies and meta-analysis, it is important to know how best to compare the strength of spatial genetic structure between studies and species. Moran’s Eigenvector Maps are a promising method that does not make an assumption of isolation-by-distance in a homogeneous environment but can discern cryptic structure that may result from multiple processes operating in heterogeneous landscapes. MEMgene uses spatial filters from Moran’s Eigenvector Maps as predictor variables to explain variation in a genetic distance matrix, and it returns adjusted R2 as a measure of the amount of genetic variation that is spatially structured. However, it is unclear whether, and under which conditions, this value can be used to compare the degree of spatial genetic structure (effect size) between studies. This study addresses the fundamental question of comparability at two levels: between independent studies (meta-analysis mode) and between species sampled at the same locations (comparative mode). We used published datasets containing 9,900 haploid, biallelic, neutral loci simulated on a quasi-continuous, square landscape under four demographic scenarios (island model, isolation-by-distance, expansion from one or two refugia). We varied the genetic resolution (number of individuals and loci) and the number of random sampling locations. We considered two measures of effect size, the MEMgene adjusted R2 and multivariate Moran’s I, which is related to Moran’s Eigenvector Maps. Both metrics were highly sensitive to the number of locations, even when using standardized effect sizes, SES, and the number of individuals sampled per location, but not to the number of loci. In comparative mode, using the same Moran Eigenvector Maps for all species, even those with missing values at some sampling locations, reduced bias due to the number of locations under isolation-by-distance (stationary process) but increased it under expansion from one or two refugia (non-stationary process). More robust measures of effect size need to be developed before the strength of spatial genetic structure can be accurately compared, either in a meta-analysis of independent empirical studies or within a comparative, multispecies landscape genetic study.

Genetic adaptation of Tibetan poplar (Populus szechuanica var. tibetica) to high altitudes on the Qinghai-Tibetan Plateau

Plant adaptation to high altitudes has long been a substantial focus of ecological and evolutionary research. However, the genetic mechanisms underlying such adaptation remain poorly understood. Here, we address this issue by sampling, genotyping, and comparing populations of Tibetan poplar, Populus szechuanica var. tibetica, distributed from low (~2,000 m) to high altitudes (~3,000 m) of Sejila Mountain on the Qinghai-Tibet Plateau. Population structure analyses allow clear classification of two groups according to their altitudinal distributions. However, in contrast to the genetic variation within each population, differences between the two populations only explain a small portion of the total genetic variation (3.64%). We identified asymmetrical gene flow from high- to low-altitude populations. Integrating population genomic and landscape genomic analyses, we detected two hotspot regions, one containing four genes associated with altitudinal variation, and the other containing ten genes associated with response to solar radiation. These genes participate in abiotic stress resistance and regulation of reproductive processes. Our results provide insight into the genetic mechanisms underlying high-altitude adaptation in Tibetan poplar.

Fine-scale landscape genetics unveiling contemporary asymmetric movement of red panda (Ailurus fulgens) in Kangchenjunga landscape, India

Wildlife management in rapid changing landscapes requires critical planning through cross cutting networks, and understanding of landscape features, often affected by the anthropogenic activities. The present study demonstrates fine-scale spatial patterns of genetic variation and contemporary gene flow of red panda (Ailurus fulgens) populations with respect to landscape connectivity in Kangchenjunga Landscape (KL), India. The study found about 1,309.54 km2 area suitable for red panda in KL-India, of which 62.21% area fell under the Protected Area network. We identified 24 unique individuals from 234 feces collected at nine microsatellite loci. The spatially explicit and non-explicit Bayesian clustering algorithms evident to exhibit population structuring and supported red panda populations to exist in meta-population frame work. In concurrence to the habitat suitability and landscape connectivity models, gene flow results supported a contemporary asymmetric movement of red panda by connecting KL-India in a crescent arc. We demonstrate the structural-operational connectivity of corridors in KL-India that facilitated red panda movement in the past. We also seek for cooperation in Nepal, Bhutan and China to aid in preparing for a comprehensive monitoring plan for the long-term conservation and management of red panda in trans-boundary landscapes.

Disease swamps molecular signatures of genetic-environmental associations to abiotic factors in Tasmanian devil (Sarcophilus harrisii) populations

Landscape genomics studies focus on identifying candidate genes under selection via spatial variation in abiotic environmental variables, but rarely by biotic factors (i.e., disease). The Tasmanian devil (Sarcophilus harrisii) is found only on the environmentally heterogeneous island of Tasmania and is threatened with extinction by a transmissible cancer, devil facial tumor disease (DFTD). Devils persist in regions of long-term infection despite epidemiological model predictions of species' extinction, suggesting possible adaptation to DFTD. Here, we test the extent to which spatial variation and genetic diversity are associated with the abiotic environment (i.e., climatic variables, elevation, vegetation cover) and/or DFTD. We employ genetic-environment association analyses using 6886 SNPs from 3287 individuals sampled pre- and post-disease arrival across the devil's geographic range. Pre-disease, we find significant correlations of allele frequencies with environmental variables, including 365 unique loci linked to 71 genes, suggesting local adaptation to abiotic environment. The majority of candidate loci detected pre-DFTD are not detected post-DFTD arrival. Several post-DFTD candidate loci are associated with disease prevalence and were in linkage disequilibrium with genes involved in tumor suppression and immune response. Loss of apparent signal of abiotic local adaptation post-disease suggests swamping by strong selection resulting from the rapid onset of DFTD.

Whole Genome Scan Reveals Molecular Signatures of Divergence and Selection Related to Important Traits in Durum Wheat Germplasm

The first breeding program in the world for durum wheat was conceived in Italy in the early 1900s. Over the decades, pressure exerted by natural and artificial selection could have progressively reduced the genetic diversity of the durum wheat germplasm. In the present study, a large panel of Italian durum wheat accessions that includes landraces, old and modern cultivars was subjected to genotyping using the Illumina iSelect 15K wheat SNP array. The aim was to assess the impact that selection has in shaping Italian durum wheat genetic diversity and to exploit the patterns of genetic diversity between populations to identify molecular signatures of divergence and selection. Relatively small differences in genetic diversity have been observed among accessions, which have been selected and cultivated in Italy over the past 150 years. Indeed, directional selection combined with that operated by farmers/breeders resulted in the increase of linkage disequilibrium (LD) and in changes of the allelic frequencies in DNA regions that control important agronomic traits. Results from this study also show that major well-known genes and/or QTLs affecting plant height (RHT), earliness (VRN, PPD) and grain quality (GLU, PSY, PSD, LYC, PPO, LOX3) co-localized with outlier SNP loci. Interestingly, many of these SNPs fall in genomic regions where genes involved in nitrogen metabolism are. This finding highlights the key role these genes have played in the transition from landraces to modern cultivars. Finally, our study remarks on the need to fully exploit the genetic diversity of Italian landraces by intense pre-breeding activities aimed at introducing a new source of adaptability and resistance in the genetic background of modern cultivars, to contrast the effect of climate change. The list of divergent loci and loci under selection associated with useful agronomic traits represents an invaluable resource to detect new allelic variants for target genes and for guiding new genomic selection programs in durum wheat.

Current State and Perspectives in Population Genomics of the Common Bean

* Correspondence: r [...].

Knowledge status and sampling strategies to maximize cost-benefit ratio of studies in landscape genomics of wild plants

To avoid local extinction due to the changes in their natural ecosystems, introduced by anthropogenic activities, species undergo local adaptation. Landscape genomics approach, through genome–environment association studies, has helped evaluate the local adaptation in natural populations. Landscape genomics, is still a developing discipline, requiring refinement of guidelines in sampling design, especially for studies conducted in the backdrop of stark socioeconomic realities of the rainforest ecologies, which are global biodiversity hotspots. In this study we aimed to devise strategies to improve the cost-benefit ratio of landscape genomics studies by surveying sampling designs and genome sequencing strategies used in existing studies. We conducted meta-analyses to evaluate the importance of sampling designs, in terms of (i) number of populations sampled, (ii) number of individuals sampled per population, (iii) total number of individuals sampled, and (iv) number of SNPs used in different studies, in discerning the molecular mechanisms underlying local adaptation of wild plant species. Using the linear mixed effects model, we demonstrated that the total number of individuals sampled and the number of SNPs used, significantly influenced the detection of loci underlying the local adaptation. Thus, based on our findings, in order to optimize the cost-benefit ratio of landscape genomics studies, we suggest focusing on increasing the total number of individuals sampled and using a targeted (e.g. sequencing capture) Pool-Seq approach and/or a random (e.g. RAD-Seq) Pool-Seq approach to detect SNPs and identify SNPs under selection for a given environmental cline. We also found that the existing molecular evidences are inadequate in predicting the local adaptations to climate change in tropical forest ecosystems.

Sampling strategy optimization to increase statistical power in landscape genomics: A simulation-based approach

An increasing number of studies are using landscape genomics to investigate local adaptation in wild and domestic populations. Implementation of this approach requires the sampling phase to consider the complexity of environmental settings and the burden of logistical constraints. These important aspects are often underestimated in the literature dedicated to sampling strategies. In this study, we computed simulated genomic data sets to run against actual environmental data in order to trial landscape genomics experiments under distinct sampling strategies. These strategies differed by design approach (to enhance environmental and/or geographical representativeness at study sites), number of sampling locations and sample sizes. We then evaluated how these elements affected statistical performances (power and false discoveries) under two antithetical demographic scenarios. Our results highlight the importance of selecting an appropriate sample size, which should be modified based on the demographic characteristics of the studied population. For species with limited dispersal, sample sizes above 200 units are generally sufficient to detect most adaptive signals, while in random mating populations this threshold should be increased to 400 units. Furthermore, we describe a design approach that maximizes both environmental and geographical representativeness of sampling sites and show how it systematically outperforms random or regular sampling schemes. Finally, we show that although having more sampling locations (between 40 and 50 sites) increase statistical power and reduce false discovery rate, similar results can be achieved with a moderate number of sites (20 sites). Overall, this study provides valuable guidelines for optimizing sampling strategies for landscape genomics experiments.

Selection at behavioural, developmental and metabolic genes is associated with the northward expansion of a successful tropical colonizer

What makes a species able to colonize novel environments? This question is key to understand the dynamics of adaptive radiations and ecological niche shifts, but the mechanisms that underlie expansion into novel habitats remain poorly understood at a genomic scale. Lizards from the genus Anolis are typically tropical, and the green anole (Anolis carolinensis) constitutes an exception since it expanded into temperate North America from subtropical Florida. Thus, we used the green anole as a model to investigate signatures of selection associated with colonization of a new environment, namely temperate North America. To this end, we analysed 29 whole-genome sequences, covering the entire native range of the species. We used a combination of recent methods to quantify both positive and balancing selection in northern populations, including F_ST outlier methods, machine learning and ancestral recombination graphs. We naively scanned for genes of interest and assessed the overlap between multiple tests. Strikingly, we identified many genes involved in behaviour, suggesting that the recent successful colonization of northern environments may have been linked to behavioural shifts as well as physiological adaptation. Using a candidate genes strategy, we determined that genes involved in response to cold or behaviour displayed more frequently signals of selection, while controlling for local recombination rate, gene clustering and gene length. In addition, we found signatures of balancing selection at immune genes in all investigated genetic groups, but also at genes involved in neuronal and anatomical development.

Geographic separation and genetic differentiation of populations are not coupled with niche differentiation in threatened Kaiser's spotted newt (Neurergus kaiseri)

The combination of niche modelling and landscape genetics (genomics) helps to disentangle processes that have shaped population structure in the evolutionary past and presence of species. Herein, we integrate a comprehensive genomic dataset with ecological parameters and niche modelling for the threatened Kaiser's newt, a newt species adapted to mountain spring-ponds in Iran. Genomic analysis suggests the existence of two highly differentiated clades North and South of the Dez River. Genetic variation between the two clades (76.62%) was much greater than within clades (16.25%), suggesting that the Dez River prevented gene flow. River disconnectivity, followed by geographic distance, contributed mostly to genetic differentiation between populations. Environmental niche and landscape resistance had no significant influence. Though a significant difference between climatic niches occupied by each clade at the landscape-scale, habitat niches at the local-scale were equivalent. 'Niche similarity analysis' supported niche conservatism between the two clades despite the southward shift in the climatic niche of the Southern clade. Accordingly, populations of different clades may occupy different climatic niches within their ancestral niche. Our results indicate that the change of climatic conditions of geographically and genetically separated populations does not necessarily result in the shift of an ecological niche.

Genome-wide signatures of local adaptation among seven stoneflies species along a nationwide latitudinal gradient in Japan

Background

Environmental heterogeneity continuously produces a selective pressure that results in genomic variation among organisms; understanding this relationship remains a challenge in evolutionary biology. Here, we evaluated the degree of genome-environmental association of seven stonefly species across a wide geographic area in Japan and additionally identified putative environmental drivers and their effect on co-existing multiple stonefly species. Double-digest restriction-associated DNA (ddRAD) libraries were independently sequenced for 219 individuals from 23 sites across four geographical regions along a nationwide latitudinal gradient in Japan.

Results

A total of 4251 candidate single nucleotide polymorphisms (SNPs) strongly associated with local adaptation were discovered using Latent mixed models; of these, 294 SNPs showed strong correlation with environmental variables, specifically precipitation and altitude, using distance-based redundancy analysis. Genome–genome comparison among the seven species revealed a high sequence similarity of candidate SNPs within a geographical region, suggesting the occurrence of a parallel evolution process.

Conclusions

Our results revealed genomic signatures of local adaptation and their influence on multiple, co-occurring species. These results can be potentially applied for future studies on river management and climatic stressor impacts.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5453-3) contains supplementary material, which is available to authorized users.