Detection of outlier loci and their utility for fisheries management

Genetic investigation of population structure in Atlantic chub mackerel, Scomber colias Gmelin, 1789 along the West African coast

Sustainable management of transboundary fish stocks hinges on accurate delineation of population structure. Genetic analysis offers a powerful tool to identify potential subpopulations within a seemingly homogenous stock, facilitating the development of effective, coordinated management strategies across international borders. Along the West African coast, the Atlantic chub mackerel (Scomber colias) is a commercially important and ecologically significant species, yet little is known about its genetic population structure and connectivity. Currently, the stock is managed as a single unit in West African waters despite new research suggesting morphological and adaptive differences. Here, eight microsatellite loci were genotyped on 1,169 individuals distributed across 33 sampling sites from Morocco (27.39°N) to Namibia (22.21°S). Bayesian clustering analysis depicts one homogeneous population across the studied area with null overall differentiation (F ST = 0.0001ns), which suggests panmixia and aligns with the migratory potential of this species. This finding has significant implications for the effective conservation and management of S. colias within a wide scope of its distribution across West African waters from the South of Morocco to the North-Centre of Namibia and underscores the need for increased regional cooperation in fisheries management and conservation.

Genetic analyses reveal panmixia in Indian waters and population subdivision across Indian Ocean and Indo-Malay Archipelago for Decapterus russelli

The Indian Scad, Decapterus russelli is an important pelagic carangid widely distributed throughout the Indian Ocean and the Indo-West Pacific. Despite being widely distributed in the Indian Ocean, the information regarding genetic structuring and diversity has been lacking compared to its Indo Malay Archipelago counterparts. The present study was conducted to investigate the genetic stock structure of D. russelli based on mitochondrial (Cyt b) and nuclear (DrAldoB1) markers along Indian waters. The results indicated the presence of a single panmictic stock across the Indian Ocean region. High haplotype diversity associated with low nucleotide diversity suggested a population bottleneck followed by rapid population growth. Phylogenetic analysis revealed the absence of geographical clustering of lineages with the most common haplotype distributed globally. The pelagic life style, migratory capabilities, and larval dispersal may be the contributing factors to the observed spatial homogeneity of D. russelli. However, significant genetic differentiation was observed between the populations from Indian Ocean and Indo-Malay Archipelago. Hierarchical molecular variance analysis (AMOVA), pairwise FST comparisons and SAMOVA showed existence of two distinct genetic stocks of D. russelli in the Indian Ocean and IMA. The observed interpopulation genetic variation was high. A plausible explanation for the genetic differentiation observed between the Indo-Malay Archipelago and the Indian Ocean regions suggest the influence of historic isolation, ocean surface currents and biotic and abiotic features of the ocean. Also, there was a significant relationship between genetic distance and geographical distance between population pairs in a manner consistent with isolation-by-distance. These resulted in the evolution of a phylogeographic break for this species between these regions. The findings of these results suggest that D. russelli from the Indian Ocean shall be managed in its entire area of distribution as a single stock. Further, the Indian Ocean and Indo-Malayan stocks can be managed separately.

Aquaculture-driven evolution of the salmon louse mtDNA genome

Resistance toward the antiparasitic pyrethroid, deltamethrin, is reported in the Atlantic salmon louse (Lepeophtheirus salmonis salmonis), a persistent ectoparasite of farmed and wild salmonids. The resistance mechanism is linked to mitochondrial DNA (mtDNA), where genetic markers for resistance have been identified. Here, we investigated how widespread pyrethroid use in aquaculture may have influenced mtDNA variation in lice, and the dispersion of resistant haplotypes across the North Atlantic, using historical (2000-2002 "pre-resistance") and contemporary (2014-2017 "post-resistance") samples. To study this, we sequenced ATPase 6 and cytochrome b, genotyped two genetic markers for deltamethrin resistance, and genotyped microsatellites as "neutral" controls of potential population bottlenecks. Overall, we observed a modest reduction in mtDNA diversity in the period 2000-2017, but no reduction in microsatellite variation was observed. The reduction in mtDNA variation was especially distinct in two of the contemporary samples, fixed for one and two haplotypes, respectively. By contrast, all historical samples consisted of close to one mtDNA haplotype per individual. No population genetic structure was detected among the historical samples for mtDNA nor microsatellites. By contrast, significant population genetic differentiation was observed for mtDNA among some of the contemporary samples. However, the observed population genetic structure was tightly linked with the pattern of deltamethrin resistance, and we therefore conclude that it primarily reflects the transient mosaic of pyrethroid usage in time and space. Two historically undetected mtDNA haplotypes dominated in the contemporary samples, both of which were linked to deltamethrin resistance, demonstrating primarily two origins of deltamethrin resistance in the North Atlantic. Collectively, these data demonstrate that the widespread use of pyrethroids in commercial aquaculture has substantially altered the patterns of mtDNA diversity in lice across the North Atlantic, and that long-distance dispersion of resistance is rapid due to high level of genetic connectivity that is observed in this species.

The paradox of adaptive trait clines with nonclinal patterns in the underlying genes

Multivariate climate change presents an urgent need to understand how species adapt to complex environments. Population genetic theory predicts that loci under selection will form monotonic allele frequency clines with their selective environment, which has led to the wide use of genotype-environment associations (GEAs). This study used a set of simulations to elucidate the conditions under which allele frequency clines are more or less likely to evolve as multiple quantitative traits adapt to multivariate environments. Phenotypic clines evolved with nonmonotonic (i.e., nonclinal) patterns in allele frequencies under conditions that promoted unique combinations of mutations to achieve the multivariate optimum in different parts of the landscape. Such conditions resulted from interactions among landscape, demography, pleiotropy, and genetic architecture. GEA methods failed to accurately infer the genetic basis of adaptation under a range of scenarios due to first principles (clinal patterns did not evolve) or statistical issues (clinal patterns evolved but were not detected due to overcorrection for structure). Despite the limitations of GEAs, this study shows that a back-transformation of multivariate ordination can accurately predict individual multivariate traits from genotype and environmental data regardless of whether inference from GEAs was accurate. In addition, frameworks are introduced that can be used by empiricists to quantify the importance of clinal alleles in adaptation. This research highlights that multivariate trait prediction from genotype and environmental data can lead to accurate inference regardless of whether the underlying loci display clinal or nonmonotonic patterns.

Genotyping-in-Thousands by sequencing panel development and application to inform kokanee salmon (Oncorhynchus nerka) fisheries management at multiple scales

The ability to differentiate life history variants is vital for estimating fisheries management parameters, yet traditional survey methods can be inaccurate in mixed-stock fisheries. Such is the case for kokanee, the freshwater resident form of sockeye salmon (Oncorhynchus nerka), which exhibits various reproductive ecotypes (stream-, shore-, deep-spawning) that co-occur with each other and/or anadromous O. nerka in some systems across their pan-Pacific distribution. Here, we developed a multi-purpose Genotyping-in-Thousands by sequencing (GT-seq) panel of 288 targeted single nucleotide polymorphisms (SNPs) to enable accurate kokanee stock identification by geographic basin, migratory form, and reproductive ecotype across British Columbia, Canada. The GT-seq panel exhibited high self-assignment accuracy (93.3%) and perfect assignment of individuals not included in the baseline to their geographic basin, migratory form, and reproductive ecotype of origin. The GT-seq panel was subsequently applied to Wood Lake, a valuable mixed-stock fishery, revealing high concordance (>98%) with previous assignments to ecotype using microsatellites and TaqMan® SNP genotyping assays, while improving resolution, extending a long-term time-series, and demonstrating the scalability of this approach for this system and others.

Genome-Wide Investigation of the Multiple Origins Hypothesis for Deep-Spawning Kokanee Salmon (Oncorhynchus nerka) across its Pan-Pacific Distribution

Salmonids have emerged as important study systems for investigating molecular processes underlying parallel evolution given their tremendous life history variation. Kokanee, the resident form of anadromous sockeye salmon (Oncorhynchus nerka), have evolved multiple times across the species' pan-Pacific distribution, exhibiting multiple reproductive ecotypes including those that spawn in streams, on lake-shores, and at lake depths >50 m. The latter has only been detected in 5 locations in Japan and British Columbia, Canada. Here, we investigated the multiple origins hypothesis for deep-spawning kokanee, using 9721 single nucleotide polymorphisms distributed across the genome analyzed for the vast majority of known populations in Japan (Saiko Lake) and Canada (Anderson, Seton, East Barrière Lakes) relative to stream-spawning populations in both regions. We detected 397 outlier loci, none of which were robustly identified in paired-ecotype comparisons in Japan and Canada independently. Bayesian clustering and principal components analyses based on neutral loci revealed 6 distinct clusters, largely associated with geography or translocation history, rather than ecotype. Moreover, a high level of divergence between Canadian and Japanese populations, and between deep- and stream-spawning populations regionally, suggests the deep-spawning ecotype independently evolved on the 2 continents. On a finer level, Japanese kokanee populations exhibited low estimates of heterozygosity, significant levels of inbreeding, and reduced effective population sizes relative to Canadian populations, likely associated with transplantation history. Along with preliminary evidence for hybridization between deep- and stream-spawning ecotypes in Saiko Lake, these findings should be considered within the context of on-going kokanee fisheries management in Japan.

Genome-wide analysis reveals demographic and life-history patterns associated with habitat modification in landlocked, deep-spawning sockeye salmon (Oncorhynchus nerka)

Human-mediated habitat fragmentation in freshwater ecosystems can negatively impact genetic diversity, demography, and life history of native biota, while disrupting the behavior of species that are dependent on spatial connectivity to complete their life cycles. In the Alouette River system (British Columbia, Canada), dam construction in 1928 impacted passage of anadromous sockeye salmon (Oncorhynchus nerka), with the last records of migrants occurring in the 1930s. Since that time, O. nerka persisted as a resident population in Alouette Reservoir until experimental water releases beginning in 2005 created conditions for migration; two years later, returning migrants were observed for the first time in ~70 years, raising important basic and applied questions regarding life-history variation and population structure in this system. Here, we investigated the genetic distinctiveness and population history of Alouette Reservoir O. nerka using genome-wide SNP data (n = 7,709 loci) collected for resident and migrant individuals, as well as for neighboring anadromous sockeye salmon and resident kokanee populations within the Fraser River drainage (n = 312 individuals). Bayesian clustering and principal components analyses based on neutral loci revealed five distinct clusters, largely associated with geography, and clearly demonstrated that Alouette Reservoir resident and migrant individuals are genetically distinct from other O. nerka populations in the Fraser River drainage. At a finer level, there was no clear evidence for differentiation between Alouette Reservoir residents and migrants; although we detected eight high-confidence outlier loci, they all mapped to sex chromosomes suggesting that differences were likely due to uneven sex ratios rather than life history. Taken together, these data suggest that contemporary Alouette Reservoir O. nerka represents a landlocked sockeye salmon population, constituting the first reported instance of deep-water spawning behavior associated with this life-history form. This finding punctuates the need for reassessment of conservation status and supports ongoing fisheries management activities in Alouette Reservoir.

Weak spatial-genetic structure in a native invasive, the southern pine beetle (Dendroctonus frontalis), across the eastern United States

The southern pine beetle, Dendroctonus frontalis, is a native pest of pine trees that has recently expanded its range into the northeastern United States. Understanding its colonization, dispersal, and connectivity will be critical for mitigating negative economic and ecological impacts in the newly invaded areas. Characterization of spatial-genetic structure can contribute to this; however, previous studies have reached different conclusions about regional population genetic structure, with one study reporting a weak east-west pattern, and the most recent reporting an absence of structure. Here we systematically assessed several explanations for the absence of spatial-genetic structure. To do this, we developed nine new microsatellite markers and combined them with an existing 24-locus data matrix for the same individuals. We then reanalyzed this full dataset alongside datasets in which certain loci were omitted with the goal of creating more favorable signal to noise ratios. We also partitioned the data based on the sex of D. frontalis individuals, and then employed a broad suite of genotypic clustering and isolation-by-distance (IBD) analyses. We found that neither inadequate information content in the molecular marker set, nor unfavorable signal-to-noise ratio, nor insensitivity of the analytical approaches could explain the absence of structure. Regardless of dataset composition, there was little evidence for clusters (i.e., distinct geo-genetic groups) or clines (i.e., gradients of increasing allele frequency differences over larger geographic distances), with one exception: significant IBD was repeatedly detected using an individual-based measure of relatedness whenever datasets included males (but not for female-only datasets). This is strongly indicative of broad-scale female-biased dispersal, which has not previously been reported for D. frontalis, in part owing to logistical limitations of direct approaches (e.g., capture-mark-recapture). Weak spatial-genetic structure suggests long-distance connectivity and that gene flow is high, but additional research is needed to understand range expansion and outbreak dynamics in this species using alternate approaches.

The use of non-lethal sampling for transcriptomics to assess the physiological status of wild fishes

Fishes respond to different abiotic and biotic stressors through changes in gene expression as a part of an integrated physiological response. Transcriptomics approaches have been used to quantify gene expression patterns as a reductionist approach to understand responses to environmental stressors in animal physiology and have become more commonly used to study wild fishes. We argue that non-lethal sampling for transcriptomics should become the norm for assessing the physiological status of wild fishes, especially when there are conservation implications. Processes at the level of the transcriptome provide a "snapshot" of the cellular conditions at a given time; however, by using a non-lethal sampling protocol, researchers can connect the transcriptome profile with fitness-relevant ecological endpoints such as reproduction, movement patterns and survival. Furthermore, telemetry is a widely used approach in fisheries to understand movement patterns in the wild, and when combined with transcriptional profiling, provides arguably the most powerful use of non-lethal sampling for transcriptomics in wild fishes. In this review, we discuss the different tissues that can be successfully incorporated into non-lethal sampling strategies, which is particularly useful in the context of the emerging field of conservation transcriptomics. We briefly describe different methods for transcriptional profiling in fishes from high-throughput qPCR to whole transcriptome approaches. Further, we discuss strategies and the limitations of using transcriptomics for non-lethally studying fishes. Lastly, as 'omics' technology continues to advance, transcriptomics paired with different omics approaches to study wild fishes will provide insight into the factors that regulate phenotypic variation and the physiological responses to changing environmental conditions in the future.

Population Genetics of Odontarrhena (Brassicaceae) from Albania: The Effects of Anthropic Habitat Disturbance, Soil, and Altitude on a Ni-Hyperaccumulator Plant Group from a Major Serpentine Hotspot

Albanian taxa and populations of the genus Odontarrhena are most promising candidates for research on metal tolerance and Ni-agromining, but their genetic structure remains unknown. We investigated phylogenetic relationships and genetic differentiation in relation to distribution and ploidy of the taxa, anthropic site disturbance, elevation, soil type, and trace metals at each population site. After performing DNA sequencing of selected accessions, we applied DNA-fingerprinting to analyze the genetic structure of 32 populations from ultramafic and non-ultramafic outcrops across Albania. Low sequence divergence resulted in poorly resolved phylograms, but supported affinity between the two diploid serpentine endemics O. moravensis and O. rigida. Analysis of molecular variance (AMOVA) revealed significant population differentiation, but no isolation by distance. Among-population variation was higher in polyploids than in diploids, in which genetic distances were lower. Genetic admixing at population and individual level occurred especially in the polyploids O. chalcidica, O. decipiens, and O. smolikana. Admixing increased with site disturbance. Outlier loci were higher in serpentine populations but decreased along altitude with lower drought and heat stress. Genetic variability gained by gene flow and hybridization at contact zones with "resident" species of primary ultramafic habitats promoted expansion of the tetraploid O. chalcidica across anthropogenic sites.

Genomic analysis reveals neutral and adaptive patterns that challenge the current management regime for East Atlantic cod Gadus morhua L

Challenging long-held perceptions of fish management units can help to protect vulnerable stocks. When a fishery consisting of multiple genetic stocks is managed as a single unit, overexploitation and depletion of minor genetic units can occur. Atlantic cod (Gadus morhua) is an economically and ecologically important marine species across the North Atlantic. The application of new genomic resources, including SNP arrays, allows us to detect and explore novel structure within specific cod management units. In Norwegian waters, coastal cod (i.e. those not undertaking extensive migrations) are divided into two arbitrary management units defined by ICES: one between 62° and 70°N (Norwegian coastal cod; NCC) and one between 58° and 62°N (Norwegian coastal south; NCS). Together, these capture a fishery area of >25,000 km2 containing many spawning grounds. To assess whether these geographic units correctly represent genetic stocks, we analysed spawning cod of NCC and NCS for more than 8,000 SNPs along with samples of Russian White Sea cod, north-east Arctic cod (NEAC: the largest Atlantic stock), and outgroup samples representing the Irish and Faroe Sea's. Our analyses revealed large differences in spatial patterns of genetic differentiation across the genome and revealed a complex biological structure within NCC and NCS. Haplotype maps from four chromosome sets show regional specific SNP indicating a complex genetic structure. The current management plan dividing the coastal cod into only two management units does not accurately reflect the genetic units and needs to be revised. Coastal cod in Norway, while highly heterogenous, is also genetically distinct from neighbouring stocks in the north (NEAC), west (Faroe Island) and the south. The White Sea cod are highly divergent from other cod, possibly yielding support to the earlier notion of subspecies rank.

Genome wide analysis reveals genetic divergence between Goldsinny wrasse populations

Background

Marine fish populations are often characterized by high levels of gene flow and correspondingly low genetic divergence. This presents a challenge to define management units. Goldsinny wrasse (Ctenolabrus rupestris) is a heavily exploited species due to its importance as a cleaner-fish in commercial salmonid aquaculture. However, at the present, the population genetic structure of this species is still largely unresolved. Here, full-genome sequencing was used to produce the first genomic reference for this species, to study population-genomic divergence among four geographically distinct populations, and, to identify informative SNP markers for future studies.

Results

After construction of a de novo assembly, the genome was estimated to be highly polymorphic and of ~600Mbp in size. 33,235 SNPs were thereafter selected to assess genomic diversity and differentiation among four populations collected from Scandinavia, Scotland, and Spain. Global F_ST among these populations was 0.015–0.092. Approximately 4% of the investigated loci were identified as putative global outliers, and ~ 1% within Scandinavia. SNPs showing large divergence (F_ST > 0.15) were picked as candidate diagnostic markers for population assignment. One hundred seventy-three of the most diagnostic SNPs between the two Scandinavian populations were validated by genotyping 47 individuals from each end of the species’ Scandinavian distribution range. Sixty-nine of these SNPs were significantly (p < 0.05) differentiated (mean F_{ST_173_loci} = 0.065, F_{ST_69_loci} = 0.140). Using these validated SNPs, individuals were assigned with high probability (≥ 94%) to their populations of origin.

Conclusions

Goldsinny wrasse displays a highly polymorphic genome, and substantial population genomic structure. Diversifying selection likely affects population structuring globally and within Scandinavia. The diagnostic loci identified now provide a promising and cost-efficient tool to investigate goldsinny wrasse populations further.

Genetic analysis redraws the management boundaries for the European sprat

Sustainable fisheries management requires detailed knowledge of population genetic structure. The European sprat is an important commercial fish distributed from Morocco to the Arctic circle, Baltic, Mediterranean, and Black seas. Prior to 2018, annual catch advice on sprat from the International Council for the Exploration of the Sea (ICES) was based on five putative stocks: (a) North Sea, (b) Kattegat-Skagerrak and Norwegian fjords, (c) Baltic Sea, (d) West of Scotland-southern Celtic Seas, and (e) English Channel. However, there were concerns that the sprat advice on stock size estimates management plan inadequately reflected the underlying biological units. Here, we used ddRAD sequencing to develop 91 SNPs that were thereafter used to genotype approximately 2,500 fish from 40 locations. Three highly distinct and relatively homogenous genetic groups were identified: (a) Norwegian fjords; (b) Northeast Atlantic including the North Sea, Kattegat-Skagerrak, Celtic Sea, and Bay of Biscay; and (c) Baltic Sea. Evidence of genetic admixture and possibly physical mixing was detected in samples collected from the transition zone between the North and Baltic seas, but not between any of the other groups. These results have already been implemented by ICES with the decision to merge the North Sea and the Kattegat-Skagerrak sprat to be assessed as a single unit, thus demonstrating that genetic data can be rapidly absorbed to align harvest regimes and biological units.

Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye (Sander vitreus) in a large freshwater ecosystem

RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.

Transforming Ocean Conservation: Applying the Genetic Rescue Toolkit

Although oceans provide critical ecosystem services and support the most abundant populations on earth, the extent of damage impacting oceans and the diversity of strategies to protect them is disconcertingly, and disproportionately, understudied. While conventional modes of conservation have made strides in mitigating impacts of human activities on ocean ecosystems, those strategies alone cannot completely stem the tide of mounting threats. Biotechnology and genomic research should be harnessed and developed within conservation frameworks to foster the persistence of viable ocean ecosystems. This document distills the results of a targeted survey, the Ocean Genomics Horizon Scan, which assessed opportunities to bring novel genetic rescue tools to marine conservation. From this Horizon Scan, we have identified how novel approaches from synthetic biology and genomics can alleviate major marine threats. While ethical frameworks for biotechnological interventions are necessary for effective and responsible practice, here we primarily assessed technological and social factors directly affecting technical development and deployment of biotechnology interventions for marine conservation. Genetic insight can greatly enhance established conservation methods, but the severity of many threats may demand genomic intervention. While intervention is controversial, for many marine areas the cost of inaction is too high to allow controversy to be a barrier to conserving viable ecosystems. Here, we offer a set of recommendations for engagement and program development to deploy genetic rescue safely and responsibly.

Population Genomics of an Anadromous Hilsa Shad Tenualosa ilisha Species across Its Diverse Migratory Habitats: Discrimination by Fine-Scale Local Adaptation

The migration of anadromous fish in heterogenic environments unceasingly imposes a selective pressure that results in genetic variation for local adaptation. However, discrimination of anadromous fish populations by fine-scale local adaptation is challenging because of their high rate of gene flow, highly connected divergent population, and large population size. Recent advances in next-generation sequencing (NGS) have expanded the prospects of defining the weakly structured population of anadromous fish. Therefore, we used NGS-based restriction site-associated DNA (NextRAD) techniques on 300 individuals of an anadromous Hilsa shad (Tenualosa ilisha) species, collected from nine strategic habitats, across their diverse migratory habitats, which include sea, estuary, and different freshwater rivers. The NextRAD technique successfully identified 15,453 single nucleotide polymorphism (SNP) loci. Outlier tests using the F_ST OutFLANK and pcadapt approaches identified 74 and 449 SNPs (49 SNPs being common), respectively, as putative adaptive loci under a divergent selection process. Our results, based on the different cluster analyses of these putatively adaptive loci, suggested that local adaptation has divided the Hilsa shad population into two genetically structured clusters, in which marine and estuarine collection sites were dominated by individuals of one genetic cluster and different riverine collection sites were dominated by individuals of another genetic cluster. The phylogenetic analysis revealed that all the riverine populations of Hilsa shad were further subdivided into the north-western riverine (turbid freshwater) and the north-eastern riverine (clear freshwater) ecotypes. Among all of the putatively adaptive loci, only 36 loci were observed to be in the coding region, and the encoded genes might be associated with important biological functions related to the local adaptation of Hilsa shad. In summary, our study provides both neutral and adaptive contexts for the observed genetic divergence of Hilsa shad and, consequently, resolves the previous inconclusive findings on their population genetic structure across their diverse migratory habitats. Moreover, the study has clearly demonstrated that NextRAD sequencing is an innovative approach to explore how dispersal and local adaptation can shape genetic divergence of non-model anadromous fish that intersect diverse migratory habitats during their life-history stages.

Population genetics and historical demographic inferences of the blue crab Callinectes sapidus in the US based on microsatellites

The native range of the blue crab Callinectes sapidus spans Nova Scotia to northern Argentina. In the US, it constitutes a keystone species in estuarine habitats of the Atlantic coast and Gulf of Mexico (GOM), serving as both predator and prey to other species, and also has historically represented a multi-billion dollar fishery. Knowledge relevant to effective management and monitoring of this ecologically and economically important species, such as levels of population genetic differentiation and genetic diversity, is necessary. Although several population genetics studies have attempted to address these questions in one or more parts of its distribution, conflicting results and potential problems with the markers used, as well as other issues, have obscured our understanding on them. In this study, we examined large-scale genetic connectivity of the blue crab in the US, using 16 microsatellites, and genotyped individuals from Chesapeake Bay, in the US Atlantic, and from nine localities along the US GOM coast. Consistent with the high long-distance dispersal potential of this species, very low levels of genetic differentiation were detected for the blue crab among the ten US localities examined, suggesting it constitutes a large panmictic population within this region. Estimations of genetic diversity for the blue crab appear to be high in the US, and provide a baseline for monitoring temporal changes in this species. Demographic analyses indicate a recent range expansion of the US population, probably during the Holocene. In addition, capitalizing on published microsatellite data from southern Brazil, our analyses detected high genetic differentiation between localities in the US and Brazil. These results point to the need for examination of genetic diversity and differentiation along the area spanning the US to southern Brazil.

Genotyping-in-Thousands by sequencing (GT-seq) panel development and application to minimally invasive DNA samples to support studies in molecular ecology

Minimally invasive sampling (MIS) is widespread in wildlife studies; however, its utility for massively parallel DNA sequencing (MPS) is limited. Poor sample quality and contamination by exogenous DNA can make MIS challenging to use with modern genotyping-by-sequencing approaches, which have been traditionally developed for high-quality DNA sources. Given that MIS is often more appropriate in many contexts, there is a need to make such samples practical for harnessing MPS. Here, we test the ability for Genotyping-in-Thousands by sequencing (GT-seq), a multiplex amplicon sequencing approach, to effectively genotype minimally invasive cloacal DNA samples collected from the Western Rattlesnake (Crotalus oreganus), a threatened species in British Columbia, Canada. As there was no previous genetic information for this species, an optimized panel of 362 SNPs was selected for use with GT-seq from a de novo restriction site-associated DNA sequencing (RADseq) assembly. Comparisons of genotypes generated within and among RADseq and GT-seq for the same individuals found low rates of genotyping error (GT-seq: 0.50%; RADseq: 0.80%) and discordance (2.57%), the latter likely due to the different genotype calling models employed. GT-seq mean genotype discordance between blood and cloacal swab samples collected from the same individuals was also minimal (1.37%). Estimates of population diversity parameters were similar across GT-seq and RADseq data sets, as were inferred patterns of population structure. Overall, GT-seq can be effectively applied to low-quality DNA samples, minimizing the inefficiencies presented by exogenous DNA typically found in minimally invasive samples and continuing the expansion of molecular ecology and conservation genetics in the genomics era.

Can we rely on selected genetic markers for population identification? Evidence from coastal Atlantic cod

The use of genetic markers under putative selection in population studies carries the potential for erroneous identification of populations and misassignment of individuals to population of origin. Selected markers are nevertheless attractive, especially in marine organisms that are characterized by weak population structure at neutral loci. Highly fecund species may tolerate the cost of strong selective mortality during early life stages, potentially leading to a shift in offspring genotypes away from the parental proportions. In Atlantic cod, recent genetic studies have uncovered different genotype clusters apparently representing phenotypically cryptic populations that coexist in coastal waters. Here, we tested if a high-graded SNP panel specifically designed to classify individual cod to population of origin may be unreliable because of natural selection acting on the SNPs or their linked background. Temporal samples of cod were collected from two fjords, starting at the earliest life stage (pelagic eggs) and carried on until late autumn (bottom-settled juveniles), covering the period during summer of high natural mortality. Despite the potential for selective mortality during the study period, we found no evidence for selection, as both cod types occurred throughout the season, already in the earliest egg samples, and there was no evidence for a shift during the season in the proportions of one or the other type. We conclude that high-graded marker panels under putative natural selection represent a valid and useful tool for identifying biological population structure in this highly fecund species and presumably in others.

The population genomics of yellowfin tuna (Thunnus albacares) at global geographic scale challenges current stock delineation

Yellowfin tuna, Thunnus albacares, is one of the most important seafood commodities in the world. Despite its great biological and economic importance, conflicting evidence arises from classical genetic and tagging studies concerning the yellowfin tuna population structure at local and global oceanic scales. Access to more powerful and cost effective genetic tools would represent the first step towards resolving the population structure of yellowfin tuna across its distribution range. Using a panel of 939 neutral Single Nucleotide Polymorphisms (SNPs), and the most comprehensive data set of yellowfin samples available so far, we found genetic differentiation among the Atlantic, Indian and Pacific oceans. The genetic stock structure analysis carried out with 33 outlier SNPs, putatively under selection, identified discrete populations within the Pacific Ocean and, for the first time, also within the Atlantic Ocean. Stock assessment approaches that consider genetic differences at neutral and adaptive genomic loci should be routinely implemented to check the status of the yellowfin tuna, prevent illegal trade, and develop more sustainable management measures.

Genomic Changes Associated with Reproductive and Migratory Ecotypes in Sockeye Salmon (Oncorhynchus nerka)

Mechanisms underlying adaptive evolution can best be explored using paired populations displaying similar phenotypic divergence, illuminating the genomic changes associated with specific life history traits. Here, we used paired migratory [anadromous vs. resident (kokanee)] and reproductive [shore- vs. stream-spawning] ecotypes of sockeye salmon (Oncorhynchus nerka) sampled from seven lakes and two rivers spanning three catchments (Columbia, Fraser, and Skeena) in British Columbia, Canada to investigate the patterns and processes underlying their divergence. Restriction-site associated DNA sequencing was used to genotype this sampling at 7,347 single nucleotide polymorphisms, 334 of which were identified as outlier loci and candidates for divergent selection within at least one ecotype comparison. Sixty-eight of these outliers were present in two or more comparisons, with 33 detected across multiple catchments. Of particular note, one locus was detected as the most significant outlier between shore and stream-spawning ecotypes in multiple comparisons and across catchments (Columbia, Fraser, and Snake). We also detected several genomic islands of divergence, some shared among comparisons, potentially showing linked signals of differential selection. The single nucleotide polymorphisms and genomic regions identified in our study offer a range of mechanistic hypotheses associated with the genetic basis of O. nerka life history variation and provide novel tools for informing fisheries management.

The molecular basis of differential morphology and bleaching thresholds in two morphs of the coral Pocillopora acuta

Processes of cnidarian evolution, including hybridization and phenotypic plasticity, have complicated the clear diagnosis of species boundaries within the phylum. Pocillopora acuta, a species of scleractinian coral that was recently split from the widespread Pocillopora damicornis species complex, occurs in at least two distinct morphs on the Great Barrier Reef. Contrasting morphology combined with evidence of differential bleaching thresholds among sympatrically distributed colonies suggest that the taxonomy of this recently described species is not fully resolved and may represent its own species complex. To examine the basis of sympatric differentiation between the two morphs, we combined analyses of micro- and macro-skeletal morphology with genome wide sequencing of the coral host, as well as ITS2 genotyping of the associated Symbiodinium communities. We found consistent differences between morphs on both the macro- and micro-skeletal scale. In addition, we identified 18 candidate functional genes that relate to skeletal formation and morphology that may explain how the two morphs regulate growth to achieve their distinct growth forms. With inconclusive results in endosymbiotic algal community diversity between the two morphs, we propose that colony morphology may be linked to bleaching susceptibility. We conclude that cryptic speciation may be in the early stages within the species P. acuta.

Genomic signatures of local adaptation reveal source-sink dynamics in a high gene flow fish species

Understanding source-sink dynamics is important for conservation management, particularly when climatic events alter species’ distributions. Following a 2011 ‘marine heatwave’ in Western Australia, we observed high recruitment of the endemic fisheries target species Choerodon rubescens, towards the cooler (southern) end of its distribution. Here, we use a genome wide set of 14 559 single-nucleotide polymorphisms (SNPs) to identify the likely source population for this recruitment event. Most loci (76%) showed low genetic divergence across the species’ range, indicating high levels of gene flow and confirming previous findings using neutral microsatellite markers. However, a small proportion of loci showed strong patterns of differentiation and exhibited patterns of population structure consistent with local adaptation. Clustering analyses based on these outlier loci indicated that recruits at the southern end of C. rubescens’ range originated 400 km to the north, at the centre of the species’ range, where average temperatures are up to 3 °C warmer. Survival of these recruits may be low because they carry alleles adapted to an environment different to the one they now reside in, but their survival is key to establishing locally adapted populations at and beyond the range edge as water temperatures increase with climate change.

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.

From population genomics to conservation and management: a workflow for targeted analysis of markers identified using genome-wide approaches in Atlantic salmon Salmo salar

A genotyping assay for the Ion Torrent Ion PGM platform was developed for fast and cost-effective targeted genotyping of key single nucleotide polymorphisms (SNPs) earlier identified using a genome-wide SNP array in Atlantic salmon Salmo salar. The method comprised a simple primer design step for multiplex-polymerase chain reaction (PCR), followed by two rounds of Ion Torrent Ion PGM sequencing to empirically evaluate marker efficiency in large multiplexes and to optimise or exclude them when necessary. Of 282 primer pairs initially tested, 217 were successfully amplified, indicating good amplification success (>75%). These markers included the sdy partial gene product to determine genetic sex, as well as three additional modules comprising SNPs for assessing neutral genetic variation (N_SNP = 150), examining functional genetic variation associated with sea age at maturity (N_SNP = 5), and for performing genetic subpopulation assignment (N_SNP = 61). The assay was primarily developed to monitor long-term genetic changes in S. salar from the Teno River, but modules are likely suitable for application in a wide range of S. salar populations. Furthermore, the fast and versatile assay development pipeline offers a strategy for developing targeted sequencing assays in any species.

Single nucleotide polymorphisms to discriminate different classes of hybrid between wild Atlantic salmon and aquaculture escapees

Many wild Atlantic salmon (Salmo salar) populations are threatened by introgressive hybridization from domesticated fish that have escaped from aquaculture facilities. A detailed understanding of the hybridization dynamics between wild salmon and aquaculture escapees requires discrimination of different hybrid classes; however, markers currently available to discriminate the two types of parental genome have limited power to do this. Using a high‐density Atlantic salmon single nucleotide polymorphism (SNP) array, in combination with pooled‐sample allelotyping and an F_st outlier approach, we identified 200 SNPs that differentiated an important Atlantic salmon stock from the escapees potentially hybridizing with it. By simulating multiple generations of wild–escapee hybridization, involving wild populations in two major phylogeographic lineages and a genetically diverse set of escapees, we showed that both the complete set of SNPs and smaller subsets could reliably assign individuals to different hybrid classes up to the third hybrid (F3) generation. This set of markers will be a useful tool for investigating the genetic interactions between native wild fish and aquaculture escapees in many Atlantic salmon populations.

Adaptive and neutral markers both show continent-wide population structure of mountain pine beetle (Dendroctonus ponderosae)

Assessments of population genetic structure and demographic history have traditionally been based on neutral markers while explicitly excluding adaptive markers. In this study, we compared the utility of putatively adaptive and neutral single‐nucleotide polymorphisms (SNPs) for inferring mountain pine beetle population structure across its geographic range. Both adaptive and neutral SNPs, and their combination, allowed range‐wide structure to be distinguished and delimited a population that has recently undergone range expansion across northern British Columbia and Alberta. Using an equal number of both adaptive and neutral SNPs revealed that adaptive SNPs resulted in a stronger correlation between sampled populations and inferred clustering. Our results suggest that adaptive SNPs should not be excluded prior to analysis from neutral SNPs as a combination of both marker sets resulted in better resolution of genetic differentiation between populations than either marker set alone. These results demonstrate the utility of adaptive loci for resolving population genetic structure in a nonmodel organism.

Genomic evidence for ecological divergence against a background of population homogeneity in the marine snail Chlorostoma funebralis

The balance between natural selection, gene flow and genetic drift is difficult to resolve in marine invertebrates with extensive dispersal and fluctuating population sizes. The intertidal snail Chlorostoma funebralis has planktonic larvae and previous work using mtDNA polymorphism reported no genetic population structure. Nevertheless, recent studies have documented differences in thermal tolerance and transcriptomic responses to heat stress between northern and southern California, USA, populations. To gain insight into the dynamics influencing adaptive divergence, we used double-digest restriction site-associated DNA (ddRAD) sequencing to identify 1861 genomewide, quality-filtered single-nucleotide polymorphism (SNP) loci for C. funebralis collected from three northern and three southern California sites (15 individuals per population). Considering all SNPs, there was no evidence for genetic differentiation among populations or regions (average FST  = 0.0042). However, outlier tests revealed 34 loci putatively under divergent selection between northern and southern populations, and structure and SNP tree analyses based on these outliers show clear genetic differentiation between geographic regions. Three of these outliers are known or hypothesized to be involved in stress granule formation, a response to environmental stress such as heat. Combined with previous work that found thermally tolerant southern populations show high baseline expression of stress response genes, these results further suggest that thermal stress is a strong selective pressure across C. funebralis populations. Overall, this study increases our understanding of the factors constraining local adaptation in marine organisms, while suggesting that ecologically driven, strong differentiation can occur at relevant loci in a species with planktonic larvae.

Genome-Wide SNP Discovery, Genotyping and Their Preliminary Applications for Population Genetic Inference in Spotted Sea Bass (Lateolabrax maculatus)

Next-generation sequencing and the collection of genome-wide single-nucleotide polymorphisms (SNPs) allow identifying fine-scale population genetic structure and genomic regions under selection. The spotted sea bass (Lateolabrax maculatus) is a non-model species of ecological and commercial importance and widely distributed in northwestern Pacific. A total of 22 648 SNPs was discovered across the genome of L. maculatus by paired-end sequencing of restriction-site associated DNA (RAD-PE) for 30 individuals from two populations. The nucleotide diversity (π) for each population was 0.0028±0.0001 in Dandong and 0.0018±0.0001 in Beihai, respectively. Shallow but significant genetic differentiation was detected between the two populations analyzed by using both the whole data set (F_ST = 0.0550, P < 0.001) and the putatively neutral SNPs (F_ST = 0.0347, P < 0.001). However, the two populations were highly differentiated based on the putatively adaptive SNPs (F_ST = 0.6929, P < 0.001). Moreover, a total of 356 SNPs representing 298 unique loci were detected as outliers putatively under divergent selection by F_ST-based outlier tests as implemented in BAYESCAN and LOSITAN. Functional annotation of the contigs containing putatively adaptive SNPs yielded hits for 22 of 55 (40%) significant BLASTX matches. Candidate genes for local selection constituted a wide array of functions, including binding, catalytic and metabolic activities, etc. The analyses with the SNPs developed in the present study highlighted the importance of genome-wide genetic variation for inference of population structure and local adaptation in L. maculatus.

Population Structure and Adaptive Divergence in a High Gene Flow Marine Fish: The Small Yellow Croaker (Larimichthys polyactis)

The spatial distribution of genetic diversity has been long considered as a key component of policy development for management and conservation of marine fishes. However, unraveling the population genetic structure of migratory fish species is challenging due to high potential for gene flow. Despite the shallow population differentiation revealed by putatively neutral loci, the higher genetic differentiation with panels of putatively adaptive loci could provide greater resolution for stock identification. Here, patterns of population differentiation of small yellow croaker (Larimichthys polyactis) were investigated by genotyping 15 highly polymorphic microsatellites in 337 individuals of 15 geographic populations collected from both spawning and overwintering grounds. Outlier analyses indicated that the locus Lpol03 might be under directional selection, which showed a strong homology with Grid2 gene encoding the glutamate receptor δ2 protein (GluRδ2). Based on Lpol03, two distinct clusters were identified by both STRUCTURE and PCoA analyses, suggesting that there were two overwintering aggregations of L. polyactis. A novel migration pattern was suggested for L. polyactis, which was inconsistent with results of previous studies based on historical fishing yield statistics. These results provided new perspectives on the population genetic structure and migratory routes of L. polyactis, which could have significant implications for sustainable management and utilization of this important fishery resource.

Biocomplexity in Populations of European Anchovy in the Adriatic Sea

The sustained exploitation of marine populations requires an understanding of a species' adaptive seascape so that populations can track environmental changes from short- and long-term climate cycles and from human development. The analysis of the distributions of genetic markers among populations, together with correlates of life-history and environmental variability, can provide insights into the extent of adaptive variation. Here, we examined genetic variability among populations of mature European anchovies (n = 531) in the Adriatic (13 samples) and Tyrrhenian seas (2 samples) with neutral and putative non-neutral microsatellite loci. These genetic markers failed to confirm the occurrence of two anchovy species in the Adriatic Sea, as previously postulated. However, we found fine-scale population structure in the Adriatic, especially in northern areas, that was associated with four of the 13 environmental variables tested. Geographic gradients in sea temperature, salinity and dissolved oxygen appear to drive adaptive differences in spawning time and early larval development among populations. Resolving adaptive seascapes in Adriatic anchovies provides a means to understand mechanisms underpinning local adaptation and a basis for optimizing exploitation strategies for sustainable harvests.

Establishing the evolutionary compatibility of potential sources of colonizers for overfished stocks: a population genomics approach

Identifying fish stock structure is fundamental to pinpoint stocks that might contribute colonizers to overfished stocks. However, a stock's potential to contribute to rebuilding hinges on demographic connectivity, a challenging parameter to measure. With genomics as a new tool, fisheries managers can detect signatures of natural selection and thus identify fishing areas likely to contribute evolutionarily compatible colonizers to an overfished area (i.e. colonizers that are not at a fitness disadvantage in the overfished area and able to reproduce at optimal rates). Identifying evolutionarily compatible stocks would help narrow the focus on establishing demographic connectivity where it matters. Here, we genotype 4723 SNPs in 616 orange roughy (Hoplostethus atlanticus) across five fishing areas off the Tasmanian coast in Australia. We ask whether these areas form a single genetic unit, and test for signatures of local adaptation. Results from amova, structure, discriminant analysis of principal components, BAYESASS and isolation by distance suggest that sampled locations are subjected to geneflow amounts that are above what is needed to establish 'drift connectivity'. However, it remains unclear whether there is a single panmictic population or several highly connected populations. Most importantly, we did not find any evidence of local adaptation, suggesting that the examined orange roughy stocks are evolutionarily compatible. The data have helped test an assumption of the orange roughy management programme and to formulate hypotheses regarding stock demographic connectivity. Overall, our results demonstrate the potential of genomics to inform fisheries management, even when evidence for stock structure is sparse.

Genetic evidence for ecological divergence in kokanee salmon

The evolution of locally adapted phenotypes among populations that experience divergent selective pressures is a central mechanism for generating and maintaining biodiversity. Recently, the advent of high-throughput DNA sequencing technology has provided tools for investigating the genetic basis of this process in natural populations of nonmodel organisms. Kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), occurs as two reproductive ecotypes, which differ in spawning habitat (tributaries vs. shorelines); however, outside of the spawning season the two ecotypes co-occur in many lakes and lack diagnostic morphological characteristics. We used restriction site-associated DNA (RAD) sequencing to identify 6145 SNPs and genotype kokanee from multiple spawning sites in Okanagan Lake (British Columbia, Canada). Outlier tests revealed 18 loci putatively under divergent selection between ecotypes, all of which exhibited temporally stable allele frequencies within ecotypes. Six outliers were annotated to sequences in the NCBI database, two of which matched genes associated with early development. There was no evidence for neutral genetic differentiation; however, outlier loci demonstrated significant structure with respect to ecotype and had high assignment accuracy in mixed composition simulations. The absence of neutral structure combined with a small number of highly divergent outlier loci is consistent with theoretical predictions for the early stages of ecological divergence. These outlier loci were then applied to a realistic fisheries scenario in which additional RAD sequencing was used to genotype kokanee collected by trawl in Okanagan Lake, providing preliminary evidence that this approach may be an effective tool for conservation and management.

Diversity and relative abundance of the bacterial pathogen, Flavobacterium spp., infecting reproductive ecotypes of kokanee salmon

Background

Understanding the distribution and abundance of pathogens can provide insight into the evolution and ecology of their host species. Previous research in kokanee, the freshwater form of sockeye salmon (Oncorhynchus nerka), found evidence that populations spawning in streams may experience a greater pathogen load compared with populations that spawn on beaches. In this study we tested for differences in the abundance and diversity of the gram-negative bacteria, Flavobacterium spp., infecting tissues of kokanee in both of these spawning habitats (streams and beaches). Molecular assays were carried out using primers designed to amplify a ~200 nucleotide region of the gene encoding the ATP synthase alpha subunit (AtpA) within the genus Flavobacterium. Using a combination of DNA sequencing and quantitative PCR (qPCR) we compared the diversity and relative abundance of Flavobacterium AtpA amplicons present in DNA extracted from tissue samples of kokanee collected from each spawning habitat.

Results

We identified 10 Flavobacterium AtpA haplotypes among the tissues of stream-spawning kokanee and seven haplotypes among the tissues of beach-spawning kokanee, with only two haplotypes shared between spawning habitats. Haplotypes occurring in the same clade as F. psychrophilum were the most prevalent (92% of all reads, 60% of all haplotypes), and occurred in kokanee from both spawning habitats (streams and beaches). Subsequent qPCR assays did not find any significant difference in the relative abundance of Flavobacterium AtpA amplicons between samples from the different spawning habitats.

Conclusions

We confirmed the presence of Flavobacterium spp. in both spawning habitats and found weak evidence for increased Flavobacterium diversity in kokanee sampled from stream-spawning sites. However, the quantity of Flavobacterium DNA did not differ between spawning habitats. We recommend further study aimed at quantifying pathogen diversity and abundance in population-level samples of kokanee combined with environmental sampling to better understand the ecology of pathogen infection in this species.

Electronic supplementary material

The online version of this article (doi:10.1186/1756-0500-7-778) contains supplementary material, which is available to authorized users.

Time scale matters: genetic analysis does not support adaptation-by-time as the mechanism for adaptive seasonal declines in kokanee reproductive life span

Seasonal declines of fitness-related traits are often attributed to environmental effects or individual-level decisions about reproductive timing and effort, but genetic variation may also play a role. In populations of Pacific salmon (Oncorhynchus spp.), seasonal declines in reproductive life span have been attributed to adaptation-by-time, in which divergent selection for different traits occurs among reproductively isolated temporal components of a population. We evaluated this hypothesis in kokanee (freshwater obligate Oncorhynchus nerka) by testing for temporal genetic structure in neutral and circadian-linked loci. We detected no genetic differences in presumably neutral loci among kokanee with different arrival and maturation dates within a spawning season. Similarly, we detected no temporal genetic structure in OtsClock1b, Omy1009uw, or OmyFbxw11, candidate loci associated with circadian function. The genetic evidence from this study and others indicates a lack of support for adaptation-by-time as an important evolutionary mechanism underlying seasonal declines in reproductive life span and a need for greater consideration of other mechanisms such as time-dependent, adaptive adjustment of reproductive effort.

Genotyping by sequencing resolves shallow population structure to inform conservation of Chinook salmon (Oncorhynchus tshawytscha)

Recent advances in population genomics have made it possible to detect previously unidentified structure, obtain more accurate estimates of demographic parameters, and explore adaptive divergence, potentially revolutionizing the way genetic data are used to manage wild populations. Here, we identified 10 944 single-nucleotide polymorphisms using restriction-site-associated DNA (RAD) sequencing to explore population structure, demography, and adaptive divergence in five populations of Chinook salmon (Oncorhynchus tshawytscha) from western Alaska. Patterns of population structure were similar to those of past studies, but our ability to assign individuals back to their region of origin was greatly improved (>90% accuracy for all populations). We also calculated effective size with and without removing physically linked loci identified from a linkage map, a novel method for nonmodel organisms. Estimates of effective size were generally above 1000 and were biased downward when physically linked loci were not removed. Outlier tests based on genetic differentiation identified 733 loci and three genomic regions under putative selection. These markers and genomic regions are excellent candidates for future research and can be used to create high-resolution panels for genetic monitoring and population assignment. This work demonstrates the utility of genomic data to inform conservation in highly exploited species with shallow population structure.

Finding markers that make a difference: DNA pooling and SNP-arrays identify population informative markers for genetic stock identification

Genetic stock identification (GSI) using molecular markers is an important tool for management of migratory species. Here, we tested a cost-effective alternative to individual genotyping, known as allelotyping, for identification of highly informative SNPs for accurate genetic stock identification. We estimated allele frequencies of 2880 SNPs from DNA pools of 23 Atlantic salmon populations using Illumina SNP-chip. We evaluated the performance of four common strategies (global F ST, pairwise F ST, Delta and outlier approach) for selection of the most informative set of SNPs and tested their effectiveness for GSI compared to random sets of SNP and microsatellite markers. For the majority of cases, SNPs selected using the outlier approach performed best followed by pairwise F ST and Delta methods. Overall, the selection procedure reduced the number of SNPs required for accurate GSI by up to 53% compared with randomly chosen SNPs. However, GSI accuracy was more affected by populations in the ascertainment group rather than the ranking method itself. We demonstrated for the first time the compatibility of different large-scale SNP datasets by compiling the largest population genetic dataset for Atlantic salmon to date. Finally, we showed an excellent performance of our top SNPs on an independent set of populations covering the main European distribution range of Atlantic salmon. Taken together, we demonstrate how combination of DNA pooling and SNP arrays can be applied for conservation and management of salmonids as well as other species.

High degree of genetic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus)

Populations of widespread marine organisms are typically characterized by a low degree of genetic differentiation in neutral genetic markers, but much less is known about differentiation in genes whose functional roles are associated with specific selection regimes. To uncover possible adaptive population divergence and heterogeneous genomic differentiation in marine three-spined sticklebacks (Gasterosteus aculeatus), we used a candidate gene-based genome-scan approach to analyse variability in 138 microsatellite loci located within/close to (<6 kb) functionally important genes in samples collected from ten geographic locations. The degree of genetic differentiation in markers classified as neutral or under balancing selection-as determined with several outlier detection methods-was low (F(ST) = 0.033 or 0.011, respectively), whereas average FST for directionally selected markers was significantly higher (F(ST) = 0.097). Clustering analyses provided support for genomic and geographic heterogeneity in selection: six genetic clusters were identified based on allele frequency differences in the directionally selected loci, whereas four were identified with the neutral loci. Allelic variation in several loci exhibited significant associations with environmental variables, supporting the conjecture that temperature and salinity, but not optic conditions, are important drivers of adaptive divergence among populations. In general, these results suggest that in spite of the high degree of physical connectivity and gene flow as inferred from neutral marker genes, marine stickleback populations are strongly genetically structured in loci associated with functionally relevant genes.

Multi-genetic marker approach and spatio-temporal analysis suggest there is a single panmictic population of swordfish Xiphias gladius in the Indian Ocean

Genetic population structure of swordfish Xiphias gladius was examined based on 2231 individual samples, collected mainly between 2009 and 2010, among three major sampling areas within the Indian Ocean (IO; twelve distinct sites), Atlantic (two sites) and Pacific (one site) Oceans using analysis of nineteen microsatellite loci (n = 2146) and mitochondrial ND2 sequences (n = 2001) data. Sample collection was stratified in time and space in order to investigate the stability of the genetic structure observed with a special focus on the South West Indian Ocean. Significant AMOVA variance was observed for both markers indicating genetic population subdivision was present between oceans. Overall value of F-statistics for ND2 sequences confirmed that Atlantic and Indian Oceans swordfish represent two distinct genetic stocks. Indo-Pacific differentiation was also significant but lower than that observed between Atlantic and Indian Oceans. However, microsatellite F-statistics failed to reveal structure even at the inter-oceanic scale, indicating that resolving power of our microsatellite loci was insufficient for detecting population subdivision. At the scale of the Indian Ocean, results obtained from both markers are consistent with swordfish belonging to a single unique panmictic population. Analyses partitioned by sampling area, season, or sex also failed to identify any clear structure within this ocean. Such large spatial and temporal homogeneity of genetic structure, observed for such a large highly mobile pelagic species, suggests as satisfactory to consider swordfish as a single panmictic population in the Indian Ocean.

Transcriptome-wide comparison of sequence variation in divergent ecotypes of kokanee salmon

Background

High throughput next-generation sequencing technology has enabled the collection of genome-wide sequence data and revolutionized single nucleotide polymorphism (SNP) discovery in a broad range of species. When analyzed within a population genomics framework, SNP-based genotypic data may be used to investigate questions of evolutionary, ecological, and conservation significance in natural populations of non-model organisms. Kokanee salmon are recently diverged freshwater populations of sockeye salmon (Oncorhynchus nerka) that exhibit reproductive ecotypes (stream-spawning and shore-spawning) in lakes throughout western North America and northeast Asia. Current conservation and management strategies may treat these ecotypes as discrete stocks, however their recent divergence and low levels of gene flow make in-season genetic stock identification a challenge. The development of genome-wide SNP markers is an essential step towards fine-scale stock identification, and may enable a direct investigation of the genetic basis of ecotype divergence.

Results

We used pooled cDNA samples from both ecotypes of kokanee to generate 750 million base pairs of transcriptome sequence data. These raw data were assembled into 11,074 high coverage contigs from which we identified 32,699 novel single nucleotide polymorphisms. A subset of these putative SNPs was validated using high-resolution melt analysis and Sanger resequencing to genotype independent samples of kokanee and anadromous sockeye salmon. We also identified a number of contigs that were composed entirely of reads from a single ecotype, which may indicate regions of differential gene expression between the two reproductive ecotypes. In addition, we found some evidence for greater pathogen load among the kokanee sampled in stream-spawning habitats, suggesting a possible evolutionary advantage to shore-spawning that warrants further study.

Conclusions

This study provides novel genomic resources to support population genetic and genomic studies of both kokanee and anadromous sockeye salmon, and has the potential to produce markers capable of fine-scale stock assessment. While this RNAseq approach was successful at identifying a large number of new SNP loci, we found that the frequency of alleles present in the pooled transcriptome data was not an accurate predictor of population allele frequencies.