Genetic analysis of a documented population bottleneck: introduced Bennett's wallabies (Macropus rufogriseus rufogriseus) in New Zealand

Few bottlenecks of wild populations are sufficiently well-documented to constitute models for testing theories about the impact of bottlenecks on genetic variation, and subsequent population persistence. Relevant details of the Bennett's wallaby (Macropus rufogriseus rufogriseus) introduction into New Zealand were recorded (founder number, source and approximate bottleneck duration) and suggest this may provide a rare opportunity to examine the efficacy of tests designed to detect recent bottlenecks in wild populations. We first assessed the accuracy of historic accounts of the introduction using genetic diversity detected in mitochondrial DNA (mtDNA) and at five microsatellite loci. Phylogenetic analyses of mtDNA D-loop sequence haplotypes were consistent with the reported origin of the founders as Tasmania, rather than one of the Bass Strait islands in which Bennett's wallabies are also found. Microsatellite allele frequencies from the Tasmanian source population were then used to seed bottleneck simulations encompassing varying sizes and numbers of generations, in order to assess the severity of bottleneck consistent with diversity observed in the New Zealand population. The results suggested that the founder number was unlikely to have been as small as the three animals suggested by the account of the introduction. Nonetheless, the bottleneck was probably severe; in the range of three to five pairs of wallabies for one to three generations. It resulted in significantly reduced levels of allelic diversity and heterozygosity relative to the source population. This bottleneck is only detectable under the infinite allele model (IAM) and not under the stepwise mutation model (SMM) or the two-phase model (TPM), and possible explanations for this are discussed.     

Genetic consequences of a single-founder population bottleneck in Trifolium amoenum (Fabaceae)

We investigated the genetic consequences of a single-founder bottleneck in a population of showy Indian clover (Trifolium amoenum), a species presumed to be extinct until rediscovered near Occidental, California, in 1993. Electrophoretic variation was evaluated in the bottlenecked population and in a larger population (Dillon Beach) discovered during the course of this study, as well as in populations of two closely related species, T. albopurpureum var. dichotomum and T. macraei. We found a surprisingly high amount of polymorphism in the single-founder T. amoenum population from Occidental (15% of loci polymorphic; an average of 1.1 alleles per locus). However, this represents a 53% reduction in number of polymorphic loci and a 20% reduction in average number of alleles per locus compared to three Trifolium populations with putatively similar mating systems (the Dillon Beach T. amoenum population and both populations of T. albopurpureum var. dichotomum). Expanding the genetic base of the Occidental T. amoenum population is a priority due to concerns about loss of evolutionary potential and the possibility of deleterious effects associated with inbreeding. However, using seed from the Dillon Beach T. amoenum population may not be beneficial due to distinct, presumably adaptive differences between plants from the two populations and concerns about outbreeding depression.     

Reliability of genetic bottleneck tests for detecting recent population declines

The identification of population bottlenecks is critical in conservation because populations that have experienced significant reductions in abundance are subject to a variety of genetic and demographic processes that can hasten extinction. Genetic bottleneck tests constitute an appealing and popular approach for determining if a population decline has occurred because they only require sampling at a single point in time, yet reflect demographic history over multiple generations. However, a review of the published literature indicates that, as typically applied, microsatellite-based bottleneck tests often do not detect bottlenecks in vertebrate populations known to have experienced declines. This observation was supported by simulations that revealed that bottleneck tests can have limited statistical power to detect bottlenecks largely as a result of limited sample sizes typically used in published studies. Moreover, commonly assumed values for mutation model parameters do not appear to encompass variation in microsatellite evolution observed in vertebrates and, on average, the proportion of multi-step mutations is underestimated by a factor of approximately two. As a result, bottleneck tests can have a higher probability of 'detecting' bottlenecks in stable populations than expected based on the nominal significance level. We provide recommendations that could add rigor to inferences drawn from future bottleneck tests and highlight new directions for the characterization of demographic history.     

Fragmented habitat compensates for the adverse effects of genetic bottleneck

1. Download : Download high-res image (256KB)
2. Download : Download full-size image

     

Low genetic variability of the koala Phascolarctos cinereus in south-eastern Australia following a severe population bottleneck

Genotyping of koalas at CA-repeat microsatellite loci has revealed significant differences in the levels of allelic diversity ( A ) and expected heterozygosity ( H¯ _E) between populations from north-eastern and south-eastern Australia. In the 10 populations studied, allelic diversity ranged from 8.0 in the Nowendoc population to 1.7 in the Kangaroo Is. population, and values of H¯ _E ranged from 0.831 in the Nowendoc population to 0.331 in the Kangaroo Is. population. Data from pooled populations revealed koalas from the north-eastern region had significantly higher levels of allelic diversity ( A = 11.5 ± 1.4) than those from south-eastern Australia ( A = 5.3 ± 1.0). Furthermore significantly higher heterozygosity levels were found in the north-eastern ( H¯ _E= 0.851) vs. the south-eastern ( H¯ _E= 0.436) regions of Australia. Following a near-extinction bottleneck in the 1920s, mainland Victorian and Kangaroo Is. koalas have been involved in an extensive program of relocations. The source populations of the relocated animals were islands in Westernport Bay, which were founded by very few individuals in the late 1800s and early 1900s. The significantly lower levels of variation between south-eastern Australian populations suggests that human intervention has had a severe effect on levels of genetic diversity in this region, and this may have long-term genetic consequences.     

Genome-wide signatures of population bottlenecks and diversifying selection in European wolves

Genomic resources developed for domesticated species provide powerful tools for studying the evolutionary history of their wild relatives. Here we use 61K single-nucleotide polymorphisms (SNPs) evenly spaced throughout the canine nuclear genome to analyse evolutionary relationships among the three largest European populations of grey wolves in comparison with other populations worldwide, and investigate genome-wide effects of demographic bottlenecks and signatures of selection. European wolves have a discontinuous range, with large and connected populations in Eastern Europe and relatively smaller, isolated populations in Italy and the Iberian Peninsula. Our results suggest a continuous decline in wolf numbers in Europe since the Late Pleistocene, and long-term isolation and bottlenecks in the Italian and Iberian populations following their divergence from the Eastern European population. The Italian and Iberian populations have low genetic variability and high linkage disequilibrium, but relatively few autozygous segments across the genome. This last characteristic clearly distinguishes them from populations that underwent recent drastic demographic declines or founder events, and implies long-term bottlenecks in these two populations. Although genetic drift due to spatial isolation and bottlenecks seems to be a major evolutionary force diversifying the European populations, we detected 35 loci that are putatively under diversifying selection. Two of these loci flank the canine platelet-derived growth factor gene, which affects bone growth and may influence differences in body size between wolf populations. This study demonstrates the power of population genomics for identifying genetic signals of demographic bottlenecks and detecting signatures of directional selection in bottlenecked populations, despite their low background variability.     

Genetic diversity and bottleneck studies in the Marwari horse breed

Genetic diversity within the Marwari breed of horses was evaluated using 26 different microsatellite pairs with 48 DNA samples from unrelated horses. This molecular characterisation was undertaken to evaluate the problem of genetic bottlenecks also, if any, in this breed. The estimated mean (± s.e.) allelic diversity was 5.9 (± 2.24), with a total of 133 alleles. A high level of genetic variability within this breed was observed in terms of high values of mean (±s.e.) effective number of alleles (3.3 ± 1.27), observed heterozygosity (0.5306 ± 0.22), expected Levene’s heterozygosity (0.6612 ± 0.15), expected Nei’s heterozygosity (0.6535 ± 0.14), and polymorphism information content (0.6120 ± 0.03). Low values of Wright’s fixation index, F_IS (0.2433 ± 0.05) indicated low levels of inbreeding. This basic study indicated the existence of substantial genetic diversity in the Marwari horse population. No significant genotypic linkage disequilibrium was detected across the population, suggesting no evidence of linkage between loci. A normal ‘L’ shaped distribution of mode-shift test, non-significant heterozygote excess on the basis of different models, as revealed from Sign, Standardized differences and Wilcoxon sign rank tests as well as non-significantM ratio value suggested that there was no recent bottleneck in the existing Marwari breed population, which is important information for equine breeders. This study also revealed that the Marwari breed can be differentiated from some other exotic breeds of horses on the basis of three microsatellite primers.     

Past population history versus recent population decline – founder effects in island species and their genetic signatures

Island populations are mostly characterized by low genetic diversity if compared with their mainland conspecifics. This is often explained as a consequence of founder effects in the wake of island colonization and concomitant bottlenecks. In a recent contribution, Stuessy et al. (Journal of Biogeography, 2012, 39, 1565–1566) point out that the genetic imprint of past founder effects is no longer visible today, as most island colonizations occurred millions of generations ago. The authors argue that low genetic diversity detectable today is mainly caused by recent environmental factors such as anthropogenic habitat destruction. This scenario should be complemented by the influence of long‐term isolation and small habitat size, which often lead to strong population fluctuations and repeated bottlenecks. In consequence, inbreeding and genetic drift, coupled with the potential effects of purging in small populations, may also result in genetic diversity remaining low for a long time after colonization.     

Genetic diversity, population structure, effective population size and demographic history of the Finnish wolf population

The Finnish wolf population (Canis lupus) was sampled during three different periods (1996-1998, 1999-2001 and 2002-2004), and 118 individuals were genotyped with 10 microsatellite markers. Large genetic variation was found in the population despite a recent demographic bottleneck. No spatial population subdivision was found even though a significant negative relationship between genetic relatedness and geographic distance suggested isolation by distance. Very few individuals did not belong to the local wolf population as determined by assignment analyses, suggesting a low level of immigration in the population. We used the temporal approach and several statistical methods to estimate the variance effective size of the population. All methods gave similar estimates of effective population size, approximately 40 wolves. These estimates were slightly larger than the estimated census size of breeding individuals. A Bayesian model based on Markov chain Monte Carlo simulations indicated strong evidence for a long-term population decline. These results suggest that the contemporary wolf population size is roughly 8% of its historical size, and that the population decline dates back to late 19th century or early 20th century. Despite an increase of over 50% in the census size of the population during the whole study period, there was only weak evidence that the effective population size during the last period was higher than during the first. This may be caused by increased inbreeding, diminished dispersal within the population, and decreased immigration to the population during the last study period.     

Comparing population trend estimates of migratory birds from breeding censuses and capture data at a spring migration bottleneck

Europe has a well‐established network of breeding bird monitoring that is used to produce supranational indices of population trends for many species. However, a comparison of breeding bird censuses with other methods may be beneficial to confirm the validity of such indices. The aim of this study was to assess the value of standardized capture data of migratory birds at migration bottlenecks as an indicator of the effective breeding populations. One limitation to this method is that several populations are co‐occurring at these bottlenecks and their catchment areas need to be clearly identified to allow extrapolation of population indices. Here, we used standardized trends in capture numbers of 30 species on the island of Ponza, a migration bottleneck in the central Mediterranean, and compared them to population trends estimated in the putative catchment breeding areas between 2005 and 2016. The catchment areas were identified through the analysis of ring recoveries during the breeding season of birds passing through Ponza. Our results show an agreement between the population trends observed on Ponza and those in the breeding areas in 15 out of 30 species. The correlations were strongest in species with a more robust definition of the catchment areas, that is, species with more than 10 recoveries, and for which the recoveries were most likely of breeding birds. The main reason for disagreement between the two indices in the remaining species might be related to different intensity of sampling in different areas. This issue can be solved by further developing monitoring projects in underrepresented countries, as well as by intensifying monitoring through ringing, both in the breeding grounds and at migration bottlenecks. These results show that spring migration monitoring at bottlenecks has the potential to provide a valuable complement and an independent control of breeding bird surveys, allowing raising early warnings of population declines and contributing to their conservation.     

Genetics of the Shimokita macaque population suggest an ancient bottleneck

The macaque population of the Shimokita Peninsula represents the northernmost distribution of this species and is isolated from other populations in the Tohoku region of Japan. A previous protein-based study revealed a high level of genetic variability in this population and considerable differentiation from other populations. In order to reassess the genetic features of the Shimokita macaques, we examined 11 autosomal microsatellite loci and three Y chromosomal microsatellite loci. We observed considerable differentiation from other Japanese populations of macaques, but in contrast to the previous results, we observed significantly lower genetic variability in this population. There was a weak indication of a population bottleneck, suggesting a decay over time from an excess of heterozygotes that might be expected in the initial stages of a bottleneck. This may indicate that an ancient bottleneck occurred during the warm period after the last glacial period rather than a recent bottleneck due to hunting in modern times. The frequencies of private alleles were exceptionally high in the Shimokita population, suggesting that the difference in variability as determined in various studies was due to accidental sampling of marker loci with low power to resolve genetic variations in the protein-based studies. The assessments of interpopulation differentiation as determined using autosomal and Y chromosomal markers were highly correlated, and using both types of markers the Shimokita population was found to be the most differentiated of the study populations, probably due to infrequent gene flow with surrounding populations.     

Genetic diversity and population structure of Tasmanian devils, the largest marsupial carnivore

Genetic diversity and population structure were investigated across the core range of Tasmanian devils (Sarcophilus laniarius; Dasyuridae), a wide-ranging marsupial carnivore restricted to the island of Tasmania. Heterozygosity (0.386-0.467) and allelic diversity (2.7-3.3) were low in all subpopulations and allelic size ranges were small and almost continuous, consistent with a founder effect. Island effects and repeated periods of low population density may also have contributed to the low variation. Within continuous habitat, gene flow appears extensive up to 50 km (high assignment rates to source or close neighbour populations; nonsignificant values of pairwise FST), in agreement with movement data. At larger scales (150-250 km), gene flow is reduced (significant pairwise FST) but there is no evidence for isolation by distance. The most substantial genetic structuring was observed for comparisons spanning unsuitable habitat, implying limited dispersal of devils between the well-connected, eastern populations and a smaller northwestern population. The genetic distinctiveness of the northwestern population was reflected in all analyses: unique alleles; multivariate analyses of gene frequency (multidimensional scaling, minimum spanning tree, nearest neighbour); high self-assignment (95%); two distinct populations for Tasmania were detected in isolation by distance and in Bayesian model-based clustering analyses. Marsupial carnivores appear to have stronger population subdivisions than their placental counterparts.     

Genomic signatures of population bottleneck and recovery in Northwest Atlantic pinnipeds

Population increases over the past several decades provide natural settings in which to study the evolutionary processes that occur during bottleneck, growth, and spatial expansion. We used parallel natural experiments of historical decline and subsequent recovery in two sympatric pinniped species in the Northwest Atlantic, the gray seal (Halichoerus grypus atlantica) and harbor seal (Phoca vitulina vitulina), to study the impact of recent demographic change in genomic diversity. Using restriction site‐associated DNA sequencing, we assessed genomic diversity at over 8,700 polymorphic gray seal loci and 3,700 polymorphic harbor seal loci in samples from multiple cohorts collected throughout recovery over the past half‐century. Despite significant differences in the degree of genetic diversity assessed in the two species, we found signatures of historical bottlenecks in the contemporary genomes of both gray and harbor seals. We evaluated temporal trends in diversity across cohorts, as well as compared samples from sites at both the center and edge of a recent gray seal range expansion, but found no significant change in genomewide diversity following recovery. We did, however, find that the variance and degree of allele frequency change measured over the past several decades were significantly different from neutral expectations of drift under population growth. These two cases of well‐described demographic history provide opportunities for critical evaluation of current approaches to simulating and understanding the genetic effects of historical demographic change in natural populations.     

Genetic variation and population structure in Scandinavian wolverine (Gulo gulo) populations

Wolverine (Gulo gulo) numbers in Scandinavia were significantly reduced during the early part of the century as a result of predator removal programmes and hunting. Protective legislation in both Sweden and Norway in the 1960s and 1970s has now resulted in increased wolverine densities in Scandinavia. We report here the development of 15 polymorphic microsatellite markers in wolverine and their use to examine the population sub-structure and genetic variability in free-ranging Scandinavian wolverine populations as well as in a sample of individuals collected before 1970. Significant subdivision between extant populations was discovered, in particular for the small and isolated population of southern Norway, which represents a recent recolonization. Overall genetic variability was found to be lower than previously reported for other mustelids, with only two to five alleles per locus and observed heterozygosities (H(O)) ranging from 0.269 to 0.376 across the examined populations, being lowest in southern Norway. Analysis of the mitochondrial DNA control region revealed no variation throughout the surveyed populations. As the historical sample did not show higher levels of genetic variability, our results are consistent with a reduction in the genetic variation in Scandinavian wolverines that pre-dates the demographic bottleneck observed during the last century. The observed subdivision between populations calls for management caution when issuing harvest quotas, especially for the geographically isolated south Norwegian population.     

Genomic inference accurately predicts the timing and severity of a recent bottleneck in a nonmodel insect population

The analysis of molecular data from natural populations has allowed researchers to answer diverse ecological questions that were previously intractable. In particular, ecologists are often interested in the demographic history of populations, information that is rarely available from historical records. Methods have been developed to infer demographic parameters from genomic data, but it is not well understood how inferred parameters compare to true population history or depend on aspects of experimental design. Here, we present and evaluate a method of SNP discovery using RNA sequencing and demographic inference using the program δaδi, which uses a diffusion approximation to the allele frequency spectrum to fit demographic models. We test these methods in a population of the checkerspot butterfly Euphydryas gillettii. This population was intentionally introduced to Gothic, Colorado in 1977 and has as experienced extreme fluctuations including bottlenecks of fewer than 25 adults, as documented by nearly annual field surveys. Using RNA sequencing of eight individuals from Colorado and eight individuals from a native population in Wyoming, we generate the first genomic resources for this system. While demographic inference is commonly used to examine ancient demography, our study demonstrates that our inexpensive, all‐in‐one approach to marker discovery and genotyping provides sufficient data to accurately infer the timing of a recent bottleneck. This demographic scenario is relevant for many species of conservation concern, few of which have sequenced genomes. Our results are remarkably insensitive to sample size or number of genomic markers, which has important implications for applying this method to other nonmodel systems.     

Genetic diversity within a declining natural population of Ferocactus histrix (DC) Lindsay

Ferocactus histrix is a barrel cactus that is widespread in Mexico. A population located in Llanos de Ojuelos, a semiarid zone representative of many disturbed regions in north‐central Mexico, was studied. Over a period of 10 years (1997 to 2007), the average number of individuals decreased from 21.95 to 3.53 plants per 300 m². A change in population size structure was also registered over this period of time. In 2008, a plot selected on the basis of plant abundance was established within the population and a genetic analysis was conducted with ISTR and ISSR markers. This analysis revealed low levels of genetic diversity (expected heterozygosity (H_E) = 0.073, Shannon index (I) = 0.113 and H_E = 0.178, I = 0.271, respectively) compared with those of most studied cacti species. The genetic diversity between the different life stages was also evaluated, and a gradual decrease in levels of genetic variation was observed from adults to juveniles and seedlings (H_E = 0.130, I = 0.192 to H_E = 0.103, I = 0.157). These differences, however, were not significant. Loci fixation and a decrease in the frequency of rare alleles were observed in seedling and juvenile classes. The decline in genetic variation may be associated with recent bottlenecks experienced by the population of F. histrix. If the sizes of local populations of F. histrix continue to decrease, genetic variation will be gradually lost, and the risk of extinction will increase.     

Impact of a population bottleneck on symmetry and genetic diversity in the northern elephant seal

Abstract The northern elephant seal (NES) suffered a severe population bottleneck towards the end of the nineteenth century. Theoretical expectations for the impact of population bottlenecks include the loss of genetic diversity and a loss of fitness (e.g. through a disruption of developmental stability); however, there are few direct demonstrations in natural populations. Here, we report on the comparison of archive samples collected prior to and following the NES population bottleneck. Measures of genetic diversity show a loss of variation consistent with expectations and suggest a strong disruption in the pattern of allele frequencies following the bottleneck. Measures of bilateral characters show an increase in fluctuating asymmetry.     

Microsatellite analysis of a population crash and bottleneck in the Mauna Kea silversword, Argyroxiphium sandwicense ssp. sandwicense (Asteraceae), and its implications for reintroduction

The Mauna Kea silversword, Argyroxiphium sandwicense ssp. sandwicense, has experienced both a severe population crash associated with an increase in alien ungulate populations on Mauna Kea, and a population bottleneck associated with reintroduction. In this paper, we address the genetic consequences of both demographic events using eight microsatellite loci. The population crash was not accompanied by a significant reduction in number of alleles or heterozygosity. However, the population bottleneck was accompanied by significant reductions in observed number of alleles, effective number of alleles, and expected heterozygosity, though not in observed heterozygosity. The effective size of the population bottleneck was calculated using both observed heterozygosities and allele frequency variances. Both methods corroborated the historical census size of the population bottleneck of at most three individuals. The results suggest that: (i) small populations, even those that result from severe reductions in historical population size and extent, are not necessarily genetically depauperate; and (ii) species reintroduction plans need to be conceived and implemented carefully, with due consideration to the genetic impact of sampling for reintroduction.