Are low but statistically significant levels of genetic differentiation in marine fishes 'biologically meaningful'? A case study of coastal Atlantic cod

Patterns of persistence: Genetic and behavioral population complexity of winter flounder amid population declines

Winter flounder Pseudopleuronectes americanus (Walbaum 1792) are a coastal flatfish species of economic and cultural importance that have dwindled to <15, % of their historic abundance in the southern New England/Mid-Atlantic region of the United States, with evidence indicating near-extirpation of certain local populations. This species exhibits intricate behaviors in spawning and migration that contribute to population complexity and resilience. These behaviors encompass full or partial philopatry to natal estuaries, the generation of multiple pulses of larval delivery, and partial migration. The patterns of genetic diversity within and among estuaries and cohorts presented here carry important implications in understanding the susceptibility to demographic shocks, even if the full extent of genetic diversity within and among winter flounder stocks on the US East Coast remains unresolved. Our findings reveal connectivity between estuaries in Long Island, New York, suggesting the potential for genetic rescue of depleted subpopulations. Family reconstruction and relatedness analysis indicate that split cohorts and migration contingents are not the result of genetically distinct lineages. We found no evidence for genetic structure separating these groups, and in some instances, we were able to detect closely related individuals that belonged to different migratory contingents or cohorts. Characterizing the spatial and behavioral organization of this species at the population level is crucial for comprehending its potential for recovery, not only in terms of biomass but also in reinstating the complex population structure that supports resilience. The search for generality in winter flounder spawning and migration behavior remains elusive, but perhaps the lack of generalities within this species is what has allowed it to persist in the face of decades of environmental and anthropogenic stressors.

The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

Genetic divergence at species boundaries of the dolphinfish (Coryphaena hippurus) in the Tropical Eastern Pacific

Background

Marine species constitute commercially important resources, and knowledge about mechanisms that shape phylogeographic patterns and genetic structure provides valuable information for conservation. The dolphinfish, Coryphaena hippurus, is one of the most important species caught in the Tropical Eastern Pacific (TEP). However, the lack of consensus about the existence of genetically differentiated populations in the area has hindered the adoption of management strategies to ensure its viability.

Methods

We assessed genetic variation and phylogeographic structure using two mitochondrial genes and 14 nuclear DNA microsatellite loci. Population genetic tools were used to characterize the spatial distribution of genetic variation of C. hippurus in the TEP, evaluate the extent of connectivity between dolphinfish populations, infer potential barriers to gene flow, and test for signals of contemporary and historical demographic expansions.

Results

Mitochondrial DNA sequences showed genetic homogeneity across locations in the TEP, as well as a strong signal of population expansion dated to the late Pleistocene. In contrast, nuclear microsatellite markers resolved four genetically distinct groups with a remarked genetic differentiation between the most distant locations, at the northern and southern boundaries of the species' range. High mean genetic diversity was found at all localities (Hs = 0.66-0.81). Notwithstanding, positive F _IS and low effective population size (Ne = 77.9-496.4) were also recorded.

Conclusions

The distribution of genetic variation could be related to expansion-contraction cycles following seasonal temperature changes at transitional areas, promoting population subdivisions. However, we cannot rule out the effect of oceanographic dynamics to the observed patterns. Although this marine species remains highly abundant despite commercial exploitation, the low Ne values are of conservation concern and must be considered in fishery management plans.

Eco-evolutionary dynamics of Atlantic cod spatial behavior maintained after the implementation of a marine reserve

The effects of marine reserves on the life history and demography of the protected populations are well-established, typically increasing population density and body size. However, little is known about how marine reserves may alter the behavior of the populations that are the target of protection. In theory, marine reserves can relax selection on spatial behavioral phenotypes that were previously targeted by the fishery and also drive selection in favor of less mobile individuals. In this study, we used acoustic telemetry to monitor the individual spatial behavior of 566 Atlantic cod (Gadus morhua Linnaeus, 1758) moving within a marine reserve and a control site in southern Norway, starting 1 year before the implementation of the marine reserve and lasting up to 9 years after. Following a before-after-control-impact approach, we investigated changes in (1) survival, (2) selection acting on behavioral traits, and (3) mean behavioral phenotypes, after the implementation of the marine reserve. We focused on three behavioral traits commonly used to describe the mobility of aquatic animals: home range size, depth position, and diel vertical migration range. Survival increased after reserve implementation, but contrary to our expectations, it subsequently decreased to preprotection levels after just 3 years. Further, we found no significance in selection patterns acting on any of the three behavioral traits after reserve implementation. Although some changes related to water column use (the tendency to occupy deeper waters) were observed in the marine reserve after 9 years, they cannot unequivocally be attributed to protection. Our results show that survival and behavioral responses to marine reserves in some cases may be more complex than previously anticipated and highlight the need for appropriately scaled management experiments and more integrated approaches to understand the effects of marine protected areas on harvested aquatic species.

Genetic Population Structure of Summer Flounder Paralichthys dentatus using Microsatellite DNA Analysis

Summer flounder Paralichthys dentatus supports one of the most valuable commercial and recreational fisheries along the Atlantic Coast of the U.S. However, in recent decades the management of this species has proven to be one of the most contentious for any exploited marine resource in the region. A coastwide catch quota is imposed annually for summer flounder of which 60% is allocated to the commercial fishery and 40% to the recreational fishery. The allocation is further divided among the individual coastal states from North Carolina to Massachusetts based on their landings in the 1980s. This process, based on political jurisdictions, does not consider the species' biological stock structure. Previous genetic studies (allozyme, mtDNA, and SNPs) provided contradictory results regarding the possible population structure of summer. To address this issue, we used DNA microsatellite analysis at 9 loci to define the coastwide population structure of summer flounder. In total, 1,182 specimens were analyzed from 18 collection sites. Most collections were from the continental shelf during the fall-winter spawning season. These were supplemented with additional samples from inshore waters from North Carolina to Florida, and inshore sites which support significant recreational fisheries at Nantucket Shoals, Massachusetts and Fire Island, New York. The overall level of genetic differentiation in pairwise comparison between collections was very low, mean F ST = 0.001. There was no evidence of genetic differentiation between collections from north and south of Cape Hatteras. Our microsatellite results are consistent with an earlier SNP study which failed to find significant allelic heterogeneity among coastwide collections of summer flounder. However, a subset of pairwise F ST comparisons between some collections proved statistically significant. Furthermore, in STRUCTURE analysis we found evidence of two genetic clusters within the species' northern landings area, however, this finding was not supported by DPAC analysis. We conclude that summer flounder most likely constitute a single population along their entire Atlantic Coast distribution.

Genomic Signatures of Domestication Selection in the Australasian Snapper (Chrysophrys auratus)

Domestication of teleost fish is a recent development, and in most cases started less than 50 years ago. Shedding light on the genomic changes in key economic traits during the domestication process can provide crucial insights into the evolutionary processes involved and help inform selective breeding programmes. Here we report on the recent domestication of a native marine teleost species in New Zealand, the Australasian snapper (Chrysophrys auratus). Specifically, we use genome-wide data from a three-generation pedigree of this species to uncover genetic signatures of domestication selection for growth. Genotyping-By-Sequencing (GBS) was used to generate genome-wide SNP data from a three-generation pedigree to calculate generation-wide averages of FST between every generation pair. The level of differentiation between generations was further investigated using ADMIXTURE analysis and Principal Component Analysis (PCA). After that, genome scans using Bayescan, LFMM and XP-EHH were applied to identify SNP variants under putative selection following selection for growth. Finally, genes near candidate SNP variants were annotated to gain functional insights. Analysis showed that between generations FST values slightly increased as generational time increased. The extent of these changes was small, and both ADMIXTURE analysis and PCA were unable to form clear clusters. Genome scans revealed a number of SNP outliers, indicative of selection, of which a small number overlapped across analyses methods and populations. Genes of interest within proximity of putative selective SNPs were related to biological functions, and revealed an association with growth, immunity, neural development and behaviour, and tumour repression. Even though few genes overlapped between outlier SNP methods, gene functionalities showed greater overlap between methods. While the genetic changes observed were small in most cases, a number of outlier SNPs could be identified, of which some were found by more than one method. Multiple outlier SNPs appeared to be predominately linked to gene functionalities that modulate growth and survival. Ultimately, the results help to shed light on the genomic changes occurring during the early stages of domestication selection in teleost fish species such as snapper, and will provide useful candidates for the ongoing selective breeding in the future of this and related species.

Historical Demographic Processes Dominate Genetic Variation in Ancient Atlantic Cod Mitogenomes

Ancient DNA (aDNA) approaches have been successfully used to infer the long-term impacts of climate change, domestication, and human exploitation in a range of terrestrial species. Nonetheless, studies investigating such impacts using aDNA in marine species are rare. Atlantic cod (Gadus morhua), is an economically important species that has experienced dramatic census population declines during the last century. Here, we investigated 48 ancient mitogenomes from historical specimens obtained from a range of archeological excavations in northern Europe dated up to 6,500 BCE. We compare these mitogenomes to those of 496 modern conspecifics sampled across the North Atlantic Ocean and adjacent seas. Our results confirm earlier observations of high levels of mitogenomic variation and a lack of mutation-drift equilibrium—suggestive of population expansion. Furthermore, our temporal comparison yields no evidence of measurable mitogenomic changes through time. Instead, our results indicate that mitogenomic variation in Atlantic cod reflects past demographic processes driven by major historical events (such as oscillations in sea level) and subsequent gene flow rather than contemporary fluctuations in stock abundance. Our results indicate that historical and contemporaneous anthropogenic pressures such as commercial fisheries have had little impact on mitogenomic diversity in a wide-spread marine species with high gene flow such as Atlantic cod. These observations do not contradict evidence that overfishing has had negative consequences for the abundance of Atlantic cod and the importance of genetic variation in implementing conservation strategies. Instead, these observations imply that any measures toward the demographic recovery of Atlantic cod in the eastern Atlantic, will not be constrained by recent loss of historical mitogenomic variation.

Depletion of coastal predatory fish sub-stocks coincided with the largest sea urchin grazing event observed in the NE Atlantic

In this contribution, we propose fishery driven predator release as the cause for the largest grazing event ever observed in the NE Atlantic. Based on the evolving appreciation of limits to population connectivity, published and previously unpublished data, we discuss whether overfishing caused a grazer bloom of the sea urchin (Strongylocentrotus droebachiensis) resulting in overgrazing of more than 2000 km² kelp (Laminaria hyperborea) forest along Norwegian and Russian coasts during the 1970 s. We show that coastal fisheries likely depleted predatory coastal fish stocks through modernization of fishing methods and fleet. These fish were important predators on urchins and the reduction coincided with the urchin bloom. From this circumstantial evidence, we hypothesize that coastal predatory fish were important in regulating sea urchins, and that a local population dynamics perspective is necessary in management of coastal ecosystems.

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.

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.

Genetic analysis provides insights into species distribution and population structure in East Atlantic horse mackerel (Trachurus trachurus and T. capensis)

Two sister species of horse mackerel (Trachurus trachurus and T. capensis) are described that are intensively harvested in East Atlantic waters. To address long-standing uncertainties as to their respective geographical ranges, overlap and intraspecific population structure this study combined genetic (mitochondrial DNA and microsatellite) analysis and targeted sampling of the hitherto understudied West African coast. mtDNA revealed two reciprocally monophyletic clades corresponding to each species with interspecies nuclear differentiation supported by FST values. The T. trachurus clade was found across the north-east Atlantic down to Ghana but was absent from Angolan and South African samples. The T. capensis clade was found only in South Africa, Angola and a single Ghanaian individual. This pattern suggests that both species may overlap in the waters around Ghana. The potential for cryptic hybridization and/or indiscriminate harvesting of both species in the region is discussed. For T. capensis mtDNA supports high gene flow across the Benguela upwelling system, which fits with the species' ecology. The data add to evidence of a lack of significant genetic structure throughout the range of T. trachurus though the assumption of demographic panmixia is cautioned against. For both species, resolution of stock recruitment heterogeneity relevant to fishery management, as well as potential hybridization, will require more powerful genomic analyses.

Population genetic structure and phylogenetic analysis of gray's pipefish, Halicampus grayi in the South China Sea

BACKGROUND

The evolution of male pregnancy is the most distinctive characteristic of syngnathids, and they were recognized as flagship species for marine conservation. Genus Halicampus is an important branch of syngnathid fishes that has not received the attention it deserves.

OBJECTIVE

To sequence the mitochondrial genome of Halicampus grayi, and investage the genetic structure of its populations.

METHODS

Degenerate primers were designed to amplify the entire mitochondrial genome of H. grayi. The phylogenetic relationship between H. grayi and other syngnathids were conducted using maximum-likelihood method. Population genetic structure of three geographic population of H. grayi were determined using median-joining haplotype network based on COI and Cytb sequences.

RESULTS

The complete mitochondrial genome of Halicampus grayi was assembled into a 17,059 bp circular sequence, which contains 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and 1 D-loop region. The overall base composition of H. grayi is 29.93% A, 29.31% T, 16.23% G and 24.54% C, with a slight A + T rich feature (59.24%). Phylogenetic analysis indicated that H. grayi has a close relationship with Trachyrhamphus serratus. Population genetic analysis revealed a relatively high genetic diversity across different geographic populations of H. grayi, and the results of median-joining haplotype network indicated a lack of structure in populations of H. grayi.

CONCLUSION

The mitogenome of H. grayi will provided important information about the origin and evolution issues of syngnathid fishes, and the high-level genetic diversity detected in their populations will provide insight into the gene flow pattern of marine fishes.

Sex matters: Otolith shape and genomic variation in deacon rockfish (Sebastes diaconus)

Little is known about intraspecific variation within the deacon rockfish (Sebastes diaconus), a recently described species found in the northeast Pacific Ocean. We investigated population structure among fish sampled from two nearshore reefs (Siletz Reef and Seal Rock) and one offshore site (Stonewall Bank) within a <50-km2 area off the Oregon coast. Fish from the three sample sites exhibited small but statistically significant differences based on genetic variation at >15,000 neutral loci, whether analyzed independently or classified into nearshore and offshore groups. Male and females were readily distinguished using genetic data and 92 outlier loci were associated with sex, potentially indicating differential selection between males and females. Morphometric results indicated that there was significant secondary sexual dimorphism in otolith shape, but further sampling is required to disentangle potential confounding influence of age. This study is the first step toward understanding intraspecific variation within the deacon rockfish and the potential management implications. Since differentiation among the three sample sites was small, we consider the results to be suggestive of a single stock. However, future studies should evaluate how the stock is affected by differences in sex, age, and gene flow between the nearshore and offshore environments.

Indications of strong adaptive population genetic structure in albacore tuna (Thunnus alalunga) in the southwest and central Pacific Ocean

Albacore tuna (Thunnus alalunga) has a distinctly complex life history in which juveniles and adults separate geographically but at times inhabit the same spaces sequentially. The species also migrates long distances and presumably experiences varied regimes of physical stress over a lifetime. There are, therefore, many opportunities for population structure to arise based on stochastic differences or environmental factors that promote local adaptation. However, with the extent of mobility consistently demonstrated by tagged individuals, there is also a strong argument for panmixia within an ocean basin. It is important to confirm such assumptions from a population genetics standpoint for this species in particular because albacore is one of the principal market tuna species that sustains massive global fisheries and yet is also a slow-growing temperate tuna. Consequently, we used 1,837 neutral SNP loci and 89 loci under potential selection to analyze population genetic structure among five sample groups collected from the western and central South Pacific. We found no evidence to challenge panmixia at neutral loci, but strong indications of structuring at adaptive loci. One population sample, from French Polynesia in 2004, was particularly differentiated. Unfortunately, the current study cannot infer whether the divergence is geographic or temporal, or possibly caused by sample distribution. We encourage future studies to include potentially adaptive loci and to continue fine-scale observations within an ocean basin, and not to assume genome-wide panmixia.

Possible adverse impact of contaminants on Atlantic cod population dynamics in coastal ecosystems

While many in-laboratory ecotoxicological studies have shown the adverse impact of pollutants to the fitness of an individual, direct evidence from the field on the population dynamics of wildlife animals has been lacking. Here, we provide empirical support for a negative effect of pollution on Atlantic cod (Gadus morhua) population dynamics in coastal waters of Norway by combining unique time series of juvenile cod abundance, body size, environmental concentration of toxic contaminants and a spatially structured population dynamics model. The study shows that mercury concentration might have decreased the reproductive potential of cod in the region despite the general decline in the environmental concentration of mercury, cadmium and hexachlorobenzene since the implementation of national environmental laws. However, some cod populations appeared to be more resistant to mercury pollution than others, and the strength and shape of mercury effect on cod reproductive potential was fjord-specific. Additionally, cod growth rate changed at scales smaller than fjords with a gradient related to the exposure to the open ocean and offshore cod. These spatial differences in life-history traits emphasize the importance of local adaptation in shaping the dynamics of local wildlife populations. Finally, this study highlights the possibility to mitigate pollution effects on natural populations by reducing the overall pollution level, but also reveals that pollution reduction alone is not enough to rebuild local cod populations. Cod population recovery probably requires complementary efforts on fishing regulation and habitat restoration.

Intraspecific DNA contamination distorts subtle population structure in a marine fish: Decontamination of herring samples before restriction-site associated sequencing and its effects on population genetic statistics

Wild specimens are often collected in challenging field conditions, where samples may be contaminated with the DNA of conspecific individuals. This contamination can result in false genotype calls, which are difficult to detect, but may also cause inaccurate estimates of heterozygosity, allele frequencies and genetic differentiation. Marine broadcast spawners are especially problematic, because population genetic differentiation is low and samples are often collected in bulk and sometimes from active spawning aggregations. Here, we used contaminated and clean Pacific herring (Clupea pallasi) samples to test (a) the efficacy of bleach decontamination, (b) the effect of decontamination on RAD genotypes and (c) the consequences of contaminated samples on population genetic analyses. We collected fin tissue samples from actively spawning (and thus contaminated) wild herring and nonspawning (uncontaminated) herring. Samples were soaked for 10 min in bleach or left untreated, and extracted DNA was used to prepare DNA libraries using a restriction site-associated DNA (RAD) approach. Our results demonstrate that intraspecific DNA contamination affects patterns of individual and population variability, causes an excess of heterozygotes and biases estimates of population structure. Bleach decontamination was effective at removing intraspecific DNA contamination and compatible with RAD sequencing, producing high-quality sequences, reproducible genotypes and low levels of missing data. Although sperm contamination may be specific to broadcast spawners, intraspecific contamination of samples may be common and difficult to detect from high-throughput sequencing data and can impact downstream analyses.

Adaptive genetic variation underlies biocomplexity of Atlantic Cod in the Gulf of Maine and on Georges Bank

Atlantic cod (Gadus morhua) populations in the Gulf of Maine (GoM) are at a fraction of their historical abundance, creating economic hardships for fishermen and putting at risk the genetic diversity of the remaining populations. An understanding of the biocomplexity among GoM populations will allow for adaptive genetic diversity to be conserved to maximize the evolutionary potential and resilience of the fishery in a rapidly changing environment. We used restriction-site-associated DNA sequencing (RADseq) to characterize the population structure and adaptive genetic diversity of five spawning aggregations from the western GoM and Georges Bank. We also analyzed cod caught in the eastern GoM, an under-sampled area where spawning aggregations have been extirpated. Using 3,128 single nucleotide polymorphisms (SNPs), we confirmed the existence of three genetically separable spawning groups: (1) winter spawning cod from the western GoM, (2) spring spawning cod, also from the western GoM, and (3) Georges Bank cod. Non-spawning cod from the eastern GoM could not be decisively linked to either of the three spawning groups and may represent a unique component of the resource, a mixed sample, or cod from other unsampled source populations. The genetic differentiation among the three major spawning groups was primarily driven by loci putatively under selection, particularly loci in regions known to contain genomic inversions on linkage groups (LG) 7 and 12. These LGs have been found to be linked to thermal regime in cod across the Atlantic, and so it is possible that variation in timing of spawning in western GoM cod has resulted in temperature-driven adaptive divergence. This complex population structure and adaptive genetic differentiation could be crucial to ensuring the long-term productivity and resilience of the cod fishery, and so it should be considered in future management plans.

Disentangling structural genomic and behavioural barriers in a sea of connectivity

Genetic divergence among populations arises through natural selection or drift and is counteracted by connectivity and gene flow. In sympatric populations, isolating mechanisms are thus needed to limit the homogenizing effects of gene flow to allow for adaptation and speciation. Chromosomal inversions act as an important mechanism maintaining isolating barriers, yet their role in sympatric populations and divergence with gene flow is not entirely understood. Here, we revisit the question of whether inversions play a role in the divergence of connected populations of the marine fish Atlantic cod (Gadus morhua), by exploring a unique data set combining whole‐genome sequencing data and behavioural data obtained with acoustic telemetry. Within a confined fjord environment, we find three genetically differentiated Atlantic cod types belonging to the oceanic North Sea population, the western Baltic population and a local fjord‐type cod. Continuous behavioural tracking over 4 year revealed temporally stable sympatry of these types within the fjord. Despite overall weak genetic differentiation consistent with high levels of gene flow, we detected significant frequency shifts of three previously identified inversions, indicating an adaptive barrier to gene flow. In addition, behavioural data indicated that North Sea cod and individuals homozygous for the LG12 inversion had lower fitness in the fjord environment. However, North Sea and fjord‐type cod also occupy different depths, possibly contributing to prezygotic reproductive isolation and representing a behavioural barrier to gene flow. Our results provide the first insights into a complex interplay of genomic and behavioural isolating barriers in Atlantic cod and establish a new model system towards an understanding of the role of genomic structural variants in adaptation and diversification.

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.

Fine-scale population differences in Atlantic cod reproductive success: A potential mechanism for ecological speciation in a marine fish

Successful resource-management and conservation outcomes ideally depend on matching the spatial scales of population demography, local adaptation, and threat mitigation. For marine fish with high dispersal capabilities, this remains a fundamental challenge. Based on daily parentage assignments of more than 4,000 offspring, we document fine-scaled temporal differences in individual reproductive success for two spatially adjacent (<10 km) populations of a broadcast-spawning marine fish. Distinguished by differences in genetics and life history, Atlantic cod (Gadus morhua) from inner- and outer-fjord populations were allowed to compete for mating and reproductive opportunities. After accounting for phenotypic variability in several traits, reproductive success of outer-fjord cod was significantly lower than that of inner-fjord cod. This finding, given that genomically different cod ecotypes inhabit inner- and outer-fjord waters, raises the intriguing hypothesis that the populations might be diverging because of ecological speciation. Individual reproductive success, skewed within both sexes (more so among males), was positively affected by body size, which also influenced the timing of reproduction, larger individuals spawning later among females but earlier among males. Our work suggests that spatial mismatches between management and biological units exist in marine fishes and that studies of reproductive interactions between putative populations or ecotypes can provide an informative basis on which determination of the scale of local adaptation can be ascertained.

Stable coexistence of genetically divergent Atlantic cod ecotypes at multiple spatial scales

Coexistence in the same habitat of closely related yet genetically different populations is a phenomenon that challenges our understanding of local population structure and adaptation. Identifying the underlying mechanisms for such coexistence can yield new insight into adaptive evolution, diversification and the potential for organisms to adapt and persist in response to a changing environment. Recent studies have documented cryptic, sympatric populations of Atlantic cod (Gadus morhua) in coastal areas. We analysed genetic origin of 6,483 individual cod sampled annually over 14 years from 125 locations along the Norwegian Skagerrak coast and document stable coexistence of two genetically divergent Atlantic cod ecotypes throughout the study area and study period. A "fjord" ecotype dominated in numbers deep inside fjords while a "North Sea" ecotype was the only type found in offshore North Sea. Both ecotypes coexisted in similar proportions throughout coastal habitats at all spatial scales. The size-at-age of the North Sea ecotype on average exceeded that of the fjord ecotype by 20% in length and 80% in weight across all habitats. Different growth and size among individuals of the two types might be one of several ecologically significant variables that allow for stable coexistence of closely related populations within the same habitat. Management plans, biodiversity initiatives and other mitigation strategies that do not account for the mixture of species ecotypes are unlikely to meet objectives related to the sustainability of fish and fisheries.

Population assignment and local adaptation along an isolation-by-distance gradient in Pacific cod (Gadus macrocephalus)

The discernment of populations as management units is a fundamental prerequisite for sustainable exploitation of species. A lack of clear stock boundaries complicates not only the identification of spatial management units, but also the assessment of mixed fisheries by population assignment and mixed stock analysis. Many marine species, such as Pacific cod, are characterized by isolation by distance, showing significant differentiation but no clear stock boundaries. Here, we used restriction-site-associated DNA (RAD) sequencing to investigate population structure and assess power to genetically assign Pacific cod to putative populations of origin. Samples were collected across the species range in the eastern Pacific Ocean, from the Salish Sea to the Aleutian Islands. A total of 6,425 putative biallelic single nucleotide polymorphisms were identified from 276 individuals. We found a strong isolation-by-distance signal along coastlines that mirrored previous microsatellite results and pronounced genetic differentiation between coastal samples and those from the inland waters of the Salish Sea, with no evidence for hybridization between these two populations. Individual assignment success based on two methods was high overall (≥84%) but decreased from south to north. Assignment to geographic location of origin also was successful, with average distance between capture and assignment location of 220 km. Outlier analyses identified more loci potentially under selection along the coast than between Salish Sea and coastal samples, suggesting more diverse adaptation to latitudinal environmental factors than inshore vs. offshore environments. Our results confirm previous observations of sharp genetic differentiation of the Salish Sea population and isolation by distance along the coast, but also highlight the feasibility of using modern genomic techniques to inform stock boundaries and fisheries management in a low FST marine species.

Weak population structure of the Spot-tail shark Carcharhinus sorrah and the Blacktip shark C. limbatus along the coasts of the Arabian Peninsula, Pakistan, and South Africa

The increase in demand for shark meat and fins has placed shark populations worldwide under high fishing pressure. In the Arabian region, the spot-tail shark Carcharhinus sorrah and the Blacktip shark Carcharhinus limbatus are among the most exploited species. In this study, we investigated the population genetic structure of C. sorrah (n = 327) along the coasts of the Arabian Peninsula and of C. limbatus (n = 525) along the Arabian coasts, Pakistan, and KwaZulu-Natal, South Africa, using microsatellite markers (15 and 11 loci, respectively). Our findings support weak population structure in both species. Carcharhinus sorrah exhibited a fine structure, subdividing the area into three groups. The first group comprises all samples from Bahrain, the second from the UAE and Yemen, and the third from Oman. Similarly, C. limbatus exhibited population subdivision into three groups. The first group, comprising samples from Bahrain and Kuwait, was highly differentiated from the second and third groups, comprising samples from Oman, Pakistan, the UAE, and Yemen; and South Africa and the Saudi Arabian Red Sea, respectively. Population divisions were supported by pairwise F ST values and discriminant analysis of principal components (DAPC), but not by STRUCTURE. We suggest that the mostly low but significant pairwise F ST values in our study are suggestive of fine population structure, which is possibly attributable to behavioral traits such as residency in C. sorrah and site fidelity and philopatry in C. limbatus. However, for all samples obtained from the northern parts of the Gulf (Bahrain and/or Kuwait) in both species, the higher but significant pairwise F ST values could possibly be a result of founder effects during the Tethys Sea closure. Based on DAPC and F ST results, we suggest each population to be treated as independent management unit, as conservation concerns emerge.

Contemporary ancestor? Adaptive divergence from standing genetic variation in Pacific marine threespine stickleback

Background

Populations that have repeatedly colonized novel environments are useful for studying the role of ecology in adaptive divergence – particularly if some individuals persist in the ancestral habitat. Such “contemporary ancestors” can be used to demonstrate the effects of selection by comparing phenotypic and genetic divergence between the derived population and their extant ancestors. However, evolution and demography in these “contemporary ancestors” can complicate inferences about the source (standing genetic variation, de novo mutation) and pace of adaptive divergence. Marine threespine stickleback (Gasterosteus aculeatus) have colonized freshwater environments along the Pacific coast of North America, but have also persisted in the marine environment. To what extent are marine stickleback good proxies of the ancestral condition?

Results

We sequenced > 5800 variant loci in over 250 marine stickleback from eight locations extending from Alaska to California, and phenotyped them for platedness and body shape. Pairwise F_ST varied from 0.02 to 0.18. Stickleback were divided into five genetic clusters, with a single cluster comprising stickleback from Washington to Alaska. Plate number, Eda, body shape, and candidate loci showed evidence of being under selection in the marine environment. Comparisons to a freshwater population demonstrated that candidate loci for freshwater adaptation varied depending on the choice of marine populations.

Conclusions

Marine stickleback are structured into phenotypically and genetically distinct populations that have been evolving as freshwater stickleback evolved. This variation complicates their usefulness as proxies of the ancestors of freshwater populations. Lessons from stickleback may be applied to other “contemporary ancestor”-derived population studies.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1228-8) contains supplementary material, which is available to authorized users.

Female-biased dispersal and non-random gene flow of MC1R variants do not result in a migration load in barn owls

Non-random gene flow is a widely neglected force in evolution and ecology. This genotype-dependent dispersal is difficult to assess, yet can impact the genetic variation of natural populations and their fitness. In this work, we demonstrate a high immigration rate of barn owls (Tyto alba) inside a Swiss population surveyed during 15 years. Using ten microsatellite loci as an indirect method to characterize dispersal, two-third of the genetic tests failed to detect a female-biased dispersal, and Monte Carlo simulations confirmed a low statistical power to detect sex-biased dispersal in case of high dispersal rate of both sexes. The capture-recapture data revealed a female-biased dispersal associated with an excess of heterozygote for the melanocortin-1 receptor gene (MC1R), which is responsible for their ventral rufous coloration. Thus, female homozygotes for the MC1R_WHITE allele might be negatively selected during dispersal. Despite the higher immigration of females that are heterozygote at MC1R, non-random gene flow should not lead to a migration load regarding this gene because we did not detect an effect of MC1R on survival and reproductive success in our local population. The present study highlights the usefulness of using multiple methods to correctly decrypt dispersal and gene flow. Moreover, despite theoretical expectations, we show that non-random dispersal of particular genotypes does not necessarily lead to migration load in recipient populations.

Analysis of coastal cod (Gadus morhua L.) sampled on spawning sites reveals a genetic gradient throughout Norway's coastline

BACKGROUND

Atlantic cod (Gadus morhua L.) has formed the basis of many economically significant fisheries in the North Atlantic, and is one of the best studied marine fishes, but a legacy of overexploitation has depleted populations and collapsed fisheries in several regions. Previous studies have identified considerable population genetic structure for Atlantic cod. However, within Norway, which is the country with the largest remaining catch in the Atlantic, the population genetic structure of coastal cod (NCC) along the entire coastline has not yet been investigated. We sampled > 4000 cod from 55 spawning sites. All fish were genotyped with 6 microsatellite markers and Pan I (Dataset 1). A sub-set of the samples (1295 fish from 17 locations) were also genotyped with an additional 9 microsatellites (Dataset 2). Otoliths were read in order to exclude North East Arctic Cod (NEAC) from the analyses, as and where appropriate.

RESULTS

We found no difference in genetic diversity, measured as number of alleles, allelic richness, heterozygosity nor effective population sizes, in the north-south gradient. In both data sets, weak but significant population genetic structure was revealed (Dataset 1: global FST = 0.008, P < 0.0001. Dataset 2: global FST = 0.004, P < 0.0001). While no clear genetic groups were identified, genetic differentiation increased among geographically-distinct samples. Although the locus Gmo132 was identified as a candidate for positive selection, possibly through linkage with a genomic region under selection, overall trends remained when this locus was excluded from the analyses. The most common allele in loci Gmo132 and Gmo34 showed a marked frequency change in the north-south gradient, increasing towards the frequency observed in NEAC in the north.

CONCLUSION

We conclude that Norwegian coastal cod displays significant population genetic structure throughout its entire range, that follows a trend of isolation by distance. Furthermore, we suggest that a gradient of genetic introgression between NEAC and NCC contributes to the observed population genetic structure. The current management regime for coastal cod in Norway, dividing it into two stocks at 62°N, represents a simplification of the level of genetic connectivity among coastal cod in Norway, and needs revision.

Inferring genetic connectivity in real populations, exemplified by coastal and oceanic Atlantic cod

Estimates of migration are important for understanding the dynamics of natural populations. A statistic known as F_ST is often used to measure levels of genetic differentiation among natural populations. Equations that translate F_ST into estimates of migration are based on “ideal” populations, which are subject to many simplifying assumptions compared with real populations. Therefore, theoretical estimates of migration might not be realistic. We modeled populations of Atlantic cod in the North Sea and the adjacent Skagerrak region to compare how migration is related to the complexities of real populations, and how actual migration compares with predictions based on theory. Results are intended to help apply population genetic theory to practical situations.

Genetic data are commonly used to estimate connectivity between putative populations, but translating them to demographic dispersal rates is complicated. Theoretical equations that infer a migration rate based on the genetic estimator F_ST, such as Wright’s equation, F_ST ≈ 1/(4N_em + 1), make assumptions that do not apply to most real populations. How complexities inherent to real populations affect migration was exemplified by Atlantic cod in the North Sea and Skagerrak and was examined within an age-structured model that incorporated genetic markers. Migration was determined under various scenarios by varying the number of simulated migrants until the mean simulated level of genetic differentiation matched a fixed level of genetic differentiation equal to empirical estimates. Parameters that decreased the N_e/N_t ratio (where N_e is the effective and N_t is the total population size), such as high fishing mortality and high fishing gear selectivity, increased the number of migrants required to achieve empirical levels of genetic differentiation. Higher maturity-at-age and lower selectivity increased N_e/N_t and decreased migration when genetic differentiation was fixed. Changes in natural mortality, fishing gear selectivity, and maturity-at-age within expected limits had a moderate effect on migration when genetic differentiation was held constant. Changes in population size had the greatest effect on the number of migrants to achieve fixed levels of F_ST, particularly when genetic differentiation was low, F_ST ≈ 10⁻³. Highly variable migration patterns, compared with constant migration, resulted in higher variance in genetic differentiation and higher extreme values. Results are compared with and provide insight into the use of theoretical equations to estimate migration among real populations.

Cryptic species and parallel genetic structuring in Lethrinid fish: Implications for conservation and management in the southwest Indian Ocean

Analysis of genetic variation can provide insights into ecological and evolutionary diversification which, for commercially harvested species, can also be relevant to the implementation of spatial management strategies and sustainability. In comparison with other marine biodiversity hot spots, there has been less genetic research on the fauna of the southwest Indian Ocean (SWIO). This is epitomized by the lack of information for lethrinid fish, which support socioeconomically important fisheries in the region. This study combines comparative phylogeographic and population genetic analyses with ecological niche modeling to investigate historical and contemporary population dynamics of two species of emperor fish (Lethrinus mahsena and Lethrinus harak) across the SWIO. Both species shared similarly shallow phylogeographic patterns and modeled historical (LGM) habitat occupancies. For both species, allele frequency and kinship analyses of microsatellite variation revealed highly significant structure with no clear geographical pattern and nonrandom genetic relatedness among individuals within samples. The genetic patterns for both species indicate recurrent processes within the region that prevent genetic mixing, at least on timescales of interest to fishery managers, and the potential roles of recruitment variability and population isolation are discussed in light of biological and environmental information. This consistency in both historical and recurrent population processes indicates that the use of model species may be valuable in management initiatives with finite resources to predict population structure, at least in cases wherein biogeographic and ecological differences between taxa are minimized. Paradoxically, mtDNA sequencing and microsatellite analysis of samples from the Seychelles revealed a potential cryptic species occurring in sympatry with, and seemingly morphologically identical to, L. mahsena. BLAST results point to the likely misidentification of species and incongruence between voucher specimens, DNA barcodes, and taxonomy within the group, which highlights the utility and necessity of genetic approaches to characterize baseline biodiversity in the region before such model-based methods are employed.

Ancient DNA reveals the Arctic origin of Viking Age cod from Haithabu, Germany

Knowledge of the range and chronology of historic trade and long-distance transport of natural resources is essential for determining the impacts of past human activities on marine environments. However, the specific biological sources of imported fauna are often difficult to identify, in particular if species have a wide spatial distribution and lack clear osteological or isotopic differentiation between populations. Here, we report that ancient fish-bone remains, despite being porous, brittle, and light, provide an excellent source of endogenous DNA (15-46%) of sufficient quality for whole-genome reconstruction. By comparing ancient sequence data to that of modern specimens, we determine the biological origin of 15 Viking Age (800-1066 CE) and subsequent medieval (1066-1280 CE) Atlantic cod (Gadus morhua) specimens from excavation sites in Germany, Norway, and the United Kingdom. Archaeological context indicates that one of these sites was a fishing settlement for the procurement of local catches, whereas the other localities were centers of trade. Fish from the trade sites show a mixed ancestry and are statistically differentiated from local fish populations. Moreover, Viking Age samples from Haithabu, Germany, are traced back to the North East Arctic Atlantic cod population that has supported the Lofoten fisheries of Norway for centuries. Our results resolve a long-standing controversial hypothesis and indicate that the marine resources of the North Atlantic Ocean were used to sustain an international demand for protein as far back as the Viking Age.

Reduced Single Nucleotide Polymorphism Panels for Assigning Atlantic Albacore and Bay of Biscay Anchovy Individuals to Their Geographic Origin: Toward Sustainable Fishery Management

There is an increasing trend upon adding a detailed description of the origin of seafood products driven by a general interest in the implementation of sustainable fishery management plans for the conservation of marine ecosystems. North Atlantic albacore ("Bonito del Norte con Eusko Label") and Bay of Biscay anchovy ("Anchoa del Cantábrico") are two commercially important fish populations with high economical value and vulnerable to commercial fraud. This fact, together with the overexploited situation of these two populations, makes it necessary to develop a tool to identify individual origin and to detect commercial fraud. In the present study, we have developed and validated a traceability tool consisting of reduced panels of gene-associated single nucleotide polymorphisms (SNPs) suitable for assigning individuals of two species to their origin with unprecedented accuracy levels. Only 48 SNPs are necessary to assign 81.1% albacore and 93.4% anchovy individuals with 100% accuracy to their geographic origin. The total accuracy of the results demonstrates how gene-associated SNPs can revolutionize food traceability. Gene-associated SNP panels are not of mere commercial interest, but they also can result in a positive impact on sustainability of marine ecosystems through conservation of fish populations through establishing a more effective and sustainable fishery management framework and contributing to the prevention of falsified labeling.

Fine-scale population dynamics in a marine fish species inferred from dynamic state-space models

Identifying the spatial scale of population structuring is critical for the conservation of natural populations and for drawing accurate ecological inferences. However, population studies often use spatially aggregated data to draw inferences about population trends and drivers, potentially masking ecologically relevant population sub-structure and dynamics. The goals of this study were to investigate how population dynamics models with and without spatial structure affect inferences on population trends and the identification of intrinsic drivers of population dynamics (e.g. density dependence). Specifically, we developed dynamic, age-structured, state-space models to test different hypotheses regarding the spatial structure of a population complex of coastal Atlantic cod (Gadus morhua). Data were from a 93-year survey of juvenile (age 0 and 1) cod sampled along >200 km of the Norwegian Skagerrak coast. We compared two models: one which assumes all sampled cod belong to one larger population, and a second which assumes that each fjord contains a unique population with locally determined dynamics. Using the best supported model, we then reconstructed the historical spatial and temporal dynamics of Skagerrak coastal cod. Cross-validation showed that the spatially structured model with local dynamics had better predictive ability. Furthermore, posterior predictive checks showed that a model which assumes one homogeneous population failed to capture the spatial correlation pattern present in the survey data. The spatially structured model indicated that population trends differed markedly among fjords, as did estimates of population parameters including density-dependent survival. Recent biomass was estimated to be at a near-record low all along the coast, but the finer scale model indicated that the decline occurred at different times in different regions. Warm temperatures were associated with poor recruitment, but local changes in habitat and fishing pressure may have played a role in driving local dynamics. More generally, we demonstrated how state-space models can be used to test evidence for population spatial structure based on survey time-series data. Our study shows the importance of considering spatially structured dynamics, as the inferences from such an approach can lead to a different ecological understanding of the drivers of population declines, and fundamentally different management actions to restore populations.

Swept away: ocean currents and seascape features influence genetic structure across the 18,000 Km Indo-Pacific distribution of a marine invertebrate, the black-lip pearl oyster Pinctada margaritifera

BACKGROUND

Genetic structure in many widely-distributed broadcast spawning marine invertebrates remains poorly understood, posing substantial challenges for their fishery management, conservation and aquaculture. Under the Core-Periphery Hypothesis (CPH), genetic diversity is expected to be highest at the centre of a species' distribution, progressively decreasing with increased differentiation towards outer range limits, as populations become increasingly isolated, fragmented and locally adapted. The unique life history characteristics of many marine invertebrates such as high dispersal rates, stochastic survival and variable recruitment are also likely to influence how populations are organised. To examine the microevolutionary forces influencing population structure, connectivity and adaptive variation in a highly-dispersive bivalve, populations of the black-lip pearl oyster Pinctada margaritifera were examined across its ~18,000 km Indo-Pacific distribution.

RESULTS

Analyses utilising 9,624 genome-wide SNPs and 580 oysters, discovered differing patterns of significant and substantial broad-scale genetic structure between the Indian and Pacific Ocean basins. Indian Ocean populations were markedly divergent (F st = 0.2534-0.4177, p < 0.001), compared to Pacific Ocean oysters, where basin-wide gene flow was much higher (F st = 0.0007-0.1090, p < 0.001). Partitioning of genetic diversity (hierarchical AMOVA) attributed 18.1% of variance between ocean basins, whereas greater proportions were resolved within samples and populations (45.8% and 35.7% respectively). Visualisation of population structure at selectively neutral loci resolved three and five discrete genetic clusters for the Indian and Pacific Oceans respectively. Evaluation of genetic structure at adaptive loci for Pacific populations (89 SNPs under directional selection; F st = 0.1012-0.4371, FDR = 0.05), revealed five clusters identical to those detected at neutral SNPs, suggesting environmental heterogeneity within the Pacific. Patterns of structure and connectivity were supported by Mantel tests of isolation by distance (IBD) and independent hydrodynamic particle dispersal simulations.

CONCLUSIONS

It is evident that genetic structure and connectivity across the natural range of P. margaritifera is highly complex, and produced by the interaction of ocean currents, IBD and seascape features at a broad scale, together with habitat geomorphology and local adaptation at regional levels. Overall population organisation is far more elaborate than generalised CPH predictions, however valuable insights for regional fishery management, and a greater understanding of range-wide genetic structure in a highly-dispersive marine invertebrate have been gained.

Is the ballan wrasse (Labrus bergylta) two species? Genetic analysis reveals within-species divergence associated with plain and spotted morphotype frequencies

The ballan wrasse (Labrus bergylta) is a marine fish belonging to the family Labridae characterized by 2 main morphotypes that occur in sympatry: spotty and plain. Previous studies have revealed differences in their life‐history traits, such as growth and maturation; however, the genetic relationship between forms is presently unknown. Using 20 recently developed microsatellite markers, we conducted a genetic analysis of 41 and 48 spotty and plain ballan wrasse collected in Galicia (northwest Spain). The 2 morphotypes displayed highly significant genetic differences to each other (F_ST = 0.018, P < 0.0001). A similar degree of genetic differentiation (F_ST = 0.025, P < 0.0001) was shown using the STRUCTURE clustering approach with no priors at K = 2. In this case, the frequency of spotty and plain morphotypes was significantly different (χ² = 9.46, P = 0.002). It is concluded that there is significant genetic heterogeneity within this species, which appears to be highly associated with the spotty and plain forms, but not completely explained by them. Given the previously demonstrated biological differences between morphotypes, and the present genetic analyses, we speculate about the convenience of a taxonomic re‐evaluation of this species.

A Parallel Population Genomic and Hydrodynamic Approach to Fishery Management of Highly-Dispersive Marine Invertebrates: The Case of the Fijian Black-Lip Pearl Oyster Pinctada margaritifera

Fishery management and conservation of marine species increasingly relies on genetic data to delineate biologically relevant stock boundaries. Unfortunately for high gene flow species which may display low, but statistically significant population structure, there is no clear consensus on the level of differentiation required to resolve distinct stocks. The use of fine-scale neutral and adaptive variation, considered together with environmental data can offer additional insights to this problem. Genome-wide genetic data (4,123 SNPs), together with an independent hydrodynamic particle dispersal model were used to inform farm and fishery management in the Fijian black-lip pearl oyster Pinctada margaritifera, where comprehensive fishery management is lacking, and the sustainability of exploitation uncertain. Weak fine-scale patterns of population structure were detected, indicative of broad-scale panmixia among wild oysters, while a hatchery-sourced farmed population exhibited a higher degree of genetic divergence (Fst = 0.0850-0.102). This hatchery-produced population had also experienced a bottleneck (NeLD = 5.1; 95% C.I. = [5.1-5.3]); compared to infinite NeLD estimates for all wild oysters. Simulation of larval transport pathways confirmed the existence of broad-scale mixture by surface ocean currents, correlating well with fine-scale patterns of population structuring. Fst outlier tests failed to detect large numbers of loci supportive of selection, with 2-5 directional outlier SNPs identified (average Fst = 0.116). The lack of biologically significant population genetic structure, absence of evidence for local adaptation and larval dispersal simulation, all indicate the existence of a single genetic stock of P. margaritifera in the Fiji Islands. This approach using independent genomic and oceanographic tools has allowed fundamental insights into stock structure in this species, with transferability to other highly-dispersive marine taxa for their conservation and management.

Using dense locality sampling resolves the subtle genetic population structure of the dispersive fish species Plecoglossus altivelis

In dispersive species with continuous distributions, genetic differentiation between local populations is often absent or subtle and thus difficult to detect. To incorporate such subtle differentiation into management plans, it may be essential to analyse many samples from many localities using adequate numbers of high-resolution genetic markers. Here, we evaluated the usefulness of dense locality sampling in resolving genetic population structure in the ayu (Plecoglossus altivelis), a dispersive fish important in Japanese inland fisheries. Genetic variability in, and differentiation between, ayu populations around the Japan-Ryukyu Archipelago were investigated in 4746 individuals collected from 120 localities by genotyping 12 microsatellite markers. These individuals represented the two subspecies of ayu, namely the Ryukyuan subspecies (Plecoglossus altivelis ryukyuensis) and both amphidromous and landlocked forms of the nominotypical subspecies (P. a. altivelis) along the archipelago. We successfully detected an absence of genetic differentiation within the landlocked form and subtle but significant differentiation and clear geographic patterns of genetic variation among populations of the amphidromous form, which had been considered genetically homogeneous. This suggests that dense locality sampling effectively resolves subtle differences in genetic population structure, reducing stochastic deviation in the detection of genetic differentiation and geographic patterns in local populations of this dispersive species. Resampling analyses based on empirical data sets clearly demonstrate the effectiveness of increasing the number of locality samples for stable and reliable estimations of genetic fixation indices. The genetic population structure observed within the amphidromous form provides useful information for identifying management or conservation units in ayu.

"Islands of Divergence" in the Atlantic Cod Genome Represent Polymorphic Chromosomal Rearrangements

In several species genetic differentiation across environmental gradients or between geographically separate populations has been reported to center at "genomic islands of divergence," resulting in heterogeneous differentiation patterns across genomes. Here, genomic regions of elevated divergence were observed on three chromosomes of the highly mobile fish Atlantic cod (Gadus morhua) within geographically fine-scaled coastal areas. The "genomic islands" extended at least 5, 9.5, and 13 megabases on linkage groups 2, 7, and 12, respectively, and coincided with large blocks of linkage disequilibrium. For each of these three chromosomes, pairs of segregating, highly divergent alleles were identified, with little or no gene exchange between them. These patterns of recombination and divergence mirror genomic signatures previously described for large polymorphic inversions, which have been shown to repress recombination across extensive chromosomal segments. The lack of genetic exchange permits divergence between noninverted and inverted chromosomes in spite of gene flow. For the rearrangements on linkage groups 2 and 12, allelic frequency shifts between coastal and oceanic environments suggest a role in ecological adaptation, in agreement with recently reported associations between molecular variation within these genomic regions and temperature, oxygen, and salinity levels. Elevated genetic differentiation in these genomic regions has previously been described on both sides of the Atlantic Ocean, and we therefore suggest that these polymorphisms are involved in adaptive divergence across the species distributional range.

Harvest Pressure on Coastal Atlantic Cod (Gadus morhua) from Recreational Fishing Relative to Commercial Fishing Assessed from Tag-Recovery Data

Marine recreational fishing is a popular outdoor activity. However, knowledge about the magnitude of recreational catches relative to commercial catches in coastal fisheries is generally sparse. Coastal Atlantic cod (Gadus morhua) is a target species for recreational fishers in the North Atlantic. In Norway, recreational fishers are allowed to use a variety of traps and nets as well as long-line and rod and line when fishing for cod. From 2005 to 2013, 9729 cod (mean size: 40 cm, range: 15–93 cm) were tagged and released in coastal Skagerrak, southeast Norway. Both high-reward (NOK 500) and low-reward tags (NOK 50) were used in this study. Because some harvested fish (even those posting high-reward tags) may go unreported by fishers, reporting rates were estimated from mark-recovery models that incorporate detection parameters in their structure, in addition to survival and mortality estimates. During 2005 to 2013, a total of 1707 tagged cod were recovered and reported by fishers. We estimate the overall annual survival to be 33% (SE 1.5). Recreational rod and line fishing were responsible for 33.7% (SE 2.4) of total mortality, followed by commercial fisheries (15.1% SE 0.8) and recreational fixed gear (6.8% SE 0.4). Natural mortality was 44.4% (SE 2.5) of total mortality. Our findings suggest that recreational fishing—rod and line fishing in particular—is responsible for a substantial part of fishing mortality exerted on coastal cod in southern Norway.

Population structure in Atlantic cod in the eastern North Sea-Skagerrak-Kattegat: early life stage dispersal and adult migration

Background

In marine fish species, where pelagic egg and larvae drift with ocean currents, population structure has been suggested to be maintained by larval retention due to hydrographic structuring and by homing of adult fish to natal areas. Whilst natal homing of adults has been demonstrated for anadromous and coral reef fishes, there are few documented examples of philopatric migration in temperate marine fish species.

Results

Here, we demonstrate temporally stable genetic differentiation among spawning populations of Atlantic cod (Gadus morhua L.), and present genetic and behavioural evidence for larval drift and philopatric migration in the eastern North Sea-Skagerrak-Kattegat area. We show that juvenile cod collected in the eastern Skagerrak and central Kattegat are genetically similar to cod from offshore spawning areas in the eastern North Sea. Genetic assignment of individual 2–5 year old fish indicates that cod residing at, or migrating towards, spawning areas in Kattegat and the North Sea display philopatric behaviours.

Conclusions

Together these findings suggest a loop between spawning, larval drift and adult return-migrations to spawning areas and underlines that both oceanographic processes and migratory behaviour in the adult phase may be important for stock separation and integrity in marine temperate fishes such as Atlantic cod.

Electronic supplementary material

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

The Story of a Hitchhiker: Population Genetic Patterns in the Invasive Barnacle Balanus(Amphibalanus) improvisus Darwin 1854

Understanding the ecological and evolutionary forces that determine the genetic structure and spread of invasive species is a key component of invasion biology. The bay barnacle, Balanus improvisus (= Amphibalanus improvisus), is one of the most successful aquatic invaders worldwide, and is characterised by broad environmental tolerance. Although the species can spread through natural larval dispersal, human-mediated transport through (primarily) shipping has almost certainly contributed to the current global distribution of this species. Despite its worldwide distribution, little is known about the phylogeography of this species. Here, we characterize the population genetic structure and model dispersal dynamics of the barnacle B. improvisus, and describe how human-mediated spreading via shipping as well as natural larval dispersal may have contributed to observed genetic variation. We used both mitochondrial DNA (cytochrome c oxidase subunit I: COI) and nuclear microsatellites to characterize the genetic structure in 14 populations of B. improvisus on a global and regional scale (Baltic Sea). Genetic diversity was high in most populations, and many haplotypes were shared among populations on a global scale, indicating that long-distance dispersal (presumably through shipping and other anthropogenic activities) has played an important role in shaping the population genetic structure of this cosmopolitan species. We could not clearly confirm prior claims that B. improvisus originates from the western margins of the Atlantic coasts; although there were indications that Argentina could be part of a native region. In addition to dispersal via shipping, we show that natural larval dispersal may play an important role for further colonisation following initial introduction.