Phylogenomics of the adaptive radiation of Triturus newts supports gradual ecological niche expansion towards an incrementally aquatic lifestyle

A staging table of Balkan crested newt embryonic development to serve as a baseline in evolutionary developmental studies

There is an increased interest in the evolution and development of newts from the genus Triturus because: (1) morphological differentiation among the nine constituent species largely corresponds to different ecological preferences, (2) hybridization between different species pairs has various evolutionary outcomes in terms of life history traits and morphology, and (3) the genus expresses a balanced lethal system that causes arrested growth and death of half of the embryos. These features provide natural experimental settings for molecular, morphological, and life-history studies. Therefore, we produce a staging table for the Balkan crested newt (T. ivanbureschi). We provide detailed descriptions of 34 embryonic stages based on easily observable and interpretable external morphological characters, to ensure reproducibility. Compared with previous staging tables for Triturus, we include a vastly increased sample size and provide high-resolution photographs in lateral, ventral, and dorsal view, complemented by videos of specific developmental periods, and accompanied by detailed explanations on how to delineate the specific stages. Our staging table will serve as a baseline in comparative studies on Triturus newts: an emerging model system in evolutionary and developmental studies.

Spatial genetic structure in European marbled newts revealed with target enrichment by sequence capture

European marbled newts come in two species that have abutting ranges. The northern species, Triturus marmoratus, is found in France and the northern part of the Iberian Peninsula, whereas the southern species, T. pygmaeus, is found in the southwestern corner of the Iberian Peninsula. We study the intraspecific genetic differentiation of the group because morphological data show geographical variation and because the Iberian Peninsula is a recognized center of speciation and intraspecific genetic diversity for all kinds of organisms, amphibians included. We use target enrichment by sequence capture to generate c. 7 k nuclear DNA markers. We observe limited genetic exchange between the species, which confirms their distinctiveness. Both species show substantial genetic structuring that is only in part mirrored by morphological variation. Genetically differentiated groups are found in the south (T. marmoratus) and west (T. pygmaeus) of the species ranges. Our observations highlight the position of the Iberian Peninsula as a hotspot for genetic differentiation.

Piecing the barcoding puzzle of Palearctic water frogs (Pelophylax) sheds light on amphibian biogeography and global invasions

Palearctic water frogs (genus Pelophylax) are an outstanding model in ecology and evolution, being widespread, speciose, either threatened or threatening to other species through biological invasions, and capable of siring hybrid offspring that escape the rules of sexual reproduction. Despite half a century of genetic research and hundreds of publications, the diversity, systematics and biogeography of Pelophylax still remain highly confusing, in no small part due to a lack of correspondence between studies. To provide a comprehensive overview, we gathered >13,000 sequences of barcoding genes from >1700 native and introduced localities and built multigene mitochondrial (~17 kb) and nuclear (~10 kb) phylogenies. We mapped all currently recognized taxa and their phylogeographic lineages (>40) to get a grasp on taxonomic issues, cyto-nuclear discordances, the genetic makeup of hybridogenetic hybrids, and the origins of introduced populations. Competing hypotheses for the molecular calibration were evaluated through plausibility tests, implementing a new approach relying on predictions from the anuran speciation continuum. Based on our timetree, we propose a new biogeographic paradigm for the Palearctic since the Paleogene, notably by attributing a prominent role to the dynamics of the Paratethys, a vast paleo-sea that extended over most of Europe. Furthermore, our results show that distinct marsh frog lineages from Eastern Europe, the Balkans, the Near East, and Central Asia (P. ridibundus ssp.) are naturally capable of inducing hybridogenesis with pool frogs (P. lessonae). We identified 14 alien lineages (mostly of P. ridibundus) over ~20 areas of invasions, especially in Western Europe, with genetic signatures disproportionally pointing to the Balkans and Anatolia as the regions of origins, in line with exporting records of the frog leg industry and the stocks of pet sellers. Pelophylax thus emerges as one of the most invasive amphibians worldwide, and deserves much higher conservation concern than currently given by the authorities fighting biological invasions.

Phylogenomics of the narrowly endemic Eurycheilichthys (Siluriformes: Loricariidae): Sympatric species with non-sister relationships suggest mainly allopatric speciation

Armored catfishes of the genus Eurycheilichthys are endemic to Southern Brazil and Misiones (Argentina) comprising nine species of small size, with a high degree of sympatry and species diversity distributed in two river basins. Here we use new genome-wide data to infer a species phylogeny and test species boundaries for this poorly known group. We estimate 1) the phylogenetic relationships of the species of Eurycheilichthys based on 29,350 loci in 65 individuals of nine species plus outgroups, and 2) the population structure and differentiation based on 43,712 loci and 62 individuals to estimate how geography may have acted on speciation and formation of the sympatric species groups. Analyses support the monophyly of the genus and suggest two species-inclusive clades (East and West) with high support and very recently diverged species. Western clade contains E. limulus (from upper Jacuí River basin) that is sister to Western species of the Taquari-Antas basin plus E. paucidens. The Eastern clade contains E. pantherinus (from Uruguay River basin) sister to the Eastern species of the Taquari-Antas basin E. coryphaenus, plus the central-distributed species E. planus and E. vacariensis, and the more widely-distributed species E. luisae. Eurycheilichthys luisae is not monophyletic and may contain one or more cryptic species or hybrid individuals. A stronger diversity on structure of lineages on the Taquari-Antas, when compared to upper Uruguay and Jacuí River basins, and the fact that most of the sympatrically distributed taxa have non-sister relationships suggest a scenario of mainly allopatric speciation and may indicate a more dynamic landscape with headwater capture events among these tributaries.

A synoptic review of the Amphibians of Iran: bibliography, taxonomy, synonymy, distribution, conservation status, and identification key to the eggs, larvae, and adults

This study provides an illustrated account, a comprehensive update of the systematics, and a bibliography of the 15 species of anurans in five families, eight genera; and of the six species of urodeles in two families, four genera in Iran. Bufonidae, with eight species, is the most diverse family; Salamandridae has five species and Ranidae has four species. This study also presents updated identification keys for the eggs, larvae, and metamorphosed amphibians of Iran. We designated specimen NMW 19855.1 as neotype of Pelophylax persicus (Schneider, 1799) comb. nov.. Along with distribution maps obtained from all the reliable localities and museum specimens known at this time, the modelled habitat of species, and for the first time, the National Red List of amphibians based on the IUCN red list categories and criteria. Based on our evaluation we propose to categorize Bufo eichwaldi, Paradactylodon persicus, Neurergus derjugini, and N. kaiseri as Vulnerable at National Red List, and to move Bufotes (Calliopersa) luristanicus, B. (C.) surdus, Firouzophrynus olivaceus, and Rana pseudodalmatina from the category of Least Concern (LC) to Near Threatened (NT). The National Red List of amphibians that we propose has significant implications for endangered species management and conservation. Forty-one percent of amphibian species in Iran are endemic to the country, and more than forty percent of the Iranian amphibians are at risk of extinction. Zagros Mountain forest and Hyrcaniain forests have more than 80% (i.e. 18 species) of the diversity of Iranian amphibians. A considerable amount of scientific literature published on Iranian amphibians in Persian language is not easily accessible to researchers outside Iran. This monograph attempts to remedy the situation and provides broader access to international herpetology. We recognize that taxonomy is always in a state of flux, and the names and synonymies used here reflect our current view.

Histological changes of the skin during postembryonic development of the crested newt Triturus ivanbureschi (Urodela, Salamandridae)

BACKGROUND

Amphibian skin has been studied for many decades, especially the metamorphic changes in the skin of frogs. Less attention has been paid to salamander skin. Here, we describe changes in the skin structure during postembryonic development in a salamandrid species, the Balkan crested newt Triturus ivanbureschi.

METHOD

Using traditional histological techniques we examined the skin in the trunk region of three premetamorphic larval stages (hatchling, mid larval and late larval) and two postmetamorphic stages (juvenile, just after metamorphosis, and adult).

RESULTS

In larval stages, skin consists only of the epidermis, which gradually develops from the single epithelial cell layer in hatchlings, to a stratified epidermis with gland nests and characteristic Leydig cells at the late larval stage. During metamorphosis, Leydig cells disappear, and the dermal layer develops. In postmetamorphic stages, skin is differentiated on stratified epidermis and the dermis with well-developed glands. Three types of glands were observed in the skin of the postmetamorphic stages: mucous, granular and mixed. Gland composition appears to be stage- and sex-specific, with juveniles and adult female being more similar to each other. In juveniles and adult female, there are a similar proportion of glands in both dorsal and ventral skin, whereas in adult male granular glands dominated the dorsal skin, while mixed glands dominated the ventral skin.

CONCLUSION

Our results provide a baseline for future comparative research of skin anatomy in salamanders.

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.

Higher temperature induces oxidative stress in hybrids but not in parental species: A case study of crested newts

Ectotherms are particularly sensitive to global warming due to their limited capacity to thermoregulate, which can impact their performance and fitness. From a physiological standpoint, higher temperatures often enhance biological processes that can induce the production of reactive oxygen species and result in a state of cellular oxidative stress. Temperature alters interspecific interactions, including species hybridization. Hybridization under different thermal conditions could amplify parental (genetic) incompatibilities, thus affecting a hybrid's development and distribution. Understanding the impact of global warming on the physiology of hybrids and particularly their oxidative status could help in predicting future scenarios in ecosystems and in hybrids. In the present study, we investigated the effect of water temperature on the development, growth and oxidative stress of two crested newt species and their reciprocal hybrids. Larvae of Triturus macedonicus and T. ivanbureschi, and their T. macedonicus-mothered and T. ivanbureschi-mothered hybrids were exposed for 30 days to temperatures of 19°C and 24°C. Under the higher temperature, the hybrids experienced increases in both growth and developmental rates, while parental species exhibited accelerated growth (T. macedonicus) or development (T. ivanbureschi). Warm conditions also had different effects on the oxidative status of hybrid and parental species. Parental species had enhanced antioxidant responses (catalase, glutathione peroxidase, glutathione S-transferase and SH groups), which allowed them to alleviate temperature-induced stress (revealed by the absence of oxidative damage). However, warming induced an antioxidant response in the hybrids, including oxidative damage in the form of lipid peroxidation. These findings point to a greater disruption of redox regulation and metabolic machinery in hybrid newts, which can be interpreted as the cost of hybridization that is likely linked to parental incompatibilities expressed under a higher temperature. Our study aims to improve mechanistic understanding of the resilience and distribution of hybrid species that cope with climate-driven changes.

Changes in thyroid histomorphology and thyroglobulin immunostaining upon exposure to thiourea in Triturus newts

Amphibians are useful bioindicators for monitoring aquatic health and the influence of xenobiotics such as endocrine disrupting chemicals. Because aquatic ecosystems experience the majority of global pollution, aquatic organisms are most exposed and vulnerable to endocrine disruptors. Furthermore, penetration of endocrine disruptors into aquatic organisms especially in amphibians is even easier because of more permeable skin, resulting in high bioavailability and bioaccumulation of chemicals. One of the most potent endocrine disruptors is thiourea, which chemically blocks the synthesis of thyroid hormones and prevents metamorphosis in amphibians. We investigated the influence of thiourea on histomorphology of the thyroid gland in Triturus newts at the metamorphic stage, when thyroid hormone concentrations should reach their maximum level. Chronic exposure to thiourea induced hypertrophy and hyperplasia of follicular cells as well as a significant reduction of interstitial tissue. The intensity of the thyroglobulin immunostaining signal significantly decreases upon chronic exposure to thiourea. Successful cross-reactivity of human primary antibody in immunochemical detection of thyroglobulin in Urodela confirms potential homology in thyroglobulin structure throughout the vertebrates.

Salamanders reveal novel trajectories of amphibian MHC evolution

Genes of the major histocompatibility complex (MHC) code for immune proteins that are crucial for pathogen recognition in vertebrates. MHC research in nonmodel taxa has long been hampered by its genomic complexity that makes the locus-specific genotyping challenging. The recent progress in sequencing and genotyping methodologies allows an extensive phylogenetic coverage in studies of MHC evolution. Here, we analyzed the peptide-binding region of MHC class I (MHC-I) in 30 species of salamanders from six families representative of Urodela phylogeny. This extensive dataset revealed an extreme diversity of MHC-I in salamanders, both in terms of sequence diversity (about 3000 variants) and architecture (2-22 gene copies per species). The signal of positive selection was moderate and consistent between both peptide-binding domains, but varied greatly between genera. Positions of positively selected sites mostly coincided with human peptide-binding sites, suggesting similar structural properties of MHC-I molecules across distant vertebrate lineages. Finally, we provided evidence for the common intraexonic recombination at MHC-I and for the role of life history traits in the processes of MHC-I expansion/contraction. Our study revealed novel evolutionary trajectories of amphibian MHC and it contributes to the understanding of the mechanisms that generated extraordinary MHC diversity throughout vertebrate evolution.

The Reproductive Success of Triturus ivanbureschi × T. macedonicus F1 Hybrid Females (Amphibia: Salamandridae)

Simple Summary

Two moderately related large-bodied newt species endemic to the Balkan Peninsula, the Balkan crested newt (Triturus ivanbureschi) and the Macedonian crested newt (T. macedonicus), coexist and hybridize in central Serbia. Many generations of mutual hybrid crossings and backcrossings with parental species shaped the genetic composition of hybrid populations. Natural populations have admixed nuclear DNA (nuDNA) of parental species and T. ivanbureschi mitochondrial DNA (mtDNA), which is usually maternally inherited. The mechanisms that direct gene flow and shape the first generations of hybrids could explain the formation of hybrid zones and their maintenance in nature. We followed and compared life history traits related to reproduction of the first generation of reciprocal hybrids obtained by experimental crossing. Our results suggested that possible incompatibilities between mitochondrial and nuclear genomes, which could lead to the exclusion of T. macedonicus mtDNA in natural populations, most likely act at later stages of development or subsequent hybrid generations. Results from this study add to the growing knowledge of Triturus hybrid biology and ecology, which is the baseline for conservation programs necessary to protect these highly endangered amphibians.

Abstract

Two large-bodied newt species, Triturus ivanbureschi and T. macedonicus, hybridize in nature across the Balkan Peninsula. Consequences of hybridization upon secondary contact of two species include species displacement and asymmetrical introgression of T. ivanbureschi mtDNA. We set an experimental reciprocal cross of parental species and obtained two genotypes of F₁ hybrids (with T. ivanbureschi or T. macedonicus mtDNA). When hybrids attained sexual maturity, they were engaged in mutual crossings and backcrossing with parental species. We followed reproductive traits over two successive years. Our main aim was to explore the reproductive success of F₁ females carrying different parental mtDNA. Additionally, we tested for differences in reproductive success within female genotypes depending on the crossing with various male genotypes (hybrids or parental species). Both female genotypes had similar oviposition periods, number of laid eggs and hatched larvae but different body and egg sizes. Overall reproductive success (percentage of egg-laying females and viability of embryos) was similar for both genotypes. The type of crossing led to some differences in reproductive success within female genotypes. The obtained results suggest that processes that led to exclusion of T. macedonicus mtDNA in natural populations may be related to the survival at postembryonic stages of F₂ generation or reproductive barriers that emerged in subsequent hybrid generations.

Consequences of hybridization on life history and growth in postmetamorphic Triturus macedonicus

The Macedonian crested newt T. macedonicus forms a complex hybrid zone with the Balkan crested newt T. ivanbureschi in the central part of the Balkan Peninsula. In this study, we compared life history parameters (sex ratio, survival rate and age of sexual maturation) and growth parameters (body length and mass) of T. macedonicus and T. macedonicus × T. ivanbureschi hybrids over the first three postmetamorphic years. Survival rates were high and similar for both genotypes. Sexual dimorphism in body size (length and mass) was evident at the beginning of the third postmetamorphic year, after the first breeding year, for both T. macedonicus and hybrids. We did not find clear and consistent difference in any of the analysed traits suggesting that hybridization with T. ivabureschi does not affect postmetamorphic growth and survival of T. macedonicus.

Next-generation phylogeography of the banded newts (Ommatotriton): A phylogenetic hypothesis for three ancient species with geographically restricted interspecific gene flow and deep intraspecific genetic structure

Technological developments now make it possible to employ many markers for many individuals in a phylogeographic setting, even for taxa with large and complex genomes such as salamanders. The banded newt (genus Ommatotriton) from the Near East has been proposed to contain three species (O. nesterovi, O. ophryticus and O. vittatus) with unclear phylogenetic relationships, apparently limited interspecific gene flow and deep intraspecific geographic mtDNA structure. We use parallel tagged amplicon sequencing to obtain 177 nuclear DNA markers for 35 banded newts sampled throughout the range. We determine population structure (with Bayesian clustering and principal component analysis), interspecific gene flow (by determining the distribution of species-diagnostic alleles) and phylogenetic relationships (by maximum likelihood inference of concatenated sequence data and based on a summary-coalescent approach). We confirm that the three proposed species are genetically distinct. A sister relationship between O. nesterovi and O. ophryticus is suggested. We find evidence for introgression between O. nesterovi and O. ophryticus, but this is geographically limited. Intraspecific structuring is extensive, with the only recognized banded newt subspecies, O. vittatus cilicensis, representing the most distinct lineage below the species level. While mtDNA mostly mirrors the pattern observed in nuclear DNA, all banded newt species show mito-nuclear discordance as well.

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.

Effects of thiourea on the skull of Triturus newts during ontogeny

Background

In amphibians, thyroid hormone (TH) has a profound role in cranial development, especially in ossification of the late-appearing bones and remodeling of the skull. In the present study, we explored the influence of TH deficiency on bone ossification and resulting skull shape during the ontogeny of Triturus newt hybrid larvae obtained from interspecific crosses between T. ivanbureschi and T. macedonicus.

Methods

Larvae were treated with two concentrations of thiourea (an endocrine disruptor that chemically inhibits synthesis of TH) during the midlarval and late larval periods. Morphological differences of the cranium were assessed at the end of the midlarval period (ontogenetic stage 62) and the metamorphic stage after treatment during the late larval period.

Results

There was no difference in the ossification level and shape of the skull between the experimental groups (control and two treatment concentrations) at stage 62. During the late larval period and metamorphosis, TH deficit had a significant impact on the level of bone ossification and skull shape with no differences between the two treatment concentrations of thiourea. The most pronounced differences in bone development were: the palatopterygoid failed to disintegrate into the palatal and pterygoid portions, retardation was observed in development of the maxilla, nasal and prefrontal bones and larval organization of the vomer was retained in thiourea-treated larvae.

Conclusions

This implies that deficiency of TH caused retardation in development and arrested metamorphic cranium skeletal reorganization, which resulted in divergent cranial shape compared to the control group. Our results confirmed that skull remodeling and ossification of late-appearing bones is TH–dependent, as in other studied Urodela species. Also, our results indicate that TH plays an important role in the establishment of skull shape during the ontogeny of Triturus newts, especially during the late larval period and metamorphosis, when TH concentrations reach their maximum.

Phylogenomic resolution of the Ceratitis FARQ complex (Diptera: Tephritidae)

The Ceratitis FARQ complex (formerly FAR complex) includes four frugivorous tephritids, Ceratitis fasciventris, C. anonae, C. rosa and C. quilicii, the latter two causing important agricultural losses in Africa. Although FARQ species can be identified on the basis of subtle morphological differences, they cannot be resolved as monophyletic when trying phylogenetic tree reconstructions based on mitochondrial or nuclear gene fragments except for microsatellites. In this study, we used mitogenome and genome-wide SNPs to investigate the phylogenetic relationship within the complex as well as between all four Ceratitis subgenera. The analysis of 13 species supported the monophyly of the Ceratitis subgenera Ceratitis, Ceratalaspis, Pardalaspis, and recovered Pterandrus as paraphyletic but could not properly resolve species within the FARQ complex. Conversely, gene and species tree reconstructions based on 785,484 genome-wide SNPs could consistently resolve the FARQ taxa and provide insights into their phylogenetic relationships. Gene flow was detected by TreeMix analysis from C. quilicii to C. fasciventris, suggesting the existence of introgression events in the FARQ complex. Our results suggest that genome-wide SNPs represent a suitable tool for the molecular diagnosis of FARQ species and could possibly be used to develop rapid diagnostic methods or to trace the origins of intercepted samples.

Molecular Evolution of Antigen-Processing Genes in Salamanders: Do They Coevolve with MHC Class I Genes?

Proteins encoded by antigen-processing genes (APGs) prepare antigens for presentation by the major histocompatibility complex class I (MHC I) molecules. Coevolution between APGs and MHC I genes has been proposed as the ancestral gnathostome condition. The hypothesis predicts a single highly expressed MHC I gene and tight linkage between APGs and MHC I. In addition, APGs should evolve under positive selection, a consequence of the adaptive evolution in MHC I. The presence of multiple highly expressed MHC I genes in some teleosts, birds, and urodeles appears incompatible with the coevolution hypothesis. Here, we use urodele amphibians to test two key expectations derived from the coevolution hypothesis: 1) the linkage between APGs and MHC I was studied in Lissotriton newts and 2) the evidence for adaptive evolution in APGs was assessed using 42 urodele species comprising 21 genera from seven families. We demonstrated that five APGs (PSMB8, PSMB9, TAP1, TAP2, and TAPBP) are tightly linked (<0.5 cM) to MHC I. Although all APGs showed some codons under episodic positive selection, we did not find a pervasive signal of positive selection expected under the coevolution hypothesis. Gene duplications, putative gene losses, and divergent allelic lineages detected in some APGs demonstrate considerable evolutionary dynamics of APGs in salamanders. Overall, our results indicate that if coevolution between APGs and MHC I occurred in urodeles, it would be more complex than envisaged in the original formulation of the hypothesis.

Tracing species replacement in Iberian marbled newts

Secondary contact between closely related species can lead to the formation of hybrid zones, allowing for interspecific gene flow. Hybrid zone movement can take place if one of the species possesses a competitive advantage over the other, ultimately resulting in species replacement. Such hybrid zone displacement is predicted to leave a genomic footprint across the landscape in the form of asymmetric gene flow (or introgression) of selectively neutral alleles from the displaced to the advancing species. Hybrid zone movement has been suggested for marbled newts in the Iberian Peninsula, supported by asymmetric gene flow and a distribution relict (i.e., an enclave) of Triturus marmoratus in the range of T. pygmaeus. We developed a panel of nuclear and mitochondrial SNP markers to test for the presence of a T. marmoratus genomic footprint in the Lisbon peninsula, south of the enclave. We found no additional populations of T. marmoratus. Analysis with the software Structure showed no genetic traces of T. marmoratus in T. pygmaeus. A principal component analysis showed some variation within the local T. pygmaeus, but it is unclear if this represents introgression from T. marmoratus. The results may be explained by (a) species replacement without introgressive hybridization and (b) displacement with hybridization followed by the near-complete erosion of the footprint by purifying selection. We predict that testing for a genomic footprint north of the reported enclave would confirm that species replacement in these marbled newts occurred with hybridization.

Genetic Divergence Across Glacial Refugia Despite Interglacial Gene Flow in a Crested Newt

MtDNA-based phylogeography has illuminated the impact of the Pleistocene Ice Age on species distribution dynamics and the build-up of genetic divergence. The well-known shortcomings of mtDNA in biogeographical inference can be compensated by integrating multilocus data and species distribution modelling into phylogeography. We re-visit the phylogeography of the Italian crested newt (Triturus carnifex), a species distributed in two of Europe’s main glacial refugia, the Balkan and Italian Peninsulas. While a new 51 nuclear DNA marker dataset supports the existence of three lineages previously suggested by mtDNA (Balkan, northern Italy and southern Italy), the nuclear DNA dataset also provides improved resolution where these lineages have obtained secondary contact. We observe geographically restricted admixture at the contact between the Balkan and northern Italy gene pools and identify a potential mtDNA ghost lineage here. At the contact between the northern and southern Italy gene pools we find admixture over a broader area, as well as asymmetric mtDNA introgression. Our species distribution model is in agreement with a distribution restricted to distinct refugia during Pleistocene glacial cycles and postglacial expansion with secondary contact. Our study supports: (1) the relevance of the north-western Balkan Peninsula as a discrete glacial refugium; (2) the importance of north-eastern Italy and the northern Apennine as suture zones; and (3) the applicability of a refugia-within-refugia scenario within the Italian Peninsula.

Phylotranscriptomic evidence for pervasive ancient hybridization among Old World salamanders

Hybridization can leave genealogical signatures in an organism's genome, originating from the parental lineages and persisting over time. This potentially confounds phylogenetic inference methods that aim to represent evolution as a strictly bifurcating tree. We apply a phylotranscriptomic approach to study the evolutionary history of, and test for inter-lineage introgression in the Salamandridae, a Holarctic salamanders group of interest in studies of toxicity and aposematism, courtship behavior, and molecular evolution. Although the relationships between the 21 currently recognized salamandrid genera have been the subject of numerous molecular phylogenetic studies, some branches have remained controversial and sometimes affected by discordances between mitochondrial vs. nuclear trees. To resolve the phylogeny of this family, and understand the source of mito-nuclear discordance, we generated new transcriptomic (RNAseq) data for 20 salamandrids and used these along with published data, including 28 mitochondrial genomes, to obtain a comprehensive nuclear and mitochondrial perspective on salamandrid evolution. Our final phylotranscriptomic data set included 5455 gene alignments for 40 species representing 17 of the 21 salamandrid genera. Using concatenation and species-tree phylogenetic methods, we find (1) Salamandrina sister to the clade of the "True Salamanders" (consisting of Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), (2) Ichthyosaura sister to the Near Eastern genera Neurergus and Ommatotriton, (3) Triturus sister to Lissotriton, and (4) Cynops paraphyletic with respect to Paramesotriton and Pachytriton. Combining introgression tests and phylogenetic networks, we find evidence for introgression among taxa within the clades of "Modern Asian Newts" and "Modern European Newts". However, we could not unambiguously identify the number, position, and direction of introgressive events. Combining evidence from nuclear gene analysis with the observed mito-nuclear phylogenetic discordances, we hypothesize a scenario with hybridization and mitochondrial capture among ancestral lineages of (1) Lissotriton into Ichthyosaura and (2) Triturus into Calotriton, plus introgression of nuclear genes from Triturus into Lissotriton. Furthermore, both mitochondrial capture and nuclear introgression may have occurred among lineages assigned to Cynops. More comprehensive genomic data will, in the future, allow testing this against alternative scenarios involving hybridization with other, extinct lineages of newts.

Phylogenomics reveals novel relationships among Neotropical crocodiles (Crocodylus spp.)

Extant species in the order Crocodylia are remnants of an ancient lineage of large-bodied archosaur reptiles. Despite decades of systematic studies, phylogenetic relationships among members of the genus Crocodylus (true crocodiles) in the Neotropics are poorly understood. Here we estimated phylogenomic relationships among the four extant Crocodylus species in the Americas. Species-tree reconstructions using genotypic data from 17,538 SNPs collected for 33 individuals spanning six Crocodylus species (four ingroup and two outgroup) revealed novel relationships for all Neotropical species. For the first time, C. acutus, the American crocodile, was recovered as monophyletic when individuals from Antillean and continental populations were analyzed together. Our results also contradict previous inferences based on mitochondrial DNA data and a limited number of nuclear markers by robustly grouping Morelet's crocodile (C. moreletii) as the sister species to C. acutus, suggesting a novel phylogeographic hypothesis for the group. The present study punctuates the importance of using nuclear genome-wide information and representative sampling for resolving phylogenetic relationships, especially in broadly distributed species and those with complex evolutionary histories.

Extensive cytonuclear discordance in a crested newt from the Balkan Peninsula glacial refugium

Integration of multilocus data and species distribution modelling into phylogeography allows mitochondrial DNA (mtDNA)-based scenarios to be fine-tuned. We address the question of whether extensive mtDNA substructuring in the crested newt Triturus macedonicus from the Balkan Peninsula is matched in the nuclear genome. We determine the intraspecific population structure based on 52 nuclear DNA markers and project a species distribution model on climate layers for the Last Glacial Maximum. We show that T. macedonicus accumulated nuclear DNA population structure in an area predicted to have been climatically stable during the Pleistocene, with four nuclear DNA groups in the western part of the species range. The distribution of these nuclear DNA groups shows little agreement with that of mtDNA structuring, which shows three highly distinct species-specific clades and a fourth one introgressed from another crested newt species. This cytonuclear discordance conveys that historical biogeographical scenarios based on mtDNA exclusively should be interpreted with caution. Our findings further highlight the important role the Balkan Peninsula has played in the evolution and preservation of European biodiversity.

The transcriptome of the newt Cynops orientalis provides new insights into evolution and function of sexual gene networks in sarcopterygians

Amphibians evolved in the Devonian period about 400 Mya and represent a transition step in tetrapod evolution. Among amphibians, high-throughput sequencing data are very limited for Caudata, due to their largest genome sizes among terrestrial vertebrates. In this paper we present the transcriptome from the fire bellied newt Cynops orientalis. Data here presented display a high level of completeness, comparable to the fully sequenced genomes available from other amphibians. Moreover, this work focused on genes involved in gametogenesis and sexual development. Surprisingly, the gsdf gene was identified for the first time in a tetrapod species, so far known only from bony fish and basal sarcopterygians. Our analysis failed to isolate fgf24 and foxl3, supporting the possible loss of both genes in the common ancestor of Rhipidistians. In Cynops, the expression analysis of genes described to be sex-related in vertebrates singled out an expected functional role for some genes, while others displayed an unforeseen behavior, confirming the high variability of the sex-related pathway in vertebrates.

Integrating hybrid zone analyses in species delimitation: lessons from two anuran radiations of the Western Mediterranean

Molecular ecologists often rely on phylogenetic evidence for assessing the species-level systematics of newly discovered lineages. Alternatively, the extent of introgression at phylogeographic transitions can provide a more direct test to assign candidate taxa into subspecies or species categories. Here, we compared phylogenetic versus hybrid zone approaches of species delimitation in two groups of frogs from the Western Mediterranean region (Discoglossus and Pelodytes), by using genomic data (ddRAD). In both genera, coalescent analyses recovered almost all nominal taxa as "species". However, the least-diverged pairs D. g. galganoi/jeanneae and P. punctatus/hespericus admix over hundreds of kilometers, suggesting that they have not yet developed strong reproductive isolation and should be treated as conspecifics. In contrast, the comparatively older D. scovazzi/pictus and P. atlanticus/ibericus form narrow contact zones, consistent with species distinctiveness. Due to their complementarity, we recommend taxonomists to combine phylogenomics with hybrid zone analyses to scale the gray zone of speciation, i.e., the evolutionary window separating widely admixing lineages versus nascent reproductively isolated species. The radically different transitions documented here conform to the view that genetic incompatibilities accumulating with divergence generate a weak barrier to gene flow for long periods of time, until their effects multiply and the speciation process then advances rapidly. Given the variability of the gray zone among taxonomic groups, at least from our current abilities to measure it, we recommend to customize divergence thresholds within radiations to categorize lineages for which no direct test of speciation is possible.

Species-specific habitat preferences do not shape the structure of a crested newt hybrid zone (Triturus cristatus x T. carnifex)

Reproductive isolation barriers maintain the integrity of species by preventing interspecific gene flow. They involve temporal, habitat or behavioral isolation acting before fertilization, and postzygotic isolation manifested as hybrid mortality or sterility. One of the approaches of how to study reproductive isolation barriers is through the analysis of hybrid zones. In this paper, we describe the structure of a hybrid zone between two crested newt species (Triturus cristatus and T. carnifex) in the southern part of the Czech Republic using morphological, microsatellite, and mitochondrial (mtDNA) markers. Specifically, we tested the hypothesis that the structure of the hybrid zone is maintained by species-specific habitat preferences. Comparing the genetic structure of populations with geographical and ecological parameters, we found that the hybrid zone was structured primarily geographically, with T. cristatus-like populations occurring in the northeast and T. carnifex-like populations in the southwest. Despite T. cristatus tending to occur in deeper ponds and T. carnifex on localities with more shading, the effect of both ecological parameters on the structure of the zone was minimal. Next, we corroborated that T. carnifex individuals and some hybrids possess mtDNA of T. dobrogicus, whose nuclear background was not detected in the studied hybrid zone. Hybridization between T. carnifex and T. dobrogicus (resulting in unidirectional mtDNA introgression) had to predate subsequent formation of the hybrid zone between T. cristatus and T. carnifex. Populations of crested newts in the southern part of the Czech Republic thus represent a genetic mosaic of nuclear and mitochondrial genomes of three species.

MHC structuring and divergent allele advantage in a urodele amphibian: a hierarchical multi-scale approach

Proteins encoded by extraordinarily polymorphic major histocompatibility complex (MHC) genes are involved in the adaptive immune response. Balancing selection is believed to maintain MHC polymorphism in the long term, although neutral processes also play a role in shaping MHC diversity. However, the relative contribution of these processes is poorly understood. Here we characterized MHC class II variation of a low-dispersal, pond-breeding newt (Triturus carnifex) over a restricted, geographically structured area. We aimed to (1) evaluate the contribution of selection and neutral processes to shaping MHC diversity at two geographic scales, and (2) test for signatures of divergent allele advantage (DAA), which is a potentially important mechanism of balancing selection. The dominant role of selection in shaping MHC variation was suggested by the lack of correlation between MHC and neutral (microsatellite) variation. Although most variation occurred within populations for both types of markers, they differed in the extent of structuring at the two spatial scales. MHC structuring was more pronounced at local scales, suggesting the role of local selection, while structuring was not detectable at a larger scale, possibly due to the effect of balancing selection. Microsatellites showed the opposite pattern. As expected under DAA, the observed genotypes combined more sequence diversity than expected under a random association of alleles. Thus, DAA may contribute to maintaining MHC polymorphism, which is ancient, as supported by signatures of historical positive selection and trans-species polymorphism. Our results point to the importance of a multi-scale approach in studying MHC variation, especially in low-dispersal taxa, which are genetically structured at fine spatial scales.

Phylogeography of a cryptic speciation continuum in Eurasian spadefoot toads (Pelobates)

Cryptic phylogeographic diversifications provide unique models to examine the role of phylogenetic divergence on the evolution of reproductive isolation, without extrinsic factors such as ecological and behavioural differentiation. Yet, to date very few comparative studies have been attempted within such radiations. Here, we characterize a new speciation continuum in a group of widespread Eurasian amphibians, the Pelobates spadefoot toads, by conducting multilocus (restriction site associated DNA sequencing and mitochondrial DNA) phylogenetic, phylogeographic and hybrid zone analyses. Within the P. syriacus complex, we discovered species-level cryptic divergences (>5 million years ago [My]) between populations distributed in the Near-East (hereafter P. syriacus sensu stricto [s.s.]) and southeastern Europe (hereafter P. balcanicus), each featuring deep intraspecific lineages. Altogether, we could scale hybridizability to divergence time along six different stages, spanning from sympatry without gene flow (P. fuscus and P. balcanicus, >10 My), parapatry with highly restricted hybridization (P. balcanicus and P. syriacus s.s., >5 My), narrow hybrid zones (~15 km) consistent with partial reproductive isolation (P. fuscus and P. vespertinus, ~3 My), to extensive admixture between Pleistocene and refugial lineages (≤2 My). This full spectrum empirically supports a gradual build up of reproductive barriers through time, reversible up until a threshold that we estimate at ~3 My. Hence, cryptic phylogeographic lineages may fade away or become reproductively isolated species simply depending on the time they persist in allopatry, and without definite ecomorphological divergence.

Testing the evolutionary constraints of metamorphosis: The ontogeny of head shape in Triturus newts

In vertebrates with complex, biphasic, life cycles, larvae have a distinct morphology and ecological preferences compared to metamorphosed juveniles and adults. In amphibians, abrupt and rapid metamorphic changes transform aquatic larvae to terrestrial juveniles. The main aim of this study is to test whether, relative to larval stages, metamorphosis (1) resets the pattern of variation between ontogenetic stages and species, (2) constrains intraspecific morphological variability, and (3) similar to the "hour-glass" model reduces morphological disparity. We explore postembryonic ontogenetic trajectories of head shape (from hatching to completed metamorphosis) of two well-defined, morphologically distinct Triturus newts species and their F1 hybrids. Variation in head shape is quantified and compared on two levels: dynamic (across ontogenetic stages) and static (at a particular stage). Our results show that the ontogenetic trajectories diverge early during development and continue to diverge throughout larval stages and metamorphosis. The high within-group variance and the largest disparity level (between-group variance) characterize the metamorphosed stage. Hence, our results indicate that metamorphosis does not canalize head shape variation generated during larval development and that metamorphosed phenotype is not more constrained relative to larval ones. Therefore, metamorphosis cannot be regarded as a developmental constraint, at least not for salamander head shape.