What remains from a 454 run: estimation of success rates of microsatellite loci development in selected newt species (Calotriton asper, Lissotriton helveticus, and Triturus cristatus) and comparison with Illumina-based approaches

Context-dependent dispersal determines relatedness and genetic structure in a patchy amphibian population

Dispersal is a central process in ecology and evolution with far reaching consequences for the dynamics and genetics of spatially structured populations (SSPs). Individuals can adjust their decisions to disperse according to local fitness prospects, resulting in context-dependent dispersal. By determining dispersal rate, distance and direction, these individual-level decisions further modulate the demography, relatedness and genetic structure of SSPs. Here, we examined how context-dependent dispersal influences the dynamics and genetics of a great crested newt (Triturus cristatus) SSP. We collected capture-recapture data of 5564 individuals and genetic data of 950 individuals across an SSP in northern Germany. We added genetic data from six sites outside this SSP to assess genetic structure and gene flow at a regional level. Dispersal rates within the SSP were high but dispersal distances were short. Dispersal was context-dependent: individuals preferentially immigrated into high-quality ponds where breeding probabilities were higher. The studied SSP behaved like a patchy population, where subpopulations at each pond were demographically interdependent. High context-dependent dispersal led to weak but significant spatial genetic structure and relatedness within the SSP. At the regional level, a strong hierarchical genetic structure with very few first-generation migrants as well as low effective dispersal rates suggest the presence of independent demographic units. Overall, our study highlights the importance of habitat quality for driving context-dependent dispersal and therefore demography and genetic structure in SSPs. Limited capacity for long-distance dispersal seems to increase genetic structure within a population and leads to demographic isolation in anthropogenic landscapes.

Isolation and characterization of microsatellites from a cicada, Yezoterpnosia nigricosta (Hemiptera: Cicadidae), distributed in subarctic and cool temperate forests

The cicada Yezoterpnosia nigricosta (Hemiptera: Cicadidae) is distributed in subarctic and cool temperate forests in Japan, China and the Russian Far East. Due to its limited and isolated distribution at higher altitudes in cool temperate forests on the main island of Japan, especially in the central and southern areas, this species is listed as a threatened or near-threatened species on the Red List of 12 prefectures in Japan. Moreover, there are concerns about the impacts of climate change on the species' distribution and population demography. In this study, seventeen microsatellite markers were developed for Y. nigricosta, and marker suitability was evaluated using 32 individuals from two populations in Nagano prefecture (central Japan) and Hokkaido, a northern island of Japan. The number of alleles, expected heterozygosity and fixation index at each locus were 1-15 (mean = 4.294), 0.000-0.914 (mean = 0.519) and -0.225-0.456 (mean = 0.108), respectively. Furthermore, there was moderate genetic differentiation between the two populations (F_ST = 0.111, F'_ST = 0.237). These markers will be useful to evaluate the genetic structure and to infer population demographic history of Y. nigricosta populations, which can contribute to population genetics studies of this species.

Multiple glacial refugia and contemporary dispersal shape the genetic structure of an endemic amphibian from the Pyrenees

Historical factors (colonization scenarios, demographic oscillations) and contemporary processes (population connectivity, current population size) largely contribute to shaping species' present-day genetic diversity and structure. In this study, we use a combination of mitochondrial and nuclear DNA markers to understand the role of Quaternary climatic oscillations and present-day gene flow dynamics in determining the genetic diversity and structure of the newt Calotriton asper (Al. Dugès, 1852), endemic to the Pyrenees. Mitochondrial DNA did not show a clear phylogeographic pattern and presented low levels of variation. In contrast, microsatellites revealed five major genetic lineages with admixture patterns at their boundaries. Approximate Bayesian computation analyses and linear models indicated that the five lineages likely underwent separate evolutionary histories and can be tracked back to distinct glacial refugia. Lineage differentiation started around the Last Glacial Maximum at three focal areas (western, central and eastern Pyrenees) and extended through the end of the Last Glacial Period in the central Pyrenees, where it led to the formation of two more lineages. Our data revealed no evidence of recent dispersal between lineages, whereas borders likely represent zones of secondary contact following expansion from multiple refugia. Finally, we did not find genetic evidence of sex-biased dispersal. This work highlights the importance of integrating past evolutionary processes and present-day gene flow and dispersal dynamics, together with multilocus approaches, to gain insights into what shaped the current genetic attributes of amphibians living in montane habitats.

Landscape genetics of northern crested newt Triturus cristatus populations in a contrasting natural and human-impacted boreal forest

Among vertebrates, amphibians currently have the highest proportion of threatened species worldwide, mainly through loss of habitat, leading to increased population isolation. Smaller amphibian populations may lose more genetic diversity, and become more dependent on immigration for survival. Investigations of landscape factors and patterns mediating migration and population genetic differentiation are fundamental for knowledge-based conservation. The pond-breeding northern crested newt (Triturus cristatus) populations are decreasing throughout Europe, and are a conservation concern. Using microsatellites, we studied the genetic structure of the northern crested newt in a boreal forest ecosystem containing two contrasting landscapes, one subject to recent change and habitat loss by clear-cutting and roadbuilding, and one with little anthropogenic disturbance. Newts from 12 breeding ponds were analyzed for 13 microsatellites and 7 landscape and spatial variables. With a Maximum-likelihood population-effects model we investigated important landscape factors potentially explaining genetic patterns. Results indicate that intervening landscape factors between breeding ponds, explain the genetic differentiation in addition to an isolation-by-distance effect. Geographic distance, gravel roads, and south/south-west facing slopes reduced landscape permeability and increased genetic differentiation for these newts. The effect was opposite for streams, presumably being more favorable for newt dispersal. Populations within or bordering on old growth forest had a higher allelic richness than populations in managed forest outside these areas. Old growth forest areas may be important source habitats in the conservation of northern crested newt populations.

Microsatellite development from genome skimming and transcriptome sequencing: comparison of strategies and lessons from frog species

Background

Even though microsatellite loci frequently have been isolated using recently developed next-generation sequencing (NGS) techniques, this task is still difficult because of the subsequent polymorphism screening requires a substantial amount of time. Selecting appropriate polymorphic microsatellites is a critical issue for ecological and evolutionary studies. However, the extent to which assembly strategy, read length, sequencing depth, and library layout produce a measurable effect on microsatellite marker development remains unclear. Here, we use six frog species for genome skimming and two frog species for transcriptome sequencing to develop microsatellite markers, and investigate the effect of different isolation strategies on the yield of microsatellites.

Results

The results revealed that the number of isolated microsatellites increases with increased data quantity and read length. Assembly strategy could influence the yield and the polymorphism of microsatellite development. Larger k-mer sizes produced fewer total number of microsatellite loci, but these loci had a longer repeat length, suggesting greater polymorphism. However, the proportion of each type of nucleotide repeats was not affected; dinucleotide repeats were always the dominant type. Finally, the transcriptomic microsatellites displayed lower levels of polymorphisms and were less abundant than genomic microsatellites, but more likely to be functionally linked loci.

Conclusions

These observations provide deep insight into the evolution and distribution of microsatellites and how different isolation strategies affect microsatellite development using NGS.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5329-y) contains supplementary material, which is available to authorized users.

Jailed in the mountains: Genetic diversity and structure of an endemic newt species across the Pyrenees

The Pyrenees represent a natural laboratory for biogeographic, evolutionary and ecological research of mountain fauna as a result of the high variety of habitats and the profound effect of the glacial and interglacial periods. There is a paucity of studies providing a detailed insight into genetic processes and better knowledge on the patterns of genetic diversity and how they are maintained under high altitude conditions. This is of particular interest when considering the course of past climate conditions and glaciations in a species which is considered site tenacious, with long generation times. Here we analyzed the genetic patterns of diversity and structure of the endemic Pyrenean brook newt (Calotriton asper) along its distribution range, with special emphasis on the distinct habitat types (caves, streams, and lakes), and the altitudinal and geographical ranges, using a total set of 900 individuals from 44 different localities across the Pyrenean mountain range genotyped for 19 microsatellite loci. We found evidence for a negative longitudinal and positive altitudinal gradient of genetic diversity in C. asper populations. The fact that genetic diversity was markedly higher westwards is in accordance with other Pyrenean species. However, the impact of altitudinal gradient on the genetic diversity seems to differ from other species, and mostly from other amphibians. We found that lower altitudes can act as a barrier probably because the lowlands do not provide a suitable habitat for C. asper. Regarding the distinct habitat types, caves had significantly lower values of genetic diversity compared to streams or lakes. The mean F_ST value was relatively high (0.304) with maximum values as high as 0.771, suggesting a highly structured total population. Indeed, populations were grouped into five subclusters, the eastern populations (cluster 1) remained grouped into two subclusters and the central-western Pyrenees (cluster 2) into three subclusters. The increase of isolation with geographical distance is consistent with the population structure detected. In conclusion, C. asper seems to be adapted to high altitude mountain habitats, and its genetic diversity is higher in the western Pyrenees. In terms of conservation priority, we consider more relevant the populations that represent a reservoir of genetic diversity.

Getting off to a good start? Genetic evaluation of the ex situ conservation project of the Critically Endangered Montseny brook newt (Calotriton arnoldi)

Ex situ management strategies play an important role in the conservation of threatened species when the wild survival of the species cannot be ensured. Molecular markers have become an outstanding tool for the evaluation and management of captive breeding programs. Two main genetic objectives should be prioritized when planning breeding programs: the maintenance of maximum neutral genetic diversity, and to obtain “self-sustaining” captive populations. In this study, we use 24 microsatellite loci to analyze and evaluate the genetic representativity of the initial phases of the captive breeding program of the Montseny brook newt, Calotriton arnoldi, an Iberian endemic listed as Critically Endangered. The results show that the initial captive stock has 74–78% of the alleles present in the wild populations, and captures roughly 93–95% of their total genetic diversity as observed in a previous study on wild newts, although it does not reach the desired 97.5%. Moreover, the percentage of unrelatedness among individuals does not exceed 95%. Therefore, we conclude that the genetic diversity of the captive stock should be improved by incorporating genetic material from unrelated wild newts. In recognition of the previously described significant genetic and morphological differentiation between eastern and western wild populations of C. arnoldi, we suggest maintaining two distinct breeding lines, and we do not recommend outbreeding between these lines. Our comparisons of genetic diversity estimates between real and distinct sample-sized simulated populations corroborated that a minimum of 20 individuals are needed for each captive population, in order to match the level of genetic diversity present in the wild populations. Thus, the current initial stock should be reinforced by adding wild specimens. The captive stock and subsequent cohorts should be monitored in order to preserve genetic variation. In order to avoid genetic adaptation to captivity, occasionally incorporating previously genotyped individuals from the wild into the captive populations is recommended.

Islands within an island: Population genetic structure of the endemic Sardinian newt, Euproctus platycephalus

The identification of historic and contemporary barriers to dispersal is central to the conservation of endangered amphibians, but may be hindered by their complex life history and elusive nature. The complementary information generated by mitochondrial (mtDNA) and microsatellite markers generates a valuable tool in elucidating population structure and the impact of habitat fragmentation. We applied this approach to the study of an endangered montane newt, Euproctus platycephalus. Endemic to the Mediterranean island of Sardinia, it is threatened by anthropogenic activity, disease, and climate change. We have demonstrated a clear hierarchy of structure across genetically divergent and spatially distinct subpopulations. Divergence between three main mountain regions dominated genetic partitioning with both markers. Mitochondrial phylogeography revealed a deep division dating to ca. 1 million years ago (Mya), isolating the northern region, and further differentiation between the central and southern regions ca. 0.5 Mya, suggesting an association with Pleistocene severe glacial oscillations. Our findings are consistent with a model of southward range expansion during glacial periods, with postglacial range retraction to montane habitat and subsequent genetic isolation. Microsatellite markers revealed further strong population structure, demonstrating significant divergence within the central region, and partial differentiation within the south. The northern population showed reduced genetic diversity. Discordance between mitochondrial and microsatellite markers at this scale indicated a further complexity of population structure, in keeping with male-biased dispersal and female philopatry. Our study underscores the need to elucidate cryptic population structure in the ecology and conservation strategies for endangered island-restricted amphibians, especially in the context of disease and climate change.

No signs of inbreeding despite long-term isolation and habitat fragmentation in the critically endangered Montseny brook newt (Calotriton arnoldi)

Endemic species with restricted geographic ranges potentially suffer the highest risk of extinction. If these species are further fragmented into genetically isolated subpopulations, the risk of extinction is elevated. Habitat fragmentation is generally considered to have negative effects on species survival, despite some evidence for neutral or even positive effects. Typically, non-negative effects are ignored by conservation biology. The Montseny brook newt (Calotriton arnoldi) has one of the smallest distribution ranges of any European amphibian (8 km²) and is considered critically endangered by the International Union for Conservation of Nature. Here we apply molecular markers to analyze its population structure and find that habitat fragmentation owing to a natural barrier has resulted in strong genetic division of populations into two sectors, with no detectable migration between sites. Although effective population size estimates suggest low values for all populations, we found low levels of inbreeding and relatedness between individuals within populations. Moreover, C. arnoldi displays similar levels of genetic diversity to its sister species Calotriton asper, from which it separated around 1.5 million years ago and which has a much larger distribution range. Our extensive study shows that natural habitat fragmentation does not result in negative genetic effects, such as the loss of genetic diversity and inbreeding on an evolutionary timescale. We hypothesize that species in such conditions may evolve strategies (for example, special mating preferences) to mitigate the effects of small population sizes. However, it should be stressed that the influence of natural habitat fragmentation on an evolutionary timescale should not be conflated with anthropogenic habitat loss or degradation when considering conservation strategies.

High gene flow between alternative morphs and the evolutionary persistence of facultative paedomorphosis

Paedomorphosis and metamorphosis are two major developmental processes that characterize the evolution of complex life cycles in many lineages. Whereas these processes were fixed in some taxa, they remained facultative in others, with alternative phenotypes expressed in the same populations. From a genetic perspective, it is still unknown whether such phenotypes form a single population or whether they show some patterns of isolation in syntopy. This has deep implications for understanding the evolution of the phenotypes, i.e. towards their persistence or their fixation and speciation. Newts and salamanders are excellent models to test this hypothesis because they exhibit both developmental processes in their populations: the aquatic paedomorphs retain gills, whereas the metamorphs are able to colonize land. Using microsatellite data of coexisting paedomorphic and metamorphic palmate newts (Lissotriton helveticus), we found that they formed a panmictic population, which evidences sexual compatibility between the two phenotypes. The high gene flow could be understood as an adaptation to unstable habitats in which phenotypic plasticity is favored over the fixation of developmental alternatives. This makes then possible the persistence of a polyphenism: only metamorphosis could be maintained in case of occasional drying whereas paedomorphosis could offer specific advantages in organisms remaining in water.

Using kernels and ecological niche modeling to delineate conservation areas in an endangered patch-breeding phenotype

Efficient delineation of conservation areas is a great challenge in maintaining biodiversity. Kernel density estimators (KDEs) are a powerful tool in this perspective, but they have not been applied at the population level on patch-distributed organisms. This would be particularly worthy for species that need broad habitats beyond those where they can be sampled; such as terrestrial lands for pond-breeding amphibians. The aim of this study was to compare different approaches for the identification of suitable areas for conservation: KDE, ecological niche modelling, and a combination of KDE and niche models. Paedomorphosis was chosen as a model system because this is an important form of intraspecific variation that is present in numerous taxa, but geographically localized within species and globally endangered. 277 ponds were sampled in one of the hotspots of paedomorphosis to determine the abundance and distribution of paedomorphs (i.e., individuals retaining gills at the adult stage) of the palmate newt (Lissotriton helveticus), with emphasis on the connections between the most valuable populations. KDEs gave insights into the surface areas required to balance the maintenance of certain number of connected ponds and the respective number of disjoint areas in which the whole population is divided. The inclusion of barriers in the models helped in accurately designing the limits of the areas to protect. Alone, habitat models were not able to successfully delineate the area to protect, but the integration between terrestrial suitable areas or barriers and KDE allowed an objective identification of areas required for conservation. Overall, the best performance was observed by the KDE integrating ecological barriers, and by the combination between KDE and niche modelling. In a broader perspective, KDEs are thus a pertinent tool in providing quantitative spatial measurements to delineate conservation areas based on patch-abundance data with a specific focus to connectivity.

Isolation and characterization of new microsatellite markers for application in population genetic studies of Caiman latirostris and related species

Wild populations ofCaiman latirostrisare subject to sustainable use programs in Argentina, becoming a species with important impact in the regional economy, based in their skin and meat. Genetic studies are fundamental to acquire information on important parameters for conservation and management, which may be obtained from analysis of molecular markers. Some microsatellites have been previously isolated in this species, but due to some difficulties in using them, we obtained new ones using Next Generation Sequencing approach. This study reports eight new microsatellites forC. latirostrisand tests their utility in a related species,Caiman yacare, with successful application in population genetics and mating systems studies. In addition, we shared data about a novel and fast bioinformatics tool to find microsatellites and to design their corresponding primers.

Genetic fingerprinting proves cross-correlated automatic photo-identification of individuals as highly efficient in large capture-mark-recapture studies

Capture-mark-recapture (CMR) approaches are the backbone of many studies in population ecology to gain insight on the life cycle, migration, habitat use, and demography of target species. The reliable and repeatable recognition of an individual throughout its lifetime is the basic requirement of a CMR study. Although invasive techniques are available to mark individuals permanently, noninvasive methods for individual recognition mainly rest on photographic identification of external body markings, which are unique at the individual level. The re-identification of an individual based on comparing shape patterns of photographs by eye is commonly used. Automated processes for photographic re-identification have been recently established, but their performance in large datasets (i.e., > 1000 individuals) has rarely been tested thoroughly. Here, we evaluated the performance of the program AMPHIDENT, an automatic algorithm to identify individuals on the basis of ventral spot patterns in the great crested newt (Triturus cristatus) versus the genotypic fingerprint of individuals based on highly polymorphic microsatellite loci using GENECAP. Between 2008 and 2010, we captured, sampled and photographed adult newts and calculated for 1648 samples/photographs recapture rates for both approaches. Recapture rates differed slightly with 8.34% for GENECAP and 9.83% for AMPHIDENT. With an estimated rate of 2% false rejections (FRR) and 0.00% false acceptances (FAR), AMPHIDENT proved to be a highly reliable algorithm for CMR studies of large datasets. We conclude that the application of automatic recognition software of individual photographs can be a rather powerful and reliable tool in noninvasive CMR studies for a large number of individuals. Because the cross-correlation of standardized shape patterns is generally applicable to any pattern that provides enough information, this algorithm is capable of becoming a single application with broad use in CMR studies for many species.

Stronger transferability but lower variability in transcriptomic- than in anonymous microsatellites: evidence from Hylid frogs

A simple way to quickly optimize microsatellites in nonmodel organisms is to reuse loci available in closely related taxa; however, this approach can be limited by the stochastic and low cross-amplification success experienced in some groups (e.g. amphibians). An efficient alternative is to develop loci from transcriptome sequences. Transcriptomic microsatellites have been found to vary in their levels of cross-species amplification and variability, but this has to date never been tested in amphibians. Here, we compare the patterns of cross-amplification and levels of polymorphism of 18 published anonymous microsatellites isolated from genomic DNA vs. 17 loci derived from a transcriptome, across nine species of tree frogs (Hyla arborea and Hyla cinerea group). We established a clear negative relationship between divergence time and amplification success, which was much steeper for anonymous than transcriptomic markers, with half-lives (time at which 50% of the markers still amplify) of 1.1 and 37 My, respectively. Transcriptomic markers are significantly less polymorphic than anonymous loci, but remain variable across diverged taxa. We conclude that the exploitation of amphibian transcriptomes for developing microsatellites seems an optimal approach for multispecies surveys (e.g. analyses of hybrid zones, comparative linkage mapping), whereas anonymous microsatellites may be more informative for fine-scale analyses of intraspecific variation. Moreover, our results confirm the pattern that microsatellite cross-amplification is greatly variable among amphibians and should be assessed independently within target lineages. Finally, we provide a bank of microsatellites for Palaearctic tree frogs (so far only available for H. arborea), which will be useful for conservation and evolutionary studies in this radiation.