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.

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.

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.

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.

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.

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.

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.

The balanced lethal system of crested newts: a ghost of sex chromosomes past?

Balanced lethal systems are more than biological curiosities: as theory predicts, they should quickly be eliminated through the joint forces of recombination and selection. That such systems might become fixed in natural populations poses a challenge to evolutionary theory. Here we address the case of a balanced lethal system fixed in crested newts and related species, which makes 50% of offspring die early in development. All adults are heteromorphic for chromosome pair 1. The two homologues (1A and 1B) have different recessive deleterious alleles fixed on a nonrecombining segment, so that heterozygotes are viable, while homozygotes are lethal. Given such a strong segregation load, how could autosomes stop recombining? We propose a role for a sex-chromosome turnover from pair 1 (putative ancestral sex chromosome) to pair 4 (currently active sex chromosome). Accordingly, 1A and 1B represent two variants (Y(A) and Y(B)) of the Y chromosome from an ancestral male-heterogametic system. We formalize a scenario in which turnovers are driven by sex ratio selection stemming from gene-environment interactions on sex determination. Individual-based simulations show that a balanced lethal system can be fixed with significant likelihood, provided the masculinizing allele on chromosome 4 appears after the elimination of the feminizing allele on chromosome 1. Our study illustrates how strikingly maladaptive traits might evolve through natural selection.

Monitoring endangered freshwater biodiversity using environmental DNA.

Freshwater ecosystems are among the most endangered habitats on Earth, with thousands of animal species known to be threatened or already extinct. Reliable monitoring of threatened organisms is crucial for data-driven conservation actions but remains a challenge owing to nonstandardized methods that depend on practical and taxonomic expertise, which is rapidly declining. Here, we show that a diversity of rare and threatened freshwater animals--representing amphibians, fish, mammals, insects and crustaceans--can be detected and quantified based on DNA obtained directly from small water samples of lakes, ponds and streams. We successfully validate our findings in a controlled mesocosm experiment and show that DNA becomes undetectable within 2 weeks after removal of animals, indicating that DNA traces are near contemporary with presence of the species. We further demonstrate that entire faunas of amphibians and fish can be detected by high-throughput sequencing of DNA extracted from pond water. Our findings underpin the ubiquitous nature of DNA traces in the environment and establish environmental DNA as a tool for monitoring rare and threatened species across a wide range of taxonomic groups.

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.

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.

Unraveling the rapid radiation of crested newts (Triturus cristatus superspecies) using complete mitogenomic sequences.

BACKGROUND: The rapid radiation of crested newts (Triturus cristatus superspecies) comprises four morphotypes: 1) the T. karelinii group, 2) T. carnifex - T. macedonicus, 3) T. cristatus and 4) T. dobrogicus. These vary in body build and the number of rib-bearing pre-sacral vertebrae (NRBV). The phylogenetic relationships of the morphotypes have not yet been settled, despite several previous attempts, employing a variety of molecular markers. We here resolve the crested newt phylogeny by using complete mitochondrial genome sequences. RESULTS: Bayesian inference based on the mitogenomic data yields a fully bifurcating, significantly supported tree, though Maximum Likelihood inference yields low support values. The internal branches connecting the morphotypes are short relative to the terminal branches. Seen from the root of Triturus (NRBV = 13), a basal dichotomy separates the T. karelinii group (NRBV = 13) from the remaining crested newts. The next split divides the latter assortment into T. carnifex - T. macedonicus (NRBV = 14) versus T. cristatus (NRBV = 15) and T. dobrogicus (NRBV = 16 or 17). CONCLUSIONS: We argue that the Bayesian full mitochondrial DNA phylogeny is superior to previous attempts aiming to recover the crested newt species tree. Furthermore, our new phylogeny involves a maximally parsimonious interpretation of NRBV evolution. Calibrating the phylogeny allows us to evaluate potential drivers for crested newt cladogenesis. The split between the T. karelinii group and the three other morphotypes, at ca. 10.4 Ma, is associated with the separation of the Balkan and Anatolian landmasses (12-9 Ma). No currently known vicariant events can be ascribed to the other two splits, first at ca. 9.3 Ma, separating T. carnifex - T. macedonicus, and second at ca. 8.8 Ma, splitting T. cristatus and T. dobrogicus. The crested newt morphotypes differ in the duration of their annual aquatic period. We speculate on the role that this ecological differentiation could have played during speciation.

Pituitary adenylate cyclase-activating polypeptide and its receptor PAC1 in the testis of Triturus carnifex and Podarcis sicula.

The pituitary adenylate cyclase-activating polypeptide (PACAP) is a member of the glucagon-related family that occurs in two amidated forms with 38 (PACAP38) and 27 (PACAP27) amino acids. First discovered in the brain, it was then localized in several peripheral tissues of mammals, including the testis. However, current knowledge of the expression and function of PACAP and its receptor PAC(1) in the reproductive system of non-mammalian vertebrates, and particularly in the testis, is still limited. The aim of this work was to study the presence of PACAP and its receptor PAC(1) in the testis of two non-mammalian vertebrates during the breeding season: the crested newt Triturus carnifex and the wall lizard Podarcis sicula. The expression and distribution of this neuropeptide and its receptor PAC(1) were investigated by using in situ hybridization and immunohistochemistry techniques. Our results demonstrated that PACAP and its receptor PAC(1) were highly represented in the testis of these two species. In particular, we showed that they are present within some germ cells and that PACAP, unlike in mammals, is expressed also in the somatic cells (Sertoli and Leydig cells) of the testis of these two non-mammalian vertebrates, suggesting that this neuropeptide is involved in the hormonal control of spermatogenesis and steroidogenesis.

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.

Triturus newts.

Ben Wielstra introduces Eurasian newts of the genus 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.

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.

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.