The complete mitochondrial genome of the Critically Endangered Saba Green Iguana, Iguana iguana (Squamata: Iguanidae).

The populations of native iguanas in the Caribbean Lesser Antilles are threatened by the wide occurrence and spread of non-native iguanas. Until recently, competitive hybridization was not believed to threaten the Saba Green Iguana, a subpopulation of Iguana iguana (Linnaeus, 1758) from the island of Saba. However, the arrival of non-native iguanas has put the native population at risk, leading to a change in the conservation status of the Saba Green Iguana to Critically Endangered, according to guidelines from the International Union for the Conservation of Nature. Here, we generated the complete mitogenome of the Saba Green Iguana using Oxford Nanopore long-read technology. The mitogenome is 16,626 bp long and has 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes, and a control region (1194 bp). Noteworthy, this is only the second published mitogenome for the Iguana iguana species complex, despite the known high intraspecific genetic variation.

Phylogenomic inference of species and subspecies diversity in the Palearctic salamander genus Salamandra.

The salamander genus Salamandra is widespread across Europe, North Africa, and the Near East and is renowned for its conspicuous and polymorphic colouration and diversity of reproductive modes. The phylogenetic relationships within the genus, and especially in the highly polymorphic species S. salamandra, have been very challenging to elucidate, leaving its real evolutionary history and classification at species and subspecies levels a topic of debate and contention. However, the distribution of diversity and species delimitation within the genus are critically important for identifying evolutionarily significant units for conservation and management, especially in light of threats posed by the pathogenic chytrid fungus Batrachochytrium salamandrivorans that is causing massive declines of S. salamandra populations in central Europe. Here, we conducted a phylogenomic analysis from across the taxonomic and geographic breadth of the genus Salamandra in its entire range. Bayesian, maximum likelihood and network-based phylogenetic analyses of up to 4905 ddRADseq-loci (294,300 nucleotides of sequence) supported the distinctiveness of all currently recognised species (Salamandra algira, S. atra, S. corsica, S. infraimmaculata, S. lanzai, and S. salamandra), and all five species for which we have multiple exemplars were confirmed as monophyletic. Within S. salamandra, two main clades can be distinguished: one clade with the Apenninic subspecies S. s. gigliolii nested within the Iberian S. s. bernardezi/fastuosa; and a second clade comprising all other Iberian, Central and East European subspecies. Our analyses revealed that some of the currently recognized subspecies of S. salamandra are paraphyletic and may require taxonomic revision, with the Central- and Eastern-European subspecies all being poorly differentiated at the analysed genomic markers. Salamandra s. longirostris - sometimes considered a separate species - was nested within S. salamandra, consistent with its subspecies status. The relationships identified within and between Salamandra species provide valuable context for future systematic and biogeographic studies, and help elucidate critical evolutionary units for conservation and taxonomy.

Hidden in plain sight: a new frog species of the genus Blommersia from the oceanic island of Mayotte, Comoros archipelago.

The amphibian fauna of the western Indian ocean volcanic island of Mayotte is currently constituted by two species belonging to two genera of the anuran family Mantellidae: Blommersia transmarina and Boophis nauticus. These were recently described after intense fieldwork on the herpetofauna of the island. We here describe a third new species of frog from Mayotte, based on morphological and molecular data, that occurs in sympatry with the others and was utterly unnoticed until now. Genetic analyses of the16S rRNA gene, including all described and several undescribed species of the genus Blommersia from Madagascar and Mayotte, confirms that the new species is the sister species of Blommersia transmarina. Both species show apparent morphological differences as well as different life histories, ecology and genetics that confirm Blommersia nataliae sp. nov. as a new species. We propose an IUCN Red List status of Critically Endangered for B. nataliae sp. nov.

ACDC, a global database of amphibian cytochrome-b sequences using reproducible curation for GenBank records.

Genetic data are a crucial and exponentially growing resource across all biological sciences, yet curated databases are scarce. The widespread occurrence of sequence and (meta)data errors in public repositories calls for comprehensive improvements of curation protocols leading to robust research and downstream analyses. We collated and curated all available GenBank cytochrome-b sequences for amphibians, a benchmark marker in this globally declining vertebrate clade. The Amphibia's Curated Database of Cytochrome-b (ACDC) consists of 36,514 sequences representing 2,309 species from 398 genera (median = 2 with 50% interquartile ranges of 1-7 species/genus). We updated the taxonomic identity of >4,800 sequences (ca. 13%) and found 2,359 (6%) conflicting sequences with 84% of the errors originating from taxonomic misidentifications. The database (accessible at https://doi.org/10.6084/m9.figshare.9944759 ) also includes an R script to replicate our study for other loci and taxonomic groups. We provide recommendations to improve genetic-data quality in public repositories and flag species for which there is a need for taxonomic refinement in the face of increased rate of amphibian extinctions in the Anthropocene.

Functional colour genes and signals of selection in colour-polymorphic salamanders.

Coloration has been associated with multiple biologically relevant traits that drive adaptation and diversification in many taxa. However, despite the great diversity of colour patterns present in amphibians the underlying molecular basis is largely unknown. Here, we use insight from a highly colour-variable lineage of the European fire salamander (Salamandra salamandra bernardezi) to identify functional associations with striking variation in colour morph and pattern. The three focal colour morphs-ancestral black-yellow striped, fully yellow and fully brown-differed in pattern, visible coloration and cellular composition. From population genomic analyses of up to 4,702 loci, we found no correlations of neutral population genetic structure with colour morph. However, we identified 21 loci with genotype-phenotype associations, several of which relate to known colour genes. Furthermore, we inferred response to selection at up to 142 loci between the colour morphs, again including several that relate to coloration genes. By transcriptomic analysis across all different combinations, we found 196 differentially expressed genes between yellow, brown and black skin, 63 of which are candidate genes involved in animal coloration. The concordance across different statistical approaches and 'omic data sets provide several lines of evidence for loci linked to functional differences between colour morphs, including TYR, CAMK1 and PMEL. We found little association between colour morph and the metabolomic profile of its toxic compounds from the skin secretions. Our research suggests that current ecological and evolutionary hypotheses for the origins and maintenance of these striking colour morphs may need to be revisited.

Towards completion of the species inventory of small-sized leaf-tailed geckos: two new species of Uroplatus from northern Madagascar.

The northern part of Madagascar is well known for its high species diversity and endemism. Exceptional species richness is related to the existence of large forest blocks and mountain complexes. These areas shelter a diverse variety of habitats occupied by a wide diversity of species, including leaf-tailed geckos of the genus Uroplatus. Based on morphological and molecular evidence, we here formally name two evolutionary lineages as new species that previously had been considered as candidate species (Uroplatus spp. Ca3 and Ca4), both small-sized species of the Uroplatus ebenaui group. Genetically, both new species are related to U. finiavana with a genetic divergence (uncorrected pairwise distance) of 10.3-12.8% in a fragment of the mitochondrial 16S rRNA gene, and separated from each other by 10.3-10.7%. Uroplatus fangorn sp. nov. is described from the Marojejy massif and is also known from Andrevorevo, Lohanandroranga and Sorata; it is similar to U. ebenaui and U. fetsy but is distinguishable by its rather short tail and only partially black pigmented oral mucosa. Uroplatus fivehy sp. nov. is described from the Sorata Massif and is wider-ranging, occurring in an area from Marotandrano and Makira to Sorata, comprising Marojejy, Anjanaharibe-Sud, Ankitsika and Betaolana; it is morphologically similar to U. finiavana but distinguishable by tail size and shape. Both newly described species are found in rainforest from mid to high elevation and range respectively from 840-1417 m for U. fivehy sp nov. and 1300-1800 m for U. fangorn sp. nov. According to their respective geographical distributions, we propose to classify the two new species as Vulnerable under the IUCN Red List criteria, due to their occurrence in some forests outside of the protected area network plus the continuing decline of forest patches in the north of Madagascar.

Challenges in estimating ancestral state reconstructions: the evolution of migration in Sylvia warblers as a study case.

Our current understanding of how species have evolved is mainly based on comparative phylogenetic methods, which use phylogenies to infer the evolution of traits. The development of ancestral state reconstruction (ASR) methods has provided the tools to reconstruct trait evolution, which are widely used in fields like evolutionary biology, macroecology and paleontology. As there are different elements involved in those analyses, with different levels of uncertainty (i.e. relating to branch length estimation, trait coding, statistical framework, taxon sampling or software), the various combinations of these elements likely have a strong impact on the reconstruction of the evolution of traits, potentially leading to opposite conclusions. To assess the impact of these different elements in ASR, we performed a set of analyses, including all possible combinations of such elements and using the evolution of migratory behavior in Sylvia warblers as a case study, which was coded as a continuous or as a discrete character. Our results show that taxon sampling, character coding, tree shape, statistical framework and software all significantly affect ASR, both individually and in combination. Not all reconstructed tree nodes show discrepancies, but in the critical ones most pairwise comparisons between analyses lead to conflicting and unexpectedly antagonistic results (zero migration vs fully migratory), thus challenging interpretations of trait evolution. We propose some possible solutions to partly inform decisions, involving the method selection and the incorporation of biological or fossil evidence regarding how traits evolve, but our results demand serious rethinking about how the research community currently uses ASR.

Intraspecific variation in lizard heat tolerance alters estimates of climate impact.

Research addressing the effects of global warming on the distribution and persistence of species generally assumes that population variation in thermal tolerance is spatially constant or overridden by interspecific variation. Typically, this rationale is implicit in sourcing one critical thermal maximum (CTmax ) population estimate per species to model spatiotemporal cross-taxa variation in heat tolerance. Theory suggests that such an approach could result in biased or imprecise estimates and forecasts of impact from climate warming, but limited empirical evidence in support of those expectations exists. We experimentally quantify the magnitude of intraspecific variation in CTmax among lizard populations, and the extent to which incorporating such variability can alter estimates of climate impact through a biophysical model. To do so, we measured CTmax from 59 populations of 15 Iberian lizard species (304 individuals). The overall median CTmax across all individuals from all species was 42.8°C and ranged from 40.5 to 48.3°C, with species medians decreasing through xeric, climate-generalist and mesic taxa. We found strong statistical support for intraspecific differentiation in CTmax by up to a median of 3°C among populations. We show that annual restricted activity (operative temperature > CTmax ) over the Iberian distribution of our study species differs by a median of >80 hr per 25-km2 grid cell based on different population-level CTmax estimates. This discrepancy leads to predictions of spatial variation in annual restricted activity to change by more than 20 days for six of the study species. Considering that during restriction periods, reptiles should be unable to feed and reproduce, current projections of climate-change impacts on the fitness of ectotherm fauna could be under- or over-estimated depending on which population is chosen to represent the physiological spectra of the species in question. Mapping heat tolerance over the full geographical ranges of single species is thus critical to address cross-taxa patterns and drivers of heat tolerance in a biologically comprehensive way.

Two new species of leaf-tailed geckos (Uroplatus) from the Tsaratanana mountain massif in northern Madagascar.

The Tsaratanana Massif is the highest mountain massif of Madagascar and is characterized by a high species-level endemism of its biota. Here we describe two new small-sized species of leaf-tailed geckos of the Uroplatus ebenaui group from this region. Named in a preliminary way as confirmed candidate species Uroplatus ebenaui [Ca1] and [Ca2] in previous studies, we here provide detailed data on their molecular and morphological differentiation and describe them as Uroplatus fotsivava sp. nov. and Uroplatus kelirambo sp. nov. Their closest relatives are U. fiera and U. finiavana, respectively, as revealed by a phylogeny based on DNA sequences of four mitochondrial genes. However, integration of various lines of evidence confirms that the new species are independent evolutionary lineages, differing from other Uroplatus by high mitochondrial divergences, lack of haplotype sharing in the nuclear CMOS gene, and shape and relative size of the tail. While U. fotsivava is known from multiple sites at elevations between 1100-1538 m above sea level, U. kelirambo was found only at high elevations between 2000-2200 m a.s.l.; it therefore is the only known forest-dwelling nocturnal gecko reaching such high elevations in Madagascar, and might be narrowly endemic to montane forests of the Tsaratanana Massif.

Spatial Biodiversity Patterns of Madagascar's Amphibians and Reptiles.

Madagascar has become a model region for testing hypotheses of species diversification and biogeography, and many studies have focused on its diverse and highly endemic herpetofauna. Here we combine species distribution models of a near-complete set of species of reptiles and amphibians known from the island with body size data and a tabulation of herpetofaunal communities from field surveys, compiled up to 2008. Though taxonomic revisions and novel distributional records arose since compilation, we are confident that the data are appropriate for inferring and comparing biogeographic patterns among these groups of organisms. We observed species richness of both amphibians and reptiles was highest in the humid rainforest biome of eastern Madagascar, but reptiles also show areas of high richness in the dry and subarid western biomes. In several amphibian subclades, especially within the Mantellidae, species richness peaks in the central eastern geographic regions while in reptiles different subclades differ distinctly in their richness centers. A high proportion of clades and subclades of both amphibians and reptiles have a peak of local endemism in the topographically and bioclimatically diverse northern geographic regions. This northern area is roughly delimited by a diagonal spanning from 15.5°S on the east coast to ca. 15.0°S on the west coast. Amphibian diversity is highest at altitudes between 800-1200 m above sea-level whereas reptiles have their highest richness at low elevations, probably reflecting the comparatively large number of species specialized to the extended low-elevation areas in the dry and subarid biomes. We found that the range sizes of both amphibians and reptiles strongly correlated with body size, and differences between the two groups are explained by the larger body sizes of reptiles. However, snakes have larger range sizes than lizards which cannot be readily explained by their larger body sizes alone. Range filling, i.e., the amount of suitable habitat occupied by a species, is less expressed in amphibians than in reptiles, possibly reflecting their lower dispersal capacity. Taxonomic composition of communities assessed by field surveys is largely explained by bioclimatic regions, with communities from the dry and especially subarid biomes distinctly differing from humid and subhumid biomes.

The complete mitochondrial genome of the Endangered European brown frog Rana pyrenaica through RNAseq.

We sequenced the complete mitogenome of the Pyrenean frog Rana pyrenaica, which was determined from an Illumina Hi-seq RNAseq run (Illumina Inc., San Diego, CA). The genome is 17,213 bp in size, including 13 protein-coding genes, 21 transfer RNAs, two ribosomal RNAs and a control region. It shows the typical gene order of previously available frog mitogenomes, although it lacks the tRNAPhe . This is the first complete mitogenome described for a Western Palearctic brown frog species.

Individual and Geographic Variation of Skin Alkaloids in Three Swamp-Forest Species of Madagascan Poison Frogs (Mantella).

Seventy skins of three mantellid frog species from Madagascan swamp-forest habitats, Mantella aurantiaca, M. crocea, and M. milotympanum, were individually examined for skin alkaloids using GC/MS. These poison frogs were found to differ significantly in their alkaloid composition from species of Mantella originating from non-flooded rainforest in eastern Madagascar, which were examined in earlier work. Only 16 of the previously detected 106 alkaloids were represented among the 60 alkaloids from the swamp-forest frogs of the present study. We hypothesize this difference is related mainly to habitat but cannot exclude a phylogenetic component as the three swamp-forest species are a closely related monophyletic group. The paucity of alkaloids with unbranched-carbon skeletons (ant-derived) and the commonness of alkaloids with branched-carbon skeletons (mite-derived) indicate that oribatid mites are a major source of alkaloids in these species of mantellids. Furthermore, most of the alkaloids have an oxygen atom in their formulae. Differences in alkaloids were observed among species, populations of the same species, and habitats. In M. aurantiaca, small geographic distances among populations were associated with differences in alkaloid profiles, with a remote third site illustrating even greater differences. The present study and an earlier study of three other mantellid species suggest that oribatid mites, and not ants, are the major source of alkaloids in the species of mantellids examined thus far.

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.

Seabird aggregative patterns: a new tool for offshore wind energy risk assessment.

The emerging development of offshore wind energy has raised public concern over its impact on seabird communities. There is a need for an adequate methodology to determine its potential impacts on seabirds. Environmental Impact Assessments (EIAs) are mostly relying on a succession of plain density maps without integrated interpretation of seabird spatio-temporal variability. Using Taylor's power law coupled with mixed effect models, the spatio-temporal variability of species' distributions can be synthesized in a measure of the aggregation levels of individuals over time and space. Applying the method to a seabird aerial survey in the Ebro Delta, NW Mediterranean Sea, we were able to make an explicit distinction between transitional and feeding areas to define and map the potential impacts of an offshore wind farm project. We use the Ebro Delta study case to discuss the advantages of potential impacts maps over density maps, as well as to illustrate how these potential impact maps can be applied to inform on concern levels, optimal EIA design and monitoring in the assessment of local offshore wind energy projects.

Vertebrate time-tree elucidates the biogeographic pattern of a major biotic change around the K-T boundary in Madagascar.

The geographic and temporal origins of Madagascar's biota have long been in the center of debate. We reconstructed a time-tree including nearly all native nonflying and nonmarine vertebrate clades present on the island, from DNA sequences of two single-copy protein-coding nuclear genes (BDNF and RAG1) and a set of congruent time constraints. Reconstructions calculated with autocorrelated or independent substitution rates over clades agreed in placing the origins of the 31 included clades in Cretaceous to Cenozoic times. The two clades with sister groups in South America were the oldest, followed by those of a putative Asian ancestry that were significantly older than the prevalent clades of African ancestry. No colonizations from Asia occurred after the Eocene, suggesting that dispersal and vicariance of Asian/Indian groups were favored over a comparatively short period during, and shortly after, the separation of India and Madagascar. Species richness of clades correlates with their age but those clades that have a large proportion of species diversity in rainforests are significantly more species-rich. This finding suggests an underlying pattern of continuous speciation through time in Madagascar's vertebrates, with accelerated episodes of adaptive diversification in those clades that succeeded radiating into the rainforests.

Discordant patterns of nuclear and mitochondrial introgression in Iberian populations of the European common frog (Rana temporaria).

Amphibians often show complex histories of intraspecific and interspecific genetic introgression, which might differ in mitochondrial and nuclear genes. In our study of the genetic differentiation of the European common frog, Rana temporaria (159 specimens from 23 populations were analyzed for 24 presumptive allozyme loci; 82 specimens were sequenced for a 540-bp fragment of the mitochondrial 16S rRNA gene), multilocus correspondence analysis (CA) and Bayesian assignment tests of the nuclear data were concordant in identifying 2 population groups corresponding to 1) the Pyrenees in the east and 2) the Galicia and Asturias regions in the west, the latter corresponding to the subspecies R. temporaria parvipalmata. Geographically intermediate populations were genetically intermediate in the allozyme CA and, less clearly in the Bayesian assignment, with mitochondrial haplotypes exclusively belonging to the parvipalmata group. This indicates different degrees of introgression in the mitochondrial and nuclear genomes. Although Pyrenean high-altitude populations are morphologically distinct from low-altitude populations, these 2 groups were not separate clusters in any analysis. This suggests that the morphological differences may be due to fast adaptation to elevational gradients, likely under maintenance of gene flow, and that the underlying genetic changes are not detectable by the analyzed markers. We argue that a parsimonious explanation for the observed pattern along the east-west axis in northern Spain may be competition between invading and resident populations, with no need to invoke selection. However, in order to conclusively rule out selective processes, additional and finer scale data are required to test for asymmetric mating preference/behaviour, sex-biased gene flow, or sex-biased survival of potential hybrids.

Molecular phylogeny and biogeography of Malagasy frogs of the genus Gephyromantis.

The genus Gephyromantis is a clade within the Malagasy-Comoroan family Mantellidae composed of rainforest frogs that live and breed to varying degrees independently from water. Based on DNA sequences of five mitochondrial and five nuclear genes we inferred the phylogeny of these frogs with full taxon coverage at the species level. Our preferred consensus tree from a partitioned Bayesian analysis of 5843 base pairs of 51 nominal and candidate species supports various major clades within the genus although the basal relationships among these remain unresolved. The data provide strong evidence for the monophyly of the subgenera Gephyromantis (after exclusion of Gephyromantis klemmeri), Laurentomantis, Vatomantis, and Phylacomantis. Species assigned to the subgenus Duboimantis belong to two strongly supported clades of uncertain relationships. G. klemmeri, previously in the subgenus Gephyromantis, was placed with high support sister to the Laurentomantis clade, and the Laurentomantis+G. klemmeri clade was sister to Vatomantis. A reconstruction of ancestral distribution areas indicates a diversification of several subgenera in the northern biogeographic regions of Madagascar and the dispersal out of northern Madagascar for several clades.

Speciation in little: the role of range and body size in the diversification of Malagasy mantellid frogs.

BACKGROUND: The rate and mode of lineage diversification might be shaped by clade-specific traits. In Madagascar, many groups of organisms are characterized by tiny distribution ranges and small body sizes, and this high degree of microendemism and miniaturization parallels a high species diversity in some of these groups. We here investigate the geographic patterns characterizing the radiation of the frog family Mantellidae that is virtually endemic to Madagascar. We integrate a newly reconstructed near-complete species-level timetree of the Mantellidae with georeferenced distribution records and maximum male body size data to infer the influence of these life-history traits on each other and on mantellid diversification. RESULTS: We reconstructed a molecular phylogeny based on nuclear and mitochondrial DNA for 257 species and candidate species of the mantellid frog radiation. Based on this phylogeny we identified 53 well-supported pairs of sister species that we used for phylogenetic comparative analyses, along with whole tree-based phylogenetic comparative methods. Sister species within the Mantellidae diverged at 0.2-14.4 million years ago and more recently diverged sister species had geographical range centroids more proximate to each other, independently of their current sympatric or allopatric occurrence. The largest number of sister species pairs had non-overlapping ranges, but several examples of young microendemic sister species occurring in full sympatry suggest the possibility of non-allopatric speciation. Range sizes of species included in the sister species comparisons increased with evolutionary age, as did range size differences between sister species, which rejects peripatric speciation. For the majority of mantellid sister species and the whole mantellid radiation, range and body sizes were associated with each other and small body sizes were linked to higher mitochondrial nucleotide substitution rates and higher clade diversity. In contrast, small range sizes were unexpectedly associated with a slow-down of mitochondrial substitution rates. CONCLUSIONS: Based on these results we define a testable hypothesis under which small body sizes result in limited dispersal capabilities and low physiological tolerances, causing smaller and more strongly fragmented ranges. This can be thought to facilitate reproductive isolation and thus favor speciation. Contrary to the expectation of the faster speciation of such microendemic phenotype species, we only found small body sizes of mantellid frogs to be linked to higher diversification and substitution rates, but not small range sizes. A joint analysis of various species-rich regional anuran radiations might provide enough species with all combinations of range and body sizes for a more conclusive test of this hypothesis.

A multigenic perspective on phylogenetic relationships in the largest family of salamanders, the Plethodontidae.

Despite several recent studies, the phylogeny of plethodontid salamanders is not yet fully resolved and the phylogenetic positions of several key genera, especially Aneides, Hemidactylium, Hydromantes and Karsenia, are contentious. Here we present a combined dataset of complete mitochondrial genomes and three nuclear loci for 20 species (16 genera) of plethodontids, representing all major clades in the family. The combined dataset without mitochondrial third codon positions provides a fully resolved, statistically well-supported tree. In this topology two major clades are recovered. A northern clade includes Aneides, Desmognathus, Ensatina, Hydromantes, Karsenia, Phaeognathus and Plethodon, with Plethodon being the sister taxon to the rest of the clade. Hydromantes and Karsenia are sister taxa, and Aneides is recovered as the sister taxon to Ensatina. Desmognathus+Phaeognathus form the sister taxon to Aneides+Ensatina. An eastern/southern clade comprises two subclades. One subclade, the spelerpines (Eurycea, Gyrinophilus, Pseudotriton, Stereochilus, Urspelerpes) is the sister taxon to a subclade comprising Hemidactylium, Batrachoseps and the tropical plethodontids (represented by Bolitoglossa, Nototriton and Thorius). In this topology Hemidactylium is well-supported as the sister taxon to Batrachoseps. Only when mitochondrial third codon positions are included using maximum likelihood analysis is Hemidactylium recovered as the sister taxon to Batrachoseps+tropical genera. Hypothesis testing of alternative topologies supports these conclusions. On the basis of these results we propose a conservative taxonomy for Plethodontidae.

Deep genealogical lineages in the widely distributed African helmeted terrapin: evidence from mitochondrial and nuclear DNA (Testudines: Pelomedusidae: Pelomedusa subrufa).

We investigated the phylogeographic differentiation of the widely distributed African helmeted terrapin Pelomedusa subrufa based on 1503 base pairs of mitochondrial DNA (partial cyt b and ND4 genes with adjacent tRNAs) and 1937 bp of nuclear DNA (partial Rag1, Rag2, R35 genes). Congruent among different analyses, nine strongly divergent mitochondrial clades were found, representing three major geographical groupings: (1) A northern group which includes clades I from Cameroon, II from Ghana and Ivory Coast, III from Benin, Burkina Faso and Niger, IV from the Central African Republic, and V from Kenya, (2) a northeastern group consisting of clades VI from Somalia, and VII from Saudi Arabia and Yemen, and (3) a southern group comprising clade VIII from Botswana, the Democratic Republic of Congo, Madagascar and Malawi, and clade IX from South Africa. Malagasy and continental African populations were not clearly differentiated, indicating very recent arrival or introduction of Pelomedusa in Madagascar. The southern group was in some phylogenetic analyses sister to Pelusios, rendering Pelomedusa paraphyletic with respect to that genus. However, using partitioned Bayesian analyses and sequence data of the three nuclear genes, Pelomedusa was monophyletic, suggesting that its mitochondrial paraphyly is due to either ancient introgressive hybridization or phylogenetic noise. Otherwise, nuclear sequence data recovered a lower level of divergence, but corroborated the general differentiation pattern of Pelomedusa as revealed by mtDNA. This, and the depth of the divergences between clades, indicates ancient differentiation. The divergences observed fall within, and in part exceed considerably, the differentiation typically occurring among chelonian species. To test whether Pelomedusa is best considered a single species composed of deep genealogical lineages, or a complex of up to nine distinct species, we suggest a future taxonomic revision that should (1) extend the geographical sampling of molecular data, specifically focusing on contact zones and the possible sympatric occurrence of lineages without admixture, and (2) evaluate the morphology of the various genealogical lineages using the type specimens or topotypical material of the numerous junior synonyms of P. subrufa.