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.

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.

Data concatenation, Bayesian concordance and coalescent-based analyses of the species tree for the rapid radiation of Triturus newts.

The phylogenetic relationships for rapid species radiations are difficult to disentangle. Here we study one such case, namely the genus Triturus, which is composed of the marbled and crested newts. We analyze data for 38 genetic markers, positioned in 3-prime untranslated regions of protein-coding genes, obtained with 454 sequencing. Our dataset includes twenty Triturus newts and represents all nine species. Bayesian analysis of population structure allocates all individuals to their respective species. The branching patterns obtained by data concatenation, Bayesian concordance analysis and coalescent-based estimations of the species tree differ from one another. The data concatenation based species tree shows high branch support but branching order is considerably affected by allele choice in the case of heterozygotes in the concatenation process. Bayesian concordance analysis expresses the conflict between individual gene trees for part of the Triturus species tree as low concordance factors. The coalescent-based species tree is relatively similar to a previously published species tree based upon morphology and full mtDNA and any conflicting internal branches are not highly supported. Our findings reflect high gene tree discordance due to incomplete lineage sorting (possibly aggravated by hybridization) in combination with low information content of the markers employed (as can be expected for relatively recent species radiations). This case study highlights the complexity of resolving rapid radiations and we acknowledge that to convincingly resolve the Triturus species tree even more genes will have to be consulted.

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.

Exploring the effect of asymmetric mitochondrial DNA introgression on estimating niche divergence in morphologically cryptic species.

If potential morphologically cryptic species, identified based on differentiated mitochondrial DNA, express ecological divergence, this increases support for their treatment as distinct species. However, mitochondrial DNA introgression hampers the correct estimation of ecological divergence. We test the hypothesis that estimated niche divergence differs when considering nuclear DNA composition or mitochondrial DNA type as representing the true species range. We use empirical data of two crested newt species (Amphibia: Triturus) which possess introgressed mitochondrial DNA from a third species in part of their ranges. We analyze the data in environmental space by determining Fisher distances in a principal component analysis and in geographical space by determining geographical overlap of species distribution models. We find that under mtDNA guidance in one of the two study cases niche divergence is overestimated, whereas in the other it is underestimated. In the light of our results we discuss the role of estimated niche divergence in species delineation.

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.

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.

Postglacial species displacement in Triturus newts deduced from asymmetrically introgressed mitochondrial DNA and ecological niche models.

BACKGROUND: If the geographical displacement of one species by another is accompanied by hybridization, mitochondrial DNA can introgress asymmetrically, from the outcompeted species into the invading species, over a large area. We explore this phenomenon using the two parapatric crested newt species, Triturus macedonicus and T. karelinii, distributed on the Balkan Peninsula in south-eastern Europe, as a model. RESULTS: We first delimit a ca. 54,000 km(2) area in which T. macedonicus contains T. karelinii mitochondrial DNA. This introgression zone bisects the range of T. karelinii, cutting off a T. karelinii enclave. The high similarity of introgressed mitochondrial DNA haplotypes with those found in T. karelinii suggests a recent transfer across the species boundary. We then use ecological niche modeling to explore habitat suitability of the location of the present day introgression zone under current, mid-Holocene and Last Glacial Maximum conditions. This area was inhospitable during the Last Glacial Maximum for both species, but would have been habitable at the mid-Holocene. Since the mid-Holocene, habitat suitability generally increased for T. macedonicus, whereas it decreased for T. karelinii. CONCLUSION: The presence of a T. karelinii enclave suggests that T. karelinii was the first to colonize the area where the present day introgression zone is positioned after the Last Glacial Maximum. Subsequently, we propose T. karelinii was outcompeted by T. macedonicus, which captured T. karelinii mitochondrial DNA via introgressive hybridization in the process. Ecological niche modeling suggests that this replacement was likely facilitated by a shift in climate since the mid-Holocene. We suggest that the northwestern part of the current introgression zone was probably never inhabited by T. karelinii itself, and that T. karelinii mitochondrial DNA spread there through T. macedonicus exclusively. Considering the spatial distribution of the introgressed mitochondrial DNA and the signal derived from ecological niche modeling, we do not favor the hypothesis that foreign mitochondrial DNA was pulled into the T. macedonicus range by natural selection.

Is dispersal guided by the environment? A comparison of interspecific gene flow estimates among differentiated regions of a newt hybrid zone.

How frequently genes pass through a hybrid zone may be influenced by the environment. Accordingly, in long hybrid zones that span more than one environmental setting, different patterns may emerge. The varied conditions allow testing of hypotheses on dispersal as a function of the environment. We reconstruct the amount and direction of gene flow across a heterogeneous hybrid zone of two species of marbled newts (Triturus marmoratus and Triturus pygmaeus), in four widely separated areas of the Iberian Peninsula from one mitochondrial and three nuclear genes. The main variables associated with the position of the contact zone are precipitation, rivers, altitude and relief. In some sections of the contact zone, however, its position is not correlated with any environmental factor and is instead determined by the shortest geographical distance between fixed positions at either side (mountains in the East and river in the West). In areas where the position of the zone is stable, gene flow was bidirectional. External data show that T. pygmaeus has superseded T. marmoratus over a large area and here gene flow was unidirectional. The prediction that a major river would reduce gene flow was not confirmed.