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.     

Tyrosine phosphorylated proteins in the reproductive tract of the viviparous lizard Eulamprus tympanum and the oviparous lizard Lampropholis guichenoti

Plastic changes occur in the morphology of the uterus at various stages of the reproductive cycle in both oviparous and viviparous lizards and these may be influenced by estrogen. Estrogen driven phosphorylation of effector proteins on tyrosine residues plays a major role in the plastic modulation of uterine anatomy and physiology in vertebrates. We used electrophoresis and Western blotting to characterize the phosphotyrosine protein profiles at various stages of the reproductive pathway in an oviparous lizard Lampropholis guichenoti and a viviparous lizard Eulamprus tympanum. L. guichenoti displayed major bands in the 200-35 kDa range and a triplet of bands of molecular masses 61 kDa, 52 kDa and 48 kDa in 50% of specimens and a 38 kDa band in all specimens. In contrast, E. tympanum samples all displayed a single major band at 40 kDa, which was significantly elevated at the early pregnancy stage. Somewhat paradoxically, the viviparous species, which has the more complex uterine epithelial changes during pregnancy, has the fewest phosphotyrosine bands, so how tyrosine phosphorylation is affected during the evolution of viviparity is not clear.     

HoxA10-like proteins in the reproductive tract of the viviparous lizard Eulamprus tympanum and the oviparous lizard Lampropholis guichenoti

The gene HoxA10 and its protein product are essential for the formation of the extensions of the plasma membrane called uterodomes or pinopods in mammalian uterine epithelia. In mice, the presence of the HoxA10 protein and uterodomes is needed for uterine receptivity to blastocyst implantation. The viviparous lizard Eulamprus tympanum displays uterodomes whereas the oviparous lizard Lampropholis guichenoti does not. To explore the theory that HoxA10 is involved in the formation of uterodomes we investigated whether HoxA10 immunoreactive proteins were present in both species during their reproductive cycles. Oviduct proteins from vitellogenic, gravid or non-reproductive L. guichenoti (n=19) and E. tympanum (n=28) were separated by electrophoresis and analysed by Western blot and specific antibodies to HoxA10. E. tympanum displayed HoxA10 immunoreactive bands at 59 and 63 kDa in 20 out of the 28 samples. All of the L. guichenoti samples displayed HoxA10 immunoreactive bands, 18 had bands at 59 and 64 kDa and 1 animal had a single band at 59 kDa. There were no significant differences in the level of HoxA10 immunoreactivity between the different stages of reproductive cycle in either species. The different molecular mass of the larger band in L. guichenoti (64 kDa) compared to E. tympanum (63 kDa) indicates that the two lizards express different isoforms of the HoxA10-like proteins and it will be interesting in future studies to determine whether there are differences in the biological activity of the proteins that regulate different physiological functions in the uterus of viviparous and oviparous lizards.     

Genetic analyses reveal further cryptic lineages within the Myotis nattereri species complex

In recent years, new cryptic mammalian species have been discovered in Europe, many of which belong to the order Chiroptera. Within this order, some species such as Myotis nattereri contain several cryptic lineages/species, especially in the Mediterranean region. Here we present genetic, phylogenetic and morphological analysis on the Myotis nattereri species complex, focusing on France which is thought to be a contact zone for the different lineages/species. We sequenced the full Cytochrome b gene from individuals from 23 localities and investigated diagnostic morphological characteristics. Our results reveal new phylogenetic relationships within the Myotis nattereri species complex and among closely related species. We discuss morphological characters which are synapomorphic to each lineage and can be used to differentiate them. Our results also demonstrate the presence of a new lineage within the Myotis nattereri species complex. This lineage, endemic to Corsica, possibly represents a new cryptic species for which we present preliminary ecological data. We further identify the presence of three lineages/species in France and detail their distribution with potential contact zones.     

Reproductive biology of viviparous and oviparous species of the leaf beetle genus Oreina

In five species of the genus Oreina Chevrolat (Coleoptera, Chrysomelidae) we compared the size of offspring, the fecundity of the females, the timing of offspring production and female investment over the season. Two of the species, O. elongata and O. luctuosa, laid eggs, while O. cacaliae, O. gloriosa and O. variabilis gave birth to larvae. Offspring size corrected for female size was similar in the two oviparous species and in the viviparous O. cacaliae. In the two other viviparous species the larvae were two to three times bigger in relation to the female. The greater size of the offspring was not traded off for lower fecundity in these latter two species, yet the production of bigger larvae was associated with a longer laying period and thereby a spreading of reproductive investment over the season. The prediction of life history theory that higher investment in individual offspring should be traded off for lower fecundity could not be confirmed. The investigation of egg and larval development showed that in one of the oviparous species, O. luctuosa, the length of the egg stage was more variable. This corroborates the view that in this species the eggs can be retained for varying times before being laid. Greater size at birth does not necessarily lead to shortened developmental times: the larval periods of O. cacaliae, O. elongata, O. gloriosa and O. variabilis were all comparable although the larvae of the first two species were relatively smaller when laid; only the small larvae of O. luctuosa needed significantly longer for their development. For all growth parameters examined the differences between species were larger than the differences between populations. A comparison of larval growth of the oligophagous species O. cacaliae on three plant genera showed that larval growth rate is influenced by the food plant. However, the plant on which the larvae grew worst is apparently not chosen for oviposition in the field. A comparison with a phylogeny of the species based on allozymes suggests that species with similar reproductive parameters are closely related, yet that viviparity evolved independently in O. cacaliae on one hand and O. variabilis and O. gloriosa on the other.     

Development and characterization of 79 nuclear markers amplifying in viviparous and oviparous clades of the European common lizard

The European common lizard (Zootoca vivipara) is a widely distributed species across Europe and Asia exhibiting two reproductive modes (oviparity/viviparity), six major lineages and several sublineages. It has been used to tackle a large variety of research questions, nevertheless, few nuclear DNA sequence markers have been developed for this species. Here we developed 79 new nuclear DNA sequence markers using a clonation protocol. These markers were amplified in several oviparous and viviparous specimens including samples of all extant clades, to test the amplification success and their diversity. 49.4% of the markers were polymorphic and of those, 51.3% amplified in all and 94.9% amplified in 5–7 of the extant Z. vivipara clades. These new markers will be very useful for the study of the population structure, population dynamics, and micro/macro evolution of Z. vivipara. Cross-species amplification in four lizard species (Psammodromus edwardsianus, Podarcis muralis, Lacerta bilineata, and Takydromus sexlineatus) was positive in several of the markers, and six makers amplified in all five species. The large genetic distance between P. edwardsianus and Z. vivipara further suggests that these markers may as well be employed in many other species.     

Interleukin 1 in oviductal tissues of viviparous, oviparous, and ovuliparous species of amphibians

In previous reports, we have shown that interleukin 1 (IL1), a cytokine associated with implantation in mice, is also expressed in reproductive tissues of viviparous squamate reptiles and cartilaginous fishes. In the present study, we investigated the expression of IL1B and its functional membrane receptor type I (IL1R1) in amphibians, a class of vertebrates that is characterized by different reproductive modes, including internal and external fertilization. In particular, we investigated the oviductal tissues of the aplacental viviparous Salamandra lanzai, the oviparous Triturus carnifex, and the ovuliparous Bufo bufo. In immunohistochemistry with anti-human IL1B and IL1R1 polyclonal antibodies we found that in S. lanzai, most cells in the uterine mucosa were immunoreactive for IL1B and IL1R1. In T. carnifex, IL1B and IL1R1 were present in ciliated luminal cells, and there was evidence of IL1B in glandular cells. In B. bufo, the expression of IL1B and IL1R1 was limited to the apical cytoplasm of the ciliated oviductal cells. Western blot analysis showed that a putative mature form of IL1B, similar to that seen in mammals, was present in the oviductal tissues of S. lanzai, whereas different forms, which probably correspond to an inactive pro-IL1B protein, were found in T. carnifex and B. bufo. A band that corresponded to the predicted 80-kDa human IL1R1 was found in S. lanzai and T. carnifex. Although the present study shows that IL1B and IL1R1 expression occurs in all reproductive modes, the differential expression patterns noted between ovuliparity and oviparity and viviparity may reflect the different roles of IL1 in the various reproductive modes.     

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.     

Lipids of the eggs and neonates of oviparous and viviparous lizards

The purpose of this article is to collate the compositional data for the lipids of the eggs and neonates of ten species of lizards displaying a range of parity modes, to highlight emergent trends and to identify some of the physiological changes central to the evolution of viviparity. The eggs of oviparous species and of viviparous species with a simple (type I) placenta are characterised by very high proportions of triacylglycerol which forms over 80% (wt. /wt.) of the total yolk lipid. The eggs of viviparous species with complex (types II and III) placentae contain lower proportions of triacylglycerol (about 70% of total yolk lipid) and commensurately greater proportions of phospholipid, cholesteryl ester and free cholesterol. The fatty acid compositions of the yolk lipids are very similar for all the lizard species, irrespective of parity mode; in particular, the proportions of docosahexaenoic acid are consistently low. For all the species, the proportions of both docosahexaenoic and arachidonic acids are higher in the phospholipid of the neonate compared with the egg. The difference between the lipid contents of the eggs and the neonates indicates that, in species of Pseudemoia which have a complex (type III) placenta, more than 50% of the total lipid supplied to the embryo is derived from placental transport.     

Combining molecular, morphological and ecological data to infer species boundaries in a cryptic tropical pitviper

Few operational methods exist for delimiting species boundaries, and these usually require sampling strategies that are unrealistic for widespread organisms that occur at low densities. Here we apply molecular, morphological and ecological species delimitation criteria to a wide-ranging, fragmented group of Asian green pitvipers, the Popeia popeiorum complex. A mitochondrial DNA phylogeny for the group indicates two well-differentiated clades, corresponding mainly to northern and southern parts of its range. Strong phylogeographical structure within each clade suggests isolation in forest refugia during the Pliocene and a southward colonization of the Sunda islands during the Pleistocene. Multivariate analysis of morphological characters reveals a generally conserved pattern of geographical variation, incongruent with the recovered phylogenetic history. We compare groups delineated by mtDNA variation to morphological and ecological divisions in the complex, and discuss the implications of these for the taxonomy of the group. Discordance between species boundaries inferred from different criteria suggests that combining independent sources of data provides the most reliable estimation of species boundaries in organisms that are difficult to sample in large numbers.  © 2006 The Linnean Society of London, Biological Journal of the Linnean Society , 2006, 87 , 343–364.     

Genetic monandry in 6 viviparous species of true sea snakes

Using a suite of 10 highly variable microsatellite loci, we conducted genetic paternity analyses for 76 embryos in the broods of 12 pregnant females representing 6 viviparous species of true sea snakes (Hydrophis clade) in the family Elapidae. To our surprise, we uncovered no evidence for multiple paternity within any of the clutches despite the fact that the genetic markers showed high intraspecific heterozygosities and as many as 20 conspecific alleles per locus. This outcome stands in sharp contrast to the rather high (but also variable) frequency of multiple paternity previously reported in many other reptilian species. However, because our study appears to be the first assessment of genetic parentage for any sea snake species (and indeed for any member of the elapid clade), we can only speculate as to whether this apparent monandry by females is a broader phylogenetic characteristic of elapid snakes or whether it might relate somehow to the sea snakes' peculiar lifestyle that uniquely combines viviparity with a marine existence.     

Testicular androgens and spermatogenetic cycles in the viviparous lizard

Androstenedione, testosterone and dihydrotestosterone levels were measured in the testis of 36 males of the viviparous lizard throughout a period (from end of May to end of July) characterized by the transition between two spermatogenetic cycles and by very low levels of plasma testosterone. The sudden rise of testicular testosterone and androstenedione in June is concomitant with a degeneration of the seminiferous epithelium. It coincides with a transient appearance of testicular dihydrotestosterone. During the next decline in the levels of testosterone and androstenedione, it occurs a restoration of the seminiferous tubules which resume spermatogenesis (proliferation of spermatogenia and prophase of first meiotic division). The part played by some testicular steroids in the control of spermatogenesis is discussed.     

Genetic analysis of an endemic archipelagic lizard reveals sympatric cryptic lineages and taxonomic discordance

The importance of genetic data in biodiversity conservation is well established, and knowledge of standing genetic variation within and between populations is important for designing conservation strategies. We investigated partitioning of genetic diversity in an endemic lizard (Leiocephalus psammodromus) distributed in the Turks and Caicos archipelago using mtDNA and AFLP data from 259 individuals sampled across 13 islands. Current taxonomy identifies six or more subspecies of L. psammodromus within the archipelago, several of which have undergone recent drastic reductions in range due to extirpation. However, our results indicate the presence of two independent lineages, one on each of the Turks and Caicos banks, and a third sympatric cryptic lineage on both banks. These lineages do not correspond to current taxonomy and alter our understanding of diversity and conservation of this species. Gross morphological data (mass and snout-vent length) indicate some variation in female size among lineages, indicating the possibility of cryptic morphological variation. Instead of initiating separate conservation measures for nominate subspecies, we recommend a more thorough investigation of the morphology and genetics of this group and a more inclusive conservation program. Our surprising results indicate that other endemic squamates in the Bahamas Archipelago might also exhibit sympatric cryptic diversity that does not correspond to current taxonomic understanding and could have significant impacts on our approach to conservation in this region.     

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.     

Geographically widespread mitochondrial lineages of the African saw-wings inconsistent with species boundaries

The African saw-wings (genus Psalidoprocne) are a group of swallows endemic to sub-Saharan Africa. Although currently described as five species with several subspecies, the taxonomy of the saw-wings is unclear. We sequenced two mitochondrial genes (ND2 and cytb; 1 717 sites) from 20 individuals in order to investigate phylogenetic relationships across the range of Psalidoprocne. We obtained similar results from maximum likelihood (RAxML) and Bayesian (BEAST) phylogenetic analyses, which reveal multiple well-supported mitochondrial clades in the genus. The geographic overlap of distinct lineages suggests reproductive isolation, although these clades do not correspond entirely to current species designations.     

Genetic analysis of a contact zone between two lineages of the ocellated lizard (Lacerta lepida Daudin 1802) in south‐eastern Iberia reveal a steep and narrow hybrid zone

Measuring the diffusion of genes between diverging taxa through zones of secondary contact is an essential step to understand the extent and nature of the reproductive isolation that has been achieved. Previous studies have shown that the ocellated lizard (Lacerta lepida Daudin, 1802) has endured repeated range fragmentation associated with the climatic oscillations of the Plio‐Pleistocene that promoted diversification of many different evolutionary units within the species. However, the oldest divergence within the group is estimated to have occurred much earlier, during the Miocene, around 9 Ma and corresponds to the split between the subspecies Lacerta lepida nevadensis Buchholz (1963) and Lacerta lepida lepida Daudin (1802). Although these two evolutionary units have documented genetic and morphological differentiation, most probably accumulated during periods of allopatry, little is known about patterns of gene flow between them. In this study, we performed a population genetic analysis of a putative area of secondary contact between these two taxa, using mtDNA and microsatellite data. We assessed levels of gene flow across the contact zone to clarify to what extent gene flow may be occurring. Hybridization between the subspecies was observed by the presence of genetically introgressed individuals. However, the overall coincidence of mitochondrial and multilocus nuclear clines and generally steep clines support the idea that this contact zone is acting as a barrier to gene flow. Taken together, these results suggest that L. l. lepida and L. l. nevadensis are in independent evolutionary trajectories and should be considered as two different species.     

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.     

Divergent male and female mate preferences do not explain incipient speciation between lizard lineages

Diversification in sexual signals is often taken as evidence for the importance of sexual selection in speciation. However, in order for sexual selection to generate reproductive isolation between populations, both signals and mate preferences must diverge together. Furthermore, assortative mating may result from multiple behavioral mechanisms, including female mate preferences, male mate preferences, and male–male competition; yet their relative contributions are rarely evaluated. Here, we explored the role of mate preferences and male competitive ability as potential barriers to gene flow between 2 divergent lineages of the tawny dragon lizard, Ctenophorus decresii, which differ in male throat coloration. We found stronger behavioral barriers to pairings between southern lineage males and northern lineage females than between northern males and southern females, indicating incomplete and asymmetric behavioral isolating barriers. These results were driven by both male and female mate preferences rather than lineage differences in male competitive ability. Intrasexual selection is therefore unlikely to drive the outcome of secondary contact in C. decresii, despite its widely acknowledged importance in lizards. Our results are consistent with the emerging view that although both male and female mate preferences can diverge alongside sexual signals, speciation is rarely driven by divergent sexual selection alone.