Buccal swabs allow efficient and reliable microsatellite genotyping in amphibians

Optimising recovery of DNA from minimally invasive sampling methods: Efficacy of buccal swabs, preservation strategy and DNA extraction approaches for amphibian studies

Studies in evolution, ecology and conservation are increasingly based on genetic and genomic data. With increased focus on molecular approaches, ethical concerns about destructive or more invasive techniques need to be considered, with a push for minimally invasive sampling to be optimised. Buccal swabs have been increasingly used to collect DNA in a number of taxa, including amphibians. However, DNA yield and purity from swabs are often low, limiting its use. In this study, we compare different types of swabs, preservation method and storage, and DNA extraction techniques in three case studies to assess the optimal approach for recovering DNA in anurans. Out of the five different types of swabs that we tested, Isohelix MS-02 and Rapidry swabs generated higher DNA yields than other swabs. When comparing storage buffers, ethanol is a better preservative than a non-alcoholic alternative. Dried samples resulted in similar or better final DNA yields compared to ethanol-fixed samples if kept cool. DNA extraction via a Qiagen™ DNeasy Blood and Tissue Kit and McHale's salting-out extraction method resulted in similar DNA yields but the Qiagen™ kit extracts contained less contamination. We also found that samples have better DNA recovery if they are frozen as soon as possible after collection. We provide recommendations for sample collection and extraction under different conditions, including budgetary considerations, size of individual animal sampled, access to cold storage facilities and DNA extraction methodology. Maximising efficacy of all of these factors for better DNA recovery will allow buccal swabs to be used for genetic and genomic studies in a range of vertebrates.

Comparing skin swabs, buccal swabs, and toe clips for amphibian genetic sampling, a case study with a small anuran (Acris blanchardi)

Multiple methods for collecting genetic samples from amphibians exist, each with their own implications for study design, animal welfare, and costs. Toe clipping is one common method, but there is ongoing debate regarding its potential detriment. Less invasive methods should be implemented, if efficacious, as amphibians are a particularly vulnerable vertebrate group. Skin and buccal swabbing are less invasive methods for genetic sampling, but the potential for contamination and a lower yield of DNA may exist. To compare these methods, we gathered skin swabs, buccal swabs, and toe clips from the same individuals of a relatively small anuran species, Blanchard's Cricket Frog (Acris blanchardi). We then compared DNA yield, DNA purity, amplification success rate, and genotypic data quality among sample types. We found toe clips and buccal swabs generated similar DNA yield and purity, with skin swabs yielding significantly less DNA of significantly lower purity than the other sample types. Amplification success rate was significantly higher using toe clips compared to the other sample types, though buccal swab samples amplified more readily than skin swabs. Genotypic data from toe clips and buccal swabs did not differ significantly in quality, but skin swab data quality was significantly lowest among sample types. Thus, skin swabbing could produce erroneous data in some situations, but buccal swabbing is likely an effective substitute to toe clipping, even for small species. Our results can help future researchers select which genetic sampling method might best suit their research needs.

Non-Lethal Detection of Frog Virus 3-Like (RUK13) and Common Midwife Toad Virus-Like (PDE18) Ranaviruses in Two UK-Native Amphibian Species

Ranaviruses have been involved in amphibian mass mortality events worldwide. Effective screening to control this pathogen is essential; however, current sampling methods are unsuitable for the detection of subclinical infections. Non-lethal screening is needed to prevent both further spread of ranavirus and losses of at-risk species. To assess non-lethal sampling methods, we conducted two experiments: bath exposing common frogs to RUK13 ranavirus at three concentrations, and exposing common toads to RUK13 or PDE18. Non-lethal sampling included buccal, digit, body and tank swabs, along with toe clips and stool taken across three time-points post-exposure. The presence/load of ranavirus was examined using quantitative PCR in 11 different tissues obtained from the same euthanised animals (incl. liver, gastro-intestinal tract and kidney). Buccal swab screening had the highest virus detection rate in both species (62% frogs; 71% toads) and produced consistently high virus levels compared to other non-lethal assays. The buccal swab was effective across multiple stages of infection and differing infection intensities, though low levels of infection were more difficult to detect. Buccal swab assays competed with, and even outperformed, lethal sampling in frogs and toads, respectively. Successful virus detection in the absence of clinical signs was observed (33% frogs; 50% toads); we found no difference in detectability for RUK13 and PDE18. Our results suggest that buccal swabbing could replace lethal sampling for screening and be introduced as standard practice for ranavirus surveillance.

Morphometry of two cryptic tree frog species at their hybrid zone reveals neither intermediate nor transgressive morphotypes

Under incomplete reproductive isolation, secondary contact of diverged allopatric lineages may lead to the formation of hybrid zones that allow to study recombinants over several generations as excellent systems of genomic interactions resulting from the evolutionary forces acting on certain genes and phenotypes. Hybrid phenotypes may either exhibit intermediacy or, alternatively, transgressive traits, which exceed the extremes of their parents due to epistasis and segregation of complementary alleles. While transgressive morphotypes have been examined in fish, reptiles, birds, and mammals, studies in amphibians are rare. Here, we associate microsatellite-based genotypes with morphometrics-based morphotypes of two tree frog species of the Hyla arborea group, sampled across a hybrid zone in Poland, to understand whether the genetically differentiated parental species also differ in morphology between each other and their hybrids and whether secondary contact leads to the evolution of intermediate or transgressive morphotypes. Using univariate approaches, explorative multivariate methods (principal component analyses) as well as techniques with prior grouping (discriminant function analyses), we find that morphotypes of both parental species and hybrids differ from each other. Importantly, hybrid morphotypes are neither intermediate nor transgressive but found to be more similar to H. orientalis than to H. arborea.

Differential gene expression of the healthy conjunctiva during the day

PURPOSE

To determine if there is diurnal variation in gene expression in normal healthy conjunctival cells.

METHODS

Bulbar conjunctival swab samples were collected from four healthy subjects in the morning and evening of the same day. The two swab samples were taken from one eye of each participant, with a minimum of five hours gap between the two samples. RNA was extracted and analysed using RNA sequencing (RNA-Seq).

RESULTS

A total of 121 genes were differentially expressed between the morning and the evening conjunctival samples, of which 94 genes were upregulated in the morning, and 27 genes were upregulated in the evening. Many of the genes that were upregulated in the morning were involved in defence, cell turnover and regulation of gene expression, while the genes upregulated in the evening were involved in signalling and mucin production.

CONCLUSIONS

This study has identified several genes whose expression changes over the course of the day. Knowledge of diurnal variations of conjunctival gene expression provides an insight into the regulatory status of the healthy eye and provides a baseline for examining changes during ocular surface disease.

Impact of the Miocene orogenesis on Kaloula spp. radiation and implication of local refugia on genetic diversification

The phylogeography of the Kaloula genus in East Asia is still poorly understood. One of the difficulties is the absence of fossils to corroborate molecular dating estimates. Here, we examined the mitochondrial structure of Kaloula spp. in East Asia and focused on the impact of glaciations on the northernmost species: Kaloula borealis. We determined the phylogenetic relationships, molecular dating, and genetic connectivity assessments within the genus from 1211 bp of concatenated mitochondrial 12S and 16S. The relaxed clock analyses reveal the emergence of Kaloula spp. common ancestor in East and Southeast Asia between the Eocene and Oligocene, c. 38.47 Ma (24.69-53.65). The genetic diversification of lineages then increased on the East Asian Mainland during the Lower Miocene, c. 20.10 (8.73-30.65), most likely originating from the vicariance and radiation triggered by the orogeny of the Qinghai-Tibetan Plateau. Later, the dispersal towards the North East Asian Mainland during the Upper Miocene drove the population diversification of K. borealis c. 9.01 Ma (3.66-15.29). Finally, the central mainland population became isolated following orogenesis events and diverged into K. rugifera during the Pliocene, c. 3.06 Ma (0.02-10.90). The combination of population genetic and barrier analyses revealed a significant genetic isolation between populations of Kaloula spp. matching with the massive Qinling-Daba Mountain chain located in south-central China. Finally, we highlight a young divergence within the Eastern Mainland population of K. borealis, possibly attributed to refugia in south eastern China from which populations later expanded.

Phylogeographic analysis and species distribution modelling of the wood frog Batrachyla leptopus (Batrachylidae) reveal interglacial diversification in south western Patagonia

Background

The evolutionary history of southern South American organisms has been strongly influenced by Pleistocene climate oscillations. Amphibians are good models to evaluate hypotheses about the influence of these climate cycles on population structure and diversification of the biota, because they are sensitive to environmental changes and have restricted dispersal capabilities. We test hypotheses regarding putative forest refugia and expansion events associated with past climatic changes in the wood frog Batrachyla leptopus distributed along ∼1,000 km of length including glaciated and non-glaciated areas in southwestern Patagonia.

Methods

Using three mitochondrial regions (D-loop, cyt b, and coI) and two nuclear loci (pomc and crybA1), we conducted multilocus phylogeographic analyses and species distribution modelling to gain insights of the evolutionary history of this species. Intraspecific genealogy was explored with maximum likelihood, Bayesian, and phylogenetic network approaches. Diversification time was assessed using molecular clock models in a Bayesian framework, and demographic scenarios were evaluated using approximate Bayesian computation (ABC) and extended Bayesian skyline plot (EBSP). Species distribution models (SDM) were reconstructed using climatic and geographic data.

Results

Population structure and genealogical analyses support the existence of four lineages distributed north to south, with moderate to high phylogenetic support (Bootstrap > 70%; BPP > 0.92). The diversification time of B. leptopus’ populations began at ∼0.107 mya. The divergence between A and B lineages would have occurred by the late Pleistocene, approximately 0.068 mya, and divergence between C and D lineages was approximately 0.065 mya. The ABC simulations indicate that lineages coalesced at two different time periods, suggesting the presence of at least two glacial refugia and a postglacial colonization route that may have generated two southern lineages (p = 0.93, type I error: <0.094, type II error: 0.134). EBSP, mismatch distribution and neutrality indexes suggest sudden population expansion at ∼0.02 mya for all lineages. SDM infers fragmented distributions of B. leptopus associated with Pleistocene glaciations. Although the present populations of B. leptopus are found in zones affected by the last glacial maximum (∼0.023 mya), our analyses recover an older history of interglacial diversification (0.107–0.019 mya). In addition, we hypothesize two glacial refugia and three interglacial colonization routes, one of which gave rise to two expanding lineages in the south.

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.

Large-Scale Hybridisation as an Extinction Threat to the Suweon Treefrog (Hylidae: Dryophytes suweonensis)

Simple Summary

A large number of amphibian species are now endangered, mostly because of human activities. An example is land modification, which may bring species that were previously isolated in contact, and allows them to hybridise. Here, we assessed the presence of hybrid individuals between the endangered Suweon treefrog (Dryophytes suweonensis) and the widespread Japanese treefrog (Dryophytes japonicus). We found hybrids to be relatively widespread and present at all populations where the Suweon treefrog occurred. This is important, as it results in an additional threat to the Suweon treefrog.

Abstract

Amphibians are in the midst of a sixth mass extinction, and human activities play a major role in pushing species towards extinction. Landscape anthropisation has impacts that indirectly threaten species, in addition to the obvious destruction of natural habitats. For instance, land modification may bring human-commensal species in contact with sister-clades from which they were previously isolated. The species in these new contact zones are then able to hybridise to the point of reaching lineage fusion, through which the gene pool of the two species merges and one of the parental lineages becomes extirpated. Here, we documented the patterns of hybridisation between the spatially restricted D. suweonensis and the widespread D. japonicus. On the basis of the analysis of Cytochrome c oxidase subunit I mitochondrial DNA sequences (404 individuals from 35 sites) and six polymorphic microsatellites (381 individuals from 34 sites), we revealed a generalised, bi-directional, and geographically widespread hybridisation between the two species. Evidence of fertile back-crosses is provided by relatively high numbers of individuals in cyto-nuclear disequilibrium, as well as the presence of hybrid individuals further south than the species distribution limit, determined on the basis of call properties. Hybridisation is an additional threat to the endangered D. suweonensis.

No evidence that Y-chromosome differentiation affects male fitness in a Swiss population of common frogs

The canonical model of sex-chromosome evolution assigns a key role to sexually antagonistic (SA) genes on the arrest of recombination and ensuing degeneration of Y chromosomes. This assumption cannot be tested in organisms with highly differentiated sex chromosomes, such as mammals or birds, owing to the lack of polymorphism. Fixation of SA alleles, furthermore, might be the consequence rather than the cause of recombination arrest. Here we focus on a population of common frogs (Rana temporaria) where XY males with genetically differentiated Y chromosomes (nonrecombinant Y haplotypes) coexist with both XY° males with proto-Y chromosomes (only differentiated from X chromosomes in the immediate vicinity of the candidate sex-determining locus Dmrt1) and XX males with undifferentiated sex chromosomes (genetically identical to XX females). Our study finds no effect of sex-chromosome differentiation on male phenotype, mating success or fathering success. Our conclusions rejoin genomic studies that found no differences in gene expression between XY, XY° and XX males. Sexual dimorphism in common frogs might result more from the differential expression of autosomal genes than from sex-linked SA genes. Among-male variance in sex-chromosome differentiation seems better explained by a polymorphism in the penetrance of alleles at the sex locus, resulting in variable levels of sex reversal (and thus of X-Y recombination in XY females), independent of sex-linked SA genes.

Phylogeny of spiny frogs Nanorana (Anura: Dicroglossidae) supports a Tibetan origin of a Himalayan species group

Recent advances in the understanding of the evolution of the Asian continent challenge the long-held belief of a faunal immigration into the Himalaya. Spiny frogs of the genus Nanorana are a characteristic faunal group of the Himalaya-Tibet orogen (HTO). We examine the phylogeny of these frogs to explore alternative biogeographic scenarios for their origin in the Greater Himalaya, namely, immigration, South Tibetan origin, strict vicariance. We sequenced 150 Nanorana samples from 62 localities for three mitochondrial (1,524 bp) and three nuclear markers (2,043 bp) and complemented the data with sequence data available from GenBank. We reconstructed a gene tree, phylogenetic networks, and ancestral areas. Based on the nuDNA, we also generated a time-calibrated species tree. The results revealed two major clades (Nanorana and Quasipaa), which originated in the Lower Miocene from eastern China and subsequently spread into the HTO (Nanorana). Five well-supported subclades are found within Nanorana: from the East, Central, and Northwest Himalaya, the Tibetan Plateau, and the southeastern Plateau margin. The latter subclade represents the most basal group (subgenus Chaparana), the Plateau group (Nanorana) represents the sister clade to all species of the Greater Himalaya (Paa). We found no evidence for an east-west range expansion of Paa along the Himalaya, nor clear support for a strict vicariance model. Diversification in each of the three Himalayan subclades has probably occurred in distinct areas. Specimens from the NW Himalaya are placed basally relative to the highly diverse Central Himalayan group, while the lineage from the Tibetan Plateau is placed within a more terminal clade. Our data indicate a Tibetan origin of Himalayan Nanorana and support a previous hypothesis, which implies that a significant part of the Himalayan biodiversity results from primary diversification of the species groups in South Tibet before this part of the HTO was uplifted to its recent heights.

Using Noninvasive Genetic Sampling to Survey Rare Butterfly Populations

Advances in nondestructive genetic sampling techniques continue to offer new opportunities for studying organisms, particularly those of conservation concern where more traditional invasive sampling methods are often not available. As part of a proof-of-concept, we investigated the effectiveness of using the chorion from residual butterfly egg debris as a source of viable genetic material for analysis. Laboratory material from a captive breeding population of the federally endangered Miami blue butterfly (Cyclargus thomasi bethunebakeri) was used to test efficacy and refine the methodology. The resulting best practices were subsequently evaluated using field-collected material from extant north Florida populations of the at-risk frosted elfin butterfly (Callophyrs irus). Our results demonstrated that it is possible to extract DNA of sufficiently high quantity and quality for successful gene sequencing. We additionally describe a simple, low-cost, and reliable method of collecting and storing egg debris samples that can be consistently adopted for field or laboratory work as well as deployed with projects that have a larger geographic scope and/or involve citizen scientists. Potential limitations related to field sample collection are discussed as well as needs for future evaluation.

Impacts of the synthetic androgen Trenbolone on gonad differentiation and development - comparisons between three deeply diverged anuran families

Using a recently developed approach for testing endocrine disruptive chemicals (EDCs) in amphibians, comprising synchronized tadpole exposure plus genetic and histological sexing of metamorphs in a flow-through-system, we tested the effects of 17β-Trenbolone (Tb), a widely used growth promoter in cattle farming, in three deeply diverged anuran families: the amphibian model species Xenopus laevis (Pipidae) and the non-models Bufo(tes) viridis (Bufonidae) and Hyla arborea (Hylidae). Trenbolone was applied in three environmentally and/or physiologically relevant concentrations (0.027 µg/L (10-10 M), 0.27 µg/L (10-9 M), 2.7 µg/L (10-8 M)). In none of the species, Tb caused sex reversals or masculinization of gonads but had negative species-specific impacts on gonad morphology and differentiation after the completion of metamorphosis, independently of genetic sex. In H. arborea and B. viridis, mounting Tb-concentration correlated positively with anatomical abnormalities at 27 µg/L (10-9 M) and 2.7 µg/L (10-8 M), occurring in X. laevis only at the highest Tb concentration. Despite anatomical aberrations, histologically all gonadal tissues differentiated seemingly normally when examined at the histological level but at various rates. Tb-concentration caused various species-specific mortalities (low in Xenopus, uncertain in Bufo). Our data suggest that deep phylogenetic divergence modifies EDC-vulnerability, as previously demonstrated for Bisphenol A (BPA) and Ethinylestradiol (EE2).

Speciation in a biodiversity hotspot: Phylogenetic relationships, species delimitation, and divergence times of Patagonian ground frogs from the Eupsophus roseus group (Alsodidae)

The alsodid ground frogs of the Eupsophus genus are divided into two groups, the roseus (2n = 30) and vertebralis (2n = 28), which are distributed throughout the temperate Nothofagus forests of South America. Currently, the roseus group is composed by four species, while the vertebralis group consists of two. Phylogenetic relationships and species delimitation within each group are controversial. In fact, previous analyses considered that the roseus group was composed of between four to nine species. In this work, we evaluated phylogenetic relationships, diversification times, and species delimitation within the roseus group using a multi-locus dataset. For this purpose, mitochondrial (D-loop, Cyt b, and COI) and nuclear (POMC and CRYBA1) partial sequences from 164 individuals were amplified, representing all species. Maximum Likelihood (ML) and Bayesian approaches were used to reconstruct phylogenetic relationships. Species tree was estimated using BEAST and singular value decomposition scores for species quartets (SVDquartets). Species limits were evaluated with six coalescent approaches. Diversification times were estimated using mitochondrial and nuclear rates with LogNormal relaxed clock in BEAST. Nine well-supported monophyletic lineages were recovered in Bayesian, ML, and SVDquartets, including eight named species and a lineage composed by specimens from the Villarrica population (Bootstrap:>70, PP:> 0.99). Single-locus species delimitation analyses overestimated the species number in E. migueli, E. calcaratus, and E. roseus lineages, while multi-locus analyses recovered as species the nine lineages observed in phylogenetic analyses (Ctax = 0.69). It is hypothesized that Eupsophus diversification occurred during Mid-Pleistocene (0.42-0.14 Mya), with most species having originated after the Last Southern Patagonian Glaciation (0.18 Mya). Our results revitalize the hypothesis that the E. roseus group is composed of eight species and support the Villarrica lineage as a new putative species.

Population genetic analysis of the recently rediscovered Hula painted frog (Latonia nigriventer) reveals high genetic diversity and low inbreeding

After its recent rediscovery, the Hula painted frog (Latonia nigriventer) has remained one of the world's rarest and least understood amphibian species. Together with its apparently low dispersal capability and highly disturbed niche, the low abundance of this living fossil calls for urgent conservation measures. We used 18 newly developed microsatellite loci and four different models to calculate the effective population size (N_e) of a total of 125 Hula painted frog individuals sampled at a single location. We compare the N_e estimates to the estimates of potentially reproducing adults in this population (N_ad) determined through a capture-recapture study on 118 adult Hula painted frogs captured at the same site. Surprisingly, our data suggests that, despite N_ad estimates of only ~234-244 and N_e estimates of ~16.6-35.8, the species appears to maintain a very high genetic diversity (H_O = 0.771) and low inbreeding coefficient (F_IS = -0.018). This puzzling outcome could perhaps be explained by the hypotheses of either genetic rescue from one or more unknown Hula painted frog populations nearby or by recent admixture of genetically divergent subpopulations. Independent of which scenario is correct, the original locations of these populations still remain to be determined.

Testing skin swabbing for DNA sampling in dendrobatid frogs

Skin swabbing, a minimally invasive DNA sampling method recently proposed for adult amphibians, was tested on the dendrobatid frog Allobates femoralis. I compared DNA yield from skin swabs and toe clips by evaluating obtained DNA concentrations and purity of extracts, as well as amplification success using eleven polymorphic microsatellite loci. I also tested whether storing skin swabs for two months at −20°C affected the properties of the extract or microsatellite analysis. Results show that skin swabs of adult A. femoralis suffered from high contamination and yielded significantly lower DNA quality and quantity, resulting in insufficient genotyping success, than DNA obtained from toe clips. The relatively dry skin in dendrobatid frogs may have impeded the collection of sufficient viable cells, and the presence of skin alkaloids and microbiota in the frog mucus may lead to high contamination load of skin swabs.

When new human-modified habitats favour the expansion of an amphibian pioneer species: Evolutionary history of the natterjack toad (Bufo calamita) in a coal basin

Human activities affect microevolutionary dynamics by inducing environmental changes. In particular, land cover conversion and loss of native habitats decrease genetic diversity and jeopardize the adaptive ability of populations. Nonetheless, new anthropogenic habitats can also promote the successful establishment of emblematic pioneer species. We investigated this issue by examining the population genetic features and evolutionary history of the natterjack toad (Bufo [Epidalea] calamita) in northern France, where populations can be found in native coastal habitats and coalfield habitats shaped by European industrial history, along with an additional set of European populations located outside this focal area. We predicted contrasting patterns of genetic structure, with newly settled coalfield populations departing from migration-drift equilibrium. As expected, coalfield populations showed a mosaic of genetically divergent populations with short-range patterns of gene flow, and native coastal populations indicated an equilibrium state with an isolation-by-distance pattern suggestive of postglacial range expansion. However, coalfield populations exhibited (i) high levels of genetic diversity, (ii) no evidence of local inbreeding or reduced effective population size and (iii) multiple maternal mitochondrial lineages, a genetic footprint depicting independent colonization events. Furthermore, approximate Bayesian computations suggested several evolutionary trajectories from ancient isolation in glacial refugia during the Pleistocene, with biogeographical signatures of recent expansion probably confounded by human-mediated mixing of different lineages. From an evolutionary and conservation perspective, this study highlights the ecological value of industrial areas, provided that ongoing regional gene flow is ensured within the existing lineage boundaries.

Impaired gonadal and somatic development corroborate vulnerability differences to the synthetic estrogen ethinylestradiol among deeply diverged anuran lineages

Amphibians are undergoing a global decline. One poorly investigated reason could be the pollution of aquatic habitats by endocrine disrupting compounds (EDCs). We tested the susceptibility to the synthetically stabilized estrogen 17α-ethinylestradiol (EE2) in three deeply diverged anuran species, differing in sex determination systems, types of gonadogenesis and larval ecologies. To understand whether data from the amphibian model Xenopus laevis (Pipidae) are analogous and applicable to only distantly related non-model amphibians, tadpoles of X. laevis, Hyla arborea (Hylidae) and Bufo viridis (Bufonidae) were simultaneously exposed to 50, 500 and 5000ng/L EE2 from hatching until completion of metamorphosis, using a flow-through-system under identical experimental conditions. Comparing molecularly established genetic with histologically assessed phenotypic sex in all species, we have recently shown that EE2 provoked numerous genetic-male-to-phenotypic-female sex reversals and mixed sex individuals, confirming overall its expected feminizing effect. In the present study, we focus on the influence of EE2 on gonadal and somatic development. Anatomy and histology revealed several species-specific effects. In both non-model species, H. arborea and B. viridis, high numbers of anatomically impaired gonads were observed. In H. arborea, exposed to 5000ng/L EE2, numerous underdeveloped gonads were detected. Whereas EE2 did not alter snout-to-vent length and body weight of X. laevis metamorphs, H. arborea showed a treatment-dependent decrease, while B. viridis exhibited an increase in body weight and snout-to-vent length. Apart from a concentration-dependent occurrence of yellowish skin color in several H. arborea, no organ-specific effects were detected. Since EE2 ubiquitously occurs in many aquatic ecosystems and affects sexual and somatic development, among EDCs, it may indeed contribute to amphibian decline. The inter-species variation in developmental EE2-effects corroborates species-specific vulnerability differences towards EDCs between deeply diverged amphibian groups.

The genetic contribution to sex determination and number of sex chromosomes vary among populations of common frogs (Rana temporaria)

The patterns of sex determination and sex differentiation have been shown to differ among geographic populations of common frogs. Notably, the association between phenotypic sex and linkage group 2 (LG2) has been found to be perfect in a northern Swedish population, but weak and variable among families in a southern one. By analyzing these populations with markers from other linkage groups, we bring two new insights: (1) the variance in phenotypic sex not accounted for by LG2 in the southern population could not be assigned to genetic factors on other linkage groups, suggesting an epigenetic component to sex determination; (2) a second linkage group (LG7) was found to co-segregate with sex and LG2 in the northern population. Given the very short timeframe since post-glacial colonization (in the order of 1000 generations) and its seemingly localized distribution, this neo-sex chromosome system might be the youngest one described so far. It does not result from a fusion, but more likely from a reciprocal translocation between the original Y chromosome (LG2) and an autosome (LG7), causing their co-segregation during male meiosis. By generating a strict linkage between several important genes from the sex-determination cascade (Dmrt1, Amh and Amhr2), this neo-sex chromosome possibly contributes to the 'differentiated sex race' syndrome (strictly genetic sex determination and early gonadal development) that characterizes this northern population.

Sex-chromosome differentiation and 'sex races' in the common frog (Rana temporaria)

Sex-chromosome differentiation was recently shown to vary among common frog populations in Fennoscandia, suggesting a trend of increased differentiation with latitude. By rearing families from two contrasted populations (respectively, from northern and southern Sweden), we show this disparity to stem from differences in sex-determination mechanisms rather than in XY-recombination patterns. Offspring from the northern population display equal sex ratios at metamorphosis, with phenotypic sexes that correlate strongly with paternal LG2 haplotypes (the sex chromosome); accordingly, Y haplotypes are markedly differentiated, with male-specific alleles and depressed diversity testifying to their smaller effective population size. In the southern population, by contrast, a majority of juveniles present ovaries at metamorphosis; only later in development do sex ratios return to equilibrium. Even at these later stages, phenotypic sexes correlate only mildly with paternal LG2 haplotypes; accordingly, there are no recognizable Y haplotypes. These distinct patterns of gonadal development fit the concept of 'sex races' proposed in the 1930s, with our two populations assigned to the 'differentiated' and 'semi-differentiated' races, respectively. Our results support the suggestion that 'sex races' differ in the genetic versus epigenetic components of sex determination. Analysing populations from the 'undifferentiated race' with high-density genetic maps should help to further test this hypothesis.

Sex reversal assessments reveal different vulnerability to endocrine disruption between deeply diverged anuran lineages

Multiple anthropogenic stressors cause worldwide amphibian declines. Among several poorly investigated causes is global pollution of aquatic ecosystems with endocrine disrupting compounds (EDCs). These substances interfere with the endocrine system and can affect the sexual development of vertebrates including amphibians. We test the susceptibility to an environmentally relevant contraceptive, the artificial estrogen 17α-ethinylestradiol (EE2), simultaneously in three deeply divergent systematic anuran families, a model-species, Xenopus laevis (Pipidae), and two non-models, Hyla arborea (Hylidae) and Bufo viridis (Bufonidae). Our new approach combines synchronized tadpole exposure to three EE2-concentrations (50, 500, 5,000 ng/L) in a flow-through-system and pioneers genetic and histological sexing of metamorphs in non-model anurans for EDC-studies. This novel methodology reveals striking quantitative differences in genetic-male-to-phenotypic-female sex reversal in non-model vs. model species. Our findings qualify molecular sexing in EDC-analyses as requirement to identify sex reversals and state-of-the-art approaches as mandatory to detect species-specific vulnerabilities to EDCs in amphibians.

Empirical evidence for large X-effects in animals with undifferentiated sex chromosomes

Reproductive isolation is crucial for the process of speciation to progress. Sex chromosomes have been assigned a key role in driving reproductive isolation but empirical evidence from natural population processes has been restricted to organisms with degenerated sex chromosomes such as mammals and birds. Here we report restricted introgression at sex-linked compared to autosomal markers in a hybrid zone between two incipient species of European tree frog, Hyla arborea and H. orientalis, whose homologous X and Y sex chromosomes are undifferentiated. This large X-effect cannot result from the dominance or faster-X aspects of Haldane's rule, which are specific to degenerated sex chromosomes, but rather supports a role for faster-heterogametic-sex or faster-male evolutionary processes. Our data suggest a prominent contribution of undifferentiated sex chromosomes to speciation.

Timeframe of speciation inferred from secondary contact zones in the European tree frog radiation (Hyla arborea group)

Background

Hybridization between incipient species is expected to become progressively limited as their genetic divergence increases and reproductive isolation proceeds. Amphibian radiations and their secondary contact zones are useful models to infer the timeframes of speciation, but empirical data from natural systems remains extremely scarce. Here we follow this approach in the European radiation of tree frogs (Hyla arborea group). We investigated a natural hybrid zone between two lineages (Hyla arborea and Hyla orientalis) of Mio-Pliocene divergence (~5 My) for comparison with other hybrid systems from this group.

Results

We found concordant geographic distributions of nuclear and mitochondrial gene pools, and replicated narrow transitions (~30 km) across two independent transects, indicating an advanced state of reproductive isolation and potential local barriers to dispersal. This result parallels the situation between H. arborea and H. intermedia, which share the same amount of divergence with H. orientalis. In contrast, younger lineages show much stronger admixture at secondary contacts.

Conclusions

Our findings corroborate the negative relationship between hybridizability and divergence time in European tree frogs, where 5 My are necessary to achieve almost complete reproductive isolation. Speciation seems to progress homogeneously in this radiation, and might thus be driven by gradual genome-wide changes rather than single speciation genes. However, the timescale differs greatly from that of other well-studied amphibians. General assumptions on the time necessary for speciation based on evidence from unrelated taxa may thus be unreliable. In contrast, comparative hybrid zone analyses within single radiations such as our case study are useful to appreciate the advance of speciation in space and time.

Electronic supplementary material

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

Effects of Tail Clipping on Larval Performance and Tail Regeneration Rates in the Near Eastern Fire Salamander, Salamandra infraimmaculata

Tail-tip clipping is a common technique for collecting tissue samples from amphibian larvae and adults. Surprisingly, studies of this invasive sampling procedure or of natural tail clipping--i.e., bites inflicted by predators including conspecifics--on the performance and fitness of aquatic larval stages of urodeles are scarce. We conducted two studies in which we assessed the effects of posterior tail clipping (~30 percent of tail) on Near Eastern fire salamander (Salamandra infraimmaculata) larvae. In a laboratory study, we checked regeneration rates of posterior tail-tip clipping at different ages. Regeneration rates were hump-shaped, peaking at the age of ~30 days and then decreasing. This variation in tail regeneration rates suggests tradeoffs in resource allocation between regeneration and somatic growth during early and advanced development. In an outdoor artificial pond experiment, under constant larval densities, we assessed how tail clipping of newborn larvae affects survival to, time to, and size at metamorphosis. Repeated measures ANOVA on mean larval survival per pond revealed no effect of tail clipping. Tail clipping had correspondingly no effect on larval growth and development expressed in size (mass and snout-vent length) at, and time to, metamorphosis. We conclude that despite the given variation in tail regeneration rates throughout larval ontogeny, clipping of 30% percent of the posterior tail area seems to have no adverse effects on larval fitness and survival. We suggest that future use of this imperative tool for the study of amphibian should take into account larval developmental stage during the time of application and not just the relative size of the clipped tail sample.

Sex-Chromosome Homomorphy in Palearctic Tree Frogs Results from Both Turnovers and X-Y Recombination

Contrasting with birds and mammals, poikilothermic vertebrates often have homomorphic sex chromosomes, possibly resulting from high rates of sex-chromosome turnovers and/or occasional X-Y recombination. Strong support for the latter mechanism was provided by four species of European tree frogs, which inherited from a common ancestor (∼ 5 Ma) the same pair of homomorphic sex chromosomes (linkage group 1, LG1), harboring the candidate sex-determining gene Dmrt1. Here, we test sex linkage of LG1 across six additional species of the Eurasian Hyla radiation with divergence times ranging from 6 to 40 Ma. LG1 turns out to be sex linked in six of nine resolved cases. Mapping the patterns of sex linkage to the Hyla phylogeny reveals several transitions in sex-determination systems within the last 10 My, including one switch in heterogamety. Phylogenetic trees of DNA sequences along LG1 are consistent with occasional X-Y recombination in all species where LG1 is sex linked. These patterns argue against one of the main potential causes for turnovers, namely the accumulation of deleterious mutations on nonrecombining chromosomes. Sibship analyses show that LG1 recombination is strongly reduced in males from most species investigated, including some in which it is autosomal. Intrinsically low male recombination might facilitate the evolution of male heterogamety, and the presence of important genes from the sex-determination cascade might predispose LG1 to become a sex chromosome.

Origin and genome evolution of polyploid green toads in Central Asia: evidence from microsatellite markers

Polyploidization, which is expected to trigger major genomic reorganizations, occurs much less commonly in animals than in plants, possibly because of constraints imposed by sex-determination systems. We investigated the origins and consequences of allopolyploidization in Palearctic green toads (Bufo viridis subgroup) from Central Asia, with three ploidy levels and different modes of genome transmission (sexual versus clonal), to (i) establish a topology for the reticulate phylogeny in a species-rich radiation involving several closely related lineages and (ii) explore processes of genomic reorganization that may follow polyploidization. Sibship analyses based on 30 cross-amplifying microsatellite markers substantiated the maternal origins and revealed the paternal origins and relationships of subgenomes in allopolyploids. Analyses of the synteny of linkage groups identified three markers affected by translocation events, which occurred only within the paternally inherited subgenomes of allopolyploid toads and exclusively affected the linkage group that determines sex in several diploid species of the green toad radiation. Recombination rates did not differ between diploid and polyploid toad species, and were overall much reduced in males, independent of linkage group and ploidy levels. Clonally transmitted subgenomes in allotriploid toads provided support for strong genetic drift, presumably resulting from recombination arrest. The Palearctic green toad radiation seems to offer unique opportunities to investigate the consequences of polyploidization and clonal transmission on the dynamics of genomes in vertebrates.

Beam Brook revisited: a molecular study of a historically introduced non-native amphibian (Triturus carnifex) and its limited introgression into native UK Triturus cristatus populations

Human mediated introductions of non-native species can pose a major threat to global biodiversity on several accounts i.e. through competition, the introduction of novel pathogens, and genetic pollution. Where hybridization occurs between two closely related species the F1 offspring are usually phenotypically discernible whereas F2 hybrid, backcrossed or admixed individuals become more difficult to separate. At this point the utilization of molecular methods is required in conservation efforts to differentiate and manage populations. This study demonstrated how a possible threat of hybridization from an introduced non-native (T. carnifex) with a protected native newt species (T. cristatus) could be investigated with molecular tools, and examined the current extent of its genetic introgression over an 80 years period. The results confirmed that hybridization had taken place at the site of introduction (and continues to do so), and that historically limited local dispersal of both non-natives and/or hybrids had occurred sometime in the past. However, the data suggests that although dispersal of hybrids into a local satellite site may still be occuring, hybridization with native species appears limited.

Sex-chromosome differentiation parallels postglacial range expansion in European tree frogs (Hyla arborea)

Occasional XY recombination is a proposed explanation for the sex-chromosome homomorphy in European tree frogs. Numerous laboratory crosses, however, failed to detect any event of male recombination, and a detailed survey of NW-European Hyla arborea populations identified male-specific alleles at sex-linked loci, pointing to the absence of XY recombination in their recent history. Here, we address this paradox in a phylogeographic framework by genotyping sex-linked microsatellite markers in populations and sibships from the entire species range. Contrasting with postglacial populations of NW Europe, which display complete absence of XY recombination and strong sex-chromosome differentiation, refugial populations of the southern Balkans and Adriatic coast show limited XY recombination and large overlaps in allele frequencies. Geographically and historically intermediate populations of the Pannonian Basin show intermediate patterns of XY differentiation. Even in populations where X and Y occasionally recombine, the genetic diversity of Y haplotypes is reduced below the levels expected from the fourfold drop in copy numbers. This study is the first in which X and Y haplotypes could be phased over the distribution range in a species with homomorphic sex chromosomes; it shows that XY-recombination patterns may differ strikingly between conspecific populations, and that recombination arrest may evolve rapidly (<5000 generations).

Phenotypic divergence of the common toad (Bufo bufo) along an altitudinal gradient: evidence for local adaptation

Variation in the environment can induce different patterns of genetic and phenotypic differentiation among populations. Both neutral processes and selection can influence phenotypic differentiation. Altitudinal phenotypic variation is of particular interest in disentangling the interplay between neutral processes and selection in the dynamics of local adaptation processes but remains little explored. We conducted a common garden experiment to study the phenotypic divergence in larval life-history traits among nine populations of the common toad (Bufo bufo) along an altitudinal gradient in France. We further used correlation among population pairwise estimates of quantitative trait (QST) and neutral genetic divergence (FST from neutral microsatellite markers), as well as altitudinal difference, to estimate the relative role of divergent selection and neutral genetic processes in phenotypic divergence. We provided evidence for a neutral genetic differentiation resulting from both isolation by distance and difference in altitude. We found evidence for phenotypic divergence along the altitudinal gradient (faster development, lower growth rate and smaller metamorphic size). The correlation between pairwise QSTs-FSTs and altitude differences suggested that this phenotypic differentiation was most likely driven by altitude-mediated selection rather than by neutral genetic processes. Moreover, we found different divergence patterns for larval traits, suggesting that different selective agents may act on these traits and/or selection on one trait may constrain the evolution on another through genetic correlation. Our study highlighted the need to design more integrative studies on the common toad to unravel the underlying processes of phenotypic divergence and its selective agents in the context of environmental clines.

Optimizing the trade-off between spatial and genetic sampling efforts in patchy populations: towards a better assessment of functional connectivity using an individual-based sampling scheme

Genetic data are increasingly used in landscape ecology for the indirect assessment of functional connectivity, that is, the permeability of landscape to movements of organisms. Among available tools, matrix correlation analyses (e.g. Mantel tests or mixed models) are commonly used to test for the relationship between pairwise genetic distances and movement costs incurred by dispersing individuals. When organisms are spatially clustered, a population-based sampling scheme (PSS) is usually performed, so that a large number of genotypes can be used to compute pairwise genetic distances on the basis of allelic frequencies. Because of financial constraints, this kind of sampling scheme implies a drastic reduction in the number of sampled aggregates, thereby reducing sampling coverage at the landscape level. We used matrix correlation analyses on simulated and empirical genetic data sets to investigate the efficiency of an individual-based sampling scheme (ISS) in detecting isolation-by-distance and isolation-by-barrier patterns. Provided that pseudo-replication issues are taken into account (e.g. through restricted permutations in Mantel tests), we showed that the use of interindividual measures of genotypic dissimilarity may efficiently replace interpopulation measures of genetic differentiation: the sampling of only three or four individuals per aggregate may be sufficient to efficiently detect specific genetic patterns in most situations. The ISS proved to be a promising methodological alternative to the more conventional PSS, offering much flexibility in the spatial design of sampling schemes and ensuring an optimal representativeness of landscape heterogeneity in data, with few aggregates left unsampled. Each strategy offering specific advantages, a combined use of both sampling schemes is discussed.

Isolation and characterization of microsatellite loci for the European tree frog (Hyla arborea)

We developed 11 new microsatellite markers for the European tree frog (Hyla arborea), and tested patterns of polymorphism in 54 adults (27 males and 27 females) from two ponds close to Lausanne (Western Switzerland). One marker was sex linked and two pairs displayed linkage disequilibrium. Comparisons of allele numbers with heterozygosity values support a stepwise-mutation model at neutral equilibrium, with mutation rates spanning nearly two orders of magnitude. These markers will prove useful for population genetic studies and fine-scale investigations requiring genetic assignment techniques.