Genetic diversity and effective size of the Atlantic salmon Salmo salar L. inhabiting the River Eo (Spain) following a stock collapse

Historic angling records suggest the occurrence of a drastic decline in the River Eo (Spain) Atlantic salmon population size during the past two decades, as a result of overexploitation and habitat deterioration. In recent years, the population has been apparently recovering, and the present study is aimed to report information on the level of genetic diversity and the effective size of the current population as these may have immediate consequences for its conservation. Eighty-six salmon from two temporal groups (1998–1999 and 2004–2006), representing three generations, were genotyped using a panel of eight microsatellites. Inspite of the recent decline in census numbers and the detection of the signs of a population bottleneck, the population exhibits a high level of genetic diversity, similar to that from other populations, and almost unchanged during the period of study [average allelic richness ( A ) = 14·0 and 13·9, and average heterozygosity ( H _e) = 0·843 and 0·851 in 1998–1999 and 2004–2006, respectively]. The effective population size ( N _e) estimated by two different temporal methods showed a consistent value around 80 salmon, whereas the estimates from the linkage disequilibrium (LD) method provided a value around 165 individuals for either sample. The recent growing number of salmon, as indicated by fisheries records, the relatively large estimates of the ratio N _e/ N (with range 0·23–0·44 for the temporal estimates and 0·31–0·59 for the LD estimates) and the high levels of diversity found suggest that the population has not been greatly affected by the historical census declines and can be expected to recover in the future.     

Loss of genetic variation in hatchery-reared Indian major carp, Catla catla

The hypothesis that effective population sizes are low in hatchery-reared catla ( Catla catla ) from Bangladesh, possibly leading to inbreeding and loss of variation, was tested. The study was based on analysis of seven microsatellite loci in three samples of hatchery-reared catla and four samples representing wild populations. Pair-wise estimates of genetic differentiation between samples were low between wild samples (θ ranging from 0·012 to 0·034), but high between hatchery samples (θ ranging from 0·153 to 0·185), suggesting strong genetic drift in hatcheries. Genetic variation, both in terms of expected heterozygosity and allelic richness, was significantly lower in hatchery samples than in samples of wild catla. Application of a method for reconstructing families among offspring without parental genetic data showed that the hatchery samples consisted of very few half- and full-sib families, whereas the wild samples consisted of a high number of families, suggesting that most individuals were unrelated. Finally, estimation of the effective number of parents ( N _b) in the largest sample of hatchery fish confirmed that effective population size was low ( N _b= 14·9 for multiallelic loci and N _b= 10·6 if alleles were pooled into two composite alleles). The results show that low effective population sizes leading to loss of variation and possibly inbreeding depression should be a matter of serious concern in aquaculture production of catla.     

Temporal effective size estimates of a managed walleye Sander vitreus population and implications for genetic-based management

The goal of this research was to use the long-term fishery data set and DNA from archived scales of walleye Sander vitreus in Escanaba Lake, WI, U.S.A., to improve the understanding of the underlying mechanism(s) influencing genetic diversity in naturally recruiting populations. The introduced population of S. vitreus in Escanaba Lake has a low mean effective population size ( N _E) between 124·6 and 185·5 despite a mean census size ( N _C) of 4659 ( N _E/ N _C c. 0·04), suggesting an accelerated rate of genetic drift between 1952 and 2002. These values are smaller than the median N _E range of several studies suggesting typical N _E/ N _C ratios of 0·11–0·16 in a wide range of taxa. N _E increased steadily during the past two sampled decades (1992 and 2002) and was consistent with a lowering of the variance in S. vitreus reproductive success, possibly linked to a large, sustained exploitation (mean 28%) rate. Variance in reproductive success is one of the most important factors influencing N _E in species, like S. vitreus , which have a potential for large fecundities and large juvenile mortalities (type III survivorship). The N _B estimates across six sequential cohorts (age classes of S. vitreus , assayed from 1994 to 1999) was consistent with estimates of N _E reported for 1992–2002. These results, coupled with in-depth census and exploitation data, show that the genetic characteristics of Escanaba Lake S. vitreus have changed substantially and that management activities, such as supplemental stocking and harvest practices, have profoundly influenced the genetic dynamics of S. vitreus in this lake.     

Loss of genetic variation and effective population size of Kirikuchi charr: implications for the management of small, isolated salmonid populations

Small, isolated populations may face extinction due to a combination of inbreeding depression and other threats. Effective population size ( N _e) is one comprehensive measure that allows us to evaluate the genetic status of a population, and to make management decisions regarding genetic viability. We simulated loss of genetic variation and estimated N _e for two small, isolated populations of Kirikuchi charr Salvelinus leucomaenis japonicus , the endangered, southernmost local populations of the genus Salvelinus in the world, using VORTEX, an individual-based stochastic PVA model. Approximately half of the genetic variation was lost over 200 years regardless of census population size and demographic parameters, and N _e estimates were roughly 50 in each of the two populations, suggesting the possibility of inbreeding depression. The target population size of N _e>500, by securing long-term viability, is several times that of the present size of each of the populations studied, and no local habitats maintaining such a target number are considered to exist. The results strongly indicate a need for recovering natural connections and potential gene flow among local populations. However, the impending threat to these populations from non-native charr widely distributed throughout the drainage has prevented the recovery of the connections. Given the small N _e of the two populations, it would be necessary to retain gene flow artificially within or across local populations. This will be true of many other salmonid populations that have been isolated or fragmented recently.     

Usefulness of molecular markers for detecting population bottlenecks via monitoring genetic change

It is important to detect population bottlenecks in threatened and managed species because bottlenecks can increase the risk of population extinction. Early detection is critical and can be facilitated by statistically powerful monitoring programs for detecting bottleneck-induced genetic change. We used Monte Carlo computer simulations to evaluate the power of the following tests for detecting genetic changes caused by a severe reduction in a population's effective size ( N _e): a test for loss of heterozygosity, two tests for loss of alleles, two tests for change in the distribution of allele frequencies, and a test for small N _e based on variance in allele frequencies (the 'variance test'). The variance test was most powerful; it provided an 85% probability of detecting a bottleneck of size N _e = 10 when monitoring five microsatellite loci and sampling 30 individuals both before and one generation after the bottleneck. The variance test was almost 10-times more powerful than a commonly used test for loss of heterozygosity, and it allowed for detection of bottlenecks before 5% of a population's heterozygosity had been lost. The second most powerful tests were generally the tests for loss of alleles. However, these tests had reduced power for detecting genetic bottlenecks caused by skewed sex ratios. We provide guidelines for the number of loci and individuals needed to achieve high-power tests when monitoring via the variance test. We also illustrate how the variance test performs when monitoring loci that have widely different allele frequency distributions as observed in five wild populations of mountain sheep ( Ovis canadensis ).     

A preliminary investigation of allozyme genetic variation and population geographical structure in Aphanius fasciatus from Italian brackish-water habitats

A total of 150 individuals from Aphanius fasciatus from coastal brackish-water habitats was analysed by allozyme electrophoresis to collect data on its genetic variation. From 22 enzymes, 43 putative enzyme-coding loci were resolved, 12 of which were polymorphic at P_0·99 level. Only one of the 31 probability tests showed a significant departure from the Hardy-Weinberg equilibrium. Aphanius fasciatus showed low levels of genetic polymorphism, with expected heterozygosity values ranging from 0·027 (S.E.=0·013) to 0·064 ( s.e. =0·023). Nei's genetic distances between populations ranged from 0·002 to 0·042. Weir & Cockerham F -statistics showed high levels of genetic heterogeneity among populations (jackknifed θ=0·302, s.e. =0·045) and estimates of N _m were < 1, indicating restricted gene flow. Significant positive correlation between genetic distance and geographical distance matrices, detected by Mantel's test ( g =1·941; P 0·001), is consistent with the prediction that the species is genetically structured by isolation-by-distance.     

Stability of genetic structure and effective population size inferred from temporal changes of microsatellite DNA polymorphisms in the land snail Helix aspersa (Gastropoda: Helicidae)

Temporal evolution of genetic variability may have far-reaching consequences for a diverse array of evolutionary processes. Within the polders of the Bay of Mont-Saint-Michel (France), populations of the land snail Helix aspersa are characterized by a metapopulation structure with occasional extinction processes resulting from farming practices. A temporal survey of genetic structure in H . aspersa was carried out using variability at four microsatellite loci, in ten populations sampled two years apart. Levels of within-population genetic variation, as measured by allelic richness, H _e or F _is , did not change over time and similar levels of population differentiation were demonstrated for both sampling years. The extent of genetic differentiation between temporal samples of the same population established (i) a stable structure for six populations, and (ii) substantial genetic changes for four populations. Using classical F -statistics and a maximum likelihood method, estimates of the effective population size ( N _e) illustrated a mixture of stable populations with high N _e, and unstable populations characterized by very small N _e estimates (of 5–11 individuals). Owing to human disturbances, intermittent gene flow and genetic drift are likely to be the predominant evolutionary processes shaping the observed genetic structure. However, the practice of multiple matings and sperm storage is likely to provide a reservoir of variability, minimizing the eroding genetic effects of population size reduction and increasing the effective population size.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 89–102.     

Genetic estimates of contemporary effective population size: what can they tell us about the importance of genetic stochasticity for wild population persistence?

Genetic stochasticity due to small population size contributes to population extinction, especially when population fragmentation disrupts gene flow. Estimates of effective population size ( N _e) can therefore be informative about population persistence, but there is a need for an assessment of their consistency and informative relevance. Here we review the body of empirical estimates of N _e for wild populations obtained with the temporal genetic method and published since Frankham's (1995 ) review. Theoretical considerations have identified important sources of bias for this analytical approach, and we use empirical data to investigate the extent of these biases. We find that particularly model selection and sampling require more attention in future studies.
We report a median unbiased N _e estimate of 260 (among 83 studies) and find that this median estimate tends to be smaller for populations of conservation concern, which may therefore be more sensitive to genetic stochasticity. Furthermore, we report a median N _e/ N ratio of 0.14, and find that this ratio may actually be higher for small populations, suggesting changes in biological interactions at low population abundances. We confirm the role of gene flow in countering genetic stochasticity by finding that N _e correlates strongest with neutral genetic metrics when populations can be considered isolated. This underlines the importance of gene flow for the estimation of N _e, and of population connectivity for conservation in general. Reductions in contemporary gene flow due to ongoing habitat fragmentation will likely increase the prevalence of genetic stochasticity, which should therefore remain a focal point in the conservation of biodiversity.     

Large variations in the ratio of effective breeding and census population sizes between two species of pond-breeding anurans

The viability of wild populations is frequently assessed by monitoring adult census sizes ( N _c). This approach is particularly useful for pond-breeding amphibians, because assemblages during the breeding season are relatively easy to detect and count. However, it is the genetic effective population size ( N _e) or surrogates such as effective breeding population size ( N _b) that are of primary importance for long-term viability. Although N _c estimates of one anuran amphibian ( Bufo bufo ) in Britain were much larger than those of another ( Rana temporaria ) at the same sites, the ratios of N _b to N _c were much smaller in B. bufo than in R. temporaria. These differences were sufficiently great as to reverse the effective size order at one site, such that N _b for R. temporaria was larger than that for B. bufo. Differences in adult sex ratios at breeding sites probably contributed to lower N _b values in B. bufo populations compared with those of R. temporaria . The relationship of N _b to N _c can therefore vary dramatically even between similar species, to the extent that just monitoring N _c can give misleading impressions of relative effective breeding sizes and thus of population viability. It will be increasingly important to estimate N _e or N _b in wildlife populations for assessment of conservation priorities.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 89 , 365–372.     

Effective number of breeding adults in Bufo bufo estimated from age-specific variation at minisatellite loci

Estimates of the effective number of breeding adults were derived for three semi-isolated populations of the common toad Bufo bufo based on temporal (i.e. adult-progeny) variance in allele frequency for three highly polymorphic minisatellite loci. Estimates of spatial variance in allele frequency among populations and of age-specific measures of genetic variability are also described. Each population was characterized by a low effective adult breeding number ( N _b) based on a large age-specific variance in mini-satellite allele frequency. Estimates of N _b (range 21–46 for population means across three loci) were ≊ 55–230-fold lower than estimates of total adult census size. The implications of low effective breeding numbers for long-term maintenance of genetic variability and population viability are discussed relative to the species' reproductive ecology, current land-use practices, and present and historical habitat modification and loss. The utility of indirect measures of population parameters such as N _b and N _e based on time-series data of minisatellite allele frequencies is discussed relative to similar measures estimated from commonly used genetic markers such as protein allozymes.     

Allozymic variation and relationships of the endangered cyprinodontid genus Valencia and its implications for conservation

Allozyme analysis was used to determine patterns of genetic variation and relationships within the genus Valencia , a group with two allopatric species, V. letourneuxi inhabiting the Balkan Peninsula, and V. hispanica which is endemic to the Iberian Peninsula. Products of 25 gene loci were analysed with AK *, FH *, LDH -1* and LDH -2* monomorphic when Fundulus heteroclitus macrolepidotus was not taken into account. Nine loci were diagnostic for both species of the genus Valencia . Levels of genetic variation of V. hispanica ( P =0·416, H ^e=0·118) and V. letourneuxi ( P =0·160, H _e=0·040) were higher than values reported previously. High population subdivision ( F _st=0·321) among V. hispanica populations indicates a clear interruption of genetic exchange among populations. High genetic variation and differentiation exhibited by populations of V. hispanica suggest that a recovery program should be carried out with natural stocks from the same localities.     

Genetic population structure of the catadromous Perciform: Macquaria novemaculeata (Percichthyidae)

Genetic population structure in the catadromous Australian bass Macquaria novemaculeata was investigated using samples from four locations spanning 600 km along the eastern Australian coastline. Both allozymes and mtDNA control region sequences were examined. Population subdivision estimates based on allozymes revealed low levels of population structuring ( G_st =0·043, P <0·05). However, mtDNA indicated moderate levels of geographic population structure ( G_st =0·146, P <0·0l). Phylogenetic analysis of mtDNA control region sequences (mean sequence divergence 1·9%) indicated little phylogeographic structuring. Results suggested that genotypic variation within each river population, while being affected primarily by genetic drift, was also prevented from more significant divergence by homogenizing levels of gene flow—synonymous with a one-dimensional stepping-stone model of population structure. The catadromous life history of Macquaria novemaculeata was considered to be influential on the pattern of population structure displayed. Results were compared to the few population genetic studies involving catadromous fishes, indicating that catadromy alone is unlikely to be a good predictor of population structure. A more comprehensive suite of biological characteristics than simple life-history traits must be considered fully to allow reliable predictive models of population structure to be formulated.     

Differential population history in the migratory catfishes Brachyplatystoma flavicans and Pseudoplatystoma fasciatum(Pimelodidae) from the Bolivian Amazon assessed with nuclear and mitochondrial DNA markers

The catfishes Brachyplatystoma flavicans ( n  = 49) and Pseudoplatystoma fasciatum ( n  = 69) showed comparable low allozyme diversities ( H _e = 0·012 and 0·009–0·028, respectively), but contrasting PCR‐RFLP restriction site mitochondrial DNA diversities (three haplotypes: π = 0·034–0·092 and five haplotypes: π = 0·001–0·023, respectively) in the Rio Ichilo and Beni (Bolivia). Genetic homogeneity between samples was high for B. flavicans and lower for P. fasciatum . Based on mitochondrial diversity, both species probably experienced a historic population reduction but at different time scales.     

Population genetic structure of the butterfly Melitaea didyma (Nymphalidae) along a northern distribution range border

The population genetic structure of the butterfly Melitaea didyma was studied along the northern distribution range border in Central Germany by means of allozyme electrophoresis. Individuals were sampled from a total of 21 habitat patches from four regions, and two provinces. Sampling was designed to estimate local vs. regional differentiation. High levels of variability were found, H _e= 0.14–0.21. The mean expected sample heterozygosity from one region, Mosel, was significantly lower than from the Hammelburg region, H _e= 0.17 and 0.19, respectively. Two hierarchical levels of genetic differentiation were found. Within regions individuals sampled from different patches behaved as belonging to one population with high levels of gene flow (Hammelburg F _ST= 0.015, Mosel F _ST= 0.044), though local isolation barriers did create a substructuring of these populations. The inbreeding coefficients, F _IS, were constant over all sample levels, suggesting a similar distribution of habitat patches within regions. Between regions gene flow was limited. An isolation by distance analysis indicated that the hierarchical structure, at the provincial level, may be breaking down due to isolation of regional populations. A more general observation was that the sampling design may greatly have influenced the estimation of genetic differentiation. Depending on which samples were included, overall F _ST estimates ranged from 0.059–0.090.     

Genetics of the peripheral populations of the alpine bullhead, Cottus poecilopus (Scorpaeniformes, Cottidae) in Poland

The alpine bullhead ( Cottus poecilopus ) in Poland is highly polymorphic ( P for examined populations in the range of 0.0–0.45, mean P  = 0.11, P for species 0.64) in seven out of 11 loci studied ( Pgd-1, Mdh-1, Ck-1, Gpd-1, Gpd-2, Gpi-1, Gpi-2 ). However, they reveal very low heterozygosity (mean H _o = 0.007, mean H _e = 0.019), which may result from population subdivisions, small genetically effective size ( N _e) of some populations and frequent inbreeding, the latter being probably mediated through polygyny. The index of genetic differentiation ( G _ST) for the Polish metapopulation equals 0.48, which is relatively high. By far, the strongest are genetic differences between the Sudeten-Carpathian group and the Hańcza population. Gene flow between the Sudeten-Carpathian group and the Hańcza population is very limited and their divergence may have started as early as the last (Eemian) interglacial period.     

Sibship reconstruction demonstrates the extremely low effective population size of striped bass Morone saxatilis in the Santee–Cooper system, South Carolina, USA

For organisms with great fecundity and high mortality in early life stages, such as shellfish or fishes, the need to match reproductive activity with environmental conditions conducive to spawning, fertilization, larval development and recruitment may result in extreme variance in reproductive success among individuals. The main objective of this study was to investigate evidence of large variance in the reproductive success of the striped bass Morone saxatilis in the Santee–Cooper system, South Carolina, USA. Seven microsatellite loci were analysed in 603 recruits representing three yearly cohorts from 1992 to 1994, and a group analysis was performed to identify full-sib families. Large variance in reproductive success was detected, with a few large, full-sib families contributing disproportionately to each of the cohorts. The severity of sweepstakes reproductive success varied among cohorts depending on environmentally imposed mortality. Estimations of the effective number of breeders in these long-lived fish ranged from 24 in 1992 to 44 in 1994. Furthermore, the estimated genetic effective population size ( N _e = 93) is approximately four orders of magnitude lower than estimates of adult census size ( N  =   362 000). Furthermore, the presence of large full-sib families indicates that striped bass engage in pair mating in the wild. Heterogeneity in genetic composition was also observed among cohorts, suggesting that genetically different adults contribute to different cohorts and that chance rather than fitness variation determines reproductive success.     

Population structure of lake‐type and river‐type sockeye salmon in transboundary rivers of northern British Columbia

The population structure of 'lake‐type' and 'river‐type' sockeye salmon Oncorhynchus nerka , primarily in transboundary rivers in northern British Columbia, was examined with a survey of microsatellite variation. Variation at 14 microsatellite loci was surveyed from c . 3000 lake‐type and 3200 river‐type sockeye salmon from 47 populations in six river drainages in British Columbia. The mean F _ST for the 14 microsatellite loci and 47 populations was 0·068, and 0·034 over all river‐type populations. River‐type sockeye salmon were more genetically diverse than lake‐type sockeye salmon, with expected heterozygosity of river‐type sockeye salmon 0·72 and with an average 12·7 alleles observed per locus, whereas expected heterozygosity of lake‐type sockeye salmon was 0·65 with and average 10·5 alleles observed per locus. River drainage of origin was a significant unit of population structure. There was clear evidence of genetic differentiation among river‐type populations of sockeye salmon from different drainages over a broad geographic range in British Columbia.     

A comparison of genetic diversity levels in marine, freshwater, and anadromous fishes

Electrophoretic data were analysed from 49 species of freshwater fish, 57 species of marine fish, and seven anadromous species. For each species, at least 15 individuals had been assayed for at least 15 loci in two or more subpopulations. The results showed that while average total heterozygosity ( _T ) was approximately equal in freshwater and marine species (0·062 and 0·064 respectively), subpopulation heterozygosity ( _s ) was significantly less in the former group (0·046 and 0·059 respectively). Consequently the average degree of genetic subpopulation differentiation ( _ST ) was significantly greater for freshwater species (0·222 v. 0·062). On average, it is likely that marine subpopulations exchange between 10 and 100 times more migrants per generation than freshwater subpopulations, presumably because of the relative absence of barriers to dispersal in the marine environment. The reduced values of H_s in freshwater species are likely to reflect reduced effective subpopulation sizes relative to marine species. The few andromous species that have been analysed show intermediate levels of G_ST .     

Genetic and morphological variation in a common European cyprinid, Leuciscus cephalus within and across Central European drainages

Genetic variability and differentiation of the common European cyprinid Leusciscus cephalus was investigated within and across the drainages of the rivers Rhine, Elbe and Danube using starch gel electrophoresis and morphological characters. Eleven of 24 presumptive enzyme loci showed some allozymic variation. Compared to other freshwater fish, genetic variability within all sites was high ( H _e=0·074–0·113). This was explained by the wide geographical and ecological range of the chub. In contrast, genetic and morphological differentiation between sites and drainages was low (F_ST=0·15), which may be due to the high dispersal ability of the chub and translocations of spawn by waterfowl. The phylogeographic structure of Leuciscus cephalus in north-eastern Bavaria showed low congruence with Cottus gobio populations sampled in the same area.     

The influence of thermal parameters on the acclimation responses of pinfish Lagodon rhomboides exposed to static and decreasing low temperatures

Pinfish Lagodon rhomboides acclimation rates were determined by modelling changes in critical thermal minimum ( T _{crit min}, ° C) estimates at set intervals following a temperature decrease of 3–4° C. The results showed that pinfish gained a total of 3·7° C of cold tolerance over a range of acclimation temperatures ( T _acc, ° C) from (23–12° C), that cold tolerance increased with exposure time to the reduced temperature at all T _acc, but that the rate of cold tolerance accruement (mean 0·14° C day⁻¹) was independent of T _acc. A highly significant ( P < 0·001) multivariate predictive model was generated that described the acclimation rates and thermal tolerance of pinfish exposed to reduction in water temperature: log₁₀ T _{crit min}= 0·41597 − 0·01704 T _acc+ 0·04320 T _plunge− 0·08376[log₁₀ ( t + 1)], where T _plunge is plunge temperature (° C) and t is the time (days). A comparison of the present data, with acclimation rate data for other species, suggests that factors such as latitude or geographic range may play a more important role than ambient temperature in determining cold acclimation rates in fishes.