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.

Coevolution between MHC Class I and Antigen-Processing Genes in Salamanders.

Proteins encoded by antigen-processing genes (APGs) provide major histocompatibility complex (MHC) class I (MHC-I) with antigenic peptides. In mammals, polymorphic multigenic MHC-I family is served by monomorphic APGs, whereas in certain nonmammalian species both MHC-I and APGs are polymorphic and coevolve within stable haplotypes. Coevolution was suggested as an ancestral gnathostome feature, presumably enabling only a single highly expressed classical MHC-I gene. In this view coevolution, while optimizing some aspects of adaptive immunity, would also limit its flexibility by preventing the expansion of classical MHC-I into a multigene family. However, some nonmammalian taxa, such as salamanders, have multiple highly expressed MHC-I genes, suggesting either that coevolution is relaxed or that it does not prevent the establishment of multigene MHC-I. To distinguish between these two alternatives, we use salamanders (30 species from 16 genera representing six families) to test, within a comparative framework, a major prediction of the coevolution hypothesis: the positive correlation between MHC-I and APG diversity. We found that MHC-I diversity explained both within-individual and species-wide diversity of two APGs, TAP1 and TAP2, supporting their coevolution with MHC-I, whereas no consistent effect was detected for the other three APGs (PSMB8, PSMB9, and TAPBP). Our results imply that although coevolution occurs in salamanders, it does not preclude the expansion of the MHC-I gene family. Contrary to the previous suggestions, nonmammalian vertebrates thus may be able to accommodate diverse selection pressures with flexibility granted by rapid expansion or contraction of the MHC-I family, while retaining the benefits of coevolution between MHC-I and TAPs.

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.

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.

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.

Hybridization and speciation.

Hybridization has many and varied impacts on the process of speciation. Hybridization may slow or reverse differentiation by allowing gene flow and recombination. It may accelerate speciation via adaptive introgression or cause near-instantaneous speciation by allopolyploidization. It may have multiple effects at different stages and in different spatial contexts within a single speciation event. We offer a perspective on the context and evolutionary significance of hybridization during speciation, highlighting issues of current interest and debate. In secondary contact zones, it is uncertain if barriers to gene flow will be strengthened or broken down due to recombination and gene flow. Theory and empirical evidence suggest the latter is more likely, except within and around strongly selected genomic regions. Hybridization may contribute to speciation through the formation of new hybrid taxa, whereas introgression of a few loci may promote adaptive divergence and so facilitate speciation. Gene regulatory networks, epigenetic effects and the evolution of selfish genetic material in the genome suggest that the Dobzhansky-Muller model of hybrid incompatibilities requires a broader interpretation. Finally, although the incidence of reinforcement remains uncertain, this and other interactions in areas of sympatry may have knock-on effects on speciation both within and outside regions of hybridization.

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.

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.

Multilocus species tree analyses resolve the radiation of the widespread Bufo bufo species group (Anura, Bufonidae).

New analytical methods are improving our ability to reconstruct robust species trees from multilocus datasets, despite difficulties in phylogenetic reconstruction associated with recent, rapid divergence, incomplete lineage sorting and/or introgression. In this study, we applied these methods to resolve the radiation of toads in the Bufo bufo (Anura, Bufonidae) species group, ranging from the Iberian Peninsula and North Africa to Siberia, based on sequences from two mitochondrial and four nuclear DNA regions (3490 base pairs). We obtained a fully-resolved topology, with the recently described Bufo eichwaldi from the Talysh Mountains in south Azerbaijan and Iran as the sister taxon to a clade including: (1) north African, Iberian, and most French populations, referred herein to Bufo spinosus based on the implied inclusion of populations from its type locality and (2) a second clade, sister to B. spinosus, including two sister subclades: one with all samples of Bufo verrucosissimus from the Caucasus and another one with samples of B. bufo from northern France to Russia, including the Apennine and Balkan peninsulas and most of Anatolia. Coalescent-based estimations of time to most recent common ancestors for each species and selected subclades allowed historical reconstruction of the diversification of the species group in the context of Mediterranean paleogeography and indicated a long evolutionary history in this region. Finally, we used our data to delimit the ranges of the four species, particularly the more widespread and historically confused B. spinosus and B. bufo, and identify potential contact zones, some of which show striking parallels with other co-distributed species.

Crest evolution in newts: implications for reconstruction methods, sexual selection, phenotypic plasticity and the origin of novelties.

The dorsal crest of newts (Salamandridae) is a novel, phenotypically plastic, sexually selected trait that may evolve in association with complex courtship behaviours. We estimated a near-comprehensive, time-calibrated phylogeny for salamandrids and analysed the evolution of their crests and display behaviour. Different models give conflicting reconstructions for crest evolution, showing that likelihood can estimate incorrect ancestral states with strong statistical support. The best-fitting model suggests that crests evolved once and were lost repeatedly, supporting the hypothesis that sexually selected traits may be frequently lost. We demonstrate the correlated evolution of crests and courtship behaviour and show that species with larger numbers of crest-related traits have larger repertoires of behaviours. We also show that phenotypically plastic morphological traits can be maintained over long macroevolutionary timescales (∼25-48 Myr). Finally, we use salamandrids to address how novel structures may arise, and support a model involving the expansion and subdivision of pre-existing structures.

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.

Long-term survival of a urodele amphibian despite depleted major histocompatibility complex variation.

Depletion of polymorphism at major histocompatibility complex (MHC) genes has been hypothesized to limit the ability of populations to respond to emerging pathogens, thus putting their survival at risk. As pathogens contribute substantially to the global amphibian decline, assessing patterns of MHC variation is important in devising conservation strategies. Here, we directly compare levels of MHC class II and neutral variation between multiple populations of the great crested newt (Triturus cristatus) from refugial (REF: Romania) and postglacial expansion (PGE: Germany, Poland and UK) areas. REF populations harboured high levels of adaptive variation (24 expressed alleles), exhibiting clear signatures of historical positive selection, which points to the overall importance of MHC class II variation in this species. On the other hand, PGE populations were extremely depauperate (two alleles) but nevertheless have survived for c. 10,000 years, since the postglacial expansion. Variation in putative MHC class II pseudogenes, microsatellites and allozymes also showed a significant southern richness-northern purity pattern. The populations in the postglacial expansion area thus provide the clearest example to date of the long-term survival of populations in which MHC variation, historically under positive selection, has been depleted.

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.

Heterogeneous responses of temperate-zone amphibian populations to climate change complicates conservation planning.

The pervasive and unabated nature of global amphibian declines suggests common demographic responses to a given driver, and quantification of major drivers and responses could inform broad-scale conservation actions. We explored the influence of climate on demographic parameters (i.e., changes in the probabilities of survival and recruitment) using 31 datasets from temperate zone amphibian populations (North America and Europe) with more than a decade of observations each. There was evidence for an influence of climate on population demographic rates, but the direction and magnitude of responses to climate drivers was highly variable among taxa and among populations within taxa. These results reveal that climate drivers interact with variation in life-history traits and population-specific attributes resulting in a diversity of responses. This heterogeneity complicates the identification of conservation 'rules of thumb' for these taxa, and supports the notion of local focus as the most effective approach to overcome global-scale conservation challenges.

From descriptive to predictive distribution models: a working example with Iberian amphibians and reptiles.

BACKGROUND: Aim of the study was to identify the conditions under which spatial-environmental models can be used for the improved understanding of species distributions, under the explicit criterion of model predictive performance. I constructed distribution models for 17 amphibian and 21 reptile species in Portugal from atlas data and 13 selected ecological variables with stepwise logistic regression and a geographic information system. Models constructed for Portugal were extrapolated over Spain and tested against range maps and atlas data. RESULTS: Descriptive model precision ranged from 'fair' to 'very good' for 12 species showing a range border inside Portugal ('edge species', kappa (k) 0.35-0.89, average 0.57) and was at best 'moderate' for 26 species with a countrywide Portuguese distribution ('non-edge species', k = 0.03-0.54, average 0.29). The accuracy of the prediction for Spain was significantly related to the precision of the descriptive model for the group of edge species and not for the countrywide species. In the latter group data were consistently better captured with the single variable search-effort than by the panel of environmental data. CONCLUSION: Atlas data in presence-absence format are often inadequate to model the distribution of species if the considered area does not include part of the range border. Conversely, distribution models for edge-species, especially those displaying high precision, may help in the correct identification of parameters underlying the species range and assist with the informed choice of conservation measures.

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.

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 combination of techniques proves useful in the development of nuclear markers in the newt genus Triturus.

To increase the number of markers available for study of phylogeny and phylogeography in the newt genus Triturus, we developed and tested 59 primer pairs using three different techniques. Primers were obtained from published sources, by designing exon-primed intron-crossing primers and from randomly cloned anonymous nuclear DNA fragments. Successful polymerase chain reaction products were cloned and sequenced. Five fragments were successfully amplified and sequenced for six species of Triturus: intron 7 of the ß-fibrinogen gene (ßfibint7), third intron of the calreticulin gene (CalintC), the 11th intron of the α-subunit of the platelet derived growth factor receptor (PDGFRα) and two anonymous markers (Cri1 and Cri4). The average percentage species divergence across all the markers is low (c. 3%), compared to what has been found in mitochondrial DNA (25-30%).

A Bayesian approach on molecules and behavior: reconsidering phylogenetic and evolutionary patterns of the Salamandridae with emphasis on Triturus newts.

The monophyly of European newts of the genus Triturus within the family Salamandridae has for decades rested on presumably homologous behavioral and morphological characters. Molecular data challenge this hypothesis, but the phylogenetic position of Triturus within the Salamandridae has not yet been convincingly resolved. We addressed this issue and the temporal divergence of Triturus within the Salamandridae with novel Bayesian approaches applied to DNA sequence data from three mitochondrial genes (12S, 16S and cytb). We included 38 salamandrid species comprising all 13 recognized species of Triturus and 16 out of 17 salamandrid genera. A clade comprising all the "Newts" can be separated from the "True Salamanders" and Salamandrina clades. Within the "Newts" well-supported clades are: Tylototriton-Pleurodeles, the "New World Newts" (Notophthalmus-Taricha), and the "Modern Eurasian Newts" (Cynops, Pachytriton, Paramesotriton=together the "Modern Asian Newts", Calotriton, Euproctus, Neurergus and Triturus species). We found that Triturus is a non-monophyletic species assemblage, which includes four groups that are themselves monophyletic: (i) the "Large-Bodied Triturus" (six species), (ii) the "Small-Bodied Triturus" (five species), (iii) T. alpestris and (iv) T. vittatus. We estimated that the last common ancestor of Triturus existed around 64 million years ago (mya) while the root of the Salamandridae dates back to 95 mya. This was estimated using a fossil-based molecular dating approach and an explicit framework to select calibration points that least underestimated their corresponding nodes. Using the molecular phylogeny we mapped the evolution of life history and courtship traits in Triturus and found that several Triturus-specific courtship traits evolved independently.