Long-term survival despite low genetic diversity in the critically endangered Madagascar fish-eagle

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.

Evaluation of the Impact of Population Management on the Genetic Parameters of Selected Spiral-Horned Antelopes

The rapid loss of biodiversity and the associated reduction and fragmentation of habitats means that ex situ populations have become an important part of species conservation. These populations, which are often established from a small number of founders, require careful management to avoid the negative effects of genetic drift and inbreeding. Although the inclusion of molecular data is recommended, their availability for captive breeding management remains limited. The aim of this study was to evaluate the relationship between the levels of genetic diversity in six spiral-horned antelope taxa bred under human care and their respective management strategies, conservation status, demography, and geographic origin, using 10 nuclear DNA microsatellite loci and mitochondrial control region DNA sequences. Our findings include associations between genetic diversity and management intensity but also with the diversity and contribution of wild populations to captive founders, with some populations apparently composed of animals from divergent wild lineages elevating captive genetic diversity. When population sizes are large, the potential advantages of maximizing genetic diversity in widely outcrossed populations may need careful consideration with respect to the potential disruption of adaptive diversity. Genetic data serve as a robust tool for managing captive populations, yet their interpretation necessitates a comprehensive understanding of species biology and history.

Whole-genome survey reveals extensive variation in genetic diversity and inbreeding levels among peregrine falcon subspecies

In efforts to prevent extinction, resource managers are often tasked with increasing genetic diversity in a population of concern to prevent inbreeding depression or improve adaptive potential in a changing environment. The assumption that all small populations require measures to increase their genetic diversity may be unwarranted, and limited resources for conservation may be better utilized elsewhere. We test this assumption in a case study focused on the peregrine falcon (Falco peregrinus), a cosmopolitan circumpolar species with 19 named subspecies. We used whole-genome resequencing to generate over two million single nucleotide polymorphisms (SNPs) from multiple individuals of all peregrine falcon subspecies. Our analyses revealed extensive variation among subspecies, with many island-restricted and nonmigratory populations possessing lower overall genomic diversity, elevated inbreeding coefficients (F _ROH)-among the highest reported, and extensive runs of homozygosity (ROH) compared to mainland and migratory populations. Similarly, the majority of subspecies that are either nonmigratory or restricted to islands show a much longer history of low effective population size (N _e). While mutational load analyses indicated an increased proportion of homozygous-derived deleterious variants (i.e., drift load) among nonmigrant and island populations compared to those that are migrant or reside on the mainland, no significant differences in the proportion of heterozygous deleterious variants (i.e., inbreeding load) was observed. Our results provide evidence that high levels of inbreeding may not be an existential threat for some populations or taxa. Additional factors such as the timing and severity of population declines are important to consider in management decisions about extinction potential.

Impact of model assumptions on demographic inferences: the case study of two sympatric mouse lemurs in northwestern Madagascar

BACKGROUND

Quaternary climate fluctuations have been acknowledged as major drivers of the geographical distribution of the extraordinary biodiversity observed in tropical biomes, including Madagascar. The main existing framework for Pleistocene Malagasy diversification assumes that forest cover was strongly shaped by warmer Interglacials (leading to forest expansion) and by cooler and arid glacials (leading to forest contraction), but predictions derived from this scenario for forest-dwelling animals have rarely been tested with genomic datasets.

RESULTS

We generated genomic data and applied three complementary demographic approaches (Stairway Plot, PSMC and IICR-simulations) to infer population size and connectivity changes for two forest-dependent primate species (Microcebus murinus and M. ravelobensis) in northwestern Madagascar. The analyses suggested major demographic changes in both species that could be interpreted in two ways, depending on underlying model assumptions (i.e., panmixia or population structure). Under panmixia, the two species exhibited larger population sizes across the Last Glacial Maximum (LGM) and towards the African Humid Period (AHP). This peak was followed by a population decline in M. ravelobensis until the present, while M. murinus may have experienced a second population expansion that was followed by a sharp decline starting 3000 years ago. In contrast, simulations under population structure suggested decreasing population connectivity between the Last Interglacial and the LGM for both species, but increased connectivity during the AHP exclusively for M. murinus.

CONCLUSION

Our study shows that closely related species may differ in their responses to climatic events. Assuming that Pleistocene climatic conditions in the lowlands were similar to those in the Malagasy highlands, some demographic dynamics would be better explained by changes in population connectivity than in population size. However, changes in connectivity alone cannot be easily reconciled with a founder effect that was shown for M. murinus during its colonization of the northwestern Madagascar in the late Pleistocene. To decide between the two alternative models, more knowledge about historic forest dynamics in lowland habitats is necessary. Altogether, our study stresses that demographic inferences strongly depend on the underlying model assumptions. Final conclusions should therefore be based on a comparative evaluation of multiple approaches.

Past environmental changes affected lemur population dynamics prior to human impact in Madagascar

Quaternary climatic changes have been invoked as important drivers of species diversification worldwide. However, the impact of such changes on vegetation and animal population dynamics in tropical regions remains debated. To overcome this uncertainty, we integrated high-resolution paleoenvironmental reconstructions from a sedimentary record covering the past 25,000 years with demographic inferences of a forest-dwelling primate species (Microcebus arnholdi), in northern Madagascar. Result comparisons suggest that climate changes through the African Humid Period (15.2 - 5.5 kyr) strongly affected the demographic dynamics of M. arnholdi. We further inferred a population decline in the last millennium which was likely shaped by the combination of climatic and anthropogenic impacts. Our findings demonstrate that population fluctuations in Malagasy wildlife were substantial prior to a significant human impact. This provides a critical knowledge of climatically driven, environmental and ecological changes in the past, which is essential to better understand the dynamics and resilience of current biodiversity.

The tale of springs and streams: how different aquatic ecosystems impacted the mtDNA population structure of two riffle beetles in the Western Carpathians

The Western Carpathians are a particularly interesting part of the Carpathian Arc. According to recent molecular data upon aquatic and terrestrial taxa, this mountain area is an important biodiversity hotspot of Europe. Moreover, the W Carpathians include rich systems of karst springs inhabited by specific fauna, where molecular diversity and phylogeographic patterns are yet to be fully explored. Our study aims to compare population genetic structure and molecular diversity of two related and commonly co-occurring riffle beetles, Elmis aenea (PWJ Müller, 1806) and Limnius perrisi (Dufour, 1843) in the springs and streams of the W Carpathians using the mitochondrial DNA barcoding fragment of the cytochrome c oxidase subunit I gene (COI). The relatively stable thermal and chemical conditions of springs throughout unfavourable climatic settings make these highly specific lotic systems potentially ideal for a long-term survival of some aquatic biota. Populations of both elmid species were relatively homogeneous genetically, with a single dominant haplotype. However, we revealed that E. aenea significantly dominated in the springs, while L. perrisi preferred streams. Relative isolation of the springs and their stable conditions were reflected in significantly higher molecular diversity of the E. aenea population in comparison to L. perrisi. The results of Bayesian Skyline Plot analysis also indicated the exceptional position of springs regarding maintaining the population size of E. aenea. On the other hand, it seems that streams in the W Carpathians provide more effective dispersal channels for L. perrisi, whose population expanded much earlier compared to E. aenea. Present study points out that different demographic histories of these two closely related elmid species are manifested by their different habitat preference and molecular diversity.

High genetic diversity and low differentiation retained in the European fragmented and declining Greater Spotted Eagle (Clanga clanga) population

Characterising genetic diversity and structure of populations is essential for effective conservation of threatened species. The Greater Spotted Eagle (Clanga clanga), a large and globally vulnerable raptor, is extinct or in severe decline in most of its previous range in Europe. We assessed whether the remnants of European population are genetically impoverished, and isolated from each other. We evaluated levels of genetic diversity and population structuring by sequencing mitochondrial pseudo-control region and 10 introns from various nuclear genes, and estimated length diversity in 23 microsatellite markers. The European population has expanded since the late Pleistocene, and does not exhibit signs of a recent population bottleneck. The global genetic diversity in Europe was rather similar to that detected in other similar species. Microsatellites suggested shallow but significant differentiation between the four extant populations in Estonia, Poland, Belarus and Russia (Upper Volga region) populations, but introns and mtDNA showed that only the Estonian population differed from the others. Mitochondrial diversity was highest in the northernmost Estonian population, introns suggested lower diversity in Upper Volga, microsatellites indicated equal diversity among populations. A recent bottleneck was detected in Poland, which is consistent with the observed repopulation of the region. We conclude that significant gene flow and high genetic diversity are retained in the fragmented Greater Spotted Eagle populations; there is currently no need for genetic augmentation in Europe.

Understanding population structure and historical demography of Litsea auriculata (Lauraceae), an endangered species in east China

Detecting how historical and contemporary factors contribute to genetic divergence and genetic structure is a central question in ecology and evolution. We examine this question by intergrating population genetics with ecological niche modelling of Litsea auriculata (Lauraceae), which is endangered and native to east China. Geographical and environmental factors including climatic fluctuations since the last glacial maximum (LGM) have also contribute to population demography and patterns of genetic structure. L. auriculata populations underwent expansion after divergence and dramatically decreased to the current small size with relative population bottlenecks due to climate changes. Populations separated by physical geographical barrier including geographic distance and Yangtze River, as a result contemporary gene flow among L. auriculata populations showed drastic declines in comparison with historical gene flow, resulting in a high level of population divergence. Thus, patterns of genetic structure of L. auriculata can result from both geographic and environmental factors including climate changes. This information is helpful in forming conservation strategies for L. auriculata in China.

Population structure and historical demography of Dipteronia dyeriana (Sapindaceae), an extremely narrow palaeoendemic plant from China: implications for conservation in a biodiversity hot spot

Inferring past demography is a central question in evolutionary and conservation biology. It is, however, sometimes challenging to disentangle their roles of contemporary versus historical processes in shaping the current patterns of genetic variation in endangered species. In this study, we used both chloroplast microsatellite (cpSSR) loci and nuclear microsatellite (nSSR) loci to assess the levels of genetic differentiation, genetic effective population size, contemporary/historical levels of gene flow and demographic history for five populations sampled across the range of Dipteronia dyeriana, an endangered palaeoendemism from Southwestern China. We found that D. dyeriana had a mixed pattern of moderate genetic diversity and high inbreeding. Bayesian clustering divided D. dyeriana populations into two nSSR genetic clusters. Coalescent-based approximate Bayesian computation analyses suggest the western and eastern groups of D. dyeriana likely persisted in a long-term refuge in Southern China since the beginning of the last glacial period, whereas increasingly colder and arid climates at the onset of the last glacial maximum might have fostered the fragmentation of D. dyeriana within refugia. Following their divergence, the western group kept relatively stable effective population size, whereas the eastern group had experienced 500-fold population expansion during the Holocene. Although clear loss of genetic diversity by human activities was not suggested, recent habitat fragmentation has led to a reduction of population connectivity and increased genetic differentiation by ongoing genetic drift in isolated populations, possibly owing to decreased population size in recent dozen years. Finally, we discussed the implications of these results on conservation policies.

Determining causes of genetic isolation in a large carnivore (Ursus americanus) population to direct contemporary conservation measures

The processes leading to genetic isolation influence a population’s local extinction risk, and should thus be identified before conservation actions are implemented. Natural or human-induced circumstances can result in historical or contemporary barriers to gene flow and/or demographic bottlenecks. Distinguishing between these hypotheses can be achieved by comparing genetic diversity and differentiation in isolated vs. continuous neighboring populations. In Ontario, American black bears (Ursus americanus) are continuously distributed, genetically diverse, and exhibit an isolation-by-distance structuring pattern, except on the Bruce Peninsula (BP). To identify the processes that led to the genetic isolation of BP black bears, we modelled various levels of historical and contemporary migration and population size reductions using forward simulations. We compared simulation results with empirical genetic indices from Ontario black bear populations under different levels of geographic isolation, and conducted additional simulations to determine if translocations could help achieve genetic restoration. From a genetic standpoint, conservation concerns for BP black bears are warranted because our results show that: i) a recent demographic bottleneck associated with recently reduced migration best explains the low genetic diversity on the BP; and ii) under sustained isolation, BP black bears could lose between 70% and 80% of their rare alleles within 100 years. Although restoring migration corridors would be the most effective method to enhance long-term genetic diversity and prevent inbreeding, it is unrealistic to expect connectivity to be re-established. Current levels of genetic diversity could be maintained by successfully translocating 10 bears onto the peninsula every 5 years. Such regular translocations may be more practical than landscape restoration, because areas connecting the peninsula to nearby mainland black bear populations have been irreversibly modified by humans, and form strong barriers to movement.

Reduction of Genetic Diversity of the Harpy Eagle in Brazilian Tropical Forests

Habitat loss and fragmentation intensify the effects of genetic drift and endogamy, reducing genetic variability of populations with serious consequences for wildlife conservation. The Harpy Eagle (Harpia harpyja) is a forest dwelling species that is considered near threatened and suffers from habitat loss in the forests of the Neotropical region. In this study, 72 historical and current samples were assessed using eight autosomal microsatellite markers to investigate the distribution of genetic diversity of the Harpy Eagle of the Amazonian and Atlantic forests in Brazil. The results showed that the genetic diversity of Harpy Eagle decreased in the regions where deforestation is intense in the southern Amazon and Atlantic Forest.

Effects of Forest Fragmentation on Genetic Diversity of the Critically Endangered Primate, the Pied Tamarin (Saguinus bicolor): Implications for Conservation

We analyzed DNA at 9 microsatellite loci from hair samples of 73 pied tamarins (Saguinus bicolor) located in 3 urban forest fragments and a biological reserve in the city of Manaus, Amazonas, Brazil. The forest fragments had become isolated from the continuous forest 6-15 years prior to the time of sampling. Tests for reduction in population size showed that all groups from the urban forest fragments had undergone genetic bottlenecks. Pied tamarins in this region historically formed one biological population, and the fragments were connected by high levels of gene flow. These results indicate the need to implement a conservation plan that allows for connectivity between the urban fragments, as well as protection from further constriction. Such connectivity could be achieved via the creation and protection of corridors. In addition to the current population trends explained by anthropogenic actions, the species also shows a trend of long-term demographic decline that has resulted in approximately an order of magnitude decrease and began 13 thousand years ago.

Island survivors: population genetic structure and demography of the critically endangered giant lizard of La Gomera, Gallotia bravoana

Background

The giant lizard of La Gomera (Gallotia bravoana), is an endemic lacertid of this Canary Island that lives confined to a very restricted area of occupancy in a steep cliff, and is catalogued as Critically Endangered by IUCN. We present the first population genetic analysis of the wild population as well as of captive-born individuals (for which paternity data are available) from a recovery center. Current genetic variability, and inferred past demographic changes were determined in order to discern the relative contribution of natural versus human-mediated effects on the observed decline in population size.

Results

Genetic analyses indicate that the only known natural population of the species shows low genetic diversity and acts as a single evolutionary unit. Demographic analyses inferred a prolonged decline of the species for at least 230 generations. Depending on the assumed generation time, the onset of the decline was dated between 1200–13000 years ago. Pedigree analyses of captive individuals suggest that reproductive behavior of the giant lizard of La Gomera may include polyandry, multiple paternity and female long-term sperm retention.

Conclusions

The current low genetic diversity of G. bravoana is the result of a long-term gradual decline. Because generation time is unknown in this lizard and estimates had large credibility intervals, it is not possible to determine the relative contribution of humans in the collapse of the population. Shorter generation times would favor a stronger influence of human pressure whereas longer generation times would favor a climate-induced origin of the decline. In any case, our analyses show that the wild population has survived for a long period of time with low levels of genetic diversity and a small effective population size. Reproductive behavior may have acted as an important inbreeding avoidance mechanism allowing the species to elude extinction. Overall, our results suggest that the species retains its adaptive potential and could restore its ancient genetic diversity under favorable conditions. Therefore, management of the giant lizard of La Gomera should concentrate efforts on enhancing population growth rates through captive breeding of the species as well as on restoring the carrying capacity of its natural habitat.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-014-0121-8) contains supplementary material, which is available to authorized users.

Genetic status and timing of a weevil introduction to Santa Cruz Island, Galapagos

Successful invasive species can overcome or circumvent the potential genetic loss caused by an introduction bottleneck through a rapid population expansion and admixture from multiple introductions. We explore the genetic makeup and the timing of a species introduction to Santa Cruz Island in the Galápagos archipelago. We investigate the presence of processes that can maintain genetic diversity in populations of the broad-nosed weevil Galapaganus howdenae howdenae. Analyses of combined genotypes for 8 microsatellite loci showed evidence of past population size reductions through moment and likelihood-based estimators. No evidence of admixture through multiple introductions was found, but substantial current population sizes (N0 298, 95% credible limits 50-2300), genetic diversity comparable with long-established endemics (Mean number of alleles = 3.875), and lack of genetic structure across the introduced range (F ST = 0.01359) could suggest that foundations are in place for populations to rapidly recover any loss of genetic variability. The time estimates for the introduction into Santa Cruz support an accidental transfer during the colonization period (1832-1959) predating the spurt in human population growth. Our evaluation of the genetic status of G. h. howdenae suggests potential for population growth in addition to our field observations of a concurrent expansion in range and feeding preferences towards protected areas and endemic host plants.

Genetic consequences of a century of protection: serial founder events and survival of the little spotted kiwi (Apteryx owenii)

We present the outcome of a century of post-bottleneck isolation of a long-lived species, the little spotted kiwi (Apteryx owenii, LSK) and demonstrate that profound genetic consequences can result from protecting few individuals in isolation. LSK were saved from extinction by translocation of five birds from South Island, New Zealand to Kapiti Island 100 years ago. The Kapiti population now numbers some 1200 birds and provides founders for new populations. We used 15 microsatellite loci to compare genetic variation among Kapiti LSK and the populations of Red Mercury, Tiritiri Matangi and Long Islands that were founded with birds from Kapiti. Two LSK native to D'Urville Island were also placed on Long Island. We found extremely low genetic variation and signatures of acute and recent genetic bottleneck effects in all four populations, indicating that LSK have survived multiple genetic bottlenecks. The Long Island population appears to have arisen from a single mating pair from Kapiti, suggesting there is no genetic contribution from D'Urville birds among extant LSK. The Ne/NC ratio of Kapiti Island LSK (0.03) is exceptionally low for terrestrial vertebrates and suggests that genetic diversity might still be eroding in this population, despite its large census size.

Speciation on oceanic islands: rapid adaptive divergence vs. cryptic speciation in a Guadalupe Island songbird (Aves: Junco)

The evolutionary divergence of island populations, and in particular the tempo and relative importance of neutral and selective factors, is of central interest to the study of speciation. The rate of phenotypic evolution upon island colonization can vary greatly among taxa, and cases of convergent evolution can further confound the inference of correct evolutionary histories. Given the potential lability of phenotypic characters, molecular dating of insular lineages analyzed in a phylogenetic framework provides a critical tool to test hypotheses of phenotypic divergence since colonization. The Guadalupe junco is the only insular form of the polymorphic dark-eyed junco (Junco hyemalis), and shares eye and plumage color with continental morphs, yet presents an enlarged bill and reduced body size. Here we use variation in mtDNA sequence, morphological traits and song variables to test whether the Guadalupe junco evolved rapidly following a recent colonization by a mainland form of the dark-eyed junco, or instead represents a well-differentiated "cryptic" lineage adapted to the insular environment through long-term isolation, with plumage coloration a result of evolutionary convergence. We found high mtDNA divergence of the island lineage with respect to both continental J. hyemalis and J. phaeonotus, representing a history of isolation of about 600,000 years. The island lineage was also significantly differentiated in morphological and male song variables. Moreover, and contrary to predictions regarding diversity loss on small oceanic islands, we document relatively high levels of both haplotypic and song-unit diversity on Guadalupe Island despite long-term isolation in a very small geographic area. In contrast to prevailing taxonomy, the Guadalupe junco is an old, well-differentiated evolutionary lineage, whose similarity to mainland juncos in plumage and eye color is due to evolutionary convergence. Our findings confirm the role of remote islands in driving divergence and speciation, but also their potential role as repositories of ancestral diversity.

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.

Molecular identification of traces from the White-tailed Sea Eagle

PURPOSE

Over the preceeding decades, after periods of dramatic decline and extinction in many parts of Europe, the White-tailed Sea Eagle (Haliaeetus albicilla) has re-colonized traditional breeding areas. However, this large apex predator remains threatened, not only by the bioaccumulation of environmental pollutants, but also by targeted poisoning and poaching. In connection with a forensic case, a novel PCR assay was developed for the sensitive and specific detection of sea eagle DNA traces in questioned samples of unknown origin.

METHODS

The assay amplifies a fragment of the popular phylogenetic marker gene cytochrome b. Primers were designed to bind sites with relatively high variability between homologous sequences from H. albicilla and other related European birds of prey.

RESULTS

Assay sensitivity was sufficient for single cell analysis. Specificity was tested in vitro and the primers did not cross-detect DNA from humans, chicken and the following raptors: Common Buzzard (Buteo buteo), Northern Goshawk (Accipiter gentilis), Red Kite (Milvus milvus) and Black Kite (Milvus migrans). Applicability for the analysis of poor quality samples was demonstrated with extracts from field-collected small molted down feathers that did not contain detectable amounts of sea eagle nuclear DNA. Amplicons of the expected size were generated, purified and sequenced. Sequence data were subjected to Basic Local Alignment Search Tool analysis and affiliated with cytochrome b from H. albicilla.

CONCLUSIONS

The novel PCR primers allowed for the correct assignment of traces from H. albicilla, even in mixed samples and in cases with limited and degraded biological material.

Signature of a pre-human population decline in the critically endangered Reunion Island endemic forest bird Coracina newtoni

The exceptional biodiversity of Reunion Island is threatened by anthropogenic landscape changes that took place during the 350 years of human colonization. During this period the human population size increased dramatically from 250 to 800,000. The arrival of humans together with the development of agriculture, invasive species such as rats and cats, and deforestation has lead to the extinction of more than half of the original vertebrate species of the island. For the remaining species, significant work is being carried out to identify threats and conservation status, but little genetic work has been carried on some of the most endangered species. In the last decade theoretical studies have shown the ability of neutral genetic markers to infer the demographic history of endangered species and identify and date past population size changes (expansions or bottlenecks). In this study we provide the first genetic data on the critically endangered species the Reunion cuckoo-shrike Coracina newtoni. The Reunion cuckoo-shrike is a rare endemic forest bird surviving in a restricted 12-km(2) area of forested uplands and mountains. The total known population consists of less than one hundred individuals out of which 45 were genotyped using seventeen polymorphic microsatellite loci. We found a limited level of genetic variability and weak population structure, probably due to the limited geographic distribution. Using Bayesian methods, we identified a strong decline in population size during the Holocene, most likely caused by an ancient climatic or volcanic event around 5000 years ago. This result was surprising as it appeared in apparent contradiction with the accepted theory of recent population collapse due to deforestation and predator introduction. These results suggest that new methods allowing for more complex demographic models are necessary to reconstruct the demographic history of populations.

Extremely low genetic diversity indicating the endangered status of Ranodon sibiricus (Amphibia: Caudata) and implications for phylogeography

BACKGROUND

The Siberian salamander (Ranodon sibiricus), distributed in geographically isolated areas of Central Asia, is an ideal alpine species for studies of conservation and phylogeography. However, there are few data regarding the genetic diversity in R. sibiricus populations.

METHODOLOGY/PRINCIPAL FINDINGS

We used two genetic markers (mtDNA and microsatellites) to survey all six populations of R. sibiricus in China. Both of the markers revealed extreme genetic uniformity among these populations. There were only three haplotypes in the mtDNA, and the overall nucleotide diversity in the mtDNA was 0.00064, ranging from 0.00000 to 0.00091 for the six populations. Although we recovered 70 sequences containing microsatellite repeats, there were only two loci that displayed polymorphism. We used the approximate Bayesian computation (ABC) method to study the demographic history of the populations. This analysis suggested that the extant populations diverged from the ancestral population approximately 120 years ago and that the historical population size was much larger than the present population size; i.e., R. sibiricus has experienced dramatic population declines.

CONCLUSION/SIGNIFICANCE

Our findings suggest that the genetic diversity in the R. sibiricus populations is the lowest among all investigated amphibians. We conclude that the isolation of R. sibiricus populations occurred recently and was a result of recent human activity and/or climatic changes. The Pleistocene glaciation oscillations may have facilitated intraspecies genetic homogeneity rather than enhanced divergence. A low genomic evolutionary rate and elevated inbreeding frequency may have also contributed to the low genetic variation observed in this species. Our findings indicate the urgency of implementing a protection plan for this endangered species.

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.