The draft genomes of soft-shell turtle and green sea turtle yield insights into the development and evolution of the turtle-specific body plan

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

slc26a12-A novel member of the slc26 family, is located in tandem with slc26a2 in coelacanths, amphibians, reptiles, and birds

Solute carrier family 26 (Slc26) is a family of anion exchangers with 11 members in mammals (named Slc26a1-a11). Here, we identified a novel member of the slc26 family, slc26a12, located in tandem with slc26a2 in the genomes of several vertebrate lineages. BLAST and synteny analyses of various jawed vertebrate genome databases revealed that slc26a12 is present in coelacanths, amphibians, reptiles, and birds but not in cartilaginous fishes, lungfish, mammals, or ray-finned fishes. In some avian and reptilian lineages such as owls, penguins, egrets, and ducks, and most turtles examined, slc26a12 was lost or pseudogenized. Phylogenetic analysis showed that Slc26a12 formed an independent branch with the other Slc26 members and Slc26a12, Slc26a1 and Slc26a2 formed a single branch, suggesting that these three members formed a subfamily in Slc26. In jawless fish, hagfish have two genes homologous to slc26a2 and slc26a12, whereas lamprey has a single gene homologous to slc26a2. African clawed frogs express slc26a12 in larval gills, skin, and fins. These results show that slc26a12 was present at least before the separation of lobe-finned fish and tetrapods; the name slc26a12 is appropriate because the gene duplication occurred in the distant past.

Multiple horizontal transfer events of a DNA transposon into turtles, fishes, and a frog

Horizontal transfer of transposable elements (HTT) has been reported across many species and the impact of such events on genome structure and function has been well described. However, few studies have focused on reptilian genomes, especially HTT events in Testudines (turtles). Here, as a consequence of investigating the repetitive content of Malaclemys terrapin terrapin (Diamondback turtle) we found a high similarity DNA transposon, annotated in RepBase as hAT-6_XT, shared between other turtle species, ray-finned fishes, and a frog. hAT-6_XT was notably absent in reptilian taxa closely related to turtles, such as crocodiles and birds. Successful invasion of DNA transposons into new genomes requires the conservation of specific residues in the encoded transposase, and through structural analysis, these residues were identified indicating some retention of functional transposition activity. We document six recent independent HTT events of a DNA transposon in turtles, which are known to have a low genomic evolutionary rate and ancient repeats.

When colors mislead: Genomics and bioacoustics prompt re-classification of Asian flycatcher radiation (Aves: Niltavinae)

Traditional classification of many animals, including birds, has been highly dependent on external morphological characters like plumage coloration. However, both bioacoustics and genetic or genomic data have revolutionized our understanding of the relationships of certain lineages and led to sweeping taxonomic re-organizations. In this study, we present a case of erroneous delimitation of genus boundaries in the species-rich flycatcher subfamily Niltavinae. Genera within this subfamily have historically been delineated based on blue versus brown male body plumage until recent studies based on a few mitochondrial and nuclear loci unearthed several cases of generic misclassification. Here we use extensive bioacoustic data from 43 species and genomic data from 28 species for a fundamental reclassification of species in the Niltavinae. Our study reveals that song is an important trait to classify these birds even at the genus level, whereas plumage traits exhibit ample convergence and have led to numerous historic misattributions. Our taxonomic re-organization leads to new biogeographic limits of major genera, such that the genus Cyornis now only extends as far east as the islands of Sulawesi, Sula, and Banggai, whereas Eumyias is redefined to extend far beyond Wallace's Line to the islands of Seram and Timor. Our conclusions advise against an over-reliance on morphological traits and underscore the importance of integrative datasets.

A chromosome-level genome assembly of the pig-nosed turtle (Carettochelys insculpta)

The pig-nosed turtle (Carettochelys insculpta) represents the only extant species within the Carettochelyidae family, is a unique Trionychia member fully adapted to aquatic life and currently facing endangerment. To enhance our understanding of this species and contribute to its conservation efforts, we employed high-fidelity (HiFi) and Hi-C sequencing technology to generate its genome assembly at the chromosome level. The assembly result spans 2.18 Gb, with a contig N50 of 126 Mb, encompassing 34 chromosomes that account for 99.6% of the genome. The assembly has a BUSCO score above 95% with different databases and strong collinearity with Yangtze giant softshell turtles (Rafetus swinhoei), indicating its completeness and continuity. A total of 19,175 genes and 46.86% repetitive sequences were annotated. The availability of this chromosome-scale genome represents a valuable resource for the pig-nosed turtle, providing insights into its aquatic adaptation and serving as a foundation for future turtle research.

Developmental hourglass: Verification by numerical evolution and elucidation by dynamical-systems theory

Determining the general laws between evolution and development is a fundamental biological challenge. Developmental hourglasses have attracted increased attention as candidates for such laws, but the necessity of their emergence remains elusive. We conducted evolutionary simulations of developmental processes to confirm the emergence of the developmental hourglass and unveiled its establishment. We considered organisms consisting of cells containing identical gene networks that control morphogenesis and evolved them under selection pressure to induce more cell types. By computing the similarity between the spatial patterns of gene expression of two species that evolved from a common ancestor, a developmental hourglass was observed, that is, there was a correlation peak in the intermediate stage of development. The fraction of pleiotropic genes increased, whereas the variance in individuals decreased, consistent with previous experimental reports. Reduction of the unavoidable variance by initial or developmental noise, essential for survival, was achieved up to the hourglass bottleneck stage, followed by diversification in developmental processes, whose timing is controlled by the slow expression dynamics conserved among organisms sharing the hourglass. This study suggests why developmental hourglasses are observed within a certain phylogenetic range of species.

Identification of Sex-Specific Markers and Candidate Genes Using WGS Sequencing Reveals a ZW-Type Sex-Determination System in the Chinese Soft-Shell Turtle (Pelodiscus sinensis)

Male and female Chinese soft-shelled turtles (Pelodiscus sinensis) have sex-dimorphic growth patterns, and males have higher commercial value because of their larger size and thicker calipash. Thus, developing sex-specific markers is beneficial to studies on all-male breeding in P. sinensis. Here, we developed an accurate and efficient workflow for the screening of sex-specific sequences with ZW or XY sex determination systems. Based on this workflow, female and male P. sinensis reference genomes of 2.23 Gb and 2.26 Gb were obtained using de novo assembly. After aligning and filtering, 4.01 Mb female-specific sequences were finally identified. Subsequently, the seven developed sex-specific primer pairs were 100% accurate in preliminary, population, and embryonic validation. The presence and absence of bands for the primers of P44, P45, P66, P67, P68, and P69, as well as two and one bands for the PB1 primer, indicate that the embryos are genetically female and male, respectively. NR and functional annotations identified several sex-determining candidate genes and related pathways, including Ran, Eif4et, and Crkl genes, and the insulin signaling pathway and the cAMP signaling pathway, respectively. Collectively, our results reveal that a ZW-type sex-determination system is present in P. sinensis and provide novel insights for the screening of sex-specific markers, sex-control breeding, and the studies of the sex determination mechanism of P. sinensis.

In situ hybridization analysis of odorant receptor expression in the olfactory organ of the pig-nosed turtle Carettochelys insculpta

The turtle olfactory organ consists of upper (UCE) and lower (LCE) chamber epithelium, which send axons to the ventral and dorsal portions of the olfactory bulbs, respectively. Generally, the UCE is associated with glands and contains ciliated olfactory receptor neurons (ORNs), while the LCE is devoid of glands and contains microvillous ORNs. However, the olfactory organ of the pig-nosed turtle Carettochelys insculpta appears to be a single olfactory system morphologically: there are no associated glands; ciliated ORNs are distributed throughout the olfactory organ; and the olfactory bulb is not divided into ventral and dorsal portions. In this study, we analyzed the expression of odorant receptors (ORs), the major olfactory receptors in turtles, in the pig-nosed turtle olfactory organ, via in situ hybridization. Of 690 ORs, 375 were classified as class I and 315 as class II. Some class II ORs were expressed predominantly in the posterior dorsomedial walls of the nasal cavity, while other class II ORs and all class I ORs examined were expressed in the remaining region. These results suggest that the pig-nosed turtle olfactory organ can be divided into two regions according to the expression of ORs.

The first high-quality chromosome-level genome of Eretmochelys imbricata using HiFi and Hi-C data

Eretmochelys imbricata, a critically endangered sea turtle inhabiting tropical oceans and protected across the world, had an unknown genome sequence until now. In this study, we used HiFi reads and Hi-C technology to assemble a high-quality, chromosome-level genome of E. imbricata. The genome size was 2,138.26 Mb, with contig N50 length of 123.49 Mb and scaffold N50 of 137.21 Mb. Approximately 97.52% of the genome sequence was anchored to 28 chromosomes. A total of 20,206 protein-coding genes were predicted. We also analyzed the evolutionary relationships, gene family expansions, and positive selection of E. imbricata. Our results revealed that E. imbricata diverged from Chelonia mydas 38 million years ago and had enriched olfactory receptors and aging-related genes. Our genome will be useful for studying E. imbricata and its conservation.

Histological features and Gαolf expression patterns in the nasal cavity of sea turtles

Sea turtles use olfaction to detect volatile and water-soluble substances. The nasal cavity of green turtles (Chelonia mydas) comprises morphologically defined the anterodorsal, anteroventral, and posterodorsal diverticula, as well as a single posteroventral fossa. Here, we detailed the histological features of the nasal cavity of a mature female green turtle. The posterodorsal diverticulum contained spongy-like venous sinuses and a wave-shaped sensory epithelium that favored ventilation. Secretory structures that were significant in sensory and non-sensory epithelia were probably involved in protection against seawater. These findings suggested that green turtles efficiently intake airborne substances and dissolve water-soluble substances in mucous, while suppressing the effects of salts. In addition, positive staining of Gαs/olf that couples with olfactory, but not vomeronasal, receptors was predominant in all three types of sensory epithelium in the nasal cavity. Both of airborne and water-soluble odorants seemed to be detected in cells expressing Gαolf and olfactory receptors.

The State of Squamate Genomics: Past, Present, and Future of Genome Research in the Most Speciose Terrestrial Vertebrate Order

Squamates include more than 11,000 extant species of lizards, snakes, and amphisbaenians, and display a dazzling diversity of phenotypes across their over 200-million-year evolutionary history on Earth. Here, we introduce and define squamates (Order Squamata) and review the history and promise of genomic investigations into the patterns and processes governing squamate evolution, given recent technological advances in DNA sequencing, genome assembly, and evolutionary analysis. We survey the most recently available whole genome assemblies for squamates, including the taxonomic distribution of available squamate genomes, and assess their quality metrics and usefulness for research. We then focus on disagreements in squamate phylogenetic inference, how methods of high-throughput phylogenomics affect these inferences, and demonstrate the promise of whole genomes to settle or sustain persistent phylogenetic arguments for squamates. We review the role transposable elements play in vertebrate evolution, methods of transposable element annotation and analysis, and further demonstrate that through the understanding of the diversity, abundance, and activity of transposable elements in squamate genomes, squamates can be an ideal model for the evolution of genome size and structure in vertebrates. We discuss how squamate genomes can contribute to other areas of biological research such as venom systems, studies of phenotypic evolution, and sex determination. Because they represent more than 30% of the living species of amniote, squamates deserve a genome consortium on par with recent efforts for other amniotes (i.e., mammals and birds) that aim to sequence most of the extant families in a clade.

Genomic insights into the genetic basis of eagle-beak jaw, large head, and long tail in the big-headed turtle

The big-headed turtle (Platysternon megacephalum) is an endemic chelonian species in Asia. Unlike most other turtles in the world, P. megacephalum is characterized with eagle-beak jaw, large head, and long tail. Although these unique characteristics are well recognized, the underlying genetic basis remains largely elusive. Here, we performed comparative genomic analysis between P. megacephalum and other representative species, aiming to reveal the genetic basis of the unique morphological features. Our results revealed that the eagle-beak jaw is most likely enabled by combined effects of expansion of SFRP5, extraction of FGF11, and mutation of both ZFYVE16 and PAX6. Large head is supported by mutations of SETD2 and FGRF2 and copy number variations of six head circumference modulation-related genes (TGFBR2, Twist2, Rdh10, Gas1, Chst11, and SNAP25). The long tail is probably involved in a genetic network comprising Gdf11, Lin 28, and HoxC12, two of which showed a consistent expression pattern with a model organism (mice). These findings suggest that expansion, extraction, and mutation of those genes may have profound effects on unique phenotypes of P. megacephalum.

Chromosome-Level Analysis of the Pelochelys cantorii Genome Provides Insights to Its Immunity, Growth and Longevity

The Asian giant soft-shelled turtle, Pelochelys cantorii (Trionychidae), is one of the largest aquatic turtles in China and was designated as a First-Grade Protected Animal in China in 1989. Previous investigation based on a combination of Illumina short-read, PacBio long-read and Hi-C scaffolding technologies acquired a high-quality chromosome-level genome of Pc. cantorii. In this study, comparative genomic analysis between Pc. cantorii and 16 other vertebrate genomes indicated that turtles separated from the ancestor of archosaurians approximately 256.6 (95% highest posterior density interval, 263.6-251.9) million years ago (Mya) (Upper Permian to Triassic) and that Pc. cantorii separated from the ancestor of Pd. sinensis and R. swinhoei approximately 59.3 (95% highest posterior density interval, 64.3-54.3) Mya. Moreover, several candidate genes, such as VWA5A, ABCG2, A2M and IGSF1, associated with tumor suppression, growth and age were expanded, implicating their potential roles in the exceptional longevity of turtles. This new chromosome-level assembly has important scientific value in the study of conservation of Pc. cantorii and also enriches the evolutionary investigation of turtle species.

The genome sequence of the Loggerhead sea turtle, Caretta caretta Linnaeus 1758

We present a genome assembly of Caretta caretta (the Loggerhead sea turtle; Chordata, Testudines, Cheloniidae), generated from genomic data from two unrelated females. The genome sequence is 2.13 gigabases in size. The assembly has a busco completion score of 96.1% and N50 of 130.95 Mb. The majority of the assembly is scaffolded into 28 chromosomal representations with a remaining 2% of the assembly being excluded from these.

In situ hybridization analysis of olfactory receptor expression in the sea turtle olfactory organ

The olfactory organ of turtles consists of an upper chamber epithelium (UCE) with associated glands, and a lower chamber epithelium (LCE) devoid of glands. The UCE and LCE are referred to as the air-nose and the water-nose, respectively, because the UCE is thought to detect airborne odorants, while the LCE detects waterborne odorants. However, it is not clear how the two are used in the olfactory organ. Odorant receptors (ORs) are the major olfactory receptors in turtles; they are classified as class I and II ORs, distinguished by their primary structure. Class I ORs are suggested to be receptive to water-soluble ligands and class II ORs to volatile ligands. This study analyzed the expression of class I and II ORs in hatchlings of the green sea turtle, Chelonia mydas, through in situ hybridization, to determine the localization of OR-expressing cells in the olfactory organ. Class I OR-expressing cells were distributed mainly in the LCE, implying that the LCE is receptive to waterborne odorants. Class II OR-expressing cells were distributed in both the UCE and LCE, implying that the entire olfactory organ is receptive to airborne odorants. The widespread expression of class II ORs may increase opportunities for sea turtles to sense airborne odorants.

Could theropod dinosaurs have evolved to a human level of intelligence?

Noting that some theropod dinosaurs had large brains, large grasping hands, and likely binocular vision, paleontologist Dale Russell suggested that a branch of these dinosaurs might have evolved to a human intelligence level, had dinosaurs not become extinct. I offer reasons why the likely pallial organization in dinosaurs would have made this improbable, based on four assumptions. First, it is assumed that achieving human intelligence requires evolving an equivalent of the about 200 functionally specialized cortical areas characteristic of humans. Second, it is assumed that dinosaurs had an avian nuclear type of pallial organization, in contrast to the mammalian cortical organization. Third, it is assumed that the interactions between the different neuron types making up an information processing unit within pallium are critical to its role in analyzing information. Finally, it is assumed that increasing axonal length between the neuron sets carrying out this operation impairs its efficacy. Based on these assumptions, I present two main reasons why dinosaur pallium might have been unable to add the equivalent of 200 efficiently functioning cortical areas. First, a nuclear pattern of pallial organization would require increasing distances between the neuron groups corresponding to the separate layers of any given mammalian cortical area, as more sets of nuclei equivalent to a cortical area are interposed between the existing sets, increasing axon length and thereby impairing processing efficiency. Second, because of its nuclear organization, dinosaur pallium could not reduce axon length by folding to bring adjacent areas closer together, as occurs in cerebral cortex.

Stability in gene expression and body-plan development leads to evolutionary conservation

BACKGROUND

Phenotypic evolution is mainly explained by selection for phenotypic variation arising from factors including mutation and environmental noise. Recent theoretical and experimental studies have suggested that phenotypes with greater developmental stability tend to have a constant phenotype and gene expression level within a particular genetic and environmental condition, and this positively correlates with stronger evolutionary conservation, even after the accumulation of genetic changes. This could reflect a novel mechanism that contributes to evolutionary conservation; however, it remains unclear whether developmental stability is the cause, or whether at least it contributes to their evolutionary conservation. Here, using Japanese medaka lines, we tested experimentally whether developmental stages and gene expression levels with greater stability led to their evolutionary conservation.

RESULTS

We first measured the stability of each gene expression level and developmental stage (defined here as the whole embryonic transcriptome) in the inbred F0 medaka population. We then measured their evolutionary conservation in the F3 generation by crossing the F0 line with the distantly related Japanese medaka line (Teradomori), followed by two rounds of intra-generational crossings. The results indicated that the genes and developmental stages that had smaller variations in the F0 generation showed lower diversity in the hybrid F3 generation, which implies a causal relationship between stability and evolutionary conservation.

CONCLUSIONS

These findings suggest that the stability in phenotypes, including the developmental stages and gene expression levels, leads to their evolutionary conservation; this most likely occurs due to their low potential to generate phenotypic variation. In addition, since the highly stable developmental stages match with the body-plan-establishment stage, it also implies that the developmental stability potentially contributed to the strict conservation of animal body plan.

Divergent sensory and immune gene evolution in sea turtles with contrasting demographic and life histories

Sea turtles represent an ancient lineage of marine vertebrates that evolved from terrestrial ancestors over 100 Mya. The genomic basis of the unique physiological and ecological traits enabling these species to thrive in diverse marine habitats remains largely unknown. Additionally, many populations have drastically declined due to anthropogenic activities over the past two centuries, and their recovery is a high global conservation priority. We generated and analyzed high-quality reference genomes for the leatherback (Dermochelys coriacea) and green (Chelonia mydas) turtles, representing the two extant sea turtle families. These genomes are highly syntenic and homologous, but localized regions of noncollinearity were associated with higher copy numbers of immune, zinc-finger, and olfactory receptor (OR) genes in green turtles, with ORs related to waterborne odorants greatly expanded in green turtles. Our findings suggest that divergent evolution of these key gene families may underlie immunological and sensory adaptations assisting navigation, occupancy of neritic versus pelagic environments, and diet specialization. Reduced collinearity was especially prevalent in microchromosomes, with greater gene content, heterozygosity, and genetic distances between species, supporting their critical role in vertebrate evolutionary adaptation. Finally, diversity and demographic histories starkly contrasted between species, indicating that leatherback turtles have had a low yet stable effective population size, exhibit extremely low diversity compared with other reptiles, and harbor a higher genetic load compared with green turtles, reinforcing concern over their persistence under future climate scenarios. These genomes provide invaluable resources for advancing our understanding of evolution and conservation best practices in an imperiled vertebrate lineage.

The Story of the Finest Armor: Developmental Aspects of Reptile Skin

The reptile skin is a barrier against water loss and pathogens and an armor for mechanical damages. The integument of reptiles consists of two main layers: the epidermis and the dermis. The epidermis, the hard cover of the body which has an armor-like role, varies among extant reptiles in terms of structural aspects such as thickness, hardness or the kinds of appendages it constitutes. The reptile epithelial cells of the epidermis (keratinocytes) are composed of two main proteins: intermediate filament keratins (IFKs) and corneous beta proteins (CBPs). The outer horny layer of the epidermis, stratum corneum, is constituted of keratinocytes by means of terminal differentiation or cornification which is a result of the protein interactions where CBPs associate with and coat the initial scaffold of IFKs. Reptiles were able to colonize the terrestrial environment due to the changes in these epidermal structures, which led to various cornified epidermal appendages such as scales and scutes, a beak, claws or setae. Developmental and structural aspects of the epidermal CBPs as well as their shared chromosomal locus (EDC) indicate an ancestral origin that gave rise to the finest armor of reptilians.

Creating morphological diversity in reptilian temporal skull region: A review of potential developmental mechanisms

Reptilian skull morphology is highly diverse and broadly categorized into three categories based on the number and position of the temporal fenestrations: anapsid, synapsid, and diapsid. According to recent phylogenetic analysis, temporal fenestrations evolved twice independently in amniotes, once in Synapsida and once in Diapsida. Although functional aspects underlying the evolution of tetrapod temporal fenestrations have been well investigated, few studies have investigated the developmental mechanisms responsible for differences in the pattern of temporal skull region. To determine what these mechanisms might be, we first examined how the five temporal bones develop by comparing embryonic cranial osteogenesis between representative extant reptilian species. The pattern of temporal skull region may depend on differences in temporal bone growth rate and growth direction during ontogeny. Next, we compared the histogenesis patterns and the expression of two key osteogenic genes, Runx2 and Msx2, in the temporal region of the representative reptilian embryos. Our comparative analyses suggest that the embryonic histological condition of the domain where temporal fenestrations would form predicts temporal skull morphology in adults and regulatory modifications of Runx2 and Msx2 expression in osteogenic mesenchymal precursor cells are likely involved in generating morphological diversity in the temporal skull region of reptiles.

Reptilian Innate Immunology and Ecoimmunology: What Do We Know and Where Are We Going?

Reptiles, the only ectothermic amniotes, employ a wide variety of physiological adaptations to adjust to their environments but remain vastly understudied in the field of immunology and ecoimmunology in comparison to other vertebrate taxa. To address this knowledge gap, we assessed the current state of research on reptilian innate immunology by conducting an extensive literature search of peer-reviewed articles published across the four orders of Reptilia (Crocodilia, Testudines, Squamata, and Rhynchocephalia). Using our compiled dataset, we investigated common techniques, characterization of immune components, differences in findings and type of research among the four orders, and immune responses to ecological and life-history variables. We found that there are differences in the types of questions asked and approaches used for each of these reptilian orders. The different conceptual frameworks applied to each group has led to a lack of unified understanding of reptilian immunological strategies, which, in turn, have resulted in large conceptual gaps in the field of ecoimmunology as a whole. To apply ecoimmunological concepts and techniques most effectively to reptiles, we must combine traditional immunological studies with ecoimmunological studies to continue to identify, characterize, and describe the reptilian immune components and responses. This review highlights the advances and gaps that remain to help identify targeted and cohesive approaches for future research in reptilian ecoimmunological studies.

Single-cell transcriptomes from turtle livers reveal sensitivity of hepatic immune cells to bacteria-infection

The liver is important in the synthesis, metabolism and storage of nutrients, detoxification and immune response of the body, and the liver immune response against exogenous pathogens from the intestinal tract plays a key role in the immune activities. However, the cellular composition of the liver immune atlas remains sparsely studied in reptiles. We used single-cell RNA sequencing to identify the cellular profile of the liver of the Chinese soft-shelled turtle (Pelodiscus sinensis). We obtained the transcriptional landscape based on 9938 cells from the fractionation of fresh hepatic tissues from two individuals, uninfected and infected with bacteria (Aeromonas hydrophila). We identified seven hepatic immune cell subsets, including plasma, erythroid, T/NK, B, endothelial, dendritic and Kupffer cells. Bacteria-infection altered the number of liver immune cells, as revealed by the fact that the infected turtle had more plasma, endothelial and Kupffer cells and fewer T/NK, dendritic and erythroid cells than did the uninfected turtle. Our study is the first to provide a comprehensive view of the hepatic immune landscape of P. sinensis at the single-cell resolution that outlines the characteristics of immune cells in the turtle liver and provides a liver transcriptome baseline for turtle immunology.

Comparative genomic and transcriptomic analyses of trans-kingdom pathogen Fusarium solani species complex reveal degrees of compartmentalization

BACKGROUND

The Fusarium solani species complex (FSSC) comprises fungal pathogens responsible for mortality in a diverse range of animals and plants, but their genome diversity and transcriptome responses in animal pathogenicity remain to be elucidated. We sequenced, assembled and annotated six chromosome-level FSSC clade 3 genomes of aquatic animal and plant host origins. We established a pathosystem and investigated the expression data of F. falciforme and F. keratoplasticum in Chinese softshell turtle (Pelodiscus sinensis) host.

RESULTS

Comparative analyses between the FSSC genomes revealed a spectrum of conservation patterns in chromosomes categorised into three compartments: core, fast-core (FC), and lineage-specific (LS). LS chromosomes contribute to variations in genomes size, with up to 42.2% of variations between F. vanettenii strains. Each chromosome compartment varied in structural architectures, with FC and LS chromosomes contain higher proportions of repetitive elements with genes enriched in functions related to pathogenicity and niche expansion. We identified differences in both selection in the coding sequences and DNA methylation levels between genome features and chromosome compartments which suggest a multi-speed evolution that can be traced back to the last common ancestor of Fusarium. We further demonstrated that F. falciforme and F. keratoplasticum are opportunistic pathogens by inoculating P. sinensis eggs and identified differentially expressed genes also associated with plant pathogenicity. These included the most upregulated genes encoding the CFEM (Common in Fungal Extracellular Membrane) domain.

CONCLUSIONS

The high-quality genome assemblies provided new insights into the evolution of FSSC chromosomes, which also serve as a resource for studies of fungal genome evolution and pathogenesis. This study also establishes an animal model for fungal pathogens of trans-kingdom hosts.

Genome of the endangered Guatemalan Beaded Lizard, Heloderma charlesbogerti, reveals evolutionary relationships of squamates and declines in effective population sizes

Many lizard species face extinction due to worldwide climate change. The Guatemalan Beaded Lizard, Heloderma charlesbogerti, is a member of the Family Helodermatidae that may be particularly imperiled; fewer than 600 mature individuals are believed to persist in the wild. In addition, H. charlesbogerti lizards are phenotypically remarkable. They are large in size, charismatically patterned, and possess a venomous bite. Here, we report the draft genome of the Guatemalan Beaded Lizard using DNA from a wild-caught individual. The assembled genome totals 2.31 Gb in length, similar in size to the genomes of related species. Single-copy orthologs were used to produce a novel molecular phylogeny, revealing that the Guatemalan Beaded Lizard falls into a clade with the Asian Glass Lizard (Anguidae) and in close association with the Komodo Dragon (Varanidae) and the Chinese Crocodile Lizard (Shinisauridae). In addition, we identified 31,411 protein-coding genes within the genome. Of the genes identified, we found 504 that evolved with a differential constraint on the branch leading to the Guatemalan Beaded Lizard. Lastly, we identified a decline in the effective population size of the Guatemalan Beaded Lizard approximately 400,000 years ago, followed by a stabilization before starting to dwindle again 60,000 years ago. The results presented here provide important information regarding a highly endangered, venomous reptile that can be used in future conservation, functional genetic, and phylogenetic analyses.

Gene Regulation during Carapacial Ridge Development of Mauremys reevesii: The Development of Carapacial Ridge, Ribs and Scutes

The unique topological structure of a turtle shell, including the special ribs-scapula relationship, is an evolutionarily novelty of amniotes. The carapacial ridge is a key embryonic tissue for inducing turtle carapace morphologenesis. However, the gene expression profiles and molecular regulatory mechanisms that occur during carapacial ridge development, including the regulation mechanism of rib axis arrest, the development mechanism of the carapacial ridge, and the differentiation between soft-shell turtles and hard-shell turtles, are not fully understood. In this study, we obtained genome-wide gene expression profiles during the carapacial ridge development of Mauremys reevesii using RNA-sequencing by using carapacial ridge tissues from stage 14, 15 and 16 turtle embryos. In addition, a differentially expressed genes (DEGs) analysis and a gene set enrichment analysis (GSEA) of three comparison groups were performed. Furthermore, a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis was used to analyze the pathway enrichment of the differentially expressed genes of the three comparative groups. The result displayed that the Wnt signaling pathway was substantially enriched in the CrTK14 vs. the CrTK15 comparison group, while the Hedgehog signaling pathway was significantly enriched in the CrTK15 vs. the CrTK16 group. Moreover, the regulatory network of the Wnt signaling pathway showed that Wnt signaling pathways might interact with Fgfs, Bmps, and Shh to form a regulatory network to regulate the carapacial ridge development. Next, WGCNA was used to cluster and analyze the expression genes during the carapacial ridge development of M. reevesii and P. sinensis. Further, a KEGG functional enrichment analysis of the carapacial ridge correlation gene modules was performed. Interesting, these results indicated that the Wnt signaling pathway and the MAPK signaling pathway were significantly enriched in the gene modules that were highly correlated with the stage 14 and stage 15 carapacial ridge samples of the two species. The Hedgehog signaling pathway was significantly enriched in the modules that were strongly correlated with the stage 16 carapacial ridge samples of M. reevesii, however, the PI3K-Akt signaling and the TGF-β signaling pathways were significantly enriched in the modules that were strongly correlated with the stage 16 carapacial ridge samples of P. sinensis. Furthermore, we found that those modules that were strongly correlated with the stage 14 carapacial ridge samples of M. reevesii and P. sinensis contained Wnts and Lef1. While the navajo white 3 module which was strongly correlated with the stage 16 carapacial ridge samples of M. reevesii contained Shh and Ptchs. The dark green module strongly correlated with the stage 16 carapacial ridge samples of P. sinensis which contained Col1a1, Col1a2, and Itga8. Consequently, this study systematically revealed the signaling pathways and genes that regulate the carapacial ridge development of M. reevesii and P. sinensis, which provides new insights for revealing the molecular mechanism that is underlying the turtle's body structure.

Functional complexity in the chorioallantoic membrane of an oviparous snake: Specializations for calcium uptake from the eggshell

The chorioallantoic membrane of oviparous reptiles forms a vascular interface with the eggshell. The eggshell contains calcium, primarily as calcium carbonate. Extraction and mobilization of this calcium by the chorioallantoic membrane contributes importantly to embryonic nutrition. Development of the chorioallantoic membrane is primarily known from studies of squamates and birds. Although there are pronounced differences in eggshell structure, squamate and bird embryos each mobilize calcium from eggshells. Specialized cells in the chicken chorionic epithelium transport calcium from the eggshell aided by a second population of cells that secrete protons generated by the enzyme carbonic anhydrase. Calcium transporting cells also are present in the chorioallantoic membrane of corn snakes, although these cells function differently than those of chickens. We used histology and immunohistology to characterize the morphology and functional attributes of the chorioallantoic membrane of corn snakes. We identified two populations of cells in the outer layer of the chorionic epithelium. Calbindin-D_28K , a cellular marker for calcium transport expressed in squamate chorioallantoic membranes, is localized in large, flattened cells that predominate in the chorionic epithelium. Smaller cells, interspersed among the large cells, express carbonic anhydrase 2, an enzyme not previously localized in the chorionic epithelium of an oviparous squamate. These findings indicate that differentiation of chorionic epithelial cells contributes to extraction and transport of calcium from the eggshell. The presence of specializations of chorioallantoic membranes for calcium uptake from eggshells in chickens and corn snakes suggests that eggshell calcium was a source of embryonic nutrition early in the evolution of Sauropsida.

Gene synteny, evolution and antiviral activity of type I IFNs in a reptile species, the Chinese soft-shelled turtle Pelodiscus sinensis

Type I interferons (IFNs) are critical cytokines for the establishment of antiviral status in fish, amphibian, avian and mammal, but the knowledge on type I IFNs is rather limited in reptile. In this study, seven type I IFN genes, designed as IFN1 to IFN7, were identified from a reptile species, the Chinese soft-shelled turtle (Pelodiscus sinensis). These identified type I IFNs have relatively low protein identity, when compared with those in human and chicken; but they possess conserved cysteines, predicted multi-helix structure and N-terminal signal peptide. The Chinese soft-shelled turtle IFN1 to IFN5 have two exons and one intron, but IFN6 and IFN7 are the single-exon genes. Chinese soft-shelled turtle type I IFNs are located respectively on the two conserved reptile-bird loci, named as Locus a and Locus c, and are clustered into the four of the five reptile-bird groups (named as Groups I-V) based on phylogenetic evidence, due to the lack of IFNK in the turtle. Moreover, the Chinese soft-shelled turtle type I IFNs can be induced by soft-shelled turtle iridovirus (STIV) infection and show antiviral activity in soft-shelled turtle artery (STA) cells, except IFN6. In addition, due to the difference in genome organizations, such as the number of exons and introns of type I IFN genes from fish to mammal, the definition and evolution of 'intronless' type I IFN genes were discussed in lineages of vertebrates. Thus, the finding of type I IFNs on two different loci in P. sinensis sheds light on the evolution of type I IFN genes in vertebrates.

Computational analysis on two putative mitochondrial protein-coding genes from the Emydura subglobosa genome: A functional annotation approach

Rapid advancements in automated genomic technologies have uncovered many unique findings about the turtle genome and its associated features including olfactory gene expansions and duplications of toll-like receptors. However, despite the advent of large-scale sequencing, assembly, and annotation, about 40-50% of genes in eukaryotic genomes are left without functional annotation, severely limiting our knowledge of the biological information of genes. Additionally, these automated processes are prone to errors since draft genomes consist of several disconnected scaffolds whose order is unknown; erroneous draft assemblies may also be contaminated with foreign sequences and propagate to cause errors in annotation. Many of these automated annotations are thus incomplete and inaccurate, highlighting the need for functional annotation to link gene sequences to biological identity. In this study, we have functionally annotated two genes of the red-bellied short-neck turtle (Emydura subglobosa), a member of the relatively understudied pleurodire lineage of turtles. We improved upon initial ab initio gene predictions through homology-based evidence and generated refined consensus gene models. Through functional, localization, and structural analyses of the predicted proteins, we discovered conserved putative genes encoding mitochondrial proteins that play a role in C21-steroid hormone biosynthetic processes and fatty acid catabolism-both of which are distantly related by the tricarboxylic acid (TCA) cycle and share similar metabolic pathways. Overall, these findings further our knowledge about the genetic features underlying turtle physiology, morphology, and longevity, which have important implications for the treatment of human diseases and evolutionary studies.

Glutamatergic pathways in the brains of turtles: A comparative perspective among reptiles, birds, and mammals

Glutamate acts as the main excitatory neurotransmitter in the brain and plays a vital role in physiological and pathological neuronal functions. In mammals, glutamate can cause detrimental excitotoxic effects under anoxic conditions. In contrast, Trachemys scripta, a freshwater turtle, is one of the most anoxia-tolerant animals, being able to survive up to months without oxygen. Therefore, turtles have been investigated to assess the molecular mechanisms of neuroprotective strategies used by them in anoxic conditions, such as maintaining low levels of glutamate, increasing adenosine and GABA, upregulating heat shock proteins, and downregulating KATP channels. These mechanisms of anoxia tolerance of the turtle brain may be applied to finding therapeutics for human glutamatergic neurological disorders such as brain injury or cerebral stroke due to ischemia. Despite the importance of glutamate as a neurotransmitter and of the turtle as an ideal research model, the glutamatergic circuits in the turtle brain remain less described whereas they have been well studied in mammalian and avian brains. In reptiles, particularly in the turtle brain, glutamatergic neurons have been identified by examining the expression of vesicular glutamate transporters (VGLUTs). In certain areas of the brain, some ionotropic glutamate receptors (GluRs) have been immunohistochemically studied, implying that there are glutamatergic target areas. Based on the expression patterns of these glutamate-related molecules and fiber connection data of the turtle brain that is available in the literature, many candidate glutamatergic circuits could be clarified, such as the olfactory circuit, hippocampal–septal pathway, corticostriatal pathway, visual pathway, auditory pathway, and granule cell–Purkinje cell pathway. This review summarizes the probable glutamatergic pathways and the distribution of glutamatergic neurons in the pallium of the turtle brain and compares them with those of avian and mammalian brains. The integrated knowledge of glutamatergic pathways serves as the fundamental basis for further functional studies in the turtle brain, which would provide insights on physiological and pathological mechanisms of glutamate regulation as well as neural circuits in different species.

Chromosome-level genome assembly of Asian yellow pond turtle (Mauremys mutica) with temperature-dependent sex determination system

Knowledge of sex determination has important implications in physiology, ecology and genetics, but the evolutionary mechanisms of sex determination systems in turtles have not been fully elucidated, due to a lack of reference genomes. Here, we generate a high-quality genome assembly of Asian yellow pond turtle (Mauremys mutica) using continuous long-read (PacBio platform), Illumina, and high-throughput chromatin conformation capture (Hi-C) technologies. The M. mutica haplotype has a genome size of 2.23 Gb with a contig N50 of 8.53 Mb and scaffold N50 of 141.98 Mb. 99.98% sequences of the total assembly are anchored to 26 pseudochromosomes. Comparative genomics analysis indicated that the lizard-snake-tuatara clade diverged from the bird-crocodilian-turtle clade at approximately 267.0-312.3 Mya. Intriguingly, positive selected genes are mostly enriched in the calcium signaling pathway and neuroactive ligand-receptor interaction, which are involved in the process of temperature-dependent sex determination. These findings provide important evolutionary insights into temperature-dependent sex determination system.

Potential contribution of intrinsic developmental stability toward body plan conservation

BACKGROUND

Despite the morphological diversity of animals, their basic anatomical patterns-the body plans in each animal phylum-have remained highly conserved over hundreds of millions of evolutionary years. This is attributed to conservation of the body plan-establishing developmental period (the phylotypic period) in each lineage. However, the evolutionary mechanism behind this phylotypic period conservation remains under debate. A variety of hypotheses based on the concept of modern synthesis have been proposed, such as negative selection in the phylotypic period through its vulnerability to embryonic lethality. Here we tested a new hypothesis that the phylotypic period is developmentally stable; it has less potential to produce phenotypic variations than the other stages, and this has most likely led to the evolutionary conservation of body plans.

RESULTS

By analyzing the embryos of inbred Japanese medaka embryos raised under the same laboratory conditions and measuring the whole embryonic transcriptome as a phenotype, we found that the phylotypic period has greater developmental stability than other stages. Comparison of phenotypic differences between two wild medaka populations indicated that the phylotypic period and its genes in this period remained less variational, even after environmental and mutational modifications accumulated during intraspecies evolution. Genes with stable expression levels were enriched with those involved in cell-cell signalling and morphological specification such as Wnt and Hox, implying possible involvement in body plan development of these genes.

CONCLUSIONS

This study demonstrated the correspondence between the developmental stage with low potential to produce phenotypic variations and that with low diversity in micro- and macroevolution, namely the phylotypic period. Whereas modern synthesis explains evolution as a process of shaping of phenotypic variations caused by mutations, our results highlight the possibility that phenotypic variations are readily limited by the intrinsic nature of organisms, namely developmental stability, thus biasing evolutionary outcomes.

A Novel Ex Vivo Approach Based on Proteomics and Biomarkers to Evaluate the Effects of Chrysene, MEHP, and PBDE-47 on Loggerhead Sea Turtles (Caretta caretta)

The principal aim of the present study was to develop and apply novel ex vivo tests as an alternative to cell cultures able to evaluate the possible effects of emerging and legacy contaminants in Caretta caretta. To this end, we performed ex vivo experiments on non-invasively collected whole-blood and skin-biopsy slices treated with chrysene, MEHP, or PBDE-47. Blood samples were tested by oxidative stress (TAS), immune system (respiratory burst, lysozyme, and complement system), and genotoxicity (ENA assay) biomarkers, and genotoxic and immune system effects were observed. Skin slices were analyzed by applying a 2D-PAGE/MS proteomic approach, and specific contaminant signatures were delineated on the skin proteomic profile. These reflect biochemical effects induced by each treatment and allowed to identify glutathione S-transferase P, peptidyl-prolyl cis-trans isomerase A, mimecan, and protein S100-A6 as potential biomarkers of the health-threatening impact the texted toxicants have on C. caretta. Obtained results confirm the suitability of the ex vivo system and indicate the potential risk the loggerhead sea turtle is undergoing in the natural environment. In conclusion, this work proved the relevance that the applied ex vivo models may have in testing the toxicity of other compounds and mixtures and in biomarker discovery.

Detection and population genomics of sea turtle species via non-invasive environmental DNA analysis of nesting beach sand tracks and oceanic water

Elusive aquatic wildlife, such as endangered sea turtles, are difficult to monitor and conserve. As novel molecular and genetic technologies develop, it is possible to adapt and optimize them for wildlife conservation. One such technology is environmental (e)DNA - the detection of DNA shed from organisms into their surrounding environments. We developed species-specific green (Chelonia mydas) and loggerhead (Caretta caretta) sea turtle probe-based qPCR assays, which can detect and quantify sea turtle eDNA in controlled (captive tank water and sand samples) and free ranging (oceanic water samples and nesting beach sand) settings. eDNA detection complemented traditional in-water sea turtle monitoring by enabling detection even when turtles were not visually observed. Furthermore, we report that high throughput shotgun sequencing of eDNA sand samples enabled sea turtle population genetic studies and pathogen monitoring, demonstrating that non-invasive eDNA techniques are viable and efficient alternatives to biological sampling (e.g. biopsies and blood draws). Genetic information was obtained from sand many hours after nesting events, without having to observe or interact with the target individual. This greatly reduces the sampling stress experienced by nesting mothers and emerging hatchlings, and avoids sacrificing viable eggs for genetic analysis. The detection of pathogens from sand indicates significant potential for increased wildlife disease monitoring capacity and viral variant surveillance. Together, these results demonstrate the potential of eDNA approaches to ultimately help understand and conserve threatened species such as sea turtles.

Distinguishing Evolutionary Conservation from Derivedness

While the concept of “evolutionary conservation” has enabled biologists to explain many ancestral features and traits, it has also frequently been misused to evaluate the degree of changes from a common ancestor, or “derivedness”. We propose that the distinction of these two concepts allows us to properly understand phenotypic and organismal evolution. From a methodological aspect, “conservation” mainly considers genes or traits which species have in common, while “derivedness” additionally covers those that are not commonly shared, such as novel or lost traits and genes to evaluate changes from the time of divergence from a common ancestor. Due to these differences, while conservation-oriented methods are effective in identifying ancestral features, they may be prone to underestimating the overall changes accumulated during the evolution of certain lineages. Herein, we describe our recently developed method, “transcriptomic derivedness index”, for estimating the phenotypic derivedness of embryos with a molecular approach using the whole-embryonic transcriptome as a phenotype. Although echinoderms are often considered as highly derived species, our analyses with this method showed that their embryos, at least at the transcriptomic level, may not be much more derived than those of chordates. We anticipate that the future development of derivedness-oriented methods could provide quantitative indicators for finding highly/lowly evolvable traits.

A Genomic Perspective on the Evolutionary Diversification of Turtles

To examine phylogenetic heterogeneity in turtle evolution, we collected thousands of high-confidence single-copy orthologs from 19 genome assemblies representative of extant turtle diversity and estimated a phylogeny with multispecies coalescent and concatenated partitioned methods. We also collected next-generation sequences from 26 turtle species and assembled millions of biallelic markers to reconstruct phylogenies based on annotated regions from the western painted turtle (Chrysemys picta bellii) genome (coding regions, introns, untranslated regions, intergenic, and others). We then measured gene tree-species tree discordance, as well as gene and site heterogeneity at each node in the inferred trees, and tested for temporal patterns in phylogenomic conflict across turtle evolution. We found strong and consistent support for all bifurcations in the inferred turtle species phylogenies. However, a number of genes, sites, and genomic features supported alternate relationships between turtle taxa. Our results suggest that gene tree-species tree discordance in these datasets is likely driven by population-level processes such as incomplete lineage sorting. We found very little effect of substitutional saturation on species tree topologies, and no clear phylogenetic patterns in codon usage bias and compositional heterogeneity. There was no correlation between gene and site concordance, node age, and DNA substitution rate across most annotated genomic regions. Our study demonstrates that heterogeneity is to be expected even in well resolved clades such as turtles, and that future phylogenomic studies should aim to sample as much of the genome as possible in order to obtain accurate phylogenies for assessing conservation priorities in turtles.

How the evolution of air breathing shaped hippocampal function

To make maps from airborne odours requires dynamic respiratory patterns. I propose that this constraint explains the modulation of memory by nasal respiration in mammals, including murine rodents (e.g. laboratory mouse, laboratory rat) and humans. My prior theories of limbic system evolution offer a framework to understand why this occurs. The answer begins with the evolution of nasal respiration in Devonian lobe-finned fishes. This evolutionary innovation led to adaptive radiations in chemosensory systems, including the emergence of the vomeronasal system and a specialization of the main olfactory system for spatial orientation. As mammals continued to radiate into environments hostile to spatial olfaction (air, water), there was a loss of hippocampal structure and function in lineages that evolved sensory modalities adapted to these new environments. Hence the independent evolution of echolocation in bats and toothed whales was accompanied by a loss of hippocampal structure (whales) and an absence of hippocampal theta oscillations during navigation (bats). In conclusion, models of hippocampal function that are divorced from considerations of ecology and evolution fall short of explaining hippocampal diversity across mammals and even hippocampal function in humans. This article is part of the theme issue 'Systems neuroscience through the lens of evolutionary theory'.

Gut microbiome-immune system interaction in reptiles

Reptiles are ectothermic amniotes in a world dominated by endotherms. Reptiles originated more than 300 million years ago and they often dwell in polluted environments which may expose them to pathogenic micro-organisms, radiation and/or heavy metals. Reptiles also possess greater longevity and may live much longer than similar-sized land mammals, for example, turtles, tortoises, crocodiles and tuatara are long-lived reptiles living up to 100 years or more. Many recent studies have emphasized the pivotal role of the gut microbiome on its host; thus, we postulated that reptilian gut microbiome and/or its metabolites and the interplay with their robust immune system may contribute to their longevity and overall hardiness. Herein, we discuss the composition of the reptilian gut microbiome, immune system-gut microbiome cross-talk, antimicrobial peptides, reptilian resistance to infectious diseases and cancer, ageing, as well the current knowledge of the genome and epigenome of these remarkable species. Preliminary studies have demonstrated that microbial gut flora of reptiles such as crocodiles, tortoises, water monitor lizard and python exhibit remarkable anticancer and antibacterial properties, as well as comprise novel gut bacterial metabolites and antimicrobial peptides. The underlying mechanisms between the gut microbiome and the immune system may hold clues to developing new therapies overall for health, and possible extrapolation to exploit the ancient defence systems of reptiles for Homo sapiens benefit.

Novel pegiviruses infecting wild birds and rodents

Pegivirus (family Flaviviridae) is a genus of small enveloped RNA viruses that mainly causes blood infections in various mammals including human. Herein, we carried out an extensive survey of pegiviruses from a wide range of wild animals mainly sampled in the Qinghai-Tibet Plateau of China. Three novel pegiviruses, namely Passer montanus pegivirus, Leucosticte brandti pegivirus and Montifringilla taczanowskii pegivirus, were identified from different wild birds, and one new rodent pegivirus, namely Phaiomys leucurus pegivirus, was identified from Blyth's vole. Interestingly, the pegiviruses of non-mammalian origin discovered in this study substantially broaden the host range of Pegivirus to avian species. Co-evolutionary analysis showed virus-host co-divergence over long evolutionary timescales, and indicated that pegiviruses largely followed a virus-host co-divergence relationship. Overall, this work extends the biodiversity of the Pegivirus genus to those infecting wild birds and hence revises the host range and evolutionary history of genus Pegivirus.

•

Novel pegiviruses were identified from wild-life animals in Qinghai-Tibet Plateau.

•

The three divergent species of bird pegiviruses substantially broaden the host range of Pegivirus.

•

A long-term evolutionary relationship was established between pegiviruses and their vertebrate hosts.

Dynamical systems approach to evolution-development congruence: Revisiting Haeckel's recapitulation theory

It is acknowledged that embryonic development has a tendency to proceed from common toward specific. Ernst Haeckel raised the question of why that tendency prevailed through evolution, and the question remains unsolved. Here, we revisit Haeckel's recapitulation theory, that is, the parallelism between evolution and development through numerical evolution and dynamical systems theory. By using intracellular gene expression dynamics with cell-to-cell interaction over spatially aligned cells to represent the developmental process, gene regulation networks (GRN) that govern these dynamics evolve under the selection pressure to achieve a prescribed spatial gene expression pattern. For most numerical evolutionary experiments, the evolutionary pattern changes over generations, as well as the developmental pattern changes governed by the evolved GRN exhibit remarkable similarity. Changes in both patterns consisted of several epochs where stripes are formed in a short time, whereas for other temporal regimes, the pattern hardly changes. In evolution, these quasi-stationary generations are needed to achieve relevant mutations, whereas, in development, they are due to some gene expressions that vary slowly and control the pattern change. These successive epochal changes in development and evolution are represented as common bifurcations in dynamical systems theory, regulating working network structure from feedforward subnetwork to those containing feedback loops. The congruence is the correspondence between successive acquisitions of subnetworks through evolution and changes in working subnetworks in development. Consistency of the theory with the segmentation gene-expression dynamics is discussed. Novel outlook on recapitulation and heterochrony are provided, testable experimentally by the transcriptome and network analysis.

Derivedness Index for Estimating Degree of Phenotypic Evolution of Embryos: A Study of Comparative Transcriptomic Analyses of Chordates and Echinoderms

Species retaining ancestral features, such as species called living fossils, are often regarded as less derived than their sister groups, but such discussions are usually based on qualitative enumeration of conserved traits. This approach creates a major barrier, especially when quantifying the degree of phenotypic evolution or degree of derivedness, since it focuses only on commonly shared traits, and newly acquired or lost traits are often overlooked. To provide a potential solution to this problem, especially for inter-species comparison of gene expression profiles, we propose a new method named "derivedness index" to quantify the degree of derivedness. In contrast to the conservation-based approach, which deals with expressions of commonly shared genes among species being compared, the derivedness index also considers those that were potentially lost or duplicated during evolution. By applying our method, we found that the gene expression profiles of penta-radial phases in echinoderm tended to be more highly derived than those of the bilateral phase. However, our results suggest that echinoderms may not have experienced much larger modifications to their developmental systems than chordates, at least at the transcriptomic level. In vertebrates, we found that the mid-embryonic and organogenesis stages were generally less derived than the earlier or later stages, indicating that the conserved phylotypic period is also less derived. We also found genes that potentially explain less derivedness, such as Hox genes. Finally, we highlight technical concerns that may influence the measured transcriptomic derivedness, such as read depth and library preparation protocols, for further improvement of our method through future studies. We anticipate that this index will serve as a quantitative guide in the search for constrained developmental phases or processes.

Demographic changes in Pleistocene sea turtles were driven by past sea level fluctuations affecting feeding habitat availability

Pleistocene environmental changes are generally assumed to have dramatically affected species’ demography via changes in habitat availability, but this is challenging to investigate due to our limited knowledge of how Pleistocene ecosystems changed through time. Here, we tracked changes in shallow marine habitat availability resulting from Pleistocene sea level fluctuations throughout the last glacial cycle (120–14 thousand years ago; kya) and assessed correlations with past changes in genetic diversity inferred from genome‐wide SNPs, obtained via ddRAD sequencing, in Caribbean hawksbill turtles, which feed in coral reefs commonly found in shallow tropical waters. We found sea level regression resulted in an average 75% reduction in shallow marine habitat availability during the last glacial cycle. Changes in shallow marine habitat availability correlated strongly with past changes in hawksbill turtle genetic diversity, which gradually declined to ~1/4th of present‐day levels during the Last Glacial Maximum (LGM; 26–19 kya). Shallow marine habitat availability and genetic diversity rapidly increased after the LGM, signifying a population expansion in response to warming environmental conditions. Our results suggest a positive correlation between Pleistocene environmental changes, habitat availability and species’ demography, and that demographic changes in hawksbill turtles were potentially driven by feeding habitat availability. However, we also identified challenges associated with disentangling the potential environmental drivers of past demographic changes, which highlights the need for integrative approaches. Our conclusions underline the role of habitat availability on species’ demography and biodiversity, and that the consequences of ongoing habitat loss should not be underestimated.

The population genomic structure of green turtles (Chelonia mydas) suggests a warm-water corridor for tropical marine fauna between the Atlantic and Indian oceans during the last interglacial

The occasional westward transport of warm water of the Agulhas Current, "Agulhas leakage", around southern Africa has been suggested to facilitate tropical marine connectivity between the Atlantic and Indian oceans, but the "Agulhas leakage" hypothesis does not explain the signatures of eastward gene flow observed in many tropical marine fauna. We investigated an alternative hypothesis: the establishment of a warm-water corridor during comparatively warm interglacial periods. The "warm-water corridor" hypothesis was investigated by studying the population genomic structure of Atlantic and Southwest Indian Ocean green turtles (N = 27) using 12,035 genome-wide single nucleotide polymorphisms (SNPs) obtained via ddRAD sequencing. Model-based and multivariate clustering suggested a hierarchical population structure with two main Atlantic and Southwest Indian Ocean clusters, and a Caribbean and East Atlantic sub-cluster nested within the Atlantic cluster. Coalescent-based model selection supported a model where Southwest Indian Ocean and Caribbean populations diverged from the East Atlantic population during the transition from the last interglacial period (130-115 thousand years ago; kya) to the last glacial period (115-90 kya). The onset of the last glaciation appeared to isolate Atlantic and Southwest Indian Ocean green turtles into three refugia, which subsequently came into secondary contact in the Caribbean and Southwest Indian Ocean when global temperatures increased after the Last Glacial Maximum. Our findings support the establishment of a warm-water corridor facilitating tropical marine connectivity between the Atlantic and Southwest Indian Ocean during warm interglacials.

Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species

Reconstructing past events of hybridization and population size changes are required to understand speciation mechanisms and current patterns of genetic diversity, and ultimately contribute to species' conservation. Sea turtles are ancient species currently facing anthropogenic threats including climate change, fisheries, and illegal hunting. Five of the seven extant sea turtle species are known to currently hybridize, especially along the Brazilian coast where some populations can have ~32%-42% of hybrids. Although frequently observed today, it is not clear what role hybridization plays in the evolutionary diversification of this group of reptiles. In this study, we generated whole genome resequencing data of the five globally distributed sea turtle species to estimate a calibrated phylogeny and the population size dynamics, and to understand the role of hybridization in shaping the genomes of these ancient species. Our results reveal discordant species divergence dates between mitochondrial and nuclear genomes, with a high frequency of conflicting trees throughout the nuclear genome suggesting that some sea turtle species frequently hybridized in the past. The reconstruction of the species' demography showed a general decline in effective population sizes with no signs of recovery, except for the leatherback sea turtle. Furthermore, we discuss the influence of reference bias in our estimates. We show long-lasting ancestral gene flow events within Chelonioidea that continued for millions of years after initial divergence. Speciation with gene flow is a common pattern in marine species, and it raises questions whether current hybridization events should be considered as a part of these species' evolutionary history or a conservation issue.

Chromosome-level genome assembly of the Chinese three-keeled pond turtle (Mauremys reevesii) provides insights into freshwater adaptation

Mauremys reevesii is an endangered freshwater turtle that symbolizes longevity in Chinese culture. Despite its importance, genetic studies of this species remain limited, with no genomic sequence reported to date. Here, we report a high-quality, chromosome-level genomic sequence of M. reevesii obtained using a combination of Nanopore and Hi-C sequencing technologies. The 2.37 Gb M. reevesii genome was assembled from a total of ~226.80 Gb of Nanopore sequencing data. The M. reevesii genome contig N50 is 34.73 Mb, the highest value in published turtle genomes. In total, 18,238 genes were functionally annotated. The contigs were clustered and ordered onto 27 pseudochromosomes covering ~96.55% of the genome assembled with Hi-C data. To explore genome evolution, synteny analysis was performed between M. reevesii (freshwater turtle) and Gopherus evgoodei (terrestrial turtle) genomes. In general, each chromosome of M. reevesii corresponded to one chromosome of Gopherus evgoodei, but some interchromosomal rearrangements occurred between the two species based on the assembled genomes. These interchromosomal rearrangements were further confirmed by mapping of the long-read nanopore data to the assembly. The reconstructed demographic history showed varied effective population size among freshwater, marine and terrestrial turtles. We also discovered expansion of genes related to the innate immune system in M. reevesii that may provide defence against freshwater pathogens. The high-quality genomic sequence provides a valuable genetic resource for further studies of genetics and genome evolution in turtles.

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.

Global Analysis of Transcriptome and Translatome Revealed That Coordinated WNT and FGF Regulate the Carapacial Ridge Development of Chinese Soft-Shell Turtle

The turtle carapace is composed of severely deformed fused dorsal vertebrae, ribs, and bone plates. In particular, the lateral growth in the superficial layer of turtle ribs in the dorsal trunk causes an encapsulation of the scapula and pelvis. The recent study suggested that the carapacial ridge (CR) is a new model of epithelial–mesenchymal transition which is essential for the arrangement of the ribs. Therefore, it is necessary to explore the regulatory mechanism of carapacial ridge development to analyze the formation of the turtle shell. However, the current understanding of the regulatory network underlying turtle carapacial ridge development is poor due to the lack of both systematic gene screening at different carapacial ridge development stages and gene function verification studies. In this study, we obtained genome-wide gene transcription and gene translation profiles using RNA sequencing and ribosome nascent-chain complex mRNA sequencing from carapacial ridge tissues of Chinese soft-shell turtle at different development stages. A correlation analysis of the transcriptome and translatome revealed that there were 129, 670, and 135 codifferentially expressed genes, including homodirection and opposite-direction differentially expressed genes, among three comparison groups, respectively. The pathway enrichment analysis of codifferentially expressed genes from the Kyoto Encyclopedia of Genes and Genomes showed dynamic changes in signaling pathways involved in carapacial ridge development. Especially, the results revealed that the Wnt signaling pathway and MAPK signaling pathway may play important roles in turtle carapacial ridge development. In addition, Wnt and Fgf were expressed during the carapacial ridge development. Furthermore, we discovered that Wnt5a regulated carapacial ridge development through the Wnt5a/JNK pathway. Therefore, our studies uncover that the morphogenesis of the turtle carapace might function through the co-operation between conserved WNT and FGF signaling pathways. Consequently, our findings revealed the dynamic signaling pathways acting on the carapacial ridge development of Chinese soft-shell turtle and provided new insights into uncover the molecular mechanism underlying turtle shell morphogenesis.