Temporal changes in allele frequencies but stable genetic diversity over the past 40 years in the Irish Sea population of thornback ray, Raja clavata

Population structure and demographic history for year cohort dynamics of landlocked ayu Plecoglossus altivelis altivelis in dam reservoir of Japan

Ayu Plecoglossus altivelis altivelis is a valuable osmeroid species for inland fishery in Japan. It is classified into two ecological forms of amphidromous migrating between rivers and sea and landlocked migrating between rivers and lakes or dam reservoirs. The number of dams and their reservoirs has remarkably increased in the twenty-first century under climate change, because of their respective roles in hydropower generation with negligible carbon emissions and in flood control. Dam reservoirs therefore become increasingly important as inland nursery grounds of ayu. In this study, we investigated the reproduction status of landlocked ayu migrating in the Haidzuka Dam reservoir and the Tabusa River in western Japan by molecular phylogenetic analysis based on population structure and demographic history for year cohort dynamics. A total of 849 individuals were collected monthly from October 2018 to September 2021 according to an annual life cycle of ayu. Nucleotide sequences of the partial mitochondrial DNA control region yielded 31 haplotypes, consisting of 4 shared haplotypes among the 2019, 2020 and 2021 cohorts and 27 unique haplotypes. The overall haplotype diversity and nucleotide diversity were calculated to be relatively low at 0.3503 ± 0.0206 and 0.0077 ± 0.0045, respectively, suggesting a founder event by dominant haplotypes. Star-shaped radiational haplotypes from dominant shared haplotypes on the median-joining network likely support a founder event. Although pairwise ФST values were determined to be very low among the year cohorts, only the 2019 cohort was found to have a significant difference from the 2020 and 2021 cohorts, for both of which Tajima's D values were also statistically significant. For the overall population, multimodal mismatch distribution and negative Tajima's D and Fu's Fs values in the neutrality test suggested population expansion or population subdivision. The native riverine population in the Tabusa River suffered habitat fragmentation and population bottleneck from dam construction, and therefore severe founder effect remained behind the artificially landlocked population with a low level of genetic diversity in the Haidzuka Dam reservoir.

Retrospective genomics highlights changes in genetic composition of tiger sharks (Galeocerdo cuvier) and potential loss of a south-eastern Australia population

Over the last century, many shark populations have declined, primarily due to overexploitation in commercial, artisanal and recreational fisheries. In addition, in some locations the use of shark control programs also has had an impact on shark numbers. Still, there is a general perception that populations of large ocean predators cover wide areas and therefore their diversity is less susceptible to local anthropogenic disturbance. Here we report on temporal genomic analyses of tiger shark (Galeocerdo cuvier) DNA samples that were collected from eastern Australia over the past century. Using Single Nucleotide Polymorphism (SNP) loci, we documented a significant change in genetic composition of tiger sharks born between ~1939 and 2015. The change was most likely due to a shift over time in the relative contribution of two well-differentiated, but hitherto cryptic populations. Our data strongly indicate a dramatic shift in the relative contribution of these two populations to the overall tiger shark abundance on the east coast of Australia, possibly associated with differences in direct or indirect exploitation rates.

Time matters: genetic composition and evaluation of effective population size in temperate coastal fish species

Background

Extensive knowledge on the genetic characterization of marine organisms has been assembled, mainly concerning the spatial distribution and structuring of populations. Temporal monitoring assesses not only the stability in genetic composition but also its trajectory over time, providing critical information for the accurate forecast of changes in genetic diversity of marine populations, particularly important for both fisheries and endangered species management. We assessed fluctuations in genetic composition among different sampling periods in the western Portuguese shore in three fish species.

Methods

White seabream Diplodus sargus, sand smelt Atherina presbyter and shanny Lipophrys pholis were chosen, because of their genetic patterns in distinct ecological environments, insight into historical and contemporary factors influencing population effective size (N_e), and degree of commercial exploitation. Samples were obtained near Lisbon between 2003 and 2014 and screened for genetic variation with mitochondrial and nuclear markers. Analyses included genealogies, genetic diversities, temporal structures and contemporary N_e.

Results

For mtDNA no temporal structure was detected, while for nDNA significant differences were recorded between some sampling periods for the shanny and the sand smelt. Haplotype networks revealed deep genealogies, with various levels of diversification. The shanny revealed a smaller N_e/generation when compared to the other species, which, in turn, revealed no evidence of genetic drift for most study periods. These results highlight the fact that temporal variations in genetic pool composition should be considered when evaluating the population structure of fish species with long distance dispersal, which are more vulnerable to recruitment fluctuations.

Estimating effective population size using RADseq: Effects of SNP selection and sample size

Effective population size (Ne ) is a key parameter of population genetics. However, N e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site-associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To evaluate the potential impact of SNP selection criteria on Ne estimates (Linkage Disequilibrium method) we used RADseq data for a nonmodel species, the thornback ray. In this data set, the inbreeding coefficient F IS was positively correlated with the amount of missing data, implying data were missing nonrandomly. The precision of Ne estimates decreased with the number of SNPs. Mean Ne estimates (averaged across 50 random data sets with2000 SNPs) ranged between 237 and 1784. Increasing the percentage of missing data from 25% to 50% increased Ne estimates between 82% and 120%, while increasing the minor allele frequency (MAF) threshold from 0.01 to 0.1 decreased estimates between 71% and 75%. Considering these effects is important when interpreting RADseq data-derived estimates of effective population size in empirical studies.

Climate change in the oceans: evolutionary versus phenotypically plastic responses of marine animals and plants

I summarize marine studies on plastic versus adaptive responses to global change. Due to the lack of time series, this review focuses largely on the potential for adaptive evolution in marine animals and plants. The approaches were mainly synchronic comparisons of phenotypically divergent populations, substituting spatial contrasts in temperature or CO2 environments for temporal changes, or in assessments of adaptive genetic diversity within populations for traits important under global change. The available literature is biased towards gastropods, crustaceans, cnidarians and macroalgae. Focal traits were mostly environmental tolerances, which correspond to phenotypic buffering, a plasticity type that maintains a functional phenotype despite external disturbance. Almost all studies address coastal species that are already today exposed to fluctuations in temperature, pH and oxygen levels. Recommendations for future research include (i) initiation and analyses of observational and experimental temporal studies encompassing diverse phenotypic traits (including diapausing cues, dispersal traits, reproductive timing, morphology) (ii) quantification of nongenetic trans-generational effects along with components of additive genetic variance (iii) adaptive changes in microbe-host associations under the holobiont model in response to global change (iv) evolution of plasticity patterns under increasingly fluctuating environments and extreme conditions and (v) joint consideration of demography and evolutionary adaptation in evolutionary rescue approaches.

A review of the application of molecular genetics for fisheries management and conservation of sharks and rays

Since the first investigation 25 years ago, the application of genetic tools to address ecological and evolutionary questions in elasmobranch studies has greatly expanded. Major developments in genetic theory as well as in the availability, cost effectiveness and resolution of genetic markers were instrumental for particularly rapid progress over the last 10 years. Genetic studies of elasmobranchs are of direct importance and have application to fisheries management and conservation issues such as the definition of management units and identification of species from fins. In the future, increased application of the most recent and emerging technologies will enable accelerated genetic data production and the development of new markers at reduced costs, paving the way for a paradigm shift from gene to genome-scale research, and more focus on adaptive rather than just neutral variation. Current literature is reviewed in six fields of elasmobranch molecular genetics relevant to fisheries and conservation management (species identification, phylogeography, philopatry, genetic effective population size, molecular evolutionary rate and emerging methods). Where possible, examples from the Indo-Pacific region, which has been underrepresented in previous reviews, are emphasized within a global perspective.

Genetic monitoring and complex population dynamics: insights from a 12-year study of the Rio Grande silvery minnow

The endangered Rio Grande silvery minnow persists as a remnant population in a highly fragmented and regulated arid-land river system. The species is subject to dramatic annual fluctuations in density. Since 2003, the wild population has been supplemented by hatchery-reared fish. We report on a 12-year (1999–2010) monitoring study of genetic diversity and effective population size (N_e) of wild and hatchery stocks. Our goals were to evaluate how genetic metrics responded to changes in wild fish density and whether they corresponded to the number and levels of diversity of hatchery-reared repatriates. Genetic diversity and all measures of N_e in the wild population did not correlate with wild fish density until hatchery supplementation began in earnest. Estimates of variance and inbreeding effective size were not correlated. Our results suggest source–sink dynamics where captive stocks form a genetically diverse source and the wild population behaves as a sink. Nevertheless, overall genetic diversity of silvery minnow has been maintained over the last decade, and we attribute this to a well-designed and executed propagation management plan. When multiple factors like environmental fluctuation and hatchery supplementation act simultaneously on a population, interpretation of genetic monitoring data may be equally complex and require considerable ecological data.

Temporal genetic stability and high effective population size despite fisheries-induced life-history trait evolution in the North Sea sole

Heavy fishing and other anthropogenic influences can have profound impact on a species' resilience to harvesting. Besides the decrease in the census and effective population size, strong declines in mature adults and recruiting individuals may lead to almost irreversible genetic changes in life-history traits. Here, we investigated the evolution of genetic diversity and effective population size in the heavily exploited sole (Solea solea), through the analysis of historical DNA from a collection of 1379 sole otoliths dating back from 1957. Despite documented shifts in life-history traits, neutral genetic diversity inferred from 11 microsatellite markers showed a remarkable stability over a period of 50 years of heavy fishing. Using simulations and corrections for fisheries induced demographic variation, both single-sample estimates and temporal estimates of effective population size (N(e) ) were always higher than 1000, suggesting that despite the severe census size decrease over a 50-year period of harvesting, genetic drift is probably not strong enough to significantly decrease the neutral diversity of this species in the North Sea. However, the inferred ratio of effective population size to the census size (N(e) /N(c) ) appears very small (10(-5) ), suggesting that overall only a low proportion of adults contribute to the next generation. The high N(e) level together with the low N(e) /N(c) ratio is probably caused by a combination of an equalized reproductive output of younger cohorts, a decrease in generation time and a large variance in reproductive success typical for marine species. Because strong evolutionary changes in age and size at first maturation have been observed for sole, changes in adaptive genetic variation should be further monitored to detect the evolutionary consequences of human-induced selection.

Demographic processes underlying subtle patterns of population structure in the scalloped hammerhead shark, Sphyrna lewini

Genetic diversity (θ), effective population size (N(e)), and contemporary levels of gene flow are important parameters to estimate for species of conservation concern, such as the globally endangered scalloped hammerhead shark, Sphyrna lewini. Therefore, we have reconstructed the demographic history of S. lewini across its Eastern Pacific (EP) range by applying classical and coalescent population genetic methods to a combination of 15 microsatellite loci and mtDNA control region sequences. In addition to significant population genetic structure and isolation-by-distance among seven coastal sites between central Mexico and Ecuador, the analyses revealed that all populations have experienced a bottleneck and that all current values of θ are at least an order of magnitude smaller than ancestral θ, indicating large decreases in N(e) (θ = 4N(e)μ), where μ is the mutation rate. Application of the isolation-with-migration (IM) model showed modest but significant genetic connectivity between most sampled sites (point estimates of Nm = 0.1-16.7), with divergence times (t) among all populations significantly greater than zero. Using a conservative (i.e., slow) fossil-based taxon-specific phylogenetic calibration for mtDNA mutation rates, posterior probability distributions (PPDs) for the onset of the decline in N(e) predate modern fishing in this region. The cause of decline over the last several thousand years is unknown but is highly atypical as a post-glacial demographic history. Regardless of the cause, our data and analyses suggest that S. lewini was far more abundant throughout the EP in the past than at present.

Repeated sampling detects gene flow in a flightless ground beetle in a fragmented landscape

Secondary clines level down in the course of time if the gene flow is not interrupted. Temporally repeated sampling of populations in a cline allows the investigation not only of its occurrence but also of the estimation of the amount of ongoing gene flow. We reinvestigated an allozyme gradient in Carabus auronitens populations in the Westphalian Lowlands (northwestern Germany) 15 to 20 years after it had originally been recorded. A total of 977 individuals of this flightless woodland species from 29 sample sites were genotyped at the diallelic Est-1 locus in 2005-2006 and compared to former findings, collected in 1985-1994 from the same populations. Both data sets showed clinal variation. Pairwise differences between the samples of both data sets indicated significant decrease in the steepness of the cline during the past 15 to 20 years. The estimated average gene flow per generation is 0.6% of each beetle population. Ongoing gene flow in the flightless ground beetle C. auronitens led to a less pronounced cline despite a stable degree of fragmentation (and connectivity) of the landscape. Migration and gene flow were obviously enabled by the numerous hedgerows. The corridors are seen to be a prerequisite for migration between populations and for possible future range shifts of forest insect species.

Temporally stable genetic structure of heavily exploited Atlantic herring (Clupea harengus) in Swedish waters

Information on the temporal stability of genetic structures is important to permit detection of changes that can constitute threats to biological resources. Large-scale harvesting operations are known to potentially alter the composition and reduce the variability of populations, and Atlantic herring (Clupea harengus) has a long history of heavy exploitation. In the Baltic Sea and Skagerrak waters, the census population sizes have declined by 35-50% over the last three decades. We compared the genetic structure of Atlantic herring in these waters sampled at least two different times between 1979 and 2003 by assaying 11 allozyme and nine microsatellite loci. We cannot detect any changes in the amount of genetic variation or spatial structure, and differentiation is weak with overall F(ST)=0.003 among localities for the older samples and F(ST)=0.002 for the newer ones. There are indications of temporal allele frequency changes, particularly in one of five sampling localities that is reflected in a relatively small local N(e) estimate of c. 400. The previously identified influence of selection at the microsatellite locus Cpa112 remains stable over the 24-year period studied here. Despite little genetic differentiation, migration among localities appears small enough to permit demographic independence between populations.

Remarkably low mtDNA control region diversity in an abundant demersal fish

Cape hake, Merluccius paradoxus, is a valuable commercially exploited demersal species. Using the 5' mtDNA control region we show that 96% of 1013 fishes sampled over a three-year period share one of two dominant haplotypes; 19 haplotypes were recovered in total, suggesting a genetically homogenous population of fish. Accordingly, haplotype and nucleotide diversities are low (h = 0.53, pi = 0.0014); an asymptotic haplotype accumulation curve suggests that few additional haplotypes exist. Comparing h and pi with other fish species shows that M. paradoxus and other southern African fish species have remarkably low genetic diversity values compared with other global marine fishes. Despite low genetic variability, frequency differences among M. paradoxus haplotypes suggest weakly structured populations between Namibia and South Africa. However, given the remarkably homogeneous mtDNA population genetic structure between fishes sampled along 1800 km, it is clear that faster evolving markers such a microsatellites are also needed before inferences can be made regarding stock identification and management of this species.

Unexpected panmixia in a long-lived, deep-sea fish with well-defined spawning habitat and relatively low fecundity

The marine environment presents particular challenges for our understanding of the factors that determine gene flow and consequent population structure. For marine fish, various aspects of life history have been considered important in an environment with few physical barriers, but dominated by current patterns, often varying with depth. These factors include the abundance and longevity of larval stages, typically more susceptible to movement along current paths. It also includes adult body size, fecundity and longevity with 'r-selected' species typically thought capable of greater gene flow and consequent panmixia. Here we investigate the population genetics of the orange roughy (Hoplostethus atlanticus), a clearly 'K-selected' species with habitat dependence on sea mounts for spawning, relatively large body size, a brief larval stage and relatively low fecundity. We used 14 polymorphic microsatellite loci to test the hypothesis that these characteristics will result in philopatry and genetic structure in the Atlantic Ocean. We discuss possible evolutionary mechanisms that could explain the results, which show the opposite pattern, with effective panmixia across thousands of kilometres in the North Atlantic.