High degree of cryptic population differentiation in the Baltic Sea herring Clupea harengus

Genome-wide SNPs reveal novel patterns of spatial genetic structure in Aedes albopictus (Diptera Culicidae) population in China

Introduction

Since the second half of the 20th century, Aedes albopictus, a vector for more than 20 arboviruses, has spread worldwide. Aedes albopictus is the main vector of infectious diseases transmitted by Aedes mosquitoes in China, and it has caused concerns regarding public health. A comprehensive understanding of the spatial genetic structure of this vector species at a genomic level is essential for effective vector control and the prevention of vector-borne diseases.

Methods

During 2016-2018, adult female Ae. albopictus mosquitoes were collected from eight different geographical locations across China. Restriction site-associated DNA sequencing (RAD-seq) was used for high-throughput identification of single nucleotide polymorphisms (SNPs) and genotyping of the Ae. albopictus population. The spatial genetic structure was analyzed and compared to those exhibited by mitochondrial cytochrome c oxidase subunit 1 (cox1) and microsatellites in the Ae. albopictus population.

Results

A total of 9,103 genome-wide SNP loci in 101 specimens and 32 haplotypes of cox1 in 231 specimens were identified in the samples from eight locations in China. Principal component analysis revealed that samples from Lingshui and Zhanjiang were more genetically different than those from the other locations. The SNPs provided a better resolution and stronger signals for novel spatial population genetic structures than those from the cox1 data and a set of previously genotyped microsatellites. The fixation indexes from the SNP dataset showed shallow but significant genetic differentiation in the population. The Mantel test indicated a positive correlation between genetic distance and geographical distance. However, the asymmetric gene flow was detected among the populations, and it was higher from south to north and west to east than in the opposite directions.

Conclusions

The genome-wide SNPs revealed seven gene pools and fine spatial genetic structure of the Ae. albopictus population in China. The RAD-seq approach has great potential to increase our understanding of the spatial dynamics of population spread and establishment, which will help us to design new strategies for controlling vectors and mosquito-borne diseases.

Atlantic herring (Clupea harengus) population structure in the Northeast Atlantic Ocean

The Atlantic herring Clupea harengus L has a vast geographical distribution and a complex population structure with a few very large migratory units and many small local populations. Each population has its own spawning ground and/or time, thereby maintaining their genetic integrity. Several herring populations migrate between common feeding grounds and over-wintering areas resulting in frequent mixing of populations. Thus, many herring fisheries are based on mixed populations of different demographic status. In order to avoid over-exploitation of weak populations and to conserve biodiversity, understanding the population structure and population mixing is important for maintaining biologically sustainable herring fisheries. The aim of this study was to investigate the genetic population structure of herring in the Faroese and surrounding waters, and to develop genetic markers for distinguishing between four herring management units (often called stocks), namely the Norwegian spring-spawning herring (NSSH), Icelandic summer-spawning herring (ISSH), North Sea autumn-spawning herring (NSAH), and Faroese autumn-spawning herring (FASH). Herring from the four stocks were sequenced at low coverage, and single nucleotide polymorphisms (SNPs) were called and used for population structure analysis and individual assignment. An ancestry-informative SNP panel with 118 SNPs was developed and tested on 240 individuals. The results showed that all four stocks appeared to be genetically differentiated populations, but at lower levels of differentiation between FASH and ISSH than the other two populations. Overall assignment rate with the SNP panel was 80.7%, and agreement between the genetic and traditional visual assignment was 75.5%. The NSAH and NSSH samples had the highest assignment rate (100% and 98.3%, respectively) and highest agreement between traditional and genetic assignment methods (96.6% and 94.9%, respectively). The FASH and ISSH samples had substantially lower assignment rates (72.9% and 51.7%, respectively) and agreement between traditional and genetic methods (39.5% and 48.4%, respectively).

Genome-wide patterns of genetic diversity, population structure and demographic history in mānuka (Leptospermum scoparium) growing on indigenous Māori land

Leptospermum scoparium J. R. Forst et G. Forst, known as mānuka by Māori, the indigenous people of Aotearoa (New Zealand), is a culturally and economically significant shrub species, native to New Zealand and Australia. Chemical, morphological and phylogenetic studies have indicated geographical variation of mānuka across its range in New Zealand, and genetic differentiation between New Zealand and Australia. We used pooled whole genome re-sequencing of 76 L. scoparium and outgroup populations from New Zealand and Australia to compile a dataset totalling ~2.5 million SNPs. We explored the genetic structure and relatedness of L. scoparium across New Zealand, and between populations in New Zealand and Australia, as well as the complex demographic history of this species. Our population genomic investigation suggests there are five geographically distinct mānuka gene pools within New Zealand, with evidence of gene flow occurring between these pools. Demographic modelling suggests three of these gene pools have undergone expansion events, whilst the evolutionary histories of the remaining two have been subjected to contractions. Furthermore, mānuka populations in New Zealand are genetically distinct from populations in Australia, with coalescent modelling suggesting these two clades diverged ~9-12 million years ago. We discuss the evolutionary history of this species and the benefits of using pool-seq for such studies. Our research will support the management and conservation of mānuka by landowners, particularly Māori, and the development of a provenance story for the branding of mānuka based products.

Population structure, genomic diversity and demographic history of Komodo dragons inferred from whole-genome sequencing

Population and conservation genetics studies have greatly benefited from the development of new techniques and bioinformatic tools associated with next-generation sequencing. Analysis of extensive data sets from whole-genome sequencing of even a few individuals allows the detection of patterns of fine-scale population structure and detailed reconstruction of demographic dynamics through time. In this study, we investigated the population structure, genomic diversity and demographic history of the Komodo dragon (Varanus komodoensis), the world's largest lizard, by sequencing the whole genomes of 24 individuals from the five main Indonesian islands comprising the entire range of the species. Three main genomic groups were observed. The populations of the Island of Komodo and the northern coast of Flores, in particular, were identified as two distinct conservation units. Degrees of genomic divergence among island populations were interpreted as a result of changes in sea level affecting connectivity across islands. Demographic inference suggested that Komodo dragons probably experienced a relatively steep population decline over the last million years, reaching a relatively stable N_e during the Saalian glacial cycle (400-150 thousand years ago) followed by a rapid N_e decrease. Genomic diversity of Komodo dragons was similar to that found in endangered or already extinct reptile species. Overall, this study provides an example of how whole-genome analysis of a few individuals per population can help define population structure and intraspecific demographic dynamics. This is particularly important when applying population genomics data to conservation of rare or elusive endangered species.

Re-examination of population structure in Arctic ringed seals using DArTseq genotyping

Although Arctic ringed seals Phoca hispida hispida are currently abundant and broadly distributed, their numbers are projected to decline substantially by the year 2100 due to climate warming. While understanding population structure could provide insight into the impact of environmental changes on this subspecies, detecting demographically important levels of exchange can be difficult in taxa with high abundance. We used a next-generation sequencing approach (DArTseq) to genotype ~5700 single nucleotide polymorphisms in 79 seals from 4 Pacific Arctic regions. Comparison of the 2 most geographically separated strata (eastern Bering vs. northeastern Chukchi-Beaufort Seas) revealed a statistically significant level of genetic differentiation (FST = 0.001, p = 0.005) that, while small, was 1 to 2 orders of magnitude greater than expected based on divergence estimated for similarly sized populations connected by low (1% yr-1) dispersal. A relatively high proportion (72 to 88%) of individuals within these strata could be genetically assigned to their stratum of origin. These results indicate that demographically important structure may be present among Arctic ringed seals breeding in different areas, increasing the risk that declines in the number of seals breeding in areas most negatively affected by environmental warming could occur.

Mycelial biomass and concentration of loline alkaloids driven by complex population structure in Epichloë uncinata and meadow fescue (Schedonorus pratensis)

Many efforts have been made to select and isolate naturally occurring animal-friendly Epichloë strains for later reinfection into elite cultivars. Often this process involves large-scale screening of Epichloë-infected wild grass populations where strains are characterized and alkaloids measured. Here, we describe for the first time the use of genotyping-by-sequencing (GBS) on a collection of 217 Epichloë-infected grasses (7 S. arundinaceum, 4 L. perenne, and 206 S. pratensis). This genotyping strategy is cheaper than complete genome sequencing, is suitable for a large number of individuals, and, when applied to endophyte-infected grasses, conveniently genotypes both organisms. In total, 6273 single nucleotide polymorphisms (SNPs) in the endophyte data set and 38 323 SNPs in the host data set were obtained. Our findings reveal a composite structure with three distinct endophyte clusters unrelated to the three main S. pratensis gene pools that have most likely spread from different glacial refugia in Eurasia. All three gene pools can establish symbiosis with E. uncinata. A comparison of the endophyte clusters with microsatellite-based fingerprinting of the same samples allows a quick test to discriminate between these clusters using two simple sequence repeats (SSRs). Concentrations of loline alkaloids and mycelial biomass are correlated and differ significantly among the plant and endophyte subpopulations; one endophyte strain has higher levels of lolines than others, and one specific host genotype is particularly suitable to host E. uncinata. These findings pave the way for targeted artificial inoculations of specific host-endophyte combinations to boost loline production in the symbiota and for genome association studies with the aim of isolating genes involved in the compatibility between meadow fescue and E. uncinata.

RAD genotyping reveals fine-scale population structure and provides evidence for adaptive divergence in a commercially important fish from the northwestern Pacific Ocean

Exploring factors shaping genetic structure of marine fish is challenging due to fewer barriers to gene flow in the ocean. However, genome-wide sequence data can greatly enhance our ability to delineate previously unidentified population structure as well as potential adaptive divergence. The small yellow croaker (Larimichthys polyactis) is a commercially important fish species with high gene flow and its overwintering populations experience heterogeneous environment, suggesting possible population differentiation and adaptive divergence. To delineate patterns of population structure as well as test for signatures of local adaptation, a total of 68,666 quality filtered SNP markers were identified for 80 individuals from four overwintering populations by using restriction site-associated DNA sequencing (RAD-seq). Significant genetic differentiation among overwintering populations from the Central Yellow Sea, the South Yellow Sea and the North East China Sea were detected (Pair-wise F_ST: 0.00036–0.00390), which were consistent with population division of overwintering groups inferred from traditional ecological approaches. In addition, a total of 126 unique SNPs were detected to be significantly associated with environmental parameters (temperature, salinity and turbidity). These candidate SNPs were involved in multiple pathways such as energy metabolism and phagocytosis, suggesting they may play key roles in growth and innate immunity. Our results suggested the existence of hitherto unrecognized cryptic population structure and local adaptation in this high gene flow marine fish and thus gain new insights into the design of management strategies.

Geo-referenced population-specific microsatellite data across American continents, the MacroPopGen Database

Population genetic data from nuclear DNA has yet to be synthesized to allow broad scale comparisons of intraspecific diversity versus species diversity. The MacroPopGen database collates and geo-references vertebrate population genetic data across the Americas from 1,308 nuclear microsatellite DNA studies, 897 species, and 9,090 genetically distinct populations where genetic differentiation (FST) was measured. Caribbean populations were particularly distinguished from North, Central, and South American populations, in having higher differentiation (FST = 0.12 vs. 0.07-0.09) and lower mean numbers of alleles (MNA = 4.11 vs. 4.84-5.54). While mammalian populations had lower MNA (4.86) than anadromous fish, reptiles, amphibians, freshwater fish, and birds (5.34-7.81), mean heterozygosity was largely similar across groups (0.57-0.63). Mean FST was consistently lowest in anadromous fishes (0.06) and birds (0.05) relative to all other groups (0.09-0.11). Significant differences in Family/Genera variance among continental regions or taxonomic groups were also observed. MacroPopGen can be used in many future applications including latitudinal analyses, spatial analyses (e.g. central-margin), taxonomic comparisons, regional assessments of anthropogenic impacts on biodiversity, and conservation of wild populations.

First genomic study on Lake Tanganyika sprat Stolothrissa tanganicae: a lack of population structure calls for integrated management of this important fisheries target species

BACKGROUND

Clupeid fisheries in Lake Tanganyika (East Africa) provide food for millions of people in one of the world's poorest regions. Due to climate change and overfishing, the clupeid stocks of Lake Tanganyika are declining. We investigate the population structure of the Lake Tanganyika sprat Stolothrissa tanganicae, using for the first time a genomic approach on this species. This is an important step towards knowing if the species should be managed separately or as a single stock. Population structure is important for fisheries management, yet understudied for many African freshwater species. We hypothesize that distinct stocks of S. tanganicae could be present due to the large size of the lake (isolation by distance), limnological variation (adaptive evolution), or past separation of the lake (historical subdivision). On the other hand, high mobility of the species and lack of obvious migration barriers might have resulted in a homogenous population.

RESULTS

We performed a population genetic study on wild-caught S. tanganicae through a combination of mitochondrial genotyping (96 individuals) and RAD sequencing (83 individuals). Samples were collected at five locations along a north-south axis of Lake Tanganyika. The mtDNA data had low global FST and, visualised in a haplotype network, did not show phylogeographic structure. RAD sequencing yielded a panel of 3504 SNPs, with low genetic differentiation (FST = 0.0054; 95% CI: 0.0046-0.0066). PCoA, fineRADstructure and global FST suggest a near-panmictic population. Two distinct groups are apparent in these analyses (FST = 0.1338 95% CI: 0.1239,0.1445), which do not correspond to sampling locations. Autocorrelation analysis showed a slight increase in genetic difference with increasing distance. No outlier loci were detected in the RADseq data.

CONCLUSION

Our results show at most very weak geographical structuring of the stock and do not provide evidence for genetic adaptation to historical or environmental differences over a north-south axis. Based on these results, we advise to manage the stock as one population, integrating one management strategy over the four riparian countries. These results are a first comprehensive study on the population structure of these important fisheries target species, and can guide fisheries management.

Temporal changes in genetic diversity and forage yield of perennial ryegrass in monoculture and in combination with red clover in swards

Agricultural grasslands are often cultivated as mixtures of grasses and legumes, and an extensive body of literature is available regarding interspecific interactions, and how these relate to yield and agronomic performance. However, knowledge of the impact of intraspecific diversity on grassland functioning is scarce. We investigated these effects during a 4-year field trial established with perennial ryegrass (Lolium perenne) and red clover (Trifolium pratense). We simulated different levels of intraspecific functional diversity by sowing single cultivars or by combining cultivars with contrasting growth habits, in monospecific or bispecific settings (i.e. perennial ryegrass whether or not in combination with red clover). Replicate field plots were established for seven seed compositions. We determined yield parameters and monitored differences in genetic diversity in the ryegrass component among seed compositions, and temporal changes in the genetic composition and genetic diversity at the within plot level. The composition of cultivars of both species affected the yield and species abundance. In general, the presence of clover had a positive effect on the yield. The cultivar composition of the ryegrass component had a significant effect on the yield, both in monoculture, and in combination with clover. For the genetic analyses, we validated empirically that genotyping-by-sequencing of pooled samples (pool-GBS) is a suitable method for accurate measurement of population allele frequencies, and obtained a dataset of 22,324 SNPs with complete data. We present a method to investigate the temporal dynamics of cultivars in seed mixtures grown under field conditions, and show how cultivar abundances vary during subsequent years. We screened the SNP panel for outlier loci, putatively under selection during the cultivation period, but none were detected.

Mind the gut: genomic insights to population divergence and gut microbial composition of two marine keystone species

BACKGROUND

Deciphering the mechanisms governing population genetic divergence and local adaptation across heterogeneous environments is a central theme in marine ecology and conservation. While population divergence and ecological adaptive potential are classically viewed at the genetic level, it has recently been argued that their microbiomes may also contribute to population genetic divergence. We explored whether this might be plausible along the well-described environmental gradient of the Baltic Sea in two species of sand lance (Ammodytes tobianus and Hyperoplus lanceolatus). Specifically, we assessed both their population genetic and gut microbial composition variation and investigated not only which environmental parameters correlate with the observed variation, but whether host genome also correlates with microbiome variation.

RESULTS

We found a clear genetic structure separating the high-salinity North Sea from the low-salinity Baltic Sea sand lances. The observed genetic divergence was not simply a function of isolation by distance, but correlated with environmental parameters, such as salinity, sea surface temperature, and, in the case of A. tobianus, possibly water microbiota. Furthermore, we detected two distinct genetic groups in Baltic A. tobianus that might represent sympatric spawning types. Investigation of possible drivers of gut microbiome composition variation revealed that host species identity was significantly correlated with the microbial community composition of the gut. A potential influence of host genetic factors on gut microbiome composition was further confirmed by the results of a constrained analysis of principal coordinates. The host genetic component was among the parameters that best explain observed variation in gut microbiome composition.

CONCLUSIONS

Our findings have relevance for the population structure of two commercial species but also provide insights into potentially relevant genomic and microbial factors with regards to sand lance adaptation across the North Sea-Baltic Sea environmental gradient. Furthermore, our findings support the hypothesis that host genetics may play a role in regulating the gut microbiome at both the interspecific and intraspecific levels. As sequencing costs continue to drop, we anticipate that future studies that include full genome and microbiome sequencing will be able to explore the full relationship and its potential adaptive implications for these species.

Assessing the genetic diversity of the critically endangered Chinese sturgeon Acipenser sinensis using mitochondrial markers and genome-wide single-nucleotide polymorphisms from RAD-seq

As a living fossil, the endangered Chinese sturgeon (Acipenser sinensis) has been considered a national treasure in China. Here, the famous Gezhouba Dam and Three Gorges Dam on the Yangtze River were built in 1988 and 2006, for economic purposes. The natural population of Chinese sturgeon has declined since then, as these dams block its migratory route to the original spawning grounds in the middle reaches of the Yangtze River. In 2013 and 2014, there was an absence of spawning where it typically happened near the Gezhouba Dam. Nevertheless, from April to June in 2015, over 1,000 larvae with different body lengths (10-35 cm) were detected along the Shanghai Yangtze Estuary; but only little is currently known about the population genetic structure of the Chinese sturgeon. Herein, we inferred population genetic parameters from 462 available Chinese sturgeon specimens based on a 421-bp fragment of the mitochondrial DNA (mtDNA) D-loop region and 1,481,620 SNPs (single-nucleotide polymorphisms) generated by restriction site-associated DNA sequencing (RAD-seq). For the D-loop dataset, 15 haplotypes were determined. Randomly picked 23 individuals, representing the 15 D-loop haplotype groups, were subsequently used for further RAD-seq validation. The average nucleotide diversity calculated from the mtDNA and RAD datasets was 0.0086 and 0.000478, respectively. The overall effective female population size was calculated to be 1,255 to 2,607, and the long-term effective population size was estimated to range from 11,950 to 119,500. We observed that the genetic variability and the effective female population size of the current population in the Yangtze River are severely low, which are similar to the data reported over 10 years ago. The deduced relatively small effective population of female fish, limiting the genetic connectivity among Chinese sturgeon, should be considered a serious threat to this endangered species.

Oceanographic variation influences spatial genomic structure in the sea scallop, Placopecten magellanicus

Environmental factors can influence diversity and population structure in marine species and accurate understanding of this influence can both improve fisheries management and help predict responses to environmental change. We used 7163 SNPs derived from restriction site-associated DNA sequencing genotyped in 245 individuals of the economically important sea scallop, Placopecten magellanicus, to evaluate the correlations between oceanographic variation and a previously identified latitudinal genomic cline. Sea scallops span a broad latitudinal area (>10 degrees), and we hypothesized that climatic variation significantly drives clinal trends in allele frequency. Using a large environmental dataset, including temperature, salinity, chlorophyll a, and nutrient concentrations, we identified a suite of SNPs (285-621, depending on analysis and environmental dataset) potentially under selection through correlations with environmental variation. Principal components analysis of different outlier SNPs and environmental datasets revealed similar northern and southern clusters, with significant associations between the first axes of each (R2adj = .66-.79). Multivariate redundancy analysis of outlier SNPs and the environmental principal components indicated that environmental factors explained more than 32% of the variance. Similarly, multiple linear regressions and random-forest analysis identified winter average and minimum ocean temperatures as significant parameters in the link between genetic and environmental variation. This work indicates that oceanographic variation is associated with the observed genomic cline in this species and that seasonal periods of extreme cold may restrict gene flow along a latitudinal gradient in this marine benthic bivalve. Incorporating this finding into management may improve accuracy of management strategies and future predictions.

Genomic associations with bill length and disease reveal drift and selection across island bird populations

Island species provide excellent models for investigating how selection and drift operate in wild populations, and for determining how these processes act to influence local adaptation and speciation. Here, we examine the role of selection and drift in shaping genomic and phenotypic variation across recently separated populations of Berthelot's pipit (Anthus berthelotii), a passerine bird endemic to three archipelagos in the Atlantic. We first characterized genetic diversity and population structuring that supported previous inferences of a history of recent colonizations and bottlenecks. We then tested for regions of the genome associated with the ecologically important traits of bill length and malaria infection, both of which vary substantially across populations in this species. We identified a SNP associated with variation in bill length among individuals, islands, and archipelagos; patterns of variation at this SNP suggest that both phenotypic and genotypic variation in bill length is largely shaped by founder effects. Malaria was associated with SNPs near/within genes involved in the immune response, but this relationship was not consistent among archipelagos, supporting the view that disease resistance is complex and rapidly evolving. Although we found little evidence for divergent selection at candidate loci for bill length and malaria resistance, genome scan analyses pointed to several genes related to immunity and metabolism as having important roles in divergence and adaptation. Our findings highlight the utility and challenges involved with combining association mapping and population genetic analysis in nonequilibrium populations, to disentangle the effects of drift and selection on shaping genotypes and phenotypes.

Genetic Diversity and Connectivity in Maurolicus muelleri in the Bay of Biscay Inferred from Thousands of SNP Markers

Mesopelagic fish are largely abundant poorly studied fish that are still intact, but which, due to their potentially great added value, will be imminently exploited by humans. Therefore, studies that provide information to anticipate the anthropogenic impact on this important resource are urgently needed. In particular, knowledge about their connectivity, potential adaptation and resilience are needed. This information can be obtained through the analysis of genome-wide markers which are now relatively easily and cost-efficiently discovered thanks to high-throughput sequencing technologies. Here, we have generated thousands of SNP markers in Maurolicus muelleri, based on the restriction-site associated DNA sequencing method, and preformed population connectivity and genetic diversity analyses in a subset of samples collected from the Bay of Biscay. Our study proves the method valid for obtaining genome-wide markers in this species and provides the first insights into the population genomics of M. muelleri. Importantly, the genomic resources developed here are made available for future studies and set the basics for additional endeavors on this issue.

Population genetics analysis of the Nujiang catfish Creteuchiloglanis macropterus through a genome-wide single nucleotide polymorphisms resource generated by RAD-seq

Advances in genome scanning using high-throughput sequencing technologies has led to a revolution in studies of non-model organisms. The glyptosternoid fish Creteuchiloglanis macropterus, is widely distributed in the main stem and tributaries of the Nujiang River basin. Here, we analyzed IIB restriction-site-associated DNA (2b-RAD) sequences and mitochondrial DNA sequences, to assess the genomic signature of adaptation by detecting and estimating the degree of genetic differentiation among ten Creteuchiloglanis macropterus populations from the Nujiang River. The analyses revealed significant population differentiation among the up-tributaries, main stem, mid-tributary and low-tributary. Annotation of contigs containing outlier SNPs revealed that the candidate genes showed significant enrichment in several important biological process terms between up-tributaries and low-tributary, and exhibited prominent enrichment in the term macromolecular metabolic process between all tributaries and the main stem. Population dynamics analyses indicated that the Late Pleistocene glaciations strongly influenced the demographic history of C. macropterus. Our results provide strong evidence for the utility of RAD-seq in population genetics studies, and our generated SNP resource should provide a valuable tool for population genomics studies of C. macropterus in the future.

Genomic evidence for local adaptation in the ovoviviparous marine fish Sebastiscus marmoratus with a background of population homogeneity

Advances in next-generation sequencing techniques have allowed for the generation of genome-wide sequence data, to gain insight into the dynamics influencing genetic structure and the local adaptation of marine fish. Here, using genotyping-by-sequencing (GBS) technique, we identified 31,119 single nucleotide polymorphisms (SNPs) for Sebastiscus marmoratus in 59 individuals from three populations in Chinese coastal waters. Based on all SNPs, there was little evidence of genetic differentiation among populations. However, outlier tests revealed 329 SNPs putatively under divergent selection across populations. Structural and phylogenetic topology analyses based on the outliers showed clear genetic differentiation among populations. Gene Ontology (GO) annotation results revealed that most of these outliers are known or hypothesized to be involved in metabolic process. Together with previous work using mitochondrial cytochrome b sequences, the present results further suggest that the population structure is strongly influenced by locally adaptive pressure. Overall, adaptive evolution in a heterogeneous environment plays an important role in inducing genetic differentiation among local populations. This study increases understanding of the factors (including gene flow and local adaptation) promoting and constraining population genetic differentiation in marine organisms.

Using a butterflyfish genome as a general tool for RAD-Seq studies in specialized reef fish

Data from a large-scale restriction site-associated DNA sequencing (RAD-Seq) study of nine butterflyfish species in the Red Sea and Arabian Sea provided a means to test the utility of a recently published draft genome (Chaetodon austriacus) and assess apparent bias in this method of isolating nuclear loci. We here processed double-digest restriction site-associated DNA (ddRAD) sequencing data to identify single nucleotide polymorphism (SNP) markers and their associated function with and without our reference genome to see whether it improves the quality of RAD-Seq. Our analyses indicate (i) a modest gap between the number of nonannotated versus annotated SNPs across all species, (ii) an advantage of using genomic resources for closely related but not distantly related butterflyfish species based on the ability to assign putative gene function to SNPs and (iii) an enrichment of genes among sister butterflyfish taxa related to calcium transmembrane transport and binding. The latter result highlights the potential for this approach to reveal insights into adaptive mechanisms in populations inhabiting challenging coral reef environments such as the Red Sea, Arabian Sea and Arabian Gulf with further study.

Fine-scale genetic structure due to adaptive divergence among microhabitats

It has been suggested that adaptive evolution on ecological timescales shapes communities. However, adaptation among environments relies on isolation or large selection coefficients that exceed migration effects. This reliance is tempered if adaptation is polygenic-does not depend on one allele completely replacing another but instead requires small allele frequency changes at many loci. Thus, whether individuals can evolve adaptation to fine-scale habitat variation (for example, microhabitats) is not resolved. Here we analyze the genetic divergence of the teleost fish, Fundulus heteroclitus, among microhabitats that are <200 m apart in three separate saltmarshes using 4741 single-nucleotide polymorphisms (SNPs). Among these SNPs, 1.3-2.3% have large and highly significant differences among microhabitats (mean F_ST=0.15; false discovery rate ⩽1%). The divergence among microhabitats for these outlier SNPs is larger than that among populations, exceeds neutral expectation and indicates surprising population structure among microhabitats. Thus, we suggest that polygenic selection is surprisingly effective in altering allele frequencies among many different SNPs that share similar biological functions in response to environmental and ecological differences over very small geographic distances. We acknowledge the evolutionary difficulty of large genetic divergence among well-connected habitats. Therefore, these studies are only the first step to discern whether natural selection is responsible and capable of effecting genetic divergence on such a fine scale.

Absence of population structure across elevational gradients despite large phenotypic variation in mountain chickadees (Poecile gambeli)

Montane habitats are characterized by predictably rapid heterogeneity along elevational gradients and are useful for investigating the consequences of environmental heterogeneity for local adaptation and population genetic structure. Food-caching mountain chickadees inhabit a continuous elevation gradient in the Sierra Nevada, and birds living at harsher, high elevations have better spatial memory ability and exhibit differences in male song structure and female mate preference compared to birds inhabiting milder, low elevations. While high elevation birds breed, on average, two weeks later than low elevation birds, the extent of gene flow between elevations is unknown. Despite phenotypic variation and indirect evidence for local adaptation, population genetic analyses based on 18 073 single nucleotide polymorphisms across three transects of high and low elevation populations provided no evidence for genetic differentiation. Analyses based on individual genotypes revealed no patterns of clustering, pairwise estimates of genetic differentiation (F_ST, Nei's D) were very low, and AMOVA revealed no evidence for genetic variation structured by transect or by low and high elevation sites within transects. In addition, we found no consistent evidence for strong parallel allele frequency divergence between low and high elevation sites within the three transects. Large elevation-related phenotypic variation may be maintained by strong selection despite gene flow and future work should focus on the mechanisms underlying such variation.

Identifying patterns of dispersal, connectivity and selection in the sea scallop, Placopecten magellanicus, using RADseq-derived SNPs

Understanding patterns of dispersal and connectivity among marine populations can directly inform fisheries conservation and management. Advances in high-throughput sequencing offer new opportunities for estimating marine connectivity. We used restriction-site-associated DNA sequencing to examine dispersal and realized connectivity in the sea scallop Placopecten magellanicus, an economically important marine bivalve. Based on 245 individuals sampled rangewide at 12 locations from Newfoundland to the Mid-Atlantic Bight, we identified and genotyped 7163 single nucleotide polymorphisms; 112 (1.6%) were identified as outliers potentially under directional selection. Bayesian clustering revealed a discontinuity between northern and southern samples, and latitudinal clines in allele frequencies were observed in 42.9% of the outlier loci and in 24.6% of neutral loci. Dispersal estimates derived using these clines and estimates of linkage disequilibrium imply limited dispersal; 373.1 ± 407.0 km (mean ± SD) for outlier loci and 641.0 ± 544.6 km (mean ± SD) for neutral loci. Our analysis suggests restricted dispersal compared to the species range (>2000 km) and that dispersal and effective connectivity differ. These observations support the hypothesis that limited effective dispersal structures scallop populations along eastern North America. These findings can help refine the appropriate scale of management and conservation in this commercially valuable species.

Spatiotemporal heterogeneity in precipitation patterns explain population-level germination strategies in an edaphic specialist

BACKGROUND AND AIMS

Many locally endemic species in biodiversity hotspots are restricted to edaphic conditions that are fixed in the landscape, limiting their potential to track climate change through dispersal. Instead, such species experience strong selection for germination strategies that can track suitable conditions through time. Germination strategies were compared among populations across the geographic range of a California vernal pool annual, Lasthenia fremontii Local germination strategies were tested to determine the associations with geographic variation in precipitation patterns.

METHODS

This study evaluated patterns of seed germination, dormancy and mortality in response to simulated variation in the timing, amount and duration of the first autumn precipitation event using seeds from six populations that span a geographic gradient in precipitation. Next, it was tested whether the germination strategies of different populations can be predicted by historical precipitation patterns that characterize each site.

KEY RESULTS

A significant positive relationship was observed between the historical variability in autumn precipitation and the extent of dormancy in a population. Marginal populations, with histories of the most extreme but constant autumn precipitation levels, expressed the lowest dormancy levels. Populations from sites with historically higher levels of autumn precipitation tended to germinate faster, but this tendency was not statistically significant.

CONCLUSIONS

Germination in L. fremontii is cued by the onset of the first rains that characterize the beginning of winter in California's Great Central Valley. However, populations differ in how fast they germinate and the fraction of seeds that remain dormant when germination cues occur. The results suggest that seed dormancy may be a key trait for populations to track increasingly drier climates predicted by climate change models. However, the low dormancy and high mortality levels observed among seeds of the southernmost, driest populations make them most vulnerable to local extinction.

On the maintenance of genetic variation and adaptation to environmental change: considerations from population genomics in fishes

The first goal of this paper was to overview modern approaches to local adaptation, with a focus on the use of population genomics data to detect signals of natural selection in fishes. Several mechanisms are discussed that may enhance the maintenance of genetic variation and evolutionary potential, which have been overlooked and should be considered in future theoretical development and predictive models: the prevalence of soft sweeps, polygenic basis of adaptation, balancing selection and transient polymorphisms, parallel evolution, as well as epigenetic variation. Research on fish population genomics has provided ample evidence for local adaptation at the genome level. Pervasive adaptive evolution, however, seems to almost never involve the fixation of beneficial alleles. Instead, adaptation apparently proceeds most commonly by soft sweeps entailing shifts in frequencies of alleles being shared between differentially adapted populations. One obvious factor contributing to the maintenance of standing genetic variation in the face of selective pressures is that adaptive phenotypic traits are most often highly polygenic, and consequently the response to selection should derive mostly from allelic co-variances among causative loci rather than pronounced allele frequency changes. Balancing selection in its various forms may also play an important role in maintaining adaptive genetic variation and the evolutionary potential of species to cope with environmental change. A large body of literature on fishes also shows that repeated evolution of adaptive phenotypes is a ubiquitous evolutionary phenomenon that seems to occur most often via different genetic solutions, further adding to the potential options of species to cope with a changing environment. Moreover, a paradox is emerging from recent fish studies whereby populations of highly reduced effective population sizes and impoverished genetic diversity can apparently retain their adaptive potential in some circumstances. Although more empirical support is needed, several recent studies suggest that epigenetic variation could account for this apparent paradox. Therefore, epigenetic variation should be fully integrated with considerations pertaining to role of soft sweeps, polygenic and balancing selection, as well as repeated adaptation involving different genetic basis towards improving models predicting the evolutionary potential of species to cope with a changing world.

Large-Scale Genotyping-by-Sequencing Indicates High Levels of Gene Flow in the Deep-Sea Octocoral Swiftia simplex (Nutting 1909) on the West Coast of the United States

Deep-sea corals are a critical component of habitat in the deep-sea, existing as regional hotspots for biodiversity, and are associated with increased assemblages of fish, including commercially important species. Because sampling these species is so difficult, little is known about the connectivity and life history of deep-sea octocoral populations. This study evaluates the genetic connectivity among 23 individuals of the deep-sea octocoral Swiftia simplex collected from Eastern Pacific waters along the west coast of the United States. We utilized high-throughput restriction-site associated DNA (RAD)-tag sequencing to develop the first molecular genetic resource for the deep-sea octocoral, Swiftia simplex. Using this technique we discovered thousands of putative genome-wide SNPs in this species, and after quality control, successfully genotyped 1,145 SNPs across individuals sampled from California to Washington. These SNPs were used to assess putative population structure across the region. A STRUCTURE analysis as well as a principal coordinates analysis both failed to detect any population differentiation across all geographic areas in these collections. Additionally, after assigning individuals to putative population groups geographically, no significant F_ST values could be detected (F_ST for the full data set 0.0056), and no significant isolation by distance could be detected (p = 0.999). Taken together, these results indicate a high degree of connectivity and potential panmixia in S. simplex along this portion of the continental shelf.

Genomic Methods Take the Plunge: Recent Advances in High-Throughput Sequencing of Marine Mammals

The dramatic increase in the application of genomic techniques to non-model organisms (NMOs) over the past decade has yielded numerous valuable contributions to evolutionary biology and ecology, many of which would not have been possible with traditional genetic markers. We review this recent progression with a particular focus on genomic studies of marine mammals, a group of taxa that represent key macroevolutionary transitions from terrestrial to marine environments and for which available genomic resources have recently undergone notable rapid growth. Genomic studies of NMOs utilize an expanding range of approaches, including whole genome sequencing, restriction site-associated DNA sequencing, array-based sequencing of single nucleotide polymorphisms and target sequence probes (e.g., exomes), and transcriptome sequencing. These approaches generate different types and quantities of data, and many can be applied with limited or no prior genomic resources, thus overcoming one traditional limitation of research on NMOs. Within marine mammals, such studies have thus far yielded significant contributions to the fields of phylogenomics and comparative genomics, as well as enabled investigations of fitness, demography, and population structure. Here we review the primary options for generating genomic data, introduce several emerging techniques, and discuss the suitability of each approach for different applications in the study of NMOs.

Fine-scale genetic structure among greater sage-grouse leks in central Nevada

Background

Mating systems that reduce dispersal and lead to non-random mating might increase the potential for genetic structure to arise at fine geographic scales. Greater sage-grouse (Centrocercus urophasianus) have a lek-based mating system and exhibit high site fidelity and skewed mating ratios. We quantified population structure by analyzing variation at 27,866 single-nucleotide polymorphisms in 140 males from ten leks (within five lek complexes) occurring in a small geographic region in central Nevada.

Results

Lek complexes, and to a lesser extent individual leks, formed statistically identifiable clusters in ordination analyses, providing evidence for fine-scale geographic genetic differentiation. Lek geography predicted genetic differentiation even at a small geographic scale, which could be sharpened by strong site fidelity. Relatedness was also higher among individuals within lek complexes (and leks), suggesting that reproductive skew, where few males participate in most of the successful matings, could also potentially contribute to genetic differentiation. Models incorporating a habitat resistance surface as a proxy for potentially reduced movement due to landscape features indicated that both geographic distance and habitat suitability (i.e. preferred habitat) predicted genetic structure, with no significant effect of man-made barriers to movement (i.e. power lines and roads). Finally, we illustrate how data sets containing fewer loci (<4000) had less statistical precision and failed to detect the full degree of genetic structure.

Conclusion

Our results suggest that habitat features and lek site geography of sage-grouse shape fine scale genetic structure, and highlight how larger data sets can have increased precision and accuracy for quantifying ecologically relevant genetic structure over small geographic scales.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0702-4) contains supplementary material, which is available to authorized users.

Outlier Loci Detect Intraspecific Biodiversity amongst Spring and Autumn Spawning Herring across Local Scales

Herring, Clupea harengus, is one of the ecologically and commercially most important species in European northern seas, where two distinct ecotypes have been described based on spawning time; spring and autumn. To date, it is unknown if these spring and autumn spawning herring constitute genetically distinct units. We assessed levels of genetic divergence between spring and autumn spawning herring in the Baltic Sea using two types of DNA markers, microsatellites and Single Nucleotide Polymorphisms, and compared the results with data for autumn spawning North Sea herring. Temporally replicated analyses reveal clear genetic differences between ecotypes and hence support reproductive isolation. Loci showing non-neutral behaviour, so-called outlier loci, show convergence between autumn spawning herring from demographically disjoint populations, potentially reflecting selective processes associated with autumn spawning ecotypes. The abundance and exploitation of the two ecotypes have varied strongly over space and time in the Baltic Sea, where autumn spawners have faced strong depression for decades. The results therefore have practical implications by highlighting the need for specific management of these co-occurring ecotypes to meet requirements for sustainable exploitation and ensure optimal livelihood for coastal communities.

Three chromosomal rearrangements promote genomic divergence between migratory and stationary ecotypes of Atlantic cod

Identification of genome-wide patterns of divergence provides insight on how genomes are influenced by selection and can reveal the potential for local adaptation in spatially structured populations. In Atlantic cod – historically a major marine resource – Northeast-Arctic- and Norwegian coastal cod are recognized by fundamental differences in migratory and non-migratory behavior, respectively. However, the genomic architecture underlying such behavioral ecotypes is unclear. Here, we have analyzed more than 8.000 polymorphic SNPs distributed throughout all 23 linkage groups and show that loci putatively under selection are localized within three distinct genomic regions, each of several megabases long, covering approximately 4% of the Atlantic cod genome. These regions likely represent genomic inversions. The frequency of these distinct regions differ markedly between the ecotypes, spawning in the vicinity of each other, which contrasts with the low level of divergence in the rest of the genome. The observed patterns strongly suggest that these chromosomal rearrangements are instrumental in local adaptation and separation of Atlantic cod populations, leaving footprints of large genomic regions under selection. Our findings demonstrate the power of using genomic information in further understanding the population dynamics and defining management units in one of the world’s most economically important marine resources.

Development and preliminary evaluation of a genomewide single nucleotide polymorphisms resource generated by RAD-seq for the small yellow croaker (Larimichthys polyactis)

Recent advances in high-throughput sequencing technologies have offered the possibility to generate genomewide sequence data to delineate previously unidentified genetic structure, obtain more accurate estimates of demographic parameters and to evaluate potential adaptive divergence. Here, we identified 27 556 single nucleotide polymorphisms for the small yellow croaker (Larimichthys polyactis) using restriction-site-associated DNA (RAD) sequencing of 24 individuals from two populations. Significant sources of genetic variation were identified, with an average nucleotide diversity (π) of 0.00105 ± 0.000425 across individuals, and long-term effective population size was thus estimated to range between 26 172 and 261 716. According to the results, no differentiation between the two populations was detected based on the SNP data set of top quality score per contig or neutral loci. However, the two analysed populations were highly differentiated based on SNP data set of both top FST value per contig and the outlier SNPs. Moreover, local adaptation was highlighted by an FST -based outlier tests implemented in LOSITAN and a total of 538 potentially locally selected SNPs were identified. blast2go annotation of contigs containing the outlier SNPs yielded hits for 37 (66%) of 56 significant blastx matches. Candidate genes for local adaptation constituted a wide array of biological functions, including cellular response to oxidative stress, actin filament binding, ion transmembrane transport and synapse assembly. The generated SNP resources in this study provided a valuable tool for future population genetics and genomics studies of L. polyactis.

Validation of Pooled Whole-Genome Re-Sequencing in Arabidopsis lyrata

Sequencing pooled DNA of multiple individuals from a population instead of sequencing individuals separately has become popular due to its cost-effectiveness and simple wet-lab protocol, although some criticism of this approach remains. Here we validated a protocol for pooled whole-genome re-sequencing (Pool-seq) of Arabidopsis lyrata libraries prepared with low amounts of DNA (1.6 ng per individual). The validation was based on comparing single nucleotide polymorphism (SNP) frequencies obtained by pooling with those obtained by individual-based Genotyping By Sequencing (GBS). Furthermore, we investigated the effect of sample number, sequencing depth per individual and variant caller on population SNP frequency estimates. For Pool-seq data, we compared frequency estimates from two SNP callers, VarScan and Snape; the former employs a frequentist SNP calling approach while the latter uses a Bayesian approach. Results revealed concordance correlation coefficients well above 0.8, confirming that Pool-seq is a valid method for acquiring population-level SNP frequency data. Higher accuracy was achieved by pooling more samples (25 compared to 14) and working with higher sequencing depth (4.1× per individual compared to 1.4× per individual), which increased the concordance correlation coefficient to 0.955. The Bayesian-based SNP caller produced somewhat higher concordance correlation coefficients, particularly at low sequencing depth. We recommend pooling at least 25 individuals combined with sequencing at a depth of 100× to produce satisfactory frequency estimates for common SNPs (minor allele frequency above 0.05).

Exploring the utility of cross-laboratory RAD-sequencing datasets for phylogenetic analysis

BACKGROUND

Restriction site-Associated DNA sequencing (RAD-Seq) is widely applied to generate genome-wide sequence and genetic marker datasets. RAD-Seq has been extensively utilised, both at the population level and across species, for example in the construction of phylogenetic trees. However, the consistency of RAD-Seq data generated in different laboratories, and the potential use of cross-species orthologous RAD loci in the estimation of genetic relationships, have not been widely investigated. This study describes the use of SbfI RAD-Seq data for the estimation of evolutionary relationships amongst ten teleost fish species, using previously established phylogeny as a benchmark.

RESULTS

The number of orthologous SbfI RAD loci identified decreased with increasing evolutionary distance between the species, with several thousand loci conserved across five salmonid species (divergence ~50 MY), and several hundred conserved across the more distantly related teleost species (divergence ~100-360 MY). The majority (>70%) of loci identified between the more distantly related species were genic in origin, suggesting that the bias of SbfI towards genic regions is useful for identifying distant orthologs. Interspecific single nucleotide variants at each orthologous RAD locus were identified. Evolutionary relationships estimated using concatenated sequences of interspecific variants were congruent with previously published phylogenies, even for distantly (divergence up to ~360 MY) related species.

CONCLUSION

Overall, this study has demonstrated that orthologous SbfI RAD loci can be identified across closely and distantly related species. This has positive implications for the repeatability of SbfI RAD-Seq and its potential to address research questions beyond the scope of the original studies. Furthermore, the concordance in tree topologies and relationships estimated in this study with published teleost phylogenies suggests that similar meta-datasets could be utilised in the prediction of evolutionary relationships across populations and species with readily available RAD-Seq datasets, but for which relationships remain uncharacterised.

Latitudinal Gradient in Otolith Shape among Local Populations of Atlantic Herring (Clupea harengus L.) in Norway

Otolith shape analysis of Atlantic herring (Clupea harengus) in Norwegian waters shows significant differentiation among fjords and a latitudinal gradient along the coast where neighbouring populations are more similar to each other than to those sampled at larger distances. The otolith shape was obtained using quantitative shape analysis, the outlines were transformed with Wavelet and analysed with multivariate methods. The observed morphological differences are likely to reflect environmental differences but indicate low dispersal among the local herring populations. Otolith shape variation suggests also limited exchange between the local populations and their oceanic counterparts, which could be due to differences in spawning behaviour. Herring from the most northerly location (69°N) in Balsfjord, which is genetically more similar to Pacific herring (Clupea pallasii), differed in otolith shape from all the other populations. Our results suggest that the semi-enclosed systems, where the local populations live and breed, are efficient barriers for dispersal. Otolith shape can thus serve as a marker to identify the origin of herring along the coast of Norway.

Adaptation to Low Salinity Promotes Genomic Divergence in Atlantic Cod (Gadus morhua L.)

How genomic selection enables species to adapt to divergent environments is a fundamental question in ecology and evolution. We investigated the genomic signatures of local adaptation in Atlantic cod (Gadus morhua L.) along a natural salinity gradient, ranging from 35‰ in the North Sea to 7‰ within the Baltic Sea. By utilizing a 12 K SNPchip, we simultaneously assessed neutral and adaptive genetic divergence across the Atlantic cod genome. Combining outlier analyses with a landscape genomic approach, we identified a set of directionally selected loci that are strongly correlated with habitat differences in salinity, oxygen, and temperature. Our results show that discrete regions within the Atlantic cod genome are subject to directional selection and associated with adaptation to the local environmental conditions in the Baltic- and the North Sea, indicating divergence hitchhiking and the presence of genomic islands of divergence. We report a suite of outlier single nucleotide polymorphisms within or closely located to genes associated with osmoregulation, as well as genes known to play important roles in the hydration and development of oocytes. These genes are likely to have key functions within a general osmoregulatory framework and are important for the survival of eggs and larvae, contributing to the buildup of reproductive isolation between the low-salinity adapted Baltic cod and the adjacent cod populations. Hence, our data suggest that adaptive responses to the environmental conditions in the Baltic Sea may contribute to a strong and effective reproductive barrier, and that Baltic cod can be viewed as an example of ongoing speciation.

Otolith shape: a population marker for Atlantic herring Clupea harengus

Otolith shape variation of seven Atlantic herring Clupea harengus populations from Canada, the Faroe Islands, Iceland, Ireland, Norway and Scotland, U.K., covering a large area of the species' distribution, was studied in order to see if otolith shape can be used to discriminate between populations. The otolith shape was obtained using quantitative shape analysis, transformed with Wavelet and analysed with multivariate methods. Significant differences were detected among the seven populations, which could be traced to three morphological structures in the otoliths. The differentiation in otolith shape between populations was not only correlated with their spawning time, indicating a strong environmental effect, but could also be due to differing life-history strategies. A model based on the shape differences discriminates with 94% accuracy between Icelandic summer spawners and Norwegian spring spawners, which are known to mix at feeding grounds. This study shows that otolith shape could become an accurate marker for C. harengus population discrimination.

Genetic structure of capelin (Mallotus villosus) in the northwest Atlantic Ocean

Capelin (Mallotus villosus) is a commercially exploited, key forage-fish species found in the boreal waters of the North Pacific and North Atlantic Oceans. We examined the population structure of capelin throughout their range in the Canadian northwest Atlantic Ocean using genetic-based methods. Capelin collected at ten beach and five demersal spawning locations over the period 2002 through 2008 (N = 3,433 fish) were genotyped using six polymorphic microsatellite loci. Temporally distinct samples were identified at three beach spawning locations: Chance Cove, Little Lawn and Straitsview, Newfoundland. Four capelin stocks are assumed for fisheries management in the northwest Atlantic Ocean based on meristics, morphometrics, tag returns, and seasonal distribution patterns. Our results suggested groupings that were somewhat different than the assumed structure, and indicate at least seven genetically defined populations arising from two ancestral populations. The spatial mosaic of capelin from each of the two basal cluster groups explains much of the observed geographic variability amongst neighbouring samples. The genetic-defined populations were resolved at Jost’s D_est ≥ 0.01 and were composed of fish collected 1) in the Gulf of St. Lawrence, 2) along the south and east coasts of Newfoundland, 3) along coastal northern Newfoundland and southern Labrador, 4) along coastal northern Labrador, 5) near the Saguenay River, and at two nearshore demersal spawning sites, 6) one at Grebes Nest off Bellevue Beach on the east coast of Newfoundland, and 7) one off the coast of Labrador at Domino Run. Moreover, the offshore demersal spawners on the Scotian Shelf and Southeast Shoal appeared to be related to the inshore demersal spawners at Grebes Nest and in Domino Run and to beach spawners from the Gulf of St. Lawrence.

Population genomic evidence for adaptive differentiation in Baltic Sea three-spined sticklebacks

Background

The degree of genetic differentiation among populations experiencing high levels of gene flow is expected to be low for neutral genomic sites, but substantial divergence can occur in sites subject to directional selection. Studies of highly mobile marine fish populations provide an opportunity to investigate this kind of heterogeneous genomic differentiation, but most studies to this effect have focused on a relatively low number of genetic markers and/or few populations. Hence, the patterns and extent of genomic divergence in high-gene-flow marine fish populations remain poorly understood.

Results

We here investigated genome-wide patterns of genetic variability and differentiation in ten marine populations of three-spined stickleback (Gasterosteus aculeatus) distributed across a steep salinity and temperature gradient in the Baltic Sea, by utilizing >30,000 single nucleotide polymorphisms obtained with a pooled RAD-seq approach. We found that genetic diversity and differentiation varied widely across the genome, and identified numerous fairly narrow genomic regions exhibiting signatures of both divergent and balancing selection. Evidence was uncovered for substantial genetic differentiation associated with both salinity and temperature gradients, and many candidate genes associated with local adaptation in the Baltic Sea were identified.

Conclusions

The patterns of genetic diversity and differentiation, as well as candidate genes associated with adaptation, in Baltic Sea sticklebacks were similar to those observed in earlier comparisons between marine and freshwater populations, suggesting that similar processes may be driving adaptation to brackish and freshwater environments. Taken together, our results provide strong evidence for heterogenic genomic divergence driven by local adaptation in the face of gene flow along an environmental gradient in the post-glacially formed Baltic Sea.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0130-8) contains supplementary material, which is available to authorized users.