Genetic bottlenecks driven by population disconnection

Drivers of genomic diversity and phenotypic development in early phases of domestication in Hermetia illucens

The black soldier fly (BSF), Hermetia illucens, has the ability to efficiently bioremediate organic waste into usable bio-compounds. Understanding the impact of domestication and mass rearing on fitness and production traits is therefore important for sustainable production. This study aimed to assess patterns of genomic diversity and its association to phenotypic development across early generations of mass rearing under two selection strategies: selection for greater larval mass (SEL lines) and no direct artificial selection (NS lines). Genome-wide single nucleotide polymorphism (SNP) data were generated using 2bRAD sequencing, while phenotypic traits relating to production and population fitness were measured. Declining patterns of genomic diversity were observed across three generations of captive breeding, with the lowest diversity recorded for the F3 generation of both selection lines, most likely due to founder effects. The SEL cohort displayed statistically significantly greater larval weight com the NS lines with pronounced genetic and phenotypic directional changes across generations. Furthermore, lower genetic and phenotypic diversity, particularly for fitness traits, were evident for SEL lines, illustrating the trade-off between selecting for mass and the resulting decline in population fitness. SNP-based heritability was significant for growth, but was low or non-significant for fitness traits. Genotype-phenotype correlations were observed for traits, but individual locus effect sizes where small and very few of these loci demonstrated a signature for selection. Pronounced genetic drift, due to small effective population sizes, is likely overshadowing the impacts of selection on genomic diversity and consequently phenotypic development. The results hold particular relevance for genetic management and selective breeding for BSF in future.

Ancient DNA indicates a century of overhunting did not reduce genetic diversity in Pacific Walruses (Odobenus rosmarus divergens)

Pacific Walruses (Odobenus rosmarus divergens [Illiger 1815]) are gregarious marine mammals considered to be sentinels of the Arctic because of their dependence on sea ice for feeding, molting, and parturition. Like many other marine mammal species, their population sizes were decimated by historical overhunting in the nineteenth and twentieth centuries. Although they have since been protected from nearly all commercial hunting pressure, they now face rapidly accelerating habitat loss as global warming reduces the extent of summer sea ice in the Arctic. To investigate how genetic variation was impacted by overhunting, we obtained mitochondrial DNA sequences from historic Pacific Walrus samples in Alaska that predate the period of overhunting, as well as from extant populations. We found that genetic variation was unchanged over this period, suggesting Pacific Walruses are resilient to genetic attrition in response to reduced population size, and that this may be related to their high vagility and lack of population structure. Although Pacific Walruses will almost certainly continue to decline in number as the planet warms and summer sea ice is further reduced, they may be less susceptible to the ratcheting effects of inbreeding that typically accompany shrinking populations.

A renewed glance at the Palearctic golden eagle: Genetic variation in space and time

Anthropogenic pressures on nature have been causing population declines for centuries. Intensified persecution of apex predators, like the golden eagle, resulted in population bottlenecks during the 19th and 20th centuries. To study population genetics and demographic history of the golden eagle throughout its distribution, we collected museum samples from previously underrepresented regions, such as Russia and Central Asia. We used 12 microsatellite loci and a fragment of the mitochondrial DNA control region to re-evaluate phylogeography of Eurasian golden eagles and study the impacts of the population bottleneck. Our results revealed a north-south genetic gradient, expressed by the difference between Mediterranean and Holarctic lineages, as well as genetically distinct Northern Europe and Central Asia and Caucasus regions. Furthermore, Northern Europe exhibited the lowest, whereas Central Asia and Caucasus had the highest genetic diversity. Although golden eagles maintained relatively high genetic diversity, we detected genetic signatures of the recent bottleneck, including reduced genetic diversity and a decline in the effective female population size around the year 1975. Our study improves the knowledge of the genetic composition of Eurasian golden eagles and highlights the importance of understanding their historical population dynamics in the face of ongoing and future conservation efforts.

Genetic diversity, structure, and effective population size of an endangered, endemic hoary bat, 'ōpe'ape'a, across the Hawaiian Islands

Island bat species are disproportionately at risk of extinction, and Hawai'i's only native terrestrial land mammal, the Hawaiian hoary bat (Lasiurus semotus) locally known as 'ōpe'ape'a, is no exception. To effectively manage this bat species with an archipelago-wide distribution, it is important to determine the population size on each island and connectivity between islands. We used 18 nuclear microsatellite loci and one mitochondrial gene from 339 individuals collected from 1988-2020 to evaluate genetic diversity, population structure and estimate effective population size on the Islands of Hawai'i, Maui, O'ahu, and Kaua'i. Genetic differentiation occurred between Hawai'i and Maui, both of which were differentiated from O'ahu and Kaua'i. The population on Maui presents the greatest per-island genetic diversity, consistent with their hypothesized status as the original founding population. A signature of isolation by distance was detected between islands, with contemporary migration analyses indicating limited gene flow in recent generations, and male-biased sex dispersal within Maui. Historical and long-term estimates of genetic effective population sizes were generally larger than contemporary estimates, although estimates of contemporary genetic effective population size lacked upper bounds in confidence intervals for Hawai'i and Kaua'i. Contemporary genetic effective population sizes were smaller on O'ahu and Maui. We also detected evidence of past bottlenecks on all islands with the exception of Hawai'i. Our study provides population-level estimates for the genetic diversity and geographic structure of 'ōpe'ape'a, that could be used by agencies tasked with wildlife conservation in Hawai'i.

Searching for genetic evidence of demographic decline in an arctic seabird: beware of overlapping generations

Genetic data are useful for detecting sudden population declines in species that are difficult to study in the field. Yet this indirect approach has its own drawbacks, including population structure, mutation patterns, and generation overlap. The ivory gull (Pagophila eburnea), a long-lived Arctic seabird, is currently suffering from rapid alteration of its primary habitat (i.e., sea ice), and dramatic climatic events affecting reproduction and recruitment. However, ivory gulls live in remote areas, and it is difficult to assess the population trend of the species across its distribution. Here we present complementary microsatellite- and SNP-based genetic analyses to test a recent bottleneck genetic signal in ivory gulls over a large portion of their distribution. With attention to the potential effects of population structure, mutation patterns, and sample size, we found no significant signatures of population decline worldwide. At a finer scale, we found a significant bottleneck signal at one location in Canada. These results were compared with predictions from simulations showing how generation time and generation overlap can delay and reduce the bottleneck microsatellite heterozygosity excess signal. The consistency of the results obtained with independent methods strongly indicates that the species shows no genetic evidence of an overall decline in population size. However, drawing conclusions related to the species' population trends will require a better understanding of the effect of age structure in long-lived species. In addition, estimates of the effective global population size of ivory gulls were surprisingly low (~1000 ind.), suggesting that the evolutionary potential of the species is not assured.

Genetic diversity in a long-lived mammal is explained by the past's demographic shadow and current connectivity

Within-species genetic diversity is crucial for the persistence and integrity of populations and ecosystems. Conservation actions require an understanding of factors influencing genetic diversity, especially in the context of global change. Both population size and connectivity are factors greatly influencing genetic diversity; the relative importance of these factors can, however, change through time. Hence, quantifying the degree to which population size or genetic connectivity are shaping genetic diversity, and at which ecological time scale (past or present), is challenging, yet essential for the development of efficient conservation strategies. In this study, we estimated the genetic diversity of 42 colonies of Rhinolophus hipposideros, a long-lived mammal vulnerable to global change, sampling locations spanning its continental northern range. Here, we present an integrative approach that disentangles and quantifies the contribution of different connectivity measures in addition to contemporary colony size and historic bottlenecks in shaping genetic diversity. In our study, the best model explained 64% of the variation in genetic diversity. It included historic bottlenecks, contemporary colony size, connectivity and a negative interaction between the latter two. Contemporary connectivity explained most genetic diversity when considering a 65 km radius around the focal colonies, emphasizing the large geographic scale at which the positive impact of connectivity on genetic diversity is most profound and hence, the minimum scale at which conservation should be planned. Our results highlight that the relative importance of the two main factors shaping genetic diversity varies through time, emphasizing the relevance of disentangling them to ensure appropriate conservation strategies.

Varying genetic imprints of road networks and human density in North American mammal populations

Road networks and human density are major factors contributing to habitat fragmentation and loss, isolation of wildlife populations, and reduced genetic diversity. Terrestrial mammals are particularly sensitive to road networks and encroachment by human populations. However, there are limited assessments of the impacts of road networks and human density on population-specific nuclear genetic diversity, and it remains unclear how these impacts are modulated by life-history traits. Using generalized linear mixed models and microsatellite data from 1444 North American terrestrial mammal populations, we show that taxa with large home range sizes, dense populations, and large body sizes had reduced nuclear genetic diversity with increasing road impacts and human density, but the overall influence of life-history traits was generally weak. Instead, we observed a high degree of genus-specific variation in genetic responses to road impacts and human density. Human density negatively affected allelic diversity or heterozygosity more than road networks (13 vs. 5-7 of 25 assessed genera, respectively); increased road networks and human density also positively affected allelic diversity and heterozygosity in 15 and 6-9 genera, respectively. Large-bodied, human-averse species were generally more negatively impacted than small, urban-adapted species. Genus-specific responses to habitat fragmentation by ongoing road development and human encroachment likely depend on the specific capability to (i) navigate roads as either barriers or movement corridors, and (ii) exploit resource-rich urban environments. The nonuniform genetic response to roads and human density highlights the need to implement efforts to mitigate the risk of vehicular collisions, while also facilitating gene flow between populations of particularly vulnerable taxa.

Conservation genetics of the yellow-bellied toad (Bombina variegata): population structure, genetic diversity and landscape effects in an endangered amphibian

Revealing patterns of genetic diversity and barriers for gene flow are key points for successful conservation in endangered species. Methods based on molecular markers are also often used to delineate conservation units such as evolutionary significant units and management units. Here we combine phylo-geographic analyses (based on mtDNA) with population and landscape genetic analyses (based on microsatellites) for the endangered yellow-bellied toad Bombina variegata over a wide distribution range in Germany. Our analyses show that two genetic clusters are present in the study area, a northern and a southern/central one, but that these clusters are not deeply divergent. The genetic data suggest high fragmentation among toad occurrences and consequently low genetic diversity. Genetic diversity and genetic connectivity showed a negative relationship with road densities and urban areas surrounding toad occurrences, indicating that these landscape features act as barriers to gene flow. To preserve a maximum of genetic diversity, we recommend considering both genetic clusters as management units, and to increase gene flow among toad occurrences with the aim of restoring and protecting functional meta-populations within each of the clusters. Several isolated populations with especially low genetic diversity and signs of inbreeding need particular short-term conservation attention to avoid extinction. We also recommend to allow natural gene flow between both clusters but not to use individuals from one cluster for translocation or reintroduction into the other. Our results underscore the utility of molecular tools for species conservation, highlight outcomes of habitat fragmentation onto the genetic structure of an endangered amphibian and reveal particularly threatened populations in need for urgent conservation efforts.

Harnessing paleo-environmental modeling and genetic data to predict intraspecific genetic structure

Spatially explicit simulations of gene flow within complex landscapes could help forecast the responses of populations to global and anthropological changes. Simulating how past climate change shaped intraspecific genetic variation can provide a validation of models in anticipation of their use to predict future changes. We review simulation models that provide inferences on population genetic structure. Existing simulation models generally integrate complex demographic and genetic processes but are less focused on the landscape dynamics. In contrast to previous approaches integrating detailed demographic and genetic processes and only secondarily landscape dynamics, we present a model based on parsimonious biological mechanisms combining habitat suitability and cellular processes, applicable to complex landscapes. The simulation model takes as input (a) the species dispersal capacities as the main biological parameter, (b) the species habitat suitability, and (c) the landscape structure, modulating dispersal. Our model emphasizes the role of landscape features and their temporal dynamics in generating genetic differentiation among populations within species. We illustrate our model on caribou/reindeer populations sampled across the entire species distribution range in the Northern Hemisphere. We show that simulations over the past 21 kyr predict a population genetic structure that matches empirical data. This approach looking at the impact of historical landscape dynamics on intraspecific structure can be used to forecast population structure under climate change scenarios and evaluate how species range shifts might induce erosion of genetic variation within species.

Large-Scale Hybridisation as an Extinction Threat to the Suweon Treefrog (Hylidae: Dryophytes suweonensis)

Simple Summary

A large number of amphibian species are now endangered, mostly because of human activities. An example is land modification, which may bring species that were previously isolated in contact, and allows them to hybridise. Here, we assessed the presence of hybrid individuals between the endangered Suweon treefrog (Dryophytes suweonensis) and the widespread Japanese treefrog (Dryophytes japonicus). We found hybrids to be relatively widespread and present at all populations where the Suweon treefrog occurred. This is important, as it results in an additional threat to the Suweon treefrog.

Abstract

Amphibians are in the midst of a sixth mass extinction, and human activities play a major role in pushing species towards extinction. Landscape anthropisation has impacts that indirectly threaten species, in addition to the obvious destruction of natural habitats. For instance, land modification may bring human-commensal species in contact with sister-clades from which they were previously isolated. The species in these new contact zones are then able to hybridise to the point of reaching lineage fusion, through which the gene pool of the two species merges and one of the parental lineages becomes extirpated. Here, we documented the patterns of hybridisation between the spatially restricted D. suweonensis and the widespread D. japonicus. On the basis of the analysis of Cytochrome c oxidase subunit I mitochondrial DNA sequences (404 individuals from 35 sites) and six polymorphic microsatellites (381 individuals from 34 sites), we revealed a generalised, bi-directional, and geographically widespread hybridisation between the two species. Evidence of fertile back-crosses is provided by relatively high numbers of individuals in cyto-nuclear disequilibrium, as well as the presence of hybrid individuals further south than the species distribution limit, determined on the basis of call properties. Hybridisation is an additional threat to the endangered D. suweonensis.

Spatial genetic structure and diversity of natural populations of Aesculus hippocastanum L. in Greece

Horse-chestnut (Aesculus hippocastanum L.) is an endemic and relict species from the Mediterranean biodiversity hotspot and a popular ornamental tree. Knowledge about the evolutionary history of this species remains scarce. Here, we ask what historical and ecological factors shaped the pattern of genetic diversity and differentiation of this species. We genotyped 717 individuals from nine natural populations using microsatellite markers. The influence of distance, topography and habitat variables on spatial genetic structure was tested within the approaches of isolation-by-distance and isolation-by-ecology. Species niche modeling was used to project the species theoretical range through time and space. The species showed high genetic diversity and moderate differentiation for which topography, progressive range contraction through the species’ history and long-term persistence in stable climatic refugia are likely responsible. A strong geographic component was revealed among five genetic clusters that are connected with very limited gene flow. The environmental variables were a significant factor in the spatial genetic structure. Modeling results indicated that future reduction of the species range may affect its survival. The possible impact of climate changes and high need of in situ conservation are discussed.

Genetic decline, restoration and rescue of an isolated ungulate population

Isolation of small populations is expected to reduce fitness through inbreeding and loss of genetic variation, impeding population growth and compromising population persistence. Species with long generation time are the least likely to be rescued by evolution alone. Management interventions that maintain or restore genetic variation to assure population viability are consequently of significant importance. We investigated, over 27 years, the genetic and demographic consequences of a demographic bottleneck followed by artificial supplementation in an isolated population of bighorn sheep (Ovis canadensis). Based on a long-term pedigree and individual monitoring, we documented the genetic decline, restoration and rescue of the population. Microsatellite analyses revealed that the demographic bottleneck reduced expected heterozygosity and allelic diversity by 6.2% and 11.3%, respectively, over two generations. Following supplementation, first-generation admixed lambs were 6.4% heavier at weaning and had 28.3% higher survival to 1 year compared to lambs of endemic ancestry. Expected heterozygosity and allelic diversity increased by 4.6% and 14.3% after two generations through new alleles contributed by translocated individuals. We found no evidence for outbreeding depression and did not see immediate evidence of swamping of local genes. Rapid intervention following the demographic bottleneck allowed the genetic restoration and rescue of this bighorn sheep population, likely preventing further losses at both the genetic and demographic levels. Our results provide further empirical evidence that translocation can be used to reduce inbreeding depression in nature and has the potential to mitigate the effect of human-driven environmental changes on wild populations.

Challenges of next-generation sequencing in conservation management: Insights from long-term monitoring of corridor effects on the genetic diversity of mouse lemurs in a fragmented landscape

Long-term genetic monitoring of populations is essential for efforts aimed at preserving genetic diversity of endangered species. Here, we employ a framework of long-term genetic monitoring to evaluate the effects of fragmentation and the effectiveness of the establishment of corridors in restoring population connectivity and genetic diversity of mouse lemurs Microcebus ganzhorni. To this end, we supplement estimates of neutral genetic diversity with the assessment of adaptive genetic variability of the major histocompatibility complex (MHC). In addition, we address the challenges of long-term genetic monitoring of functional diversity by comparing the genotyping performance and estimates of MHC variability generated by single-stranded conformation polymorphism (SSCP)/Sanger sequencing with those obtained by high-throughput sequencing (next-generation sequencing [NGS], Illumina), an issue that is particularly relevant when previous work serves as a baseline for planning management strategies that aim to ensure the viability of a population. We report that SSCP greatly underestimates individual diversity and that discrepancies in estimates of MHC diversity attributable to the comparisons of traditional and NGS genotyping techniques can influence the conclusions drawn from conservation management scenarios. Evidence of migration among fragments in Mandena suggests that mouse lemurs are robust to the process of fragmentation and that the effect of corridors is masked by ongoing gene flow. Nonetheless, results based on a larger number of shared private alleles at neutral loci between fragment pairs found after the establishment of corridors in Mandena suggest that gene flow is augmented as a result of enhanced connectivity. Our data point out that despite low effective population size, M. ganzhorni maintains high individual heterozygosity at neutral loci and at MHC II DRB gene and that selection plays a predominant role in maintaining MHC diversity. These findings highlight the importance of long-term genetic monitoring in order to disentangle between the processes of drift and selection maintaining adaptive genetic diversity in small populations.

Pleistocene population expansions of shade-tolerant trees indicate fragmentation of the African rainforest during the Ice Ages

The fossil record in tropical Africa suggests that dry conditions during the Ice Ages caused expansion of savannahs and contraction of the rainforest. Forest refugia have been proposed to be located in areas of Central Africa that currently harbour high rates of endemic species. However, to what extent the forest was fragmented remains unknown. Nuclear microsatellites and plastid sequences of 732 trees of two species occurring in the same habitat-mature lowland evergreen rainforests-but with remarkably different dispersal capacities-animal versus gravity-were analysed. Geographical information system tools revealed intraspecific lineages partially congruent across the two species, suggesting common past barriers to gene flow in Central Africa. According to approximate Bayesian computation, the intraspecific genetic clusters diverged during the Pleistocene (less than 2 Ma), so that intraspecific differentiation is the appropriate scale to test the aridification effect of the Ice Ages on tree populations. Demographic tests revealed clear genetic signals of population expansion in both taxa, possibly following bottleneck events after forest fragmentation, with stronger evidence of expansion after the Penultimate rather than after the Last Glacial Maximum. The differential dispersal capacity may have modulated the particular response of each species to climate change, as revealed by the stronger evidence of expansion found in the animal-dispersed species than in the gravity-dispersed one.

Strong fluctuations in aboveground population size do not limit genetic diversity in populations of an endangered biennial species

Assessing genetic diversity within populations of rare species and understanding its determinants are crucial for effective species protection. While a lot is known about the relationships between genetic diversity, fitness, and current population size, very few studies explored the effects of past population size. Knowledge of past population size may, however, improve our ability to predict future population fates. We studied Gentianella praecox subsp. bohemica, a biennial species with extensive seed bank. We tested the effect of current, past minimal and maximal population size, and harmonic mean of population sizes within the last 15 years on genetic diversity and fitness. Maximum population size over the last 15 years was the best predictor of expected heterozygosity of the populations and was significantly related to current population size and management. Plant fitness was significantly related to current as well as maximum population size and expected heterozygosity. The results suggested that information on past population size may improve our understanding of contemporary genetic diversity across populations. They demonstrated that despite the strong fluctuations in population size, large reductions in population size do not result in immediate loss of genetic diversity and reduction of fitness within the populations. This is likely due to the seed bank of the species serving as reservoir of the genetic diversity of the populations. From a conservation point of view, this suggests that the restoration of small populations of short-lived species with permanent seed bank is possible as these populations may still be genetically diverse.

Loss of connectivity among island-dwelling Peary caribou following sea ice decline

Global warming threatens to reduce population connectivity for terrestrial wildlife through significant and rapid changes to sea ice. Using genetic fingerprinting, we contrasted extant connectivity in island-dwelling Peary caribou in northern Canada with continental-migratory caribou. We next examined if sea-ice contractions in the last decades modulated population connectivity and explored the possible impact of future climate change on long-term connectivity among island caribou. We found a strong correlation between genetic and geodesic distances for both continental and Peary caribou, even after accounting for the possible effect of sea surface. Sea ice has thus been an effective corridor for Peary caribou, promoting inter-island connectivity and population mixing. Using a time series of remote sensing sea-ice data, we show that landscape resistance in the Canadian Arctic Archipelago has increased by approximately 15% since 1979 and may further increase by 20-77% by 2086 under a high-emission scenario (RCP8.5). Under the persistent increase in greenhouse gas concentrations, reduced connectivity may isolate island-dwelling caribou with potentially significant consequences for population viability.

Genetic diversity of calcareous grassland plant species depends on historical landscape configuration

Background

Habitat fragmentation is considered to be a main reason for decreasing genetic diversity of plant species. However, the results of many fragmentation studies are inconsistent. This may be due to the influence of habitat conditions, having an indirect effect on genetic variation via reproduction. Consequently we took a comparative approach to analyse the impact of habitat fragmentation and habitat conditions on the genetic diversity of calcareous grassland species in this study. We selected five typical grassland species (Primula veris, Dianthus carthusianorum, Medicago falcata, Polygala comosa and Salvia pratensis) occurring in 18 fragments of calcareous grasslands in south eastern Germany. We sampled 1286 individuals in 87 populations and analysed genetic diversity using amplified fragment length polymorphisms. Additionally, we collected data concerning habitat fragmentation (historical and present landscape structure) and habitat conditions (vegetation structure, soil conditions) of the selected study sites. The whole data set was analysed using Bayesian multiple regressions.

Results

Our investigation indicated a habitat loss of nearly 80% and increasing isolation between grasslands since 1830. Bayesian analysis revealed a significant impact of the historical landscape structure, whereas habitat conditions played no important role for the present-day genetic variation of the studied plant species.

Conclusions

Our study indicates that the historical landscape structure may be more important for genetic diversity than present habitat conditions. Populations persisting in abandoned grassland fragments may contribute significantly to the species’ variability even under deteriorating habitat conditions. Therefore, these populations should be included in approaches to preserve the genetic variation of calcareous grassland species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12898-017-0129-9) contains supplementary material, which is available to authorized users.

A longitudinal genetic survey identifies temporal shifts in the population structure of Dutch house sparrows

Dutch house sparrow (Passer domesticus) densities dropped by nearly 50% since the early 1980s, and similar collapses in population sizes have been reported across Europe. Whether, and to what extent, such relatively recent demographic changes are accompanied by concomitant shifts in the genetic population structure of this species needs further investigation. Therefore, we here explore temporal shifts in genetic diversity, genetic structure and effective sizes of seven Dutch house sparrow populations. To allow the most powerful statistical inference, historical populations were resampled at identical locations and each individual bird was genotyped using nine polymorphic microsatellites. Although the demographic history was not reflected by a reduction in genetic diversity, levels of genetic differentiation increased over time, and the original, panmictic population (inferred from the museum samples) diverged into two distinct genetic clusters. Reductions in census size were supported by a substantial reduction in effective population size, although to a smaller extent. As most studies of contemporary house sparrow populations have been unable to identify genetic signatures of recent population declines, results of this study underpin the importance of longitudinal genetic surveys to unravel cryptic genetic patterns.

Considering the Influence of Nonadaptive Evolution on Primate Color Vision

Color vision in primates is variable across species, and it represents a rare trait in which the genetic mechanisms underlying phenotypic variation are fairly well-understood. Research on primate color vision has largely focused on adaptive explanations for observed variation, but it remains unclear why some species have trichromatic or polymorphic color vision while others are red-green color blind. Lemurs, in particular, are highly variable. While some species are polymorphic, many closely-related species are strictly dichromatic. We provide the first characterization of color vision in a wild population of red-bellied lemurs (Eulemur rubriventer, Ranomafana National Park, Madagascar) with a sample size (87 individuals; N_{X chromosomes} = 134) large enough to detect even rare variants (0.95 probability of detection at ≥ 3% frequency). By sequencing exon 5 of the X-linked opsin gene we identified opsin spectral sensitivity based on known diagnostic sites and found this population to be dichromatic and monomorphic for a long wavelength allele. Apparent fixation of this long allele is in contrast to previously published accounts of Eulemur species, which exhibit either polymorphic color vision or only the medium wavelength opsin. This unexpected result may represent loss of color vision variation, which could occur through selective processes and/or genetic drift (e.g., genetic bottleneck). To indirectly assess the latter scenario, we genotyped 55 adult red-bellied lemurs at seven variable microsatellite loci and used heterozygosity excess and M-ratio tests to assess if this population may have experienced a recent genetic bottleneck. Results of heterozygosity excess but not M-ratio tests suggest a bottleneck might have occurred in this red-bellied lemur population. Therefore, while selection may also play a role, the unique color vision observed in this population might have been influenced by a recent genetic bottleneck. These results emphasize the need to consider adaptive and nonadaptive mechanisms of color vision evolution in primates.

On the importance of being structured: instantaneous coalescence rates and human evolution--lessons for ancestral population size inference?

Most species are structured and influenced by processes that either increased or reduced gene flow between populations. However, most population genetic inference methods assume panmixia and reconstruct a history characterized by population size changes. This is potentially problematic as population structure can generate spurious signals of population size change through time. Moreover, when the model assumed for demographic inference is misspecified, genomic data will likely increase the precision of misleading if not meaningless parameters. For instance, if data were generated under an n-island model (characterized by the number of islands and migrants exchanged) inference based on a model of population size change would produce precise estimates of a bottleneck that would be meaningless. In addition, archaeological or climatic events around the bottleneck's timing might provide a reasonable but potentially misleading scenario. In a context of model uncertainty (panmixia versus structure) genomic data may thus not necessarily lead to improved statistical inference. We consider two haploid genomes and develop a theory that explains why any demographic model with structure will necessarily be interpreted as a series of changes in population size by inference methods ignoring structure. We formalize a parameter, the inverse instantaneous coalescence rate, and show that it is equivalent to a population size only in panmictic models, and is mostly misleading for structured models. We argue that this issue affects all population genetics methods ignoring population structure which may thus infer population size changes that never took place. We apply our approach to human genomic data.

Demographic inference using genetic data from a single individual: Separating population size variation from population structure

The rapid development of sequencing technologies represents new opportunities for population genetics research. It is expected that genomic data will increase our ability to reconstruct the history of populations. While this increase in genetic information will likely help biologists and anthropologists to reconstruct the demographic history of populations, it also represents new challenges. Recent work has shown that structured populations generate signals of population size change. As a consequence it is often difficult to determine whether demographic events such as expansions or contractions (bottlenecks) inferred from genetic data are real or due to the fact that populations are structured in nature. Given that few inferential methods allow us to account for that structure, and that genomic data will necessarily increase the precision of parameter estimates, it is important to develop new approaches. In the present study we analyze two demographic models. The first is a model of instantaneous population size change whereas the second is the classical symmetric island model. We (i) re-derive the distribution of coalescence times under the two models for a sample of size two, (ii) use a maximum likelihood approach to estimate the parameters of these models (iii) validate this estimation procedure under a wide array of parameter combinations, (iv) implement and validate a model rejection procedure by using a Kolmogorov-Smirnov test, and a model choice procedure based on the AIC, and (v) derive the explicit distribution for the number of differences between two non-recombining sequences. Altogether we show that it is possible to estimate parameters under several models and perform efficient model choice using genetic data from a single diploid individual.

Ancient and contemporary DNA reveal a pre-human decline but no population bottleneck associated with recent human persecution in the kea (Nestor notabilis)

The impact of population bottlenecks is an important factor to consider when assessing species survival. Population declines can considerably limit the evolutionary potential of species and make them more susceptible to stochastic events. New Zealand has a well documented history of decline of endemic avifauna related to human colonization. Here, we investigate the genetic effects of a recent population decline in the endangered kea (Nestor notabilis). Kea have undergone a long-lasting persecution between the late 1800s to 1970s where an estimated 150,000 kea were culled under a governmental bounty scheme. Kea now number 1,000–5,000 individuals in the wild and it is likely that the recent population decline may have reduced the genetic diversity of the species. Comparison of contemporary (n = 410), historical (n = 15) and fossil samples (n = 4) showed a loss of mitochondrial diversity since the end of the last glaciation (Otiran Glacial) but no loss of overall genetic diversity associated with the cull. Microsatellite data indicated a recent bottleneck for only one population and a range-wide decline in N_e dating back some 300 – 6,000 years ago, a period predating European arrival in NZ. These results suggest that despite a recent human persecution, kea might have experienced a large population decline before stabilizing in numbers prior to human settlement of New Zealand in response to Holocene changes in habitat distribution. Our study therefore highlights the need to understand the respective effects of climate change and human activities on endangered species dynamics when proposing conservation guidelines.

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.

Effect of rotational shepherding on demographic and genetic connectivity of calcareous grassland plants

Response to habitat fragmentation may not be generalized among species, in particular for plant communities with a variety of dispersal traits. Calcareous grasslands are one of the most species-rich habitats in Central Europe, but abandonment of traditional management has caused a dramatic decline of calcareous grassland species. In the Southern Franconian Alb in Germany, reintroduction of rotational shepherding in previously abandoned grasslands has restored species diversity, and it has been suggested that sheep support seed dispersal among grasslands. We tested the effect of rotational shepherding on demographic and genetic connectivity of calcareous grassland specialist plants and whether the response of plant populations to shepherding was limited to species dispersed by animals (zoochory). Specifically, we tested competing dispersal models and source and focal patch properties to explain landscape connectivity with patch-occupancy data of 31 species. We fitted the same connectivity models to patch occupancy and nuclear microsatellite data for the herb Dianthus carthusianorum (Carthusian pink). For 27 species, patch connectivity was explained by dispersal by rotational shepherding regardless of adaptations to zoochory, whereas population size (16% species) and patch area (0% species) of source patches were not important predictors of patch occupancy in most species. [Correction made after online publication, February 25, 2014: Population size and patch area percentages were mistakenly inverted, and have now been fixed.] Microsite diversity of focal patches significantly increased the model variance explained by patch occupancy in 90% of the species. For D. carthusianorum, patch connectivity through rotational shepherding explained both patch occupancy and population genetic diversity. Our results suggest shepherding provides dispersal for multiple plant species regardless of their dispersal adaptations and thus offers a useful approach to restore plant diversity in fragmented calcareous grasslands.

Stronger transferability but lower variability in transcriptomic- than in anonymous microsatellites: evidence from Hylid frogs

A simple way to quickly optimize microsatellites in nonmodel organisms is to reuse loci available in closely related taxa; however, this approach can be limited by the stochastic and low cross-amplification success experienced in some groups (e.g. amphibians). An efficient alternative is to develop loci from transcriptome sequences. Transcriptomic microsatellites have been found to vary in their levels of cross-species amplification and variability, but this has to date never been tested in amphibians. Here, we compare the patterns of cross-amplification and levels of polymorphism of 18 published anonymous microsatellites isolated from genomic DNA vs. 17 loci derived from a transcriptome, across nine species of tree frogs (Hyla arborea and Hyla cinerea group). We established a clear negative relationship between divergence time and amplification success, which was much steeper for anonymous than transcriptomic markers, with half-lives (time at which 50% of the markers still amplify) of 1.1 and 37 My, respectively. Transcriptomic markers are significantly less polymorphic than anonymous loci, but remain variable across diverged taxa. We conclude that the exploitation of amphibian transcriptomes for developing microsatellites seems an optimal approach for multispecies surveys (e.g. analyses of hybrid zones, comparative linkage mapping), whereas anonymous microsatellites may be more informative for fine-scale analyses of intraspecific variation. Moreover, our results confirm the pattern that microsatellite cross-amplification is greatly variable among amphibians and should be assessed independently within target lineages. Finally, we provide a bank of microsatellites for Palaearctic tree frogs (so far only available for H. arborea), which will be useful for conservation and evolutionary studies in this radiation.

Combining genetic and demographic data for prioritizing conservation actions: insights from a threatened fish species

Prioritizing and making efficient conservation plans for threatened populations requires information at both evolutionary and ecological timescales. Nevertheless, few studies integrate multidisciplinary approaches, mainly because of the difficulty for conservationists to assess simultaneously the evolutionary and ecological status of populations. Here, we sought to demonstrate how combining genetic and demographic analyses allows prioritizing and initiating conservation plans. To do so, we combined snapshot microsatellite data and a 30-year-long demographic survey on a threatened freshwater fish species (Parachondrostoma toxostoma) at the river basin scale. Our results revealed low levels of genetic diversity and weak effective population sizes (<63 individuals) in all populations. We further detected severe bottlenecks dating back to the last centuries (200–800 years ago), which may explain the differentiation of certain populations. The demographic survey revealed a general decrease in the spatial distribution and abundance of P. toxostoma over the last three decades. We conclude that demo-genetic approaches are essential for (1) identifying populations for which both evolutionary and ecological extinction risks are high; and (2) proposing conservation plans targeted toward these at risk populations, and accounting for the evolutionary history of populations. We suggest that demo-genetic approaches should be the norm in conservation practices.

We combined genetic and demographic data from a threatened freshwater fish species (Parachondrostoma toxostoma) at the river basin scale for conservation purposes. Genetic diversity and effective population sizes are very low, probably due to the strong genetic bottlenecks detected in this study. The species spatial distribution and abundance also decreased during the last decades.

The demographic history of populations experiencing asymmetric gene flow: combining simulated and empirical data

Population structure can significantly affect genetic-based demographic inferences, generating spurious bottleneck-like signals. Previous studies have typically assumed island or stepping-stone models, which are characterized by symmetric gene flow. However, many organisms are characterized by asymmetric gene flow. Here, we combined simulated and empirical data to test whether asymmetric gene flow affects the inference of past demographic changes. Through the analysis of simulated genetic data with three methods (i.e. bottleneck, M-ratio and msvar), we demonstrated that asymmetric gene flow biases past demographic changes. Most biases were towards spurious signals of expansion, albeit their strength depended on values of effective population size and migration rate. It is noteworthy that the spurious signals of demographic changes also depended on the statistical approach underlying each of the three methods. For one of the three methods, biases induced by asymmetric gene flow were confirmed in an empirical multispecific data set involving four freshwater fish species (Squalius cephalus, Leuciscus burdigalensis, Gobio gobio and Phoxinus phoxinus). However, for the two other methods, strong signals of bottlenecks were detected for all species and across two rivers. This suggests that, although potentially biased by asymmetric gene flow, some of these methods were able to bypass this bias when a bottleneck actually occurred. Our results show that population structure and dispersal patterns have to be considered for proper inference of demographic changes from genetic data.

Peak and persistent excess of genetic diversity following an abrupt migration increase

Genetic diversity is essential for population survival and adaptation to changing environments. Demographic processes (e.g., bottleneck and expansion) and spatial structure (e.g., migration, number, and size of populations) are known to shape the patterns of the genetic diversity of populations. However, the impact of temporal changes in migration on genetic diversity has seldom been considered, although such events might be the norm. Indeed, during the millions of years of a species' lifetime, repeated isolation and reconnection of populations occur. Geological and climatic events alternately isolate and reconnect habitats. We analytically document the dynamics of genetic diversity after an abrupt change in migration given the mutation rate and the number and sizes of the populations. We demonstrate that during transient dynamics, genetic diversity can reach unexpectedly high values that can be maintained over thousands of generations. We discuss the consequences of such processes for the evolution of species based on standing genetic variation and how they can affect the reconstruction of a population's demographic and evolutionary history from genetic data. Our results also provide guidelines for the use of genetic data for the conservation of natural populations.

Historical and contemporary DNA indicate fisher decline and isolation occurred prior to the European settlement of California

Establishing if species contractions were the result of natural phenomena or human induced landscape changes is essential for managing natural populations. Fishers (Martes pennanti) in California occur in two geographically and genetically isolated populations in the northwestern mountains and southern Sierra Nevada. Their isolation is hypothesized to have resulted from a decline in abundance and distribution associated with European settlement in the 1800s. However, there is little evidence to establish that fisher occupied the area between the two extant populations at that time. We analyzed 10 microsatellite loci from 275 contemporary and 21 historical fisher samples (1880-1920) to evaluate the demographic history of fisher in California. We did not find any evidence of a recent (post-European) bottleneck in the northwestern population. In the southern Sierra Nevada, genetic subdivision within the population strongly influenced bottleneck tests. After accounting for genetic subdivision, we found a bottleneck signal only in the northern and central portions of the southern Sierra Nevada, indicating that the southernmost tip of these mountains may have acted as a refugium for fisher during the anthropogenic changes of the late 19(th) and early 20(th) centuries. Using a coalescent-based Bayesian analysis, we detected a 90% decline in effective population size and dated the time of decline to over a thousand years ago. We hypothesize that fisher distribution in California contracted to the two current population areas pre-European settlement, and that portions of the southern Sierra Nevada subsequently experienced another more recent bottleneck post-European settlement.

Within- and among-population impact of genetic erosion on adult fitness-related traits in the European tree frog Hyla arborea

Assessing in wild populations how fitness is impacted by inbreeding and genetic drift is a major goal for conservation biology. An approach to measure the detrimental effects of inbreeding on fitness is to estimate correlations between molecular variation and phenotypic performances within and among populations. Our study investigated the effect of individual multilocus heterozygosity on body size, body condition and reproductive investment of males (that is, chorus attendance) and females (that is, clutch mass and egg size) in both small fragmented and large non-fragmented populations of European tree frog (Hyla arborea). Because adult size and/or condition and reproductive investment are usually related, genetic erosion may have detrimental effects directly on reproductive investment, and also on individual body size and condition that in turn may affect reproductive investment. We confirmed that the reproductive investment was highly size-dependent for both sexes. Larger females invested more in offspring production, and larger males attended the chorus in the pond more often. Our results did not provide evidence for a decline in body size, condition and reproductive effort with decreased multilocus heterozygosity both within and among populations. We showed that the lack of heterozygosity-fitness correlations within populations probably resulted from low inbreeding levels (inferior to ca. 20% full-sib mating rate), even in the small fragmented populations. The detrimental effects of fixation load were either low in adults or hidden by environmental variation among populations. These findings will be useful to design specific management actions to improve population persistence.

Effective population size dynamics and the demographic collapse of Bornean orang-utans

Bornean orang-utans experienced a major demographic decline and local extirpations during the Pleistocene and Holocene due to climate change, the arrival of modern humans, of farmers and recent commercially-driven habitat loss and fragmentation. The recent loss of habitat and its dramatic fragmentation has affected the patterns of genetic variability and differentiation among the remaining populations and increased the extinction risk of the most isolated ones. However, the contribution of recent demographic events to such genetic patterns is still not fully clear. Indeed, it can be difficult to separate the effects of recent anthropogenic fragmentation from the genetic signature of prehistoric demographic events. Here, we investigated the genetic structure and population size dynamics of orang-utans from different sites. Altogether 126 individuals were analyzed and a full-likelihood Bayesian approach was applied. All sites exhibited clear signals of population decline. Population structure is known to generate spurious bottleneck signals and we found that it does indeed contribute to the signals observed. However, population structure alone does not easily explain the observed patterns. The dating of the population decline varied across sites but was always within the 200–2000 years period. This suggests that in some sites at least, orang-utan populations were affected by demographic events that started before the recent anthropogenic effects that occurred in Borneo. These results do not mean that the recent forest exploitation did not leave its genetic mark on orang-utans but suggests that the genetic pool of orang-utans is also impacted by more ancient events. While we cannot identify the main cause for this decline, our results suggests that the decline may be related to the arrival of the first farmers or climatic events, and that more theoretical work is needed to understand how multiple demographic events impact the genome of species and how we can assess their relative contributions.

Microsatellite-based evidences of genetic bottlenecks in the cryptic species "Andrographis paniculata Nees": a potential anticancer agent

Andrographis paniculata (AP) is a medicinal plant species introduced into Malaysia. To address the genetic structure and evolutionary connectedness of the Malaysian AP with the Indian AP, a DNA sequence analysis was conducted based on 24 microsatellite markers. Out of the 24 primer sets, seven novel microsatellite primers were designed and amplified intra-specifically according to the available Indian AP sequences at the National Centre for Biotechnology Information (NCBI), where 17 of them were amplified using the cross-species strategy by employing the primers belonging to Acanthus ilicifolius Linn (Acanthaceae) and Lumnitzera racemosa Wild (Combretaceae). The primers were then applied on the Malaysian AP accessions. Sixteen of the new microsatellite loci were amplified successfully. Analysis of these microsatellite sequences, revealed some significant differences between the Indian and Malaysian AP accessions in terms of the size and type of the repeat motifs. These findings depicted the cryptic feature of this species. Despite identifying several heterozygous alleles no polymorphism was observed in the detected loci of the selected accessions. This situation was in concordance with the presence of "fixed heterozygosity" phenomenon in the mentioned loci. Accordingly, this was fully consistent with the occurrence of the genetic bottleneck and founder effect within Malaysian AP population. Apart from the amplification of new microsatellites in this species, our observations could be in agreement with the risk of genetic depletion and consequently extinction of this precious herb in Malaysia. This issue should be taken into consideration in the future studies.

Genetic erosion in wild populations makes resistance to a pathogen more costly

Populations that have suffered from genetic erosion are expected to exhibit reduced average trait values or decreased variation in adaptive traits when experiencing periodic or emergent stressors such as infectious disease. Genetic erosion may consequentially modify the ability of a potential host population to cope with infectious disease emergence. We experimentally investigate this relationship between genetic variability and host response to exposure to an infectious agent both in terms of susceptibility to infection and indirect parasite-mediated responses that also impact fitness. We hypothesized that the deleterious consequences of exposure to the pathogen (Batrachochytrium dendrobatidis) would be more severe for tadpoles descended from European treefrog (Hyla arborea) populations lacking genetic variability. Although all exposed tadpoles lacked detectable infection, we detected this relationship for some indirect host responses, predominantly in genetically depleted animals, as well as an interaction between genetic variability and pathogen dose on life span during the postmetamorphic period. Lack of infection and a decreased mass and postmetamorphic life span in low genetic diversity tadpoles lead us to conclude that genetic erosion, while not affecting the ability to mount effective resistance strategies, also erodes the capacity to invest in resistance, increased tadpole growth rate, and metamorphosis relatively simultaneously.

Evaluating ring-tailed lemurs (Lemur catta) from southwestern Madagascar for a genetic population bottleneck

In light of historical and recent anthropogenic influences on Malagasy primate populations, in this study ring-tailed lemur (Lemur catta) samples from two sites in southwestern Madagascar, Beza Mahafaly Special Reserve (BMSR) and Tsimanampetsotsa National Park (TNP), were evaluated for the genetic signature of a population bottleneck. A total of 45 individuals (20 from BMSR and 25 from TNP) were genotyped at seven microsatellite loci. Three methods were used to evaluate these populations for evidence of a historical bottleneck: M-ratio, mode-shift, and heterozygosity excess tests. Three mutation models were used for heterozygosity excess tests: the stepwise mutation model (SMM), two-phase model (TPM), and infinite allele model (IAM). M-ratio estimations indicated a potential bottleneck in both populations under some conditions. Although mode-shift tests did not strongly indicate a population bottleneck in the recent historical past when samples from all individuals were included, a female-only analysis indicated a potential bottleneck in TNP. Heterozygosity excess was indicated under two of the three mutation models (IAM and TPM), with TNP showing stronger evidence of heterozygosity excess than BMSR. Taken together, these results suggest that a bottleneck may have occurred among L. catta in southwestern Madagascar in the recent past. Given knowledge of how current major stochastic climatic events and human-induced change can negatively impact extant lemur populations, it is reasonable that comparable events in the historical past could have caused a population bottleneck. This evaluation additionally functions to highlight the continuing environmental and anthropogenic challenges faced by lemurs in southwestern Madagascar.

Heterozygosity-fitness correlations among wild populations of European tree frog (Hyla arborea) detect fixation load

Quantifying the impacts of inbreeding and genetic drift on fitness traits in fragmented populations is becoming a major goal in conservation biology. Such impacts occur at different levels and involve different sets of loci. Genetic drift randomly fixes slightly deleterious alleles leading to different fixation load among populations. By contrast, inbreeding depression arises from highly deleterious alleles in segregation within a population and creates variation among individuals. A popular approach is to measure correlations between molecular variation and phenotypic performances. This approach has been mainly used at the individual level to detect inbreeding depression within populations and sometimes at the population level but without consideration about the genetic processes measured. For the first time, we used in this study a molecular approach considering both the interpopulation and intrapopulation level to discriminate the relative importance of inbreeding depression vs. fixation load in isolated and non-fragmented populations of European tree frog (Hyla arborea), complemented with interpopulational crosses. We demonstrated that the positive correlations observed between genetic heterozygosity and larval performances on merged data were mainly caused by co-variations in genetic diversity and fixation load among populations rather than by inbreeding depression and segregating deleterious alleles within populations. Such a method is highly relevant in a conservation perspective because, depending on how populations lose fitness (inbreeding vs. fixation load), specific management actions may be designed to improve the persistence of populations.

Temporal dynamics of genetic variability in a mountain goat (Oreamnos americanus) population

The association between population dynamics and genetic variability is of fundamental importance for both evolutionary and conservation biology. We combined long-term population monitoring and molecular genetic data from 123 offspring and their parents at 28 microsatellite loci to investigate changes in genetic diversity over 14 cohorts in a small and relatively isolated population of mountain goats (Oreamnos americanus) during a period of demographic increase. Offspring heterozygosity decreased while parental genetic similarity and inbreeding coefficients (F(IS) ) increased over the study period (1995-2008). Immigrants introduced three novel alleles into the population and matings between residents and immigrants produced more heterozygous offspring than local crosses, suggesting that immigration can increase population genetic variability. The population experienced genetic drift over the study period, reflected by a reduced allelic richness over time and an 'isolation-by-time' pattern of genetic structure. The temporal decline of individual genetic diversity despite increasing population size probably resulted from a combination of genetic drift due to small effective population size, inbreeding and insufficient counterbalancing by immigration. This study highlights the importance of long-term genetic monitoring to understand how demographic processes influence temporal changes of genetic diversity in long-lived organisms.