Ribs ‘n’ genes: Triturus hybrid zones

Hybrid zones are the regions where different species meet, mate and produce offspring. In a study just published in the Biological Journal of the Linnean Society we explore nine Triturus hybrid zones (all but one existing in nature), using a large amount of genetic markers (forty allozymes and a mitochondrial gene) and a highly informative morphological character (the number of rib-bearing vertabrae), for a huge number of newts (well over 700). Hybrid zones are narrow and individuals show less hybridity than would be expected under random mating, which shows that the different species are genetically pretty isolated. The degree of genetic isolation increases with genetic divergence. We do find mitochondrial DNA in the ‘wrong’ species, which illustrates that individual genes can be exchanged between species that for the rest manage to remain distinct. The number of ribs is a good indicator of species identity. Although there is some variation within species, most of this occurs close to hybrid zones, suggesting hybridization is to blame here. This synthesis marks the end of an era in a way: a lot of the data was already collected decades ago. It is also a strong base for future research on the Triturus hybrid zones using the new Ion Torrent protocol.

Reference: Arntzen, J.W., Wielstra, B., Wallis, G.P. (2014). The modality of nine Triturus newt hybrid zones, assessed with nuclear, mitochondrial and morphological data. Biological Journal of the Linnaean Society 113(2): 604-622.

newton-logo

I conducted this work as a Newton International Fellow.
Posted in Uncategorized | 1 Comment

The final nail in the coffin of Triturus arntzeni

During my work on the crested newt traditionally referred to as Triturus karelinii it soon became clear that more than one species is involved (although sorting out the details took a bit of time and is still not quite finished). The name arntzeni had been used before to refer to newts from the Balkan part of the range. This name was proposed in honor of my mentor Pim Arntzen. I never bought it. Not the part of Pim deserving his own newt – he most certainly does. But I did not believe that the type locality, near the small village Vrtovać in Serbia, belonged to Triturus karelinii sensu lato. It was located in such an inconvenient place, in a region where the distribution was not well understood but where I strongly suspected Triturus macedonicus would occur. We went and visit the type locality during a field trip in 2010 with Jelka Crnobrnja-Isailović. We had been joking about the identification of these newts beforehand and when we walked up to the pond judgement day had arrived. As usual there was a bit of competition of who would catch the first crested newt and as usual it was Pim who netted the first newt. I don’t remember the exact words he used when he took it out of the net but it was pretty dry and went something like “Well… that’s macedonicus…”. I have to admit I thought it all pretty funny.

Pim at the pond

Pim at the pond

3b

These sedated newts just don’t look like karelinii

Last year we published a paper in Zootaxa where we reviewed all the available data and made clear that newts from the type locality of arntzeni strongly resembled macedonicus. So basically Pim burned his own name! Spartak Litvinchuk provided some tissue from the actual type material but this was too degraded to sequence DNA. Although DNA data from fresh material pointed towards macedonicus, arguably the amount of markers used was a bit low (and note that mitochondrial DNA is pretty useless in this region due to wide scale introgression). However, with the new Triturus Ion Torrent protocol sequencing a lot of markers for poorly preserved material suddenly became possible. In a follow-up paper, again in Zootaxa, we compared high quality genetic profiles for nine newts from the arntzeni type locality (including the holotype to which, according to the rules of taxonomy, the name is officially attached) to a comprehensive sample including all crested newt species. We got an unambiguous picture of the genetic ancestry of the Vrtovać newts: as suspected they were mostly macedonicus, but there was also a considerable amount of karelinii genes present, evidencing hybridization. These finding strongly support that the name arntzeni is not valid. Instead we introduced the name ivanbureschi and this time we made sure the type locality is nowhere near a contact zone with another species of crested newt!

3c

This figure might look a bit complex but sums up the genetic identity of the Vrtovać newts. The bottom left corner represents a pure macedonicus genotype and the bottom right pure karelinii (now ivanbureschi). The top of the triangle represents pure F1 (first generation) hybrids. This plot is made with the HIest R package by Ben Fitzpatrick.
By the way, this pond is the type locality for ivanbureschi, we didn't put this picture in the original paper, not sure why actually.

By the way, this pond is the type locality for ivanbureschi. We didn’t put this picture in the original paper, not sure why actually.

Reference: Wielstra, B., Arntzen, J.W. (2014). Kicking Triturus arntzeni when it’s down: large-scale nuclear genetic data confirm that newts from the type locality are genetically admixed. Zootaxa 3802(3): 381-388.

Reference: Wielstra, B., Litvinchuk, S.N., Naumov, B., Tzankov, N., Arntzen, J.W. (2013). A revised taxonomy of crested newts in the Triturus karelinii group (Amphibia: Caudata: Salamandridae), with the description of a new species. Zootaxa 3682(3): 441-453.

newton-logo

I conducted this work as a Newton International Fellow.
Posted in Newts | 1 Comment

Triturus goes genomics, kind of

You can only learn so much about a study system if you have few genetic markers available. Particularly if that study system has an extensive history of hybridization, as is the case for Triturus. Because salamanders have massive and complex genomes it is not possible to simply sequence one. Not yet at least. Luckily we had some Triturus transcriptome data laying around. The transcriptome contains messenger DNA – the transcripts of functional genes – but not the long introns and endless repeats that torment salamander genomes. Hence it provides relatively simple, genome-wide reference data for marker design. After designing a large set of markers, I tested which of these worked for all crested newt species and multiplexed the successful ones for a large set of newts. Next the whole bunch was sequenced at Naturalis’ next-generation sequencing facility on an Ion Torrent machine and subsequently the huge amount of genetic data was sorted out with a bioinformatics pipeline and converted it to a workable format. It sounds easy, but it was quite an undertaking.

?????????????????

The top picture shows the localities sampled (circled numbers, three newts each) and the bottom picture shows the probability with which they belong to their own species (all 1, as would be expected if our method works).

Using data for a set of newts representing all species and also some putative hybrids we showed that the protocol can be used to allocate individuals to the proper species and pick out genetically admixed newts. Although the amount of markers is still relatively modest, the data provide a very detailed picture on the genetic composition of crested newts. A paper describing the methodology has recently been published in Molecular Ecology Resources. Basically we can now provide a detailed picture of the distribution of the different species and the genetic distinction of and gene flow between these species: an important next step in my research. Meanwhile I have sequenced about 1500 individuals throughout the range of Triturus as the basis for quite some papers to come. Stay tuned!

2b

These are a couple of newts from the population marked with an H on the map above.  You can click on the picture to see a larger version. The throat and belly pattern is very variable in this populations. Some animals look more like macedonicus or ivanbureschi and the ones depicted here look especially messy. Based on the Ion Torrent data these newts indeed show mixed macedonicus and ivanbureschi genetic ancestry and almost all are identified as backcrosses towards ivanbureschi (with the remainder being F2 hybrids).

By the way, many thanks to Wieslaw Babik, Michał Stuglik and Piotr Zieliński from Jagiellonian University, Krakow, Poland for helping with the design of this protocol!

Reference: Wielstra, B., Duijm, E., Lagler, P., Lammers, Y., Meilink, W., Ziermann, J.M., Arntzen, J.W. (2014). Parallel tagged amplicon sequencing of transcriptome-based genetic markers for Triturus newts with the Ion Torrent next-generation sequencing platform. Molecular Ecology Resources 14(5): 1080-1089.

newton-logo

I conducted this work as a Newton International Fellow.
Posted in Newts | 12 Comments

Mitochondrial mess-up

The different crested newt species typically have very distinct mitochondrial DNA. However, as these guys hybridize like rabbits, often the mitochondrial DNA of one species has locally been transferred (introgressed) into the wrong species. ‘Locally’ can refer to quite an extensive range actually. Species distribution modelling of distinct geographical populations, only recognizable based on DNA divergence, is a good way to test whether they represent cryptic species – multiple species erroneously classified as a single one because no one previously realized they were distinct. However, if you would use mitochondrial DNA to identify the putative species, introgression would cause misidentifications and could lead to under- or overestimation of the ecological niche of the putative species. We explore this problem in a paper just out in PLoS ONE.

?????????????????

This figure shows the setting in which we test the effect of introgressed mitochondrial DNA on niche estimation. The two areas with red dot on the wrongly colored background are causing trouble. If using mitochondrial DNA as a guidance, you would interpret these localities as beloning to the red species, while excluding them from the green or blue species.

Reference: Wielstra, B., Arntzen, J.W. (2014). Exploring the effect of asymmetric mitochondrial DNA introgression on estimating niche divergence in morphologically cryptic species. PLoS ONE 9(4): e95504.

newton-logo

I conducted this work as a Newton International Fellow.
Posted in Newts | 2 Comments

Climate and phylogenetic signal

In a paper published in Annales Zoologici Fennici we tested the influence of shared ancestry and climate on the variation in life history traits in Triturus newts. Some variation can be explained by climate, some by shared ancestry and some by both.

Reference: Vukov, T.D., Cvijanović, M., Wielstra, B., Kalezić M.L. (2014). The roles of phylogeny and climate in shaping the variation in life-history traits observed in the newt genus Triturus (Caudata, Salamandridae). Annales Zoologici Fennici 51(5): 445-456.

Posted in Uncategorized | Leave a comment

Newton International Fellowship

My application for a Newton International Fellowship has been succesful! I will conduct this two-year postdoctoral fellowship, using the Ion Torrent protocol I designed during my postdoc at Naturalis to study hybridzones in Triturus, in the lab of Professor Terry Burke at the University of Sheffield.

newton-logo

Posted in Uncategorized | Leave a comment

Cryptic crested newt species

The crested newt traditionally referred to as ‘Triturus karelinii’ comprises three distinct mitochondrial DNA clades. These clades are found in the east, the centre and the west of the range. The difference between mitochondrial DNA clades is comparable to the difference between the mitochondrial DNA of recognized crested newt species. We wanted to see if the three mitochondrial DNA clades might in fact represent different species.

BAPSThis plot shows individuals (thin bars) within populations (thick bars) roughly ordered from west to east. Based on their genotype, individuals are allocated (0-100%) to three geographical genetic groups (represented by different colors).

In a paper published in Molecular Phylogenetics and Evolution we show that, based on three nuclear genes, there are three discrete geographical groups, in line with the three species hypothesis. We suggest that these three groups should be considered distinct species, but as it is as yet unclear if they can be distinguished based on morphology, we refer to them as ‘cryptic species’ for now.

mapmpe13Here you see the three cryptic species that make up the T. karelinii-group of crested newts. An intriguing finding is that asymmetric DNA introgression from the western group into the central group (the red-green hatched area). We suggest that this pattern can be explained by the central group having expanded its range at the expense of the western group, while the two hybridized in the process. An interesting hypothesis to test in a future study!

Reference: Wielstra, B., Baird, A.B., Arntzen, J.W. (2013). A multimarker phylogeography of crested newts (Triturus cristatus superspecies) reveals cryptic species. Molecular Phylogenetics and Evolution 67(1): 167-175.

Posted in Uncategorized | 4 Comments

Mitochondrial DNA capture in smooth newts

Hybridizing species sometimes exchange genes in nature (a phenomenon called introgression). Introgression has particularly been documented for mitochondrial DNA. This might mean mitochondrial DNA is more susceptible to introgression, but it is also the case that researchers have simply studied mitochondrial DNA more. Usually introgression is restricted to close to the hybrid zone. There are however more pronounced examples. It does not get more extreme than in smooth newts, where the original mitochondrial DNA of Lissotriton montandoni has been completely replaced by that of L. vulgaris.

smooth carpathian

A male Lissotriton montandoni on the left and a male L. vulgaris on the right. Pictures by team Babik.

In a paper published in Molecular Ecology we extensively sample both mitochondrial and nuclear DNA for L. montandoni and surrounding L. vulgaris populations. We also use species distribution modelling to determine range dynamics of L. montandoni since the Last Glacial Maximum. We show that mitochondrial DNA introgression occurred several times, and at different moments, in different parts of the range of L. montandoni. This is probably related to the inferred range fragmentation under glacial conditions, and  independent range expansion from these range fragments into L. vulgaris territory when the climate ameliorated. In contrast, there is little evidence of recent nuclear gene flow between the species.

mtDNA capture

A visualization of complete mitochondrial DNA capture in response to species displacement and range reduction. We are dealing with two species here, let’s call them green and red. These species possess distinct mitochondrial DNA. Panels are ordered chronologically. In the panels, dots reflect localities, while boxes represent rough outlines of the ranges. Background shading (green or red) reflects species identity, while localities are colored according to the mitochondrial DNA type present (so dots are green or red).
In I) the ranges of the two species are geographically separated. At this stage, green mitochondrial DNA is only found in localities of the green species, and red mitochondrial DNA in localities of the red species. However, the green species is doing well for itself and, over the generations, its population increases. Surplus individuals from the right edge of the range start colonizing new localities further right. In these localities, numbers start increasing again and some offspring colonize new localities even further right, and so on. With time, the green species expands its range to the right, towards that of the red one.
In II) the ranges of the two species have come into contact. Where they meet, the two species start reproducing with one another (hybridization), resulting in offspring that are a mix of both green and red. Initially there are relatively few green individuals, but they keep pouring into the into the hybrid zone, while this is not the case for red individuals. Hybrids already present mate with these green individuals, and their offspring again mate with green individuals, and so on. Over time, individuals in the hybrid zone become ‘greener’. However, because at the initial stage of invasion by the green species most matings will concern individuals that also have red genes, some of these red genes rather than their green counterparts could locally get fixed in (introgress into) the green species by chance (and several processes that we won’t go into now might reinforce introgression). In this case mitochondrial DNA introgresses. Hence, at the right edge of the green species range, you start to see localities where individuals belong to the green species, but possess mitochondrial DNA typical of the red species (red dots in a green range).
In III) the individuals in the initial hybrid zone have become all green, except their mitochondrial DNA. Still the green species keeps expanding its range further to the right. Once again green individuals meet red ones, start hybridizing and gradually take over, and so on. In consequence, the hybrid zone between the two species moves towards the right, as the green species replaces the red species. Yet, the members of the green species at the frontier, as well as the red individuals they hybridize with, only possess red mitochondrial DNA. So the region to the right of the initial hybrid zone becomes ‘greener’, but the overturn between green and red mitochondrial DNA still aligns with that initial hybrid zone. One way to put this is that the red mitochondrial DNA ‘surfs the wave’ of the green species expansion. In effect you end up with red localities on a green background over a considerable area.
In IV) the hybrid zone between the green and the red species has stabilized at a region where the green species does not have an edge over the red species as it did before. At the same time the part of the green species’ range where it still possessed the green mitochondrial DNA type becomes unsuitable and here the green species goes extinct (black background). As a result, there are now only members of the green species left that possess red mitochondrial DNA, while the original green mitochondrial DNA has been lost.

Reference: Zieliński, P., Nadachowska-Brzyska, K., Wielstra, B., Szkotak, R., Covaciu-Marcov, S., Cogălniceanu, D., Babik, W. (2013). No evidence for nuclear introgression despite complete mtDNA replacement in the Carpathian newt (Lissotriton montandoni). Molecular Ecology 22(7): 1884-1903.

Posted in Uncategorized | Leave a comment

Crested newt skulls

The Triturus karelinii-group of crested newts comprises three mitochondrial DNA lineages, but no morphological differences are known. In a paper published in Zoologischer Anzeiger we analyse skull shape to see if there are differences between lineages. While we do not find three discrete geographical groups, the observed differences among lineages are also not smaller than between lineages and another crested newt species, T. macedonicus. We just need to keep looking for morphological differences.

Reference: Ivanović, A., Üzüm, N., Wielstra, B., Olgun, K., Litvinchuk, S.N., Kalezić, M.L., Arntzen, J.W. (2013) Is mitochondrial DNA divergence of Near Eastern crested newts (Triturus karelinii group) reflected by differentiation of skull shape? Zoologischer Anzeiger 252(2): 269-277.

Posted in Uncategorized | Leave a comment

Tracing glacial refugia of crested newts

The Quaternary Ice Age heavily influenced the distribution of species. During the colder glacial periods, species went extinct in part of their range, while during warmer interglacial periods, they could recolonize these regions again. This contraction-expansion pattern left its mark on genetic diversity across species’ ranges, with high diversity in regions where species survived continuously, and low diversity in regions where they were periodically wiped out. Past range shifts can also be visualized by projecting species distribution models, based on the environmental conditions currently experienced, on climate reconstructions of the past.

In a study published in Frontiers in Zoology we explore how all the marbled and crested newts species (so the entire genus Triturus) responded to the climate change associated with the Ice Age. We conduct a phylogeographical survey, meaning we sequence a lot of mitochondrial DNA, for many populations throughout each of the species ranges, and look at variation in genetic diversity across species ranges. Additionally, we compare species distribution models projected on current and on past climate layers (the Last Glacial Maximum, about 21,000 years ago).

Fzoolscenario

A visualization of the biogeographical scenario proposed, showing the positions of glacial refugia in dark shades and the regions postglacially colonized in light shades for each individual species. In grey areas we infer that one species displaced another as they shifted their ranges.

By combining the two independent techniques of phylogeography and species distribution modelling, we obtain a more complete understanding of the historical biogeography of the crested and marbled newts than both approaches would have provided on their own.

Reference: Wielstra, B., Crnobrnja-Isailović, J., Litvinchuk, S.N., Reijnen, B., Skidmore, A.K., Sotiropoulos, K., Toxopeus, A.G., Tzankov, N., Vukov, T., Arntzen, J.W. (2013). Tracing glacial refugia of Triturus newts based on mitochondrial DNA phylogeography and species distribution modeling. Frontiers in Zoology 10: 13.

Posted in Uncategorized | 2 Comments