Distortion of symmetrical introgression in a hybrid zone: evidence for locus-specific selection and uni-directional range expansion

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.

The Spatial Signature of Introgression After a Biological Invasion With Hybridization

The accumulation of genome-wide molecular data has emphasized the important role of hybridization in the evolution of many organisms, which may carry introgressed genomic segments resulting from past admixture events with other taxa. Despite a number of examples of hybridization occurring during biological invasions, the resulting spatial patterns of genomic introgression remain poorly understood. Preliminary simulation studies have suggested a heterogeneous spatial level of introgression for invasive taxa after range expansion. We investigated in detail the robustness of this pattern and its persistence over time for both invasive and local organisms. Using spatially explicit simulations, we explored the spatial distribution of introgression across the area of colonization of an invasive taxon hybridizing with a local taxon. The general pattern for neutral loci supported by our results is an increasing introgression of local genes into the invasive taxon with the increase in the distance from the source of the invasion and a decreasing introgression of invasive genes into the local taxon. However, we also show there is some variation in this general trend depending on the scenario investigated. Spatial heterogeneity of introgression within a given taxon is thus an expected neutral pattern in structured populations after a biological invasion with a low to moderate amount of hybridization. We further show that this pattern is consistent with published empirical observations. Using additional simulations, we argue that the spatial pattern of Neanderthal introgression in modern humans, which has been documented to be higher in Asia than in Europe, can be explained by a model of hybridization with Neanderthals in Eurasia during the range expansion of modern humans from Africa. Our results support the view that weak hybridization during range expansion may explain spatially heterogeneous introgression patterns without the need to invoke selection.

Phylogeographic analyses point to long-term survival on the spot in micro-endemic Lycian salamanders

Lycian salamanders (genus Lyciasalamandra) constitute an exceptional case of micro-endemism of an amphibian species on the Asian Minor mainland. These viviparous salamanders are confined to karstic limestone formations along the southern Anatolian coast and some islands. We here study the genetic differentiation within and among 118 populations of all seven Lyciasalamandra species across the entire genus’ distribution. Based on circa 900 base pairs of fragments of the mitochondrial 16SrDNA and ATPase genes, we analysed the spatial haplotype distribution as well as the genetic structure and demographic history of populations. We used 253 geo-referenced populations and CHELSA climate data to infer species distribution models which we projected on climatic conditions of the Last Glacial Maximum (LGM). Within all but one species, distinct phyloclades were identified, which only in parts matched current taxonomy. Most haplotypes (78%) were private to single populations. Sometimes population genetic parameters showed contradicting results, although in several cases they indicated recent population expansion of phyloclades. Climatic suitability of localities currently inhabited by salamanders was significantly lower during the LGM compared to recent climate. All data indicated a strong degree of isolation among Lyciasalamandra populations, even within phyloclades. Given the sometimes high degree of haplotype differentiation between adjacent populations, they must have survived periods of deteriorated climates during the Quaternary on the spot. However, the alternative explanation of male biased dispersal combined with a pronounced female philopatry can only be excluded if independent nuclear data confirm this result.

Past lake shore dynamics explain present pattern of unidirectional introgression across a habitat barrier

Introgression patterns between divergent lineages are often characterized by asymmetry in the direction and among-marker variation in the extent of gene flow, and therefore inform on the mechanisms involved in differentiation and speciation. In the present study, we test the hypothesis that unidirectional introgression between two phenotypically and genetically distinct lineages of the littoral, rock-dwelling cichlid fish Tropheus moorii across a wide sandy bay is linked to observed differences in mate preferences between the two lineages. This hypothesis predicts bi-directional nuclear gene flow and was rejected by congruent patterns of introgression in mtDNA, AFLP and microsatellite markers, with admixture confined to the populations west of the bay. This pattern can be explained on the basis of habitat changes in the course of lake level fluctuations, which first facilitated the development of a symmetric admixture zone including the area corresponding to the present sand bay and then shaped asymmetry by causing local extinctions and cessation of gene flow when this area became once more inhabitable. This conforms with previous assumptions that habitat dynamics are a primary determinant of population-level evolution in Tropheus. In this respect, Tropheus may be representative of species whose preferred habitat is subject to frequent re-structuring.

Hybridization and introgression during density-dependent range expansion: European wildcats as a case study

Interbreeding between historically allopatric species with incomplete reproductive barriers may result when species expand their range. The genetic consequences of such hybridization depend critically on the dynamics of the range expansion. Hybridization models during range expansion have been developed but assume dispersal to be independent from neighboring population densities. However, organisms may disperse because they are attracted by conspecifics or because they prefer depopulated areas. Here, through spatially explicit simulations, we assess the effect of various density-dependent dispersal modes on the introgression between two species. We find huge introgression from the local species into the invasive one with all dispersal modes investigated, even when the hybridization rate is relatively low. This represents a general expectation for neutral genes even if the dispersal modes differ in colonization times and amount of introgression. Invasive individuals attracted by conspecifics need more time to colonize the whole area and are more introgressed by local genes, whereas the opposite is found for solitary individuals. We applied our approach to a recent expansion of European wildcats in the Jura Mountains and the hybridization with domestic cats. We show that the simulations explained better the observed level of introgression at nuclear, mtDNA, and Y chromosome markers, when using solitary dispersal for wildcats instead of random or gregarious dispersal, in accordance with ecological knowledge. Thus, use of density-dependent dispersal models increases the predictive power of the approach.

Genomic replacement of native Cobitis lutheri with introduced C. tetralineata through a hybrid swarm following the artificial connection of river systems

River connections via artificial canals will bring about secondary contacts between previously isolated fish species. Here, we present a genetic consequence of such a secondary contact between Cobitis fish species, C. lutheri in the Dongjin River, and C. tetralineata in the Seomjin River in Korea. The construction of water canals about 80 years ago has unidirectionally introduced C. tetralineata into the native habitat of C. lutheri, and then these species have hybridized in the main stream section of the Dongjin River. According to the divergence population genetic analyses of DNA sequence data, the two species diverged about 3.3 million years ago, which is interestingly coincident with the unprecedented paleoceanographic change that caused isolations of the paleo-river systems in northeast Asia due to sea-level changes around the late Pliocene. Multilocus genotypic data of nine microsatellites and three nuclear loci revealed an extensively admixed structure in the hybrid zone with a high proportion of various post-F1 hybrids. Surprisingly, pure native C. lutheri was absent in the hybrid zone in contrast to the 7% of pure C. tetralineata. Such a biased proportion must have resulted from the dominant influence of continually introducing C. tetralineata on the native C. lutheri which has no supply of natives from other tributaries to the hybrid zone due to numerous low-head dams. In addition, mating experiments indicated that there is no discernible reproductive isolation between them. All the results suggest that the gene pool of native C. lutheri is being rapidly replaced by that of continually introducing C. tetralineata through a hybrid swarm for the last 80 years, which will ultimately lead to the genomic extinction of natives in this hybrid zone.

Phylogeography of Sardinian cave salamanders (genus Hydromantes) is mainly determined by geomorphology

Detecting the factors that determine the interruption of gene flow between populations is key to understanding how speciation occurs. In this context, caves are an excellent system for studying processes of colonization, differentiation and speciation, since they represent discrete geographical units often with known geological histories. Here, we asked whether discontinuous calcareous areas and cave systems represent major barriers to gene flow within and among the five species of Sardinian cave salamanders (genus Hydromantes) and whether intraspecific genetic structure parallels geographic distance within and among caves. We generated mitochondrial cytochrome b gene sequences from 184 individuals representing 48 populations, and used a Bayesian phylogeographic approach to infer possible areas of cladogenesis for these species and reconstruct historical and current dispersal routes among distinct populations. Our results show deep genetic divergence within and among all Sardinian cave salamander species, which can mostly be attributed to the effects of mountains and discontinuities in major calcareous areas and cave systems acting as barriers to gene flow. While these salamander species can also occur outside caves, our results indicate that there is a very poor dispersal of these species between separate cave systems.

Out of the swamp: unidirectional hybridization with weedy species may explain the prevalence of Amaranthus tuberculatus as a weed

*Amaranthus tuberculatus represents one of the most dramatic cases of weed invasion documented in the midwestern USA. The species is infamous for evolving resistance to multiple herbicides, and predicting whether these resistances may be transferred to widespread weeds of the Amaranthus hybridus aggregate is a matter of epidemiological concern. Here, we explore the patterns of genetic exchange between Amaranthus tuberculatus and A. hybridus in an effort to understand whether allele introgression occurs throughout the genome and if fecundity penalties are associated with genetic exchange. *We evaluated 192 homoploid BC(1)s at 197 amplified fragment length polymorphism (AFLP) loci, as well as two loci associated with herbicide resistance: ALS and PPO. We also assessed the fecundity of each genotype by evaluation of seed production or pollen development. *It was discovered that genetic exchange between the species is unidirectional. Whereas A. hybridus alleles transfer with little or no penalty to A. tuberculatus, the reciprocal exchange is significantly distorted and potentially of limited evolutionary consequence. *Our previous hypothesis suggesting unidirectional introgression at ALS owing to circumstantial linkage is now modified to account for the more generalized distortion of genetic exchange observed in this study.

CONTRASTING PATTERNS OF HYBRIDIZATION IN LARGE HOUSE SPIDERS (TEGENARIA ATRICA GROUP, AGELENIDAE)

The integrity of species is not fixed and may vary geographically. Here we investigate the geographic distributions and interactions of species in the Tegenaria atrica group (Araneae: Agelenidae). Detailed mapping of T. saeva and T. gigantea in England and Wales shows them to be broadly allopatric in southern England with a tightly defined, and possibly long-standing, narrow zone of parapatry in central southern England. In the north of England (Yorkshire), by contrast, the species are broadly sympatric as a result of recent range expansions. GIS techniques are used to map the species distributions and to quantify, we believe for the first time, the intimacy of interspecific interactions. The extent and nature of hybridization in these two areas is examined through regression and multivariate analyses of morphology. We show that the relative incidence of hybridization is much greater in Yorkshire than within the parapatric zone in the south. Clear patterns of asymmetric introgression are observed in both northern and southern England, with a greater impact of T. gigantea on T. saeva than vice versa. We find no sign of morphological reproductive character displacement at the zone of parapatry that might indicate reinforcement, although we cannot exclude more subtle effects, for example via cuticular pheromones. The integrity of these two species seems to be breaking down in northern England, a process that might gain momentum as the gene pools become more similar.