Contrasting evolutionary histories of the legless lizards slow worms (Anguis) shaped by the topography of the Balkan Peninsula

Background: Genetic architecture of a species is a result of historical changes in population size and extent of distribution related to climatic and environmental factors and contemporary processes of dispersal and gene flow. Population-size and range contractions, expansions and shifts have a substantial effect on genetic diversity and intraspecific divergence, which is further shaped by gene-flow limiting barriers. The Balkans, as one of the most important sources of European biodiversity, is a region where many temperate species persisted during the Pleistocene glaciations and where high topographic heterogeneity offers suitable conditions for local adaptations of populations. In this study, we investigated the phylogeographical patterns and demographic histories of four species of semifossorial slow-worm lizards (genus Anguis) present in the Balkan Peninsula, and tested the relationship between genetic diversity and topographic heterogeneity of the inhabited ranges. Results: We inferred phylogenetic relationships, compared genetic structure and historical demography of slow worms using nucleotide sequence variation of mitochondrial DNA. Four Anguis species with mostly parapatric distributions occur in the Balkan Peninsula. They show different levels of genetic diversity. A signature of population growth was detected in all four species but with various courses in particular populations. We found a strong correlation between genetic diversity of slow-worm populations and topographic ruggedness of the ranges (mountain systems) they inhabit. Areas with more rugged terrain harbour higher genetic diversity. Conclusions: Phylogeographical pattern of the genus Anguis in the Balkans is concordant with the refugiawithin- refugia model previously proposed for both several other taxa in the region and other main European Peninsulas. While slow-worm populations from the southern refugia mostly have restricted distributions and have not dispersed much from their refugial areas, populations from the extra-Mediterranean refugia in northern parts of the Balkans have colonized vast areas of eastern, central, and western Europe. Besides climatic historical events, the heterogeneous topography of the Balkans has also played an important role in shaping genetic diversity of slow worms.BMC Evolutionary Biology (2016), 16(1): 9

Evolution of retinoic acid receptors in chordates: insights from three lamprey species, Lampetra fluviatilis, Petromyzon marinus, and Lethenteron japonicum

International audienceBackground : Retinoic acid (RA) signaling controls many developmental processes in chordates, from early axis specification to late organogenesis. The functions of RA are chiefly mediated by a subfamily of nuclear hormone receptors, the retinoic acid receptors (RARs), that act as ligand-activated transcription factors. While RARs have been extensively studied in jawed vertebrates (that is, gnathostomes) and invertebrate chordates, very little is known about the repertoire and developmental roles of RARs in cyclostomes, which are extant jawless vertebrates. Here, we present the first extensive study of cyclostome RARs focusing on three different lamprey species: the European freshwater lamprey, Lampetra fluviatilis, the sea lamprey, Petromyzon marinus, and the Japanese lamprey, Lethenteron japonicum.Results : We identified four rar paralogs (rar1, rar2, rar3, and rar4) in each of the three lamprey species, and phylogenetic analyses indicate a complex evolutionary history of lamprey rar genes including the origin of rar1 and rar4 by lineage-specific duplication after the lamprey-hagfish split. We further assessed their expression patterns during embryonic development by in situ hybridization. The results show that lamprey rar genes are generally characterized by dynamic and highly specific expression domains in different embryonic tissues. In particular, lamprey rar genes exhibit combinatorial expression domains in the anterior central nervous system (CNS) and the pharyngeal region.Conclusions : Our results indicate that the genome of lampreys encodes at least four rar genes and suggest that the lamprey rar complement arose from vertebrate-specific whole genome duplications followed by a lamprey-specific duplication event. Moreover, we describe a combinatorial code of lamprey rar expression in both anterior CNS and pharynx resulting from dynamic and highly specific expression patterns during embryonic development. This ‘RAR code’ might function in regionalization and patterning of these two tissues by differentially modulating the expression of downstream effector genes during development

Bichir external gills arise via heterochronic shift that accelerates hyoid arch development.

In most vertebrates, pharyngeal arches form in a stereotypic anterior-to-posterior progression. To gain insight into the mechanisms underlying evolutionary changes in pharyngeal arch development, here we investigate embryos and larvae of bichirs. Bichirs represent the earliest diverged living group of ray-finned fishes, and possess intriguing traits otherwise typical for lobe-finned fishes such as ventral paired lungs and larval external gills. In bichir embryos, we find that the anteroposterior way of formation of cranial segments is modified by the unique acceleration of the entire hyoid arch segment, with earlier and orchestrated development of the endodermal, mesodermal, and neural crest tissues. This major heterochronic shift in the anteroposterior developmental sequence enables early appearance of the external gills that represent key breathing organs of bichir free-living embryos and early larvae. Bichirs thus stay as unique models for understanding developmental mechanisms facilitating increased breathing capacity

Slow worm, Anguis fragilis (Reptilia: Anguidae) as a species complex: genetic structure reveals deep divergences

a b s t r a c t Phylogenetic relationships of the Western Palearctic legless lizard genus Anguis were inferred based on a fragment of mitochondrial DNA and two nuclear protein-coding loci, C-mos and PRLR. A. cephallonica from the Peloponnese was confirmed as a valid species. It is the sister taxon to a clade comprising all other evolutionary lineages, which were shown to represent three distinct species: (1) A. fragilis sensu stricto occurring in Western and Central Europe, the north-western Balkans, with possibly isolated populations in the eastern Balkans, and presumably also in western Scandinavia and Italy; (2) A. colchica distributed from the eastern Czech Republic and the Baltic region eastward to northern Iran, presumably also in eastern Scandinavia, and the north-eastern Balkans; (3) A. graeca restricted to the southern Balkans, and partially sympatric with A. cephallonica. According to the more variable mitochondrial marker, A. graeca appears to be the sister species to A. colchica, and these taxa together form a sister clade to A. fragilis, whereas the less variable nuclear markers show A. colchica to be closer to A. fragilis. The C-mos gene has not provided substantial variation within this species complex, while the PRLR gene, which was used for the first time in phylogeographic study in a reptile, distinguished all species successfully. Intra-specific differentiation of A. colchica is discussed, and subspecific status of the Caucasian and Caspian populations is proposed. The uncovered genetic differences should be taken into account in all future biogeographical, morphological and ecological studies, as well as in conservation

Cartilage diversification and modularity drove the evolution of the ancestral vertebrate head skeleton

Abstract The vertebrate head skeleton has evolved a myriad of forms since their divergence from invertebrate chordates. The connection between novel gene expression and cell types is therefore of importance in this process. The transformation of the jawed vertebrate (gnathostome) head skeleton from oral cirri to jointed jaw elements required a diversity of cartilages as well as changes in the patterning of these tissues. Although lampreys are a sister clade to gnathostomes, they display skeletal diversity with distinct gene expression and histologies, a useful model for addressing joint evolution. Specifically, the lamprey tissue known as mucocartilage has noted similarities with the jointed elements of the mandibular arch in jawed vertebrates. We thus asked whether the cells in lamprey mucocartilage and gnathostome joint tissue could be considered homologous. To do this, we characterized new genes that are involved in gnathostome joint formation and characterized the histochemical properties of lamprey skeletal types. We find that most of these genes are minimally found in mucocartilage and are likely later innovations, but we do identify new activity for gdf5/6/7b in both hyaline and mucocartilage, supporting its role as a chondrogenic regulator. Contrary to previous works, our histological assays do not find any perichondrial fibroblasts surrounding mucocartilage, suggesting that mucocartilage is non-skeletogenic tissue that is partially chondrified. Interestingly, we also identify new histochemical features of the lamprey otic capsule that diverge from normal hyaline. Paired with our new insights into lamprey mucocartilage, we propose a broader framework for skeletal evolution in which an ancestral soxD/E and gdf5/6/7 network directs mesenchyme along a spectrum of cartilage-like features

Comparative Approaches in Vertebrate Cartilage Histogenesis and Regulation: Insights from Lampreys and Hagfishes

Jawed vertebrates (gnathostomes) have been the dominant lineage of deuterostomes for nearly three hundred fifty million years. Only a few lineages of jawless vertebrates remain in comparison. Composed of lampreys and hagfishes (cyclostomes), these jawless survivors are important systems for understanding the evolution of vertebrates. One focus of cyclostome research has been head skeleton development, as its evolution has been a driver of vertebrate morphological diversification. Recent work has identified hyaline-like cartilage in the oral cirri of the invertebrate chordate amphioxus, making cyclostomes critical for understanding the stepwise acquisition of vertebrate chondroid tissues. Our knowledge of cyclostome skeletogenesis, however, has lagged behind gnathostomes due to the difficulty of manipulating lamprey and hagfish embryos. In this review, we discuss and compare the regulation and histogenesis of cyclostome and gnathostome skeletal tissues. We also survey differences in skeletal morphology that we see amongst cyclostomes, as few elements can be confidently homologized between them. A recurring theme is the heterogeneity of skeletal morphology amongst living vertebrates, despite conserved genetic regulation. Based on these comparisons, we suggest a model through which these mesenchymal connective tissues acquired distinct histologies and that histological flexibility in cartilage existed in the last common ancestor of modern vertebrates

An Fgf–Shh positive feedback loop drives growth in developing unpaired fins

The origin and diversification of appendage types is a central question in vertebrate evolution. Understanding the genetic mechanisms that underlie fin and limb development can reveal relationships between different appendages. Here we demonstrate, using chemical genetics, a mutually agonistic interaction between Fgf and Shh genes in the developing dorsal fin of the channel catfish, Ictalurus punctatus. We also find that Fgf8 and Shh orthologs are expressed in the apical ectodermal ridge and zone of polarizing activity, respectively, in the median fins of representatives from other major vertebrate lineages. These findings demonstrate the importance of this feedback loop in median fins and offer developmental evidence for a median fin-first scenario for vertebrate paired appendage origins

The biogeography of Elaphe sauromates (Pallas, 1814), with a description of a new rat snake species

Background The rat snake genus Elaphe once comprised several dozens of species distributed in temperate through tropical zones of the New and Old World. Based on molecular-genetic analyses in early 2000s, the genus was split into several separate genera, leaving only 15 Palearctic and Oriental species as its members. One of the three species also occurring in Europe is Elaphe sauromates, a robust snake from the Balkans, Anatolia, Caucasus, Ponto-Caspian steppes, and Levant that has been suspected to be composed of two or more genetically diverse populations. Here, we studied the genetic structure and morphological variation of E. sauromates, aiming to better understand its inter-population relationships and biogeography, and subsequently revise its taxonomy. Methods We reconstructed the phylogeography and analyzed the genetic structure of E. sauromates populations originating from most of its geographic range using both mitochondrial (COI, ND4) and nuclear (C-MOS, MC1R, PRLR, RAG1) DNA gene fragments. We employed Maximum likelihood and Bayesian inference methods for the phylogenetic tree reconstructions, supplemented with species delimitation methods, analysis of haplotype networks, and calculation of uncorrected p-distances. Morphological variation in 15 metric and 18 meristic characters was studied using parametric univariate tests as well as multivariate general linearized models. In total, we analyzed sequences originating from 63 specimens and morphological data from 95 specimens of E. sauromates sensu lato. Results The molecular phylogeny identified two clearly divergent sister lineages within E. sauromates, with both forming a lineage sister to E. quatuorlineata. The genetic distance between them (5.80–8.24% in mtDNA) is similar to the distances among several other species of the genus Elaphe. Both lineages are also moderately morphologically differentiated and, while none of the characters are exclusively diagnostic, their combination can be used for confident lineage identification. Here, following the criteria of genetic and evolutionary species concepts, we describe the lineage from eastern Anatolia and parts of the Lesser and Great Caucasus as a new species E. urartica sp. nov. Discussion Elaphe urartica sp. nov. represents a cryptic species whose ancestors presumably diverged from their common ancestor with E. sauromates around the Miocene-Pliocene boundary. The intraspecific genetic structure indicates that the recent diversity of both species has been predominantly shaped by Pleistocene climatic oscillations, with glacial refugia mainly located in the Balkans, Crimea, and/or Anatolia in E. sauromates and Anatolia and/or the Caucasus in E. urartica sp. nov

Morphological and genetic differentiation in the anguid lizard Pseudopus apodus supports the existence of an endemic subspecies in the Levant

The Levant represents one of the most important reptile diversity hotspots and centers of endemism in the Western Palearctic. The region harbored numerous taxa in glacial refugia during the Pleistocene climatic oscillations. Due to the hostile arid conditions in the warmer periods they were not always able to spread or come into contact with populations from more distant regions. One large and conspicuous member of the Levantine herpetofauna is the legless anguid lizard Pseudopus apodus. This species is distributed from the Balkans to Central Asia with a portion of its range running along the eastern Mediterranean coast. Mitochondrial and nuclear DNA sequences, microsatellite genotypes, and morphology show that populations in this region differ from the two named subspecies and presumably had a long independent evolutionary history during the Quaternary. Here we describe the Levantine population as a new subspecies and present biogeographic scenarios for its origin and diversification. The new subspecies is genetically highly diverse, and it forms a sister lineage to Pseudopus from the remaining parts of the range according to mtDNA. It is the largest-bodied of the three subspecies, but occupies the smallest range

Multilocus phylogeny and coalescent species delimitation in Kotschy's gecko, Mediodactylus kotschyi: Hidden diversity and cryptic species

Kotschy's Gecko, Mediodactylus kotschyi, is a small gecko native to southeastern Europe and the Levant. It displays great morphological variation with a large number of morphologically recognized subspecies. However, it has been suggested that it constitutes a species complex of several yet unrecognized species. In this study, we used multilocus sequence data (three mitochondrial and three nuclear gene fragments) to estimate the phylogenetic relationships of 174 specimens from 129 sampling localities, covering a substantial part of the distribution range of the species. Our results revealed high genetic diversity of M. kotschyi populations and contributed to our knowledge about the phylogenetic relationships and the estimation of the divergence times between them. Diversification within M. kotschyi began approximately 15 million years ago (Mya) in the Middle Miocene, whereas the diversification within most of the major clades have been occurred in the last 5 Mya. Species delimitation analysis suggests there exists five species within the complex, and we propose to tentatively recognize the following taxa as full species: M. kotschyi (mainland Balkans, most of Aegean islands, and Italy), M. orientalis (Levant, Cyprus, southern Anatolia, and south-eastern Aegean islands), M. danilewskii (Black Sea region and south-western Anatolia), M. bartoni (Crete), and M. oertzeni (southern Dodecanese Islands). This newly recognized diversity underlines the complex biogeographical history of the Eastern Mediterranean region