Evolutionary history of Ichthyosaura alpestris (Caudata, Salamandridae) inferred from the combined analysis of nuclear and mitochondrial markers

Abstract

Widespread species with morphologically and ecologically differentiated populations are key to understand speciation because they allow investigating the different stages of the continuous process of population divergence. The alpine newt, Ichthyosaura alpestris, with a range that covers a large part of Central Europe as well as isolated regions in all three European Mediterranean peninsulas, and with strong ecological and life-history differences among populations, is an excellent system for such studies. We sampled individuals across most of the range of the species, and analyzed mitochondrial (1442bp) and nuclear (two nuclear genes -1554bp- and 35 allozyme loci) markers to produce a time-calibrated phylogeny and reconstruct the historical biogeography of the species. Phylogenetic analyses of mtDNA data produced a fully resolved topology, with an endemic, Balkan clade (Vlasina) which is sister to a clade comprising an eastern and a western group. Within the former, one clade (subspecies I. a. veluchiensis) is sister to a clade containing subspecies I. a. montenegrina and I. a. serdara as well as samples from southern Romania, Bosnia-Herzegovina, Serbia and Bulgaria (subspecies I. a. reiseri and part of I. a. alpestris). Within the western group, populations from the Italian peninsula (subspecies I. a. apuana and I. a. inexpectata) are sister to a clade containing samples from the Iberian Peninsula (subspecies I. a. cyreni) and the remainder of the samples from subspecies I. a. alpestris (populations from Hungary, Austria, Poland, France, Germany and the larger part of Romania). Results of ∗BEAST analyses on a combined mtDNA and nDNA dataset consistently recovered with high statistical support four lineages with unresolved inter-relationships: (1) subspecies I. a. veluchiensis; (2) subspecies I. a. apuana+I. a. inexpectata; (3) subspecies I. a. cyreni+part of subspecies I. a. alpestris (the westernmost populations, plus most Romanian populations); and (4) the remaining populations, including subspecies I. a. serdara, I. a. reiseri and I. a. montenegrina and part of subspecies I. a. alpestris, plus samples from Vlasina. Our time estimates are consistent with ages based on the fossil record and suggest a widespread distribution for the I. alpestris ancestor, with the split of the major eastern and western lineages during the Miocene, in the Tortonian. Our study provides a solid, comprehensive background on the evolutionary history of the species based on the most complete combined (mtDNA+nDNA+allozymes) dataset to date. The combination of the historical perspective provided by coalescent-based analyses of mitochondrial and nuclear DNA variation with individual-based multilocus assignment methods based on multiple nuclear markers (allozymes) also allowed identification of instances of discordance across markers that highlight the complexity and dynamism of past and ongoing evolutionary processes in the species.