Two transects reveal remarkable variation in gene flow on opposite ends of a European toad hybrid zone.

Speciation entails a reduction in gene flow between lineages. The rates at which genomic regions become isolated varies across space and time. Barrier markers are linked to putative genes involved in (processes of) reproductive isolation, and, when observed over two transects, indicate species-wide processes. In contrast, transect-specific putative barrier markers suggest local processes. We studied two widely separated transects along the 900 km hybrid zone between Bufo bufo and B. spinosus, in northern and southern France, for ~1200 RADseq markers. We used genomic and geographic cline analyses to identify barrier markers based on their restricted introgression, and found that some markers are transect-specific, while others are shared between transects. Twenty-six barrier markers were shared across both transects, of which some are clustered in the same chromosomal region, suggesting that their associated genes are involved in reduced gene flow across the entire hybrid zone. Transect-specific barrier markers were twice as numerous in the southern than in the northern transect, suggesting that the overall barrier effect is weaker in northern France. We hypothesize that this is consistent with a longer period of secondary contact in southern France. The smaller number of introgressed genes in the northern transect shows considerably more gene flow towards the southern (B. spinosus) than the northern species (B. bufo). We hypothesize that hybrid zone movement in northern France and hybrid zone stability in southern France explain this pattern. The Bufo hybrid zone provides an excellent opportunity to separate a general barrier effect from localized gene flow-reducing conditions.

Phylogenomics of the adaptive radiation of Triturus newts supports gradual ecological niche expansion towards an incrementally aquatic lifestyle.

Newts of the genus Triturus (marbled and crested newts) exhibit substantial variation in the number of trunk vertebrae (NTV) and a higher NTV corresponds to a longer annual aquatic period. Because the Triturus phylogeny has thwarted resolution to date, the evolutionary history of NTV, annual aquatic period, and their potential coevolution has remained unclear. To resolve the phylogeny of Triturus, we generated a c. 6000 transcriptome-derived marker data set using a custom target enrichment probe set, and conducted phylogenetic analyses using: (1) data concatenation with RAxML, (2) gene-tree summary with ASTRAL, and (3) species-tree estimation with SNAPP. All analyses produce the same, highly supported topology, despite cladogenesis having occurred over a short timeframe, resulting in short internal branch lengths. Our new phylogenetic hypothesis is consistent with the minimal number of inferred changes in NTV count necessary to explain the diversity in NTV observed today. Although a causal relationship between NTV, body form, and aquatic ecology has yet to be experimentally established, our phylogeny indicates that these features have evolved together, and suggest that they may underlie the adaptive radiation that characterizes Triturus.

A signature of dynamic biogeography: enclaves indicate past species replacement.

Understanding how species have replaced each other in the past is important to predicting future species turnover. While past species replacement is difficult to detect after the fact, the process may be inferred from present-day distribution patterns. Species with abutting ranges sometimes show a characteristic distribution pattern, where a section of one species range is enveloped by that of the other. Such an enclave could indicate past species replacement: when a species is partly supplanted by a competitor, but a population endures locally while the invading species moves around and past it, an enclave forms. If the two species hybridize and backcross, the receding species is predicted to leave genetic traces within the expanding one under a scenario of species replacement. By screening dozens of genes in hybridizing crested newts, we uncover genetic remnants of the ancestral species, now inhabiting an enclave, in the range of the surrounding invading species. This independent genetic evidence supports the past distribution dynamics we predicted from the enclave. We suggest that enclaves provide a valuable tool in understanding historical species replacement, which is important because a major conservation concern arising from anthropogenic climate change is increased species replacement in the future.

Description of a new species of crested newt, previously subsumed in Triturus ivanbureschi (Amphibia: Caudata: Salamandridae).

Multilocus molecular data play a pivotal role in diagnosing cryptic species (i.e. genetically distinct but morphologically similar species). A multilocus phylogeographic survey has provided compelling evidence that Triturus ivanbureschi sensu lato comprises two distinct gene pools with restricted gene flow. We conclude that this taxon had better be treated as two distinct (albeit morphologically cryptic) species. The name T. ivanbureschi should be restricted to the western species, which is distributed in western Asiatic Turkey plus the south-eastern Balkan Peninsula. No name is as yet available for the eastern species, which is distributed in northern Asiatic Turkey. We propose the name T. anatolicus sp. nov. for the eastern species and provide a formal species description.

Kicking Triturus arntzeni when it's down: large-scale nuclear genetic data confirm that newts from the type locality are genetically admixed.

We collected nuclear DNA data (52 markers) with next-generation sequencing for nine Triturus newt specimens, including the holotype and two of the paratypes of T. arntzeni, from the type locality at Vrtovac in eastern Serbia. We compare these data to a reference set composed of the four crested newt species distributed in eastern Serbia namely T. cristatus, T. dobrogicus, T. ivanbureschi and T. macedonicus to determine to which of these species the newts from the type locality of T. arntzeni should be attributed. The majority of alleles in individuals from Vrtovac is derived from T. macedonicus, but a considerable number of T. ivanbureschi alleles is also present; alleles typical for T. cristatus and T. dobrogicus are found at low frequency. Accordingly, we interpret Vrtovac as a T. macedonicus - T. ivanbureschi hybrid population, albeit not composed of F1 hybrids but of genetically admixed individuals derived through multiple generations of backcrossing. The data support the notion that the name T. arntzeni should not be applied to a species newly distinguished in T. karelinii sensu lato (to which the name T. ivanbureschi has been given). We conclude that because of the hybrid nature of the individuals from Vrtovac, the name T. arntzeni should be placed not only in the synonymy of T. macedonicus but also in the synonymy of T. ivanbureschi. In this study we demonstrate that next-generation sequencing can provide high quality data for type material with degraded DNA and therefore can play an important role in taxonomy.

Parallel tagged amplicon sequencing of transcriptome-based genetic markers for Triturus newts with the Ion Torrent next-generation sequencing platform.

Next-generation sequencing is a fast and cost-effective way to obtain sequence data for nonmodel organisms for many markers and for many individuals. We describe a protocol through which we obtain orthologous markers for the crested newts (Amphibia: Salamandridae: Triturus), suitable for analysis of interspecific hybridization. We use transcriptome data of a single Triturus species and design 96 primer pairs that amplify c. 180 bp fragments positioned in 3-prime untranslated regions. Next, these markers are tested with uniplex PCR for a set of species spanning the taxonomical width of the genus Triturus. The 52 markers that consistently show a single band of expected length at gel electrophoreses for all tested crested newt species are then amplified in five multiplex PCRs (with a plexity of ten or eleven) for 132 individual newts: a set of 84 representing the seven (candidate) species and a set of 48 from a presumed hybrid population. After pooling multiplexes per individual, unique tags are ligated to link amplicons to individuals. Subsequently, individuals are pooled equimolar and sequenced on the Ion Torrent next-generation sequencing platform. A bioinformatics pipeline identifies the alleles and recodes these to a genotypic format. Next, we test the utility of our markers. baps allocates the 84 crested newt individuals representing (candidate) species to their expected (candidate) species, confirming the markers are suitable for species delineation. newhybrids, a hybrid index and hiest confirm the 48 individuals from the presumed hybrid population to be genetically admixed, illustrating the potential of the markers to identify interspecific hybridization. We expect the set of markers we designed to provide a high resolving power for analysis of hybridization in Triturus.

No evidence for nuclear introgression despite complete mtDNA replacement in the Carpathian newt (Lissotriton montandoni).

Patterns of interspecific introgression may vary geographically, and the distribution of introgressed variants can yield insight into the historical dynamics of genetic interactions between hybridizing species. Urodele amphibians, often characterized by limited mobility, deep intraspecific genetic structuring and vulnerability to climatic changes, constitute suitable models for such historical inferences. Here, we combine an extensive survey of the mitochondrial (mtDNA) and nuclear (15 microsatellites) genomes in the Carpathian newt, Lissotriton montandoni (Lm) with species distribution modelling (SDM). Populations of the smooth newt, L. vulgaris (Lv) from the areas surrounding the Lm range were also sampled to test whether gene flow between these hybridizing species extends beyond the area of strict syntopy. The extent of introgression differs dramatically between the mitochondrial genome and the nuclear genome. While multiple, spatially and temporally distinct introgression events from Lv resulted in complete mtDNA replacement in Lm, there was little evidence of recent interspecific nuclear gene flow in the assayed markers. Microsatellite differentiation within Lm defines three units, probably derived from separate glacial refugia, located in the northern, eastern and southern part of the Carpathians. In situ survival and range fragmentation of Lm are supported by SDM, corroborating the role of the Carpathians as a major refugial area. Our results, in combination with previous reports of extensive introgression of the major histocompatibility complex (MHC) genes, emphasize the complexity of historical gene exchange between Lm and Lv.

A multimarker phylogeography of crested newts (Triturus cristatus superspecies) reveals cryptic species.

The crested newt Triturus cristatus superspecies is composed of five recognized species. One of these, T. karelinii sensu lato, comprises three geographically structured mitochondrial DNA lineages: 'eastern', 'central' and 'western T. karelinii'. Genetic divergence among these lineages is comparable to that of recognized crested newt species, but morphologically they are indistinguishable. Here, we conduct a multimarker phylogeographical survey to explore the evolutionary independence of these mitochondrial DNA lineages and we include representatives of the other species to guide our interpretation of the results. All markers show distinct patterns when analyzed singly (as a phylogeny or haplotype network) and none of them sort haplotypes fully in line with species or mitochondrial DNA lineage. A multilocus approach (BAPS and *BEAST) on the other hand shows that not only the recognized species, but also the three mitochondrial DNA lineages represent discrete nuclear DNA gene pools. A mismatch is found in the northwest of Asiatic Turkey, where several populations identified as 'central T. karelinii' based on nuclear DNA possesses 'western T. karelinii' mitochondrial DNA. We invoke asymmetric mitochondrial DNA introgression to explain this pattern and support this with a historical biogeographical scenario. The three spatial groups in T. karelinii sensu lato should be regarded as distinct species.

A revised taxonomy of crested newts in the Triturus karelinii group (Amphibia: Caudata: Salamandridae), with the description of a new species.

We present a taxonomic revision of the crested newt Triturus karelinii sensu lato. Based on the presence of discrete nuclear DNA gene pools, deep genetic divergence of mitochondrial and nuclear DNA, and no indication of gene flow, we interpret this taxon as comprising two species: one covering the southern Caspian Sea shore, the Caucasus and the Crimea, i.e. the eastern part of the total range and another covering northern Asiatic Turkey and western Asiatic Turkey plus the southeastern Balkan Peninsula, i.e. the central and western part of the total range. We acknowledge that the central/western species should likely be further subdivided into a central and a western taxon, but we prefer to await a more detailed genetic analysis of the putative contact zone, positioned in northwestern Asiatic Turkey. The name T. karelinii (Strauch, 1870) applies to the eastern species as the type locality is positioned along the coast of the Gulf of Gorgan, Iran. The name T. arntzeni has been applied to the central/western species with Vrtovac, Serbia as the type locality. We show that not T. karelinii sensu lato but T. macedonicus occurs at Vrtovac. Hence, the name T. arntzeni Litvinchuk, Borkin, Dzukic and Kalezic, 1999 (in Litvinchuk et al., 1999) is a junior synonym of T. macedonicus (Karaman, 1922) and should not be used for the central/western species. We propose the name T. ivanbureschi sp. nov. for the central/western species and provide a formal species description.

Cryptic crested newt diversity at the Eurasian transition: the mitochondrial DNA phylogeography of Near Eastern Triturus newts.

Crested newts of the Triturus karelinii group occur in a phylogeographically understudied region: the Near East. Controversy surrounds the systematic position of these newts within the complete crested newt assemblage (the Triturus cristatus superspecies). We explore the situation using mitochondrial sequence data (ND2 and ND4, approximately 1.7kb) and employing different methods of phylogenetic inference (Bayesian inference and Maximum Likelihood using mixed models) and molecular dating (r8s and BEAST). The T. karelinii group is monophyletic and constitutes one of four main lineages in the T. cristatus superspecies. The separation of the T. karelinii group from the remaining crested newts around 9Ma is related to the formation of the Mid-Aegean Trench, which separated the Balkan and Anatolian landmasses. The T. karelinii group comprises three geographically structured clades (eastern, central and western). The genetic divergence shown by these clades is comparable to that among recognized crested newt species. We suggest the uplift of the Armenian Plateau to be responsible for the separation of the eastern clade around 7Ma, and the re-establishment of a marine connection between the Black Sea and the Mediterranean at the end of the Messinian Salinity Crisis to have caused the split between the central and western clade around 5.5Ma. Genetic structuring within the three clades dates to the Quaternary Ice Age (<2.59Ma) and is associated with alternating periods of isolation and reconnection caused by periodic changes in sea level and surface runoff.

Multiple nuclear and mitochondrial genes resolve the branching order of a rapid radiation of crested newts (Triturus arntzenii) (corrected).

Newts of the genus Triturus are parapatrically distributed across Europe. Within this group, the crested newts (Triturus cristatus superspecies) radiated in a short temporal interval. Given the relatively short timespan in between branching events and to address the gene tree-species tree problem, we sequenced two mitochondrial and five nuclear genes from populations representing the distribution range of all the five crested newt species. We built gene trees using non-hierarchical Bayesian phylogenetics and phylogenetic networks, and a species tree with a recently developed method, which uses a hierarchical Bayesian approach. While the single gene trees did not provide resolution, the hierarchical Bayesian method yielded an almost fully resolved species tree, even though branching events followed one another closely. Results show a previously undetected basal dichotomy between T. karelinii and the other four species and a deep differentiation of T. karelinii in two lineages, here raised to full species status.