Interannual variation in effective number of breeders and estimation of effective population size in long-lived iteroparous lake sturgeon (Acipenser fulvescens)

Individual-based analyses reveal effects of behavioral and demographic variables associated with multi-annual reproductive success of male and female lake sturgeon

Quantifying effects of individual attributes and population demographic characteristics that affect inter- and intrasexual interactions and adult reproductive success, and the spatial and temporal contexts in which they are expressed is important to effective species management. Multi-year individual-based analyses using genetically determined parentage allowed the examination of variables associated with the reproductive success of male and female lake sturgeon (Acipenser fulvescens) in the well-studied population in Black Lake, Michigan, USA. Spawning lake sturgeon (a total of 599 individuals where many were captured more than once based on 1024 total captures) and larvae (N = 3436) were genotyped during each of seven consecutive years (2001-2007). Factors associated with individual reproductive success differed between sexes and varied among spawning groups within a year and among years depending on spawning date (higher reproductive success earlier in the season for females) and spawning locations (higher reproductive success in upstream spawning zones for females). Female reproductive success increased nonlinearly with increasing body size. Male reproductive success increased with increasing residence time in spawning areas and, to a modest degree, with increasing body size in a nonlinear fashion. Fixed effects of repeatability in spawn timing and location across years led to consistently higher or lower reproductive success for females. Results identified factors, including time spent at spawning areas by males and intersexual encounters and mate number, that contributed to higher interindividual variance in reproductive success and affected population levels of recruitment, the degree of subpopulation genetic structure (lack of isolation by time), and effective population size.

Conservation Genetics of Lake Sturgeon (Acipenser fulvescens): Nuclear Phylogeography Drives Contemporary Patterns of Genetic Structure and Diversity

Sustainable management of exploited and endangered species is facilitated by knowledge of their geographic genetic structure. Lake sturgeon (Acipenser fulvescens) epitomizes both categories, but genetic information has largely been limited to the Laurentian Great Lakes basin. We assessed the hierarchical geographic genetic structure of lake sturgeon across their Canadian range using a variation at 14 microsatellite loci. Observed patterns showed evidence of two ancestral groups which originated from Mississippian and Missourian glacial refugia. Coalescent analysis indicates the two lineages most recently shared common ancestry during the late Pleistocene and were likely isolated by the late Wisconsinan ice advance, with subsequent interpopulation divergences within each lineage reflecting their reciprocal isolation as glacial meltwaters receded. Hierarchical patterns of genetic relationships among contemporary populations largely reflect colonization histories and connections within primary and secondary watersheds. Populations in western Canada showed strong similarities based on their shared Missourian origins and colonization from glacial Lake Agassiz. By contrast, populations in the Great Lakes–St. Lawrence River drainage were largely founded from a Mississippian source. Sturgeon populations in northern parts of Ontario and Quebec showed evidence of mixed ancestry from secondary contact between the two refugial groups through Holocene meltwater lakes. Within major watersheds, the strong similarity among geographically separate populations reflects their shared ancestry during postglacial colonization. The general lack of structure within major river systems highlights historically continuous habitat (connectivity) and gene flow rather than contemporary barriers (dams). These data highlight the importance of Quaternary and prehistoric events on patterns of genetic diversity and divergence within and among contemporary populations, as well as the importance of these populations for conserving the species’ evolutionary legacy.

Aquatic insects differentially affect lake sturgeon larval phenotypes and egg surface microbial communities

Documentation of how interactions among members of different stream communities [e.g., microbial communities and aquatic insect taxa exhibiting different feeding strategies (FS)] collectively influence the growth, survival, and recruitment of stream fishes is limited. Considerable spatial overlap exists between early life stages of stream fishes, including species of conservation concern like lake sturgeon (Acipenser fulvescens), and aquatic insects and microbial taxa that abundantly occupy substrates on which spawning occurs. Habitat overlap suggests that species interactions across trophic levels may be common, but outcomes of these interactions are poorly understood. We conducted an experiment where lake sturgeon eggs were fertilized and incubated in the presence of individuals from one of four aquatic insect FS taxa including predators, facultative and obligate-scrapers, collector-filterers/facultative predators, and a control (no insects). We quantified and compared the effects of different insect taxa on the taxonomic composition and relative abundance of egg surface bacterial and lower eukaryotic communities, egg size, incubation time to hatch, free embryo body size (total length) at hatch, yolk-sac area, (a measure of resource utilization), and percent survival to hatch. Mean egg size varied significantly among insect treatments. Eggs exposed to predators had a lower mean percent survival to hatch. Eggs exposed to predators had significantly shorter incubation periods. At hatch, free embryos exposed to predators had significantly smaller yolk sacs and total length. Multivariate analyses revealed that egg bacterial and lower eukaryotic surface community composition varied significantly among insect treatments and between time periods (1 vs 4 days post-fertilization). Quantitative PCR documented significant differences in bacterial 16S copy number, and thus abundance on egg surfaces varied across insect treatments. Results indicate that lethal and non-lethal effects associated with interactions between lake sturgeon eggs and free embryos and aquatic insects, particularly predators, contributed to lake sturgeon trait variability that may affect population levels of recruitment.

Population genomic monitoring provides insight into conservation status but no correlation with demographic estimates of extinction risk in a threatened trout

The current extinction crisis requires effective assessment and monitoring tools. Genetic approaches are appealing given the relative ease of field sampling required to estimate genetic diversity characteristics assumed related to population size, evolutionary potential, and extinction risk, and to evaluate hybridization with non‐native species simultaneously. However, linkages between population genetic metrics of diversity from survey‐style field collections and demographic estimates of population size and extinction risk are still in need of empirical examples, especially for remotely distributed species of conservation concern where the approach might be most beneficial. We capitalized on an exceptional opportunity to evaluate congruence between genetic diversity metrics and demographic‐based estimates of abundance and extinction risk from a comprehensive Multiple Population Viability Analysis (MPVA) in a threatened fish, the Lahontan cutthroat trout (LCT). We sequenced non‐native trout reference samples and recently collected and archived tissue samples of most remaining populations of LCT (N = 60) and estimated common genetic assessment metrics, predicting minimal hybridization with non‐native trout, low diversity, and declining diversity over time. We further hypothesized genetic metrics would correlate positively with MPVA‐estimated abundance and negatively with extinction probability. We uncovered several instances of hybridization that pointed to immediate management needs. After removing hybridized individuals, cautious interpretation of low effective population sizes (2–63) suggested reduced evolutionary potential for many LCT populations. Other genetic metrics did not decline over time nor correlate with MPVA‐based estimates of harmonic mean abundance or 30‐year extinction probability. Our results demonstrate benefits of genetic monitoring for efficiently detecting hybridization and, though genetic results were disconnected from demographic assessment of conservation status, they suggest reduced evolutionary potential and likely a higher conservation risk than currently recognized for this threatened fish. We emphasize that genetic information provides essential complementary insight, in addition to demographic information, for evaluating species status.

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.

RAPTURE (RAD capture) panel facilitates analyses characterizing sea lamprey reproductive ecology and movement dynamics

Genomic tools are lacking for invasive and native populations of sea lamprey (Petromyzon marinus). Our objective was to discover single nucleotide polymorphism (SNP) loci to conduct pedigree analyses to quantify reproductive contributions of adult sea lampreys and dispersion of sibling larval sea lampreys of different ages in Great Lakes tributaries. Additional applications of data were explored using additional geographically expansive samples. We used restriction site-associated DNA sequencing (RAD-Seq) to discover genetic variation in Duffins Creek (DC), Ontario, Canada, and the St. Clair River (SCR), Michigan, USA. We subsequently developed RAD capture baits to genotype 3,446 RAD loci that contained 11,970 SNPs. Based on RAD capture assays, estimates of variance in SNP allele frequency among five Great Lakes tributary populations (mean F ST 0.008; range 0.00-0.018) were concordant with previous microsatellite-based studies; however, outlier loci were identified that contributed substantially to spatial population genetic structure. At finer scales within streams, simulations indicated that accuracy in genetic pedigree reconstruction was high when 200 or 500 independent loci were used, even in situations of high spawner abundance (e.g., 1,000 adults). Based on empirical collections of larval sea lamprey genotypes, we found that age-1 and age-2 families of full and half-siblings were widely but nonrandomly distributed within stream reaches sampled. Using the genomic scale set of SNP loci developed in this study, biologists can rapidly genotype sea lamprey in non-native and native ranges to investigate questions pertaining to population structuring and reproductive ecology at previously unattainable scales.

Improved genetic identification of acipenseriform embryos with application to the endangered pallid sturgeon Scaphirhynchus albus

We produced pallid sturgeon Scaphirhynchus albus embryos at five pre-hatch developmental stages and isolated and quantified genomic DNA from four of the stages using four commercial DNA isolation kits. Genomic DNA prepared using the kit that produced the largest yields and concentrations were used for microsatellite DNA analyses of 10-20 embryos at each of the five developmental stages. We attempted to genotype the hatchery-produced embryos at 19 microsatellite loci and confirmed reliable genotyping by comparing the microsatellite genotypes to those of known parents. Embryos at stages 5 and 8 did not produce reliable genotyping while those at stages 14, 24 and 33 did. We used the same DNA isolation method on 262 wild-caught acipenseriform embryos collected from the lower Yellowstone River. A total of 200 of the wild embryos were successfully identified to stages 8 to 34 and the rest could not be staged. Using a combination of single nucleotide polymorphism and microsatellite markers, 249 of the wild-caught embryos were genetically identified as paddlefish Polyodon spathula, five were identified as shovelnose sturgeon Scaphirhynchus platorynchus and eight failed to amplify. None were identified as pallid sturgeon. This study demonstrates that early-stage wild-spawned acipenseriform embryos can be genetically identified less than 24 h post-spawn. This methodology will be useful for recovery efforts for endangered pallid sturgeon and can be applied to other acipenseriform species.

Genetic evidence supports polygamous mating system in a wild population of Prochilodus lineatus (Characiformes: Prochilodontidae), a Neotropical shoal spawner fish

ABSTRACT Behavioral observations made on fish have revealed remarkably diverse reproductive strategies, including polygamy by both sexes. Still, to date, most Neotropical species remain unstudied as to whether the observed reproductive behavior in natural populations correlates with their genetic mating systems. Here, we investigated the genetic mating system of a wild population of Prochilodus lineatus settled in the Middle Uruguay River basin. By using sibship reconstruction and parental inference methods based on microsatellites’ genotypes, we inferred 45 females and 47 males as potential parents of the 87 larvae analyzed. We found evidence supporting polygamous mating in both sexes: while a high percentage of males (44.7%) fertilized the eggs of one female, 55.3% of the inferred males fertilized eggs of up to four females. Likewise, while 44.5% of the inferred females had their eggs fertilized by one only male, 55.5% of females were fertilized by multiple males. The estimated proxy of the effective population size (Nb) was 126, exhibiting moderate to high levels of genetic diversity. The genetic evidence contributed in this study complements earlier behavioral observations of formation of spawning nuclei of aggregating breeders, which may be promoting a polygamous mating strategy in this long-distance migratory fish.

Molecular diet analysis reveals predator-prey community dynamics and environmental factors affecting predation of larval lake sturgeon Acipenser fulvescens in a natural system

This study utilized molecular tools to quantify the prevalence of predation during the vulnerable drifting larval life-history stage of lake sturgeon Acipenser fulvescens. How predators, the co-distributed prey community and abiotic environmental conditions (e.g., stream substrata) affected predation levels was quantified. Nightly D-frame drift net surveys were used to estimate the biomass of A. fulvescens and co-distributed prey. Gastrointestinal diet samples (n = 1,140) from 28 species of potential fish predators were collected during electrofishing surveys. Sampling was conducted for 17 days across 2015 and 2016. Based on DNA barcode analysis using sturgeon-specific mtDNA cytochrome oxidase I primers, A. fulvescens DNA was detected in 73 of 1,140 diet samples (6.40%) from 16 of the 28 predator species examined. A logistic regression model was used to analyse the effects of biotic and abiotic variables associated with the likelihood a predator had consumed larval A. fulvescens. Increasing lunar illumination and biomass of larval A. fulvescens increased predation rates on larval A. fulvescens. Higher discharge and greater biomass and proportions of alternative prey decreased predation rates of larval A. fulvescens. Predation rates were slightly higher in habitats with sand substrata. Most predator species preyed upon larval A. fulvescens at similar rates. The study revealed considerably higher incidence of predation on larval A. fulvescens than previous studies had documented using traditional morphological diet analysis. Co-distributed prey and abiotic environmental variables that affected the predation rates of a species of regional conservation concern can inform future management actions.

Genetic substructure and admixture as important factors in linkage disequilibrium-based estimation of effective number of breeders in recovering wildlife populations

The number of effective breeders (Nb ) and effective population size (Ne ) are population parameters reflective of evolutionary potential, susceptibility to stochasticity, and viability. We have estimated these parameters using the linkage disequilibrium-based approach with LDNE through the latest phase of population recovery of the brown bears (Ursus arctos) in Finland (1993-2010; N = 621). This phase of the recovery was recently documented to be associated with major changes in genetic composition. In particular, differentiation between the northern and the southern genetic cluster declined rapidly within 1.5 generations. Based on this, we have studied effects of the changing genetic structure on Nb and Ne , by comparing estimates for whole Finland with the estimates for the two genetic clusters. We expected a potentially strong relationship between estimate sizes and genetic differentiation, which should disappear as the population recovers and clusters merge. Consistent with this, our estimates for whole Finland were lower than the sum of the estimates of the two genetic clusters and both approaches produced similar estimates in the end. Notably, we also found that admixed genotypes strongly increased the estimates. In all analyses, our estimates for Ne were larger than Nb and likely reflective for brown bears of the larger region of Finland and northwestern Russia. Conclusively, we find that neglecting genetic substructure may lead to a massive underestimation of Nb and Ne . Our results also suggest the need for further empirical analysis focusing on individuals with admixed genotypes and their potential high influence on Nb and Ne .

Rethinking the influence of hydroelectric development on gene flow in a long-lived fish, the Lake Sturgeon Acipenser fulvescens

Many hydroelectric dams have been in place for 50 - >100 years, which for most fish species means that enough generations have passed for fragmentation induced divergence to have accumulated. However, for long-lived species such as Lake Sturgeon, Acipenser fulvescens, it should be possible to discriminate between historical population structuring and contemporary gene flow and improve the broader understanding of anthropogenic influence. On the Winnipeg River, Manitoba, two hypotheses were tested: 1) Measureable quantities of former reservoir dwelling Lake Sturgeon now reside downstream of the Slave Falls Generating Station, and 2) genetically differentiated populations of Lake Sturgeon occur upstream and downstream, a result of historical structuring. Genetic methods based on ten microsatellite markers were employed, and simulations were conducted to provide context. With regards to contemporary upstream to downstream contributions, the inclusion of length-at-age data proved informative. Both pairwise relatedness and Bayesian clustering analysis substantiated that fast-growing outliers, apparently entrained after residing in the upstream reservoir for several years, accounted for ~15% of the Lake Sturgeon 525–750 mm fork length captured downstream. With regards to historical structuring, upstream and downstream populations were found to be differentiated (F_ST = 0.011, and 0.013–0.014 when fast-growing outliers were excluded), and heterozygosity metrics were higher for downstream versus upstream juveniles. Historical asymmetric (downstream) gene flow in the vicinity of the generating station was the most logical explanation for the observed genetic structuring. In this section of the Winnipeg River, construction of a major dam does not appear to have fragmented a previously panmictic Lake Sturgeon population, but alterations to habitat may be influencing upstream to downstream contributions in unexpected ways.

Highly Connected Populations and Temporal Stability in Allelic Frequencies of a Harvested Crab from the Southern Pacific Coast

For marine invertebrates with a benthic adult form and a planktonic larva phase, the connectivity among populations is mainly based on larval dispersal. While an extended larval phase will promote gene flow, other factors such as an intensive fishery and geographical barriers could lead to changes in genetic variability. In this study, the population genetic structure of the commercial crab Metacarcinus edwardsii was analyzed along 700 km of the Chilean coast. The analysis, based on eight microsatellite loci genotyped from megalopae and adult crabs, considered temporal and spatial patterns of genetic variation. The results showed no evidence of spatial patterns in genetic structure, suggesting high connectivity among the sampling sites. The temporal analysis showed no evidence of changes in allele frequencies and no evidence of a recent bottleneck. The lack of spatial structure and allele variation over time could be explained by the interaction of factors such as i) low reproductive variance due to the capability of females to store sperm in the seminal receptacle, which can be used for successive broods, ii) high larval dispersal and iii) high individual reproductive output. Using our data as priors, a genetic modelling approach coincided, predicting this temporal and spatial stability. The same analysis showed that a reduction in population size leads to the loss of genetic variability in populations, as well as of the genetic cohesiveness between populations, pointing out the importance management for species under exploitation, such as M. edwardsii.

Making sense of the relationships between Ne, Nb and Nc towards defining conservation thresholds in Atlantic salmon (Salmo salar)

Effective population size over a generation (Ne) or over a reproductive cycle (Nb) and the adult census size (Nc) are important parameters in both conservation and evolutionary biology. Ne provides information regarding the rate of loss of genetic diversity and can be tracked back in time to infer demographic history of populations, whereas Nb may often be more easily quantified than Nc for short-term abundance monitoring. In this study, we propose (1) an empirical context to Waples et al. (2014) who introduced a correction to bias due to overlapping generations, and (2) a mathematical relationship between Ne and Nb for direct application in Atlantic salmon populations in Québec, Canada. To achieve this, we investigate the relationships between Ne, Nb and Nc in 10 Atlantic salmon populations, Canada, for which we genotyped 100 randomly sampled young-of-the year individuals for 5 consecutive years. The results show a positive correlation between Ne, Nb and Nc, suggesting that Nb is an indicative parameter for tracking effective population size and abundance of Atlantic salmon. However, our model allows predicting Nc from Nb values at 27% that can be partly explained by high variance in Nb/Nc both among populations (37%) and among years (19%). This result illustrates the need for thorough calibration of Nb/Nc before using Nb in monitoring programs, as well as a full understanding of the limits of such an approach. Finally, we discuss the importance of these results for the management of wild populations.

RNA-sequencing of the sturgeon Acipenser baeri provides insights into expression dynamics of morphogenic differentiation and developmental regulatory genes in early versus late developmental stages

Background

Acipenser baeri, one of the critically endangered animals on the verge of extinction, is a key species for evolutionary, developmental, physiology and conservation studies and a standout amongst the most important food products worldwide. Though the transcriptome of the early development of A. baeri has been published recently, the transcriptome changes occurring in the transition from embryonic to late stages are still unknown. The aim of this work was to analyze the transcriptomes of embryonic and post-embryonic stages of A. baeri and identify differentially expressed genes (DEGs) and their expression patterns using mRNA collected from specimens at big yolk plug, wide neural plate and 64 day old sturgeon developmental stages for RNA-Seq.

Results

The paired-end sequencing of the transcriptome of samples of A. baeri collected at two early (big yolk plug (T1, 32 h after fertilization) and wide neural plate formation (T2, 45 h after fertilization)) and one late (T22, 64 day old sturgeon) developmental stages using Illumina Hiseq2000 platform generated 64039846, 64635214 and 75293762 clean paired-end reads for T1, T2 and T22, respectively. After quality control, the sequencing reads were de novo assembled to generate a set of 149,265 unigenes with N50 value of 1277 bp. Functional annotation indicated that a substantial number of these unigenes had significant similarity with proteins in public databases. Differential expression profiling allowed the identification of 2789, 12,819 and 10,824 DEGs from the respective T1 vs. T2, T1 vs. T22 and T2 vs. T22 comparisons. High correlation of DEGs’ features was recorded among early stages while significant divergences were observed when comparing the late stage with early stages. GO and KEGG enrichment analyses revealed the biological processes, cellular component, molecular functions and metabolic pathways associated with identified DEGs. The qRT-PCR performed for candidate genes in specimens confirmed the validity of the RNA-seq data.

Conclusions

This study presents, for the first time, an extensive overview of RNA-Seq based characterization of the early and post-embryonic developmental transcriptomes of A. baeri and provided 149,265 gene sequences that will be potentially valuable for future molecular and genetic studies in A. baeri.

Electronic supplementary material

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

Retention of Ancestral Genetic Variation Across Life-Stages of an Endangered, Long-Lived Iteroparous Fish

As with many endangered, long-lived iteroparous fishes, survival of razorback sucker depends on a management strategy that circumvents recruitment failure that results from predation by non-native fishes. In Lake Mohave, AZ-NV, management of razorback sucker centers on capture of larvae spawned in the lake, rearing them in off-channel habitats, and subsequent release ("repatriation") to the lake when adults are sufficiently large to resist predation. The effects of this strategy on genetic diversity, however, remained uncertain. After correction for differences in sample size among groups, metrics of mitochondrial DNA (mtDNA; number of haplotypes, N H , and haplotype diversity, H D ) and microsatellite (number of alleles, N A , and expected heterozygosity, H E ) diversity did not differ significantly between annual samples of repatriated adults and larval year-classes or among pooled samples of repatriated adults, larvae, and wild fish. These findings indicate that the current management program thus far maintained historical genetic variation of razorback sucker in the lake. Because effective population size, N e , is closely tied to the small census population size (N c = ~1500-3000) of razorback sucker in Lake Mohave, this population will remain at risk from genetic, as well as demographic risk of extinction unless N c is increased substantially.

Effects of overlapping generations on linkage disequilibrium estimates of effective population size

Use of single-sample genetic methods to estimate effective population size has skyrocketed in recent years. Although the underlying models assume discrete generations, they are widely applied to age-structured species. We simulated genetic data for 21 iteroparous animal and plant species to evaluate two untested hypotheses regarding performance of the single-sample method based on linkage disequilibrium (LD): (1) estimates based on single-cohort samples reflect the effective number of breeders in one reproductive cycle (Nb), and (2) mixed-age samples reflect the effective size per generation (Ne). We calculated true Ne and Nb, using the model species' vital rates, and verified these with individual-based simulations. We show that single-cohort samples should be equally influenced by Nb and Ne and confirm this with simulated results: [Formula: see text] was a linear (r(2) = 0.98) function of the harmonic mean of Ne and Nb. We provide a quantitative bias correction for raw [Formula: see text] based on the ratio Nb/Ne, which can be estimated from two or three simple life history traits. Bias-adjusted estimates were within 5% of true Nb for all 21 study species and proved robust when challenged with new data. Mixed-age adult samples produced downwardly biased estimates in all species, which we attribute to a two-locus Wahlund effect (mixture LD) caused by combining parents from different cohorts in a single sample. Results from this study will facilitate interpretation of rapidly accumulating genetic estimates in terms of both Ne (which influences long-term evolutionary processes) and Nb (which is more important for understanding eco-evolutionary dynamics and mating systems).

Time-series analysis reveals genetic responses to intensive management of razorback sucker (Xyrauchen texanus)

Time-series analysis is used widely in ecology to study complex phenomena and may have considerable potential to clarify relationships of genetic and demographic processes in natural and exploited populations. We explored the utility of this approach to evaluate population responses to management in razorback sucker, a long-lived and fecund, but declining freshwater fish species. A core population in Lake Mohave (Arizona-Nevada, USA) has experienced no natural recruitment for decades and is maintained by harvesting naturally produced larvae from the lake, rearing them in protective custody, and repatriating them at sizes less vulnerable to predation. Analyses of mtDNA and 15 microsatellites characterized for sequential larval cohorts collected over a 15-year time series revealed no changes in geographic structuring but indicated significant increase in mtDNA diversity for the entire population over time. Likewise, ratios of annual effective breeders to annual census size (N b /N a) increased significantly despite sevenfold reduction of N a. These results indicated that conservation actions diminished near-term extinction risk due to genetic factors and should now focus on increasing numbers of fish in Lake Mohave to ameliorate longer-term risks. More generally, time-series analysis permitted robust testing of trends in genetic diversity, despite low precision of some metrics.

What remains from a 454 run: estimation of success rates of microsatellite loci development in selected newt species (Calotriton asper, Lissotriton helveticus, and Triturus cristatus) and comparison with Illumina-based approaches

The development of microsatellite loci has become more efficient using next-generation sequencing (NGS) approaches, and many studies imply that the amount of applicable loci is large. However, few studies have sought to quantify the number of loci that are retained for use out of the thousands of sequence reads initially obtained. We analyzed the success rate of microsatellite loci development for three amphibian species using a 454 NGS approach on tetra-nucleotide motif-enriched species-specific libraries. The number of sequence reads obtained differed strongly between species and ranged from 19,562 for Triturus cristatus to 55,626 for Lissotriton helveticus, with 52,075 reads obtained for Calotriton asper. PHOBOS was used to identify sequences with tetra-nucleotide repeat motifs with a minimum repeat number of ten and high quality primer binding sites. Of 107 sequences for T. cristatus, 316 for C. asper and 319 for L. helveticus, we tested the amplification success, polymorphism, and degree of heterozygosity for 41 primer combinations each for C. asper and T. cristatus, and 22 for L. helveticus. We found 11 polymorphic loci for T. cristatus, 20 loci for C. asper, and 15 loci for L. helveticus. Extrapolated, the number of potentially amplifiable loci (PALs) resulted in estimated species-specific success rates of 0.15% (T. cristatus), 0.30% (C. asper), and 0.39% (L. helveticus). Compared with representative Illumina NGS approaches, our applied 454-sequencing approach on specifically enriched sublibraries proved to be quite competitive in terms of success rates and number of finally applicable loci.