Mitochondrial and nuclear DNA survey of Zootoca vivipara across the eastern Italian Alps: evolutionary relationships, historical demography and conservation implications

Viviparity is associated with larger female size and higher sexual size dimorphism in a reproductively bimodal lizard

Squamate reptiles are central for studying phenotypic correlates of evolutionary transitions from oviparity to viviparity because these transitions are numerous, with many of them being recent. Several models of life-history theory predict that viviparity is associated with increased female size, and thus more female-biased sexual size dimorphism (SSD). Yet, the corresponding empirical evidence is overall weak and inconsistent. The lizard Zootoca vivipara, which occupies a major part of Northern Eurasia and includes four viviparous and two non-sister oviparous lineages, represents an excellent model for testing these predictions. We analysed how sex-specific body size and SSD is associated with parity mode, using body length data for nearly 14,000 adult individuals from 97 geographically distinct populations, which cover almost the entire species' range and represent all six lineages. Our analyses controlled for lineage identity, climatic seasonality (the strongest predictor of geographic body size variation in previous studies of this species) and several aspects of data heterogeneity. Parity mode, lineage and seasonality are significantly associated with female size and SSD; the first two predictors accounted for 14%-26% of the total variation each, while seasonality explained 5%-7%. Viviparous populations exhibited a larger female size than oviparous populations, with no concomitant differences in male size. The variation of male size was overall low and poorly explained by our predictors. Albeit fully expected from theory, the strong female bias of the body size differences between oviparous and viviparous populations found in Z. vivipara is not evident from available data on three other lizard systems of closely related lineages differing in parity mode. We confront this pattern with the data on female reproductive traits in the considered systems and the frequencies of evolutionary changes of parity mode in the corresponding lizard families and speculate why the life-history correlates of live-bearing in Z. vivipara are distinct. Comparing conspecific populations, our study provides the most direct evidence for the predicted effect of parity mode on adult body size but also demonstrates that the revealed pattern may not be general. This might explain why across squamates, viviparity is only weakly associated with larger size.

Genetic basis and expression of ventral colour in polymorphic common lizards

Colour is an important visual cue that can correlate with sex, behaviour, life history or ecological strategies, and has evolved divergently and convergently across animal lineages. Its genetic basis in non-model organisms is rarely known, but such information is vital for determining the drivers and mechanisms of colour evolution. Leveraging genetic admixture in a rare contact zone between oviparous and viviparous common lizards (Zootoca vivipara), we show that females (N = 558) of the two otherwise morphologically indistinguishable reproductive modes differ in their ventral colouration (from pale to vibrant yellow) and intensity of melanic patterning. We find no association between female colouration and reproductive investment, and no evidence for selection on colour. Using a combination of genetic mapping and transcriptomic evidence, we identified two candidate genes associated with ventral colour differentiation, DGAT2 and PMEL. These are genes known to be involved in carotenoid metabolism and melanin synthesis respectively. Ventral melanic spots were associated with two genomic regions, including a SNP close to protein tyrosine phosphatase (PTP) genes. Using genome re-sequencing data, our results show that fixed coding mutations in the candidate genes cannot account for differences in colouration. Taken together, our findings show that the evolution of ventral colouration and its associations across common lizard lineages is variable. A potential genetic mechanism explaining the flexibility of ventral colouration may be that colouration in common lizards, but also across squamates, is predominantly driven by regulatory genetic variation.

Genomic Signatures Associated with Transitions to Viviparity in Cyprinodontiformes

The transition from oviparity to viviparity has occurred independently over 150 times across vertebrates, presenting one of the most compelling cases of phenotypic convergence. However, whether the repeated, independent evolution of viviparity is driven by redeployment of similar genetic mechanisms and whether these leave a common signature in genomic divergence remains largely unknown. Although recent investigations into the evolution of viviparity have demonstrated striking similarity among the genes and molecular pathways involved across disparate vertebrate groups, quantitative tests for genome-wide convergent have provided ambivalent answers. Here, we investigate the potential role of molecular convergence during independent transitions to viviparity across an order of ray-finned freshwater fish (Cyprinodontiformes). We assembled de novo genomes and utilized publicly available genomes of viviparous and oviparous species to test for molecular convergence across both coding and noncoding regions. We found no evidence for an excess of molecular convergence in amino acid substitutions and in rates of sequence divergence, implying independent genetic changes are associated with these transitions. However, both statistical power and biological confounds could constrain our ability to detect significant correlated evolution. We therefore identified candidate genes with potential signatures of molecular convergence in viviparous Cyprinodontiformes lineages. Motif enrichment and gene ontology analyses suggest transcriptional changes associated with early morphogenesis, brain development, and immunity occurred alongside the evolution of viviparity. Overall, however, our findings indicate that independent transitions to viviparity in these fish are not strongly associated with an excess of molecular convergence, but a few genes show convincing evidence of convergent evolution.

Understanding the evolution of viviparity using intraspecific variation in reproductive mode and transitional forms of pregnancy

How innovations such as vision, flight and pregnancy evolve is a central question in evolutionary biology. Examination of transitional (intermediate) forms of these traits can help address this question, but these intermediate phenotypes are very rare in extant species. Here we explore the biology and evolution of transitional forms of pregnancy that are midway between the ancestral state of oviparity (egg-laying) and the derived state, viviparity (live birth). Transitional forms of pregnancy occur in only three vertebrates, all of which are lizard species that also display intraspecific variation in reproductive phenotype. In these lizards (Lerista bougainvillii, Saiphos equalis, and Zootoca vivipara), geographic variation of three reproductive forms occurs within a single species: oviparity, viviparity, and a transitional form of pregnancy. This phenomenon offers the valuable prospect of watching 'evolution in action'. In these species, it is possible to conduct comparative research using different reproductive forms that are not confounded by speciation, and are of relatively recent origin. We identify major proximate and ultimate questions that can be addressed in these species, and the genetic and genomic tools that can help us understand how transitional forms of pregnancy are produced, despite predicted fitness costs. We argue that these taxa represent an excellent prospect for understanding the major evolutionary shift between egg-laying and live birth, which is a fundamental innovation in the history of animals.

The functional genetic architecture of egg-laying and live-bearing reproduction in common lizards

All amniotes reproduce either by egg-laying (oviparity), which is ancestral to vertebrates or by live-bearing (viviparity), which has evolved many times independently. However, the genetic basis of these parity modes has never been resolved and, consequently, its convergence across evolutionary scales is currently unknown. Here, we leveraged natural hybridizations between oviparous and viviparous common lizards (Zootoca vivipara) to describe the functional genes and genetic architecture of parity mode and its key traits, eggshell and gestation length, and compared our findings across vertebrates. In these lizards, parity trait genes were associated with progesterone-binding functions and enriched for tissue remodelling and immune system pathways. Viviparity involved more genes and complex gene networks than did oviparity. Angiogenesis, vascular endothelial growth and adrenoreceptor pathways were enriched in the viviparous female reproductive tissue, while pathways for transforming growth factor were enriched in the oviparous. Natural selection on these parity mode genes was evident genome-wide. Our comparison to seven independent origins of viviparity in mammals, squamates and fish showed that genes active in pregnancy were related to immunity, tissue remodelling and blood vessel generation. Therefore, our results suggest that pre-established regulatory networks are repeatedly recruited for viviparity and that these are shared at deep evolutionary scales.

Genetic diversity of major histocompatibility complex class I genes in Zootoca vivipara

The Major Histocompatibility Complex (MHC), as a family of highly polymorphic genes associated with immunity in the genome of the vertebrate, has become an important indicator for assessing the evolutionary potential of wildlife. In order to better protect Zootoca vivipara in the Greater Khingan Range and Lesser Khingan Range, to understand the genetic structure of Z. vivipara, and to explore the mechanism and phylogenetic relationship of the gene polymorphisms, the MHC molecular marker method was used to analyze Z. vivipara population. Forty-seven alleles were obtained from four populations. The four populations were highly polymorphic, rich in genetic information, and had significant genetic diversity. There were certain inbreeding phenomena. There was a high degree of genetic differentiation among populations, which was caused by genetic drift and natural selection. The sequence undergoes genetic duplication and recombination. The existence of trans-species polymorphism was found in the constructed phylogenetic tree. The present study provides a theoretical basis for species protection of Z. vivipara.

Chromosome-Level Assembly of the Common Lizard (Zootoca vivipara) Genome

Squamate reptiles exhibit high variation in their phenotypic traits and geographical distributions and are therefore fascinating taxa for evolutionary and ecological research. However, genomic resources are very limited for this group of species, consequently inhibiting research efforts. To address this gap, we assembled a high-quality genome of the common lizard, Zootoca vivipara (Lacertidae), using a combination of high coverage Illumina (shotgun and mate-pair) and PacBio sequencing data, coupled with RNAseq data and genetic linkage map generation. The 1.46-Gb genome assembly has a scaffold N50 of 11.52 Mb with N50 contig size of 220.4 kb and only 2.96% gaps. A BUSCO analysis indicates that 97.7% of the single-copy Tetrapoda orthologs were recovered in the assembly. In total, 19,829 gene models were annotated to the genome using a combination of ab initio and homology-based methods. To improve the chromosome-level assembly, we generated a high-density linkage map from wild-caught families and developed a novel analytical pipeline to accommodate multiple paternity and unknown father genotypes. We successfully anchored and oriented almost 90% of the genome on 19 linkage groups. This annotated and oriented chromosome-level reference genome represents a valuable resource to facilitate evolutionary studies in squamate reptiles.

Reversals in complex traits uncovered as reticulation events: Lessons from the evolution of parity-mode, chromosome morphology, and maternal resource transfer

Complex traits include, among many others, the evolution of eyes, wings, body forms, reproductive modes, human intelligence, social behavior, diseases, and chromosome morphology. Dollo's law states that the evolution of complex traits is irreversible. However, potential exceptions have been proposed. Here, we investigated whether reticulation, a simple and elegant means by which complex characters may be reacquired, could account for suggested reversals in the evolution of complex characters using two datasets with sufficient genetic coverage and a total of five potential reversals. Our analyses uncovered a potential reversal in the evolution of parity mode and a potential reversal in the evolution of placentotrophy of fish (Cyprinodontiformes) as reticulation events. Moreover, in a reptile that exhibits a potential reversal from viviparity to oviparity (Zootoca vivipara), reticulation provided the most parsimonious explanation for sex chromosome evolution. Therefore, three of the five studied potential reversals were unraveled as reticulation events. This constitutes the first evidence that accounting for reticulation can fundamentally influence the interpretation of the evolution of complex traits, that testing for reticulation is crucial for obtaining robust phylogenies, and that complex ancestral characters may be reacquired through hybridization with a lineage that still exhibits the trait. Hybridization, rather than reappearance of ancestral traits by means of small evolutionary steps, may thus account for suggested exceptions to Dollo's law. Consequently, ruling out reticulation is required to claim the evolutionary reversal of complex characters and potential exceptions to Dollo's rule.

Museum specimens indicate genetic erosion in an endangered lizard

Genetic variability, one of the main factors that guarantees species persistence, and species’ conservation status are generally evaluated with indices calculated at the present time. Natural history collections might help compare historical and current genetic diversity so to identify major trends. Here we analysed museum specimens of the lizard Zootoca vivipara carniolica, with a specific and stringent protocol for degraded DNA, in order to contrast its past and current genetic variability, using fragments of one mitochondrial DNA gene. Part of the distributional range of Z. v. carniolica (Po Plain, Italy), heavily impacted by human activities, was investigated. We found two previously unknown haplotypes in populations that are extinct today, suggesting the loss of these haplotypes and thus an overall shrinking of genetic variability. We argue that these results, together with the increasing threats posed by climate and land use changes, suggest that specific conservation measures for the persistence of Z. v. carniolica in Northern Italian lowlands have to be considered.

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.

Very similar, but not identical: morphological taxonomic identification to improve the resolution of fine-scale distribution of Zootoca (vivipara) carniolica

Improving the knowledge about the distribution of threatened species is a fundamental step forward for implementing effective measures for species conservation. Distribution maps of herpetological species are often inferred from morphological identification but, in some cases, morphology alone is not reliable for assessing the distribution of look-alike species. The common lizard Zootoca vivipara includes several lineages which, to date, cannot be distinguished by any conventional morphological traits. Therefore, a reliable taxonomic identification for the main Zootoca lineages is only feasible with genetic data. In this study, we identify a morphological trait that can be used to distinguish individuals of the parapatric Z. v. vivipara and Z. carniolica, without the support of molecular analyses. This method may provide a cost-effective tool to better depict the distribution of the threatened Z. carniolica and therefore prioritize the conservation of the locations where the oviparous Alpine lineage is found.

Genetic and ecological data reveal species boundaries between viviparous and oviparous lizard lineages

Identification of cryptic species is an essential aim for conservation biologists to avoid premature extinctions of 'unrecognized' species. Integrating different types of data can undoubtedly aid in resolving the issue of species delimitation. We studied here two lineages of the common lizard Zootoca vivipara that display different reproductive mode (the viviparous Z. v. vivipara and the oviparous Z. v. carniolica) and that overlap their distributional ranges in the European Alps. With the purpose of delimiting species' boundaries, we analyzed their ecological, genetic and natural history features. More than 300 samples were collected and analyzed at cytochrome b and 11 microsatellites loci for investigating genetic variation, population structure, individual relatedness and evolutionary histories of the two lineages. Additionally, we compared their ecological niches using eight ecological variables. Genetic data showed contrasting patterns of genetic structure between the two lineages, different demographic dynamics and no hybridization events. Also strong ecological differences (such as temperature) emerged between the two lineages, and niche overlap was limited. Taken together, these results indicate that Z. v. vivipara and Z. v. carniolica should be recognized as two separate species, and particular conservation consideration should be given to the oviparous lineage that tends to live in areas threatened by increasing impact of human activities. However, recent and rapid climate warming might determine an increasing risk for the persistence of the viviparous lineage, being adapted to cold environments.