Genetic estimates of contemporary effective population size in an endangered butterfly indicate a possible role for genetic compensation

Demographic resilience of brook trout populations subjected to experimental size-selective harvesting

Sustainable management of exploited populations benefits from integrating demographic and genetic considerations into assessments, as both play a role in determining harvest yields and population persistence. This is especially important in populations subject to size-selective harvest, because size selective harvesting has the potential to result in significant demographic, life-history, and genetic changes. We investigated harvest-induced changes in the effective number of breeders (N ^ b ) for introduced brook trout populations (Salvelinus fontinalis) in alpine lakes from western Canada. Three populations were subject to 3 years of size-selective harvesting, while three control populations experienced no harvest. TheN ^ c decreased consistently across all harvested populations (on average 60.8%) but fluctuated in control populations. There were no consistent changes inN ^ b between control or harvest populations, but one harvest population experienced a decrease inN ^ b of 63.2%. TheN ^ b /N ^ c ratio increased consistently across harvest lakes; however we found no evidence of genetic compensation (where variance in reproductive success decreases at lower abundance) based on changes in family evenness (FE ^ ) and the number of full-sibling families (N ^ fam ). We found no relationship betweenFE ^ andN ^ c or betweenN ^ fam /N ^ c andFE ^ . We posit that change inN ^ b was buffered by constraints on breeding habitat prior to harvest, such that the same number of breeding sites were occupied before and after harvest. These results suggest that effective size in harvested populations may be resilient to considerable changes in Nc in the short-term, but it is still important to monitor exploited populations to assess the risk of inbreeding and ensure their long-term survival.

Temporal Analysis of Effective Population Size and Mating System in a Social Wasp

Highly social species are successful because they cooperate in obligately integrated societies. We examined temporal genetic variation in the eusocial wasp Vespula maculifrons to gain a greater understanding of evolution in highly social taxa. First, we wished to test if effective population sizes of eusocial species were relatively low due to the reproductive division of labor that characterizes eusocial taxa. We thus estimated the effective population size of V. maculifrons by examining temporal changes in population allele frequencies. We sampled the genetic composition of a V. maculifrons population at 3 separate timepoints spanning a 13-year period. We found that effective population size ranged in the hundreds of individuals, which is similar to estimates in other, non-eusocial taxa. Second, we estimated levels of polyandry in V. maculifrons in different years to determine if queen mating system varied over time. We found no significant change in the number or skew of males mated to queens. In addition, mating skew was not significant within V. maculifrons colonies. Therefore, our data suggest that queen mate number may be subject to stabilizing selection in this taxon. Overall, our study provides novel insight into the selective processes operating in eusocial species by analyzing temporal genetic changes within populations.

Using Noninvasive Genetic Sampling to Survey Rare Butterfly Populations

Advances in nondestructive genetic sampling techniques continue to offer new opportunities for studying organisms, particularly those of conservation concern where more traditional invasive sampling methods are often not available. As part of a proof-of-concept, we investigated the effectiveness of using the chorion from residual butterfly egg debris as a source of viable genetic material for analysis. Laboratory material from a captive breeding population of the federally endangered Miami blue butterfly (Cyclargus thomasi bethunebakeri) was used to test efficacy and refine the methodology. The resulting best practices were subsequently evaluated using field-collected material from extant north Florida populations of the at-risk frosted elfin butterfly (Callophyrs irus). Our results demonstrated that it is possible to extract DNA of sufficiently high quantity and quality for successful gene sequencing. We additionally describe a simple, low-cost, and reliable method of collecting and storing egg debris samples that can be consistently adopted for field or laboratory work as well as deployed with projects that have a larger geographic scope and/or involve citizen scientists. Potential limitations related to field sample collection are discussed as well as needs for future evaluation.

Conservation Genomics of the Declining North American Bumblebee Bombus terricola Reveals Inbreeding and Selection on Immune Genes

The yellow-banded bumblebee Bombus terricola was common in North America but has recently declined and is now on the IUCN Red List of threatened species. The causes of B. terricola's decline are not well understood. Our objectives were to create a partial genome and then use this to estimate population data of conservation interest, and to determine whether genes showing signs of recent selection suggest a specific cause of decline. First, we generated a draft partial genome (contig set) for B. terricola, sequenced using Pacific Biosciences RS II at an average depth of 35×. Second, we sequenced the individual genomes of 22 bumblebee gynes from Ontario and Quebec using Illumina HiSeq 2500, each at an average depth of 20×, which were used to improve the PacBio genome calls and for population genetic analyses. The latter revealed that several samples had long runs of homozygosity, and individuals had high inbreeding coefficient F, consistent with low effective population size. Our data suggest that B. terricola's effective population size has decreased orders of magnitude from pre-Holocene levels. We carried out tests of selection to identify genes that may have played a role in ameliorating environmental stressors underlying B. terricola's decline. Several immune-related genes have signatures of recent positive selection, which is consistent with the pathogen-spillover hypothesis for B. terricola's decline. The new B. terricola contig set can help solve the mystery of bumblebee decline by enabling functional genomics research to directly assess the health of pollinators and identify the stressors causing declines.

Robust estimates of a high Ne/N ratio in a top marine predator, southern bluefin tuna

Genetic studies of several marine species with high fecundity have produced "tiny" estimates (≤10-3) of the ratio of effective population size (Ne) to adult census size (N), suggesting that even very large populations might be at genetic risk. A recent study using close-kin mark-recapture methods estimated adult abundance at N ≈ 2 × 106 for southern bluefin tuna (SBT), a highly fecund top predator that supports a lucrative (~$1 billion/year) fishery. We used the same genetic and life history data (almost 13,000 fish collected over 5 years) to generate genetic and demographic estimates of Ne per generation and Nb (effective number of breeders) per year and the Ne/N ratio. Demographic estimates, which accounted for age-specific vital rates, skip breeding, variation in fecundity at age, and persistent individual differences in reproductive success, suggest that Ne/N is >0.1 and perhaps about 0.5. The genetic estimates supported this conclusion. Simulations using true Ne = 5 × 105 (Ne/N = 0.25) produced results statistically consistent with the empirical genetic estimates, whereas simulations using Ne = 2 × 104 (Ne/N = 0.01) did not. Our results show that robust estimates of Ne and Ne/N can be obtained for large populations, provided sufficiently large numbers of individuals and genetic markers are used and temporal replication (here, 5 years of adult and juvenile samples) is sufficient to provide a distribution of estimates. The high estimated Ne/N ratio in SBT is encouraging and suggests that the species will not be compromised by a lack of genetic diversity in responding to environmental change and harvest.

Effective size of density-dependent two-sex populations: the effect of mating systems

Density dependence in vital rates is a key feature affecting temporal fluctuations of natural populations. This has important implications for the rate of random genetic drift. Mating systems also greatly affect effective population sizes, but knowledge of how mating system and density regulation interact to affect random genetic drift is poor. Using theoretical models and simulations, we compare N_e in short-lived, density-dependent animal populations with different mating systems. We study the impact of a fluctuating, density-dependent sex ratio and consider both a stable and a fluctuating environment. We find a negative relationship between annual N_e /N and adult population size N due to density dependence, suggesting that loss of genetic variation is reduced at small densities. The magnitude of this decrease was affected by mating system and life history. A male-biased, density-dependent sex ratio reduces the rate of genetic drift compared to an equal, density-independent sex ratio, but a stochastic change towards male bias reduces the N_e /N ratio. Environmental stochasticity amplifies temporal fluctuations in population size and is thus vital to consider in estimation of effective population sizes over longer time periods. Our results on the reduced loss of genetic variation at small densities, particularly in polygamous populations, indicate that density regulation may facilitate adaptive evolution at small population sizes.

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.

Harvest-induced evolution and effective population size

Much has been written about fishery‐induced evolution (FIE) in exploited species, but relatively little attention has been paid to the consequences for one of the most important parameters in evolutionary biology—effective population size (N_e). We use a combination of simulations of Atlantic cod populations experiencing harvest, artificial manipulation of cod life tables, and analytical methods to explore how adding harvest to natural mortality affects N_e, census size (N), and the ratio N_e/N. We show that harvest‐mediated reductions in N_e are due entirely to reductions in recruitment, because increasing adult mortality actually increases the N_e/N ratio. This means that proportional reductions in abundance caused by harvest represent an upper limit to the proportional reductions in N_e, and that in some cases N_e can even increase with increased harvest. This result is a quite general consequence of increased adult mortality and does not depend on harvest selectivity or FIE, although both of these influence the results in a quantitative way. In scenarios that allowed evolution, N_e recovered quickly after harvest ended and remained higher than in the preharvest population for well over a century, which indicates that evolution can help provide a long‐term buffer against loss of genetic variability.

Conservation genetics of a threatened butterfly: comparison of allozymes, RAPDs and microsatellites

Background

Addressing genetic issues in the management of fragmented wild populations of threatened species is one of the most important challenges in conservation biology. Nowadays, a diverse array of molecular methods exists to assess genetic diversity and differentiation of wild populations such as allozymes, dominant markers and co-dominant markers. However it remains worthwhile i) to compare the genetic estimates obtained using those several markers in order to ii) test their relative utility, reliability and relevance and iii) the impact of these results for the design of species-specific conservation measures.

Results

Following the successful isolation of 15 microsatellites loci for the cranberry fritillary butterfly, Boloria aquilonaris, we analyzed the genetic diversity and structure of eight populations located in four different landscapes, at both the regional and the landscape scales. We confront results based on microsatellites to those obtained using allozymes and RAPDs on the same samples. Genetic population analyses using different molecular markers indicate that the B. aquilonaris populations are characterized by a weak genetic variation, likely due to low effective population size and low dispersal at the regional scale. This results in inbreeding in some populations, which may have detrimental consequences on their long term viability. However, gene flow within landscape is limited but not inexistent, with some long range movements resulting in low or no isolation by distance. Spatial structuring was detected among the most isolated populations.

Conclusions

The use of allozymes and RAPD are of very limited value to determine population structuring at small spatial (i.e. landscape) scales, microsatellites giving much higher estimate resolution. The use of RAPD data is also limited for evidencing inbreeding. However, coarse-grain spatial structure (i.e. regional scale), and gene flow estimates based on RAPD and microsatellites data gave congruent results. At a time with increasing development of new molecular methods and markers, dominant markers may still be worthwhile to consider in organisms for which no genomic information is available, and for which limited resources are available.

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.

Using museum specimens to assess historical distribution and genetic diversity in an endangered butterfly

The Miami blue butterfly,Cyclargus thomasi bethunebakeri, is a state-endangered taxon in Florida and a candidate for federal listing. This once common butterfly saw a dramatic decline in population number and abundance in the 1970s and 1980s, but significant collections of individuals prior to this decline are deposited in natural history museums. Using museum specimens, we quantified the genetic diversity in a historical population present in Key Largo, Florida in 1940, 1960, and 1980. Genetic diversity was consistently high within this historical population, but diversity was observed to decrease over the decades sampled. A comparison of historical diversity from the Key Largo population with the extant populations on Bahia Honda State Park (BHSP) and Key West National Wildlife Refuge (KWNWR) revealed differences in allelic frequencies, but only minor differences in the overall number of alleles. The historical distribution of butterflies throughout the Florida Keys further suggests a metapopulation structure. This structure involved partially-isolated populations ofC. t. bethunebakerithat were loosely connected via gene flow and that underwent localized extinction and colonization events along the chain of suitable habitat in the Florida Keys. It appears that a “mini-metapopulation” currently exists on BHSP and KWNWR; structures that are similar to the historical metapopulation structure and distribution of populations on a larger scale. Knowledge of historical distribution helps to plan future reintroduction events with captive-bred butterflies. Additional populations of butterflies may represent undiscovered genetic diversity that, if appropriate, may be further incorporated into captive-breeding efforts.

Effective population size of natural populations of Drosophila buzzatii, with a comparative evaluation of nine methods of estimation

Allozyme and microsatellite data from numerous populations of Drosophila buzzatii have been used (i) to determine to what degree N(e) varies among generations within populations, and among populations, and (ii) to evaluate the congruence of four temporal and five single-sample estimators of N(e) . Effective size of different populations varied over two orders of magnitude, most populations are not temporally stable in genetic composition, and N(e) showed large variation over generations in some populations. Short-term N(e) estimates from the temporal methods were highly correlated, but the smallest estimates were the most precise for all four methods, and the most consistent across methods. Except for one population, N(e) estimates were lower when assuming gene flow than when assuming populations that were closed. However, attempts to jointly estimate N(e) and immigration rate were of little value because the source of migrants was unknown. Correlations among the estimates from the single-sample methods generally were not significant although, as for the temporal methods, estimates were most consistent when they were small. These single-sample estimates of current N(e) are generally smaller than the short-term temporal estimates. Nevertheless, population genetic variation is not being depleted, presumably because of past or ongoing migration. A clearer picture of current and short-term effective population sizes will only follow with better knowledge of migration rates between populations. Different methods are not necessarily estimating the same N(e) , they are subject to different bias, and the biology, demography and history of the population(s) may affect different estimators differently.

Changes in the genetic structure of Aedes aegypti (Diptera: Culicidae) populations in Queensland, Australia, across two seasons: implications for potential mosquito releases

Diseases transmitted by mosquitoes could be controlled if vector populations were replaced with strains that have reduced vector competency. Such a strategy is being developed for control of dengue virus which is transmitted by Aedes aegypti (L.) (Diptera: Culicidae). Mosquitoes artificially infected with the bacterium, Wolbachia pipientis Hertig, are being assessed as candidates for release at the adult stage with the aim of replacement of the wild population. Wolbachia can reduce the capacity of Ae. aegypti to transmit dengue virus and has potential to be driven through the natural population via a system of cytoplasmic incompatibility. Deployment of benign mosquito strains will be influenced by population size and structure of wild-type Ae. aegypti in proposed release areas, as well as rates of gene flow among populations in the wet and dry tropical seasons. Mosquitoes from northern Queensland were screened with genetic markers to find an optimal locality for release of a benign strain of Ae. aegypti. The inland towns of Chillagoe and Charters Towers and the coastal town of Ingham had mosquito populations that were partly genetically isolated from mosquitoes in other areas across both seasons. These locations may be suitable release sites if it is important for the released strain to be restricted during initial phases of implementation. Smaller genetic differences were also evident among other regions and were consistent over two seasons (wet and dry).

Metapopulation structure and fine-scaled genetic structuring in crop-wild hybrid weed beets

This study explores the microspatial and temporal genetic variation in crop-wild hybrid weed beets that emerged from the seed bank in a cultivated field surveyed over two successive years. We demonstrate the occurrence of demes highly genetically differentiated, kin-structured, characterized by moderate effective population sizes, differing in propensity for selfing, and arising from nonrandom genetic subsets of the seed bank. Only one deme identified in the first survey year significantly contributed to the weed beets that emerged in the second year. Spatial structuring appears to be primarily due to gravity seed dispersal and limited pollen flow among weed beet demes. Within each genetic cluster identified by Bayesian assignments and multivariate analyses, F(IS) estimates and level of biparental inbreeding--revealed by progeny analyses--dropped to non-significant values. This suggests that random mating occurs at the scale of genetically distinct demes over a very short scale. Our results highlight the need to carefully depict genetic discontinuities in weed species, when attempting to describe their local genetic neighborhoods within which genetic drift and selective processes occur.

Genetic effective size, N(e), tracks density in a small freshwater cyprinid, Pecos bluntnose shiner (Notropis simus pecosensis)

Genetic monitoring tracks changes in measures of diversity including allelic richness, heterozygosity and genetic effective size over time, and has emerged as an important tool for understanding evolutionary consequences of population management. One proposed application of genetic monitoring has been to estimate abundance and its trajectory through time. Here, genetic monitoring was conducted across five consecutive year for the Pecos bluntnose shiner, a federally threatened minnow. Temporal changes in allele frequencies at seven microsatellite DNA loci were used to estimate variance effective size (N(eV)) across adjacent years in the time series. Likewise, effective size was computed using the linkage disequilibrium method (N(eD)) for each sample. Estimates of N(e) were then compared to estimates of adult fish density obtained from traditional demographic monitoring. For Pecos bluntnose shiner, density (catch-per-unit-effort), N(eV) and N(eD) were positively associated across this time series. Results for Pecos bluntnose shiner were compared to a related and ecologically similar species, the Rio Grande silvery minnow. In this species, density and N(eV) were negatively associated, which suggested decoupling of abundance and effective size trajectories. Conversely, density and N(eD) were positively associated. For Rio Grande silvery minnow, discrepancies among estimates of N(e) and their relationships with adult fish density could be related to effects of high variance in reproductive success in the wild and/or effects of supplementation of the wild population with captive-bred and reared fish. The efficacy of N(e) as a predictor of density and abundance may depend on intrinsic population dynamics of the species and how these dynamics are influenced by the landscape features, management protocols and other factors.