High-Quality, Chromosome-Level Reference Genomes of the Viviparous Caribbean Skinks Spondylurus nitidus and S. culebrae.

New World mabuyine skinks are a diverse radiation of morphologically cryptic lizards with unique reproductive biologies. Recent studies examining population-level data (morphological, ecological, and genomic) have uncovered novel biodiversity and phenotypes, including the description of dozens of new species and insights into the evolution of their highly complex placental structures. Beyond the potential for this diverse group to serve as a model for the evolution of viviparity in lizards, much of the taxonomic diversity is concentrated in regions experiencing increasing environmental instability from climate and anthropogenic change. Consequently, a better understanding of genome structure and diversity will be an important tool in the adaptive management and conservation of this group. Skinks endemic to Caribbean islands are particularly vulnerable to global change with several species already considered likely extinct and several remaining species either endangered or threatened. Combining PacBio long-read sequencing, Hi-C, and RNAseq data, here we present the first genomic resources for this group by describing new chromosome-level reference genomes for the Puerto Rican Skink Spondylurus nitidus and the Culebra Skink S. culebrae. Results indicate two high quality genomes, both ∼1.4 Gb, assembled nearly telomere to telomere with complete mitochondrion assembly and annotation.

Insight into the function of voltage-sensing phosphatase in hindgut-derived pseudoplacenta of a viviparous teleost Xenotoca eiseni.

Nutrient absorption is essential for animal survival and development. Our previous study on zebrafish reported that nutrient absorption in lysosome-rich enterocytes (LREs) is promoted by the voltage-sensing phosphatase (VSP), which regulates phosphoinositide (PIP) homeostasis via electrical signaling in biological membranes. However, it remains unknown whether this VSP function is shared by different absorptive tissues in other species. Here, we focused on the function of VSP in a viviparous teleost Xenotoca eiseni, whose intraovarian embryos absorb nutrients from the maternal ovarian fluid through a specialized hindgut-derived pseudoplacental structure called trophotaenia. Xenotoca eiseni VSP (Xe-VSP) is expressed in trophotaenia epithelium, an absorptive tissue functionally similar to zebrafish LREs. Notably, the apical distribution of Xe-VSP in trophotaenia epithelial cells closely resembles zebrafish VSP (Dr-VSP) distribution in zebrafish LREs, suggesting a shared role for VSP in absorptive tissues between the two species. Electrophysiological analysis using a heterologous expression system revealed that Xe-VSP preserves functional voltage sensors and phosphatase activity with the leftward shifted voltage sensitivity compared with zebrafish VSP (Dr-VSP). We also identified a single amino acid variation in the S4 helix of Xe-VSP as one of the factors contributing to the leftward shifted voltage sensitivity. This study highlights the biological variation and significance of VSP in various animal species, as well as hinting at the potential role of VSP in nutrient absorption in X. eiseni trophotaenia.NEW & NOTEWORTHY We investigate the voltage-sensing phosphatase (VSP) in Xenotoca eiseni, a viviparous fish whose intraovarian embryos utilize trophotaenia for nutrient absorption. Although X. eiseni VSP (Xe-VSP) shares key features with known VSPs, its distinct voltage sensitivity arises from species-specific amino acid variation. Xe-VSP in trophotaenia epithelium suggests its involvement in nutrient absorption, similar to VSP in zebrafish enterocytes and potentially in species with similar absorptive cells. Our findings highlight the potential role of VSP across species.

Deciphering dynamic interactions between spermatozoa and the ovarian microenvironment through integrated multi-omics approaches in viviparous Sebastes schlegelii.

The ovarian microenvironment plays a crucial role in ensuring the reproductive success of viviparous teleosts. However, the molecular mechanism underlying the interaction between spermatozoa and the ovarian microenvironment has remained elusive. This study aimed to contribute to a better understanding of this process in black rockfish (Sebastes schlegelii) using integrated multi-omics approaches. The results demonstrated significant upregulation of ovarian complement-related proteins and pattern recognition receptors, along with remodeling of glycans on the surface of spermatozoa at the early spermatozoa-storage stage (1 month after mating). As spermatozoa were stored over time, ovarian complement proteins were progressively repressed by tryptophan and hippurate, indicating a remarkable adaptation of spermatozoa to the ovarian microenvironment. Before fertilization, a notable upregulation of cellular junction proteins was observed. The study revealed that spermatozoa bind to ZPB2a protein through GSTM3 and that ZPB2a promotes spermatozoa survival and movement in a GSTM3-dependent manner. These findings shed light on a key mechanism that influences the dynamics of spermatozoa in the female reproductive tract, providing valuable insights into the molecular networks regulating spermatozoa adaptation and survival in species with internal fertilization.

Convergent genomic signatures associated with vertebrate viviparity.

BACKGROUND: Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes? RESULTS: We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages. CONCLUSIONS: Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

The genetic basis of a recent transition to live-bearing in marine snails.

Key innovations are fundamental to biological diversification, but their genetic basis is poorly understood. A recent transition from egg-laying to live-bearing in marine snails (Littorina spp.) provides the opportunity to study the genetic architecture of an innovation that has evolved repeatedly across animals. Individuals do not cluster by reproductive mode in a genome-wide phylogeny, but local genealogical analysis revealed numerous small genomic regions where all live-bearers carry the same core haplotype. Candidate regions show evidence for live-bearer-specific positive selection and are enriched for genes that are differentially expressed between egg-laying and live-bearing reproductive systems. Ages of selective sweeps suggest that live-bearer-specific alleles accumulated over more than 200,000 generations. Our results suggest that new functions evolve through the recruitment of many alleles rather than in a single evolutionary step.

Ecological factors and parity mode correlate with genome size variation in squamate reptiles.

BACKGROUND: Evidence of correlation between genome size, the nuclear haploid DNA content of a cell, environmental factors and life-history traits have been reported in many animal species. Genome size, however, spans over three orders of magnitude across taxa and such a correlation does not seem to follow a universal pattern. In squamate reptiles, the second most species-rich order of vertebrates, there are currently no studies investigating drivers of genome size variability. We run a series of phylogenetic generalized least-squares models on 227 species of squamates to test for possible relationships between genome size and ecological factors including latitudinal distribution, bioclimatic variables and microhabitat use. We also tested whether genome size variation can be associated with parity mode, a highly variable life history trait in this order of reptiles. RESULTS: The best-fitting model showed that the interaction between microhabitat use and parity mode mainly accounted for genome size variation. Larger genome sizes were found in live-bearing species that live in rock/sand ecosystems and in egg-laying arboreal taxa. On the other hand, smaller genomes were found in fossorial live-bearing species. CONCLUSIONS: Environmental factors and species parity mode appear to be among the main parameters explaining genome size variation in squamates. Our results suggest that genome size may favour adaptation of some species to certain environments or could otherwise result from the interaction between environmental factors and parity mode. Integration of genome size and genome sequencing data could help understand the role of differential genome content in the evolutionary process of genome size variation in squamates.

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.

Extended embryo retention and viviparity in the first amniotes.

The amniotic egg with its complex fetal membranes was a key innovation in vertebrate evolution that enabled the great diversification of reptiles, birds and mammals. It is debated whether these fetal membranes evolved in eggs on land as an adaptation to the terrestrial environment or to control antagonistic fetal-maternal interaction in association with extended embryo retention (EER). Here we report an oviparous choristodere from the Lower Cretaceous period of northeast China. The ossification sequence of the embryo confirms that choristoderes are basal archosauromorphs. The discovery of oviparity in this assumed viviparous extinct clade, together with existing evidence, suggests that EER was the primitive reproductive mode in basal archosauromorphs. Phylogenetic comparative analyses on extant and extinct amniotes suggest that the first amniote displayed EER (including viviparity).

Placentation in the Mexican scincid lizard Plestiodon brevirostris (Squamata: Scincidae).

Viviparity is the reproductive pattern in which females gestate eggs within their reproductive tract to complete their development and give birth to live offspring. Within extant sauropsids, only the Squamata (e.g., snakes, lizards, and amphisbaenians) evolved viviparity, representing 20% of the existing species. The genus Plestiodon is represented by 43 species and is one of the most widely distributed genera of the Scincidae in Mexico. The goal of this research has been to determine the placental morphology and ontogeny during gestation in the lizard Plestiodon brevirostris. Specimens were dissected to obtain the embryonic chambers and the embryos were categorized to carry out the correlation between the development stage and the placenta development. The embryonic chambers were processed using the conventional histological technique for light microscopy. The identified embryonic stages were 4, 29, 34, 36, and 39. A thin eggshell surrounds the egg in early developmental stages; however, this structure is already absent in the embryonic hemisphere during the developmental stage 29. The results revealed that P. brevirostris is a lecithotrophic species, but a close maternal-fetal relationship is established by tissue apposition. Ontogenically, the placental types that form in the embryonic hemisphere are the chorioplacenta, choriovitelline placenta, and chorioallantoic placenta; whereas the omphaloplacenta is formed in the abembryonic hemisphere. The structure of the chorioallantoic placenta in P. brevirostris suggests that it may play a role during gas exchange between the mother and the embryo, due to the characteristics of the epithelia that comprise it. The structure of embryonic and maternal epithelia of the omphaloplacenta suggests a role in the absorption of the eggshell during gestation and possibly in the transport or diffusion of some nutrients. In general, it is evident that ontogeny and placental characteristics of P. brevirostris match those of other species of viviparous lecithotrophic scincids.

Live birth in a 47-million-year-old snake.

Viviparity is a widespread reproductive trait in snakes, although fossil evidence bearing on its evolution is extremely sparse. Here, we report an exceptional specimen of the minute booid snake Messelophis variatus recovered in the paleolake of the Messel Formation (early-middle Eocene, Germany). This gravid female contains at least two embryos located in the posterior third of the trunk region. The morphology, size, and degree of ossification of the cranial and postcranial remains indicate they correspond with late embryos. This specimen documents the first occurrence of viviparity in a fossil snake and extends the temporal distribution of this reproductive strategy in booid snakes by over 47 Ma. The evolution of viviparity in squamates has traditionally been associated with cold climates, but its presence at the dawn of the evolution of booids during early Palaeogene thermal peaks indicates that viviparity may have evolved under different selective pressures in this clade.

Increased superfetation precedes the evolution of advanced degrees of placentotrophy in viviparous fishes of the family Poeciliidae.

The causes and consequences of the evolution of placentotrophy (post-fertilization nutrition of developing embryos of viviparous organisms by means of a maternal placenta) in non-mammalian vertebrates are still not fully understood. In particular, in the fish family Poeciliidae there is an evolutionary link between placentotrophy and superfetation (ability of females to simultaneously bear embryos at distinct developmental stages), with no conclusive evidence for which of these two traits facilitates the evolution of more advanced degrees of the other. Using a robust phylogenetic comparative method based on Ornstein-Uhlenbeck models of adaptive evolution and data from 36 poeciliid species, we detected a clear causality pattern. The evolution of extensive placentotrophy has been facilitated by the preceding evolution of more simultaneous broods. Therefore, placentas became increasingly complex as an adaptive response to evolutionary increases in the degree of superfetation. This finding represents a substantial contribution to our knowledge of the factors that have shaped placental evolution in poeciliid fishes.

The eggshell-matrix protein gene OC-17 is functionally lost in the viviparous Chinese crocodile lizard.

Thickness reduction or loss of the calcareous eggshell is one of major phenotypic changes in the transition from oviparity to viviparity. Whether the reduction of eggshells in viviparous squamates is associated with specific gene losses is unknown. Taking advantage of a newly generated high-quality genome of the viviparous Chinese crocodile lizard (Shinisaurus crocodilurus), we found that ovocleidin-17 gene (OC-17), which encodes an eggshell matrix protein that is essential for calcium deposition in eggshells, is not intact in the crocodile lizard genome. Only OC-17 transcript fragments were found in the oviduct transcriptome, and no OC-17 peptides were identified in the eggshell proteome of crocodile lizards. In contrast, OC-17 was present in the eggshells of the oviparous Mongolia racerunner (Eremias argus). Although the loss of OC-17 is not common in viviparous species, viviparous squamates show fewer intact eggshell-specific proteins than oviparous squamates. Our study implies that functional loss of eggshell-matrix protein genes may be involved in the reduction of eggshells during the transition from oviparity to viviparity in the crocodile lizard.

Exceptional parallelisms characterize the evolutionary transition to live birth in phrynosomatid lizards.

Viviparity, an innovation enhancing maternal control over developing embryos, has evolved >150 times in vertebrates, and has been proposed as an adaptation to inhabit cold habitats. Yet, the behavioral, physiological, morphological, and life history features associated with live-bearing remain unclear. Here, we capitalize on repeated origins of viviparity in phrynosomatid lizards to tease apart the phenotypic patterns associated with this innovation. Using data from 125 species and phylogenetic approaches, we find that viviparous phrynosomatids repeatedly evolved a more cool-adjusted thermal physiology than their oviparous relatives. Through precise thermoregulatory behavior viviparous phrynosomatids are cool-adjusted even in warm environments, and oviparous phrynosomatids warm-adjusted even in cool environments. Convergent behavioral shifts in viviparous species reduce energetic demand during activity, which may help offset the costs of protracted gestation. Whereas dam and offspring body size are similar among both parity modes, annual fecundity repeatedly decreases in viviparous lineages. Thus, viviparity is associated with a lower energetic allocation into production. Together, our results indicate that oviparity and viviparity are on opposing ends of the fast-slow life history continuum in both warm and cool environments. In this sense, the 'cold climate hypothesis' fits into a broader range of energetic/life history trade-offs that influence transitions to viviparity.

Ecological and life-history correlates of erythrocyte size and shape in Lepidosauria.

Blood oxygen-carrying capacity is shaped both by the ambient oxygen availability as well as species-specific oxygen demand. Erythrocytes are a critical part of oxygen transport and both their size and shape can change in relation to species-specific life-history, behavioural or ecological conditions. Here, we test whether components of the environment (altitude), life history (reproductive mode, body temperature) and behaviour (diving, foraging mode) drive erythrocyte size variation in the Lepidosauria (lizards, snakes and rhynchocephalians). We collected data on erythrocyte size (area) and shape (L/W: elongation ratio) from Lepidosauria across the globe (N = 235 species). Our analyses show the importance of oxygen requirements as a driver of erythrocyte size. Smaller erythrocytes were associated with the need for faster delivery (active foragers, high-altitude species, warmer body temperatures), whereas species with greater oxygen demands (diving species, viviparous species) had larger erythrocytes. Erythrocyte size shows considerable cross-species variation, with a range of factors linked to the oxygen delivery requirements being major drivers of these differences. A key future aspect for study would include within-individual plasticity and how changing states, for example, pregnancy, perhaps alter the size and shape of erythrocytes in Lepidosaurs.

Different Genes are Recruited During Convergent Evolution of Pregnancy and the Placenta.

The repeated evolution of the same traits in distantly related groups (convergent evolution) raises a key question in evolutionary biology: do the same genes underpin convergent phenotypes? Here, we explore one such trait, viviparity (live birth), which, qualitative studies suggest, may indeed have evolved via genetic convergence. There are >150 independent origins of live birth in vertebrates, providing a uniquely powerful system to test the mechanisms underpinning convergence in morphology, physiology, and/or gene recruitment during pregnancy. We compared transcriptomic data from eight vertebrates (lizards, mammals, sharks) that gestate embryos within the uterus. Since many previous studies detected qualitative similarities in gene use during independent origins of pregnancy, we expected to find significant overlap in gene use in viviparous taxa. However, we found no more overlap in uterine gene expression associated with viviparity than we would expect by chance alone. Each viviparous lineage exhibits the same core set of uterine physiological functions. Yet, contrary to prevailing assumptions about this trait, we find that none of the same genes are differentially expressed in all viviparous lineages, or even in all viviparous amniote lineages. Therefore, across distantly related vertebrates, different genes have been recruited to support the morphological and physiological changes required for successful pregnancy. We conclude that redundancies in gene function have enabled the repeated evolution of viviparity through recruitment of different genes from genomic "toolboxes", which are uniquely constrained by the ancestries of each lineage.

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.

Independent evolutionary transitions to pueriparity across multiple timescales in the viviparous genus Salamandra.

The ability to bear live offspring, viviparity, has evolved multiple times across the tree of life and is a remarkable adaptation with profound life-history and ecological implications. Within amphibians the ancestral reproductive mode is oviparity followed by a larval life stage, but viviparity has evolved independently in all three amphibian orders. Two types of viviparous reproduction can be distinguished in amphibians; larviparity and pueriparity. Larviparous amphibians deliver larvae into nearby ponds and streams, while pueriparous amphibians deliver fully developed juveniles and thus do not require waterbodies for reproduction. Among amphibians, the salamander genus Salamandra is remarkable as it exhibits both inter- and intraspecific variation in the occurrence of larviparity and pueriparity. While the evolutionary relationships among Salamandra lineages have been the focus of several recent studies, our understanding of how often and when transitions between modes occurred is still incomplete. Furthermore, in species with intraspecific variation, the reproductive mode of a given population can only be confirmed by direct observation of births and thus the prevalence of pueriparous populations is also incompletely documented. We used sequence capture to obtain 1,326 loci from 94 individuals from across the geographic range of the genus, focusing on potential reproductive mode transition zones. We also report additional direct observations of pueriparous births for 20 new locations and multiple lineages. We identify at least five independent transitions from the ancestral mode of larviparity to pueriparity among and within species, occurring at different evolutionary timescales ranging from the Pliocene to the Holocene. Four of these transitions occurred within species. Based on a distinct set of markers and analyses, we also confirm previous findings of introgression between species and the need for taxonomic revisions in the genus. We discuss the implications of our findings with respect to the evolution of this complex trait, and the potential of using five independent convergent transitions for further studies on the ecological context in which pueriparity evolves and the genetic architecture of this specialized reproductive mode.

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.

The Evolution of Viviparity in Vertebrates.

In the vertebrate tree of life, viviparity or live birth has independently evolved many times, resulting in a rich diversity of reproductive strategies. Viviparity is believed to be a mode of reproduction that evolved from the ancestral condition of oviparity or egg laying, where most of the fetal development occurs outside the body. Today, there is not a simple model of parity transition to explain this species-specific divergence in modes of reproduction. Most evidence points to a gradual series of evolutionary adaptations that account for this phenomenon of reproduction, elegantly displayed by various viviparous squamates that exhibit placentae formed by the appositions of maternal and embryonic tissues, which share significant homology with the tissues that form the placenta in therian mammals. In an era where the genomes of many vertebrate species are becoming available, studies are now exploring the molecular basis of this transition from oviparity to viviparity, and in some rare instances its possible reversibility, such as the Australian three-toed skink (Saiphos equalis). In contrast to the parity diversity in squamates, mammals are viviparous with the notable exception of the egg-laying monotremes. Advancing computational tools coupled with increasing genome availability across species that utilize different reproductive strategies promise to reveal the molecular underpinnings of the ancestral transition of oviparity to viviparity. As a result, the dramatic changes in reproductive physiology and anatomy that accompany these parity changes can be reinterpreted. This chapter will briefly explore the vertebrate modes of reproduction using a phylogenetic framework and where possible highlight the role of potential candidate genes that may help explain the polygenic origins of live birth.

A staging table of embryonic development for a viviparous (live-bearing) lizard Eremias multiocellata (Squamata: Lacertidae).

As the only viviparous reptile in China that has both temperature-dependent sex determination (TSD) and genetic-dependent sex determination (GSD) mechanisms, Eremias multiocellata is considered as an ideal species for studying the sex determination mechanism in viviparous lizards. However, studies on embryonic stage of viviparous lizards and morphological characteristics of each stage are limited. In the present study, the embryonic development process of E. multiocellata is divided into 15 stages (stages 28-42) according to the morphology of embryos. Embryos sizes are measured and continuous dynamic variation of some key features, including limbs, genitals, eyes, pigments, and brain scales are color imaged by a stereoscopic microscope. Furthermore, based on these morphological characteristics, we compare the similarities and differences in the embryonic development of E. multiocellata with other squamate species. Our results not only identified the staging table of E. multiocellata with continuous changes of external morphological characteristics but also developed a staging scheme for an important model species that provides a necessary foundation for study of sex determination in a viviparous lizard.