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

The burden of anthropogenic changes and mutation load in a critically endangered harrier from the Reunion biodiversity hotspot, Circus maillardi

Anthropogenic impact is causing the decline of a large proportion of species worldwide and reduces their genetic diversity. Island species typically have smaller ranges than continental species. As a consequence, island species are particularly liable to undergo population bottlenecks, giving rise to conservation challenges such as inbreeding and unmasking of deleterious genetic load. Such challenges call for more detailed assessments of the genetic make-up of threatened island populations. The Mascarene islands (Indian Ocean) present many prime examples, being unusual in having been pristine until first human arrival ~400 years ago, following which anthropogenic pressure was unusually intense. A threatened harrier (Circus maillardi) endemic to the westernmost island of the archipelago is a good example of the challenges faced by species that have declined to small population size following intense anthropogenic pressure. In this study, we use an extensive set of population genomic tools to quantify variation at near-neutral and coding loci, in order to test the historical impact of human activity on this species, and evaluate the species' (mal)adaptive potential. We observed low but significant genetic differentiation between populations on the West and North-East sides of the island, echoing observations in other endemic species. Inbreeding was significant, with a substantial fraction of samples being first or second-degree relatives. Historical effective population sizes have declined from ~3000 to 300 individuals in the past 1000 years, with a more recent drop ~100 years ago consistent with human activity. Based on our simulations and comparisons with a close relative (Circus melanoleucos), this demographic history may have allowed purging of the most deleterious variants but is unlikely to have allowed the purging of mildly deleterious variants. Our study shows how using relatively affordable methods can reveal the massive impact that human activity may have on the genetic diversity and adaptive potential of island populations, and calls for urgent action to closely monitor the reproductive success of such endemic populations, in association with genetic studies.

Demographic responses of oceanic island birds to local and regional ecological disruptions revealed by whole-genome sequencing

Disentangling the effects of ecological disruptions operating at different spatial and temporal scales in shaping past species' demography is particularly important in the current context of rapid environmental changes driven by both local and regional factors. We argue that volcanic oceanic islands provide useful settings to study the influence of past ecological disruptions operating at local and regional scales on population demographic histories. We investigate potential drivers of past population dynamics for three closely related species of passerine birds from two volcanic oceanic islands, Reunion and Mauritius (Mascarene archipelago), with distinct volcanic history. Using ABC and PSMC inferences from complete genomes, we reconstructed the demographic history of the Reunion Grey White-eye (Zosterops borbonicus (Pennant, 1781)), the Reunion Olive White-eye (Z. olivaceus (Linnaeus, 1766)) and the Mauritius Grey White-eye (Z. mauritianus (Gmelin, 1789)) and searched for possible causes underlying similarities or differences between species living on the same or different islands. Both demographic inferences strongly support ancient and long-term expansions in all species. They also reveal different trajectories between species inhabiting different islands, but consistent demographic trajectories in species or populations from the same island. Species from Reunion appear to have experienced synchronous reductions in population size during the Last Glacial Maximum, a trend not seen in the Mauritian species. Overall, this study suggests that local events may have played a role in shaping population trajectories of these island species. It also highlights the potential of our conceptual framework to disentangle the effects of local and regional drivers on past species' demography and long-term population processes.

Past volcanic activity predisposes an endemic threatened seabird to negative anthropogenic impacts

Humans are regularly cited as the main driver of current biodiversity extinction, but the impact of historic volcanic activity is often overlooked. Pre-human evidence of wildlife abundance and diversity are essential for disentangling anthropogenic impacts from natural events. Réunion Island, with its intense and well-documented volcanic activity, endemic biodiversity, long history of isolation and recent human colonization, provides an opportunity to disentangle these processes. We track past demographic changes of a critically endangered seabird, the Mascarene petrel Pseudobulweria aterrima, using genome-wide SNPs. Coalescent modeling suggested that a large ancestral population underwent a substantial population decline in two distinct phases, ca. 125,000 and 37,000 years ago, coinciding with periods of major eruptions of Piton des Neiges. Subsequently, the ancestral population was fragmented into the two known colonies, ca. 1500 years ago, following eruptions of Piton de la Fournaise. In the last century, both colonies declined significantly due to anthropogenic activities, and although the species was initially considered extinct, it was rediscovered in the 1970s. Our findings suggest that the current conservation status of wildlife on volcanic islands should be firstly assessed as a legacy of historic volcanic activity, and thereafter by the increasing anthropogenic impacts, which may ultimately drive species towards extinction.

Forests Without Frugivores and Frugivores Without Forests – An Investigation Into the Causes of a Paradox in One of the Last Archipelagos Colonized by Humans

The Mascarenes are sadly famous worldwide for the massive extinction of their native vertebrates since recent human colonization. However, extinction patterns show astonishing disparities between the two main islands and between lineages of forest vertebrates. On Réunion (2,512 km2, 3,070 m) where about a third of native habitats remains, most large-bodied vertebrates, especially frugivores, collapsed by the first half of the 18th century, while several have survived longer and some still exist on Mauritius (1,865 km2, 828 m) where more than 95% of native habitats have been transformed. Considering lineages of forest vertebrates shared by both islands (23 genera, 53 species), we test the hypothesis that differing patterns of lowland suitable habitat destruction is the main cause behind this paradox. Before that, we assess the potential impact of other major drivers of extinctions since first contact with humans. Firstly, Mauritius shows earlier and more numerous introductions of mammal predators known for their devastating impact (except northern islets which have thus become important sanctuaries for several squamates). Secondly, settlers were inveterate hunters on both islands, but while Réunion was overhunted before Mauritius, the burst of human population in the latter in late 18th century has not led to the rapid extinction of all large native vertebrates. These two factors alone therefore cannot explain the observed paradox. Rather, the early destruction of lowland habitats (<400 m) on Réunion is concomitant with most extinctions of forest vertebrate, notably frugivores that rapidly lost most lowland habitats dominated by large fleshy-fruited plants. Moreover, landform-induced fragmentation has likely decreased the ability of adjacent habitats to act as effective refuges. Conversely, Mauritius retained suitable low-fragmented habitats until the late 19th which probably allowed, at least for a time, several native vertebrates to escape from multiple human-induced disturbances. Despite the almost total destruction of native habitats since then on Mauritius, conservation actions have saved several threatened vertebrate species that play a fundamental role in the functioning of native ecosystems. The fact that there are now more favorable habitats on Réunion than on Mauritius argues for the rewilding of Réunion with these extant large vertebrates.

Plant population dynamics on oceanic islands during the Late Quaternary climate changes: genetic evidence from a tree species (Coffea mauritiana) in Reunion Island

Past climatic fluctuations have played a major role in shaping the current plant biodiversity. Although harbouring an exceptional biota, oceanic islands have received little attention in studies on species demographic history and past vegetation patterns. We investigated the impact of past climatic changes on the effective population size of a tree (Coffea mauritiana) that is endemic to Reunion Island, located in the south-western Indian Ocean (SWIO). Demographic changes were inferred using summary statistics calculated from genomic data. Using ecological niche modelling and the current distribution of genetic diversity, the paleodistribution of the species was also assessed. A reduction in the effective population size of C. mauritiana during the last glaciation maximum was inferred. The distribution of the species was reduced on the western side of the island, due to low rainfall. It appeared that a major reduction in rainfall and a slight temperature decrease prevailed in the SWIO. Our findings indicated that analyses on the current patterns of intraspecific genetic variations can efficiently contribute to past climatic changes characterisation in remote islands. Identifying area with higher resilience in oceanic islands could provide guidance in forest management and conservation faced to the global climate change.

Genetic signatures of small effective population sizes and demographic declines in an endangered rattlesnake, Sistrurus catenatus

Endangered species that exist in small isolated populations are at elevated risk of losing adaptive variation due to genetic drift. Analyses that estimate short-term effective population sizes, characterize historical demographic processes, and project the trajectory of genetic variation into the future are useful for predicting how levels of genetic diversity may change. Here, we use data from two independent types of genetic markers (single nucleotide polymorphisms [SNPs] and microsatellites) to evaluate genetic diversity in 17 populations spanning the geographic range of the endangered eastern massasauga rattlesnake (Sistrurus catenatus). First, we use SNP data to confirm previous reports that these populations exhibit high levels of genetic structure (overall Fst = 0.25). Second, we show that most populations have contemporary Ne estimates <50. Heterozygosity-fitness correlations in these populations provided no evidence for a genetic cost to living in small populations, though these tests may lack power. Third, model-based demographic analyses of individual populations indicate that all have experienced declines, with the onset of many of these declines occurring over timescales consistent with anthropogenic impacts (<200 years). Finally, forward simulations of the expected loss of variation in relatively large (Ne = 50) and small (Ne = 10) populations indicate they will lose a substantial amount of their current standing neutral variation (63% and 99%, respectively) over the next 100 years. Our results argue that drift has a significant and increasing impact on levels of genetic variation in isolated populations of this snake, and efforts to assess and mitigate associated impacts on adaptive variation should be components of the management of this endangered reptile.

Spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus revealed by museum genomics

Analyzing genetic variation through time and space is important to identify key evolutionary and ecological processes in populations. However, using contemporary genetic data to infer the dynamics of genetic diversity may be at risk of a bias, as inferences are performed from a set of extant populations, setting aside unavailable, rare, or now extinct lineages. Here, we took advantage of new developments in next-generation sequencing to analyze the spatial and temporal genetic dynamics of the grasshopper Oedaleus decorus, a steppic Southwestern-Palearctic species. We applied a recently developed hybridization capture (hyRAD) protocol that allows retrieving orthologous sequences even from degraded DNA characteristic of museum specimens. We identified single nucleotide polymorphisms in 68 historical and 51 modern samples in order to (i) unravel the spatial genetic structure across part of the species distribution and (ii) assess the loss of genetic diversity over the past century in Swiss populations. Our results revealed (i) the presence of three potential glacial refugia spread across the European continent and converging spatially in the Alpine area. In addition, and despite a limited population sample size, our results indicate (ii) a loss of allelic richness in contemporary Swiss populations compared to historical populations, whereas levels of expected heterozygosities were not significantly different. This observation is compatible with an increase in the bottleneck magnitude experienced by central European populations of O. decorus following human-mediated land-use change impacting steppic habitats. Our results confirm that application of hyRAD to museum samples produces valuable information to study genetic processes across time and space.

Microevolution of the noble crayfish (Astacus astacus) in the Southern Balkan Peninsula

Background

The noble crayfish (Astacus astacus) displays a complex historical and contemporary genetic status in Europe. The species divergence has been shaped by geological events (i.e. Pleistocene glaciations) and humanly induced impacts (i.e. translocations, pollution, etc.) on its populations due to species commercial value and its niche degradation. Until now, limited genetic information has been procured for the Balkan area and especially for the southernmost distribution of this species (i.e. Greece). It is well known that the rich habitat diversity of the Balkan Peninsula offers suitable conditions for genetically diversified populations. Thus, the present manuscript revisits the phylogenetic relationships of the noble crayfish in Europe and identifies the genetic make-up and the biogeographical patterns of the species in its southern range limit.

Results

Mitochondrial markers (i.e. COI and 16S) were used in order to elucidate the genetic structure and diversity of the noble crayfish in Europe. Two of the six European haplotypic lineages, were found exclusively in Greece. These two lineages exhibited greater haplotypic richness when compared with the rest four (of “Central European” origin) while they showed high genetic diversity. Divergence time analysis identified that the majority of this divergence was captured through Pleistocene, suggesting a southern glacial refugium (Greece, southern Balkans). Furthermore, six microsatellite markers were used in order to define the factors affecting the genetic structure and demographic history of the species in Greece. The population structure analysis revealed six to nine genetic clusters and eight putative genetic barriers. Evidence of bottleneck effects in the last ~5000 years (due to climatic and geological events and human activities) is also afforded. Findings from several other research fields (e.g. life sciences, geology or even archaeology) have been utilized to perceive the genetic make-up of the noble crayfish.

Conclusions

The southernmost part of Balkans has played a major role as a glacial refugium for A. astacus. Such refugia have served as centres of expansion to northern regions. Recent history of the noble crayfish in southern Balkans reveals the influence of environmental (climate, geology and/or topology) and anthropogenic factors.

Electronic supplementary material

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

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

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

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

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

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.

Mixed fortunes: ancient expansion and recent decline in population size of a subtropical montane primate, the Arunachal macaque Macaca munzala

Quaternary glacial oscillations are known to have caused population size fluctuations in many temperate species. Species from subtropical and tropical regions are, however, considerably less studied, despite representing most of the biodiversity hotspots in the world including many highly threatened by anthropogenic activities such as hunting. These regions, consequently, pose a significant knowledge gap in terms of how their fauna have typically responded to past climatic changes. We studied an endangered primate, the Arunachal macaque Macaca munzala, from the subtropical southern edge of the Tibetan plateau, a part of the Eastern Himalaya biodiversity hotspot, also known to be highly threatened due to rampant hunting. We employed a 534 bp-long mitochondrial DNA sequence and 22 autosomal microsatellite loci to investigate the factors that have potentially shaped the demographic history of the species. Analysing the genetic data with traditional statistical methods and advance Bayesian inferential approaches, we demonstrate a limited effect of past glacial fluctuations on the demographic history of the species before the last glacial maximum, approximately 20,000 years ago. This was, however, immediately followed by a significant population expansion possibly due to warmer climatic conditions, approximately 15,000 years ago. These changes may thus represent an apparent balance between that displayed by the relatively climatically stable tropics and those of the more severe, temperate environments of the past. This study also draws attention to the possibility that a cold-tolerant species like the Arunachal macaque, which could withstand historical climate fluctuations and grow once the climate became conducive, may actually be extremely vulnerable to anthropogenic exploitation, as is perhaps indicated by its Holocene ca. 30-fold population decline, approximately 3,500 years ago. Our study thus provides a quantitative appraisal of these demographically important events, emphasising the ability to potentially infer the occurrence of two separate historical events from contemporary genetic data.

Pyrenean ptarmigans decline under climatic and human influences through the Holocene

In Europe, the Quaternary is characterized by climatic fluctuations known to have led to many cycles of contraction and expansion of species geographical ranges. In addition, during the Holocene, historical changes in human occupation such as colonization or abandonment of traditional land uses can also affect habitats. These climatically or anthropically induced geographic range changes are expected to produce considerable effective population size change, measurable in terms of genetic diversity and organization. The rock ptarmigan (Lagopus muta) is a small-bodied grouse occurring throughout Northern hemispheric arctic and alpine tundra. This species is not considered threatened at a continental scale, but the populations in the Pyrenees are of concern because of their small population size, geographical isolation and low genetic diversity. Here, we used 11 microsatellites to investigate genetic variations and differentiations and infer the overall demographic history of Pyrenean rock ptarmigan populations. The low genetic variability found in these populations has been previously thought to be the result of a bottleneck that occurred following the last glacial maximum (i.e., 10,000 years ago) or more recently (i.e., during the last 200 years). Our results clearly indicate a major bottleneck affecting the populations in the last tenth of the Holocene. We discuss how this decline can be explained by a combination of unfavorable and successive events that increased the degree of habitat fragmentation.

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

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

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

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