Genetic structure and signature of population decrease in the critically endangered freshwater cyprinid Chondrostoma lusitanicum

River capture or ancestral polymorphism: an empirical genetic test in a freshwater fish using approximate Bayesian computation

A headwater or river capture is a phenomenon commonly invoked to explain the absence of reciprocal monophyly of genetic lineages among isolated hydrographic basins in freshwater fish. Under the assumption of river capture, a secondary contact between populations previously isolated in different basins explains the observed genetic pattern. However, the absence of reciprocal monophyly could also arise under population isolation through the retention of ancestral of polymorphisms. Here, we applied an approximate Bayesian computation (ABC) framework for estimating the relative probability of scenarios with and without secondary contact. We used Cnesterodon decemmaculatus as a study model because of the multiple possible cases of river capture and the demographic parameters estimated in a previous mitochondrial DNA study that are useful for simulating scenarios to test both hypotheses using the ABC framework. Our results showed that, in general, mitochondrial DNA is useful for distinguishing between these alternative demographic scenarios with reasonable confidence, but in extreme cases (e.g. recent divergence or large population size) there is no power to discriminate between scenarios. Testing hypotheses of drainage rearrangement under a statistically rigorous framework is fundamental for understanding the evolution of freshwater fish fauna as a complement to, or in the absence of, geological evidence.

Fish diversification at the pace of geomorphological changes: evolutionary history of western Iberian Leuciscinae (Teleostei: Leuciscidae) inferred from multilocus sequence data

The evolutionary history of western Iberian Leuciscinae, obligatory freshwater fish, is directly linked to the evolution of the hydrographic network of the Iberian Peninsula after its isolation from the rest of Europe, which involved dramatic rearrangements such as the transition from endorheic lakes to open basins draining to the Atlantic. Previous phylogenetic research on western Iberian leuciscines, using mainly mitochondrial DNA and more recently one or two nuclear genes, has found contradictory results and there remain many unresolved issues regarding species relationships, taxonomy, and evolutionary history. Moreover, there is a lack of integration between phylogenetic and divergence time estimates and information on the timing of geomorphological changes and paleobasin rearrangements in the Iberian Peninsula. This study presents the first comprehensive fossil-calibrated multilocus coalescent species tree of western Iberian Leuciscinae (including 14 species of Achondrostoma, Iberochondrostoma, Pseudochondrostoma and Squalius endemic to the Iberian Peninsula, seven of which endemic to Portugal) based on seven nuclear genes, and from which we infer their biogeographic history by comparing divergence time estimates to known dated geological events. The phylogenetic pattern suggests slow-paced evolution of leuciscines during the Early-Middle Miocene endorheic phase of the main Iberian river basins, with the shift to exorheism in the late Neogene-Quaternary allowing westward dispersals that resulted in many cladogenetic events and a high rate of endemism in western Iberia. The results of this study also: (i) confirm the paraphyly of S. pyrenaicus with respect to S. carolitertii, and thus the possible presence of a new taxon in the Portuguese Tagus currently assigned to S. pyrenaicus; (ii) support the taxonomic separation of the Guadiana and Sado populations of S. pyrenaicus; (iii) show the need for further population sampling and taxonomic research to clarify the phylogenetic status of A. arcasii from the Minho basin and of the I. lusitanicum populations in the Sado and Tagus basins; and (iv) indicate that A. occidentale, I. olisiponensis and P. duriensis are the most ancient lineages within their respective genera.

Broad-scale sampling of primary freshwater fish populations reveals the role of intrinsic traits, inter-basin connectivity, drainage area and latitude on shaping contemporary patterns of genetic diversity

Background. Worldwide predictions suggest that up to 75% of the freshwater fish species occurring in rivers with reduced discharge could be extinct by 2070 due to the combined effect of climate change and water abstraction. The Mediterranean region is considered to be a hotspot of freshwater fish diversity but also one of the regions where the effects of climate change will be more severe. Iberian cyprinids are currently highly endangered, with over 68% of the species raising some level of conservation concern.

Methods. During the FISHATLAS project, the Portuguese hydrographical network was extensively covered (all the 34 river basins and 47 sub-basins) in order to contribute with valuable data on the genetic diversity distribution patterns of native cyprinid species. A total of 188 populations belonging to 16 cyprinid species of Squalius, Luciobarbus, Achondrostoma, Iberochondrostoma, Anaecypris and Pseudochondrostoma were characterized, for a total of 3,678 cytochrome b gene sequences.

Results. When the genetic diversity of these populations was mapped, it highlighted differences among populations from the same species and between species with identical distribution areas. Factors shaping the contemporary patterns of genetic diversity were explored and the results revealed the role of latitude, inter-basin connectivity, migratory behaviour, species maximum size, species range and other species intrinsic traits in determining the genetic diversity of sampled populations. Contrastingly, drainage area and hydrological regime (permanent vs. temporary) seem to have no significant effect on genetic diversity. Species intrinsic traits, maximum size attained, inter-basin connectivity and latitude explained over 30% of the haplotype diversity variance and, generally, the levels of diversity were significantly higher for smaller sized species, from connected and southerly river basins.

Discussion. Targeting multiple co-distributed species of primary freshwater fish allowed us to assess the relative role of historical versus contemporary factors affecting genetic diversity. Since different patterns were detected for species with identical distribution areas we postulate that contemporary determinants of genetic diversity (species’ intrinsic traits and landscape features) must have played a more significant role than historical factors. Implications for conservation in a context of climate change and highly disturbed habitats are detailed, namely the need to focus management and conservation actions on intraspecific genetic data and to frequently conduct combined genetic and demographic surveys.

Cenozoic tectonic and climatic events in southern Iberian Peninsula: Implications for the evolutionary history of freshwater fish of the genus Squalius (Actinopterygii, Cyprinidae)

Southern Iberian freshwater ecosystems located at the border between the European and African plates represent a tectonically complex region spanning several geological ages, from the uplifting of the Betic Mountains in the Serravalian-Tortonian periods to the present. This area has also been subjected to the influence of changing climate conditions since the Middle-Upper Pliocene when seasonal weather patterns were established. Consequently, the ichthyofauna of southern Iberia is an interesting model system for analyzing the influence of Cenozoic tectonic and climatic events on its evolutionary history. The cyprinids Squalius malacitanus and Squalius pyrenaicus are allopatrically distributed in southern Iberia and their evolutionary history may have been defined by Cenozoic tectonic and climatic events. We analyzed MT-CYB (510 specimens) and RAG1 (140 specimens) genes of both species to reconstruct phylogenetic relationships and to estimate divergence times and ancestral distribution ranges of the species and their populations. We also assessed their levels of genetic structure and diversity as well as the amount of gene flow between populations. To investigate recent paleogeographical and climatic factors in southern Iberia, we modeled changes-through-time in sea level from the LGM to the present. Phylogenetic, geographic and population structure analyses revealed two well-supported species (S. malacitanus and S. pyrenaicus) in southern Iberia and two subclades (Atlantic and Mediterranean) within S. malacitanus. The origin of S. malacitanus and the separation of its Atlantic and Mediterranean populations occurred during the Serravalian-Tortonian and Miocene-Pliocene periods, respectively. These divergence events occurred in the Middle Pliocene and Pleistocene in S. pyrenaicus. In both species, Atlantic basins possessed populations with higher genetic diversity than Mediterranean, which may be explained by the Janda Lagoon. The isolation of S. malacitanus was earlier and related to the rising of the Betic Mountains. Divergence of its Atlantic and Mediterranean populations was associated with the creation of the freshwater systems of southern Iberia close to the Gibraltar Strait. The presence of S. pyrenaicus in southern Iberia may be the result of recent colonization associated with river capture, as demonstrated our biogeographic reconstruction.

Metapopulations in temporary streams - the role of drought-flood cycles in promoting high genetic diversity in a critically endangered freshwater fish and its consequences for the future

Genetic factors have direct and indirect impacts in the viability of endangered species. Assessing their genetic diversity levels and population structure is thus fundamental for conservation and management. In this paper we use mitochondrial and nuclear markers to address phylogeographic and demographic data on the critically endangered Anaecypris hispanica, using a broad sampling set which covered its known distribution area in the Iberian Peninsula. Our results showed that the populations of A. hispanica are strongly differentiated (high and significant ФST and FST values, corroborated by the results from AMOVA and SAMOVA) and genetically diversified. We suggest that the restricted gene flow between populations may have been potentiated by ecological, hydrological and anthropogenic causes. Bayesian skyline plots revealed a signal for expansion for all populations (tMRCA between 68kya and 1.33Mya) and a genetic diversity latitudinal gradient was detected between the populations from the Upper (more diversified) and the Lower (less diversified) Guadiana river basin. We postulate a Pleistocenic westwards colonization route for A. hispanica in the Guadiana river basin, which is in agreement with the tempo and mode of paleoevolution of this drainage. The colonization of River Guadalquivir around 60kya with migrants from the Upper Guadiana, most likely by stream capture, is also suggested. This study highlights the view that critically endangered species facing range retreats (about 47% of its known populations have disappeared in the last 15years) are not necessarily small and genetically depleted. However, the extinction risk is not negligible since A. hispanica faces the combined effect of several deterministic and stochastic negative factors and, moreover, recolonization events after localized extinctions are very unlikely to occur due to the strong isolation of populations and to the patchily ecologically-conditioned distribution of fish. The inferred species distribution models highlight the significant contribution of temperature seasonality and isothermality to A. hispanica occurrence in Guadiana environments and emphasize the importance of stable climatic conditions for the preservation of this species. Given the strong population structure, high percentage of private haplotypes and virtual absence of inter-basin gene flow we suggest that each A. hispanica population should be considered as an independent Operational Conservation Unit and that ex-situ and in-situ actions should be conducted in parallel to allow for the long-term survival of the species and the preservation of the genetic integrity of its populations.

Temporal population genetics of time travelling insects: a long term study in a seed-specialized wasp

Many animal species experiencing spatial or interannual fluctuations of their environment are capable of prolonged diapause, a kind of dormancy that extends over more than one year. Such a prolonged diapause is commonly perceived as a temporal demographic refuge in stochastic environments, but empirical evidence is still lacking of its consequences on temporal population genetic structures. In this long-term study, we investigated how a particular pattern of prolonged diapause may influence the temporal population genetics of the invasive seed-specialized wasp Megastigmus schimitscheki (Hymenoptera: Torymidae) in southeastern France. We characterized the diapause strategy of M. schimitscheki using records of emergence from diapause in 97 larval cohorts, and we conducted a temporal population genetic study on a natural invasive wasp population sampled during ten consecutive years (1999-2008) using polymorphic microsatellite markers. We found that M. schimitscheki can undergo a prolonged diapause of up to five years and displays two main adult emergence peaks after two and four years of diapause. Such a bimodal and atypical pattern did not disrupt temporal gene flow between cohorts produced in even and in odd years during the period of the study. Unexpectedly, we found that this wasp population consisted of two distinct genetic sub-populations that strongly diverged in their diapause strategies, with very few admixed individuals. One of the sub-populations displayed both short and prolonged diapause (2 and 4 years respectively) in equal proportions, whereas the other sub-population displayed mainly short diapause. This study provided empirical evidence that prolonged diapause phenotypes can substantially contribute to reproduction and impact temporal genetic structures. Prolonged diapause is likely to act as both demographic and genetic refuges for insect populations living in fluctuating environments.

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.

Genetic structure of the tree peony (Paeonia rockii) and the Qinling Mountains as a geographic barrier driving the fragmentation of a large population

Background

Tree peonies are great ornamental plants associated with a rich ethnobotanical history in Chinese culture and have recently been used as an evolutionary model. The Qinling Mountains represent a significant geographic barrier in Asia, dividing mainland China into northern (temperate) and southern (semi–tropical) regions; however, their flora has not been well analyzed. In this study, the genetic differentiation and genetic structure of Paeonia rockii and the role of the Qinling Mountains as a barrier that has driven intraspecific fragmentation were evaluated using 14 microsatellite markers.

Methodology/Principal Findings

Twenty wild populations were sampled from the distributional range of P. rockii. Significant population differentiation was suggested (F_ST value of 0.302). Moderate genetic diversity at the population level (H_S of 0.516) and high population diversity at the species level (H_T of 0.749) were detected. Significant excess homozygosity (F_IS of 0.076) and recent population bottlenecks were detected in three populations. Bayesian clusters, population genetic trees and principal coordinate analysis all classified the P. rockii populations into three genetic groups and one admixed Wenxian population. An isolation-by-distance model for P. rockii was suggested by Mantel tests (r = 0.6074, P<0.001) and supported by AMOVA (P<0.001), revealing a significant molecular variance among the groups (11.32%) and their populations (21.22%). These data support the five geographic boundaries surrounding the Qinling Mountains and adjacent areas that were detected with Monmonier's maximum-difference algorithm.

Conclusions/Significance

Our data suggest that the current genetic structure of P. rockii has resulted from the fragmentation of a formerly continuously distributed large population following the restriction of gene flow between populations of this species by the Qinling Mountains. This study provides a fundamental genetic profile for the conservation and responsible exploitation of the extant germplasm of this species and for improving the genetic basis for breeding its cultivars.

Population divergence with or without admixture: selecting models using an ABC approach

Genetic data have been widely used to reconstruct the demographic history of populations, including the estimation of migration rates, divergence times and relative admixture contribution from different populations. Recently, increasing interest has been given to the ability of genetic data to distinguish alternative models. One of the issues that has plagued this kind of inference is that ancestral shared polymorphism is often difficult to separate from admixture or gene flow. Here, we applied an approximate Bayesian computation (ABC) approach to select the model that best fits microsatellite data among alternative splitting and admixture models. We performed a simulation study and showed that with reasonably large data sets (20 loci) it is possible to identify with a high level of accuracy the model that generated the data. This suggests that it is possible to distinguish genetic patterns due to past admixture events from those due to shared polymorphism (population split without admixture). We then apply this approach to microsatellite data from an endangered and endemic Iberian freshwater fish species, in which a clustering analysis suggested that one of the populations could be admixed. In contrast, our results suggest that the observed genetic patterns are better explained by a population split model without admixture.

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 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.

Karyotype and genome size of Iberochondrostoma almacai (Teleostei, Cyprinidae) and comparison with the sister-species I.lusitanicum

This study aimed to define the karyotype of the recently described Iberian endemic Iberochondrostoma almacai, to revisit the previously documented chromosome polymorphisms of its sister species I.lusitanicum using C-, Ag-/CMA₃ and RE-banding, and to compare the two species genome sizes. A 2n = 50 karyotype (with the exception of a triploid I.lusitanicum specimen) and a corresponding haploid chromosome formula of 7M:15SM:3A (FN = 94) were found. Multiple NORs were observed in both species (in two submetacentric chromosome pairs, one of them clearly homologous) and a higher intra and interpopulational variability was evidenced in I.lusitanicum. Flow cytometry measurements of nuclear DNA content showed some significant differences in genome size both between and within species: the genome of I. almacai was smaller than that of I.lusitanicum (mean values 2.61 and 2.93 pg, respectively), which presented a clear interpopulational variability (mean values ranging from 2.72 to 3.00 pg). These data allowed the distinction of both taxa and confirmed the existence of two well differentiated groups within I. lusitanicum: one that includes the populations from the right bank of the Tejo and Samarra drainages, and another that reunites the southern populations. The peculiar differences between the two species, presently listed as “Critically Endangered”, reinforced the importance of this study for future conservation plans.