Refuge areas and suture zones in the Pyrenean and Cantabrian regions: geographic variation of the female MPI sex-linked alleles among oviparous populations of the lizard Lacerta (Zootoca) vivipara

The viviparous lizard Lacerta (Zootoca) vivipara exhibits several alleles of the mannose-6-phosphate isomerase (MPI) enzyme that are carried exclusively on the female W sex-chromosome. Previous studies showed that both the oviparous and viviparous forms of L. (Zootoca) vivipara have these female sex-linked alleles. We document the existence of geographic variation of these alleles among the oviparous populations of southwestern France and northwestern Spain. Two oviparous subgroups were identified: all females from the eastern and central Pyrenees and most females from Aquitaine and from the northern slope of the western Pyrenees exhibited the fast migrating alleles MPI¹¹⁰ or MPI¹²⁰, whereas all females from the Cantabric mountains, Spanish Basque country, and from the south slope of the western Pyrenees exhibited the slow migrating allele MPI⁹⁰. Populations with both fast and slow migrating alleles occurred at some stations in the upper Ossau valley (western Pyrenees) and also at a lowland station of south Aquitaine. The hypothesis that several oviparous forms could have retreated to different places of the Pyreneo-Iberian refugia during the Quaternary glaciations could explain the conservation or the evolution of the polymorphism of the MPI alleles, and that is consistent with the phylogeographic scenario previously proposed to account for the reproductive and cytogenetical variation observed in this species.     

Further evidence of the existence of oviparous populations of Lacerta (Zootoca) vivipara in the NW of the Balkan Peninsula

The lizard Lacerta (Zootoca) vivipara, which is viviparous in the greatest part of its distribution range, has however some oviparous populations on the southern margin of its range. The present study aimed at determining the reproductive mode and the ATA (aspartate transaminase) enzyme characteristics of four populations in Slovenia and one population in Croatia. The Slovenian females studied here presented an oviparous reproductive mode which strongly resembled those observed in the oviparous populations of south-western France and north-western Spain. Our electrophoresis analyses revealed the existence of two distinct alleles, ATA-150 and ATA-200, in the oviparous populations of Slovenia. These alleles were identical to those observed in the French and Spanish oviparous group and were distinct from the allele ATA-100 characterizing the viviparous populations that we had previously studied. Although we did not study the reproductive mode of Croatian females, the allele ATA-200 observed in one population of Croatia strongly suggested that this population might also be oviparous.     

Intraspecific phylogeography of Lacerta vivipara and the evolution of viviparity

The lacertid lizard Lacerta vivipara is one of the few squamate species with two reproductive modes. We present the intraspecific phylogeny obtained from neighbor-joining and maximum-parsimony analyses of the mtDNA cytochrome b sequences for 15 individuals from Slovenian oviparous populations, 34 individuals from western oviparous populations of southern France and northern Spain, 92 specimens from European and Russian viviparous populations, and 3 specimens of the viviparous subspecies L. v. pannonica. The phylogeny indicates that the evolutionary transition from oviparity to viviparity probably occurred once in L. vivipara. The western oviparous group from Spain and southern France is phylogenetically most closely related to the viviparous clade. However, the biarmed W chromosome characterizing the western viviparous populations is an apomorphic character, whereas the uniarmed W chromosome, existing both in the western oviparous populations and in the geographically distant eastern viviparous populations, is a plesiomorphic character. This suggests an eastern origin of viviparity. Various estimates suggest that the oviparous and viviparous clades of L. vivipara split during the Pleistocene. Our results are discussed in the framework of general evolutionary models: the concept of an oviparity-viviparity continuum in squamates, the cold climate model of selection for viviparity in squamates, and the contraction-expansion of ranges in the Pleistocene resulting in allopatric differentiation.     

Evolutionary and biogeographical implications of the karyological variations in the oviparous and viviparous forms of the lizard Lacerta (Zootoca) vivipara

The lizard Lacerta ( Zootoca ) vivipara has two modes of reproduction and is variable karyologically. We describe its karyological variation from literature data and from new data on two viviparous populations from France, on two oviparous populations from the Pyrenees in south-western France and on three oviparous populations recently discovered in Slovenia. Males have 36 chromosomes, whereas females have only 35 chromosomes in all viviparous populations and in the Pyrenean oviparous populations. This karyotype has been interpreted to result from a fusion of an ancestral sexual W chromosome with an autosome from the Zl or from the Z2 pair. The karyotype formula is 32 autosomes + ZIZ2W for the female and 32 autosomes + Z1Z1Z2Z2 for the male. The karyotype of the Slovenian oviparous populations, 34 autosomes + ZW in the male and 34 autosomes +ZW in the female, represents an evolutionary stage that preceded the chromosomal fusion. There is minor karyological variation, mainly concerning the W and Z2 chromosomes, within the Pyrenean oviparous populations. This parallels the geographic variation of the W-linked alleles of the MPI enzyme and suggests that allopatric differentiation of these oviparous populations might have occurred in the vicinity of the Pyrenees during the Pleistocene.
The viviparous populations from western Europe carry a metacentric W chromosome, whereas oviparous populations from south-western Europe and eastern viviparous populations both show an acrocentrie, or a subtelocentrie. W chromosome. This suggests that the acrocentric-subtelocentric W is a primitive character and that viviparity probably arose in an eastern lineage of the species.     

Comparisons of mitochondrial DNA (mtDNA) sequences (16S rRNA gene) between oviparous and viviparous strains of Lacerta vivipara: a preliminary study

The lizard Lacerta vivipara has allopatric oviparous and viviparous populations. The mitochondrial DNA (mtDNA) gene coding for the 16S rRNA was sequenced for several viviparous lizard populations from France, Switzerland, Bulgaria, Czech Republic, The Netherlands, Sweden, and for oviparous lizard populations from the Pyrenean and Cantabric Mountains. Seven distinct groups (three oviparous and four viviparous) were identified. The net nucleotide divergence between oviparous and viviparous haplotypes was 1.3% +/- 0.5 (mean +/- standard deviation). These results on mtDNA, together with other data obtained previously, led us to formulate a biogeographical scenario that could be tested by further research.     

Intraspecific Phylogeography of Lacerta vivipara and the Evolution of Viviparity

The lacertid lizard Lacerta vivipara is one of the few squamate species with two reproductive modes. We present the intraspecific phylogeny obtained from neighbor-joining and maximum-parsimony analyses of the mtDNA cytochrome b sequences for 15 individuals from Slovenian oviparous populations, 34 individuals from western oviparous populations of southern France and northern Spain, 92 specimens from European and Russian viviparous populations, and 3 specimens of the viviparous subspecies L. v. pannonica. The phylogeny indicates that the evolutionary transition from oviparity to viviparity probably occurred once in L. vivipara. The western oviparous group from Spain and southern France is phylogenetically most closely related to the viviparous clade. However, the biarmed W chromosome characterizing the western viviparous populations is an apomorphic character, whereas the uniarmed W chromosome, existing both in the western oviparous populations and in the geographically distant eastern viviparous populations, is a plesiomorphic character. This suggests an eastern origin of viviparity. Various estimates suggest that the oviparous and viviparous clades of L. vivipara split during the Pleistocene. Our results are discussed in the framework of general evolutionary models: the concept of an oviparity–viviparity continuum in squamates, the cold climate model of selection for viviparity in squamates, and the contraction–expansion of ranges in the Pleistocene resulting in allopatric differentiation.     

Comparison of the cold hardiness capacities of the oviparous and viviparous forms of Lacerta vivipara

The lizard Lacerta vivipara has allopatric oviparous and viviparous populations. The cold hardiness strategy of L. vivipara has previously been studied in viviparous populations, but never in oviparous ones. The present study reveals that both the oviparous and viviparous individuals of this species are able to survive in a supercooled state at -3 degrees C for at least one week when kept on dry substrates. The mean crystallisation temperatures of the body, around -4 degrees C on dry substrata and -2 degrees C on wet substrata, do not differ between oviparous and viviparous individuals. All the individuals are able to tolerate up to 48-50% of their body fluid converted into ice, but only viviparous individuals were able to stabilize their body ice content at 48%, and hence were able to survive even when frozen at -3 degrees C for times of up 24 hours. Ice contents higher than 51% have been constantly found lethal for oviparous individuals. This suggests that, in L. vivipara, the evolution towards a higher degree of freezing tolerance could parallel the evolution of the viviparous reproductive mode, a feature believed to be strongly selected under cold climatic conditions. This is the first report, among reptiles, of an intraspecific variation regarding the freeze tolerance capacities.     

Mating does not influence reproductive investment, in a viviparous lizard

Mating is crucial for females that reproduce exclusively sexually and should influence their investment into reproduction. Although reproductive adjustments in response to mate quality have been tested in a wide range of species, the effect of exposure to males and mating per se has seldom been studied. Compensatory mechanisms against the absence of mating may evolve more frequently in viviparous females, which pay higher direct costs of reproduction, due to gestation, than oviparous females. To test the existence of such mechanisms in a viviparous species, we experimentally manipulated the mating opportunity of viviparous female lizard, Lacerta (Zootoca) vivipara. We assessed the effect of mating on ovulation, postpartum body condition and parturition date, as well as on changes in locomotor performances and body temperatures during the breeding cycle. Female lizards ovulated spontaneously and mating had no influence on litter size, locomotor impairment or on selected body temperature. However, offspring production induced a more pronounced locomotor impairment and physical burden than the production of undeveloped eggs. Postpartum body condition and parturition dates were not different among females. This result suggests that gestation length is not determined by an embryonic signal. In the common lizard, viviparity is not associated with facultative ovulation and a control of litter size after ovulation, in response to the absence of mating.     

Cytokines in the oviparity/viviparity transition: evidence of the interleukin-1 system in a species with reproductive bimodality, the lizard Lacerta vivipara

Placental viviparity is a reproductive strategy usually attributed to mammals. However, it is also present in other vertebrate species, e.g. in Squamate reptiles. Although the immunological mechanisms that allow the survival of the semi-allogenic embryo in maternal tissues are still largely unknown, cytokines seem to play an important role in mammalian reproduction. Previous studies in our laboratory showed that interleukin-1 (IL-1), a cytokine associated with implantation in mice, is also expressed at the materno-fetal interface of placental viviparous Squamates. In this study, we used the model of Lacerta vivipara, which exhibits reproductive bimodality, that is, the coexistence of oviparous and viviparous populations. By means of immunohistochemistry and anti-human antibodies, we showed that uterine tissues of L. vivipara (seven oviparous and six viviparous animals) expressed the two IL-1 isoforms, IL-1alpha and IL-1beta, and the type I IL-1 receptor (IL-1R tI) both at the pre-ovulatory stage and during gestation, with no significant difference between oviparous and viviparous females. In L. vivipara, as in most oviparous Squamates, an important phase of embryonic development takes place in the mother's oviduct, before egg-laying. Moreover, although thinner than in oviparous females, an eggshell membrane persists throughout gestation in viviparous females also, which develop a very simple type of placenta. The data suggest that immunological mechanisms that allow the survival of the semi-allogenic embryo in maternal tissues are independent of the timing or intimacy of contact between maternal and fetal tissues.     

Karyological investigations of populations of the lizard Zootoca vivipara (Juaquin, 1787) from the Baltic Sea basin (western region of Russia)

Populations of the lizard Zootoca vivipara (Juaquin, 1787) (Lacertidae, Sauria) from Western and Central Europe are poorly differentiated in morphology. However, they differ from each other by several karyotype markers and their reproductive mode (oviparous and viviparous types). Here we report for the first time the results of combined investigations of Z. vivipara from Eastern Europe (the Baltic Sea basin of Western Russia). The karyotype and reproductive mode of the specimens from four previously not examined populations in Kaliningrad region have been studied. It has been shown that these lizards have viviparous mode of reproduction and their karyotypes vary on morphology of W-sex chromosomes. Females possess 2n = 35 : 32 A + Z1Z2W, with W--SV and 2n = 35 : 32 A + Z1Z2W, with W--A/ST. We have identified these populations karyologically and have first revealed that they belong to two different viviparous forms of nominative subspecies Z. v. vivipara, namely to its "western" and "Russian (eastern)" forms. These data have confirmed that the Baltic Sea basin in Eastern Europe is a zone of secondary contact of these chromosomal forms of Z. v. vivipara. The results obtained do not now allow us to speak about their sympatry, parapatry and (or) hybridization events. The results suggest that during postglasial time "western" form might colonize this region from western Europe whereas "Russian (eastern)" form might do it from Southern-Eastern Europe. It has been recommended to save the biodiversity of Z. vivipara in Kaliningrad region, namely the peripheral and very rare on vast territory of Russia "western" form and unique for western region of Russia "Russian (eastern)" form.     

Experimental evidence of early costs of reproduction in conspecific viviparous and oviparous lizards

Reproduction entails costs, and disentangling the relative importance of each stage of the reproductive cycle may be important to assess the costs and benefits of different reproductive strategies. We studied the early costs of reproduction in oviparous and viviparous lizard females of the bimodal reproductive species Zootoca vivipara. Egg retention time in oviparous females is approximately one-third of the time in viviparous females. We compared the vitellogenesis and egg retention stages that are common to both reproductive modes. Precisely, we monitored the thermoregulatory behaviour, the weight gain and the immunocompetence of the females. Moreover, we injected an antigen in half of the females (immune challenge) to study the trade-offs between reproductive mode and immune performance and between different components of the immune system. Finally, we experimentally induced parturition in viviparous females at the time of egg laying in oviparous females. Oviparous and viviparous females did not show strong differences in response to the immune challenge. However, viviparous females spent more time thermoregulating while partially hidden and gained more weight than oviparous females. The greater weight gain indicates that the initial period of egg retention is less costly for viviparous than for oviparous females or that viviparous females are able to save and accumulate energy at this period. This energy may be used by viviparous females to cope with the subsequent costs of the last two-third of the gestation. Such an ability to compensate the higher costs of a longer egg retention period may account for the frequent evolution of viviparity in squamate reptiles.     

Reproductive isolation between oviparous and viviparous lineages of the Eurasian common lizard Zootoca vivipara in a contact zone

Contact zones between two evolutionary lineages are often useful for understanding the process of speciation because the observed genetic pattern reflects the history of differentiation. The Eurasian lacertid lizard Zootoca vivipara is a potentially interesting model for studying the role of reproductive mode in the speciation of squamate reptiles because it has both oviparous (Zootoca vivipara carniolica) and viviparous (Zootoca vivipara vivipara) populations that have recently been shown to be genetically distinct. We studied a newly‐discovered syntopic area of these two Zootoca subspecies in the central Italian Alps using genetic markers to investigate the level of introgression between them. Patterns of genetic differentiation in a fragment of the mitochondrial DNA cytb gene and a set of nuclear microsatellites show that the speciation process is complete in this area, with no evidence of recent introgression. Phylogenetic and genotypic divergence suggests that the two subspecies have experienced long independent evolutionary histories, during which genetic and phenotypic differences evolved. The possible roles of biogeography, reproductive mode, and cytogenetic differentiation in this speciation process are discussed. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 114 , 566–573.     

Multiple origins of viviparity,or reversal from viviparity to oviparity? The European common lizard (Zootoca vivipara,Lacertidae) and the evolution of parity

The evolution of viviparity in squamates has been the focus of much scientific attention in previous years. In particular, the possibility of the transition from viviparity back to oviparity has been the subject of a vigorous debate. Some studies have suggested this reversal is more frequent than previously thought. However, none of them provide conclusive evidence. We investigated this problem by studying the phylogenetic relationships between oviparous and viviparous lineages of the reproductively bimodal lizard species Zootoca vivipara . Our results show that viviparous populations are not monophyletic, and that several evolutionary transitions in parity mode have occurred. The most parsimonious scenario involves a single origin of viviparity followed by a reversal back to oviparity. This is the first study with a strongly supported phylogenetic framework supporting a transition from viviparity to oviparity.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 87 , 1–11.     

Length of activity season drives geographic variation in body size of a widely distributed lizard

Understanding the factors that drive geographic variation in life history is an important challenge in evolutionary ecology. Here, we analyze what predicts geographic variation in life‐history traits of the common lizard, Zootoca vivipara, which has the globally largest distribution range of all terrestrial reptile species. Variation in body size was predicted by differences in the length of activity season, while we found no effects of environmental temperature per se. Females experiencing relatively short activity season mature at a larger size and remain larger on average than females in populations with relatively long activity seasons. Interpopulation variation in fecundity was largely explained by mean body size of females and reproductive mode, with viviparous populations having larger clutch size than oviparous populations. Finally, body size‐fecundity relationship differs between viviparous and oviparous populations, with relatively lower reproductive investment for a given body size in oviparous populations. While the phylogenetic signal was weak overall, the patterns of variation showed spatial effects, perhaps reflecting genetic divergence or geographic variation in additional biotic and abiotic factors. Our findings emphasize that time constraints imposed by the environment rather than ambient temperature play a major role in shaping life histories in the common lizard. This might be attributed to the fact that lizards can attain their preferred body temperature via behavioral thermoregulation across different thermal environments. Length of activity season, defining the maximum time available for lizards to maintain optimal performance, is thus the main environmental factor constraining growth rate and annual rates of mortality. Our results suggest that this factor may partly explain variation in the extent to which different taxa follow ecogeographic rules.     

Variation in body size and sexual size dimorphism in the most widely ranging lizard: testing the effects of reproductive mode and climate

Reproductive mode, ancestry, and climate are hypothesized to determine body size variation in reptiles but their effects have rarely been estimated simultaneously, especially at the intraspecific level. The common lizard (Zootoca vivipara) occupies almost the entire Northern Eurasia and includes viviparous and oviparous lineages, thus representing an excellent model for such studies. Using body length data for >10,000 individuals from 72 geographically distinct populations over the species' range, we analyzed how sex‐specific adult body size and sexual size dimorphism (SSD) is associated with reproductive mode, lineage identity, and several climatic variables. Variation in male size was low and poorly explained by our predictors. In contrast, female size and SSD varied considerably, demonstrating significant effects of reproductive mode and particularly seasonality. Populations of the western oviparous lineage (northern Spain, south‐western France) exhibited a smaller female size and less female‐biased SSD than those of the western viviparous (France to Eastern Europe) and the eastern viviparous (Eastern Europe to Far East) lineages; this pattern persisted even after controlling for climatic effects. The phenotypic response to seasonality was complex: across the lineages, as well as within the eastern viviparous lineage, female size and SSD increase with increasing seasonality, whereas the western viviparous lineage followed the opposing trends. Altogether, viviparous populations seem to follow a saw‐tooth geographic cline, which might reflect the nonmonotonic relationship of body size at maturity in females with the length of activity season. This relationship is predicted to arise in perennial ectotherms as a response to environmental constraints caused by seasonality of growth and reproduction. The SSD allometry followed the converse of Rensch's rule, a rare pattern for amniotes. Our results provide the first evidence of opposing body size—climate relationships in intraspecific units.     

Thermal adaptive significance of ADH and EST-6 allozymes in Indian geographical populations of Drosophila melanogaster

Indian geographical populations of Drosophila melanogaster exhibit significant correlation (r 0.95) of allelic frequencies at Est -6 and Adh loci with latitude as well as altitude. Est -6^S and AdhF allozymes are well adapted to colder environments while Est -6^F and AdhS are warm adapted. The data on allozymic clines match climatic conditions on the Indian subcontinent. On the basis of multiple regression analysis, one major conclusion is that coefficient of variation of temperature ( T _CV) along latitude/altitude accounts for alterations in allelic frequency at the Adh locus while T _max and T _max explain changes at the Est -6 locus. Thus, climatic conditions lead to thermal selection of allozymes in Indian populations of D. melanogaster .     

Development of the uterine shell glands during the preovulatory and early gestation periods in oviparous and viviparous Lacerta vivipara

The evolutionary process leading to the emergence of viviparity in Squamata consists of lengthening the period of egg retention in utero coupled with marked reduction in the thickness of the eggshell. We used light microscopy and scanning electron microscopy to study uterine structure during the reproductive cycle of oviparous and viviparous females of the reproductively bimodal Lacerta vivipara. We compared the structure of the uterine shell glands, which secrete components of the eggshell, during preovulatory and early gestation phases of the reproductive cycle and also compared histochemistry of the eggshells. The uterine glands of both reproductive forms undergo considerable growth within a period of a few weeks during folliculogenesis and vitellogenesis preceding ovulation. The majority of the proteinaceous fibers of the shell membrane are secreted early in embryonic development and the uterine glands regress shortly thereafter. This supports previous observations indicating that, in Squamata, secretion of the shell membrane occurs very rapidly after ovulation. The most striking differences between reproductive modes were larger uterine glands at late vitellogenesis in oviparous females, 101 microm compared to 60 microm in viviparous females, and greater thickness of the shell membrane during early gestation in oviparous females (52-73 microm) compared to viviparous females (4-8 microm). Our intraspecific comparison supports the conclusions of previous studies that, prior to ovulation, the uterine glandular layer is less developed in viviparous than in oviparous species, and that this is the main factor accounting for differences in the thickness of the shell membrane of the two reproductive forms of squamates.     

Chemical composition of femoral secretions of oviparous and viviparous types of male common lizards Lacerta vivipara

In spite of the importance of chemoreception in intraspecific communication of lizards, only a few studies have examined chemical composition of secretions of lizards. The secretion of the femoral glands of adult male lizards Lacerta vivipara contains a relatively low number (18) of lipophilic compounds in comparison with other related lacertid lizards. These compounds were identified on the basis of mass spectra, obtained by GC-MS. Chemicals included ten steroids (mainly cholesterol) and four carboxylic acids between n-C₁₂ and n-C₁₈, and minor components such as squalene, α-tocopherol, and two waxy esters, which may contribute to avoid oxidation of other lipophilic components in the fairly humid environments occupied by this lizard. Secretions of adult males from oviparous and viviparous populations did not differ in the numbers and quality of chemical compounds, but there were some differences in the relative proportion of some compounds. Males from oviparous populations had lower proportions of hexadecanoic acid and cholestan-3-one, and higher proportions of squalene than viviparous males. These differences might be explained by either small genetic differences between types or due to different microclimatic conditions in the original populations.     

From oviparity to viviparity: a preliminary note on the morphometric differentiation between oviparous and viviparous species assigned to the genus Liolaemus (Reptilia, Squamata, Liolaemidae)

Embryonic growth requires a considerable internal space in viviparous female lizards and this need for space should be reflected in their external morphometry. External morphological differences associated with the reproductive mode in 12 viviparous and 18 oviparous species of Liolaemus lizards were identified. Size differences between viviparous and oviparous species were elucidated by axilla-groin/snout-vent relationship. Axilla-groin distance, considered a size estimator of visceral cavity, surpassed 50% of snout-vent length in viviparous females, while it is always less than 50% in oviparous females. This difference between the two reproductive modes is statistically significant.     

Energy consumption by embryos of a viviparous lizard, Eulamprus tympanum, during development

Energy consumption during development has been measured in many oviparous lizards, but not in viviparous lizards in utero. It has always been assumed that energy consumption by embryos of viviparous lizards during development is similar to that of oviparous species. Estimation of energy consumption of viviparous lizards in vivo are confounded by the possible influence of pregnancy on maternal metabolism. Here we separated maternal and embryonic metabolism in measurements of pregnant Eulamprus tympanum throughout pregnancy. Our data support the hypothesis that the energetic cost of development in viviparous lizards (19.8 kJ g(-1)) is similar to that in oviparous lizards (mean 16.2 kJ g(-1)), at least for a species with a simple placenta. An increase in maternal metabolism of 29% above that for non-pregnant E. tympanum goes to maintain pregnancy, and represents an important component of the reproductive effort in E. tympanum.