neogen: A tool to predict genetic effective population size (Ne ) for species with generational overlap and to assist empirical Ne study design

Estimation of effective number of breeders and effective population size in an abundant and heavily exploited marine teleost

Obtaining reliable estimates of the effective number of breeders (N b) and generational effective population size (N e) for fishery-important species is challenging because they are often iteroparous and highly abundant, which can lead to bias and imprecision. However, recent advances in understanding of these parameters, as well as the development of bias correction methods, have improved the capacity to generate reliable estimates. We utilized samples of both single-cohort young of the year and mixed-age adults from two geographically and genetically isolated stocks of the Australasian snapper (Chrysophrys auratus) to investigate the feasibility of generating reliable N b and N e estimates for a fishery species. Snapper is an abundant, iteroparous broadcast spawning teleost that is heavily exploited by recreational and commercial fisheries. Employing neutral genome-wide SNPs and the linkage-disequilibrium method, we determined that the most reliable N b and N e estimates could be derived by genotyping at least 200 individuals from a single cohort. Although our estimates made from the mixed-age adult samples were generally lower and less precise than those based on a single cohort, they still proved useful for understanding relative differences in genetic effective size between stocks. The correction formulas applied to adjust for biases due to physical linkage of loci and age structure resulted in substantial upward modifications of our estimates, demonstrating the importance of applying these bias corrections. Our findings provide important guidelines for estimating N b and N e for iteroparous species with large populations. This work also highlights the utility of samples originally collected for stock structure and stock assessment work for investigating genetic effective size in fishery-important species.

Genome-wide analysis reveals the genetic stock structure of hoki (Macruronus novaezelandiae)

The assessment of the genetic structuring of biodiversity is crucial for management and conservation. This is particularly critical for widely distributed and highly mobile deep-water teleosts, such as hoki (Macruronus novaezelandiae). This species is significant to Māori people and supports the largest commercial fishery in New Zealand, but uncertainty about its stock structure presents a challenge for management. Here, we apply a comprehensive genomic analysis to shed light on the demographic structure of this species by (1) assembling the genome, (2) generating a catalogue of genome-wide SNPs to infer the stock structure and (3) identifying regions of the genome under selection. The final genome assembly used short and long reads and is near complete, representing 93.8% of BUSCO genes, and consisting of 566 contigs totalling 501 Mb. Whole-genome re-sequencing of 510 hoki sampled from 14 locations around New Zealand and Australia, at a read depth greater than 10×, produced 227,490 filtered SNPs. Analyses of these SNPs were able to resolve the stock structure of hoki into two genetically and geographically distinct clusters, one including the Australian and the other one all New Zealand locations, indicating genetic exchange between these regions is limited. Location differences within New Zealand samples were much more subtle (global F ST  = 0.0006), and while small and significant differences could be detected, they did not conclusively identify additional substructures. Ten putative adaptive SNPs were detected within the New Zealand samples, but these also did not geographically partition the dataset further. Contemporary and historical N e estimation suggest the current New Zealand population of hoki is large yet declining. Overall, our study provides the first genomic resources for hoki and provides detailed insights into the fine-scale population genetic structure to inform the management of this species.

Population Genomics Training for the Next Generation of Conservation Geneticists: ConGen 2018 Workshop

The increasing availability and complexity of next-generation sequencing (NGS) data sets make ongoing training an essential component of conservation and population genetics research. A workshop entitled “ConGen 2018” was recently held to train researchers in conceptual and practical aspects of NGS data production and analysis for conservation and ecological applications. Sixteen instructors provided helpful lectures, discussions, and hands-on exercises regarding how to plan, produce, and analyze data for many important research questions. Lecture topics ranged from understanding probabilistic (e.g., Bayesian) genotype calling to the detection of local adaptation signatures from genomic, transcriptomic, and epigenomic data. We report on progress in addressing central questions of conservation genomics, advances in NGS data analysis, the potential for genomic tools to assess adaptive capacity, and strategies for training the next generation of conservation genomicists.

Big Data in Conservation Genomics: Boosting Skills, Hedging Bets, and Staying Current in the Field

A current challenge in the fields of evolutionary, ecological, and conservation genomics is balancing production of large-scale datasets with additional training often required to handle such datasets. Thus, there is an increasing need for conservation geneticists to continually learn and train to stay up-to-date through avenues such as symposia, meetings, and workshops. The ConGen meeting is a near-annual workshop that strives to guide participants in understanding population genetics principles, study design, data processing, analysis, interpretation, and applications to real-world conservation issues. Each year of ConGen gathers a diverse set of instructors, students, and resulting lectures, hands-on sessions, and discussions. Here, we summarize key lessons learned from the 2019 meeting and more recent updates to the field with a focus on big data in conservation genomics. First, we highlight classical and contemporary issues in study design that are especially relevant to working with big datasets, including the intricacies of data filtering. We next emphasize the importance of building analytical skills and simulating data, and how these skills have applications within and outside of conservation genetics careers. We also highlight recent technological advances and novel applications to conservation of wild populations. Finally, we provide data and recommendations to support ongoing efforts by ConGen organizers and instructors-and beyond-to increase participation of underrepresented minorities in conservation and eco-evolutionary sciences. The future success of conservation genetics requires both continual training in handling big data and a diverse group of people and approaches to tackle key issues, including the global biodiversity-loss crisis.

A globally threatened shark, Carcharias taurus, shows no population decline in South Africa

Knowledge about the demographic histories of natural populations helps to evaluate their conservation status, and potential impacts of natural and anthropogenic pressures. In particular, estimates of effective population size obtained through molecular data can provide useful information to guide management decisions for vulnerable populations. The spotted ragged-tooth shark, Carcharias taurus (also known as the sandtiger or grey nurse shark), is widely distributed in warm-temperate and subtropical waters, but has suffered severe population declines across much of its range as a result of overexploitation. Here, we used multilocus genotype data to investigate the demographic history of the South African C. taurus population. Using approximate Bayesian computation and likelihood-based importance sampling, we found that the population underwent a historical range expansion that may have been linked to climatic changes during the late Pleistocene. There was no evidence for a recent anthropogenic decline. Together with census data suggesting a stable population, these results support the idea that fishing pressure and other threats have so far not been detrimental to the local C. taurus population. The results reported here indicate that South Africa could possibly harbour the last remaining, relatively pristine population of this widespread but vulnerable top predator.

Genome-wide SNPs detect no evidence of genetic population structure for reef manta rays (Mobula alfredi) in southern Mozambique

Little is known about the extent of genetic connectivity along continuous coastlines in manta rays, or whether site visitation is influenced by relatedness. Such information is pertinent to defining population boundaries and understanding localized dispersal patterns and behaviour. Here, we use 3057 genome-wide single-nucleotide polymorphisms (SNPs) to evaluate population genetic structure and assess the levels of relatedness at aggregation sites of reef manta rays (Mobula alfredi) in southern Mozambique (n = 114). Contrary to indications of limited dispersal along the southern Mozambican coastline inferred from photo-identification and telemetry studies, our results show no evidence of population structure (non-significant F_ST < 0.001) for M. alfredi along this coast. We also found no evidence that individuals sampled at the same site were more related than expected by chance for males, females or across both sexes, suggesting that kinship may not influence visitation patterns at these sites. We estimated the effective population size (N_e) of this population to be 375 (95% CI = 369-380). Comparison to a distant eastern Indian Ocean site (Western Australia, n = 15) revealed strong genetic differentiation between Mozambique and Western Australia (F_ST = 0.377), identifying the Indian Ocean basin as a barrier to dispersal. Our findings show that genetic connectivity in M. alfredi extends for several hundred kilometres along continuous coastlines. We therefore recommend that the population in Mozambique be considered a discrete management unit, and future conservation plans should prioritize integrated strategies along the entire southern coastline.

AgeStrucNb: Software for Simulating and Detecting Changes in the Effective Number of Breeders (Nb)

Estimation of the effective number of breeders per reproductive event (Nb) using single sample DNA-marker-based methods has rapidly grown in recent years. However, estimating Nb is difficult in age-structured populations because the performance of estimators is influenced by the Nb / Ne ratio, which varies among species with different life histories. We provide a computer program, AgeStrucNb, to simulate age-structured populations (including life history) and also estimate Nb. The AgeStrucNb program is composed of 4 major components to simulate, subsample, estimate, and then visualize Nb time series data. AgeStrucNb allows users to also quantify the precision and accuracy of any set of loci or sample size to estimate Nb for many species and populations. AgeStrucNb allows users to conduct power analysis to evaluate sensitivity to detect changes in Nb or the power to detect a correlation between trends in Nb and environmental variables (e.g., temperature, habitat quality, predator or pathogen abundance) that could be driving changes in Nb. The software provides Nb estimates for empirical data sets using the LDNe (linkage disequilibrium) method, includes publication-quality output graphs, and outputs genotype files in Genepop format for use in other programs. AgeStrucNb will help advance the application of genetic markers for monitoring Nb, which will help biologists to detect population declines and growth, which is crucial for research and conservation of natural and managed populations.

Sweepstakes reproductive success is absent in a New Zealand snapper (Chrysophrus auratus) population protected from fishing despite "tiny" Ne /N ratios elsewhere

A landmark study published in 2002 estimated a very small N_e /N ratio (around 10^-5 ) in a population of pink snapper (Chrysophrys auratus, Forster, 1801) in the Hauraki Gulf in New Zealand. It epitomized the tiny N_e /N ratios (<10^-3 ) reported in marine species due to the hypothesized operation of sweepstakes reproductive success (SRS). Here we re-evaluate the occurrence of SRS in marine species and the potential effect of fishing on the N_e /N ratio by studying the same species in the same region, but in a population that has been protected from fishing since 1975. We combine empirical, simulation and model-based approaches to estimate N_e (and N_b ) from genotypes of 1,044 adult fish and estimate N using recapture-probabilities. The estimated N_e /N ratio was much larger (0.33, SE: 0.14) than expected. The magnitude of estimates of population-wide variance in individual lifetime reproductive success (10-18) suggested that the sweepstakes effect was negligible in the study population. After evaluating factors that could explain the contrast between studies - experimental design, life history differences, environmental effects and the influence of exploitation on the N_e /N ratio - we conclude that the low N_e of the Hauraki Gulf population is associated with demographic instability in the harvested compared to the protected population despite circumstantial evidence that the 2002 study may have underestimated N_e . This study has broad implications for the prevailing view that reproductive success in the sea is largely driven by chance, and for genetic monitoring of populations using the N_e /N ratio and N_b .