Detecting population expansion and decline using microsatellites

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.

The impact of human conflict on the genetics of Mastomys natalensis and Lassa virus in West Africa

Environmental changes have been shown to play an important role in the emergence of new human diseases of zoonotic origin. The contribution of social factors to their spread, especially conflicts followed by mass movement of populations, has not been extensively investigated. Here we reveal the effects of civil war on the phylogeography of a zoonotic emerging infectious disease by concomitantly studying the population structure, evolution and demography of Lassa virus and its natural reservoir, the rodent Mastomys natalensis, in Guinea, West Africa. Analysis of nucleoprotein gene sequences enabled us to reconstruct the evolutionary history of Lassa virus, which appeared 750 to 900 years ago in Nigeria and only recently spread across western Africa (170 years ago). Bayesian demographic inferences revealed that both the host and the virus populations have gone recently through severe genetic bottlenecks. The timing of these events matches civil war-related mass movements of refugees and accompanying environmental degradation. Forest and habitat destruction and human predation of the natural reservoir are likely explanations for the sharp decline observed in the rodent populations, the consequent virus population decline, and the coincident increased incidence of Lassa fever in these regions. Interestingly, we were also able to detect a similar pattern in Nigeria coinciding with the Biafra war. Our findings show that anthropogenic factors may profoundly impact the population genetics of a virus and its reservoir within the context of an emerging infectious disease.

Current and historical drivers of landscape genetic structure differ in core and peripheral salamander populations

With predicted decreases in genetic diversity and greater genetic differentiation at range peripheries relative to their cores, it can be difficult to distinguish between the roles of current disturbance versus historic processes in shaping contemporary genetic patterns. To address this problem, we test for differences in historic demography and landscape genetic structure of coastal giant salamanders (Dicamptodon tenebrosus) in two core regions (Washington State, United States) versus the species' northern peripheral region (British Columbia, Canada) where the species is listed as threatened. Coalescent-based demographic simulations were consistent with a pattern of post-glacial range expansion, with both ancestral and current estimates of effective population size being much larger within the core region relative to the periphery. However, contrary to predictions of recent human-induced population decline in the less genetically diverse peripheral region, there was no genetic signature of population size change. Effects of current demographic processes on genetic structure were evident using a resistance-based landscape genetics approach. Among core populations, genetic structure was best explained by length of the growing season and isolation by resistance (i.e. a ‘flat’ landscape), but at the periphery, topography (slope and elevation) had the greatest influence on genetic structure. Although reduced genetic variation at the range periphery of D. tenebrosus appears to be largely the result of biogeographical history rather than recent impacts, our analyses suggest that inherent landscape features act to alter dispersal pathways uniquely in different parts of the species' geographic range, with implications for habitat management.

Spatial and temporal patterns of neutral and adaptive genetic variation in the endangered African wild dog (Lycaon pictus)

Deciphering patterns of genetic variation within a species is essential for understanding population structure, local adaptation and differences in diversity between populations. Whilst neutrally evolving genetic markers can be used to elucidate demographic processes and genetic structure, they are not subject to selection and therefore are not informative about patterns of adaptive variation. As such, assessments of pertinent adaptive loci, such as the immunity genes of the major histocompatibility complex (MHC), are increasingly being incorporated into genetic studies. In this study, we combined neutral (microsatellite, mtDNA) and adaptive (MHC class II DLA-DRB1 locus) markers to elucidate the factors influencing patterns of genetic variation in the African wild dog (Lycaon pictus); an endangered canid that has suffered extensive declines in distribution and abundance. Our genetic analyses found all extant wild dog populations to be relatively small (N(e)  < 30). Furthermore, through coalescent modelling, we detected a genetic signature of a recent and substantial demographic decline, which correlates with human expansion, but contrasts with findings in some other African mammals. We found strong structuring of wild dog populations, indicating the negative influence of extensive habitat fragmentation and loss of gene flow between habitat patches. Across populations, we found that the spatial and temporal structure of microsatellite diversity and MHC diversity were correlated and strongly influenced by demographic stability and population size, indicating the effects of genetic drift in these small populations. Despite this correlation, we detected signatures of selection at the MHC, implying that selection has not been completely overwhelmed by genetic drift.

Population genetics of the Halys pit viper (Gloydius halys) at the northern distribution limit in Siberia

We studied the genetic structure and diversity of a vulnerable population of the Halys pit viper (Gloydius halys) living at the northern distribution limit of the species (Novosibirsk Region, West Siberia, Russia) and fairly isolated from the nearest known locations (180-200 km). 157 snakes from three locations (separated by 5-11 km) were genotyped using 8 microsatellite loci. In addition, nucleotide sequences of a mitochondrial marker gene (ND4) were obtained from 30 individuals. Two mitochondrial haplotypes differing by one synonymous substitution could be detected. Haplotype H1 had a prevalence of 96.7%. High levels of genetic diversity (mean ) were detected and no evidence of a recent bottleneck in any sampling site could be discovered. We found that G. halys exhibits a low, although significant, level of genetic differentiation on a fine geographic scale (overall FST = 0.013). The degree of differentiation is comparably low with regard to values reported for other vipers and more similar to those in Colubridae. In addition, using local based FST estimates, we detected a significant difference in the extent of genetic drift between the centrally located and marginal local populations as predicted by the stepping-stone model of migration. Thus, from a short-term perspective, primary attention should be given to the peripheral local populations, due to limited immigration and stronger influence of genetic drift. Since all habitat patches are situated along the river and have a linear order it is important to preserve each habitat site to maintain gene flow across the metapopulation.

Impact of human management on the genetic variation of wild pepper, Capsicum annuum var. glabriusculum

Management of wild peppers in Mexico has occurred for a long time without clear phenotypic signs of domestication. However, pre-domestication management could have implications for the population's genetic richness. To test this hypothesis we analysed 27 wild (W), let standing (LS) and cultivated (C) populations, plus 7 samples from local markets (LM), with nine polymorphic microsatellite markers. Two hundred and fifty two alleles were identified, averaging 28 per locus. Allele number was higher in W, and 15 and 40% less in LS and C populations, respectively. Genetic variation had a significant population structure. In W populations, structure was associated with ecological and geographic areas according to isolation by distance. When LM and C populations where included in the analysis, differentiation was no longer apparent. Most LM were related to distant populations from Sierra Madre Oriental, which represents their probable origin. Historical demography shows a recent decline in all W populations. Thus, pre-domestication human management is associated with a significant reduction of genetic diversity and with a loss of differentiation suggesting movement among regions by man. Measures to conserve wild and managed populations should be implemented to maintain the source and the architecture of genetic variation in this important crop relative.

Reconstruction of a beech population bottleneck using archival demographic information and Bayesian analysis of genetic data

Range expansion and contraction has occurred in the history of most species and can seriously impact patterns of genetic diversity. Historical data about range change are rare and generally appropriate for studies at large scales, whereas the individual pollen and seed dispersal events that form the basis of geneflow and colonization generally occur at a local scale. In this study, we investigated range change in Fagus sylvatica on Mont Ventoux, France, using historical data from 1838 to the present and approximate Bayesian computation (ABC) analyses of genetic data. From the historical data, we identified a population minimum in 1845 and located remnant populations at least 200 years old. The ABC analysis selected a demographic scenario with three populations, corresponding to two remnant populations and one area of recent expansion. It also identified expansion from a smaller ancestral population but did not find that this expansion followed a population bottleneck, as suggested by the historical data. Despite a strong support to the selected scenario for our data set, the ABC approach showed a low power to discriminate among scenarios on average and a low ability to accurately estimate effective population sizes and divergence dates, probably due to the temporal scale of the study. This study provides an unusual opportunity to test ABC analysis in a system with a well-documented demographic history and identify discrepancies between the results of historical, classical population genetic and ABC analyses. The results also provide valuable insights into genetic processes at work at a fine spatial and temporal scale in range change and colonization.

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.

The effect of recurrent floods on genetic composition of marble trout populations

A changing global climate can threaten the diversity of species and ecosystems. We explore the consequences of catastrophic disturbances in determining the evolutionary and demographic histories of secluded marble trout populations in Slovenian streams subjected to weather extremes, in particular recurrent flash floods and debris flows causing massive mortalities. Using microsatellite data, a pattern of extreme genetic differentiation was found among populations (global F(ST) of 0.716), which exceeds the highest values reported in freshwater fish. All locations showed low levels of genetic diversity as evidenced by low heterozygosities and a mean of only 2 alleles per locus, with few or no rare alleles. Many loci showed a discontinuous allele distribution, with missing alleles across the allele size range, suggestive of a population contraction. Accordingly, bottleneck episodes were inferred for all samples with a reduction in population size of 3-4 orders of magnitude. The reduced level of genetic diversity observed in all populations implies a strong impact of genetic drift, and suggests that along with limited gene flow, genetic differentiation might have been exacerbated by recurrent mortalities likely caused by flash flood and debris flows. Due to its low evolutionary potential the species might fail to cope with an intensification and altered frequency of flash flood events predicted to occur with climate change.

Demographic processes underlying subtle patterns of population structure in the scalloped hammerhead shark, Sphyrna lewini

Genetic diversity (θ), effective population size (N(e)), and contemporary levels of gene flow are important parameters to estimate for species of conservation concern, such as the globally endangered scalloped hammerhead shark, Sphyrna lewini. Therefore, we have reconstructed the demographic history of S. lewini across its Eastern Pacific (EP) range by applying classical and coalescent population genetic methods to a combination of 15 microsatellite loci and mtDNA control region sequences. In addition to significant population genetic structure and isolation-by-distance among seven coastal sites between central Mexico and Ecuador, the analyses revealed that all populations have experienced a bottleneck and that all current values of θ are at least an order of magnitude smaller than ancestral θ, indicating large decreases in N(e) (θ = 4N(e)μ), where μ is the mutation rate. Application of the isolation-with-migration (IM) model showed modest but significant genetic connectivity between most sampled sites (point estimates of Nm = 0.1-16.7), with divergence times (t) among all populations significantly greater than zero. Using a conservative (i.e., slow) fossil-based taxon-specific phylogenetic calibration for mtDNA mutation rates, posterior probability distributions (PPDs) for the onset of the decline in N(e) predate modern fishing in this region. The cause of decline over the last several thousand years is unknown but is highly atypical as a post-glacial demographic history. Regardless of the cause, our data and analyses suggest that S. lewini was far more abundant throughout the EP in the past than at present.

A century of genetic change and metapopulation dynamics in the Galápagos warbler finches (Certhidea)

Populations that are connected by immigrants play an important role in evolutionary and conservation biology, yet we have little direct evidence of how such metapopulations change genetically over evolutionary time. We compared historic (1894-1906) to modern (1988-2006) genetic variation in 11 populations of warbler finches at 14 microsatellite loci. Although several lines of evidence suggest that Darwin's finches may be in decline, we found that the genetic diversity of warbler finches has not generally declined, and broad-scale patterns of variation remained similar over time. Contrary to expectations, inferred population sizes have generally increased over time (6-8%) as have immigration rates (8-16%), which may reflect a recent increase in the frequency and intensity of El Niño events. Individual island populations showed significant declines (18-19%) and also substantial gains (18-20%) in allelic richness over time. Changes in genetic diversity were correlated with changes in immigration rates, but did not correspond to population size or human disturbance. These results reflect the expected stabilizing properties of whole metapopulations over time. However, the dramatic and unpredictable changes observed in individual populations during this short time interval suggests that care should be taken when monitoring individual population fragments with snapshots of genetic variation.

Guns, germs and dogs: On the origin of Leishmania chagasi

The evolutionary history of Leishmania chagasi, the aetiological agent of visceral leishmaniasis in South America has been widely debated. This study addresses the problem of the origin of L. chagasi, its timing and demography with fast evolving genetic markers, a suite of Bayesian clustering algorithms and coalescent modelling. Here, using 14 microsatellite markers, 450 strains from the Leishmania donovani complex, we show that the vast majority of the Central and South American L. chagasi were nested within the Portuguese Leishmania infantum clade. Moreover, L. chagasi allelic richness was half that of their Old World counterparts. The bottleneck signature was estimated to be about 500 years old and the settlement of L. chagasi in the New World, probably via infected dogs, was accompanied by a thousand-fold population decrease. Visceral leishmaniasis, lethal if untreated, is therefore one more disease that the Conquistadores brought to the New World.

Inferring population decline and expansion from microsatellite data: a simulation-based evaluation of the Msvar method

Reconstructing the demographic history of populations is a central issue in evolutionary biology. Using likelihood-based methods coupled with Monte Carlo simulations, it is now possible to reconstruct past changes in population size from genetic data. Using simulated data sets under various demographic scenarios, we evaluate the statistical performance of Msvar, a full-likelihood Bayesian method that infers past demographic change from microsatellite data. Our simulation tests show that Msvar is very efficient at detecting population declines and expansions, provided the event is neither too weak nor too recent. We further show that Msvar outperforms two moment-based methods (the M-ratio test and Bottleneck) for detecting population size changes, whatever the time and the severity of the event. The same trend emerges from a compilation of empirical studies. The latest version of Msvar provides estimates of the current and the ancestral population size and the time since the population started changing in size. We show that, in the absence of prior knowledge, Msvar provides little information on the mutation rate, which results in biased estimates and/or wide credibility intervals for each of the demographic parameters. However, scaling the population size parameters with the mutation rate and scaling the time with current population size, as coalescent theory requires, significantly improves the quality of the estimates for contraction but not for expansion scenarios. Finally, our results suggest that Msvar is robust to moderate departures from a strict stepwise mutation model.

Joint inference of microsatellite mutation models, population history and genealogies using transdimensional Markov Chain Monte Carlo

We provide a framework for Bayesian coalescent inference from microsatellite data that enables inference of population history parameters averaged over microsatellite mutation models. To achieve this we first implemented a rich family of microsatellite mutation models and related components in the software package BEAST. BEAST is a powerful tool that performs Bayesian MCMC analysis on molecular data to make coalescent and evolutionary inferences. Our implementation permits the application of existing nonparametric methods to microsatellite data. The implemented microsatellite models are based on the replication slippage mechanism and focus on three properties of microsatellite mutation: length dependency of mutation rate, mutational bias toward expansion or contraction, and number of repeat units changed in a single mutation event. We develop a new model that facilitates microsatellite model averaging and Bayesian model selection by transdimensional MCMC. With Bayesian model averaging, the posterior distributions of population history parameters are integrated across a set of microsatellite models and thus account for model uncertainty. Simulated data are used to evaluate our method in terms of accuracy and precision of θ estimation and also identification of the true mutation model. Finally we apply our method to a red colobus monkey data set as an example.

Genetic differentiation of European grayling (Thymallus thymallus) populations in Serbia, based on mitochondrial and nuclear DNA analyses

Background

The structure and diversity of grayling (Thymallus thymallus) populations have been well studied in most of its native habitat; however the southernmost populations of the Balkan Peninsula remain largely unexplored. The purpose of this study was to assess the genetic diversity of Serbian grayling populations, detect the impact of stocking and provide guidelines for conservation and management.

Methods

Eighty grayling individuals were collected from four rivers (Ibar, Lim, Drina and Rzav). The mitochondrial DNA control region (CR; 595 bp of the 3'end and 74 bp of flanking tRNA) and the ATP6 gene (630 bp fragment) were sequenced for 20 individuals (five from each locality). In addition, all individuals were genotyped with 12 microsatellite loci. The diversity and structure of the populations as well as the recent and ancient population declines were studied using specialized software.

Results

We detected three new haplotypes in the mtDNA CR and four haplotypes in the ATP6 gene of which three had not been described before. Previously, one CR haplotype and two ATP6 gene haplotypes had been identified as allochthonous, originating from Slovenia. Reconstruction of phylogenetic relations placed the remaining two CR haplotypes from the River Danube drainage of Serbia into a new clade, which is related to the previously described sister Slovenian clade. These two clades form a new Balkan clade. Microsatellite marker analysis showed that all four populations are genetically distinct from each other without any sign of intra-population structure, although stocking of the most diverse population (Drina River) was confirmed by mtDNA analysis. Recent and historical population declines of Serbian grayling do not differ from those of other European populations.

Conclusions

Our study shows that (1) the Ibar, Lim and Drina Rivers grayling populations are genetically distinct from populations outside of Serbia and thus should be managed as native populations in spite of some introgression in the Drina River population and (2) the Rzav River population is not appropriate for further stocking activities since it originates from stocked Slovenian grayling. However, the Rzav River population does not represent an immediate danger to other populations because it is physically isolated from these.

Evolution and population genetics of exotic and re-emerging pathogens: novel tools and approaches

Given human population growth and accelerated global trade, the rate of emergence of exotic plant pathogens is bound to increase. Understanding the processes that lead to the emergence of new pathogens can help manage emerging epidemics. Novel tools for analyzing population genetic variation can be used to infer the evolutionary history of populations or species, allowing for the unprecedented reconstruction of the demographic history of pathogens. Specifically, recent advances in the application of coalescent, maximum likelihood (ML), and Bayesian methods to population genetic data combined with increasing availability of affordable sequencing and parallel computing have created the opportunity to apply these methods to a broad range of questions regarding the evolution of emerging pathogens. These approaches are particularly powerful when used to test multiple competing hypotheses. We provide several examples illustrating how coalescent analysis provides critical insights into understanding migration pathways as well as processes of divergence, speciation, and recombination.

The role of humans in the diversification of a threatened island raptor

BACKGROUND

Anthropogenic habitat modifications have led to the extinction of many species and have favoured the expansion of others. Nonetheless, the possible role of humans as a diversifying force in vertebrate evolution has rarely been considered, especially for species with long generation times. We examine the influence that humans have had on the colonization and phenotypic and genetic differentiation of an insular population of a long-lived raptor species, the Egyptian vulture (Neophron percnopterus).

RESULTS

The morphological comparison between the Canarian Egyptian vultures and the main and closest population in Western Europe (Iberia) indicated that insular vultures are significantly heavier (16%) and larger (about 3%) than those from Iberia. Bayesian and standard genetic analyses also showed differentiation (FST = 0.11, p < 0.01). The inference of changes in the effective size of the Canarian deme, using two likelihood-based Bayesian approaches, suggested that the establishment of this insular population took place some 2500 years ago, matching the date of human colonization. This is consistent with the lack of earlier fossils.

CONCLUSIONS

Archaeological remains show that first colonizers were Berber people from northern Africa who imported goats. This new and abundant food source could have allowed vultures to colonize, expand and adapt to the island environment. Our results suggest that anthropogenic environmental change can induce diversification and that this process may take place on an ecological time scale (less than 200 generations), even in the case of a long-lived species.

Genetic evidence of recent population contraction in the southernmost population of giant pandas

Anthropogenic habitat loss and fragmentation have been implicated in the endangerment and extinction of many species. Here we assess genetic variation and demographic history in the southernmost population of giant pandas (Ailuropoda melanoleuca) that continues to be threatened by habitat degradation and fragmentation, using noninvasive genetic sampling, mitochondrial control region sequence and 12 microsatellite loci. Compared to other giant panda populations, this population has medium-level genetic diversity based on the measure of both mitochondrial and nuclear markers. Mitochondrial DNA-based demographic analyses revealed that no historical population expansion or contraction has occurred, indicating a relatively stable population size. However, a Bayesian-coalescent method based on the observed allele distribution and allele frequencies of microsatellite clearly did detect, quantify and date a recent decrease in population size. Overall, the results indicate that a population contraction in the order of 95-96% has taken place over the last 910-999 years and is most likely due to anthropogenic habitat loss. These findings highlight the need for a greater focus on habitat protection and restoration for the long-term survival of this giant panda population.

The number of alleles at a microsatellite defines the allele frequency spectrum and facilitates fast accurate estimation of theta

Theoretical work focused on microsatellite variation has produced a number of important results, including the expected distribution of repeat sizes and the expected squared difference in repeat size between two randomly selected samples. However, closed-form expressions for the sampling distribution and frequency spectrum of microsatellite variation have not been identified. Here, we use coalescent simulations of the stepwise mutation model to develop gamma and exponential approximations of the microsatellite allele frequency spectrum, a distribution central to the description of microsatellite variation across the genome. For both approximations, the parameter of biological relevance is the number of alleles at a locus, which we express as a function of θ, the population-scaled mutation rate, based on simulated data. Discovered relationships between θ, the number of alleles, and the frequency spectrum support the development of three new estimators of microsatellite θ. The three estimators exhibit roughly similar mean squared errors (MSEs) and all are biased. However, across a broad range of sample sizes and θ values, the MSEs of these estimators are frequently lower than all other estimators tested. The new estimators are also reasonably robust to mutation that includes step sizes greater than one. Finally, our approximation to the microsatellite allele frequency spectrum provides a null distribution of microsatellite variation. In this context, a preliminary analysis of the effects of demographic change on the frequency spectrum is performed. We suggest that simulations of the microsatellite frequency spectrum under evolutionary scenarios of interest may guide investigators to the use of relevant and sometimes novel summary statistics.

Conservation implications of drastic reductions in the smallest and most isolated populations of giant pandas

In conservation biology, understanding the causes of endangerment is a key step to devising effective conservation strategies. We used molecular evidence (coalescent simulations of population changes from microsatellite data) and historical information (habitat and human population changes) to investigate how the most-isolated populations of giant pandas (Ailuropoda melanoleuca) in the Xiaoxiangling Mountains became highly endangered. These populations experienced a strong, recent demographic reduction (60-fold), starting approximately 250 years BP. Explosion of the human population and use of non-native crop species at the peak of the Qing Empire resulted in land-use changes, deforestation, and habitat fragmentation, which are likely to have led to the drastic reduction of the most-isolated populations of giant pandas. We predict that demographic, genetic, and environmental factors will lead to extinction of giant pandas in the Xiaoxiangling Mountains in the future if the population remains isolated. Therefore, a targeted conservation action--translocation--has been proposed and is being implemented by the Chinese government.

Genetic structure of Anopheles (Nyssorhynchus) marajoara (Diptera: Culicidae) in Colombia

Five Anopheles marajoara Galvão and Damasceno populations, representing diverse ecological conditions, were sampled throughout Colombia and analyzed using nine hypervariable DNA microsatellite loci. The overall genetic diversity (H = 0.58) was lower than that determined for some Brazilian populations using the same markers. The Caquetá population (Colombia) had the lowest gene diversity (H = 0.48), and it was the only population at Hardy-Weinberg equilibrium. Hardy-Weinberg disequilibrium in the remaining four populations was probably caused by the Wahlund effect. The assignment analyses showed two incompletely isolated gene pools separated by the Eastern Andean cordillera. However, other possible geographical barriers (rivers and other mountains) did not play any role in the moderate genetic heterogeneity found among these populations (F(ST) = 0.069). These results are noteworthy, because this species is a putative malaria vector in Colombia.

Genetic signals of ancient decline in Aleppo pine populations at the species' southwestern margins in the Mediterranean Basin

Microsatellite markers were used to characterize the structure of genetic diversity in natural Moroccan Aleppo pine (Pinus halepensis Mill.) populations, the most southwesterly marginal populations of the species in the Mediterranean Basin. Twenty-two natural populations and one artificial population, located in four regions covering most of the natural range of P. halepensis in the country, were sampled. Across this range, towards the south and west (and towards high altitudes) the populations become increasingly discrete and discontinuous. The nuclear microsatellite marker analysis suggests that a large proportion of the Aleppo pines in Morocco have derived from a single genetic lineage, represented by a central group of 11 of the examined populations located in the High and Middle Atlas Mountains. In addition, two smaller groups, represented by the marginal southwestern High Atlas populations, and three still smaller north / northeastern groups of populations located in the Rif and northeast Middle Atlas Mountains, could be genetically distinguished. Further, coalescence analysis of historical demographic population patterns suggests that ancient bottlenecks occurred in all of the natural populations. However, the population differentiation and genetic diversity levels we found were good (F(st) =15.47), presumably because of the species' good potential for long-distance dispersal of seeds and high invasive capacity, which appear to have maintained a state of stable near-equilibrium, meta-population dynamics since ancient times.

The confounding effects of population structure, genetic diversity and the sampling scheme on the detection and quantification of population size changes

The idea that molecular data should contain information on the recent evolutionary history of populations is rather old. However, much of the work carried out today owes to the work of the statisticians and theoreticians who demonstrated that it was possible to detect departures from equilibrium conditions (e.g., panmictic population/mutation-drift equilibrium) and interpret them in terms of deviations from neutrality or stationarity. During the last 20 years the detection of population size changes has usually been carried out under the assumption that samples were obtained from populations that can be approximated by a Wright-Fisher model (i.e., assuming panmixia, demographic stationarity, etc.). However, natural populations are usually part of spatial networks and are interconnected through gene flow. Here we simulated genetic data at mutation and migration-drift equilibrium under an n-island and a stepping-stone model. The simulated populations were thus stationary and not subject to any population size change. We varied the level of gene flow between populations and the scaled mutation rate. We also used several sampling schemes. We then analyzed the simulated samples using the Bayesian method implemented in MSVAR, the Markov Chain Monte Carlo simulation program, to detect and quantify putative population size changes using microsatellite data. Our results show that all three factors (genetic differentiation/gene flow, genetic diversity, and the sampling scheme) play a role in generating false bottleneck signals. We also suggest an ad hoc method to counter this effect. The confounding effect of population structure and of the sampling scheme has practical implications for many conservation studies. Indeed, if population structure is creating "spurious" bottleneck signals, the interpretation of bottleneck signals from genetic data might be less straightforward than it would seem, and several studies may have overestimated or incorrectly detected bottlenecks in endangered species.

A microsatellite-based analysis for the detection of selection on BTA1 and BTA20 in northern Eurasian cattle (Bos taurus) populations

BACKGROUND

Microsatellites surrounding functionally important candidate genes or quantitative trait loci have received attention as proxy measures of polymorphism level at the candidate loci themselves. In cattle, selection for economically important traits is a long-term strategy and it has been reported that microsatellites are linked to these important loci.

METHODS

We have investigated the variation of seven microsatellites on BTA1 (Bos taurus autosome 1) and 16 on BTA20, using bovine populations of typical production types and horn status in northern Eurasia. Genetic variability of these loci and linkage disequilibrium among these loci were compared with those of 28 microsatellites on other bovine chromosomes. Four different tests were applied to detect molecular signatures of selection.

RESULTS

No marked difference in locus variability was found between microsatellites on BTA1, BTA20 and the other chromosomes in terms of different diversity indices. Average D' values of pairwise syntenic markers (0.32 and 0.28 across BTA 1 and BTA20 respectively) were significantly (P < 0.05) higher than for non-syntenic markers (0.15). The Ewens-Watterson test, the Beaumont and Nichol's modified frequentist test and the Bayesian FST-test indicated elevated or decreased genetic differentiation, at SOD1 and AGLA17 markers respectively, deviating significantly (P < 0.05) from neutral expectations. Furthermore, lnRV, lnRH and lnRtheta' statistics were used for the pairwise population comparison tests and were significantly less variable in one population relative to the other, providing additional evidence of selection signatures for two of the 51 loci. Moreover, the three Finnish native populations showed evidence of subpopulation divergence at SOD1 and AGLA17. Our data also indicate significant intergenic linkage disequilibrium around the candidate loci and suggest that hitchhiking selection has played a role in shaping the pattern of observed linkage disequilibrium.

CONCLUSION

Hitchhiking due to tight linkage with alleles at candidate genes, e.g. the POLL gene, is a possible explanation for this pattern. The potential impact of selective breeding by man on cattle populations is discussed in the context of selection effects. Our results also suggest that a practical approach to detect loci under selection is to simultaneously apply multiple neutrality tests based on different assumptions and estimations.

An examination of genetic diversity and effective population size in Atlantic salmon populations

Effective population size (Ne) is an important parameter in the conservation of genetic diversity. Comparative studies of empirical data that gauge the relative accuracy of Ne methods are limited, and a better understanding of the limitations and potential of Ne estimators is needed. This paper investigates genetic diversity and Ne in four populations of wild anadromous Atlantic salmon (Salmo salar L.) in Europe, from the Rivers Oir and Scorff (France) and Spey and Shin (Scotland). We aimed to understand present diversity and historical processes influencing current population structure. Our results showed high genetic diversity for all populations studied, despite their wide range of current effective sizes. To improve understanding of high genetic diversity observed in the populations with low effective size, we developed a model predicting present diversity as a function of past demographic history. This suggested that high genetic diversity could be explained by a bottleneck occurring within recent centuries rather than by gene flow. Previous studies have demonstrated the efficiency of coalescence models to estimate Ne. Using nine subsets from 37 microsatellite DNA markers from the four salmon populations, we compared three coalescence estimators based on single and dual samples. Comparing Ne estimates confirmed the efficiency of increasing the number and variability of microsatellite markers. This efficiency was more accentuated for the smaller populations. Analysis with low numbers of neutral markers revealed uneven distributions of allelic frequencies and overestimated short-term Ne. In addition, we found evidence of artificial stock enhancement using native and non-native origin. We propose estimates of Ne for the four populations, and their applications for salmon conservation and management are discussed.

Inference of population history by coupling exploratory and model-driven phylogeographic analyses

Understanding the nature, timing and geographic context of historical events and population processes that shaped the spatial distribution of genetic diversity is critical for addressing questions relating to speciation, selection, and applied conservation management. Cladistic analysis of gene trees has been central to phylogeography, but when coupled with approaches that make use of different components of the information carried by DNA sequences and their frequencies, the strength and resolution of these inferences can be improved. However, assessing concordance of inferences drawn using different analytical methods or genetic datasets, and integrating their outcomes, can be challenging. Here we overview the strengths and limitations of different types of genetic data, analysis methods, and approaches to historical inference. We then turn our attention to the potentially synergistic interactions among widely-used and emerging phylogeographic analyses, and discuss some of the ways that spatial and temporal concordance among inferences can be assessed. We close this review with a brief summary and outlook on future research directions.

Genetic diversity, structure, and demographic change in tanoak, Lithocarpus densiflorus (Fagaceae), the most susceptible species to sudden oak death in California

Knowledge of population genetic structure of tanoak (Lithocarpus densiflorus) is of interest to pathologists seeking natural variation in resistance to sudden oak death disease, to resource managers who need indications of conservation priorities in this species now threatened by the introduced pathogen (Phytophthora ramorum), and to biologists with interests in demographic processes that have shaped plant populations. We investigated population genetic structure using nuclear and chloroplast DNA (cpDNA) and inferred the effects of past population demographic processes and contemporary gene flow. Our cpDNA results revealed a strong pattern of differentiation of four regional groups (coastal California, southern Oregon, Klamath mountains, and Sierra Nevada). The chloroplast haplotype phylogeny suggests relatively deep divergence of Sierra Nevada and Klamath populations from those of coastal California and southern Oregon. A widespread coastal California haplotype may have resulted from multiple refugial sites during the Last Glacial Maximum or from rapid recolonization from few refugia. Analysis of nuclear microsatellites suggests two major groups: (1) central coastal California and (2) Sierra Nevada/Klamath/southern Oregon and an area of admixture in north coastal California. The low level of nuclear differentiation is likely to be due to pollen gene flow among populations during postglacial range expansion.