Shared and nonshared genomic divergence in parallel ecotypes of Littorina saxatilis at a local scale

The succession of ecological divergence and reproductive isolation in adaptive radiations

Adaptive radiation is a major source of biodiversity but the way in which known components of ecological opportunity, ecological differentiation, and reproductive isolation underpin such biodiversity patterns remains elusive. Much is known about the evolution of ecological differentiation and reproductive isolation during single speciation events, but exactly how those processes scale up to complete adaptive radiations is less understood. Do we expect complete reproductive barriers between newly formed species before the ecological differentiation continues, or does proper species formation occur much later, long after the ecological diversification? Our goal is to improve our mechanistic understanding of adaptive radiations by analyzing an individual-based model that includes a suite of mechanisms that are known to contribute to biodiversity. The model includes variable biogeographic settings, ecological opportunities, and types of mate choice, which makes several different scenarios of an adaptive radiation possible. We find that evolving clades tend to exploit ecological opportunities early whereas reproductive barriers evolve later, demonstrating a decoupling of ecological differentiation and species formation. In many cases, we also find a long-term trend where assortative mating associated with ecological traits is replaced by sexual selection of neutral display traits as the primary mechanism for reproductive isolation. Our results propose that reticulate phylogenies are likely common and stem from initially low reproductive barriers, rather than the previously suggested idea of repeated hybridization events between well-separated species.

Genetic assessment of eight zoo populations of golden snub-nosed monkey (Rhinopithecus roxellana) implication to the conservation management of captive populations

Captive breeding programs play an important role in preserving the genetic diversity of endangered species. It is of utmost importance to conduct genetic assessment for captive populations in order to develop scientific breeding plans and conservation management strategies. Here, we genotyped 10 microsatellite loci and sequenced 368 bp of mitochondrial DNA control region for the golden snub-nosed monkey (Rhinopithecus roxellana) from eight captive populations in China, and compared the genetic indices of captive populations with a wild population. Meanwhile, we performed paternity tests to verify the genealogical records and established genetic lineages. A total of 157 individuals were identified from 161 fecal samples, including 135 captive individuals (approximately 25% of captive individuals in China). Microsatellite analysis showed that the nine populations had moderate levels of genetic diversity, with polymorphism information content (PIC) ranging from 0.43 to 0.542; the genetic diversity of captive populations (average PIC: 0.503) was slightly higher than that of the wild population (PIC: 0.438). The Structure analysis indicated that individuals of the eight captive populations contained two different genetic components. We conducted either single-blind or double-blind paternity testing on 40 offspring of captive individuals and found that five offspring from two zoos (Nanjing Hongshan Forest Zoo and Shanghai Wild Animal Park) showed discrepant kinships from their pedigree records, probably due to the inaccuracies in pedigree records. By constructing genetic pedigrees, inbred offspring were found in Beijing Zoo, Shanghai Zoo, Hangzhou Zoo, and Chengdu Zoo. Analysis based on mitochondrial DNA showed a high level of genetic diversity in the eight captive populations (mean nucleotide diversity: 0.047). However, no nucleotide diversity was found in the wild population. This study conducted a genetic survey for captive golden snub-nosed monkeys and will significantly benefit the genetic conservation management for captive populations in the future.

Genomic underpinnings of head and body shape in Arctic charr ecomorph pairs

Across its Holarctic range, Arctic charr (Salvelinus alpinus) populations have diverged into distinct trophic specialists across independent replicate lakes. The major aspect of divergence between ecomorphs is in head shape and body shape, which are ecomorphological traits reflecting niche use. However, whether the genomic underpinnings of these parallel divergences are consistent across replicates was unknown but key for resolving the substrate of parallel evolution. We investigated the genomic basis of head shape and body shape morphology across four benthivore-planktivore ecomorph pairs of Arctic charr in Scotland. Through genome-wide association analyses, we found genomic regions associated with head shape (89 SNPs) or body shape (180 SNPs) separately and 50 of these SNPs were strongly associated with both body and head shape morphology. For each trait separately, only a small number of SNPs were shared across all ecomorph pairs (3 SNPs for head shape and 10 SNPs for body shape). Signs of selection on the associated genomic regions varied across pairs, consistent with evolutionary demography differing considerably across lakes. Using a comprehensive database of salmonid QTLs newly augmented and mapped to a charr genome, we found several of the head- and body-shape-associated SNPs were within or near morphology QTLs from other salmonid species, reflecting a shared genetic basis for these phenotypes across species. Overall, our results demonstrate how parallel ecotype divergences can have both population-specific and deeply shared genomic underpinnings across replicates, influenced by differences in their environments and demographic histories.

Genomic evidence for domestication selection in three hatchery populations of Chinook salmon, Oncorhynchus tshawytscha

Fish hatcheries are widely used to enhance fisheries and supplement declining wild populations. However, substantial evidence suggests that hatchery fish are subject to differential selection pressures compared to their wild counterparts. Domestication selection, or adaptation to the hatchery environment, poses a risk to wild populations if traits specific to success in the hatchery environment have a genetic component and there is subsequent introgression between hatchery and wild fish. Few studies have investigated domestication selection in hatcheries on a genomic level, and even fewer have done so in parallel across multiple hatchery-wild population pairs. In this study, we used low-coverage whole-genome sequencing to investigate signals of domestication selection in three separate hatchery populations of Chinook salmon, Oncorhynchus tshawytscha, after approximately seven generations of divergence from their corresponding wild progenitor populations. We sequenced 192 individuals from populations across Southeast Alaska and estimated genotype likelihoods at over six million loci. We discovered a total of 14 outlier peaks displaying high genetic differentiation (F ST) between hatchery-wild pairs, although no peaks were shared across the three comparisons. Peaks were small (53 kb on average) and often displayed elevated absolute genetic divergence (D xy) and linkage disequilibrium, suggesting some level of domestication selection has occurred. Our study provides evidence that domestication selection can lead to genetic differences between hatchery and wild populations in only a few generations. Additionally, our data suggest that population-specific adaptation to hatchery environments likely occurs through different genetic pathways, even for populations with similar standing genetic variation. These results highlight the need to collect paired genotype-phenotype data to understand how domestication may be affecting fitness and to identify potential management practices that may mitigate genetic risks despite multiple pathways of domestication.

Integrating Pool-seq uncertainties into demographic inference

Next-generation sequencing of pooled samples (Pool-seq) is a popular method to assess genome-wide diversity patterns in natural and experimental populations. However, Pool-seq is associated with specific sources of noise, such as unequal individual contributions. Consequently, using Pool-seq for the reconstruction of evolutionary history has remained underexplored. Here we describe a novel Approximate Bayesian Computation (ABC) method to infer demographic history, explicitly modelling Pool-seq sources of error. By jointly modelling Pool-seq data, demographic history and the effects of selection due to barrier loci, we obtain estimates of demographic history parameters accounting for technical errors associated with Pool-seq. Our ABC approach is computationally efficient as it relies on simulating subsets of loci (rather than the whole-genome) and on using relative summary statistics and relative model parameters. Our simulation study results indicate Pool-seq data allows distinction between general scenarios of ecotype formation (single versus parallel origin) and to infer relevant demographic parameters (e.g. effective sizes and split times). We exemplify the application of our method to Pool-seq data from the rocky-shore gastropod Littorina saxatilis, sampled on a narrow geographical scale at two Swedish locations where two ecotypes (Wave and Crab) are found. Our model choice and parameter estimates show that ecotypes formed before colonization of the two locations (i.e. single origin) and are maintained despite gene flow. These results indicate that demographic modelling and inference can be successful based on pool-sequencing using ABC, contributing to the development of suitable null models that allow for a better understanding of the genetic basis of divergent adaptation.

Understanding natural selection and similarity: Convergent, parallel and repeated evolution

Parallel and convergent evolution offer some of the most compelling evidence for the significance of natural selection in evolution, as the emergence of similar adaptive solutions is unlikely to occur by random chance alone. However, these terms are often employed inconsistently, leading to misinterpretation and confusion, and recently proposed definitions have unintentionally diminished the emphasis on the evolution of similar adaptive solutions. Here, I examine various conceptual frameworks and definitions related to parallel and convergent evolution and propose a consolidated framework that enhances our comprehension of these evolutionary patterns. The primary aim of this framework is to harmonize the concepts of parallel and convergent evolution together with natural selection and the idea of similarity. Both concepts involve the evolution of similar adaptive solutions as a result of environmental challenges. The distinction lies in ancestral phenotypes. Parallel evolution takes place when the ancestral phenotypes (before selection) of the lineages are similar. Convergent evolution happens when the lineages have distinct ancestral phenotypes (before selection). Because an ancestral-based distinction will inevitably lead to cases where uncertainty in the distinction may arise, the framework includes a general term, repeated evolution, which can be used as a term applying to the evolution of similar phenotypes and genotypes as well as similar responses to environmental pressures. Based on the argument that genetic similarity may frequently arise without selection, the framework posits that the similarity of genetic sequences is not of great interest unless linked to the actions of natural selection or to the origins (mutation, standing genetic variation, gene flow) and locations of the similar sequences.

Polygenic routes lead to parallel altitudinal adaptation in Heliosperma pusillum (Caryophyllaceae)

Understanding how organisms adapt to the environment is a major goal of modern biology. Parallel evolution-the independent evolution of similar phenotypes in different populations-provides a powerful framework to investigate the evolutionary potential of populations, the constraints of evolution, its repeatability and therefore its predictability. Here, we quantified the degree of gene expression and functional parallelism across replicated ecotype formation in Heliosperma pusillum (Caryophyllaceae), and gained insights into the architecture of adaptive traits. Population structure analyses and demographic modelling support a previously formulated hypothesis of parallel polytopic divergence of montane and alpine ecotypes. We detect a large proportion of differentially expressed genes (DEGs) underlying divergence within each replicate ecotype pair, with a strikingly low number of shared DEGs across pairs. Functional enrichment of DEGs reveals that the traits affected by significant expression divergence are largely consistent across ecotype pairs, in strong contrast to the nonshared genetic basis. The remarkable redundancy of differential gene expression indicates a polygenic architecture for the diverged adaptive traits. We conclude that polygenic traits appear key to opening multiple routes for adaptation, widening the adaptive potential of organisms.

Using pose estimation to identify regions and points on natural history specimens

A key challenge in mobilising growing numbers of digitised biological specimens for scientific research is finding high-throughput methods to extract phenotypic measurements on these datasets. In this paper, we test a pose estimation approach based on Deep Learning capable of accurately placing point labels to identify key locations on specimen images. We then apply the approach to two distinct challenges that each requires identification of key features in a 2D image: (i) identifying body region-specific plumage colouration on avian specimens and (ii) measuring morphometric shape variation in Littorina snail shells. For the avian dataset, 95% of images are correctly labelled and colour measurements derived from these predicted points are highly correlated with human-based measurements. For the Littorina dataset, more than 95% of landmarks were accurately placed relative to expert-labelled landmarks and predicted landmarks reliably captured shape variation between two distinct shell ecotypes ('crab' vs 'wave'). Overall, our study shows that pose estimation based on Deep Learning can generate high-quality and high-throughput point-based measurements for digitised image-based biodiversity datasets and could mark a step change in the mobilisation of such data. We also provide general guidelines for using pose estimation methods on large-scale biological datasets.

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.

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.

Divergent and non-parallel evolution of MHC IIB in the Neotropical Midas cichlid species complex

Background

Ecological diversification is the result of divergent natural selection by contrasting habitat characteristics that favours the evolution of distinct phenotypes. This process can happen in sympatry and in allopatry. Habitat-specific parasite communities have the potential to drive diversification among host populations by imposing selective pressures on their host's immune system. In particular, the hyperdiverse genes of the major histocompatibility complex (MHC) are implicated in parasite-mediated host divergence. Here, we studied the extent of divergence at MHC, and discuss how it may have contributed to the Nicaraguan Midas cichlid species complex diversification, one of the most convincing examples of rapid sympatric parallel speciation.

Results

We genotyped the MHC IIB for individuals from six sympatric Midas cichlid assemblages, each containing species that have adapted to exploit similar habitats. We recovered large allelic and functional diversity within the species complex. While most alleles were rare, functional groups of alleles (supertypes) were common, suggesting that they are key to survival and that they were maintained during colonization and subsequent radiations. We identified lake-specific and habitat-specific signatures for both allelic and functional diversity, but no clear pattern of parallel divergence among ecomorphologically similar phenotypes.

Conclusions

Colonization and demographic effects of the fish could have contributed to MHC evolution in the Midas cichlid in conjunction with habitat-specific selective pressures, such as parasites associated to alternative preys or environmental features. Additional ecological data will help evaluating the role of host–parasite interactions in the Midas cichlid radiations and aid in elucidating the potential role of non-parallel features differentiating crater lake species assemblages.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01997-9.

Extensive introgression at late stages of species formation: Insights from grasshopper hybrid zones

The process of species formation is characterised by the accumulation of multiple reproductive barriers. The evolution of hybrid male sterility, or Haldane's rule, typically characterises later stages of species formation, when reproductive isolation is strongest. Yet, understanding how quickly reproductive barriers evolve and their consequences for maintaining genetic boundaries between emerging species remains a challenging task because it requires studying taxa that hybridise in nature. Here, we address these questions using the meadow grasshopper Pseudochorthippus parallelus, where populations that show multiple reproductive barriers, including hybrid male sterility, hybridise in two natural hybrid zones. Using mitochondrial data, we infer that such populations have diverged some 100,000 years ago, at the beginning of the last glacial cycle in Europe. Nuclear data shows that contractions at multiple glacial refugia, and post-glacial expansions have facilitated genetic differentiation between lineages that today interact in hybrid zones. We find extensive introgression throughout the sampled species range, irrespective of current strength of reproductive isolation. Populations exhibiting hybrid male sterility in two hybrid zones show repeatable patterns of genomic differentiation, consistent with shared genomic constraints affecting ancestral divergence or with the role of those regions in reproductive isolation. Together, our results suggest that reproductive barriers that characterise late stages of species formation can evolve relatively quickly, particularly when associated with strong demographic changes. Moreover, we show that such barriers persist in the face of extensive gene flow, allowing future studies to identify associated genomic regions.

Comparing genome scans among species of the stickleback order reveals three different patterns of genetic diversity

Comparing genome scans among species is a powerful approach for investigating the patterns left by evolutionary processes. In particular, this offers a way to detect candidate genes that drive convergent evolution. We compared genome scan results to investigate if patterns of genetic diversity and divergence are shared among divergent species within the stickleback order (Gasterosteiformes): the threespine stickleback (Gasterosteus aculeatus), ninespine stickleback (Pungitius pungitus), and tubesnout (Aulorhynchus flavidus). Populations were sampled from the southern and northern edges of each species' range, to identify patterns associated with latitudinal changes in genetic diversity. Weak correlations in genetic diversity (F ST and expected heterozygosity) and three different patterns in the genomic landscape were found among these species. Additionally, no candidate genes for convergent evolution were detected. This is a counterexample to the growing number of studies that have shown overlapping genetic patterns, demonstrating that genome scan comparisons can be noisy due to the effects of several interacting evolutionary forces.

Parallel shape divergence between ecotypes of the limpet Nacella concinna along the Antarctic Peninsula: a new model species for parallel evolution?

Parallel phenotypic divergence is the independent differentiation between phenotypes of the same lineage or species occupying ecologically similar environments in different populations. We tested in the Antarctic limpet Nacella concinna the extent of parallel morphological divergence in littoral and sublittoral ecotypes throughout its distribution range. These ecotypes differ in morphological, behavioural and physiological characteristics. We studied the lateral and dorsal outlines of shells and the genetic variation of the mitochondrial gene Cytochrome Oxidase subunit I from both ecotypes in 17 sample sites along more than 2,000 km. The genetic data indicate that both ecotypes belong to a single evolutionary lineage. The magnitude and direction of phenotypic variation differ between ecotypes across sample sites; completely parallel ecotype-pairs (i.e., they diverge in the same magnitude and in the same direction) were detected in 84.85% of lateral and 65.15% in dorsal view comparisons. Besides, specific traits (relative shell height, position of shell apex, and elliptical/pear-shape outline variation) showed high parallelism. We observed weak morphological covariation between the two shape shell views, indicating that distinct evolutionary forces and environmental pressures could be acting on this limpet shell shape. Our results demonstrate there is a strong parallel morphological divergence pattern in N. concinna along its distribution, making this Antarctic species a suitable model for the study of different evolutionary forces shaping the shell evolution of this limpet.

Whole-genome resequencing confirms reproductive isolation between sympatric demes of brown trout (Salmo trutta) detected with allozymes

The sympatric existence of genetically distinguishable populations of the same species remains a puzzle in ecology. Coexisting salmonid fish populations are known from over 100 freshwater lakes. Most studies of sympatric populations have used limited numbers of genetic markers making it unclear if genetic divergence involves certain parts of the genome. We returned to the first reported case of salmonid sympatry, initially detected through contrasting homozygosity at a single allozyme locus (coding for lactate dehydrogenase A) in brown trout in the small Lakes Bunnersjöarna, Sweden. First, we verified the existence of the two coexisting demes using a 96-SNP fluidigm array. We then applied whole-genome resequencing of pooled DNA to explore genome-wide diversity within and between these demes; nucleotide diversity was higher in deme I than in deme II. Strong genetic divergence is observed with genome-wide F_ST  ≈ 0.2. Compared with data from populations of similar small lakes, this divergence is of similar magnitude as that between reproductively isolated populations. Individual whole-genome resequencing of two individuals per deme suggests higher inbreeding in deme II versus deme I, indicating different degree of isolation. We located two gene-copies for LDH-A and found divergence between demes in a regulatory section of one of these genes. However, we did not find a perfect fit between the sequence data and previous allozyme results, and this will require further research. Our data demonstrates genome-wide divergence governed mostly by genetic drift but also by diversifying selection in coexisting populations. This type of hidden biodiversity needs consideration in conservation management.

Rapid Parallel Adaptation to Anthropogenic Heavy Metal Pollution

The impact of human-mediated environmental change on the evolutionary trajectories of wild organisms is poorly understood. In particular, capacity of species to adapt rapidly (in hundreds of generations or less), reproducibly and predictably to extreme environmental change is unclear. Silene uniflora is predominantly a coastal species, but it has also colonized isolated, disused mines with phytotoxic, zinc-contaminated soils. To test whether rapid, parallel adaptation to anthropogenic pollution has taken place, we used reduced representation sequencing (ddRAD) to reconstruct the evolutionary history of geographically proximate mine and coastal population pairs and found largely independent colonization of mines from different coastal sites. Furthermore, our results show that parallel evolution of zinc tolerance has occurred without gene flow spreading adaptive alleles between mine populations. In genomic regions where signatures of selection were detected across multiple mine-coast pairs, we identified genes with functions linked to physiological differences between the putative ecotypes, although genetic differentiation at specific loci is only partially shared between mine populations. Our results are consistent with a complex, polygenic genetic architecture underpinning rapid adaptation. This shows that even under a scenario of strong selection and rapid adaptation, evolutionary responses to human activities (and other environmental challenges) may be idiosyncratic at the genetic level and, therefore, difficult to predict from genomic data.

Polygenic basis and the role of genome duplication in adaptation to similar selective environments

Genetic changes underlying adaptation vary greatly in terms of complexity and, within the same species, genetic responses to similar selective pressures may or may not be the same. We examine both complex (supergene) and simple (SNP) genetic variants occurring in populations of rainbow trout (Oncorhynchus mykiss) independently isolated from ocean access and compared them to each other and to an anadromous below-barrier population representing their ancestral source to search for signatures of both parallel and non-parallel adaptation. All landlocked populations displayed an increased frequency of a large inversion on chromosome Omy05, while three of the four populations exhibited elevated frequencies of another inversion located on chromosome Omy20. In addition, we identified numerous regions outside these two inversions that also show significant shifts in allele frequencies consistent with adaptive evolution. However, there was little concordance among above-barrier populations in these specific genomic regions under selection. In part, the lack of concordance appears to arise from ancestral autopolyploidy in rainbow trout that provides duplicate genomic regions of similar functional composition for selection to act upon. Thus, while selection acting on landlocked populations universally favors the resident ecotype, outside of the major chromosomal inversions, the resulting genetic changes are largely distinct among populations. Our results indicate that selection on standing genetic variation is likely the primary mode of rapid adaptation, and that both supergene complexes and individual loci contribute to adaptive evolution, further highlighting the diversity of adaptive genomic variation involved in complex phenotypic evolution.

Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing

BACKGROUND

Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS.

RESULTS

In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA.

CONCLUSIONS

Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.

Parallel adaptation in autopolyploid Arabidopsis arenosa is dominated by repeated recruitment of shared alleles

Relative contributions of pre-existing vs de novo genomic variation to adaptation are poorly understood, especially in polyploid organisms. We assess this in high resolution using autotetraploid Arabidopsis arenosa, which repeatedly adapted to toxic serpentine soils that exhibit skewed elemental profiles. Leveraging a fivefold replicated serpentine invasion, we assess selection on SNPs and structural variants (TEs) in 78 resequenced individuals and discover significant parallelism in candidate genes involved in ion homeostasis. We further model parallel selection and infer repeated sweeps on a shared pool of variants in nearly all these loci, supporting theoretical expectations. A single striking exception is represented by TWO PORE CHANNEL 1, which exhibits convergent evolution from independent de novo mutations at an identical, otherwise conserved site at the calcium channel selectivity gate. Taken together, this suggests that polyploid populations can rapidly adapt to environmental extremes, calling on both pre-existing variation and novel polymorphisms.

Relative contributions of pre-existing versus de novo genomic variation to adaptation remain unclear. Here, the authors address this problem by examining the adaptation of autotetraploid Arabidopsis arenosa to serpentine soils and find that both types of variations contribute to rapid adaptation.

Linking ecology, morphology, and metabolism: Niche differentiation in sympatric populations of closely related species of the genus Littorina (Neritrema)

Divergence of ecological niches in phylogenetically closely related species indicates the importance of ecology in speciation, especially for sympatric species are considered. Such ecological diversification provides an advantage of alleviating interspecies competition and promotes more efficient exploitation of environmental resources, thus being a basis for ecological speciation. We analyzed a group of closely related species from the subgenus Neritrema (genus Littorina, Caenogastropoda) from the gravel-bouldery shores. In two distant sites at the Barents and Norwegian Sea, we examined the patterns of snail distribution during low tide (quantitative sampling stratified by intertidal level, presence of macrophytes, macrophyte species, and position on them), shell shape and its variability (geometric morphometrics), and metabolic characteristics (metabolomic profiling). The studied species diversified microbiotopes, which imply an important role of ecological specification in the recent evolution of this group. The only exception to this trend was the species pair L. arcana / L. saxatilis, which is specifically discussed. The ecological divergence was accompanied by differences in shell shape and metabolomic characteristics. Significant differences were found between L. obtusata versus L. fabalis and L. saxatilis / L. arcana versus L. compressa both in shell morphology and in metabolomes. L. saxatilis demonstrated a clear variability depending on intertidal level which corresponds to a shift in conditions within the occupied microhabitat. Interestingly, the differences between L. arcana (inhabiting the upper intertidal level) and L. compressa (inhabiting the lower one) were analogous to those between the upper and lower fractions of L. saxatilis. No significant level-dependent changes were found between the upper and lower fractions of L. obtusata, most probably due to habitat amelioration by fucoid macroalgae. All these results are discussed in the contexts of the role of ecology in speciation, ecological niche dynamics and conservatism, and evolutionary history of the Neritrema species.

Using replicate hybrid zones to understand the genomic basis of adaptive divergence

Combining hybrid zone analysis with genomic data is a promising approach to understanding the genomic basis of adaptive divergence. It allows for the identification of genomic regions underlying barriers to gene flow. It also provides insights into spatial patterns of allele frequency change, informing about the interplay between environmental factors, dispersal and selection. However, when only a single hybrid zone is analysed, it is difficult to separate patterns generated by selection from those resulting from chance. Therefore, it is beneficial to look for repeatable patterns across replicate hybrid zones in the same system. We applied this approach to the marine snail Littorina saxatilis, which contains two ecotypes, adapted to wave-exposed rocks vs. high-predation boulder fields. The existence of numerous hybrid zones between ecotypes offered the opportunity to test for the repeatability of genomic architectures and spatial patterns of divergence. We sampled and phenotyped snails from seven replicate hybrid zones on the Swedish west coast and genotyped them for thousands of single nucleotide polymorphisms. Shell shape and size showed parallel clines across all zones. Many genomic regions showing steep clines and/or high differentiation were shared among hybrid zones, consistent with a common evolutionary history and extensive gene flow between zones, and supporting the importance of these regions for divergence. In particular, we found that several large putative inversions contribute to divergence in all locations. Additionally, we found evidence for consistent displacement of clines from the boulder-rock transition. Our results demonstrate patterns of spatial variation that would not be accessible without continuous spatial sampling, a large genomic data set and replicate hybrid zones.

Phylogeographic analysis of Saxifraga fortunei complex (Saxifragaceae) reveals multiple origins of morphological and ecological variations in the Japanese Archipelago

Phenotypic polymorphism within a species is a notable phenomenon in evolutionary biology to understand the process of adaptive speciation and other historical events. The Saxifraga fortunei complex is a widespread herb found in East Asia. It includes several ecotypic taxa corresponding to their habitat environments. The distribution of the various ecotypes in a limited area of the Japanese Archipelago makes the species a suitable model to investigate the impact of population demographic history and natural selection on lineage diversification. Here, Sanger-based sequencing was used to estimate the divergence timeframe between populations of the Eurasian continent and Japan. Genome-wide SNPs obtained by ddRAD sequencing were used to investigate the phylogeographic origins of ecotypic taxa. The phylogenetic analyses revealed the divergence of the Japanese population from the continental population in the late Miocene. Two distinct regional clades of North and South Japan were identified; phenotypic diversification was evident only in the southern clade. The South Japan clades displayed a historical distribution expansion from north to south. The phenotypic variations appeared to have generated during the expansion. The ecotypic boundaries were incongruent with the genetic grouping. We propose that morphological and ecological specialization in Japanese populations was repeatedly generated by local natural selection.

The genomics of phenotypically differentiated Asellus aquaticus cave, surface stream and lake ecotypes

Organisms well suited for the study of ecotype formation have wide distribution ranges, where they adapt to multiple drastically different habitats repeatedly over space and time. Here we study such ecotypes in a Crustacean model, Asellus aquaticus, a commonly occurring isopod found in freshwater habitats as diverse as streams, caves and lakes. Previous studies focusing on cave vs. surface ecotypes have attributed depigmentation, eye loss and prolonged antennae to several south European cave systems. Likewise, surveys across multiple Swedish lakes have identified the presence of dark-pigmented "reed" and light-pigmented "stonewort" ecotypes, which can be found within the same lake. In this study, we sequenced the first draft genome of A. aquaticus, and subsequently use this to map reads and call variants in surface stream, cave and two lake ecotypes. In addition, the draft genome was combined with a RADseq approach to perform a quantitative trait locus (QTL) mapping study using a laboratory bred F₂ and F₄ cave × surface intercross. We identified genomic regions associated with body pigmentation, antennae length and body size. Furthermore, we compared genome-wide differentiation between natural populations and found several genes potentially associated with these habitats. The assessment of the cave QTL regions in the light-dark comparison of lake populations suggests that the regions associated with cave adaptation are also involved with genomic differentiation in the lake ecotypes. These demonstrate how troglomorphic adaptations can be used as a model for related ecotype formation.

Genetic variation for adaptive traits is associated with polymorphic inversions in Littorina saxatilis

Chromosomal inversions have long been recognized for their role in local adaptation. By suppressing recombination in heterozygous individuals, they can maintain coadapted gene complexes and protect them from homogenizing effects of gene flow. However, to fully understand their importance for local adaptation we need to know their influence on phenotypes under divergent selection. For this, the marine snail Littorina saxatilis provides an ideal study system. Divergent ecotypes adapted to wave action and crab predation occur in close proximity on intertidal shores with gene flow between them. Here, we used F2 individuals obtained from crosses between the ecotypes to test for associations between genomic regions and traits distinguishing the Crab‐/Wave‐adapted ecotypes including size, shape, shell thickness, and behavior. We show that most of these traits are influenced by two previously detected inversion regions that are divergent between ecotypes. We thus gain a better understanding of one important underlying mechanism responsible for the rapid and repeated formation of ecotypes: divergent selection acting on inversions. We also found that some inversions contributed to more than one trait suggesting that they may contain several loci involved in adaptation, consistent with the hypothesis that suppression of recombination within inversions facilitates differentiation in the presence of gene flow.

Transcriptome Expression of Biomineralization Genes in Littoraria flava Gastropod in Brazilian Rocky Shore Reveals Evidence of Local Adaptation

Understanding how selection shapes population differentiation and local adaptation in marine species remains one of the greatest challenges in the field of evolutionary biology. The selection of genes in response to environment-specific factors and microenvironmental variation often results in chaotic genetic patchiness, which is commonly observed in rocky shore organisms. To identify these genes, the expression profile of the marine gastropod Littoraria flava collected from four Southeast Brazilian locations in ten rocky shore sites was analyzed. In this first L. flava transcriptome, 250,641 unigenes were generated, and 24% returned hits after functional annotation. Independent paired comparisons between 1) transects, 2) sites within transects, and 3) sites from different transects were performed for differential expression, detecting 8,622 unique differentially expressed genes. Araçá (AR) and São João (SJ) transect comparisons showed the most divergent gene products. For local adaptation, fitness-related differentially expressed genes were chosen for selection tests. Nine and 24 genes under adaptative and purifying selection, respectively, were most related to biomineralization in AR and chaperones in SJ. The biomineralization-genes perlucin and gigasin-6 were positively selected exclusively in the site toward the open ocean in AR, with sequence variants leading to pronounced protein structure changes. Despite an intense gene flow among L. flava populations due to its planktonic larva, gene expression patterns within transects may be the result of selective pressures. Our findings represent the first step in understanding how microenvironmental genetic variation is maintained in rocky shore populations and the mechanisms underlying local adaptation in marine species.

The population genomics of repeated freshwater colonizations by Gulf pipefish

How organisms adapt to the novel challenges imposed by the colonization of a new habitat has long been a central question in evolutionary biology. When multiple populations of the same species independently adapt to similar environmental challenges, the question becomes whether the populations have arrived at their adaptations through the same genetic mechanisms. In recent years, genetic techniques have been used to tackle these questions by investigating the genome-level changes underlying local adaptation. Here, we present a genomic analysis of colonization of freshwater habitats by a primarily marine fish, the Gulf pipefish (Syngnathus scovelli). We sample pipefish from four geographically distinct freshwater locations and use double-digest restriction site associated DNA sequencing to compare them to 12 previously studied saltwater populations. The two most geographically distant and isolated freshwater populations are the most genetically distinct, although demographic analysis suggests that these populations are experiencing ongoing migration with their saltwater neighbours. Additionally, outlier regions were found genome-wide, showing parallelism across ecotype pairs. We conclude that these multiple freshwater colonizations involve similar genomic regions, despite the large geographical distances and different underlying mechanisms. These similar patterns are probably facilitated by the interacting effects of intrinsic barriers, gene flow among populations and ecological selection in the Gulf pipefish.

High-quality carnivoran genomes from roadkill samples enable comparative species delineation in aardwolf and bat-eared fox

In a context of ongoing biodiversity erosion, obtaining genomic resources from wildlife is essential for conservation. The thousands of yearly mammalian roadkill provide a useful source material for genomic surveys. To illustrate the potential of this underexploited resource, we used roadkill samples to study the genomic diversity of the bat-eared fox (Otocyon megalotis) and the aardwolf (Proteles cristatus), both having subspecies with similar disjunct distributions in Eastern and Southern Africa. First, we obtained reference genomes with high contiguity and gene completeness by combining Nanopore long reads and Illumina short reads. Then, we showed that the two subspecies of aardwolf might warrant species status (P. cristatus and P. septentrionalis) by comparing their genome-wide genetic differentiation to pairs of well-defined species across Carnivora with a new Genetic Differentiation index (GDI) based on only a few resequenced individuals. Finally, we obtained a genome-scale Carnivora phylogeny including the new aardwolf species.

Performing parentage analysis for polysomic inheritances based on allelic phenotypes

Polyploidy poses several problems for parentage analysis. We present a new polysomic inheritance model for parentage analysis based on genotypes or allelic phenotypes to solve these problems. The effects of five factors are simultaneously accommodated in this model: (1) double-reduction, (2) null alleles, (3) negative amplification, (4) genotyping errors and (5) self-fertilization. To solve genotyping ambiguity (unknown allele dosage), we developed a new method to establish the likelihood formulas for allelic phenotype data and to simultaneously include the effects of our five chosen factors. We then evaluated and compared the performance of our new method with three established methods by using both simulated data and empirical data from the cultivated blueberry (Vaccinium corymbosum). We also developed and compared the performance of two additional estimators to estimate the genotyping error rate and the sample rate. We make our new methods freely available in the software package polygene, at http://github.com/huangkang1987/polygene.

The genomic signature of ecological divergence along the benthic-limnetic axis in allopatric and sympatric threespine stickleback

The repeated occurrence of similar phenotypes in independent lineages (i.e., parallel evolution) in response to similar ecological conditions can provide compelling insights into the process of adaptive evolution. An intriguing question is to what extent repeated phenotypic changes are underlain by repeated changes at the genomic level and whether patterns of genomic divergence differ with the geographic context in which populations evolve. Here, we combined genomic, morphological and ecological data sets to investigate the genomic signatures of divergence across populations of threespine stickleback (Gasterosteus aculeatus) that adapted to contrasting ecological niches (benthic or limnetic) in either sympatry or allopatry. We found that genome-wide differentiation (F_ST ) was an order of magnitude higher and substantially more repeatable for sympatric benthic and limnetic specialists compared to allopatric populations with similar levels of ecological divergence. We identified genomic regions consistently differentiated between sympatric ecotypes that were also differentiated between or associated with benthic vs. limnetic niche in allopatric populations. These candidate regions were enriched on three chromosomes known to be involved in the benthic-limnetic divergence of threespine stickleback. Some candidate regions overlapped with QTL for body shape and trophic traits such as gill raker number, traits that strongly differ between benthic and limnetic ecotypes. In summary, our study shows that magnitude and repeatability of genomic signatures of ecological divergence in threespine stickleback highly depend on the geographic context. The identified candidate regions provide starting points to identify functionally important genes for the adaptation to benthic and limnetic niches.

Parallel Alpine Differentiation in Arabidopsis arenosa

Parallel evolution provides powerful natural experiments for studying repeatability of evolution and genomic basis of adaptation. Well-documented examples from plants are, however, still rare, as are inquiries of mechanisms driving convergence in some traits while divergence in others. Arabidopsis arenosa, a predominantly foothill species with scattered morphologically distinct alpine occurrences is a promising candidate. Yet, the hypothesis of parallelism remained untested. We sampled foothill and alpine populations in all regions known to harbor the alpine ecotype and used SNP genotyping to test for repeated alpine colonization. Then, we combined field surveys and a common garden experiment to quantify phenotypic parallelism. Genetic clustering by region but not elevation and coalescent simulations demonstrated parallel origin of alpine ecotype in four mountain regions. Alpine populations exhibited parallelism in height and floral traits which persisted after two generations in cultivation. In contrast, leaf traits were distinctive only in certain region(s), reflecting a mixture of plasticity and genetically determined non-parallelism. We demonstrate varying degrees and causes of parallelism and non-parallelism across populations and traits within a plant species. Parallel divergence along a sharp elevation gradient makes A. arenosa a promising candidate for studying genomic basis of adaptation.

Putative chromosomal rearrangements are associated primarily with ecotype divergence rather than geographic separation in an intertidal, poorly dispersing snail

Littorina saxatilis is becoming a model system for understanding the genomic basis of ecological speciation. The parallel formation of crab-adapted ecotypes that exhibit partial reproductive isolation from wave-adapted ecotypes has enabled genomic investigation of conspicuous shell traits. Recent genomic studies suggest that chromosomal rearrangements may enable ecotype divergence by reducing gene flow. However, the genomic architecture of traits that are divergent between ecotypes remains poorly understood. Here, we use 11,504 single nucleotide polymorphism (SNP) markers called using the recently released L. saxatilis genome to genotype 462 crab ecotype, wave ecotype and phenotypically intermediate Spanish L. saxatilis individuals with scored phenotypes. We used redundancy analysis to study the genetic architecture of loci associated with shell shape, shape corrected for size, shell size and shell ornamentation, and to compare levels of co-association among different traits. We discovered 341 SNPs associated with shell traits. Loci associated with trait divergence between ecotypes were often located inside putative chromosomal rearrangements recently characterized in Swedish L. saxatilis. In contrast, we found that shell shape corrected for size varied primarily by geographic site rather than by ecotype and showed little association with these putative rearrangements. We conclude that genomic regions of elevated divergence inside putative rearrangements were associated with divergence of L. saxatilis ecotypes along steep environmental axes-consistent with models of adaptation with gene flow-but were not associated with divergence among the three geographical sites. Our findings support predictions from models indicating the importance of genomic regions of reduced recombination allowing co-association of loci during ecological speciation with ongoing gene flow.

Multi-model seascape genomics identifies distinct environmental drivers of selection among sympatric marine species

BACKGROUND

As global change and anthropogenic pressures continue to increase, conservation and management increasingly needs to consider species' potential to adapt to novel environmental conditions. Therefore, it is imperative to characterise the main selective forces acting on ecosystems, and how these may influence the evolutionary potential of populations and species. Using a multi-model seascape genomics approach, we compare putative environmental drivers of selection in three sympatric southern African marine invertebrates with contrasting ecology and life histories: Cape urchin (Parechinus angulosus), Common shore crab (Cyclograpsus punctatus), and Granular limpet (Scutellastra granularis).

RESULTS

Using pooled (Pool-seq), restriction-site associated DNA sequencing (RAD-seq), and seven outlier detection methods, we characterise genomic variation between populations along a strong biogeographical gradient. Of the three species, only S. granularis showed significant isolation-by-distance, and isolation-by-environment driven by sea surface temperatures (SST). In contrast, sea surface salinity (SSS) and range in air temperature correlated more strongly with genomic variation in C. punctatus and P. angulosus. Differences were also found in genomic structuring between the three species, with outlier loci contributing to two clusters in the East and West Coasts for S. granularis and P. angulosus, but not for C. punctatus.

CONCLUSION

The findings illustrate distinct evolutionary potential across species, suggesting that species-specific habitat requirements and responses to environmental stresses may be better predictors of evolutionary patterns than the strong environmental gradients within the region. We also found large discrepancies between outlier detection methodologies, and thus offer a novel multi-model approach to identifying the principal environmental selection forces acting on species. Overall, this work highlights how adding a comparative approach to seascape genomics (both with multiple models and species) can elucidate the intricate evolutionary responses of ecosystems to global change.

Tracking parallel adaptation of shell morphology through geological times in the land snail genus Pupilla (Gastropoda: Stylommatophora: Pupillidae)

Changing environmental conditions force species either to disperse or to adapt locally either genetically or via phenotypic plasticity. Although limits of plasticity can be experimentally tested, the predictability of genetic adaptation is restricted due to its stochastic nature. Nevertheless, our understanding of evolutionary adaptation has been improving in particular through studies of parallel adaptation. Based on molecular phylogenetic inferences and morphological investigations of both recent and fossil shells we tracked the morphological changes in three land snails, Pupilla alpicola, Pupilla loessica and Pupilla muscorum. These species differ in habitat requirements as well as historical and extant distributions with P. alpicola and P. loessica being more similar to each other than to P. muscorum. Therefore, we hypothesized, that the three species reacted independently and individually to the conditions changing throughout the Pleistocene, but expected that changes within P. alpicola and P. loessica would be more similar compared to P. muscorum. Indeed, intraspecific shell shape differences across time were similar in P. alpicola and P. loessica, suggesting that similar niche shifts have led to similar transformations in parallel. In contrast, extant P. muscorum populations were practically identical in shape to their ancestors. They have probably tracked their ecological niches through time.

Assortative mating, sexual selection, and their consequences for gene flow in Littorina

When divergent populations are connected by gene flow, the establishment of complete reproductive isolation usually requires the joint action of multiple barrier effects. One example where multiple barrier effects are coupled consists of a single trait that is under divergent natural selection and also mediates assortative mating. Such multiple-effect traits can strongly reduce gene flow. However, there are few cases where patterns of assortative mating have been described quantitatively and their impact on gene flow has been determined. Two ecotypes of the coastal marine snail, Littorina saxatilis, occur in North Atlantic rocky-shore habitats dominated by either crab predation or wave action. There is evidence for divergent natural selection acting on size, and size-assortative mating has previously been documented. Here, we analyze the mating pattern in L. saxatilis with respect to size in intensively sampled transects across boundaries between the habitats. We show that the mating pattern is mostly conserved between ecotypes and that it generates both assortment and directional sexual selection for small male size. Using simulations, we show that the mating pattern can contribute to reproductive isolation between ecotypes but the barrier to gene flow is likely strengthened more by sexual selection than by assortment.

Genetic Characterization of Cupped Oyster Resources in Europe Using Informative Single Nucleotide Polymorphism (SNP) Panels

The Pacific oyster, Crassostrea gigas, was voluntarily introduced from Japan and British Columbia into Europe in the early 1970s, mainly to replace the Portuguese oyster, Crassostrea angulata, in the French shellfish industry, following a severe disease outbreak. Since then, the two species have been in contact in southern Europe and, therefore, have the potential to exchange genes. Recent evolutionary genomic works have provided empirical evidence that C. gigas and C. angulata exhibit partial reproductive isolation. Although hybridization occurs in nature, the rate of interspecific gene flow varies across the genome, resulting in highly heterogeneous genome divergence. Taking this biological property into account is important to characterize genetic ancestry and population structure in oysters. Here, we identified a subset of ancestry-informative makers from the most differentiated regions of the genome using existing genomic resources. We developed two different panels in order to (i) easily differentiate C. gigas and C. angulata, and (ii) describe the genetic diversity and structure of the cupped oyster with a particular focus on French Atlantic populations. Our results confirm high genetic homogeneity among Pacific cupped oyster populations in France and reveal several cases of introgressions between Portuguese and Japanese oysters in France and Portugal.

Parallelism in eco-morphology and gene expression despite variable evolutionary and genomic backgrounds in a Holarctic fish

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds.

Comparing the Performance of Microsatellites and RADseq in Population Genetic Studies: Analysis of Data for Pike (Esox lucius) and a Synthesis of Previous Studies

Population genetic studies reveal biodiversity patterns and inform about drivers of evolutionary differentiation and adaptation, including gene flow, drift and selection. This can advance our understanding and aid decision making regarding management and conservation efforts. Microsatellites have long been used in population genetic studies. Thanks to the development of newer techniques, sequencing approaches such as restriction site associated DNA sequencing (RADseq) are on their way to replace microsatellites for some applications. However, the performance of these two marker types in population genetics have rarely been systematically compared. We utilized three neutrally and adaptively differentiated populations of anadromous pike (Esox lucius) to assess the relative performance of microsatellites and RADseq with respect to resolution and conclusiveness of estimates of population differentiation and genetic structure. To this end, the same set of individuals (N = 64) were genotyped with both RADseq and microsatellite markers. To assess effects of sample size, the same subset of 10 randomly chosen individuals from each population (N = 30 in total) were also genotyped with both methods. Comparisons of estimated genetic diversity and structure showed that both markers were able to uncover genetic structuring. The full RADseq dataset provided the clearest detection of the finer scaled genetic structuring, and the other three datasets (full and subset microsatellite, and subset RADseq) provided comparable results. A search for outlier loci performed on the full SNP dataset pointed to signs of selection potentially associated with salinity and temperature, exemplifying the utility of RADseq to inform about the importance of different environmental factors. To evaluate whether performance differences between the markers are general or context specific, the results of previous studies that have investigated population structure using both marker types were synthesized. The synthesis revealed that RADseq performed as well as, or better than microsatellites in detecting genetic structuring in the included studies. The differences in the ability to detect population structure, both in the present and the previous studies, are likely explained by the higher number of loci typically utilized in RADseq compared to microsatellite analysis, as increasing the number of markers will (regardless of the marker type) increase power and allow for clearer detection and higher resolution of genetic structure.

Phenology-dependent cold exposure and thermal performance of Ostrinia nubilalis ecotypes

Background

Understanding adaptation involves establishing connections between selective agents and beneficial population responses. However, relatively little attention has been paid to seasonal adaptation, in part, because it requires complex and integrative knowledge about seasonally fluctuating environmental factors, the effects of variable phenology on exposure to those factors, and evidence for temporal specialization. In the European corn borer moth, Ostrinia nubilalis, sympatric pheromone strains exploit the same host plant (Zea mays) but may genetically differ in phenology and be reproductively “isolated by time.” Z strain populations in eastern North America have been shown to have a prolonged larval diapause and produce one annual mating flight (July), whereas E strain populations complete an earlier (June) and a later (August) mating flight by shortening diapause duration. Here, we find evidence consistent with seasonal “adaptation by time” between these ecotypes.

Results

We use 12 years of field observation of adult seasonal abundance to estimate phenology of ecotype life cycles and to quantify life-stage specific climatic conditions. We find that the observed reduction of diapause duration in the E strain leads their non-diapausing, active life stages to experience a ~ 4 °C colder environment compared to the equivalent life stages in the Z strain. For a representative pair of populations under controlled laboratory conditions, we compare life-stage specific cold tolerance and find non-diapausing, active life stages in the E strain have as much as a 60% greater capacity to survive rapid cold shock. Enhanced cold hardiness appears unrelated to life-stage specific changes in the temperature at which tissues freeze.

Conclusions

Our results suggest that isolation by time and adaptation by time may both contribute to population divergence, and they argue for expanded study in this species of allochronic populations in nature experiencing the full spectrum of seasonal environments. Cyclical selective pressures are inherent properties of seasonal habitats. Diverse fluctuating selective agents across each year (temperature, predation, competition, precipitation, etc.) may therefore be underappreciated drivers of biological diversity.

The evolution of polymorphism in the warning coloration of the Amazonian poison frog Adelphobates galactonotus

While intraspecific variation in aposematic signals can be selected for by different predatory responses, their evolution is also contingent on other processes shaping genetic variation. We evaluate the relative contributions of selection, geographic isolation, and random genetic drift to the evolution of aposematic color polymorphism in the poison frog Adelphobates galactonotus, distributed throughout eastern Brazilian Amazonia. Dorsal coloration was measured for 111 individuals and genetic data were obtained from 220 individuals at two mitochondrial genes (mtDNA) and 7963 Single Nucleotide Polymorphisms (SNPs). Four color categories were described (brown, blue, yellow, orange) and our models of frog and bird visual systems indicated that each color was distinguishable for these taxa. Using outlier and correlative analyses we found no compelling genetic evidence for color being under divergent selection. A time-calibrated mtDNA tree suggests that the present distribution of dorsal coloration resulted from processes occurring during the Pleistocene. Separate phylogenies based on SNPs and mtDNA resolved the same well supported clades, each containing different colored populations. Ancestral character state analysis provided some evidence for evolutionary transitions in color type. Genetic structure was more strongly associated with geographic features, than color category, suggesting that the distribution of color is explained by localized processes. Evidence for geographic isolation together with estimates of low effective population size implicates drift as playing a key role in color diversification. Our results highlight the relevance of considering the neutral processes involved with the evolution of traits with important fitness consequences.

A developmentally descriptive method for quantifying shape in gastropod shells

The growth of snail shells can be described by simple mathematical rules. Variation in a few parameters can explain much of the diversity of shell shapes seen in nature. However, empirical studies of gastropod shell shape variation typically use geometric morphometric approaches, which do not capture this growth pattern. We have developed a way to infer a set of developmentally descriptive shape parameters based on three-dimensional logarithmic helicospiral growth and using landmarks from two-dimensional shell images as input. We demonstrate the utility of this approach, and compare it to the geometric morphometric approach, using a large set of Littorina saxatilis shells in which locally adapted populations differ in shape. Our method can be modified easily to make it applicable to a wide range of shell forms, which would allow for investigations of the similarities and differences between and within many different species of gastropods.

Contingent Convergence: The Ability To Detect Convergent Genomic Evolution Is Dependent on Population Size and Migration

Outlier scans, in which the genome is scanned for signatures of selection, have become a prominent tool in studies of local adaptation, and more recently studies of genetic convergence in natural populations. However, such methods have the potential to be confounded by features of demographic history, such as population size and migration, which are considerably varied across natural populations. In this study, we use forward-simulations to investigate and illustrate how several measures of genetic differentiation commonly used in outlier scans (FST, DXY and Δπ) are influenced by demographic variation across multiple sampling generations. In a factorial design with 16 treatments, we manipulate the presence/absence of founding bottlenecks (N of founding individuals), prolonged bottlenecks (proportional size of diverging population) and migration rate between two populations with ancestral and diverged phenotypic optima. Our results illustrate known constraints of individual measures associated with reduced population size and a lack of migration; but notably we demonstrate how relationships between measures are similarly dependent on these features of demography. We find that false-positive signals of convergent evolution (the same simulated outliers detected in independent treatments) are attainable as a product of similar population size and migration treatments (particularly for DXY), and that outliers across different measures (for e.g., FST and DXY) can occur with little influence of selection. Taken together, we show how underappreciated, yet quantifiable measures of demographic history can influence commonly employed methods for detecting selection.

Is embryo abortion a post-zygotic barrier to gene flow between Littorina ecotypes?

Genetic incompatibilities contribute to reproductive isolation between many diverging populations, but it is still unclear to what extent they play a role if divergence happens with gene flow. In contact zones between the "Crab" and "Wave" ecotypes of the snail Littorina saxatilis, divergent selection forms strong barriers to gene flow, while the role of post‐zygotic barriers due to selection against hybrids remains unclear. High embryo abortion rates in this species could indicate the presence of such barriers. Post‐zygotic barriers might include genetic incompatibilities (e.g. Dobzhansky–Muller incompatibilities) but also maladaptation, both expected to be most pronounced in contact zones. In addition, embryo abortion might reflect physiological stress on females and embryos independent of any genetic stress. We examined all embryos of >500 females sampled outside and inside contact zones of three populations in Sweden. Females' clutch size ranged from 0 to 1,011 embryos (mean 130 ± 123), and abortion rates varied between 0% and 100% (mean 12%). We described female genotypes by using a hybrid index based on hundreds of SNPs differentiated between ecotypes with which we characterized female genotypes. We also calculated female SNP heterozygosity and inversion karyotype. Clutch size did not vary with female hybrid index, and abortion rates were only weakly related to hybrid index in two sites but not at all in a third site. No additional variation in abortion rate was explained by female SNP heterozygosity, but increased female inversion heterozygosity added slightly to increased abortion. Our results show only weak and probably biologically insignificant post‐zygotic barriers contributing to ecotype divergence, and the high and variable abortion rates were marginally, if at all, explained by hybrid index of females.

Cryptic ecological and geographic diversification in coral-associated nudibranchs

Coral reefs are among the most biologically diverse ecosystems of the world, yet little is known about the processes creating and maintaining their diversity. Ecologically, corallivory in nudibranchs resembles phytophagy in insects- a process that for decades has served as a model for ecological speciation via host shifting. This study uses extensive field collections, DNA sequencing, and phylogenetic analyses to reconstruct the evolutionary history of coral-associated nudibranchs and assess the relative roles that host shifting and geography may have played in their diversification. We find that the number of species is three times higher than the number previously known to science, with evidence for both allopatric and ecological divergence through host shifting and host specialization. Results contribute to growing support for the importance of ecological diversification in marine environments and provide evidence for new species in the genus Tenellia.

Genomic architecture of parallel ecological divergence: Beyond a single environmental contrast

The study of parallel ecological divergence provides important clues to the operation of natural selection. Parallel divergence often occurs in heterogeneous environments with different kinds of environmental gradients in different locations, but the genomic basis underlying this process is unknown. We investigated the genomics of rapid parallel adaptation in the marine snail Littorina saxatilis in response to two independent environmental axes (crab-predation versus wave-action and low-shore versus high-shore). Using pooled whole-genome resequencing, we show that sharing of genomic regions of high differentiation between environments is generally low but increases at smaller spatial scales. We identify different shared genomic regions of divergence for each environmental axis and show that most of these regions overlap with candidate chromosomal inversions. Several inversion regions are divergent and polymorphic across many localities. We argue that chromosomal inversions could store shared variation that fuels rapid parallel adaptation to heterogeneous environments, possibly as balanced polymorphism shared by adaptive gene flow.