Transcriptomics of differences in thermal plasticity associated with selection for an exaggerated male sexual trait.

The information about the magnitude of differences in thermal plasticity both between and within populations, as well as identification of the underlying molecular mechanisms are key to understanding the evolution of thermal plasticity. In particular, genes underlying variation in the physiological response to temperature can provide raw material for selection acting on plastic traits. Using RNAseq, we investigate the transcriptional response to temperature in males and females from bulb mite populations selected for the increased frequency of one of two discrete male morphs (fighter- and scrambler-selected populations) that differ in relative fitness depending on temperature. We show that different mechanisms underlie the divergence in thermal response between fighter- and scrambler-selected populations at decreased vs. increased temperature. Temperature decrease to 18 °C was associated with higher transcriptomic plasticity of males with more elaborate armaments, as indicated by a significant selection-by-temperature interaction effect on the expression of 40 genes, 38 of which were upregulated in fighter-selected populations in response to temperature decrease. In response to 28 °C, no selection-by-temperature interaction in gene expression was detected. Hence, differences in phenotypic response to temperature increase likely depended on genes associated with their distinct morph-specific thermal tolerance. Selection of males also drove gene expression patterns in females. These patterns could be associated with temperature-dependent fitness differences between females from fighter- vs. scrambler-selected populations reported in previous studies. Our study shows that selection for divergent male sexually selected morphologies and behaviors has a potential to drive divergence in metabolic pathways underlying plastic response to temperature in both sexes.

Lizard host abundances and climatic factors explain phylogenetic diversity and prevalence of blood parasites on an oceanic island.

Host abundance might favour the maintenance of a high phylogenetic diversity of some parasites via rapid transmission rates. Blood parasites of insular lizards represent a good model to test this hypothesis because these parasites can be particularly prevalent in islands and host lizards highly abundant. We applied deep amplicon sequencing and analysed environmental predictors of blood parasite prevalence and phylogenetic diversity in the endemic lizard Gallotia galloti across 24 localities on Tenerife, an island in the Canary archipelago that has experienced increasing warming and drought in recent years. Parasite prevalence assessed by microscopy was over 94%, and a higher proportion of infected lizards was found in warmer and drier locations. A total of 33 different 18s rRNA parasite haplotypes were identified, and the phylogenetic analyses indicated that they belong to two genera of Adeleorina (Apicomplexa: Coccidia), with Karyolysus as the dominant genus. The most important predictor of between-locality variation in parasite phylogenetic diversity was the abundance of lizard hosts. We conclude that a combination of climatic and host demographic factors associated with an insular syndrome may be favouring a rapid transmission of blood parasites among lizards on Tenerife, which may favour the maintenance of a high phylogenetic diversity of parasites.

Assessing fine-scale pondscape connectivity with amphibian eyes: An integrative approach using genomic and capture-mark-recapture data.

In the face of habitat loss, preserving functional connectivity is essential to maintain genetic diversity and the demographic dynamics required for the viability of biotic communities. This requires knowledge of the dispersal behaviour of target species, which can be modelled as kernels, or probability density functions of dispersal distances at increasing geographic distances. We present an integrative approach to investigate the relationships between genetic connectivity and demographic parameters in organisms with low vagility focusing on five syntopic pond-breeding amphibians. We genotyped 1056 individuals of two anuran and three urodele species (1732-3913 SNPs per species) from populations located in a landscape comprising 64 ponds to characterize fine-scale genetic structure in a comparative framework, and combined these genetic data with information obtained in a previous 2-year capture-mark-recapture (CMR) study. Specifically, we contrasted graphs reconstructed from genomic data with connectivity graphs based on dispersal kernels and demographic information obtained from CMR data from previous studies, and assessed the effects of population size, population density, geographical distances, inverse movement probabilities and the presence of habitat patches potentially functioning as stepping stones on genetic differentiation. Our results show a significant effect of local population sizes on patterns of genetic differentiation at small spatial scales. In addition, movement records and cluster-derived kernels provide robust inferences on most likely dispersal paths that are consistent with genomic inferences on genetic connectivity. The integration of genetic and CMR data holds great potential for understanding genetic connectivity at spatial scales relevant to individual organisms, with applications for the implementation of management actions at the landscape level.

Ancient duplication, coevolution, and selection at the MHC class IIA and IIB genes of birds.

Introduction: The Major Histocompatibility Complex (MHC) of vertebrates is a dynamically evolving multigene family primarily responsible for recognizing non-self peptide antigens and triggering a pathogen-specific adaptive immune response. In birds, the MHC was previously thought to evolve via concerted evolution with high degree of gene homogenization and the rapid loss of orthology. However, the discovery of two ancient avian MHC-IIB gene lineages (DAB1 and DAB2) originating before the radiation of extant birds indicated that despite the action of concerted evolution, orthology may be detectable for long evolutionary periods. Methods: Here, we take advantage of rapidly accumulating digital genomic resources to search for the signal of an ancient duplication at the avian MHC-IIA genes, as well as to compare phylogenetic distribution and selection between MHC-IIA and IIB gene lineages. Results: The analysis of MHC sequences from over 230 species representing ca. 70 bird families revealed the presence of two ancient MHC-IIA gene lineages (DAA1 and DAA2) and showed that their phylogenetic distribution matches exactly the distribution of DAB1 and DAB2 lineages, suggesting tight coevolution. The early post-duplication divergence of DAA1 and DAA2 was driven by positive selection fixing radical amino acid differences within the membrane-proximal domain and, most probably, being functionally related to the interactions between α2 and ß2 chains of the MHC-II heterodimer. We detected no evidence for an overall (gene-wide) relaxation or intensification of selection at either DAA1/DAB1 or DAA2/DAB2, but codon-specific differences in selection signature were found at the peptide-binding sites between the two gene lineages, perhaps implying specialization to different pathogen regimes. Discussion: Our results suggest that specific pairing of MHC-II α and ß chains may have an adaptive significance, a conclusion that advances knowledge on the macroevolution of the avian MHC-II and opens exciting novel directions for future research.

Genomic analysis reveals complex population structure within the smooth newt, Lissotriton vulgaris, in Central Europe.

Species with wide-range distributions usually display high genetic variation. This variation can be partly explained by historical lineages that were temporally isolated from each other and are back into secondary reproductive contact, and partly by local adaptations. The smooth newt (Lissotriton vulgaris) is one of the most widely distributed amphibians species across Eurasia and forms a species complex with a partially overlapping distribution and morphology. In the present study, we explored the population genomic structure of smooth newt lineages in the Carpathian Basin (CB) relying on single-nucleotide polymorphisms. Our dataset included new and previously published data to study the secondary contact zone between lineages in the CB and also tested for the barrier effect of rivers to gene flow between these lineages. We confirmed the presence of the South L. v. vulgaris Lineage distributed in Transdanubia and we provided new distribution records of L. v. ampelensis inhabiting the eastern territories of the CB. High genetic diversity of smooth newts was observed, especially in the North Hungarian Mountains and at the interfluves of the main rivers in the South with four distinct lineages of L. v. vulgaris and one lineage of L. v. ampelensis showing a low level of admixture with the spatially closest L. v. vulgaris lineage. Moreover, admixture detected at the interfluve of the main rivers (i.e. Danube and Tisza) suggested a secondary contact zone in the area. Finally, we found that the river Danube has a very weak effect on population divergence, while the river Tisza is a geographical barrier limiting gene flow between smooth newt lineages. As the range boundaries of L. v. vulgaris and L. v. ampelensis in the CB coincide with the river Tisza, our study underpins the influence of rivers in lineage diversification.

MHC architecture in amphibians - ancestral reconstruction, gene rearrangements and duplication patterns.

The hypervariable major histocompatibility complex (MHC) is a crucial component of vertebrate adaptive immunity, but large-scale studies on MHC macroevolution in non-model vertebrates have long been constrained by methodological limitations. Here, we used rapidly accumulating genomic data to reconstruct macroevolution of the MHC region in amphibians. We retrieved contigs containing the MHC region from genome assemblies of 32 amphibian species and examined major structural rearrangements, duplication patterns and gene structure across the amphibian phylogeny. Based on the few available caecilian and urodele genomes we showed that the structure of ancestral MHC region in amphibians was probably relatively simple and compact, with a close physical linkage between MHC-I and MHC-II regions. This ancestral MHC architecture was generally conserved in anurans, although the evolution of class I subregion proceeded towards more extensive duplication and rapid expansion of gene copy number, providing evidence for dynamic evolutionary trajectories. Although in anurans we recorded tandems of duplicated MHC-I genes outside the core subregion, our phylogenetic analyses of MHC-I sequences provided little support for an expansion of nonclassical MHC-Ib genes across amphibian families. Finally, we found that intronic regions of amphibian classical MHC genes were much longer when compared to other tetrapod lineages (birds and mammals), which could partly be driven by the expansion of genome size. Our study reveals novel evolutionary patterns of the MHC region in amphibians and provides a comprehensive framework for further studies on the MHC macroevolution across vertebrates.

Widespread introgression of MHC genes in Iberian Podarcis lizards.

Major histocompatibility complex (MHC) genes are crucial for the adaptive immune response of jawed vertebrates. Their variation, reaching extreme levels, is driven mainly by an arms race between hosts and pathogens. One hypothesised mechanism contributing to MHC polymorphism is adaptive introgression, the exchange of genetic variants between hybridising species favoured by selection, yet its effect on MHC variation is poorly understood. Detection of adaptive MHC introgression, though challenging, may be facilitated by the analysis of species complexes forming multiple hybrid zones. Here, we investigated MHC introgression in six hybrid zones formed by seven species of Podarcis lizards inhabiting the Iberian Peninsula. To differentiate adaptive introgression from neutral introgression, we compared the patterns of gene exchange in MHC and genome-wide markers. We found elevated sharing of MHC alleles in the proximity of contact beyond the areas of detectable genome-wide admixture in most hybrid zones and, in half of them, asymmetric MHC exchange. In general, the elevated MHC allele sharing between species pairs with abutting ranges compared to geographically isolated species pairs also supports the prevalence of introgression. Collectively, our results demonstrate widespread MHC introgression in the Iberian Podarcis complex and suggest its adaptiveness. Contrary to previous results from Triturus newts, we did not observe differences in the rate of introgression between MHC classes. Our work adds support to the emerging view of adaptive introgression as a key mechanism shaping MHC diversity. It also raises questions about the effect of elevated MHC variation and factors leading to the asymmetry of adaptive introgression.

The effect of temperature and invasive alien predator on genetic and phenotypic variation in the damselfly Ischnura elegans: cross-latitude comparison.

BACKGROUND: Understanding and predicting how organisms respond to human-caused environmental changes has become a major concern in conservation biology. Here, we linked gene expression and phenotypic data to identify candidate genes underlying existing phenotypic trait differentiation under individual and combined environmental variables. For this purpose, we used the damselfly Ischnura elegans. Egg clutches from replicated high- (southern Sweden) and central-latitude (southern Poland) populations facing different degrees of seasonal time constraints were collected. Damselfly larvae were exposed to experimental treatments: current and mild warming temperatures crossed with the presence or absence of an invasive alien predator cue released by the spiny-cheek crayfish, Faxonius limosus, which is only present in Poland to date. We measured the following traits: larval development time, body size, mass and growth rate, and used the larvae for gene expression analysis by RNA-seq. Data were analysed using a multivariate approach. RESULTS: We showed latitudinal differences in coping with mild warming and predator cues. When exposed to an increased temperature and a predator cue, central-latitude individuals had the shortest development and the fastest growth compared to high-latitude individuals. There was a general effect of predator cues regarding mass and growth rate reduction independent of latitude. Transcriptome analysis revealed that metabolic pathways related to larval anatomy and development tended to be upregulated in response to mild warming but only in fast-growing central-latitude individuals. Metabolic pathways linked to oxidative stress tended to be downregulated in response to a predator cue, especially in central-latitude individuals. CONCLUSION: Different phenotypic and transcriptomic responses to environmental factors might be attributed to the variability in I. elegans life history strategies between the two latitudes caused by seasonal time constraints and to its coexistence with the invasive alien predator in nature. By providing insights into how organisms may respond to future anthropogenic changes, our results may be of particular interest in conservation biology.

Different waves of postglacial recolonisation and genomic structure of bank vole populations in NE Poland.

Previous studies indicated that in some species phylogeographic patterns obtained in the analysis of nuclear and mitochondrial DNA (mtDNA) markers can be different. Such mitonuclear discordance can have important evolutionary and ecological consequences. In the present study, we aimed to check whether there was any discordance between mtDNA and nuclear DNA in the bank vole population in the contact zone of its two mtDNA lineages. We analysed the population genetic structure of bank voles using genome-wide genetic data (SNPs) and diversity of sequenced heart transcriptomes obtained from selected individuals from three populations inhabiting areas outside the contact zone. The SNP genetic structure of the populations confirmed the presence of at least two genetic clusters, and such division was concordant with the patterns obtained in the analysis of other genetic markers and functional genes. However, genome-wide SNP analyses revealed the more detailed structure of the studied population, consistent with more than two bank vole recolonisation waves, as recognised previously in the study area. We did not find any significant differences between individuals representing two separate mtDNA lineages of the species in functional genes coding for protein-forming complexes, which are involved in the process of cell respiration in mitochondria. We concluded that the contemporary genetic structure of the populations and the width of the contact zone were shaped by climatic and environmental factors rather than by genetic barriers. The studied populations were likely isolated in separate Last Glacial Maximum refugia for insufficient amount of time to develop significant genetic differentiation.

No detectable changes in reproductive behaviour of Caenorhabditis elegans males after 97 generations under obligatory outcrossing.

In Caenorhabditis elegans, a species reproducing mostly via self-fertilization, numerous signatures of selfing syndrome are observed, including differences in reproductive behaviour compared to related obligatory outcrossing species. In this study we investigated the effect of nearly 100 generations of obligatory outcrossing on several characteristics of male reproductive behaviour. A genetically uniform ancestral population carrying a mutation changing the reproductive system to obligatory outcrossing was split into four independent populations. We predicted that the transition from the natural reproductive system, where males were extremely rare, to obligatory outcrossing, where males comprise 50% of the population and are necessary for reproduction, will increase the selection pressure on higher effectiveness of mating behaviour. Several characteristics of male mating behaviour during a 15 min interaction as well as copulation success were compared between the ancestral and evolved populations. No significant differences in male mating behaviour or fertilization success were detected between generations 1 and 97 of obligatory outcrossing populations. We found, however, that longer contact with females increased chances of successful copulation, although this effect did not differ between populations. We conclude that either selection acting on male mating behaviour has not been strong enough, or mutational input of new adaptive variants has not been sufficient to cause noticeable behavioural differences after 97 generations of evolution starting from genetically uniform population.

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.

Landscape genetics reveals contrasting patterns of connectivity in two newt species (Lissotriton montandoni and L. vulgaris).

Ecologically distinct species may respond to landscape changes in different ways. In addition to basic ecological data, the extent of the geographic range has been successfully used as an indicator of species sensitivity to anthropogenic landscapes, with widespread species usually found to be less sensitive compared to range-restricted species. In this study, we investigate connectivity patterns of two closely related but ecologically distinct newt species - the range-restricted, Lissotriton montandoni and the widespread, L. vulgaris - using genomic data, a highly replicated setting (six geographic regions per species), and tools from landscape genetics. Our results show the importance of forest for connectivity in both species, but at the same time suggest differential use of forested habitat, with L. montandoni and L. vulgaris showing the highest connectivity at forest-core and forest-edges, respectively. Anthropogenic landscapes (i.e., higher crop- or urban-cover) increased resistance in both species, but the effect was one to three orders of magnitude stronger in L. montandoni than in L. vulgaris. This result is consistent with a view of L. vulgaris as an ecological generalist. Even so, currently, the negative impact of anthropogenic landscapes is mainly seen in connectivity among L. vulgaris populations, which show significantly stronger isolation and lower effective sizes relative to L. montandoni. Overall, this study emphasizes how habitat destruction is compromising genetic connectivity not only in endemic, range-restricted species of conservation concern but also in widespread generalist species, despite their comparatively lower sensitivity to anthropogenic landscape changes.

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.

Genomic evidence that a sexually selected trait captures genome-wide variation and facilitates the purging of genetic load.

The evolution of costly traits such as deer antlers and peacock trains, which drove the formation of Darwinian sexual selection theory, has been suggested to both reflect and affect patterns of genetic variance across the genome, but direct tests are missing. Here, we used an evolve and resequence approach to reveal patterns of genome-wide diversity associated with the expression of a sexually selected weapon that is dimorphic among males of the bulb mite, Rhizoglyphus robini. Populations selected for the weapon showed reduced genome-wide diversity compared to populations selected against the weapon, particularly in terms of the number of segregating non-synonymous positions, indicating enhanced purifying selection. This increased purifying selection reduced inbreeding depression, but outbred female fitness did not improve, possibly because any benefits were offset by increased sexual antagonism. Most single nucleotide polymorphisms (SNPs) that consistently diverged in response to selection were initially rare and overrepresented in exons, and enriched in regions under balancing or relaxed selection, suggesting they are probably moderately deleterious variants. These diverged SNPs were scattered across the genome, further demonstrating that selection for or against the weapon and the associated changes to the mating system can both capture and influence genome-wide variation.

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.

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.

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.