From species divergence to population structure: a multimarker approach on the most basal lineage of Salamandridae, the spectacled salamanders (genus Salamandrina) from Italy.

The Apennine Peninsula is one of Europe's main glacial refugial areas and harbors a large number of lineages and species. Here, a pattern of higher genetic diversity in the south compared to that of the north is characteristic of most vertebrates; however, most studies that have produced these results have relied only on inferences based on mitochondrial DNA. The spectacled salamanders (genus Salamandrina) are endemic to the Apennine Peninsula and have diverged into two sibling species: S. terdigitata (in the south) and S. perspicillata (in the north), presumably in the late Miocene or early Pliocene. By sequencing one mitochondrial (cytb) and two nuclear genes (RAG1 and POMC) and genotyping 10 microsatellite loci, we traced the evolution of these sibling species from their divergence to their contemporary population structure at a fine scale. Using a multilocus coalescent-based approach, we estimated the temporal divergence of both species at approximately 2.25 mya (million years ago), which, hence, is much younger than previous estimates. The classical pattern of high genetic diversity in the south and lower diversity in the north was confirmed only for some markers, and the demographic histories of the two species differed substantially. Whereas S. perspicillata (north) expanded from a single major refugium in the center of the Apennine Peninsula, populations of S. terdigitata (south) persisted through cooler periods in multiple refugia. Further, the fine-scale population genetic structure of 16 S. perspicillata populations revealed significant genetic differentiation, even across short geographic distances. The results of our study stress that for a better understanding of phylogeographic patterns and past demographic processes, both mitochondrial and multiple nuclear loci should be analyzed to avoid gene-specific, and possibly biased results.

Nuclear and mitochondrial multilocus phylogeny and survey of alkaloid content in true salamanders of the genus Salamandra (Salamandridae).

The genus Salamandra represents a clade of six species of Palearctic salamanders of either contrasted black-yellow, or uniformly black coloration, known to contain steroidal alkaloid toxins in high concentrations in their skin secretions. This study reconstructs the phylogeny of the genus Salamandra based on DNA sequences of segments of 10 mitochondrial and 13 nuclear genes from 31 individual samples representing all Salamandra species and most of the commonly recognized subspecies. The concatenated analysis of the complete dataset produced a fully resolved tree with most nodes strongly supported, suggesting that a clade composed of the Alpine salamander (S. atra) and the Corsican fire salamander (S. corsica) is the sister taxon to a clade containing the remaining species, among which S. algira and S. salamandra are sister species. Separate analyses of mitochondrial and nuclear data partitions disagreed regarding basal nodes and in the position of the root but concordantly recovered the S. atra/S. corsica as well as the S. salamandra/S. algira relationship. A species-tree analysis suggested almost simultaneous temporal splits between these pairs of species, which we hypothesize was caused by vicariance events after the Messinian salinity crisis (from late Miocene to early Pliocene). A survey of toxins with combined gas chromatography/mass spectroscopy confirmed the presence of samandarine and/or samandarone steroidal alkaloids in all species of Salamandra as well as in representatives of their sister group, Lyciasalamandra. Samandarone was also detected in lower concentrations in other salamandrids including Calotriton, Euproctus, Lissotriton, and Triturus, suggesting that the presence and possible biosynthesis of this alkaloid is plesiomorphic within the Salamandridae.

Radically different phylogeographies and patterns of genetic variation in two European brown frogs, genus Rana.

We reconstruct range-wide phylogeographies of two widespread and largely co-occurring Western Palearctic frogs, Rana temporaria and R. dalmatina. Based on tissue or saliva samples of over 1000 individuals, we compare a variety of genetic marker systems, including mitochondrial DNA, single-copy protein-coding nuclear genes, microsatellite loci, and single nucleotide polymorphisms (SNPs) of transcriptomes of both species. The two focal species differ radically in their phylogeographic structure, with R. temporaria being strongly variable among and within populations, and R. dalmatina homogeneous across Europe with a single strongly differentiated population in southern Italy. These differences were observed across the various markers studied, including microsatellites and SNP density, but especially in protein-coding nuclear genes where R. dalmatina had extremely low heterozygosity values across its range, including potential refugial areas. On the contrary, R. temporaria had comparably high range-wide values, including many areas of probable postglacial colonization. A phylogeny of R. temporaria based on various concatenated mtDNA genes revealed that two haplotype clades endemic to Iberia form a paraphyletic group at the base of the cladogram, and all other haplotypes form a monophyletic group, in agreement with an Iberian origin of the species. Demographic analysis suggests that R. temporaria and R. dalmatina have genealogies of roughly the same time to coalescence (TMRCA ~3.5 mya for both species), but R. temporaria might have been characterized by larger ancestral and current effective population sizes than R. dalmatina. The high genetic variation in R. temporaria can therefore be explained by its early range expansion out of Iberia, with subsequent cycles of differentiation in cryptic glacial refugial areas followed by admixture, while the range expansion of R. dalmatina into central Europe is a probably more recent event.

A multigene species tree for Western Mediterranean painted frogs (Discoglossus).

Painted frogs (Discoglossus) are an anuran clade that originated in the Upper Miocene. Extant species are morphologically similar and have a circum-Mediterranean distribution. We assembled a multilocus dataset from seven nuclear and four mitochondrial genes for several individuals of all but one of the extant species and reconstructed a robust phylogeny by applying a coalescent-based species-tree method and a concatenation approach, both of which gave congruent results. The earliest phylogenetic split within Discoglossus separates D. montalentii from a clade comprising all other species. Discoglossus montalentii is monophyletic for haplotype variation at all loci and has distinct morphological, bioacoustic and karyotypic characters. We find moderate support for a sister-group relationship between the Iberian taxa and the Moroccan D. scovazzi, and high support for a D. pictus -D. sardus clade distributed around the Tyrrhenian basin. Topological discordance among gene trees during the speciation of D. galganoi, D. scovazzi, D. pictus and D. sardus is interpreted as the consequence of nearly simultaneous, vicariant diversification. The timing of these events is unclear, but possibly coincided with the final geotectonic rearrangement of the Western Mediterranean in the Middle Miocene or later during the Messinian salinity crisis. The Iberian taxa D. galganoi galganoi and D. g. jeanneae are reciprocally monophyletic in mitochondrial DNA but not in nuclear gene trees, and are therefore treated as subspecies of D. galganoi.

Mitochondrial Introgression, Color Pattern Variation, and Severe Demographic Bottlenecks in Three Species of Malagasy Poison Frogs, Genus Mantella.

Madagascar is a biodiversity hotspot particularly rich in amphibian diversity and only a few charismatic Malagasy amphibians have been investigated for their population-level differentiation. The Mantellamadagascariensis group is composed of two rainforest and three swamp forest species of poison frogs. We first confirm the monophyly of this clade using DNA sequences of three nuclear and four mitochondrial genes, and subsequently investigate the population genetic differentiation and demography of the swamp forest species using one mitochondrial, two nuclear and a set of nine microsatellite markers. Our results confirm the occurrence of two main mitochondrial lineages, one dominated by Mantellaaurantiaca (a grouping supported also by our microsatellite-based tree) and the other by Mantellacrocea + Mantellamilotympanum. These two main lineages probably reflect an older divergence in swamp Mantella. Widespread mitochondrial introgression suggests a fairly common occurrence of inter-lineage gene flow. However, nuclear admixture seems to play only a limited role in this group, and the analyses of the RAG-1 marker points to a predominant incomplete lineage sorting scenario between all five species of the group, which probably diverged relatively recently. Our demographic analyses show a common, severe and recent demographic contraction, inferred to be in temporal coincidence with the massive deforestation events that took place in the past 1000 years. Current data do not allow to conclusively delimit independent evolutionary units in these frogs, and we therefore refrain to suggest any taxonomic changes.

Turtle soup, Prohibition, and the population genetic structure of Diamondback Terrapins (Malaclemys terrapin).

Diamondback terrapins (Malaclemys terrapin) were a popular food item in early twentieth century America, and were consumed in soup with sherry. Intense market demand for terrapin meat resulted in population declines, notably along the Atlantic seaboard. Efforts to supply terrapins to markets resulted in translocation events, as individuals were moved about to stock terrapin farms. However, in 1920 the market for turtle soup buckled with the enactment of the eighteenth amendment to the United States' Constitution-which initiated the prohibition of alcoholic drinks-and many terrapin fisheries dumped their stocks into local waters. We used microsatellite data to show that patterns of genetic diversity along the terrapin's coastal range are consistent with historical accounts of translocation and cultivation activities. We identified possible instances of human-mediated dispersal by estimating gene flow over historical and contemporary timescales, Bayesian model testing, and bottleneck tests. We recovered six genotypic clusters along the Gulf and Atlantic coasts with varying degrees of admixture, including increased contemporary gene flow from Texas to South Carolina, from North Carolina to Maryland, and from North Carolina to New York. In addition, Bayesian models incorporating translocation events outperformed stepping-stone models. Finally, we were unable to detect population bottlenecks, possibly due to translocation reintroducing genetic diversity into bottlenecked populations. Our data suggest that current patterns of genetic diversity in the terrapin were altered by the demand for turtle soup followed by the enactment of alcohol prohibition. In addition, our study shows that population genetic tools can elucidate metapopulation dynamics in taxa with complex genetic histories impacted by anthropogenic activities.

A simplified molecular method for distinguishing among species and ploidy levels in European water frogs (Pelophylax).

Western Palearctic water frogs in the genus Pelophylax are a set of morphologically similar anuran species that form hybridogenetic complexes. Fully reliable identification of species and especially of hybrid ploidy depends on karyological and molecular methods. In central Europe, native water frog populations consist of the Pelophylax esculentus complex, that is, P. lessonae (LL), P. ridibundus (RR) and the hybrid form P. esculentus that can have different karyotypes (RL, LLR and RRL). We developed existing molecular methods further and propose a simple PCR method based on size-differences in the length of the serum albumin intron-1 and the RanaCR1, a non-LTR retrotransposon of the chicken repeat (CR) family. This PCR yields taxon-specific banding patterns that can easily be screened by standard agarose gel electrophoresis and correctly identify species in all of the 160 samples that had been identified to karyotype with other methods. To distinguish ploidy levels in LR, LLR and RRL specimens, we used the ratio of the peak heights of the larger (ridibundus specific) to the smaller (lessonae specific) bands of fluorescently labelled PCR products resolved on a capillary DNA sequencer and obtained a correct assignment of the karyotype in 93% of cases. Our new method will cut down time and expenses drastically for a reliable identification of water frogs of the P. esculentus complex and potentially for identification of other hybridogenetic complexes and/or taxa, and it even serves as a good indicator of the ploidy status of hybrid individuals.

Molecular characterization of MHC class II in a nonmodel anuran species, the fire-bellied toad Bombina bombina.

While the anuran Xenopus comprises one of the best characterized nonmammalian taxa regarding the major histocompatibility complex (MHC), the organization of this gene complex has never been studied in other anurans, and information on amphibian MHC (other than Xenopus) is generally very scarce. Here, we describe the characterization of the first MHC class II B cDNA sequences from a nonmodel anuran species, the European fire-bellied toad (Bombina bombina). We isolated two transcript sequences differing substantially in amino acid composition and length within the beta2 domain. To investigate the variability of the peptide binding region in this species, we sequenced a 158-bp large fragment from wild B. bombina (n = 20) and identified eight distinct alleles. All substitutions but one were nonsynonymous, and many of the highly polymorphic sites corresponded with amino acid positions known to be involved in antigen binding. The level of variation we found in B. bombina was similar compared to that previously found in a comparable sample of a wild urodelan species, Ambystoma tigrinum, and to that found in Xenopus laevis. Based on the cDNA data and the individual's allelic diversity, we conclude that Bombina possesses at least two class II B loci. With our new beta1 primers, we were able to generate sequences in other species of anurans. We provide here a first phylogenetic analysis of this gene in amphibians.

Population genetics of the diamondback terrapin (Malaclemys terrapin).

We examined the population genetic structure of the diamondback terrapins (Malaclemys terrapin), within and among estuaries. Based on mark-recapture studies, these estuarine turtles have high site fidelity that is likely to make them vulnerable to local extinctions. We tested if observed site fidelity of adults would be reflected in intraestuarine population genetic structure of six highly polymorphic microsatellite loci (five tetranucleotide and one dinucleotide). No evidence was found for population structuring within the Charleston estuary nor among three different estuaries in South Carolina. We then examined four other terrapin populations from North Carolina to New York, as well as from the Florida Keys and from Texas. With increasing geographical distance, genetic differentiation increased from South Carolina through New York, but overall values were low. The dinucleotide locus contributed significantly more to the genetic differentiation of some population comparisons than any of the other loci. Interestingly, terrapins from South Carolina to New York were much more genetically similar to those from Texas (rho = 0.154) than to those from Florida (rho = 0.357). We attribute this pattern to extensive translocations of terrapins during the early 20th century to replenish diminished populations and to provide turtle farms with stocks. Terrapins collected in Texas were especially sought for shipment to the northeastern US because of their larger size. Our study indicates no population structure within or among adjacent estuaries. Thus, the mark-recapture information from adult and subadult feeding locations is a poor predictor of population genetic structure. Additionally, it appears that past human activities may have drastically altered the genetics of current populations. Finally, our data suggest that translocation of eggs or head starting of terrapins within estuaries or among adjacent estuaries is acceptable from a genetic standpoint.