Genomic and transcriptomic investigations of the evolutionary transition from oviparity to viviparity

A synthetic review: natural history of amniote reproductive modes in light of comparative evolutionary genomics

There is a current lack of consensus on whether the ancestral parity mode was oviparity (egg-laying) or viviparity (live-birth) in amniotes and particularly in squamates (snakes, lizards, and amphisbaenids). How transitions between parity modes occur at the genomic level has primary importance for how science conceptualises the origin of amniotes, and highly variable parity modes in Squamata. Synthesising literature from medicine, poultry science, reproductive biology, and evolutionary biology, I review the genomics and physiology of five broad processes (here termed the 'Main Five') expected to change during transitions between parity modes: eggshell formation, embryonic retention, placentation, calcium transport, and maternal-fetal immune dynamics. Throughout, I offer alternative perspectives and testable hypotheses regarding proximate causes of parity mode evolution in amniotes and squamates. If viviparity did evolve early in the history of lepidosaurs, I offer the nucleation site hypothesis as a proximate explanation. The framework of this hypothesis can be extended to amniotes to infer their ancestral state. I also provide a mechanism and hypothesis on how squamates may transition from viviparity to oviparity and make predictions about the directionality of transitions in three species. After considering evidence for differing perspectives on amniote origins, I offer a framework that unifies (i) the extended embryonic retention model and (ii) the traditional model which describes the amniote egg as an adaptation to the terrestrial environment. Additionally, this review contextualises the origin of amniotes and parity mode evolution within Medawar's paradigm. Medawar posited that pregnancy could be supported by immunosuppression, inertness, evasion, or immunological barriers. I demonstrate that this does not support gestation or gravidity across most amniotes but may be an adequate paradigm to explain how the first amniote tolerated internal fertilization and delayed egg deposition. In this context, the eggshell can be thought of as an immunological barrier. If serving as a barrier underpins the origin of the amniote eggshell, there should be evidence that oviparous gravidity can be met with a lack of immunological responses in utero. Rare examples of two species that differentially express very few genes during gravidity, suggestive of an absent immunological reaction to oviparous gravidity, are two skinks Lampropholis guichenoti and Lerista bougainvillii. These species may serve as good models for the original amniote egg. Overall, this review grounds itself in the historical literature while offering a modern perspective on the origin of amniotes. I encourage the scientific community to utilise this review as a resource in evolutionary and comparative genomics studies, embrace the complexity of the system, and thoughtfully consider the frameworks proposed.

Common lizard primary oviduct cell culture: A model system for the genetic and cellular basis of oviparity and viviparity

Reproduction by egg-laying (oviparity) or live-bearing (viviparity) is a genetically determined trait fundamental to the biology of amniotes. Squamates are an emerging model for the genetics of reproductive mode yet lack cell culture models valuable for exploring molecular mechanisms. Here, we report a novel primary culture model for reproductive biology: cell cultures derived from the oviduct tissues (infundibulum, uterus and vagina) of oviparous and viviparous common lizards (Lacertidae: Zootoca vivipara). We maintained and expanded these cultures for over 100 days, including repeated subculturing and successful revival of cryopreserved cells. Immunocytochemical investigation suggested expression of both epithelial and fibroblast-like proteins, and RNA sequencing of cultured cells as compared to in vivo oviduct tissue showed changes in gene expression in response to the cell culture environment. Despite this, we confirmed the maintenance of distinct gene expression patterns in viviparous and oviparous cells after 60+ days of cell culture, finding 354 differentially expressed genes between viviparous and oviparous cells. Furthermore, we confirmed the expression of 15 viviparity-associated candidate genes in cells maintained for 60+ days in culture. Our study demonstrates the feasibility and utility of oviduct cell culture for molecular analysis of reproductive mode and provides a tool for future genetic experiments.

Transcriptomic profiling revealed immune-related signaling pathways in response to experimental infection of Leishmania donovani in two desert lizards from Northwest China

Little is known about the immune response of lizards to Leishmania parasties. In this study, we conducted the first liver transcriptome analysis of two lizards (Phrynocephalus przewalskii and Eremias multiocellata) challenged with L. donovani, endemic to the steppe desert region of northwestern China. Our results revealed that multiple biological processes and immune-related signaling pathways are closely associated with the immune response to experimental L. donovani infection in the two lizards, and that both lizards show similar changes to mammals in terms of immunity to Leishmania. However, the interspecific divergence of the two lizards leads to different transcriptomic changes. In particular, in contrast to P. przewalskii, the challenged E. mutltiocellata was characterized by the induction of down-regulation of most DEGs. These findings will contribute to the scarce resources on lizard immunity and provide a reference for further research on immune mechanisms in reptiles.

Genetically Encoded Lizard Color Divergence for Camouflage and Thermoregulation

Local adaptation is critical in speciation and evolution, yet comprehensive studies on proximate and ultimate causes of local adaptation are generally scarce. Here, we integrated field ecological experiments, genome sequencing, and genetic verification to demonstrate both driving forces and molecular mechanisms governing local adaptation of body coloration in a lizard from the Qinghai-Tibet Plateau. We found dark lizards from the cold meadow population had lower spectrum reflectance but higher melanin contents than light counterparts from the warm dune population. Additionally, the colorations of both dark and light lizards facilitated the camouflage and thermoregulation in their respective microhabitat simultaneously. More importantly, by genome resequencing analysis, we detected a novel mutation in Tyrp1 that underpinned this color adaptation. The allele frequencies at the site of SNP 459# in the gene of Tyrp1 are 22.22% G/C and 77.78% C/C in dark lizards and 100% G/G in light lizards. Model-predicted structure and catalytic activity showed that this mutation increased structure flexibility and catalytic activity in enzyme TYRP1, and thereby facilitated the generation of eumelanin in dark lizards. The function of the mutation in Tyrp1 was further verified by more melanin contents and darker coloration detected in the zebrafish injected with the genotype of Tyrp1 from dark lizards. Therefore, our study demonstrates that a novel mutation of a major melanin-generating gene underpins skin color variation co-selected by camouflage and thermoregulation in a lizard. The resulting strong selection may reinforce adaptive genetic divergence and enable the persistence of adjacent populations with distinct body coloration.

Evolutionary traits and functional roles of chemokines and their receptors in the male pregnancy of the Syngnathidae

Vertebrates have developed various modes of reproduction, some of which are found in Teleosts. Over 300 species of the Syngnathidae (seahorses, pipefishes and seadragons) exhibit male pregnancies; the males have specialized brood pouches that provide immune protection, nourishment, and oxygen regulation. Chemokines play a vital role at the mammalian maternal-fetal interface; however, their functions in fish reproduction are unclear. This study revealed the evolutionary traits and potential functions of chemokine genes in 22 oviparous, ovoviviparous, and viviparous fish species through comparative genomic analyses. Our results showed that chemokine gene copy numbers and evolutionary rates vary among species with different modes of reproduction. Syngnathidae lost cxcl13 and cxcr5, which are involved in key receptor-ligand pairs for lymphoid organ development. Notably, Syngnathidae have site-specific mutations in cxcl12b and ccl44, suggesting immune function during gestation. Moreover, transcriptome analysis revealed that chemokine gene expression varies among Syngnathidae species with different types of brood pouches, suggesting adaptive variations in chemokine functions among seahorses and their relatives. Furthermore, challenge experiments on seahorse brood pouches revealed a joint immune function of chemokine genes during male pregnancy. This study provides insights into the evolutionary diversity of chemokine genes associated with different reproductive modes in fish.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00205-x.

Live-bearing cockroach genome reveals convergent evolutionary mechanisms linked to viviparity in insects and beyond

Live birth (viviparity) has arisen repeatedly and independently among animals. We sequenced the genome and transcriptome of the viviparous Pacific beetle-mimic cockroach and performed comparative analyses with two other viviparous insect lineages, tsetse flies and aphids, to unravel the basis underlying the transition to viviparity in insects. We identified pathways undergoing adaptive evolution for insects, involved in urogenital remodeling, tracheal system, heart development, and nutrient metabolism. Transcriptomic analysis of cockroach and tsetse flies revealed that uterine remodeling and nutrient production are increased and the immune response is altered during pregnancy, facilitating structural and physiological changes to accommodate and nourish the progeny. These patterns of convergent evolution of viviparity among insects, together with similar adaptive mechanisms identified among vertebrates, highlight that the transition to viviparity requires changes in urogenital remodeling, enhanced tracheal and heart development (corresponding to angiogenesis in vertebrates), altered nutrient metabolism, and shifted immunity in animal systems.

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).

Alternative signal pathways underly fertilization and egg activation in a fish with contrasting modes of spawning

BACKGROUND

The processes of fertilization and egg activation are vital for early embryogenesis. However, while the mechanisms associated with key events during these processes differ among species and modes of spawning, the signal pathways underlying these processes are opaque for many fishes, including economically important species.

RESULTS

We investigated phenotypic traits, ultrastructure and protein expression levels in the eggs of the topmouth culter (Culter alburnus), a protected and economically important freshwater fish that exhibits two spawning modes, producing semi-buoyant eggs and adhesive eggs. Unfertilized eggs of C. alburnus were examined, as well as eggs at fertilization and 30 min post fertilization. Our results showed that in semi-buoyant eggs, energy metabolism was activated at fertilization, followed by elevated protein expression of cytoskeleton and extracellular matrix (ECM)-receptor interactions that resulted in rapid egg swelling; a recognized adaptation for lotic habitats. In contrast, in adhesive eggs fertilization initiated the process of sperm-egg fusion and blocking of polyspermy, followed by enhanced protein expression of lipid metabolism and the formation of egg envelope adhesion and hardening, which are adaptive in lentic habitats.

CONCLUSION

Our findings indicate that alternative signal pathways differ between modes of spawning and timing during the key processes of fertilization and egg activation, providing new insights into the molecular mechanisms involved in adaptive early embryonic development in teleost fishes.

A chromosome-level genome assembly of Ostrea denselamellosa provides initial insights into its evolution

The oyster Ostrea denselamellosa is a live-bearing species with a sharp decline in the natural population. Despite recent breakthroughs in long-read sequencing, high quality genomic data are very limited in O. denselamellosa. Here, we carried out the first whole genome sequencing at the chromosome-level in O. denselamellosa. Our studies yielded a 636 Mb assembly with scaffold N50 around 71.80 Mb. 608.3 Mb (95.6% of the assembly) were anchored to 10 chromosomes. A total of 26,412 protein-coding genes were predicted, of which 22,636 (85.7%) were functionally annotated. By comparative genomics, we found that long interspersed nuclear element (LINE) and short interspersed nuclear element (SINE) made up a larger proportion in O. denselamellosa genome than in other oysters'. Moreover, gene family analysis showed some initial insight into its evolution. This high-quality genome of O. denselamellosa provides a valuable genomic resource for studies of evolution, adaption and conservation in oysters.

Comparing the potential for maternal-fetal signalling in oviparous and viviparous lizards

The evolution of a placenta requires several steps including changing the timing of reproductive events, facilitating nutrient exchange, and the capacity for maternal-fetal communication. To understand the evolution of maternal-fetal communication, we used ligand-receptor gene expression as a proxy for the potential for cross-talk in a live-bearing lizard (Pseudemoia entrecasteauxii) and homologous tissues in a related egg-laying lizard (Lampropholis guichenoti). Approximately 70% of expressed ligand/receptor genes were shared by both species. Gene ontology (GO) analysis showed that there was no GO-enrichment in the fetal membranes of the egg-laying species, but live-bearing fetal tissues were significantly enriched for 50 GO-terms. Differences in enrichment suggest that the evolution of viviparity involved reinforcing specific signalling pathways, perhaps to support fetal control of placentation. One identified change was in transforming growth factor beta signalling. Using immunohistochemistry, we show the production of the signalling molecule inhibin beta B (INHBB) occurs in viviparous fetal membranes but was absent in closely related egg-laying tissues, suggesting that the evolution of viviparity may have involved changes to signalling via this pathway. We argue that maternal-fetal signalling evolved through co-opting expressed signalling molecules and recruiting new signalling molecules to support the complex developmental changes required to support a fetus in utero. This article is part of the theme issue 'Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom'.

Embryonic specializations for vertebrate placentation

The vertebrate placenta, a close association of fetal and parental tissue for physiological exchange, has evolved independently in sharks, teleost fishes, coelacanths, amphibians, squamate reptiles and mammals. This transient organ forms during pregnancy and is an important contributor to embryonic development in both viviparous and oviparous, brooding species. Placentae may be involved in transport of respiratory gases, wastes, immune molecules, hormones and nutrients. Depending on the taxon, the embryonic portion of the placenta is comprised of either extraembryonic membranes (yolk sac or chorioallantois) or temporary embryonic tissues derived via hypertrophy of pericardium, gill epithelium, gut, tails or fins. These membranes and tissues have been recruited convergently into placentae in several lineages. Here, we highlight the diversity and common features of embryonic tissues involved in vertebrate placentation and suggest future studies that will provide new knowledge about the evolution of pregnancy. This article is part of the theme issue 'Extraembryonic tissues: exploring concepts, definitions and functions across the animal kingdom'.

The evolution of reproductive modes and life cycles in amphibians

Amphibians have undergone important evolutionary transitions in reproductive modes and life-cycles. We compare large-scale macroevolutionary patterns in these transitions across the three major amphibian clades: frogs, salamanders, and caecilians. We analyse matching reproductive and phylogenetic data for 4025 species. We find that having aquatic larvae is ancestral for all three groups and is retained by many extant species (33-44%). The most frequent transitions in each group are to relatively uncommon states: live-bearing in caecilians, paedomorphosis in salamanders, and semi-terrestriality in frogs. All three groups show transitions to more terrestrial reproductive modes, but only in caecilians have these evolved sequentially from most-to-least aquatic. Diversification rates are largely independent of reproductive modes. However, in salamanders direct development accelerates diversification whereas paedomorphosis decreases it. Overall, we find a widespread retention of ancestral modes, decoupling of trait transition rates from patterns of species richness, and the general independence of reproductive modes and diversification.

Gut microbiota differs between two cold-climate lizards distributed in thermally different regions

BACKGROUND

The metabolic cold-climate adaption hypothesis predicts that animals from cold environments have relatively high metabolic rates compared with their warm-climate counterparts. However, studies testing this hypothesis are sparse. Here, we compared gut microbes between two cold-climate lizard species of the genus Phrynocephalus to see if gut microbiota could help lizards adapt to cold environments by promoting metabolism. We conducted a 2 species (P. erythrurus and P. przewalskii) × 2 temperatures (24 and 30 °C) factorial design experiment, whereby we kept lizards of two Phrynocephalus species at 24 and 30 °C for 25 d and then collected their fecal samples to analyze and compare the microbiota based on 16S rRNA gene sequencing technology.

RESULTS

The gut microbiota was mainly composed of bacteria of the phyla Proteobacteria, Firmicutes, Bacteroidetes, and Verrucomicrobia in both species (Proteobacteria > Firmicutes > Verrucomicrobiota in P. erythrurus, and Bacteroidetes > Proteobacteria > Firmicutes in P. przewalskii). Further analysis revealed that the gut microbiota promoted thermal adaptation in both lizard species, but with differences in the relative abundance of the contributory bacteria between the two species. An analysis based on the Kyoto Encyclopedia of Genes and Genomes revealed that the gut microbiota played important roles in metabolism, genetic information processing, cellular processes, and environmental information processing in both species. Furthermore, genes related to metabolism were more abundant in P. erythrurus at 24 °C than in other species ⋅ temperature combinations.

CONCLUSION

Our study provides evidence that gut microbiota promotes thermal adaptation in both species but more evidently in P. erythrurus using colder habitats than P. przewalskii all year round, thus confirming the role of gut microbiota in cold-climate adaptation in lizards.

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.

Molecular Mechanisms of Spawning Habits for the Adaptive Radiation of Endemic East Asian Cyprinid Fishes

Despite the widespread recognition of adaptive radiation as a driver of speciation, the mechanisms by which natural selection generates new species are incompletely understood. The evolutionary radiation of endemic East Asian cyprinids has been proposed as evolving through a change in spawning habits, involving a transition from semibuoyant eggs to adhesive eggs in response to crosslinked river-lake system formation. Here, we investigated the molecular mechanisms that underpin this radiation, associated with egg hydration and adhesiveness. We demonstrated that semibuoyant eggs enhance hydration by increasing the degradation of yolk protein and accumulation of Ca2+ and Mg2+ ions, while adhesive eggs improve adhesiveness and hardness of the egg envelope by producing an adhesive layer and a unique 4th layer to the egg envelope. Based on multiomics analyses and verification tests, we showed that during the process of adaptive radiation, adhesive eggs downregulated the "vitellogenin degradation pathway," "zinc metalloprotease pathway," and "ubiquitin-proteasome pathway" and the pathways of Ca2+ and Mg2+ active transport to reduce their hydration. At the same time, adhesive eggs upregulated the crosslinks of microfilament-associated proteins and adhesive-related proteins, the hardening-related proteins of the egg envelope, and the biosynthesis of glycosaminoglycan in the ovary to generate adhesiveness. These findings illustrate the novel molecular mechanisms associated with hydration and adhesiveness of freshwater fish eggs and identify critical molecular mechanisms involved in the adaptive radiation of endemic East Asian cyprinids. We propose that these key egg attributes may function as "magic traits" in this adaptive radiation.

Transcriptomic Insight into Viviparous Growth in Water Lily

Water lily is an important ornamental flower plant which is capable of viviparous plantlet development. But no study has been reported on the molecular basis of viviparity in water lily. Hence, we performed a comparative transcriptome study between viviparous water lily Nymphaea micrantha and a nonviviparous species Nymphaea colorata at four developmental stages. The higher expression of highly conserved AUX/IAA, ARF, GH3, and SAUR gene families in N. micrantha compared to N. colorata is predicted to have a major impact on the development and evolution of viviparity in water lily. Likewise, differential regulation of hormone signaling, brassinosteroid, photosynthesis, and energy-related pathways in the two species provide clues of their involvement in viviparity phenomenon. This study revealed the complex mechanism of viviparity trait in water lily. The transcriptomic signatures identified are important basis for future breeding and research of viviparity in water lily and other plant species.

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.

Genetic basis for the evolution of pelvic-fin brooding, a new mode of reproduction, in a Sulawesian fish

Modes of reproduction in animals are diverse, with different modes having evolved independently in multiple lineages across a variety of taxa. However, an understanding of the genomic change driving the transition between different modes of reproduction is limited. Several ricefishes (Adrianichthyidae) on the island of Sulawesi have a unique mode of reproduction called "pelvic-fin brooding," wherein females carry externally fertilized eggs until hatching using their pelvic fins. Phylogenomic analysis demonstrated pelvic-fin brooders to have evolved at least twice in two distant clades of the Adrianichthyidae. We investigated the genetic architecture of the evolution of this unique mode of reproduction. Morphological analyses and laboratory observations revealed that females of pelvic-fin brooders have longer pelvic fins and a deeper abdominal concavity, and that they can carry an egg clutch for longer than non-brooding adrianichthyids, suggesting that these traits play important roles in this reproductive mode. Quantitative trait locus mapping using a cross between a pelvic-fin brooder Oryzias eversi and a non-brooding O. dopingdopingensis reveals different traits involved in pelvic-fin brooding to be controlled by different loci on different chromosomes. Genomic analyses of admixture detected no signatures of introgression between two lineages with pelvic-fin brooders, indicating that introgression is unlikely to be responsible for repeated evolution of pelvic-fin brooding. These findings suggest that multiple independent mutations may have contributed to the convergent evolution of this novel mode of reproduction.

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.

Species persistence with hybridization in toad-headed lizards driven by divergent selection and low recombination

Speciation plays a central role in evolutionary studies, and particularly how reproductive isolation (RI) evolves. The origins and persistence of RI are distinct processes that require separate evaluations. Treating them separately clarifies the drivers of speciation and then it is possible to link the processes to understand large-scale patterns of diversity. Recent genomic studies have focused predominantly on how species or RI originate. However, we know little about how species persist in face of gene flow. Here, we evaluate a contact zone of two closely related toad-headed lizards (Phrynocephalus) using a chromosome-level genome assembly and population genomics. To some extent, recent asymmetric introgression from Phrynocephalus putjatai to P. vlangalii reduces their genomic differences. However, their highly divergent regions (HDRs) have heterogeneous distributions across the genomes. Functional gene annotation indicates that many genes within HDRs are involved in reproduction and RI. Compared with allopatric populations, contact areas exhibit recent divergent selection on the HDRs and a lower population recombination rate. Taken together, this implies that divergent selection and low genetic recombination help maintain RI. This study provides insights into the genomic mechanisms that drive RI and two species persistence in the face of gene flow during the late stage of speciation.

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.

Comparative Analysis Among Different Species Reveals That the Androgen Receptor Regulates Chicken Follicle Selection Through Species-Specific Genes Related to Follicle Development

The differences in reproductive processes at the molecular level between viviparous and oviparous animals are evident, and the site in the ovary that synthesizes sex hormones (androgens and oestrogens) and the trends for enriching sex hormones during follicle development in chickens are different from those in mammals, suggesting that the effect of sex hormones on follicle development in chickens is probably different from that in viviparous animals. To explore the specific role of androgen receptors (ARs) on chicken follicular development, we matched the correspondence of follicular development stages among chickens, humans, cows and identified chicken-specific genes related to follicle development (GAL-SPGs) by comparing follicle development-related genes and their biological functions among species (chickens, humans, and cows). A comparison of the core transcription factor regulatory network of granulosa cells (or ovaries) based on super-enhancers among species (chicken, human, and mouse) revealed that AR is a core transcriptional regulator specific to chickens. In vivo experiments showed that inhibition of AR significantly reduced the number of syf (selected stage follicles) in chickens and decreased the expression of GAL-SPGs in F5 follicles, while in vitro experiments showed that inhibition of AR expression in chicken granulosa cells (GCs) significantly down-regulated the expression levels of GAL-SPGs, indicating that AR could regulate follicle selection through chicken-specific genes related to follicle development. A comparison among species (77 vertebrates) of the conserved genomic regions, where chicken super-enhancers are located, revealed that the chicken AR super-enhancer region is conserved in birds, suggesting that the role of AR in follicle selection maybe widespread in birds. In summary, we found that AR can regulate follicle selection through chicken-specific genes related to follicle development, which also emphasizes the important role of AR in follicle selection in chickens and provides a new perspective for understanding the unique process of follicle development in chickens. Our study will contribute to the application of androgens to the control of egg production in chickens and suggests that researchers can delve into the mechanisms of follicle development in birds based on androgen/androgen receptors.

Geological and climatic influences on population differentiation of the Phrynocephalus vlangalii species complex (Sauria: Agamidae) in the northern Qinghai-Tibet Plateau

Extremely heterogeneous topography and complex paleoclimate history of the Qinghai-Tibet Plateau (QTP) have a key role in promoting genetic divergence and lineage/species formation. Here, we sequenced one nuclear and three mitochondrial markers of 532 individuals from the entire range of the Phrynocephalus vlangalii species complex including two species, P. putjatai and P. vlangalii, endemic to the northern QTP. We integrated multilocus phylogeny, demographic analysis and geographic barrier detection to evaluate the population structure and dynamics. We found a new mitochondrial clade (PV-I) in the Gonghe County population of P. vlangalii, partial mitochondrial DNA replacement within P. vlangalii and complete mitochondrial DNA replacement between P. putjatai and P. vlangalii. Neutrality test, mismatch distribution analysis and Extended Bayesian Skyline Plot (EBSP) analysis all supported a significant expansion of the Qaidam Basin population of P. vlangalii (PV-II-2) from 0.091 to 0.026 Ma after Penultimate Glaciation. The uplift of the Arjin and Anyemanqen Mountains during the Kunhuang Movement (∼1.2 Ma) split populations of P. vlangalii in Akesai, Qaidam Basin and source of the Yellow River. The uplift of the Elashan Mountains during the second phase of the Qingzang Movement (∼2.5 Ma) contributed to the divergence of the Gonghe County population of P. vlangalii from other conspecific populations. The third phase of the Qingzang Movement (∼1.7 Ma) contributed to the divergence of the Xinghai population of P. vlangalii from P. putjatai and to the divergence of the northern populations of P. putjatai from the southern conspecific populations. Our data support the idea that the geological and climatic changes following the orogeny of the QTP may have promoted population differentiation and shaped the current population patterns of the P. vlangalii species complex in the northeastern QTP.

Sex Chromosomes and Master Sex-Determining Genes in Turtles and Other Reptiles

Among tetrapods, the well differentiated heteromorphic sex chromosomes of birds and mammals have been highly investigated and their master sex-determining (MSD) gene, Dmrt1 and SRY, respectively, have been identified. The homomorphic sex chromosomes of reptiles have been the least studied, but the gap with birds and mammals has begun to fill. This review describes our current knowledge of reptilian sex chromosomes at the cytogenetic and molecular level. Most of it arose recently from various studies comparing male to female gene content. This includes restriction site-associated DNA sequencing (RAD-Seq) experiments in several male and female samples, RNA sequencing and identification of Z- or X-linked genes by male/female comparative transcriptome coverage, and male/female transcriptomic or transcriptome/genome substraction approaches allowing the identification of Y- or W-linked transcripts. A few putative master sex-determining (MSD) genes have been proposed, but none has been demonstrated yet. Lastly, future directions in the field of reptilian sex chromosomes and their MSD gene studies are considered.

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.

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.

Structural changes to the brood pouch of male pregnant seahorses (Hippocampus abdominalis) facilitate exchange between father and embryos

INTRODUCTION

Embryonic growth and development require efficient respiratory gas exchange. Internal incubation of developing young thus presents a significant physiological challenge, because respiratory gas diffusion to embryos is impeded by the additional barrier of parental tissue between the embryo and the environment. Therefore, live-bearing species exhibit a variety of adaptations facilitating respiratory gas exchange between the parent (usually the mother) and embryos. Syngnathid fishes are the only vertebrates to exhibit male pregnancy, allowing comparative studies of the biology and evolution of internal incubation of embryos, independent of the female reproductive tract. Here, we examine the fleshy, sealed, seahorse brood pouch, and provide the first quantification of structural changes to this gestational organ across pregnancy.

METHODS

We used histological analysis and morphometrics to quantify the surface area for exchange across the brood pouch epithelium, and the structure of the vascular bed of the brood pouch.

RESULTS

We show dramatic remodelling of gestational tissues as pregnancy progresses, including an increase in tortuosity of the gestational epithelium, an increase in capillary density, and a decrease in diffusion distance between capillaries and the pouch lumen.

DISCUSSION

These changes produce an increased surface area and expansion of the vascular bed of the placenta that likely facilitates respiratory gas exchange. These changes mirror the remodelling of gestational tissue in viviparous amniotes and elasmobranchs, and provide further evidence of the convergence of adaptations to support pregnancy in live-bearing animals.

Parallel evolution of placental calcium transfer in the lizard Mabuya and eutherian mammals

An exceptional case of parallel evolution between lizards and eutherian mammals occurs in the evolution of viviparity. In the lizard genus Mabuya, viviparity provided the environment for the evolution of yolk-reduced eggs and obligate placentotrophy. One major event that favored the evolution of placentation was the reduction of the eggshell. As with all oviparous reptiles, lizard embryos obtain calcium from both the eggshell and egg yolk. Therefore, the loss of the eggshell likely imposes a constraint for the conservation of the egg yolk, which can only be obviated by the evolution of alternative mechanisms for the transport of calcium directly from the mother. The molecular and cellular mechanisms employed to solve these constraints, in a lizard with only a rudimentary eggshell such as Mabuya, are poorly understood. Here, we used RT-qPCR on placental and uterine samples during different stages of gestation in Mabuya, and demonstrate that transcripts of the calcium transporters trpv6, cabp28k, cabp9k and pmca are expressed and gradually increase in abundance through pregnancy stages, reaching their maximum expression when bone mineralization occurs. Furthermore, CABP28K/9K proteins were studied by immunofluorescence, demonstrating expression in specific regions of the mature placenta. Our results indicate that the machinery for calcium transportation in the Mabuya placenta was co-opted from other tissues elsewhere in the vertebrate bodyplan. Thus, the calcium transportation machinery in the placenta of Mabuya evolved in parallel with the mammalian placenta by redeploying the expression of similar calcium transporter genes.

Effect of early thermal environment on the morphology and performance of a lizard species with bimodal reproduction

Early developmental environments influence life-history traits and survival of reptiles. In fact, rising global temperatures have already caused widespread extinctions among lizards. Viviparous species might be more susceptible to increasing temperatures because of their inability to meet their energetic demands following rapid environmental changes. Nevertheless, we do not yet fully understand how lizards with different reproductive modes can respond to climate change. We investigated the effect of both maternal thermal environment during pregnancy and incubation temperature on hatchling morphology and physiological performance of two populations of the lizard Saiphos equalis differing in their mode of reproduction, to test whether reproductive mode affects the ability to buffer against rising temperatures. Gravid females from both populations were subjected to current or projected end-of-century (future) thermal environments, to evaluate differences in the body size, growth rate, thermal preference, and locomotor performance of their offspring. Our results show that independently of the mode of reproduction, high temperatures accelerated gestation periods. Thermal environments did not affect hatchling thermal preference, but viviparous hatchlings consistently preferred lower temperatures. Unlike viviparous lizards, oviparous hatchlings incubated under future temperatures were smaller and had a lower growth rate compared to current-incubated hatchlings. Finally, thermal environments did not affect hatchling endurance and speed when controlling for body size. Our results show that global warming is likely to have a negative impact on S. equalis, but suggest that some of its effects may be ameliorated by maternal responses during pregnancy, particularly in viviparous populations.

Disease-Associated Circular RNAs: From Biology to Computational Identification

Circular RNAs (circRNAs) are endogenous RNAs with a covalently closed continuous loop, generated through various backsplicing events of pre-mRNA. An accumulating number of studies have shown that circRNAs are potential biomarkers for major human diseases such as cancer and Alzheimer's disease. Thus, identification and prediction of human disease-associated circRNAs are of significant importance. To this end, a computational analysis-assisted strategy is indispensable to detect, verify, and quantify circRNAs for downstream applications. In this review, we briefly introduce the biology of circRNAs, including the biogenesis, characteristics, and biological functions. In addition, we outline about 30 recent bioinformatic analysis tools that are publicly available for circRNA study. Principles for applying these computational strategies and considerations will be briefly discussed. Lastly, we give a complete survey on more than 20 key computational databases that are frequently used. To our knowledge, this is the most complete and updated summary on publicly available circRNA resources. In conclusion, this review summarizes key aspects of circRNA biology and outlines key computational strategies that will facilitate the genome-wide identification and prediction of circRNAs.

The evolution and physiology of male pregnancy in syngnathid fishes

The seahorses, pipefishes and seadragons (Syngnathidae) are among the few vertebrates in which pregnant males incubate developing embryos. Syngnathids are popular in studies of sexual selection, sex-role reversal, and reproductive trade-offs, and are now emerging as valuable comparative models for the study of the biology and evolution of reproductive complexity. These fish offer the opportunity to examine the physiology, behavioural implications, and evolutionary origins of embryo incubation, independent of the female reproductive tract and female hormonal milieu. Such studies allow us to examine flexibility in regulatory systems, by determining whether the pathways underpinning female pregnancy are also co-opted in incubating males, or whether novel pathways have evolved in response to the common challenges imposed by incubating developing embryos and releasing live young. The Syngnathidae are also ideal for studies of the evolution of reproductive complexity, because they exhibit multiple parallel origins of complex reproductive phenotypes. Here we assay the taxonomic distribution of syngnathid parity mode, examine the selective pressures that may have led to the emergence of male pregnancy, describe the biology of syngnathid reproduction, and highlight pressing areas for future research. Experimental tests of a range of hypotheses, including many generated with genomic tools, are required to inform overarching theories about the fitness implications of pregnancy and the evolution of male pregnancy. Such information will be widely applicable to our understanding of fundamental reproductive and evolutionary processes in animals.

Mother-offspring conflict for water and its mitigation in the oviparous form of the reproductively bimodal lizard, Zootoca vivipara

Parent-offspring conflicts are widespread given that resources are often limited. Recent evidence has shown that availability of water can trigger such conflict during pregnancy in viviparous squamate species (lizards and snakes) and thus questions the role of water in the evolution of reproductive modes. Here, we examined the impact of water restriction during gravidity in the oviparous form of the bimodal common lizard (Zootoca vivipara), using a protocol previously used on the viviparous form. Females were captured in early gravidity from six populations along a 600 m altitudinal gradient to investigate whether environmental conditions (altitude, water access and temperature) exacerbate responses to water restriction. Females were significantly dehydrated after water restriction, irrespective of their reproductive status (gravid vs. non-reproductive), relative reproductive effort (relative clutch mass), and treatment timing (embryonic development stage). Female dehydration, together with reproductive performance, varied with altitude, probably due to long term acclimation or local adaptation. This moderate water-based intergenerational conflict in gravid females contrasts sharply with previous findings for the viviparous form, with implications to the evolutionary reversion from viviparity to oviparity. It is likely that oviparity constitutes a water-saving reproductive mode which might help mitigate intensive temperature-driven population extinctions at low altitudes.