Interspecific introgression of MHC genes in Triturus newts: Evidence from multiple contact zones.

The major histocompatibility complex (MHC) genes are central to the adaptive immune response in vertebrates. Selection generally maintains high MHC variation because the spectrum of recognized pathogens depends on MHC polymorphism. Novel alleles favoured by selection originate by interallelic recombination or de novo mutations but may also be acquired by introgression from related species. However, the extent and prevalence of MHC introgression remain an open question. In this study, we tested for MHC introgression in six hybrid zones formed by six Triturus newt species. We sequenced and genotyped the polymorphic second exons of the MHC class I and II genes and compared their interspecific similarity at various distances from the centre of the hybrid zone. We found evidence for introgression of both MHC classes in the majority of examined hybrid zones, with support for a more substantial class I introgression. Furthermore, the overall MHC allele sharing outside of hybrid zones was elevated between pairs of Triturus species with abutting ranges, regardless of the phylogenetic distance between them. No effect of past hybrid zone movement on MHC allele sharing was found. Finally, using previously published genome-wide data, we demonstrated that MHC introgression was more extensive than genome-wide introgression, supporting its adaptive potential. Our study thus provides evidence for the prevalence of MHC introgression across multiple Triturus hybrid zones, indicating that MHC introgression between divergent hybridizing species may be widespread and adaptive.

Reliability of environmental DNA surveys to detect pond occupancy by newts at a national scale.

The distribution assessment and monitoring of species is key to reliable environmental impact assessments and conservation interventions. Considerable effort is directed towards survey and monitoring of great crested newts (Triturus cristatus) in England. Surveys are increasingly undertaken using indirect methodologies, such as environmental DNA (eDNA). We used a large data set to estimate national pond occupancy rate, as well as false negative and false positive error rates, for commercial eDNA protocols. Additionally, we explored a range of habitat, landscape and climatic variables as predictors of pond occupancy. In England, 20% of ponds were estimated to be occupied by great crested newts. Pond sample collection error rates were estimated as 5.2% false negative and 1.5% false positive. Laboratory error indicated a negligible false negative rate when 12 qPCR replicates were used. Laboratory false positive error was estimated at 2% per qPCR replicate and is therefore exaggerated by high levels of laboratory replication. Including simple habitat suitability variables into the model revealed the importance of fish, plants and shading as predictors of newt presence. However, variables traditionally considered as important for newt presence may need more precise and consistent measurement if they are to be employed as reliable predictors in modelling exercises.

Triturus newts.

Ben Wielstra introduces Eurasian newts of the genus Triturus.

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.

Investigation of axonal regeneration of Triturus ivanbureschi by using physiological and proteomic strategies.

Peripheral nerve injuries are frequently observed and successful treatment depends mainly on the injury type, location of the damage, and the elapsed time prior to treatment. The regenerative capacity is limited only to the embryonic period in many mammalian tissues, but urodele amphibians do not lose this feature during adulthood. The main purpose of this study is to define the recovery period after serious sciatic nerve damage of a urodele amphibian, Triturus ivanbureschi. Experimental transection damage was performed on the sciatic nerves of T. ivanbureschi specimens. The recovery period of sciatic nerves were investigated by walking track analysis, electrophysiological recordings, and bottom-up proteomic strategies at different time points during a 35-day period. A total of 34 proteins were identified related to the nerve regeneration process. This study showed that the expression levels of certain proteins differ between distal and proximal nerve endings during the regeneration period. In distal nerve stumps, transport proteins, growth factors, signal, and regulatory molecules are highly expressed, whereas in proximal nerve stumps, neurite elongation proteins, and cytoskeletal proteins are highly expressed.

Oxidative cost of interspecific hybridization: a case study of two Triturus species and their hybrids.

Oxidative stress has most recently been suggested as one of the possible mechanisms responsible for reduced fitness of hybrids. To explore possible oxidative cost of hybridization, we examined anti-oxidant defence system parameters (superoxide dismutase, catalase, glutathione peroxidase, glutathione S-transferase, glutathione reductase, glutathione, SH groups), their interconnectedness (index of integration) and levels of oxidative damage [concentrations of lipid peroxides, TBARS (thiobarbituric acid reactive substances)] in laboratory-reared newt species, Triturus macedonicus and T. ivanbureschi, and their hybrid. Our results showed that parental species differed in anti-oxidant defence system parameters, but not in the levels of integration of the whole system and oxidative damage. Individuals of T. ivanbureschi had higher activities of superoxide dismutase, glutathione S-transferase and concentrations of glutathione. Hybrid individuals of crested newts displayed higher levels of the anti-oxidant defence system (higher superoxide dismutase, catalase, glutathione peroxidase activities and concentrations of SH groups), and a lower overall correlation of anti-oxidant system (lower index of integration) in comparison with both parental species, suggesting that they may possess a less efficient anti-oxidant defence system and a higher investment in maintaining oxidative balance. The higher investment in the anti-oxidant system could divert limited resources away from other functions and affect further hybrid fitness. The presented findings contribute to a better understanding of the anti-oxidant defence system of crested newts and their interspecies differences, and support the hypothesis that oxidative stress is one of the costs of interspecific hybridization.

Seasonal variation in environmental DNA detection in sediment and water samples.

The use of aquatic environmental DNA (eDNA) to detect the presence of species depends on the seasonal activity of the species in the sampled habitat. eDNA may persist in sediments for longer than it does in water, and analysing sediment could potentially extend the seasonal window for species assessment. Using the great crested newt as a model, we compare how detection probability changes across the seasons in eDNA samples collected from both pond water and pond sediments. Detection of both aquatic and sedimentary eDNA varied through the year, peaking in the summer (July), with its lowest point in the winter (January): in all seasons, detection probability of eDNA from water exceeded that from sediment. Detection probability of eDNA also varied between study areas, and according to great crested newt habitat suitability and sediment type. As aquatic and sedimentary eDNA show the same seasonal fluctuations, the patterns observed in both sample types likely reflect current or recent presence of the target species. However, given the low detection probabilities found in the autumn and winter we would not recommend using either aquatic or sedimentary eDNA for year-round sampling without further refinement and testing of the methods.

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.

The detection of great crested newts year round via environmental DNA analysis.

OBJECTIVE: Analysis of environmental DNA (eDNA) is a method that has been used for the detection of various species within water bodies. The great crested newt (Triturus cristatus) has a short eDNA survey season (mid-April to June). Here we investigate whether this season could be extended into other months using the current methodology as stipulated by Natural England. RESULTS: Here we present data to show that in monthly water samples taken from two ponds (March 2014-February 2015) we were able to detect great crested newt DNA in all months in at least one of the ponds. Similar levels of great crested newt eDNA (i.e. highly positive identification) were detected through the months of March-August, suggesting it may be possible to extend the current survey window. In order to determine how applicable these observations are for ponds throughout the rest of the UK, further work in multiple other ponds over multiple seasons is suggested. Nevertheless, the current work clearly demonstrates, in two ponds, the efficacy and reproducibility of eDNA detection for determining the presence of great crested newts.

Vertebral shape and body elongation in Triturus newts.

Body elongation in vertebrates is often related to a lengthening of the vertebrae and an increase in their number. Changes in the number and shape of vertebrae are not necessarily linked. In tailed amphibians, a change in body shape is mostly associated with an increase in the number of trunk and tail vertebrae. Body elongation without a numerical change of vertebrae is rare. In Triturus aquatic salamanders body elongation is achieved by trunk elongation through an increase in the number of trunk vertebrae. We used computed microtomography and three-dimensional geometric morphometrics to document the size, shape and number of trunk vertebrae in seven Triturus species. The data suggest that body elongation has occurred more frequently than body shortening, possibly related to a more aquatic versus a more terrestrial locomotor style. Our results show that body elongation is achieved through an increase in the number of trunk vertebrae, and that interspecific differences in vertebral shape are correlated with this pattern of elongation. More gracile trunk vertebrae were found in the more elongated species. The shape differences are such that single trunk vertebrae can be used for the identification of species with a possible application in the identification of subfossil and fossil material.

Investigation of histone H4 hyperacetylation dynamics in the 5S rRNA genes family by chromatin immunoprecipitation assay.

Oogenesis is a critical event in the formation of female gamete, whose role in development is to transfer genomic information to the next generation. During this process, the gene expression pattern changes dramatically concomitant with genome remodelling, while genomic information is stably maintained. The aim of the present study was to investigate the presence of H4 acetylation of the oocyte and somatic 5S rRNA genes in Triturus cristatus, using chromatin immunoprecipitation assay (ChIP). Our findings suggest that some epigenetic mechanisms such as histone acetylation could be involved in the transcriptional regulation of 5S rRNA gene families.

Micronucleus test on Triturus carnifex as a tool for environmental biomonitoring.

The amphibian micronucleus test has been widely used during the last 30 years to test the genotoxic properties of several chemicals and as a tool for ecogenotoxic monitoring. The vast majority of these studies were performed on peripheral blood of urodelan larvae and anuran tadpoles and to a lesser extent adults were also used. In this study, we developed protocols for measuring micronuclei in adult shed skin cells and larval gill cells of the Italian crested newt (Triturus carnifex). Amphibians were collected from ponds in two protected areas in Italy that differed in their radon content. Twenty-three adult newts and 31 larvae were captured from the radon-rich pond, while 20 adults and 27 larvae were taken from the radon-free site. The animals were brought to the laboratory and the micronucleus test was performed on peripheral blood and shed skins taken from the adults and on larval gills. Samples from the radon-rich site showed micronucleus frequencies higher than those from the radon-free site and the difference was statistically significant in gill cells (P < 0.00001). Moreover, the larval gills seem to be more sensitive than the adult tissues. This method represents an easy (and noninvasive in the case of the shed skin) application of the micronucleus assay that can be useful for environmental studies in situ.

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.