Unexpected hybridization patterns in Near Eastern terrapins (Mauremys caspica,M. rivulata) indicate ancient gene flow across the Fertile Crescent

Phylogeography of the endangered "eyed" turtles (genus Sacalia) and the discovery of a lineage derived from natural interspecific hybridization

The herpetofauna of the Indomalayan bioregion of Asia suffers from severe habitat loss, unsustainable harvesting, and lack of research and conservation. Here, we investigated the range-wide phylogeography of the endangered "eyed" turtles (genus Sacalia, including the Beale's Eyed Turtle S. bealei and the Four-eyed Turtle S. quadriocellata) and discovered a natural interspecific hybrid turtle population in China. Based on phylogeny of the mitochondrial Cytochrome b gene of 101 samples in this study and public data, three major clades and six subclades were identified: S. bealei (SBE) in eastern-southern China, east S. quadriocellata in South China (northern east [SQUen] and southern east [SQUes] subclades), and west S. quadriocellata mainly in Vietnam (northern west [SQUwn], central west [SQUwc], and southern west [SQUws] subclades). We sequenced 16 nuclear DNA loci of 87 samples from SBE, SQUen, SQUes, and SQUwn subclades. Population genetic clustering analysis suggested a structure similar to the mitochondrial phylogeny, where most samples were classified into four genetic clusters corresponding to the four mtDNA subclades. However, a proportion of samples carrying SQUen mtDNA haplotypes formed an additional distinct cluster SHY. Those samples are found in the contact zone of the two species bearing mosaic and intermediate morphological characteristics. We detected an admixed ancestry in SHY from SBE and SQUen that conformed to an intrapopulation breeding scenario for at least hundreds of generations after the initial hybrid event, leading to a conclusion that SHY is a distinct and near-panmictic population derived from natural interspecific hybridization. In addition, SQUes (Hainan Island endemic) is of special concern due to significant isolation and low genetic diversity. We suggest that seven evolutionarily significant units should be recognized to facilitate appropriate conservation actions. These findings also highlight the urgent need for further herpetological research and conservation in this region.

Population structure and gene flow of the syntopic turtles Emys and Mauremys from coastal and inland regions of Anatolia (Turkey): results from mitochondrial and microsatellite data

Revealing the genetic basis of the existence of different species living together in different geographic regions provides clarification of this phylogeographic differentiation. In this study, we investigated the population genetics and evaluated the level of genetic variation of inland and coastal populations of Mauremys and Emys in Turkey. Tissue samples of 196 terrapins were studied which were collected from syntopic coastal (Gölbent-Söke/Aydın; M. rivulata and E. orbicularis) and inland populations (Bahçesaray/Aksaray; M. caspica and E. orbicularis). DNA was isolated using the InnuPREP DNA Mini Kit. Mitochondrial DNA sequences and allelic variation at 13 microsatellite loci for Mauremys and 12 microsatellite loci for Emys were examined.  Three haplotypes were found for Emys orbicularis (Im, Ip and Iw) collected from the coastal region and two haplotypes for Emys orbicularis (Ig and Im) collected from inland. Two haplotypes were identified for M. caspica (Cmt8 and Cmt9) and three haplotypes were identified for M. rivulata (Rmt3, Rmt24 and Rmt26). Using microsatellites and the software STRUCTURE the most probable value for K was revealed as two 2 for both species. The F_ST value between M. rivulata and M. caspica was 0.39, and between the coastal and inland populations of E. orbicularis 0.09. It can be concluded that Emys populations tend to evolve by somehow preserving the allelic richness they have and Mauremys populations continue to differentiate so that new species emerge in the evolutionary process to reach the ideal allelic structure.

How often do they do it? An in-depth analysis of the hybrid zone of two grass snake species (Natrix astreptophora and Natrix helvetica)

We examined the contact zone of two parapatric species of grass snake (Natrix astreptophora and Natrix helvetica) in southern France. To this end, we used comprehensive sampling, analysed mtDNA sequences and microsatellite loci, and built Species Distribution Models for current and past climatic conditions. The contact zone had established by the mid-Holocene during range expansions from glacial refuges in the Iberian Peninsula (N. astreptophora) and southern or western France (N. helvetica). The contact zone represents a narrow bimodal hybrid zone, with steep genetic transition from one taxon to the other and rare hybridization, supporting species status for N. astreptophora and N. helvetica. Our results suggest that the steepness of the clines is a more robust tool for species delimitation than cline width. In addition, we discovered in western France, beyond the hybrid zone, a remote population of N. helvetica with genetic signatures of hybridization with N. astreptophora, most likely the result of human-mediated long-distance dispersal. For N. helvetica, we identified a southern and a northern population cluster, connected by broad-scale gene flow in a unimodal hybrid zone running across France. This pattern either reflects genetic divergence caused by allopatry in two microrefuges and subsequent secondary contact or introgression of foreign alleles into the southern cluster.

Mitochondrial DNA sequence analysis reveals multiple Pleistocene glacial refugia for the Yellow-spotted mountain newt, Neurergus derjugini (Caudata: Salamandridae) in the mid-Zagros range in Iran and Iraq

Phylogeography is often used to investigate the effects of glacial cycles on current genetic structure of various plant and animal species. This approach can also identify the number and location of glacial refugia as well as the recolonization routes from those refugia to the current locations. To identify the location of glacial refugia of the Yellow‐spotted mountain newt, Neurergus derjugini, we employed phylogeography patterns and genetic variability of this species by analyzing partial ND4 sequences (867 bp) of 67 specimens from 15 sampling localities from the whole species range in Iran and Iraq. Phylogenetic trees concordant with haplotype networks showed a clear genetic structure among populations as three groups corresponding to the populations in the north, center, and south. Evolutionary ages of clades north and south ranging from 0.15 to 0.17 Myr, while the oldest clade is the central clade, corresponding to 0.32 Myr. Bayesian skyline plots of population size change through time show a relatively slight increase until about 25 kyr (around the last glacial maximum) and a decline of population size about 2.5 kyr. The presence of geographically structured clades in north, center, and south sections of the species range signifies the disjunct populations that have emerged in three different refugium. This study illustrates the importance of the effect of previous glacial cycles in shaping the genetic structure of mountain species in the Zagros range. These areas are important in terms of long‐term species persistence and therefore valuable areas for conservation of biodiversity.

Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves, Passeriformes, Passer domesticus × Passer hispaniolensis)

Under different environmental conditions, hybridization between the same species might result in different patterns of genetic admixture. Particularly, species pairs with large distribution ranges and long evolutionary history may have experienced several independent hybridization events over time in different zones of overlap. In birds, the diverse hybrid populations of the house sparrow (Passer domesticus) and the Spanish sparrow (Passer hispaniolensis) provide a striking example. Throughout their range of sympatry, these two species do not regularly interbreed; however, a stabilized hybrid form (Passer italiae) exists on the Italian Peninsula and on several Mediterranean islands. The spatial distribution pattern on the Eurasian continent strongly contrasts the situation in North Africa, where house sparrows and Spanish sparrows occur in close vicinity of phenotypically intermediate populations across a broad mosaic hybrid zone. In this study, we investigate patterns of divergence and admixture among the two parental species, stabilized and nonstabilized hybrid populations in Italy and Algeria based on a mitochondrial marker, a sex chromosomal marker, and 12 microsatellite loci. In Algeria, despite strong spatial and temporal separation of urban early-breeding house sparrows and hybrids and rural late-breeding Spanish sparrows, we found strong genetic admixture of mitochondrial and nuclear markers across all study populations and phenotypes. That pattern of admixture in the North African hybrid zone is strikingly different from i) the Iberian area of sympatry where we observed only weak asymmetrical introgression of Spanish sparrow nuclear alleles into local house sparrow populations and ii) the very homogenous Italian sparrow population where the mitogenome of one parent (P. domesticus) and the Z-chromosomal marker of the other parent (P. hispaniolensis) are fixed. The North African sparrow hybrids provide a further example of enhanced hybridization along with recent urbanization and anthropogenic land-use changes in a mosaic landscape.

Hybridization and polyploidy in the weeping lizard Liolaemus chiliensis (Squamata: Liolaemidae)

In reptiles, polyploidy is an unusual phenomenon that can originate from interspecific hybridization. In Chile, the lizard Liolaemus chiliensis provides a unique model with which to study the origin of polyploidy because it occurs in populations of diploid individuals and a few populations that also contain triploid and diploid–triploid mosaic lizards. To test whether L. chiliensis polyploids were hybrids between lineages within the species, we analysed the karyotype, the mitochondrial gene Cytb, seven microsatellite loci, and the linear and geometric morphometry of individuals throughout the species distribution. We found polyploidy at different localities throughout the distribution. Hybrids were detected with the microsatellite loci and morphological data. A clear relationship between hybridization and polyploidy could not be established because triploid and mosaic lizards were hybrids and purebreds. However, our results support the occurrence of both phenomena, and further research is needed to clarify how hybridization could be involved in the origin of polyploidy.

A river runs through it: tree frog genomics supports the Dead Sea Rift as a rare phylogeographical break

Phylogeographical breaks can be viewed as regional hotspots of diversity where the genetic integrity of incipient species is put to the test. We focus on an understudied species transition from the Middle East, namely the Dead Sea Rift in the Levant region, which presumably divided the tree frogs Hyla savignyi and H. felixarabica. Combining multilocus genetic analyses (mitochondrial DNA and RAD-sequencing) with ecological niche modelling, we test whether the rift effectively acts as a biogeographical barrier preventing this pair from admixing and merging. The answer is yes: despite weak signs of introgression, all parapatric populations were assigned to either species without cyto-nuclear discordance. Yet, the projected distributions under present and glacial conditions largely overlapped in the area, meaning their current parapatric ranges do not represent an ecological transition. Instead, we hypothesize that H. savignyi and H. felixarabica are maintained apart by limited opportunities for dispersal across the barren Jordan Valley, combined with advanced reproductive isolation. Therefore, the Dead Sea Rift may represent a rare phylogeographical break, and we encourage international efforts to assess its contribution to the rich biodiversity of the Middle East.