Can genetic estimators provide robust estimates of the effective number of breeders in small populations?

A critical assessment of estimating census population size from genetic population size (or vice versa) in three fishes

Technological and methodological advances have facilitated the use of genetic data to infer census population size (Nc) in natural populations, particularly where traditional mark-and-recapture is challenging. The effective number of breeders (Nb) describes how many adults effectively contribute to a cohort and is often correlated with Nc. Predicting Nc from Nb or vice versa in species with overlapping generations has important implications for conservation by permitting (i) estimation of the more difficult to quantify variable and (ii) inferences of Nb/Nc relationships in related species lacking data. We quantitatively synthesized Nb/Nc relationships in three salmonid fishes where sufficient data have recently accumulated. Mixed-effects models were analysed in which each variable was included as a dependent variable or predictor term (Nb from Nc and vice versa). Species-dependent Nb/Nc slope estimates were significantly positive in two of three species. Variation in species slopes was likely due to varying life histories and reinforce caution when inferring Nb/Nc from taxonomically related species. Models provided maximum probable estimates for Nb and Nc for two species. However, study, population and year effects explained substantial amounts of variation (39%-57%). Consequently, prediction intervals were wide and included or were close to zero for all population sizes and species; model predictive utility was limited. Cost-benefit trade-offs when estimating Nb and/or Nc were also discussed using a real-world system example. Our findings based on salmonids suggest that no short cuts currently exist when estimating population size and researchers should focus on quantifying the variable of interest or be aware of caveats when inferring the desired variable because of cost or logistics. We caution that the salmonid species examined share life-history traits that may obscure relationships between Nb and Nc. Sufficient data on other taxa were unavailable; additional research examining Nb/Nc relationships in species with potentially relevant life-history trait differences (e.g., differing survival curves) is needed.

Temporal trends in genetic data and effective population size support efficacy of management practices in critically endangered dusky gopher frogs (Lithobates sevosus)

Monitoring temporal changes in population genetic diversity and effective population size can provide vital information on future viability. The dusky gopher frog, Lithobates sevosus, is a critically endangered species found only in coastal Mississippi, with low genetic variability as a consequence of isolation and population size reduction. Conservation management practices have been implemented, but their efficacy has not been addressed. We genotyped individuals collected 1997-2014 to determine temporal trends in population genetic variation, structure, and effective size. Observed and expected heterozygosity and allelic richness revealed temporally stable, but low, levels of genetic variation. Positive levels of inbreeding were found in each year. There was weak genetic structure among years, which can be attributed to increased effects of genetic drift and inbreeding in small populations. L. sevosus exhibited an increase in effective population size, and currently has an estimated effective size of 33.0-58.6 individuals, which is approximately half the census size. This large ratio could possibly be explained by genetic compensation. We found that management practices have been effective at maintaining and improving effective size and genetic diversity, but that additional strategies need to be implemented to enhance viability of the species.

Repeated Habitat Disturbances by Fire Decrease Local Effective Population Size

Effective population size is a fundamental parameter in population genetics, and factors that alter effective population size will shape the genetic characteristics of populations. Habitat disturbance may have a large effect on genetic characteristics of populations by influencing immigration and gene flow, particularly in fragmented habitats. We used the Florida Sand Skink (Plestiodon reynoldsi) to investigate the effect of fire-based habitat disturbances on the effective population size in the highly threatened, severely fragmented, and fire dependent Florida scrub habitat. We screened 7 microsatellite loci in 604 individuals collected from 12 locations at Archbold Biological Station. Archbold Biological Station has an active fire management plan and detailed records of fires dating to 1967. Our objective was to determine how the timing, number, and intervals between fires affect effective population size, focusing on multiple fires in the same location. Effective population size was higher in areas that had not been burned for more than 10 years and decreased with number of fires and shorter time between fires. A similar pattern was observed in abundance: increasing abundance with time-since-fire and decreasing abundance with number of fires. The ratio of effective population size to census size was higher at sites with more recent fires and tended to decrease with time-since-last-fire. These results suggest that habitat disturbances, such as fire, may have a large effect in the genetic characteristics of local populations and that Florida Sand Skinks are well adapted to the natural fire dynamics required to maintain Florida scrub.

Use of sibling relationship reconstruction to complement traditional monitoring in fisheries management and conservation of brown trout

Declining trends in the abundance of many fish urgently call for more efficient and informative monitoring methods that would provide necessary demographic data for the evaluation of existing conservation, restoration, and management actions. We investigated how genetic sibship reconstruction from young-of-the-year brown trout (Salmo trutta L.) juveniles provides valuable, complementary demographic information that allowed us to disentangle the effects of habitat quality and number of breeders on juvenile density. We studied restored (n = 15) and control (n = 15) spawning and nursery habitats in 16 brown trout rivers and streams over 2 consecutive years to evaluate the effectiveness of habitat restoration activities. Similar juvenile densities both in restored and control spawning and nursery grounds were observed. Similarly, no differences in the effective number of breeders, Nb(SA) , were detected between habitats, indicating that brown trout readily used recently restored spawning grounds. Only a weak relationship between the Nb(SA) and juvenile density was observed, suggesting that multiple factors affect juvenile abundance. In some areas, very low estimates of Nb(SA) were found at sites with high juvenile density, indicating that a small number of breeders can produce a high number of progeny in favorable conditions. In other sites, high Nb(SA) estimates were associated with low juvenile density, suggesting low habitat quality or lack of suitable spawning substrate in relation to available breeders. Based on these results, we recommend the incorporation of genetic sibship reconstruction to ongoing and future fish evaluation and monitoring programs to gain novel insights into local demographic and evolutionary processes relevant for fisheries management, habitat restoration, and conservation.

Temporal stability of genetic variability and differentiation in the three-spined stickleback (Gasterosteus aculeatus)

Temporal variation in allele frequencies, whether caused by deterministic or stochastic forces, can inform us about interesting demographic and evolutionary phenomena occurring in wild populations. In spite of the continued surge of interest in the genetics of three-spined stickleback (Gasterosteus aculeatus) populations, little attention has been paid towards the temporal stability of allele frequency distributions, and whether there are consistent differences in effective size (Ne) of local populations. We investigated temporal stability of genetic variability and differentiation in 15 microsatellite loci within and among eight collection sites of varying habitat type, surveyed twice over a six-year time period. In addition, Nes were estimated with the expectation that they would be lowest in isolated ponds, intermediate in larger lakes and largest in open marine sites. In spite of the marked differences in genetic variability and differentiation among the study sites, the temporal differences in allele frequencies, as well as measures of genetic diversity and differentiation, were negligible. Accordingly, the Ne estimates were temporally stable, but tended to be lower in ponds than in lake or marine habitats. Hence, we conclude that allele frequencies in putatively neutral markers in three-spined sticklebacks seem to be temporally stable - at least over periods of few generations - across a wide range of habitat types differing markedly in levels of genetic variability, effective population size and gene flow.

Genetic Diversity of White Sharks, Carcharodon carcharias, in the Northwest Atlantic and Southern Africa

The white shark, Carcharodon carcharias, is both one of the largest apex predators in the world and among the most heavily protected marine fish. Population genetic diversity is in part shaped by recent demographic history and can thus provide information complementary to more traditional population assessments, which are difficult to obtain for white sharks and have at times been controversial. Here, we use the mitochondrial control region and 14 nuclear-encoded microsatellite loci to assess white shark genetic diversity in 2 regions: the Northwest Atlantic (NWA, N = 35) and southern Africa (SA, N = 131). We find that these 2 regions harbor genetically distinct white shark populations (Φ ST = 0.10, P < 0.00001; microsatellite F ST = 0.1057, P < 0.021). M-ratios were low and indicative of a genetic bottleneck in the NWA (M-ratio = 0.71, P < 0.004) but not SA (M-ratio = 0.85, P = 0.39). This is consistent with other evidence showing a steep population decline occurring in the mid to late 20th century in the NWA, whereas the SA population appears to have been relatively stable. Estimates of effective population size ranged from 22.6 to 66.3 (NWA) and 188 to 1998.3 (SA) and evidence of inbreeding was found (primarily in NWA). Overall, our findings indicate that white population dynamics within NWA and SA are determined more by intrinsic reproduction than immigration and there is genetic evidence of a population decline in the NWA, further justifying the strong domestic protective measures that have been taken for this species in this region. Our study also highlights how assessment of genetic diversity can complement other sources of information to better understand the status of threatened marine fish populations.

Challenges in analysis and interpretation of microsatellite data for population genetic studies

Advancing technologies have facilitated the ever-widening application of genetic markers such as microsatellites into new systems and research questions in biology. In light of the data and experience accumulated from several years of using microsatellites, we present here a literature review that synthesizes the limitations of microsatellites in population genetic studies. With a focus on population structure, we review the widely used fixation (F ST) statistics and Bayesian clustering algorithms and find that the former can be confusing and problematic for microsatellites and that the latter may be confounded by complex population models and lack power in certain cases. Clustering, multivariate analyses, and diversity-based statistics are increasingly being applied to infer population structure, but in some instances these methods lack formalization with microsatellites. Migration-specific methods perform well only under narrow constraints. We also examine the use of microsatellites for inferring effective population size, changes in population size, and deeper demographic history, and find that these methods are untested and/or highly context-dependent. Overall, each method possesses important weaknesses for use with microsatellites, and there are significant constraints on inferences commonly made using microsatellite markers in the areas of population structure, admixture, and effective population size. To ameliorate and better understand these constraints, researchers are encouraged to analyze simulated datasets both prior to and following data collection and analysis, the latter of which is formalized within the approximate Bayesian computation framework. We also examine trends in the literature and show that microsatellites continue to be widely used, especially in non-human subject areas. This review assists with study design and molecular marker selection, facilitates sound interpretation of microsatellite data while fostering respect for their practical limitations, and identifies lessons that could be applied toward emerging markers and high-throughput technologies in population genetics.

Alternative reproductive tactics increase effective population size and decrease inbreeding in wild Atlantic salmon

While nonanadromous males (stream-resident and/or mature male parr) contribute to reproduction in anadromous salmonids, little is known about their impacts on key population genetic parameters. Here, we evaluated the contribution of Atlantic salmon mature male parr to the effective number of breeders (Nb) using both demographic (variance in reproductive success) and genetic (linkage disequilibrium) methods, the number of alleles, and the relatedness among breeders. We used a recently published pedigree reconstruction of a wild anadromous Atlantic salmon population in which 2548 fry born in 2010 were assigned parentage to 144 anadromous female and 101 anadromous females that returned to the river to spawn in 2009 and to 462 mature male parr. Demographic and genetic methods revealed that mature male parr increased population Nb by 1.79 and 1.85 times, respectively. Moreover, mature male parr boosted the number of alleles found among progenies. Finally, mature male parr were in average less related to anadromous females than were anadromous males, likely because of asynchronous sexual maturation between mature male parr and anadromous fish of a given cohort. By increasing Nb and allelic richness, and by decreasing inbreeding, the reproductive contribution of mature male parr has important evolutionary and conservation implications for declining Atlantic salmon populations.