Combining immunological and molecular data to assess phylogenetic relations of some Greek Podarcis species

Most recent molecular studies revealed the phylogeny of Greek Podarcis species, which for years remained elusive, due to discordant data produced from various chromosomal, complement fixation and protein studies. In this report, we analyzed cellular immune responses of spleen-derived lymphocytes from six allopatric Podarcis species encountered in Greece, by assessing two-way mixed lymphocyte reaction (MLR)-induced proliferation. On the basis of stimulation indices (S.I.) as determined from cultures set up from xenogeneic splenocytes coincubated in pairs, we generated a phylogenetic tree, fully consistent with the phylogenetic relationships of Podarcis as determined by parallel analyses based on partial mitochondrial (mt) DNA sequences. Although the exact mechanisms triggering lymphocyte responses in lizard two-way xenogeneic MLR are not fully understood, our results show the potential use of cell-mediated immune responses as an additional approach to mtDNA analysis, for species delimitation within specific lizard taxa.     

Phylogenetic relationships and biogeography of Podarcis species from the Balkan Peninsula, by bayesian and maximum likelihood analyses of mitochondrial DNA sequences

Wall lizards of the genus Podarcis (Sauria, Lacertidae) comprise 17 currently recognized species in southern Europe, where they are the predominant nonavian reptile group. The taxonomy of Podarcis is complex and unstable. Based on DNA sequence data, the species of Podarcis falls into four main groups that have substantial geographic coherence (Western island group, southwestern group, Italian group, and Balkan Peninsula group). The Balkan Peninsula species are divided into two subgroups: the subgroup of P. taurica (P. taurica, P. milensis, P. gaigeae, and perhaps P. melisellensis), and the subgroup of P. erhardii (P. erhardii and P. peloponnesiaca). In the present study, the question of phylogenetic relationships among the species of Podarcis encountered in the Balkan Peninsula was addressed using partial mtDNA sequences for cytochrome b (cyt b) and 16S rRNA (16S). The data support the monophyly of Podarcis and suggest that there are three phylogenetic clades: the clade A (P. taurica, P. gaigeae, P. milensis, and P. melisellensis); the clade B (P. erhardii and P. peloponnesiaca), and the clade C (P. muralis and P. sicula). By examining intraspecific relationships it was found that extant populations of P. erhardii are paraphyletic. Furthermore, subspecies previously defined on the basis of morphological characteristics do not correspond to different molecular phylogenetic clades, suggesting that their status should be reconsidered. The distinct geographic distribution of the major clades of the phylogenetic tree and its topology suggest a spatial and temporal sequence of phylogenetic separations that coincide with some major paleogeographic separations during the geological history of the Aegean Sea. The results stress the need for a reconsideration of the evolutionary history of Balkan Podarcis species and help overcome difficulties that classical taxonomy has encountered at both the species and subspecies level.     

Phylogeography of Balkan wall lizard (Podarcis taurica) and its relatives inferred from mitochondrial DNA sequences

Wall lizards of the genus Podarcis (Sauria, Lacertidae) comprise 17 currently recognized species in southern Europe, where they are the predominant reptile group. The taxonomy of Podarcis is complex and unstable. Based on DNA sequence data the species of Podarcis falls into four main groups that have substantial geographical conherence (western island group, southwestern group, Italian group and Balkan group). The Balkan species are divided in two subgroups: the subgroup of Podarcis taurica (P. taurica, P. milensis, P. gaigeae and perhaps P. melisellensis), and the subgroup of Podarcis erhardii (P. erhardii and P. peloponnesiaca). We addressed the question of phylogenetic relations among the species of the P. taurica subgroup encountered in Greece, as they can be inferred from partial mtDNA (cyt b and 16S) sequences. Our data support the monophyly of P. taurica subgroup and suggest that P. gaigeae, P. milensis and P. melisellensis form a clade, which thereinafter connects to P. taurica. Within the previous clade, P. gaigeae is more closely related to P. milensis than to P. melisellensis. However, the specimens of P. taurica were subdivided in two different groups. The first one includes the specimens from northeastern Greece, and the other group includes the specimens from the rest of continental Greece and Ionian islands. Because the molecular clock of the cyt b and 16 rRNA genes was not rejected in our model test, it is possible to estimate times of speciation events. Based on the splitting of the island of Crete from Peloponnisos [c. 5 million years ago (Ma)], the evolutionary rate for the cyt b is 1.55% per million years (Myr) and for the 16S rRNA is 0.46% per Myr. These results suggest that the evolutionary history of P. taurica in Greece is more complex than a single evolutionary invasion. The data analysed, stress the need for a reconsideration of the evolutionary history of Greek Podarcis species and help overcome difficulties that classical taxonomy has encountered at both the species level.     

Variation in the cranium shape of wall lizards (Podarcis spp.): effects of phylogenetic constraints, allometric constraints and ecology

We used geometric morphometrics to explore the influence of phylogenetic and allometric constraints as well as ecology on variation in cranium shape in five species of monophyletic, morphologically similar Podarcis lizards (Podarcis erhardii, Podarcis melisellensis, Podarcis muralis, Podarcis sicula and Podarcis taurica). These species belong to different clades, they differ in their habitat preferences and can be classified into two distinct morphotypes: saxicolous and terrestrial. We found (i) no phylogenetic signal in cranium shape, (ii) diverging allometric slopes among species, and (iii) a significant effect of habitat on cranium shape. The saxicolous species (P. erhardii and P. muralis) had crania with elongated parietals, elongated cranium bases, shortened anterior parts of the dorsal cranium, reduced chambers of the jaw adductor muscles and larger subocular foramina. These cranial features are adaptations that compensate for a flattened cranium, dwelling on vertical surfaces and seeking refuge in crevices. The crania of the terrestrial species (P. melisellensis, P. sicula and P. taurica) tended to be more elongate and robust, with enlarged chambers of the jaw adductor muscle, reduced skull bases and shortened parietals. Terrestrial species exhibited more variation in cranium shape than saxicolous species. Our study suggests that shape variation in Podarcis sp. lizards is largely influenced by ecology, which likely affects species-specific patterns of static allometry.     

Phylogeography of the Italian wall lizard, Podarcis sicula, as revealed by mitochondrial DNA sequences

In a phylogeographical survey of the Italian wall lizard, Podarcis sicula, DNA sequence variation along an 887-bp segment of the cytochrome b gene was examined in 96 specimens from 86 localities covering the distribution range of the species. In addition, parts of the 12S rRNA and 16S rRNA genes from 12 selected specimens as representatives of more divergent cytochrome b haploclades were sequenced (together about 950 bp). Six phylogeographical main groups were found, three representing samples of the nominate subspecies Podarcis sicula sicula and closely related subspecies and the other three comprising Podarcis sicula campestris as well as all subspecies described from northern and eastern Adriatic islands. In southern Italy a population group with morphological characters of P. s. sicula but with the mitochondrial DNA features of P. s. campestris was detected indicating a probably recent hybridization zone. The present distribution patterns were interpreted as the consequence of natural events like retreats to glacial refuges and postglacial area expansions, but also as the results of multiple introductions by man.     

Comparing patterns of nuclear and mitochondrial divergence in a cryptic species complex: the case of Iberian and North African wall lizards (Podarcis, Lacertidae)

Combining different sources of information is essential for a complete understanding of the process of genetic differentiation between species. The Iberian and North African wall lizard ( Podarcis ) species complex has been the object of several studies regarding morphological and mitochondrial DNA variation but, so far, no large-scale survey of nuclear variation within this group has been accomplished. In this study, ten polymorphic allozyme loci were studied in 569 individuals collected across the Iberian Peninsula and North Africa. The obtained data were analysed using both conventional population genetic tools and recent Bayesian model-based clustering methods. Our results show that there are several well-differentiated entities corroborating the major splits observed in mtDNA analyses. These groups correspond not only to the fully recognized species Podarcis bocagei , Podarcis carbonelli , and Podarcis vaucheri but also to multiple forms within the polytypic Podarcis hispanica , all of which have a similar level of differentiation to that observed between the acknowledged species. However, relationships between forms are weakly supported both by population and individual clustering methods, suggesting a scenario of a rapid diversification that contrasts to the clear bifurcating model assumed from previous mtDNA analyses. Individual multilocus analyses report few individuals misassigned or apparently admixed, some of which are most likely explained by the persistence of high levels of ancestral polymorphism. Other admixed individuals, however, are probably the result of limited levels of gene flow between forms.  © 2007 The Linnean Society of London, Biological Journal of the Linnean Society , 2007, 91 , 121–133.     

New satellite DNA in Lacerta s. str. lizards (Sauria: Lacertidae): evolutionary pathways and phylogenetic impact

A new tandemly repeated (satellite) DNA family namely Agi160, from Lacerta agilis and Lacerta strigata (Lacerta sensu stricto (s. str.), Linnaeus 1758) have been cloned and sequenced. Agi160 is found in the above two species, as well as two other representatives of the same genus, L. viridis and L. media. DNA hybridization did not reveal it in Darevskia, Podarcis, Zootoca, Eremias, Ophisops, and Gallotia - the other genera of the family Lacertidae. The results suggest that Agi160 is a Lacerta s. str. specific family of tandem DNA repeats. However, a comparison between sequences of Agi160 and CLsat repeat units revealed 60 bp regions 62-74% identical. The latter is a satellite DNA family typical for Darevskia (syn. "L. saxicola complex") (Grechko et al., Molecular-genetic classification and phylogenetic relatedness of some species of Lacertidae lizards by taxonoprint data. Mol Biol 32:172-183, 1988.). Both Agi160 and CLsat tandem repeats share several common features (e.g., the same AT content and distribution of multiple short A-T runs, internal structure of repeated units, the presence of conservative regions). These data are indicative of their common origin and a possibly strong selective pressure upon conserving both satellites. A comparative analysis of structure, organization, and abundance of these two families of satDNA reveals evolutionary pathways that led to their formation and divergence. The data are consistent with the hypotheses of the concerted evolution of satellite DNA families. The possibility of use of Agi160 as a phylogenetic tool, defining relationships within Lacerta s. str., as well as within the whole family of Lacertidae is discussed.     

Comparative postautotomy tail activity in six Mediterranean lacertid lizard species

Tail autotomy, the self-induced tail separation from the body, is a common and effective antipredator mechanism in lizards. In this study, we examine the muscle energetics of tail shedding in six lacertid lizard species (Podarcis erhardii, Podarcis peloponnesiaca, Podarcis muralis, Podarcis gaigeae, Podarcis milensis, and Lacerta graeca) from the northeast Mediterranean region. Very long periods of postautotomy tail movement were demonstrated for all species (range=6-8 min), and differences among species were not statistically significant. Postautotomy tail movement, powered by anaerobic muscle activity, resulted in a strong increase in lactate concentrations and a concomitant depletion of muscle glycogen of exhausted tails relative to resting tails. No significant differences were found in either lactate or glycogen concentrations among the species examined. Duration of movement was negatively correlated with final lactate concentrations. The lack of differentiation in postautotomy energetic physiology in this group of species that have evolved under very different predation environments indicates that postautotomy muscle metabolism involves an overall conservative suite of characters.     

Contrasting patterns of population subdivision and historical demography in three western Mediterranean lizard species inferred from mitochondrial DNA variation

Pleistocene climatic oscillations were a major force shaping genetic variability in many taxa. We analyse the relative effects of the ice ages across a latitudinal gradient in the Western Mediterranean region, testing two main predictions: (i) species with historical distributions in northern latitudes should have experienced greater loss of suitable habitat, resulting in higher extinction of historical lineages than species distributed in southern latitudes, where the effects of the ice ages were not as drastic. This would be reflected in the observation of lower diversity and number of differentiated lineages in northern areas. (ii) a signature of demographic expansion following the climate amelioration should be obvious in northern species, whereas in the south evidence of long-term effective population size stability should be observed. We used as models three species of wall lizards (Podarcis bocagei, Podarcis carbonelli and Podarcis vaucheri) that replace each other along the study area. We investigated the patterns of mitochondrial DNA diversity and subdivision and obtained demographic parameter estimates for each species. Our results suggest that P. bocagei, the northernmost species, bears low genetic diversity, a shallow coalescent history and marks of a demographic expansion. In contrast, P. vaucheri, the species with a southernmost distribution, shows deeper coalescence events, complex geographical substructure and no evidence for population growth. The species with an intermediate distribution, P. carbonelli, shows average levels of diversity, substructure and population growth. Taken together, these results conform to our main predictions and are explained by a differential influence of the ice ages on distinct latitudes.     

Molecular phylogenetics of Iberian wall lizards (Podarcis): is Podarcis hispanica a species complex?

Phylogenetic relationships between species and morphotypes of Podarcis wall lizards from the Iberian Peninsula and north Africa were estimated using partial 12S rRNA and cytochrome b mitochondrial DNA sequences. All species except Podarcis hispanica form monophyletic units. P. hispanica is paraphyletic, although all identified morphotypes are monophyletic. These morphotypes represent highly divergent lineages showing 10-15% pairwise sequence divergence with the cytochrome b gene. The data suggest that P. hispanica is a species complex. We recommend using P. hispanica* until additional sampling delimits the number and ranges of species currently referred to P. hispanica. P. carbonelli, which has recently been raised to species status, is confirmed as a genetically distinct form. P. atrata is genetically distinct, but much more closely related to some populations of P. hispanica than previously thought.     

Molecular biogeography of the Mediterranean lizards Podarcis Wagler, 1830 and Teira Gray, 1838 (Reptilia, Lacertidae)

Aim   We discuss biogeographical hypotheses for the Mediterranean lizard species Podarcis and Teira within a phylogenetic framework based on partial mitochondrial DNA sequences.
Methods   We derived the most likely phylogenetic hypothesis from our data set (597 aligned positions from the 12S rDNA and phenyl tRNA) under parsimony, distance and maximum likelihood assumptions.
Results   The species usually included in Teira do not form a strongly monophyletic clade. In contrast, the monophyly of the genus Podarcis is rather well supported. Seven lineages are identified in the genus; in order of appearance within the tree, these are: the Balearic pityusensis and lilfordi pair, the sicula complex, a Tyrrhenian tiliguerta and raffonei pair, muralis , the Siculo-Maltese filfolensis and wagleriana pair, the Balkan group ( erhardi , peloponnesiaca , milensis , melisellensis and taurica ), and the Ibero-Maghrebian group ( bocagei , atrata , hispanica and vaucheri ).
Conclusions   The origin of the three European genera of lacertid assayed ( Lacerta , Teira and Podarcis ) is hypothesized to have occurred in the Oligocene. For Podarcis , a possible scenario of a Miocene diversification is derived from the sequence data, and the zoogeography of the lineages are discussed in relation to the palaeogeography of the Mediterranean. It is hypothesized that in the early history of the genus the main lineages separated by rapid, numerous and close events that produced a starting point very similar to a polytomy, hard to resolve by parsimony analysis of the data set.     

Spread of hemoparasites in reptiles and amphibians in Italy

Hemoparasites were harvested from 993 individuals belonging to 15 reptilian and 1 amphibian species, from various Italian localities. Hemogregarins were found in 10 reptilian species while a flagellate and microfilariae were found only in Tarentola mauritanica from Lampedusa. For each host species and place of origin the frequencies of hemogregarins are reported and discussed. Longitudinal studies with periodical thin smears were carried out on 5 Tarentola mauritanica, 4 Lacerta viridis, 26 Podarcis filfolensis, 10 Podarcis muralis, 38 Podarcis sicula, 8 Chalcides ocellatus. This material, whose study has not yet been completed, is made available by the author who strongly encourages further investigations on this subject.     

Anaerobic metabolism and activity of European lacertid lizards

1. Lactate production and physical performance during intense activity were measured in Lacerta vivipara, Podarcis muralis and Podarcis hispanica, temperate lizards from progressively lower latitudes. 2. Rates of lactate production during the initial 30 sec of bursts of activity (means = 1.75 mg g-1) and their Q10S (means = 1.36) are similar to previously obtained values for lizards. 3. Maximal anaerobic scope at activity body temperatures is lowest in Lacerta vivipara, which correlates with a tendency to less strenuous behaviour than that of the Podarcis species and may reflect the need for energy economy. 4. Q10S of maximal performance are lower in Lacerta vivipara and Podarcis muralis than Podarcis hispanica, a difference having no demonstrated metabolic basis but which is thermally adaptive for these cooler climate species.     

Molecular phylogeny and biogeography of the wall-lizard Podarcis erhardii (Squamata: Lacertidae)

Erhard's wall lizard, Podarcis erhardii (Sauria: Lacertidae), is highly diversified in Greece and especially in the southern Aegean region. Out of the 28 recognized subspecies, 27 are found in Greece from the North Sporades island-complex in the North Aegean (grossly south of the 39th parallel) to the island of Crete in the South. The species exhibits great morphological and ecological plasticity and inhabits many different habitats from rocky islets and sandy shores to mountaintops as high as 2000m. By examining intraspecific variability at a segment of the mitochondrial gene cytochrome b we have found that that extant populations of P. erhardii are paraphyletic. Furthermore, we have found that subspecies previously defined on the basis of morphological characteristics do not correspond to different molecular phylogenetic clades, so that their status should be reconsidered. The DNA based biogeographical and phylogenetic history of Podarcis in Southern Greece is congruent with available paleogeographic data of the region, which supports the view that DNA sequences may be a useful tool for the study of palaeogeography.     

Genetic admixture between the Iberian endemic lizards Podarcis bocagei and Podarcis carbonelli: evidence for limited natural hybridization and a bimodal hybrid zone

When recently diverged taxa come into contact, the extent of introgression between them is related to the degree of differentiation that they have achieved. Studying contact zones is therefore essential to understand if differentiated taxa are reproductively isolated and, ultimately, if they are likely to remain distinct. Recent work on Iberian and North African wall lizards ( Podarcis ) has documented the existence of multiple evolutionary units, diagnosable both by genetic markers and morphology, but suggests that gene flow between distinct forms has occurred. Therefore, we were interested in evaluating how species boundaries are maintained in the areas where they meet. In this work, we study the contact zone between Podarcis bocagei and Podarcis carbonelli . We sampled a transect including the only locality where these two lizards are known to occur in syntopy and analysed a battery of 15 unlinked nuclear genetic markers and mitochondrial DNA. We also conducted a preliminary analysis of morphology and fertility. Using model-based clustering approaches, we show that the two species hybridize in the population where they have direct contact, but evidences of introgression are low for nearby populations. Although a significant number of individuals show evidence of admixture, this hybrid zone is clearly bimodal, suggesting strong barriers to gene flow, of which the putative nature are discussed. Interestingly, morphological analyses do not support the existence of intermediate forms among individuals that are admixed genetically. Taken together, these results constitute further evidence validating P. bocagei and P. carbonelli as distinct species.     

Cross-species testing of 27 pre-existing microsatellites in Podarcis gaigeae and Podarcis hispanica (Squamata: Lacertidae)

We tested 27 microsatellite loci for cross-species amplification in the lacertids Podarcis gaigeae and Podarcis hispanica. We detected 11 and 15 polymorphic loci in the former and the latter species, respectively. In a larger sample of individuals from a single population of each species, the number of alleles ranged from five to 23 in 10 of the polymorphic loci in P. gaigeae, and between four and 13 in nine of polymorphic loci in P. hispanica. Two locus deviated from Hardy-Weinberg equilibrium in P. hispanica. Between 11 and 16 of the 27 loci also amplified successfully in three other Podarcis species.     

Comparative study of satellite DNA in ants of the Messor genus

The satellite DNA of ants Messor barbarus and Messor bouvieri is analysed. The results are compared with the satellite DNA data from Messor structor previously reported and with new data obtained from the genome of geographically distinct M. structor population, which have shown that this satellite DNA is highly conserved within the species. The satellite DNA is organized as tandemly repeated 79 bp monomers in all species. The sampled sequences of the three species show a high similarity and all belong to the same family of satellite DNA. Sequence comparisons suggested the occurrence of highly effective homogenization mechanism acting upon the ant genomes. In accordance with this hypothesis, putative gene conversion tracts are identified when the different monomers of the same species are compared. The highest sequence conservation in all species corresponds to a single region with inverted repeats. A CENP-B-like motif was found in this region. The possibility that it may be involved in the homogenization of satellite DNA is discussed.     

Leptin presence in plasma, liver and fat bodies in the lizard Podarcis sicula: fluctuations throughout the reproductive cycle.

Leptin is a 16 kDa peptide produced by adipocytes in response to increasing fat stores and signals to the brain to stop eating and increase energy expenditure. Recent studies point out that the role of leptin is much broader and includes the regulation of reproduction. The lizard, Podarcis sicula, is an oviparous species characterized, at this latitude, by a seasonal reproductive cycle. Since in Podarcis sicula the recovery of the gonadal function coincides with the recovery of metabolic activity, we have hypothesized that leptin might be involved in the regulation of the reproductive function in this species. We have identified an immunoreactive band in the plasma of the female of Podarcis sicula, which comigrated with recombinant mouse leptin and cross reacted with polyclonal antibodies against mouse leptin. A RIA method developed by Linco (St. Louis, MO) was utilized to measure leptin concentration in plasma, liver and fat bodies of Podarcis sicula throughout the reproductive cycle. The antibody used in this kit (Linco's Multi-Species Leptin Radioimmunoassay kit) was produced in the guinea pig against human leptin but displays crossreactivity to leptin molecules of many species. The level of leptin in the plasma of Podarcis sicula was in the same range as that of mammals. Leptin levels in plasma, liver and fat bodies fluctuated during the reproductive cycle, in a way consistent with its possible role in reproduction.     

Higher-Order Organization of Subrepeats and the Evolution of Cervid Satellite I DNA

Based on sequence analyses of 17 complete centromeric DNA monomers from ten different deer species, a model is proposed for the genesis, evolution, and genomic organization of cervid satellite I DNA. All cervid satellite I DNA arose from the initial amplification of a 31-bp DNA sequence. These 31-bp subrepeats were organized in a hierarchical fashion as 0.8-kb monomers in plesiometacarpalia deer and 1-kb monomers in telemetacarpalia deer. The higher-order repeat nature of cervid centromeric satellite DNA monomers accounts for their high intragenomic and intraspecific sequence conservation. Such high intraspecific sequence conservation validates the use of a single cervid satellite I DNA monomer from each deer species for interspecific sequence comparisons to elucidate phylogenetic relationships. Also, a specific 0.18-kb tandem duplication was observed in all 1-kb monomers, implying that 1-kb cervid satellite I DNA monomers arose from an unequal crossover event between two similar 0.8-kb ancestral DNA sequences. Received: 28 May 1996 / Accepted: 24 October 1996     

Cloning, characterization and physical mapping of three cervid satellite DNA families in the genome of the Formosan muntjac (Muntiacus reevesi micrurus)

Due to their high sequence diversity even among closely related species, satellite DNA sequences can be a useful molecular marker for phylogenetic and taxonomic analyses. To characterize the satellite DNA in the genome of a native muntjac species of Taiwan, the Formosan muntjac, satellite DNA clones representing three different cervid satellite DNA families from this species were isolated and analyzed. Genomic organization study of these satellite DNAs was also undertaken. Three Formosan muntjac satellite DNA clones were obtained and designated as FM-satI (1,391 bp), FM-satII (1,143 bp) and FM-satIV (1,103 bp), and found to share approximately 82, 81 and 98% sequence homology with the Chinese muntjac satellite I clone (C5), Indian muntjac satellite II clone (Mmv-0.7) and Chinese muntjac satellite IV clone (MR-1.0), respectively. These three satellite DNA families are organized in a pter<--FM-satII-FM-satIV-FM-satI-->qter orientation in the centromeric region with satII closely associated with the telomeric sequences. Satellite DNA sequence comparison, in combination with chromosome data concludes that the Formosan muntjac is likely a subspecies of M. reevesi, closely related to the Chinese muntjac. With the kinetochore satellite II DNA co-localizing with the telomeric sequences, the Formosan muntjac chromosomes could be truly telocentric.