Low gene flow but high genetic diversity in the threatened Mallorcan midwife toad Alytes muletensis

Genomic phylogeography illuminates deep cyto-nuclear discordances in midwife toads (Alytes)

The advent of genomic methods allows us to revisit the evolutionary history of organismal groups for which robust phylogenies are still lacking, particularly in species complexes that frequently hybridize. In this study, we conduct RAD-sequencing (RAD-seq) analyses of midwife toads (genus Alytes), an iconic group of western Mediterranean amphibians famous for their parental care behavior, but equally infamous for the difficulties to reconstruct their evolutionary history. Through admixture and phylogenetic analyses of thousands of loci, we provide the most comprehensive phylogeographic framework for the A. obstetricans complex to date, as well as the first fully resolved phylogeny for the entire genus. As part of this effort, we carefully explore the influence of different sampling schemes and data filtering thresholds on tree reconstruction, showing that several, slightly different, yet robust topologies may be retrieved with small datasets obtained by stringent SNP calling parameters, especially when admixed individuals are included. In contrast, analyses of incomplete but larger datasets converged on the same phylogeny, irrespective of the reconstruction method used or the proportion of missing data. The Alytes tree features three Miocene-diverged clades corresponding to the proposed subgenera Ammoryctis (A. cisternasii), Baleaphryne (A. maurus, A. dickhilleni and A. muletensis), and Alytes (A. obstetricans complex). The latter consists of six evolutionary lineages, grouped into three clades of Pliocene origin, and currently delimited as two species: (1) A. almogavarii almogavarii and A. a. inigoi; (2) A. obstetricans obstetricans and A. o. pertinax; (3) A. o. boscai and an undescribed taxon (A. o. cf. boscai). These results contradict the mitochondrial tree, due to past mitochondrial captures in A. a. almogavarii (central Pyrenees) and A. o. boscai (central Iberia) by A. obstetricans ancestors during the Pleistocene. Patterns of admixture between subspecies appear far more extensive than previously assumed from microsatellites, causing nomenclatural uncertainties, and even underlying the reticulate evolution of one taxon (A. o. pertinax). All Ammoryctis and Baleaphryne species form shallow clades, so their taxonomy should remain stable. Amid the prevalence of cyto-nuclear discordance among terrestrial vertebrates and the usual lack of resolution of conventional nuclear markers, our study advocates for phylogeography based on next-generation sequencing, but also encourages properly exploring parameter space and sampling schemes when building and analyzing genomic datasets.

Population structure and genetic diversity of the threatened pygmy newt Triturus pygmaeus in a network of natural and artificial ponds

Pond physical characteristics (connectivity, hydroperiod) have shown to be highly relevant in explaining species presence, reproductive success, and survival in breeding-pond amphibians. However, few studies have addressed the influence that these factors may have on the genetic variability of pond populations. We examined genetic variation at 11 microsatellite loci in Iberian endemic, the pygmy newt (Triturus pygmaeus), from 58 breeding ponds in the Doñana National Park (Andalusia), including both temporary ponds and artificially deepened ponds that remain wet during the whole year. Temporary ponds are located in the North part of the region where the surrounding habitat-wet meadows-facilitates the connectivity among populations, whereas the deepest ponds (‘zacallones’) are located in the southern edge embedded in a matrix of unsuitable habitat (thickets and dry underbrush). We investigated genetic diversity and structure within and among ponds. Our results show that both regions (Doñana-North and Doñana-South) are well-differentiated and form two main clusters. We found higher genetic diversity within ponds from the North region, which also exhibited a higher degree of genetic admixture in comparison with populations from the southern edge. Although we found an isolation-by-distance pattern within each cluster, it arose due to the effect of a few isolated ponds located on the edge of each zone, suggesting the existence of substantial gene flow between ponds in the core area. According to our findings, landscape’s permeability to movement (pond connectivity) may constitute a more important factor than hydroperiod length in determining the genetic diversity and viability of pygmy newt populations in this area. Although Doñana populations show a good state, more peripheral and isolated populations present a more worrisome condition as a result of fragmentation and thus, require conservation efforts. Our study provides key insights that could help guide management practices of this threatened and poorly-studied salamander.

Population genetic structure of the Natterjack toad (Epidalea calamita) in Ireland: implications for conservation management

Molecular methods can play a crucial role in species management and conservation. Despite the usefulness of genetic approaches, they are often not explicitly included as part of species recovery plans and conservation practises. The Natterjack toad (Epidalea calamita) is regionally Red-Listed as Endangered in Ireland. The species is declining and is now present at just seven sites within a highly restricted range. This study used 13 highly polymorphic microsatellite markers to analyse the population genetic diversity and structure. Genetic diversity was high with expected heterozygosity between 0.55 and 0.61 and allelic richness between 4.77 and 5.92. Effective population sizes were small (Ne < 100 individuals), but not abnormal for pond breeding amphibians. However, there was no evidence of historical or contemporary genetic bottlenecks or high levels of inbreeding. We identified a positive relationship between Ne and breeding pond surface area, suggesting that environmental factors are a key determinant of population size. Significant genetic structuring was detected throughout the species’ range, and we identified four genetic entities that should be considered in the species’ conservation strategies. Management should focus on preventing further population declines and future loss of genetic diversity overall and within genetic entities while maintaining adequate local effective population size through site-specific protection, human-mediated translocations and head-start programs. The apparent high levels of genetic variation give hope for the conservation of Ireland’s rarest amphibian if appropriately protected and managed.

Genetic Diversity and Population Differentiation of Pinus koraiensis in China

Pinus koraiensis is a well-known precious tree species in East Asia with high economic, ornamental and ecological value. More than fifty percent of the P. koraiensis forests in the world are distributed in northeast China, a region with abundant germplasm resources. However, these natural P. koraiensis sources are in danger of genetic erosion caused by continuous climate changes, natural disturbances such as wildfire and frequent human activity. Little work has been conducted on the population genetic structure and genetic differentiation of P. koraiensis in China because of the lack of genetic information. In this study, 480 P. koraiensis individuals from 16 natural populations were sampled and genotyped. Fifteen polymorphic expressed sequence tag-simple sequence repeat (EST-SSR) markers were used to evaluate genetic diversity, population structure and differentiation in P. koraiensis. Analysis of molecular variance (AMOVA) of the EST-SSR marker data showed that 33% of the total genetic variation was among populations and 67% was within populations. A high level of genetic diversity was found across the P. koraiensis populations, and the highest levels of genetic diversity were found in HH, ZH, LS and TL populations. Moreover, pairwise Fst values revealed significant genetic differentiation among populations (mean Fst = 0.177). According to the results of the STRUCTURE and Neighbor-joining (NJ) tree analyses and principal component analysis (PCA), the studied geographical populations cluster into two genetic clusters: cluster 1 from Xiaoxinganling Mountains and cluster 2 from Changbaishan Mountains. These results are consistent with the geographical distributions of the populations. The results provide new genetic information for future genome-wide association studies (GWAS), marker-assisted selection (MAS) and genomic selection (GS) in natural P. koraiensis breeding programs and can aid the development of conservation and management strategies for this valuable conifer species.

Urbanization reduces gene flow but not genetic diversity of stream salamander populations in the New York City metropolitan area

Natural landscape heterogeneity and barriers resulting from urbanization can reduce genetic connectivity between populations. The evolutionary, demographic, and ecological effects of reduced connectivity may lead to population isolation and ultimately extinction. Alteration to the terrestrial and aquatic environment caused by urban influence can affect gene flow, specifically for stream salamanders who depend on both landscapes for survival and reproduction. To examine how urbanization affects a relatively common stream salamander species, we compared genetic connectivity of Eurycea bislineata (northern two-lined salamander) populations within and between streams in an urban, suburban, and rural habitat around the New York City (NYC) metropolitan area. We report reduced genetic connectivity between streams within the urban landscape found to correspond with potential barriers to gene flow, that is, areas with more dense urbanization (roadways, industrial buildings, and residential housing). The suburban populations also exhibited areas of reduced connectivity correlated with areas of greater human land use and greater connectivity within a preserve protected from development. Connectivity was relatively high among neighboring rural streams, but a major roadway corresponded with genetic breaks even though the habitat contained more connected green space overall. Despite greater human disturbance across the landscape, urban and suburban salamander populations maintained comparable levels of genetic diversity to their rural counterparts. Yet small effective population size in the urban habitats yielded a high probability of loss of heterozygosity due to genetic drift in the future. In conclusion, urbanization impacted connectivity among stream salamander populations where its continual influence may eventually hinder population persistence for this native species in urban habitats.

Revisiting the genetic variability of Brazilian peridomestic populations of the Chagas disease vector Triatoma sordida () (Hemiptera, Triatominae)

Triatoma sordida is an endemic species to South America, currently considered the species most frequently found in the peridomestic environment in Brazil. This triatomine has a wide ecological tolerability that allows it to inhabit several ecotopes and use different food sources. Although the species is considered predominantly peridomestic, it is also capable of colonizing households and forming numerous intra-household colonies, leading to its inclusion among the priorities for triatomine control campaigns in Brazil. All Brazilian populations of T. sordida are considered as T. sordida sensu stricto by chromosomal analyses (which highlights their epidemiological importance), although molecular studies that characterize the genetic diversity of these populations are scarce, being restricted only to the region of Minas Gerais. Thus, several populations of this vector distributed in the states of Bahia, Goiás, Mato Grosso do Sul, and Minas Gerais were analyzed using mitochondrial markers (cyt b and nd1). Low nucleotide diversity, high haplotypic diversity, low genetic distance, and high F_ST value were observed, as well as the formation of a monophyletic clade of the Brazilian populations of T. sordida, which confirms that this species has low genetic variability, with all specimens in Brazil grouped in T. sordida sensu stricto. In addition to the genetic and evolutionary importance for the knowledge of the biology of these vectors, these results are important from an epidemiological point of view, thus being able to direct vector control programs.

Ready-to-use workflows for the implementation of genetic tools in conservation management

We present a conservation genetics tool kit, which offers two ready-to-use workflows for the routine application of genetic methods in conservation management. The workflows were optimized for work load and costs and are accompanied by an easy-to-read and richly illustrated manual with guidelines regarding sampling design, sampling of genetic material, necessary permits, laboratory methods, statistical analyses and documentation of results in a practice-oriented way. The manual also provides a detailed interpretation help for the implementation of the results in conservation management. One workflow deals with the identification of pond-breeding amphibians based on metabarcoding and environmental DNA (eDNA) from water samples. This workflow also discriminates the morphologically similar water frogs (Pelophylax sp.) and other closely related species (e.g. Triturus cristatus and T. carnifex). The second workflow studies connectivity among populations using microsatellite markers. Its statistical analyses encompass the detection of genetic groups and historical, recent and current dispersal and gene flow. Using the two workflows does not involve academic research institutes; they can be applied by environmental consultancies, laboratories from the private sector, governmental agencies or non-governmental organisations. These and additional conservation genetic workflows will hopefully foster the routine use of genetic methods in conservation management.

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.

Conservation genetics of two threatened frogs from the Mambilla highlands, Nigeria

Amphibians are the vertebrate group with the highest number of species threatened with extinction, and habitat loss and fragmentation are considered to be among the leading causes of their declines and extinctions. Little is known of the population biology of amphibian species inhabiting montane forests in Central and West Africa, where anthropogenic activities such as farming and cattle raising are major threats to native biodiversity. We used Amplified Fragment Length Polymorphisms (AFLPs) to assess the population genetic structure of two poorly known species, Cardioglossa schioetzi and Leptodactylodon bicolor (both in the Arthroleptidae), in and around Ngel Nyaki Forest Reserve on the Mambilla Plateau in eastern Nigeria. The landscape comprises continuous forest on steep slopes and small riparian forest fragments in a grassland matrix. While increased fragmentation is well documented for these and other forests in the mountains of Cameroon and Nigeria over the past century, there are no previous assessments of the impact of forest fragmentation on montane amphibian populations in this region. Our estimates of genetic diversity are similar across populations within each species with levels of heterozygosity values consistent with local population declines. Except for a pair of populations (C. schioetzi) we did not observe genetic differentiation between forest and riparian forest fragment populations, nor across sites within continuous forest (L. bicolor). Our results demonstrate recent gene flow between forest fragments and the adjacent protected forests and suggest that small forest corridors connecting these may lessen the genetic consequences of at least 30 years of intense and severe fragmentation in Ngel Nyaki.

Consequences of population topology for studying gene flow using link-based landscape genetic methods

Many landscape genetic studies aim to determine the effect of landscape on gene flow between populations. These studies frequently employ link‐based methods that relate pairwise measures of historical gene flow to measures of the landscape and the geographical distance between populations. However, apart from landscape and distance, there is a third important factor that can influence historical gene flow, that is, population topology (i.e., the arrangement of populations throughout a landscape). As the population topology is determined in part by the landscape configuration, I argue that it should play a more prominent role in landscape genetics. Making use of existing literature and theoretical examples, I discuss how population topology can influence results in landscape genetic studies and how it can be taken into account to improve the accuracy of these results. In support of my arguments, I have performed a literature review of landscape genetic studies published during the first half of 2015 as well as several computer simulations of gene flow between populations. First, I argue why one should carefully consider which population pairs should be included in link‐based analyses. Second, I discuss several ways in which the population topology can be incorporated in response and explanatory variables. Third, I outline why it is important to sample populations in such a way that a good representation of the population topology is obtained. Fourth, I discuss how statistical testing for link‐based approaches could be influenced by the population topology. I conclude the article with six recommendations geared toward better incorporating population topology in link‐based landscape genetic studies.

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.

High genetic diversity of common toad (Bufo bufo) populations under strong natural fragmentation on a Northern archipelago

The last decades have shown a surge in studies focusing on the interplay between fragmented habitats, genetic variation, and conservation. In the present study, we consider the case of a temperate pond-breeding anuran (the common toad Bufo bufo) inhabiting a naturally strongly fragmented habitat at the Northern fringe of the species' range: islands offshore the Norwegian coast. A total of 475 individuals from 19 populations (three mainland populations and 16 populations on seven adjacent islands) were genetically characterized using nine microsatellite markers. As expected for a highly fragmented habitat, genetic distances between populations were high (pairwise F st values ranging between 0.06 and 0.33), with however little differences between populations separated by ocean and populations separated by terrestrial habitat (mainland and on islands). Despite a distinct cline in genetic variation from mainland populations to peripheral islands, the study populations were characterized by overall high genetic variation, in line with effective population sizes derived from single-sample estimators which were on average about 20 individuals. Taken together, our results reinforce the notion that spatial and temporal scales of fragmentation need to be considered when studying the interplay between landscape fragmentation and genetic erosion.

High genetic diversity in the endangered and narrowly distributed amphibian species Leptobrachium leishanense

Threatened species typically have a small or declining population size, which make them highly susceptible to loss of genetic diversity through genetic drift and inbreeding. Genetic diversity determines the evolutionary potential of a species; therefore, maintaining the genetic diversity of threatened species is essential for their conservation. In this study, we assessed the genetic diversity of the adaptive major histocompatibility complex (MHC) genes in an endangered and narrowly distributed amphibian species, Leptobrachium leishanense in Southwest China. We compared the genetic variation of MHC class I genes with that observed in neutral markers (5 microsatellite loci and cytochrome b gene) to elucidate the relative roles of genetic drift and natural selection in shaping the current MHC polymorphism in this species. We found a high level of genetic diversity in this population at both MHC and neutral markers compared with other threatened amphibian species. Historical positive selection was evident in the MHC class I genes. The higher allelic richness in MHC markers compared with that of microsatellite loci suggests that selection rather than genetic drift plays a prominent role in shaping the MHC variation pattern, as drift can affect all the genome in a similar way but selection directly targets MHC genes. Although demographic analysis revealed no recent bottleneck events in L. leishanense, additional population decline will accelerate the dangerous status for this species. We suggest that the conservation management of L. leishanense should concentrate on maximizing the retention of genetic diversity through preventing their continuous population decline. Protecting their living habitats and forbidding illegal hunting are the most important measures for conservation of L. leishanense.

Genetic diversity and population structure of the critically endangered Yangtze finless porpoise (Neophocaena asiaeorientalis asiaeorientalis) as revealed by mitochondrial and microsatellite DNA

Ecological surveys have indicated that the population of the critically endangered Yangtze finless porpoise (YFP, Neophocaena asiaeorientalis asiaeorientalis) is becoming increasingly small and fragmented, and will be at high risk of extinction in the near future. Genetic conservation of this population will be an important component of the long-term conservation effort. We used a 597 base pair mitochondrial DNA (mtDNA) control region and 11 microsatellite loci to analyze the genetic diversity and population structure of the YFP. The analysis of both mtDNA and microsatellite loci suggested that the genetic diversity of the YFP will possibly decrease in the future if the population keeps declining at a rapid rate, even though these two types of markers revealed different levels of genetic diversity. In addition, mtDNA revealed strong genetic differentiation between one local population, Xingchang–Shishou (XCSS), and the other five downstream local populations; furthermore, microsatellite DNA unveiled fine but significant genetic differentiation between three of the local populations (not only XCSS but also Poyang Lake (PY) and Tongling (TL)) and the other local populations. With an increasing number of distribution gaps appearing in the Yangtze main steam, the genetic differentiation of local populations will likely intensify in the future. The YFP is becoming a genetically fragmented population. Therefore, we recommend attention should be paid to the genetic conservation of the YFP.

Population structure, historical biogeography and demographic history of the alpine toad Scutiger ningshanensis in the Tsinling Mountains of Central China

Population genetic structure, historical biogeography and historical demography of the alpine toad Scutiger ningshanensis were studied using the combined data mtDNA cytochrome b (cyt b) and the mtDNA cytochrome c oxidase subunit I (COI) as the molecular markers. This species has high genetic variation. There was a significant genetic differentiation among most populations. Three lineages were detected. The phylogenetic relationship analyses and the SAMOVA (spatial analysis of molecular variance) results showed significant phylogeographic structure. 82.15% genetic variation occurred among populations whereas differentiation within populations only contributed 17.85% to the total. Mantel test results showed a significant correlation between the pairwise calculated genetic distance and pairwise calculated geographical distance of the populations (regression coefficient  = 0.001286, correlation coefficient  = 0.77051, p (rrand≥robs)  = 0.0185<0.05), indicating the existence of isolation-by-distance pattern of genetic divergence for cyt b + COI sequence, which suggests that the distribution of genetic variation is due to geographical separation rather than natural selection. The population expansion or contraction and genetic differentiation between populations or lineages could be explained by topography and the repetitive uplifts of the Tsinling Mountains and the climatic cycles during the late Pliocene and Pleistocene. S. ningshanensis experienced a rapid population expansion about 40,000 years before present. The current decline in population size was probably caused by anthropogenic disturbance. Current populations of S. ningshanensis are from different refugia though the location of these refugia could not be determined in our study. Topography, climatic changes and repetitive population expansion/contraction together led to the high level of genetic variation in S. ningshanensis. A total of three management units (MUs) was determined, which must be considered when conservation policy is made in the future.

Molecular phylogeography and population genetic structure of an endangered species Pachyhynobius shangchengensis (hynobiid Salamander) in a fragmented habitat of southeastern China

The salamander Pachyhynobius shangchengensis (Hynobiidae) is a vulnerable species restricted to a patchy distribution associated with small mountain streams surrounded by forested slopes in the Mount Dabieshan region in southeastern China. However, molecular phylogeography and population genetic structure of P. shangchengensis remain poorly investigated. In this study, we explored the genetic structure and phylogeography of P. shangchengensis based on partial sequences of the mitochondrial DNA (mtDNA) cytochrome b and cytochrome c oxidase subunit I genes. Fifty-one haplotypes and four clades were found among 93 samples. Phylogenetic analyses revealed four deeply divergent and reciprocally monophyletic mtDNA lineages that approximately correspond to four geographic regions separated by complicated topography and long distances. The distinct geographic distributions of all lineages and the estimated divergence time suggest spatial and temporal separation coinciding with climatic changes during the Pleistocene. Analysis of molecular variance indicated that most of the observed genetic variation occurred among the four groups, implying long-term interruption of gene flow, and the possible separation of P. shangchengensis into four management units for conservation.

Phylogeographic analyses strongly suggest cryptic speciation in the giant spiny frog (Dicroglossidae: Paa spinosa) and interspecies hybridization in Paa

Species identification is one of the most basic yet crucial issues in biology with potentially far-reaching implications for fields such as conservation, population ecology, and epidemiology. The widely distributed but threatened frog Paa spinosa has been speculated to represent a complex of multiple species. In this study, 254 individuals representing species of the genus Paa were investigated along the entire range of P. spinosa: 196 specimens of P. spinosa, 8 specimens of P. jiulongensis, 5 specimens of P. boulengeri, 20 specimens of P. exilispinosa, and 25 specimens of P. shini. Approximately 1333 bp of mtDNA sequence data (genes 12S rRNA and 16S rRNA) were used. Phylogenetic analyses were conducted using maximum parsimony, maximum likelihood and Bayesian inference. BEAST was used to estimate divergence dates of major clades. Results suggest that P. spinosa can be divided into three distinct major lineages. Each major lineage totally corresponds to geographical regions, revealing the presence of three candidate cryptic species. Isolation and differentiation among lineages are further supported by the great genetic distances between the lineages. The bifurcating phylogenetic pattern also suggests an east-west dispersal trend during historic cryptic speciation. Dating analysis estimates that P. spinosa from Western China split from the remaining P. spinosa populations in the Miocene and that P. spinosa from Eastern China diverged from Central China in the Pliocene. We also found that P. exilispinosa from Mainland China and Hong Kong might have a complex of multiple species. After identifying cryptic lineages, we then determine the discrepancy between the mtDNA and the morphotypes in several individuals. This discrepancy may have been caused by introgressive hybridization between P. spinosa and P. shini.

Study of closely related species within the Physalaemus cuvieri group (Anura): contribution of microsatellite markers

Various species of the Physalaemus cuvieri group of frogs are difficult to distinguish morphologically, making molecular analysis an attractive alternative for indentifying members of this group, which is considered to be at risk because of loss of habitat. The genetic structure of natural populations of P. ephippifer and P. albonotatus species was investigated and analyzed, together with that of five previously studied populations of P. cuvieri. Nine microsatellite loci were used in the analyses. The overall G(ST) value (0.46) revealed high genetic variation among the populations, as expected for different species. Bayesian analysis implemented by the STRUCTURE software clustered the seven populations into seven groups (K = 7). All the P. albonotatus and P. ephippifer specimens were grouped into a single cluster, both species showing clear differentiation from P. cuvieri. The different grouping based on these microsatellites of some P. cuvieri individuals from Porto Nacional and from Passo Fundo suggests that they could be a new species, indicating a necessity for taxonomic reevaluation. Despite the intrinsic difficulties in analyzing closely related species, the nine microsatellite loci were found to be adequate for distinguishing these three species of the P. cuvieri group and their populations.

Mitigating amphibian disease: strategies to maintain wild populations and control chytridiomycosis

Background

Rescuing amphibian diversity is an achievable conservation challenge. Disease mitigation is one essential component of population management. Here we assess existing disease mitigation strategies, some in early experimental stages, which focus on the globally emerging chytrid fungus Batrachochytrium dendrobatidis. We discuss the precedent for each strategy in systems ranging from agriculture to human medicine, and the outlook for each strategy in terms of research needs and long-term potential.

Results

We find that the effects of exposure to Batrachochytrium dendrobatidis occur on a spectrum from transient commensal to lethal pathogen. Management priorities are divided between (1) halting pathogen spread and developing survival assurance colonies, and (2) prophylactic or remedial disease treatment. Epidemiological models of chytridiomycosis suggest that mitigation strategies can control disease without eliminating the pathogen. Ecological ethics guide wildlife disease research, but several ethical questions remain for managing disease in the field.

Conclusions

Because sustainable conservation of amphibians in nature is dependent on long-term population persistence and co-evolution with potentially lethal pathogens, we suggest that disease mitigation not focus exclusively on the elimination or containment of the pathogen, or on the captive breeding of amphibian hosts. Rather, successful disease mitigation must be context specific with epidemiologically informed strategies to manage already infected populations by decreasing pathogenicity and host susceptibility. We propose population level treatments based on three steps: first, identify mechanisms of disease suppression; second, parameterize epizootiological models of disease and population dynamics for testing under semi-natural conditions; and third, begin a process of adaptive management in field trials with natural populations.

Ecological and genetic divergence between two lineages of middle American túngara frogs Physalaemus (= Engystomops) pustulosus

BACKGROUND

Uncovering how populations of a species differ genetically and ecologically is important for understanding evolutionary processes. Here we combine population genetic methods (microsatellites) with phylogenetic information (mtDNA) to define genetic population clusters of the wide-spread Neotropical túngara frog (Physalaemus pustulosus). We measure gene flow and migration within and between population clusters and compare genetic diversity between population clusters. By applying ecological niche modeling we determine whether the two most divergent genetic groups of the túngara frog (1) inhabit different habitats, and (2) are separated geographically by unsuitable habitat across a gap in the distribution.

RESULTS

Most population structure is captured by dividing all sample localities into two allopatric genetic lineages. The Northern genetic lineage (NW Costa Rica) is genetically homogenous while the Southern lineage (SW Costa Rica and Panama) is sub-divided into three population clusters by both microsatellite and mtDNA analyses. Gene flow is higher within the Northern lineage than within the Southern lineage, perhaps due to increased landscape heterogeneity in the South. Niche modeling reveals differences in suitable habitat between the Northern and Southern lineages: the Northern lineage inhabits dry/pine-oak forests, while the Southern lineage is confined to tropical moist forests. Both lineages seem to have had little movement across the distribution gap, which persisted during the last glacial maximum. The lack of movement was more pronounced for the Southern lineage than for the Northern lineage.

CONCLUSIONS

This study confirms the finding of previous studies that túngara frogs diverged into two allopatric genetic lineages north and south of the gap in the distribution in central Costa Rica several million years ago. The allopatric distribution is attributed to unsuitable habitat and probably other unknown ecological factors present across the distribution gap. Niche conservatism possibly contributes to preventing movements across the gap and gene flow between both groups. Genetic and ecological data indicate that there is the potential for ecological divergence in allopatry between lineages. In this context we discuss whether the Northern and Southern lineages should be recognized as separate species, and we conclude that further studies of pre- and post-zygotic isolation are needed for a final assessment. Identified population clusters should motivate future behavioral and ecological research regarding within-species biodiversity and speciation mechanisms.

Genetic diversity and molecular differentiation of Chinese toad based on microsatellite markers

Genetic diversity and population structure of 9 populations of Bufo gargarizans with total 111 samples in China were assessed using seven microsatellite loci. The analysed microsatellite markers produced 161 alleles, varied from 9 to 38 alleles each locus. The number of alleles per population per locus ranged from 4.43 to 10.29. Polymorphic information content showed that all seven loci were highly informative (mean = 0.810 + or - 0.071). The average observed heterozygosity was less than the expected (0.353 + or - 0.051 and 0.828 + or - 0.067, respectively). All tested populations gave significant departures from Hardy-Weinberg equilibrium. Genetic differentiation among the populations was considerably high with the overall and pairwise F (ST) values (mean = 0.160 + or - 0.039), and showed fairly high level of inbreeding (indicated by a mean F (IS) value of 0.504 + or - 0.051) and global heterozygote deficit. In comparison to other amphibian studies; however, our results suggested that the level of genetic structuring in B. gargarizans was relatively low in the geographical scale of the study area. Interestingly, the speculated population bottleneck was found to be absent and the analyses provide only weak evidence for a recent contraction in size even though there was severe inbreeding (indicated by the F (IS) value) in the Chinese toad populations.

Assessment of genetic diversity in the critically endangered Australian corroboree frogs, Pseudophryne corroboree and Pseudophryne pengilleyi, identifies four evolutionarily significant units for conservation

The iconic and brightly coloured Australian northern corroboree frog, Pseudophryne pengilleyi, and the southern corroboree frog, Pseudophryne corroboree are critically endangered and may be extinct in the wild within 3 years. We have assembled samples that cover the current range of both species and applied hypervariable microsatellite markers and mitochondrial DNA sequences to assess the levels and patterns of genetic variation. The four loci used in the study were highly variable, the total number of alleles observed ranged from 13 to 30 and the average number of alleles per locus was 19. Expected heterozygosity of the four microsatellite loci across all populations was high and varied between 0.830 and 0.935. Bayesian clustering analyses in STRUCTURE strongly supported four genetically distinct populations, which correspond exactly to the four main allopatric geographical regions in which the frogs are currently found. Individual analyses performed on the separate regions showed that breeding sites within these four regions could not be separated into distinct populations. Twelve mtND2 haplotypes were identified from 66 individuals from throughout the four geographical regions. A statistical parsimony network of mtDNA haplotypes shows two distinct groups, which correspond to the two species of corroboree frog, but with most of the haplotype diversity distributed in P. pengilleyi. These results demonstrate an unexpectedly high level of genetic diversity in both species. Our data have important implications for how the genetic diversity is managed in the future. The four evolutionarily significant units must be protected and maintained in captive breeding programmes for as long as it is possible to do.

Postglacial range expansion from northern refugia by the wood frog, Rana sylvatica

Although the range dynamics of North American amphibians during the last glacial cycle are increasingly better understood, the recolonization history of the most northern regions and the impact of southern refugia on patterns of intraspecific genetic diversity and phenotypic variation in these regions are not well reconstructed. Here we present the phylogeographic history of a widespread and primarily northern frog, Rana sylvatica. We surveyed 551 individuals from 116 localities across the species' range for a 650-bp region of the NADH dehydrogenase subunit 2 and tRNA(TRP) mitochondrial genes. Our phylogenetic analyses revealed two distinct clades corresponding to eastern and western populations, as well as a Maritime subclade within the eastern lineage. Patterns of genetic diversity support multiple refugia. However, high-latitude refugia in the Appalachian highlands and modern-day Wisconsin appear to have had the biggest impact on northern populations. Clustering analyses based on morphology further support a distinction between eastern and western wood frogs and suggest that postglacial migration has played an important role in generating broad-scale patterns of phenotypic variation in this species.

The influence of altitude and topography on genetic structure in the long-toed salamander (Ambystoma macrodactulym)

A primary goal of molecular ecology is to understand the influence of abiotic factors on the spatial distribution of genetic variation. Features including altitudinal clines, topography and landscape characteristics affect the proportion of suitable habitat, influence dispersal patterns, and ultimately structure genetic differentiation among populations. We studied the effects of altitude and topography on genetic variation of long-toed salamanders (Ambystoma macrodactylum), a geographically widespread amphibian species throughout northwestern North America. We focused on 10 low altitude sites (< 1200 m) and 11 high-altitude sites in northwestern Montana and determined multilocus genotypes for 549 individuals using seven microsatellite loci. We tested four hypotheses: (1) gene flow is limited between high- and low-altitude sites; and, (2) gene flow is limited among high-altitude sites due to harsh habitat and extreme topographical relief between sites; (3) low-altitude sites exhibit higher among-site gene flow due to frequent flooding events and low altitudinal relief; and (4) there is a negative correlation between altitude and genetic variation. Overall F(ST) values were moderate (0.08611; P < 0.001). Pairwise F(ST) estimates between high and low populations and a population graphing method supported the hypothesis that low-altitude and high-altitude sites, taken together, are genetically differentiated from each other. Also as predicted, gene flow is more prominent among low-altitude sites than high-altitude sites; low-altitude sites had a significantly lower F(ST) (0.03995; P < 0.001) than high altitude sites (F(ST) = 0.10271; P < 0.001). Use of Bayesian analysis of population structure (BAPS) resulted in delineation of 10 genetic groups, two among low-altitude populations and eight among high-altitude populations. In addition, within high altitude populations, basin-level genetic structuring was apparent. A nonequilibrium algorithm for detecting current migration rates supported these population distinctions. Finally, we also found a significant negative correlation between genetic diversity and altitude. These results are consistent with the hypothesis that topography and altitudinal gradients shape the spatial distribution of genetic variation in a species with a broad geographical range and diverse life history. Our study sheds light on which key factors limit dispersal and ultimately species' distributions.

Genetic diversity, but not hatching success, is jointly affected by postglacial colonization and isolation in the threatened frog, Rana latastei

Both postglacial colonization and habitat fragmentation can reduce the genetic diversity of populations, which in turn can affect fitness. However, since these processes occur at different spatial and temporal scales, the consequences of either process may differ. To disentangle the relative role of isolation and postglacial colonization in determining genetic diversity and fitness, we studied microsatellite diversity of 295 individuals from 10 populations and measured the hatch rate of 218 clutches from eight populations of a threatened frog, R. latastei. The populations that were affected by fragmentation to a greater extent suffered higher embryo mortality and reduced hatch rate, while no effects of distance from glacial refugium on hatch rate were detected. Altogether, distance from glacial refugium and isolation explained > 90% of variation in genetic diversity. We found that the genetic diversity was lowest in populations both isolated and far from the glacial refugium, and that distance from refugium seems to have the primary role in determining genetic diversity. The relationship between genetic diversity and hatch rate was not significant. However, the proportion of genetic diversity lost through recent isolation had a significant, negative effect on fitness. It is possible that selection at least partially purged the negative effects of the ancestral loss of genetic diversity.