Gekko japonicus genome reveals evolution of adhesive toe pads and tail regeneration

Skin Appendage Proteins of Tetrapods: Building Blocks of Claws, Feathers, Hair and Other Cornified Epithelial Structures

Reptiles, birds, mammals and amphibians, together forming the clade tetrapods, have a large diversity of cornified skin appendages, such as scales, feathers, hair and claws. The skin appendages consist of dead epithelial cells that are tightly packed with specific structural proteins. Here, we review the molecular diversity and expression patterns of major types of skin appendage proteins, namely keratin intermediate filament proteins, keratin-associated proteins (KRTAPs) and proteins encoded by genes of the epidermal differentiation complex (EDC), including corneous beta-proteins, also known as beta-keratins. We summarize the current knowledge about the components of skin appendages with a focus on keratins and EDC proteins that have recently been identified in reptiles and birds. We discuss gaps of knowledge and suggest directions of future research.

Convergent Evolution Has Led to the Loss of Claw Proteins in Snakes and Worm Lizards

The evolution of cornified skin appendages, such as hair, feathers, and claws, is closely linked to the evolution of proteins that establish the unique mechanical stability of these epithelial structures. We hypothesized that the evolution of the limbless body anatomy of the Florida worm lizard (Rhineura floridana) and the concomitant loss of claws had led to the degeneration of genes with claw-associated functions. To test this hypothesis, we investigated the evolution of three gene families implicated in epithelial cell architecture, namely type I keratins, type II keratins, and genes of the epidermal differentiation complex in R. floridana in comparison with other squamates. We report that the orthologs of mammalian hair and nail keratins have undergone pseudogenization in R. floridana. Likewise, the epidermal differentiation complex genes tentatively named EDYM1 and EDCCs have been lost in R. floridana. The aforementioned genes are conserved in various lizards with claws, but not in snakes. Proteomic analysis of the cornified claws of the bearded dragon (Pogona vitticeps) confirmed that type I and type II hair keratin homologs, EDYM1 and EDCCs, are protein components of claws in squamates. We conclude that the convergent evolution of a limbless body was associated with the convergent loss of claw keratins and differentiation genes in squamates.

Molecular Phylogenetic Relationships Based on Mitogenomes of Spider: Insights Into Evolution and Adaptation to Extreme Environments

In this study, we performed a comparative analysis based on a total of 255 spider mitogenomes and four outgroups, of which the mitogenomes of 39 species were assembled de novo, to explore the phylogenetic relationships and the adaptive evolution of mitogenomes. Results showed that Argyroneta aquatica had the longest mitochondrial length and the most pronounced codon preference to be UUA, followed by CCU. Codon usage frequencies were similar between families and codon usage in the mitogenome of spiders was mainly influenced by natural selection pressures rather than G/C mutation bias. Our phylogenetic topology clearly explained the evolutionary relationships among the spiders, and divergence time estimates indicated that the spiders originated in the early Devonian, and that the two clades of Mesothelae and Opisthothelae separated in the late Carboniferous. Ancestral range and trait reconstruction results supported the ancestral origin of spiders to the Devonian Nearctic realm, with the trapdoor being the original trait. Selection analysis detected positive selection signals in the ATP8 gene in Desis jiaxiangi. The ND5 gene is a convergent evolutionary gene between D. jiaxiangi and A. aquatica. Positive selection signals in the ATP8 gene and convergent selection sites in the ND5 gene may facilitate metabolic adaptation to the aquatic environment in two aquatic spiders. In conclusion, our analysis contributes to a better understanding of the taxonomic status, species diversity, mitochondrial characteristics, and environmental adaptations of these spiders.

Bio-inspired designs: leveraging biological brilliance in mechanical engineering-an overview

Nature's evolutionary mastery has perfected design over the years, yielding organisms superbly adapted to their surroundings. This research delves into the promising domain of bio-inspired designs, poised to revolutionize mechanical engineering. Leveraging insights drawn from prior conversations, we categorize innovations influenced by life on land, in water, and through the air, emphasizing their pivotal contributions to mechanical properties. Our comprehensive review reveals a wealth of bio-inspired designs that have already made substantial inroads in mechanical engineering. From avian-inspired lightweight yet robust materials to hydrodynamically optimized forms borrowed from marine creatures, these innovations hold immense potential for enhancing mechanical systems. In conclusion, this study underscores the transformative potential of bio-inspired designs, offering improved mechanical characteristics and the promise of sustainability and efficiency across a broad spectrum of applications. This research envisions a future where bio-inspired designs shape the mechanical landscape, fostering a more harmonious coexistence between human technology and the natural world.

Development of α-Helical Antimicrobial Peptides with Imperfect Amphipathicity for Superior Activity and Selectivity

The advancement of antimicrobial peptides (AMPs) as therapeutic agents is hindered by their poor selectivity. Recent evidence indicates that controlled disruption of the amphipathicity of α-helical AMPs may increase the selectivity. This study investigated the role of imperfect amphipathicity in optimizing AMPs with varied sequences to enhance their activity and selectivity. Among these, the lead peptide RI-18, characterized by an imperfectly amphipathic α-helical structure, demonstrated potent and broad-spectrum antibacterial activity without inducing hemolytic or cytotoxic effects. RI-18 effectively eliminated planktonic and biofilm-associated bacteria as well as persister cells and exhibited high bacterial plasma membrane affinity, inducing rapid membrane permeabilization and rupture. Notably, RI-18 significantly reduced bacterial loads without promoting bacterial resistance, highlighting its therapeutic potential. Overall, this study identified RI-18 as a promising antimicrobial candidate. The rational strategy of tuning imperfect amphipathicity to enhance the AMP activity and selectivity may facilitate the design and development of AMPs.

Chromosome-level genome assembly and annotation of the crested gecko, Correlophus ciliatus, a lizard incapable of tail regeneration

The vast majority of gecko species are capable of tail regeneration, but singular geckos of Correlophus, Uroplatus, and Nephrurus genera are unable to regrow lost tails. Of these non-regenerative geckos, the crested gecko (Correlophus ciliatus) is distinguished by ready availability, ease of care, high productivity, and hybridization potential. These features make C. ciliatus particularly suited as a model for studying the genetic, molecular, and cellular mechanisms underlying loss of tail regeneration capabilities. We report a contiguous genome of C. ciliatus with a total size of 1.65 Gb, 152 scaffolds, L50 of 6, and N50 of 109 Mb. Repetitive content consists of 40.41% of the genome, and a total of 30,780 genes were annotated. Our assembly of the crested gecko genome provides a valuable resource for future comparative genomic studies between non-regenerative and regenerative geckos and other squamate reptiles.

Findings

We report genome sequencing, assembly, and annotation for the crested gecko, Correlophus ciliatus.

Chromosome-scale genome assembly of Astragalus membranaceus using PacBio and Hi-C technologies

Astragalus membranaceus (Fisch.) Bge (AM) is a medicinal herb plant belonging to the Leguminosae family. In this study, we present a chromosome-scale genome assembly of AM, aiming to enhance the molecular biology and functional studies of Astragali Radix. The genome size of AM is about 1.43 Gb, with a contig N50 value of 1.67 Mb. A total of 98.16% of the assembly anchored to 9 pseudochromosomes using Hi-C technology. The assembly completeness was estimated to be 97.27% using BUSCO with the long terminal repeat assembly index (LAI) of 16.22 and quality value (QV) of 48.58. Additionally, the genome contained 67.98% repetitive sequences. Genome annotation predicted 29,914 protein-coding genes, including 73 genes involved in the flavonoid biosynthetic pathway and 2,048 transcription factors. The high-quality genome assembly and gene annotation resources will greatly facilitate future functional genomic studies in Leguminosae species.

The expression pattern of Wnt6, Wnt10A, and HOXA13 during regenerating tails of Gekko Japonicus

Wnt signal is crucial to correctly regenerate tissues along the original axis in many animals. Lizards are able to regenerate their tails spontaneously, while the anterior-posterior axis information required for the successful regeneration is still elusive. In this study, we investigated the expression pattern of Wnt ligands and HOX genes during regeneration. The results of in situ hybridization revealed that Wnt6 and Wnt10A mRNA levels are higher in wound epithelium (WE) than that in blastema during regeneration. In addition, we showed that Wnt agonist positively regulated the expression of HOXA13 in cultured blastema cells, while did not show similar effect on that of HOXB13, HOXC13 and HOXD13. Finally, we found that HOXA13 showed a gradient level along the anterior-posterior axis of regenerated blastema, with higher level at the caudal end. These data proposed that Wnt6, Wnt10A and HOXA13 might play an important role in establishing distal position for regeneration.

Clinging on the brink: Whole genomes reveal human-induced population declines and severe inbreeding in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus)

Human-mediated habitat destruction has had a profound impact on increased species extinction rates and population declines worldwide. The coastal development in the United Arab Emirates (UAE) over the last two decades, serves as an example of how habitat transformation can alter the landscape of a country in just a few years. Here, we study the genomic implications of habitat transformation in the Critically Endangered Emirati Leaf-toed Gecko (Asaccus caudivolvulus), the only endemic vertebrate of the UAE. We generate a high-quality reference genome for this gecko, representing the first reference genome for the family Phyllodactylidae, and produce whole-genome resequencing data for 23 specimens from 10 different species of leaf-toed geckos. Our results show that A. caudivolvulus has consistently lower genetic diversity than any other Arabian species of Asaccus, suggesting a history of ancient population declines. However, high levels of recent inbreeding are recorded among populations in heavily developed areas, with a more than 50% increase in long runs of homozygosity within a 9-year period. Moreover, results suggest that this species does not effectively purge deleterious mutations, hence making it more vulnerable to future stochastic threats. Overall, results show that A. caudivolvulus is in urgent need of protection, and habitat preservation must be warranted to ensure the species' survival.

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.

Spatio-temporal expression patterns of glycine-rich beta proteins and cysteine-rich beta proteins in setae development of Gekko japonicus

BACKGROUND

Setae on the pad lamellae of the Japanese gecko Gekko japonicus (Schlegel, 1836), a vital epidermal derivative, are primarily composed of cornified beta-proteins (CBPs) and play a pivotal role in adhesion and climbing. The amino acid composition of CBPs might be a determining factor influencing their functional properties. However, the molecular mechanisms governed by CBP genes with diverse amino acid compositions in setae development remain unexplored.

RESULTS

Based on RNA-seq analyses, this study confirmed that all G. japonicus CBPs (GjCBPs) are involved in setae formation. Cysteine-rich CBPs encoding genes (ge-cprp-17 to ge-cprp-26) and glycine-rich CBPs encoding genes (ge-gprp-17 to ge-gprp-22) were haphazardly selected, with quantitative real-time PCR revealing their expression patterns in embryonic pad lamellae and dorsal epidermis. It is inferred that glycine-rich CBPs are integral to the formation of both dorsal scales and lamellar setae, cysteine-rich CBPs are primarily associated with setae development. Additionally, fluorescence in situ hybridization revealed spatiotemporal differences in the expression of a glycine-rich CBP encoding gene (ge-gprp-19) and a cysteine-rich CBP encoding gene (ge-cprp-17) during dorsal scales and/or lamellar development.

CONCLUSIONS

All 66 CBPs are involved in the formation of setae. Glycine-rich CBPs hold a significant role in the development of dorsal scales and lamellar setae, whereas most cysteine-rich CBPs appear to be essential components of G. japonicus setae. Even GjCBPs with similar amino acid compositions may play diverse functions. The clear spatio-temporal expression differences between the glycine-rich and cysteine-rich CBP encoding genes during epidermal scale and/or setae formation were observed. Embryonic developmental stages 39 to 42 emerged as crucial phases for setae development. These findings lay the groundwork for deeper investigation into the function of GjCBPs in the development of G. japonicus setae.

Genome-Wide Identification of Specific Genetic Loci Common to Sheep and Goat

Sheep and goat may become carriers of some zoonotic diseases. They are important livestock and experimental model animals for human beings. The fast and accurate identification of genetic materials originating from sheep and goat can prevent and inhibit the spread of some zoonotic diseases, monitor market product quality, and maintain the stability of animal husbandry and food industries. This study proposed a methodology for identifying sheep and goat common specific sites from a genome-wide perspective. A total of 150 specific sites were selected from three data sources, including the coding sequences of single copy genes from nine species (sheep, goat, cow, pig, dog, horse, human, mouse, and chicken), the dbSNPs for these species, and human 100-way alignment data. These 150 sites exhibited low intraspecific heterogeneity in the resequencing data of 1450 samples from five species (sheep, goat, cow, pig, and chicken) and high interspecific divergence in the human 100-way alignment data after quality control. The results were proven to be reliable at the data level. Using the process proposed in this study, specific sites of other species can be screened, and genome-level species identification can be performed using the screened sites.

Convergent gene losses and pseudogenizations in multiple lineages of stomachless fishes

The regressive evolution of independent lineages often results in convergent phenotypes. Several teleost groups display secondary loss of the stomach, and four gastric genes, atp4a, atp4b, pgc, and pga2 have been co-deleted in agastric (stomachless) fish. Analyses of genotypic convergence among agastric fishes showed that four genes, slc26a9, kcne2, cldn18a, and vsig1, were co-deleted or pseudogenized in most agastric fishes of the four major groups. kcne2 and vsig1 were also deleted or pseudogenized in the agastric monotreme echidna and platypus, respectively. In the stomachs of sticklebacks, these genes are expressed in gastric gland cells or surface epithelial cells. An ohnolog of cldn18 was retained in some agastric teleosts but exhibited an increased non-synonymous substitution when compared with gastric species. These results revealed novel convergent gene losses at multiple loci among the four major groups of agastric fish, as well as a single gene loss in the echidna and platypus.

Convergent mechanism underlying the acquisition of vertebrate scotopic vision

High sensitivity of scotopic vision (vision in dim light conditions) is achieved by the rods' low background noise, which is attributed to a much lower thermal activation rate (kth) of rhodopsin compared with cone pigments. Frogs and nocturnal geckos uniquely possess atypical rods containing noncanonical cone pigments that exhibit low kth, mimicking rhodopsin. Here, we investigated the convergent mechanism underlying the low kth of rhodopsins and noncanonical cone pigments. Our biochemical analysis revealed that the kth of canonical cone pigments depends on their absorption maximum (λmax). However, rhodopsin and noncanonical cone pigments showed a substantially lower kth than predicted from the λmax dependency. Given that the λmax is inversely proportional to the activation energy of the pigments in the Hinshelwood distribution-based model, our findings suggest that rhodopsin and noncanonical cone pigments have convergently acquired low frequency of spontaneous-activation attempts, including thermal fluctuations of the protein moiety, in the molecular evolutionary processes from canonical cone pigments, which contributes to highly sensitive scotopic vision.

Cytogenetic Analysis of Satellitome of Madagascar Leaf-Tailed Geckos

Satellite DNA (satDNA) consists of sequences of DNA that form tandem repetitions across the genome, and it is notorious for its diversity and fast evolutionary rate. Despite its importance, satDNA has been only sporadically studied in reptile lineages. Here, we sequenced genomic DNA and PCR-amplified microdissected W chromosomes on the Illumina platform in order to characterize the monomers of satDNA from the Henkel's leaf-tailed gecko U. henkeli and to compare their topology by in situ hybridization in the karyotypes of the closely related Günther's flat-tail gecko U. guentheri and gold dust day gecko P. laticauda. We identified seventeen different satDNAs; twelve of them seem to accumulate in centromeres, telomeres and/or the W chromosome. Notably, centromeric and telomeric regions seem to share similar types of satDNAs, and we found two that seem to accumulate at both edges of all chromosomes in all three species. We speculate that the long-term stability of all-acrocentric karyotypes in geckos might be explained from the presence of specific satDNAs at the centromeric regions that are strong meiotic drivers, a hypothesis that should be further tested.

Insights into the adaptive evolution of chromosome and essential traits through chromosome-level genome assembly of Gekko japonicus

Gekko japonicus possesses flexible climbing and detoxification abilities under insectivorous habits. Still, the evolutionary mechanisms behind these traits remain unclarified. This study presents a chromosome-level G. japonicus genome, revealing that its evolutionary breakpoint regions were enriched with specific repetitive elements and defense response genes. Gene families unique to G. japonicus and positively selected genes are mainly enriched in immune, sensory, and nervous pathways. Expansion of bitter taste receptor type 2 primarily in insectivorous species could be associated with toxin clearance. Detox cytochrome P450 in G. japonicus has undergone more birth and death processes than biosynthesis-type P450 genes. Proline, cysteine, glycine, and serine in corneous beta proteins of G. japonicus might influence flexibility and setae adhesiveness. Certain thermosensitive transient receptor potential channels under relaxed purifying selection or positive selection in G. japonicus might enhance adaptation to climate change. This genome assembly offers insights into the adaptive evolution of gekkotans.

Development-Associated Genes of the Epidermal Differentiation Complex (EDC)

The epidermal differentiation complex (EDC) is a cluster of genes that encode protein components of the outermost layers of the epidermis in mammals, reptiles and birds. The development of the stratified epidermis from a single-layered ectoderm involves an embryo-specific superficial cell layer, the periderm. An additional layer, the subperiderm, develops in crocodilians and over scutate scales of birds. Here, we review the expression of EDC genes during embryonic development. Several EDC genes are expressed predominantly or exclusively in embryo-specific cell layers, whereas others are confined to the epidermal layers that are maintained in postnatal skin. The S100 fused-type proteins scaffoldin and trichohyalin are expressed in the avian and mammalian periderm, respectively. Scaffoldin forms the so-called periderm granules, which are histological markers of the periderm in birds. Epidermal differentiation cysteine-rich protein (EDCRP) and epidermal differentiation protein containing DPCC motifs (EDDM) are expressed in the avian subperiderm where they are supposed to undergo cross-linking via disulfide bonds. Furthermore, a histidine-rich epidermal differentiation protein and feather-type corneous beta-proteins, also known as beta-keratins, are expressed in the subperiderm. The accumulating evidence for roles of EDC genes in the development of the epidermis has implications on the evolutionary diversification of the skin in amniotes.

Genome-wide characterization of the TGF-β gene family and their expression in different tissues during tail regeneration in the Schlegel's Japanese gecko Gekko japonicus

The transforming growth factor-β (TGF-β) gene family is unique to animals and is involved in various important processes including tissue regeneration. Here, we identified 52 TGF-β family genes based on genome sequences of the gecko (Gekko japonicus), compared TGF-β genes between G. japonicus and other four reptilian species, and evaluated the expression of 14 randomly selected genes in muscle, kidney, liver, heart, and brain during tail regeneration to investigate whether their expression was tissue-dependent. We detected 23 conserved domains, 13 in the TGF-β ligand subfamily, and 10 in the receptor subfamily. The pattern of higher genetic variation in the ligand subfamily than in the receptor subfamily in vertebrates might result from the precise localization of agonists and antagonists in the cell surface and intracellular compartment. TGF-β genes were unevenly distributed across 15 chromosomes in G. japonicus, presumably resulting from gene losses and gains during evolution. Genes in the TGF-β receptor subfamily (ACVR2A, ACVR2B, ACVR1, BMPR1A, ACVRL1, BMPR2 and TGFBR1) played a vital role in the TGF-β signal pathway. The expression of all 14 randomly selected TGF-β genes was tissue-specific. Our study supports the speculation that some TGF-β family genes are involved in the early stages of tail regeneration.

Comparative proteomic analysis of tail regeneration in the green anole lizard, Anolis carolinensis

As amniote vertebrates, lizards are the most closely related organisms to humans capable of appendage regeneration. Lizards can autotomize, or release their tails as a means of predator evasion, and subsequently regenerate a functional replacement. Green anoles (Anolis carolinensis) can regenerate their tails through a process that involves differential expression of hundreds of genes, which has previously been analyzed by transcriptomic and microRNA analysis. To investigate protein expression in regenerating tissue, we performed whole proteomic analysis of regenerating tail tip and base. This is the first proteomic data set available for any anole lizard. We identified a total of 2,646 proteins - 976 proteins only in the regenerating tail base, 796 only in the tail tip, and 874 in both tip and base. For over 90% of these proteins in these tissues, we were able to assign a clear orthology to gene models in either the Ensembl or NCBI databases. For 13 proteins in the tail base, 9 proteins in the tail tip, and 10 proteins in both regions, the gene model in Ensembl and NCBI matched an uncharacterized protein, confirming that these predictions are present in the proteome. Ontology and pathways analysis of proteins expressed in the regenerating tail base identified categories including actin filament-based process, ncRNA metabolism, regulation of phosphatase activity, small GTPase mediated signal transduction, and cellular component organization or biogenesis. Analysis of proteins expressed in the tail tip identified categories including regulation of organelle organization, regulation of protein localization, ubiquitin-dependent protein catabolism, small GTPase mediated signal transduction, morphogenesis of epithelium, and regulation of biological quality. These proteomic findings confirm pathways and gene families activated in tail regeneration in the green anole as well as identify uncharacterized proteins whose role in regrowth remains to be revealed. This study demonstrates the insights that are possible from the integration of proteomic and transcriptomic data in tail regrowth in the green anole, with potentially broader application to studies in other regenerative models.

Scales of non-avian reptiles and their derivatives contain corneous beta proteins coded from genes localized in the Epidermal Differentiation Complex

The evolution of modern reptiles from basic reptilian ancestors gave rise to scaled vertebrates. Scales are of different types, and their corneous layer can shed frequently during the year in lepidosaurians (lizards, snakes), 1-2 times per year in the tuatara and in some freshwater turtle, irregularly in different parts of the body in crocodilians, or simply wore superficially in marine and terrestrial turtles. Lepidosaurians possess tuberculate, non-overlapped or variably overlapped scales with inter-scale (hinge) regions. The latter are hidden underneath the outer scale surface or may be more exposed in specific body areas. Hinge regions allow stretching during growth and movement so that the skin remains mechanically functional. Crocodilian and turtles feature flat and shield scales (scutes) with narrow inter-scale regions for stretching and growth. The epidermis of non-avian reptilian hinge regions is much thinner than the exposed outer surface of scales and is less cornified. Despite the thickness of the epidermis, scales are mainly composed of variably amount of Corneous Beta Proteins (CBPs) that are coded in a gene cluster known as EDC (Epidermal Differentiation Complex). These are small proteins, 100-200 amino acid long of 8-25 kDa, rich in glycine and cysteine but also in serine, proline and valine that participate to the formation of beta-sheets in the internal part of the protein, the beta-region. This region determines the further polymerization of CBPs in filamentous proteins that, together a network of Intermediate Filament Keratins (IFKs) and other minor epidermal proteins from the EDC make the variable pliable or inflexible corneous material of reptilian scales, claws and of turtle beak. The acquisition of scales and skin derivatives with different mechanical and material properties, mainly due to the evolution of reptile CBPs, is essential for the life and different adaptations of these vertebrates.

Partial in vivo reprogramming enables injury-free intestinal regeneration via autonomous Ptgs1 induction

Tissue regeneration after injury involves the dedifferentiation of somatic cells, a natural adaptive reprogramming that leads to the emergence of injury-responsive cells with fetal-like characteristics. However, there is no direct evidence that adaptive reprogramming involves a shared molecular mechanism with direct cellular reprogramming. Here, we induced dedifferentiation of intestinal epithelial cells using OSKM (Oct4, Sox2, Klf4, and c-Myc) in vivo. The OSKM-induced forced dedifferentiation showed similar molecular features of intestinal regeneration, including a transition from homeostatic cell types to injury-responsive-like cell types. These injury-responsive-like cells, sharing gene signatures of revival stem cells and atrophy-induced villus epithelial cells, actively assisted tissue regeneration following damage. In contrast to normal intestinal regeneration involving Ptgs2 induction, the OSKM promotes autonomous production of prostaglandin E2 via epithelial Ptgs1 expression. These results indicate prostaglandin synthesis is a common mechanism for intestinal regeneration but involves a different enzyme when partial reprogramming is applied to the intestinal epithelium.

gjSOX9 Cloning, Expression, and Comparison with gjSOXs Family Members in Gekko japonicus

SOX9 plays a crucial role in the male reproductive system, brain, and kidneys. In this study, we firstly analyzed the complete cDNA sequence and expression patterns for SOX9 from Gekko japonicus SOX9 (gjSOX9), carried out bioinformatic analyses of physiochemical properties, structure, and phylogenetic evolution, and compared these with other members of the gjSOX family. The results indicate that gjSOX9 cDNA comprises 1895 bp with a 1482 bp ORF encoding 494aa. gjSOX9 was not only expressed in various adult tissues but also exhibited a special spatiotemporal expression pattern in gonad tissues. gjSOX9 was predicted to be a hydrophilic nucleoprotein with a characteristic HMG-Box harboring a newly identified unique sequence, "YKYQPRRR", only present in SOXE members. Among the 20 SOX9 orthologs, gjSOX9 shares the closest genetic relationships with Eublepharis macularius SOX9, Sphacrodactylus townsendi SOX9, and Hemicordylus capensis SOX9. gjSOX9 and gjSOX10 possessed identical physicochemical properties and subcellular locations and were tightly clustered with gjSOX8 in the SOXE group. Sixteen gjSOX family members were divided into six groups: SOXB, C, D, E, F, and H with gjSOX8, 9, and 10 in SOXE among 150 SOX homologs. Collectively, the available data in this study not only facilitate a deep exploration of the functions and molecular regulation mechanisms of the gjSOX9 and gjSOX families in G. japonicus but also contribute to basic research regarding the origin and evolution of SOX9 homologs or even sex-determination mode in reptiles.

Evolutionary Insights into the Relationship of Frogs, Salamanders, and Caecilians and Their Adaptive Traits, with an Emphasis on Salamander Regeneration and Longevity

The extant amphibians have developed uncanny abilities to adapt to their environment. I compared the genes of amphibians to those of other vertebrates to investigate the genetic changes underlying their unique traits, especially salamanders' regeneration and longevity. Using the well-supported Batrachia tree, I found that salamander genomes have undergone accelerated adaptive evolution, especially for development-related genes. The group-based comparison showed that several genes are under positive selection, rapid evolution, and unexpected parallel evolution with traits shared by distantly related species, such as the tail-regenerative lizard and the longer-lived naked mole rat. The genes, such as EEF1E1, PAFAH1B1, and OGFR, may be involved in salamander regeneration, as they are involved in the apoptotic process, blastema formation, and cell proliferation, respectively. The genes PCNA and SIRT1 may be involved in extending lifespan, as they are involved in DNA repair and histone modification, respectively. Some genes, such as PCNA and OGFR, have dual roles in regeneration and aging, which suggests that these two processes are interconnected. My experiment validated the time course differential expression pattern of SERPINI1 and OGFR, two genes that have evolved in parallel in salamanders and lizards during the regeneration process of salamander limbs. In addition, I found several candidate genes responsible for frogs' frequent vocalization and caecilians' degenerative vision. This study provides much-needed insights into the processes of regeneration and aging, and the discovery of the critical genes paves the way for further functional analysis, which could open up new avenues for exploiting the genetic potential of humans and improving human well-being.

Genetics and biology of coloration in reptiles: the curious case of the Lemon Frost geckos

Although there are more than 10,000 reptile species, and reptiles have historically contributed to our understanding of biology, genetics research into class Reptilia has lagged compared with other animals. Here, we summarize recent progress in genetics of coloration in reptiles, with a focus on the leopard gecko, Eublepharis macularius. We highlight genetic approaches that have been used to examine variation in color and pattern formation in this species as well as to provide insights into mechanisms underlying skin cancer. We propose that their long breeding history in captivity makes leopard geckos one of the most promising emerging reptilian models for genetic studies. More broadly, technological advances in genetics, genomics, and gene editing may herald a golden era for studies of reptile biology.

Protein extracts from regenerating lizard tail show an inhibitory effect on human cancer cells cultivated in-vitro

BACKGROUND

accumulating evidence indicates that during tail regeneration in lizards the initial stage of regenerative blastema is a tumor-like proliferative outgrowth that rapidly elongates into a new tail composed of fully differentiated tissues. Both oncogenes and tumor-suppressors are expressed during regeneration, and it has been hypothesized that an efficient control of cell proliferation avoids that the blastema is turned into a tumor outgrowth.

METHODS

in order to determine whether functional tumor-suppressors are present in the growing blastema we have utilized protein extracts collected from early regenerating tails of 3-5 mm that have been tested for a potential anti-tumor effect on in-vitro culture by using cancer cell lines from human mammary gland (MDA-MB-231) and prostate cancer (DU145).

RESULTS

at specific dilutions, the extract determines a reduction of viability in cancer cells after 2-4 days of culture, as supported by statistical and morphological analyses. While control cells appear viable, treated cells result damaged and produce an intense cytoplasmic granulation and degeneration.

CONCLUSIONS

this negative effect on cell viability and proliferation is absent using tissues from the original tail supporting the hypothesis that only regenerating tissues synthesize tumor-suppressor molecules. The study suggests that the regenerating tail of lizard at the stages here selected contains some molecules that determine inhibition of cell viability on the cancer cells analyzed.

Super-high procoagulant activity of gecko thrombin: A gift from sky dragon

AIMS

Gecko, the "sky dragon" named by Traditional Chinese Medicine, undergoes rapid coagulation and scarless regeneration following tail amputation in the natural ecology, providing a perfect opportunity to develop the efficient and safe drug for blood clotting. Here, gecko thrombin (gthrombin) was recombinantly prepared and comparatively studied on its procoagulant activity.

METHODS

The 3D structure of gthrombin was constructed using the homology modeling method of I-TASSER. The active gthrombin was prepared by the expression of gecko prethrombin-2 in 293 T cells, followed by purification with Ni2+ -chelating column chromatography prior to activation by snake venom-derived Ecarin. The enzymatic activities of gthrombin were assayed by hydrolysis of synthetic substrate S-2238 and the fibrinogen clotting. The vulnerable nerve cells were used to evaluate the toxicity of gthrombin at molecular and cellular levels.

RESULTS

The active recombinant gthrombin showed super-high catalytic and fibrinogenolytic efficiency than those of human under different temperatures and pH conditions. In addition, gthrombin made nontoxic effects on the central nerve cells including neurons, contrary to those of mammalian counterparts, which contribute to neuronal damage, astrogliosis, and demyelination.

CONCLUSIONS

A super-high activity but safe procoagulant candidate drug was identified from reptiles, which provided a promising perspective for clinical application in rapid blood clotting.

Expression analysis of alpha keratins and corneous beta-protein genes during embryonic development of Gekko japonicus

Epidermal appendages of birds and reptiles, including claws, feathers, scales, and setae, are primarily composed of alpha keratins (KRTs) and corneous beta-proteins (CBPs). A comprehensive and systematic knowledge of KRTs and CBPs in Schlegel's Japanese gecko (Gekko japonicus) is still lacking. In this study, 22 candidate Gecko japonicus keratin (GjKRT) family genes (12 type I genes, 10 type II genes) were identified in the G. japonicus genome. The majority of GjKRT genes across various subgroups had undergone a prolonged and highly conservative evolutionary process. Through a combination of morphological observation, RNA-seq analysis, and qRT-PCR assay, it was possible to discern the dynamic alterations in the expression of GjKRTs and Gecko japonicus corneous beta-proteins genes (GjCBPs). These findings strongly indicate that GjKRTs gradually accumulate to constitute an α-layer, which is subsequently succeeded by the formation of the corneous beta layer containing GjCBPs at late stages (40-42) of embryonic development. The epidermal appendages in G. japonicus may result from the joint accumulation of KRTs and CBPs, with stages 40-42 being critical for their development. These findings provide novel insights into KRTs and CBPs of G. japonicus and offer a foundation for investigating the functions of GjKRT and GjCBP gene families. Furthermore, this knowledge contributes to unraveling the molecular mechanisms underlying the formation of epidermal appendages in G. japonicus.

Assembly of the largest squamate reference genome to date: The western fence lizard, Sceloporus occidentalis

Spiny lizards (genus Sceloporus) have long served as important systems for studies of behavior, thermal physiology, dietary ecology, vector biology, speciation, and biogeography. The western fence lizard, Sceloporus occidentalis, is found across most of the major biogeographical regions in the western United States and northern Baja California, Mexico, inhabiting a wide range of habitats, from grassland to chaparral to open woodlands. As small ectotherms, Sceloporus lizards are particularly vulnerable to climate change, and S. occidentalis has also become an important system for studying the impacts of land use change and urbanization on small vertebrates. Here, we report a new reference genome assembly for S. occidentalis, as part of the California Conservation Genomics Project (CCGP). Consistent with the reference genomics strategy of the CCGP, we used Pacific Biosciences HiFi long reads and Hi-C chromatin-proximity sequencing technology to produce a de novo assembled genome. The assembly comprises a total of 608 scaffolds spanning 2,856 Mb, has a contig N50 of 18.9 Mb, a scaffold N50 of 98.4 Mb, and BUSCO completeness score of 98.1% based on the tetrapod gene set. This reference genome will be valuable for understanding ecological and evolutionary dynamics in S. occidentalis, the species status of the California endemic island fence lizard (S. becki), and the spectacular radiation of Sceloporus lizards.

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.

Dietary and Sexual Correlates of Gut Microbiota in the Japanese Gecko, Gekko japonicus (Schlegel, 1836)

Numerous studies have demonstrated that multiple intrinsic and extrinsic factors shape the structure and composition of gut microbiota in a host. The disorder of the gut microbiota may trigger various host diseases. Here, we collected fecal samples from wild-caught Japanese geckos (Gekko japonicus) and captive conspecifics fed with mealworms (mealworm-fed geckos) and fruit flies (fly-fed geckos), aiming to examine the dietary and sexual correlates of the gut microbiota. We used 16S rRNA gene sequencing technology to determine the composition of the gut microbiota. The dominant phyla with a mean relative abundance higher than 10% were Verrucomicrobiota, Bacteroidota, and Firmicutes. Gut microbial community richness and diversity were higher in mealworm-fed geckos than in wild geckos. Neither community evenness nor beta diversity of gut microbiota differed among wild, mealworm-fed, and fly-fed geckos. The beta rather than alpha diversity of gut microbiota was sex dependent. Based on the relative abundance of gut bacteria and their gene functions, we concluded that gut microbiota contributed more significantly to the host's metabolic and immune functions. A higher diversity of gut microbiota in mealworm-fed geckos could result from higher chitin content in insects of the order Coleoptera. This study not only provides basic information about the gut microbiota of G. japonicus but also shows that gut microbiota correlates with dietary habits and sex in the species.

Genomes of two Extinct-in-the-Wild reptiles from Christmas Island reveal distinct evolutionary histories and conservation insights

Genomics can play important roles in biodiversity conservation, especially for Extinct-in-the-Wild species where genetic factors greatly influence risk of total extinction and probability of successful reintroductions. The Christmas Island blue-tailed skink (Cryptoblepharus egeriae) and Lister's gecko (Lepidodactylus listeri) are two endemic reptile species that went extinct in the wild shortly after the introduction of a predatory snake. After a decade of management, captive populations have expanded from 66 skinks and 43 geckos to several thousand individuals; however, little is known about patterns of genetic variation in these species. Here, we use PacBio HiFi long-read and Hi-C sequencing to generate highly contiguous reference genomes for both reptiles, including the XY chromosome pair in the skink. We then analyse patterns of genetic diversity to infer ancient demography and more recent histories of inbreeding. We observe high genome-wide heterozygosity in the skink (0.007 heterozygous sites per base-pair) and gecko (0.005), consistent with large historical population sizes. However, nearly 10% of the blue-tailed skink reference genome falls within long (>1 Mb) runs of homozygosity (ROH), resulting in homozygosity at all major histocompatibility complex (MHC) loci. In contrast, we detect a single ROH in Lister's gecko. We infer from the ROH lengths that related skinks may have established the captive populations. Despite a shared recent extinction in the wild, our results suggest important differences in these species' histories and implications for management. We show how reference genomes can contribute evolutionary and conservation insights, and we provide resources for future population-level and comparative genomic studies in reptiles.

Characterization of Two Transposable Elements and an Ultra-Conserved Element Isolated in the Genome of Zootoca vivipara (Squamata, Lacertidae)

Transposable elements (TEs) constitute a considerable fraction of eukaryote genomes representing a major source of genetic variability. We describe two DNA sequences isolated in the lizard Zootoca vivipara, here named Zv516 and Zv817. Both sequences are single-copy nuclear sequences, including a truncation of two transposable elements (TEs), SINE Squam1 in Zv516 and a Tc1/Mariner-like DNA transposon in Zv817. FISH analyses with Zv516 showed the occurrence of interspersed signals of the SINE Squam1 sequence on all chromosomes of Z. vivipara and quantitative dot blot indicated that this TE is present with about 4700 copies in the Z. vivipara genome. FISH and dot blot with Zv817 did not produce clear hybridization signals. Bioinformatic analysis showed the presence of active SINE Squam 1 copies in the genome of different lacertids, in different mRNAs, and intronic and coding regions of various genes. The Tc1/Mariner-like DNA transposon occurs in all reptiles, excluding Sphenodon and Archosauria. Zv817 includes a trait of 284 bp, representing an amniote ultra-conserved element (UCE). Using amniote UCE homologous sequences from available whole genome sequences of major amniote taxonomic groups, we performed a phylogenetic analysis which retrieved Prototheria as the sister group of Metatheria and Eutheria. Within diapsids, Testudines are the sister group to Aves + Crocodylia (Archosauria), and Sphenodon is the sister group to Squamata. Furthermore, large trait regions flanking the UCE are conserved at family level.

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.

Madagascar Leaf-Tail Geckos (Uroplatus spp.) Share Independently Evolved Differentiated ZZ/ZW Sex Chromosomes

Geckos are an excellent group to study the evolution of sex determination, as they possess a remarkable variability ranging from a complete absence of sex chromosomes to highly differentiated sex chromosomes. We explored sex determination in the Madagascar leaf-tail geckos of the genus Uroplatus. The cytogenetic analyses revealed highly heterochromatic W chromosomes in all three examined species (Uroplatus henkeli, U. alluaudi, U. sikorae). The comparative gene coverage analysis between sexes in U. henkeli uncovered an extensive Z-specific region, with a gene content shared with the chicken chromosomes 8, 20, 26 and 28. The genomic region homologous to chicken chromosome 28 has been independently co-opted for the role of sex chromosomes in several vertebrate lineages, including monitors, beaded lizards and monotremes, perhaps because it contains the amh gene, whose homologs were repeatedly recruited as a sex-determining locus. We demonstrate that all tested species of leaf-tail geckos share homologous sex chromosomes despite the differences in shape and size of their W chromosomes, which are not homologous to the sex chromosomes of other closely related genera. The rather old (at least 40 million years), highly differentiated sex chromosomes of Uroplatus geckos can serve as a great system to study the convergence of sex chromosomes evolved from the same genomic region.

Spontaneous neuronal regeneration in the forebrain of the leopard gecko (Eublepharis macularius) following neurochemical lesioning

BACKGROUND

Neurogenesis is the ability to generate new neurons from resident stem/progenitor populations. Although often understood as a homeostatic process, several species of teleost fish, salamanders, and lacertid lizards are also capable of reactive neurogenesis, spontaneously replacing lost or damaged neurons. Here, we demonstrate that reactive neurogenesis also occurs in a distantly related lizard species, Eublepharis macularius, the leopard gecko.

RESULTS

To initiate reactive neurogenesis, the antimetabolite 3-acetylpyridine (3-AP) was administered. Four days following 3-AP administration there is a surge in neuronal cell death within a region of the forebrain known as the medial cortex (homolog of the mammalian hippocampal formation). Neuronal cell death is accompanied by a shift in resident microglial morphology and an increase neural stem/progenitor cell proliferation. By 30 days following 3-AP administration, the medial cortex was entirely repopulated by NeuN+ neurons. At the same time, local microglia have reverted to a resting state and cell proliferation by neural stem/progenitors has returned to levels comparable with uninjured controls.

CONCLUSIONS

Together, these data provide compelling evidence of reactive neurogenesis in leopard geckos, and indicate that the ability of lizards to spontaneously replace lost or damaged forebrain neurons is more taxonomically widespread and evolutionarily conserved than previously considered.

The evolutionary history and spectral tuning of vertebrate visual opsins

Many animals depend on the sense of vision for survival. In eumetazoans, vision requires specialized, light-sensitive cells called photoreceptors. Light reaches the photoreceptors and triggers the excitation of light-detecting proteins called opsins. Here, we describe the story of visual opsin evolution from the ancestral bilaterian to the extant vertebrate lineages. We explain the mechanisms determining color vision of extant vertebrates, focusing on opsin gene losses, duplications, and the expression regulation of vertebrate opsins. We describe the sequence variation both within and between species that has tweaked the sensitivities of opsin proteins towards different wavelengths of light. We provide an extensive resource of wavelength sensitivities and mutations that have diverged light sensitivity in many vertebrate species and predict how these mutations were accumulated in each lineage based on parsimony. We suggest possible natural and sexual selection mechanisms underlying these spectral differences. Understanding how molecular changes allow for functional adaptation of animals to different environments is a major goal in the field, and therefore identifying mutations affecting vision and their relationship to photic selection pressures is imperative. The goal of this review is to provide a comprehensive overview of our current understanding of opsin evolution in vertebrates.

Can extreme climatic events induce shifts in adaptive potential? A conceptual framework and empirical test with Anolis lizards

Multivariate adaptation to climatic shifts may be limited by trait integration that causes genetic variation to be low in the direction of selection. However, strong episodes of selection induced by extreme climatic pressures may facilitate future population-wide responses if selection reduces trait integration and increases adaptive potential (i.e., evolvability). We explain this counter-intuitive framework for extreme climatic events in which directional selection leads to increased evolvability and exemplify its use in a case study. We tested this hypothesis in two populations of the lizard Anolis scriptus that experienced hurricane-induced selection on limb traits. We surveyed populations immediately before and after the hurricane as well as the offspring of post-hurricane survivors, allowing us to estimate both selection and response to selection on key functional traits: forelimb length, hindlimb length, and toepad area. The direct selection was parallel in both islands and strong in several limb traits. Even though overall limb integration did not change after the hurricane, both populations showed a non-significant tendency toward increased evolvability after the hurricane despite the direction of selection not being aligned with the axis of most variance (i.e., body size). The population with comparably lower between-limb integration showed a less constrained response to selection. Hurricane-induced selection, not aligned with the pattern of high trait correlations, likely conflicts with selection occurring during normal ecological conditions that favours functional coordination between limb traits, and would likely need to be very strong and more persistent to elicit a greater change in trait integration and evolvability. Future tests of this hypothesis should use G-matrices in a variety of wild organisms experiencing selection due to extreme climatic events.

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.

Akt/mTOR integrate energy metabolism with Wnt signal to influence wound epithelium growth in Gekko Japonicus

The formation of wound epithelium initiates regeneration of amputated tail in Gekko japonicus. Energy metabolism is indispensable for the growth of living creatures and typically influenced by temperature. In this study, we reveal that low temperature lowers energy metabolism level and inhibits the regeneration of amputated tails of Gekko japonicus. We further find that low temperature attenuates the activation of protein kinase B (Akt) and mammalian target of rapamycin (mTOR) in regenerated tissues upon injury signals, and the inhibition of Akt hinders proliferation of the wound epithelium. Additionally, wingless/integrated (Wnt) inhibition suppresses epithelium proliferation and formation by inhibiting Akt activation. Finally, low temperature elevates the activity of adenylate-activated kinase (AMPK) pathway and in turn attenuates wound epithelium formation. Meanwhile, either mTOR downregulation or AMPK upregulation is associated with worse wound epithelium formation. Summarily, low temperature restricts wound epithelium formation by influencing energy sensory pathways including Akt/mTOR and AMPK signaling, which is also modulated by injury induced Wnt signal. Our results provide a mechanism that incorporates the injury signals with metabolic pathway to facilitate regeneration.

Low temperature inhibits the regeneration of amputated tails of Gekko japonicus by influencing the energy sensory Akt/mTOR pathway, which is also modulated by injury-induced Wnt signal.

Chromosome-Level Genome Assembly Reveals Dynamic Sex Chromosomes in Neotropical Leaf-Litter Geckos (Sphaerodactylidae: Sphaerodactylus)

Sex determination is a critical element of successful vertebrate development, suggesting that sex chromosome systems might be evolutionarily stable across lineages. For example, mammals and birds have maintained conserved sex chromosome systems over long evolutionary time periods. Other vertebrates, in contrast, have undergone frequent sex chromosome transitions, which is even more amazing considering we still know comparatively little across large swaths of their respective phylogenies. One reptile group in particular, the gecko lizards (infraorder Gekkota), shows an exceptional lability with regard to sex chromosome transitions and may possess the majority of transitions within squamates (lizards and snakes). However, detailed genomic and cytogenetic information about sex chromosomes is lacking for most gecko species, leaving large gaps in our understanding of the evolutionary processes at play. To address this, we assembled a chromosome-level genome for a gecko (Sphaerodactylidae: Sphaerodactylus) and used this assembly to search for sex chromosomes among six closely related species using a variety of genomic data, including whole-genome re-sequencing, RADseq, and RNAseq. Previous work has identified XY systems in two species of Sphaerodactylus geckos. We expand upon that work to identify between two and four sex chromosome cis-transitions (XY to a new XY) within the genus. Interestingly, we confirmed two different linkage groups as XY sex chromosome systems that were previously unknown to act as sex chromosomes in tetrapods (syntenic with Gallus chromosome 3 and Gallus chromosomes 18/30/33), further highlighting a unique and fascinating trend that most linkage groups have the potential to act as sex chromosomes in squamates.

Genomic adaptations for arboreal locomotion in Asian flying treefrogs

SignificanceTo adapt to arboreal lifestyles, treefrogs have evolved a suite of complex traits that support vertical movement and gliding, thus presenting a unique case for studying the genetic basis for traits causally linked to vertical niche expansion. Here, based on two de novo-assembled Asian treefrog genomes, we determined that genes involved in limb development and keratin cytoskeleton likely played a role in the evolution of their climbing systems. Behavioral and morphological evaluation and time-ordered gene coexpression network analysis revealed the developmental patterns and regulatory pathways of the webbed feet used for gliding in Rhacophorus kio.

Immunolocalization of Adenomatous Polyposis Coli protein (apc) in the regenerating lizard tail suggests involvement in tissue differentiation and regulation of growth

Lizard tail regeneration is likely regulated by the balanced activity of oncogenes and tumor suppressors that control cell proliferation avoiding tumorigenic degeneration. One of the main tumor suppressor genes present in the regenerating tail is the "adenomatous polyposis coli (apc)" but the localization of its coded protein (apc) is not known. This protein may be involved in regulation of apical-basal tail regeneration in lizards. The present immunohistochemical study shows that apc is localized in apical wound epidermis and regenerating ependyme, two tissues that proliferate and also express onco-genes. Apc is not present in blastema cells but localizes in differentiating cells of regenerating scales, muscles and less intensely in the non-apical ependymal epithelium and cartilage. This suggests that apc is involved in the induction of their differentiation. The apc immunolabeling is mainly nuclear in the basal epidermal layer of the apical wound epidermis where it may be involved in modulating keratinocytes proliferation, like in the forming scales. In regenerating muscle and cartilage apc is mainly cytoplasmic while sparse labeled nuclei are seen in proliferative areas of these tissues. In the regenerating spinal cord, the nuclear and cytoplasmic apc labeling is present in ependymal cells of the distal-most ependymal ampulla but the labeling fades in more proximal regions and mainly remains in the cytoplasm facing the central canal and in sparse nuclei. It is suggested that the pattern of immunolabeling for apc indicates that this tumor suppressor may contribute to tissue differentiation within the regenerating tail. This article is protected by copyright. All rights reserved.

Biomimetic fracture model of lizard tail autotomy

Lizard tail autotomy is an antipredator strategy consisting of sturdy attachment at regular times but quick detachment during need. We propose a biomimetic fracture model of lizard tail autotomy using multiscale hierarchical structures. The structures consist of uniformly distributed micropillars with nanoporous tops, which recapitulate the high-density mushroom-shaped microstructures found on the lizard tail's muscle fracture plane. The biomimetic experiments showed adhesion enhancement when combining nanoporous interfacial surfaces with flexible micropillars in tensile and peel modes. The fracture modeling identified micro- and nanostructure-based toughening mechanisms as the critical factor. Under wet conditions, capillarity-assisted energy dissipation pertaining to liquid-filled microgaps and nanopores further increased the adhesion performance. This research presents insights on lizard tail autotomy and provides new biomimetic ideas to solve adhesion problems.

A chromosome-level genome assembly and annotation of the desert horned lizard, Phrynosoma platyrhinos, provides insight into chromosomal rearrangements among reptiles

BACKGROUND

The increasing number of chromosome-level genome assemblies has advanced our knowledge and understanding of macroevolutionary processes. Here, we introduce the genome of the desert horned lizard, Phrynosoma platyrhinos, an iguanid lizard occupying extreme desert conditions of the American southwest. We conduct analysis of the chromosomal structure and composition of this species and compare these features across genomes of 12 other reptiles (5 species of lizards, 3 snakes, 3 turtles, and 1 bird).

FINDINGS

The desert horned lizard genome was sequenced using Illumina paired-end reads and assembled and scaffolded using Dovetail Genomics Hi-C and Chicago long-range contact data. The resulting genome assembly has a total length of 1,901.85 Mb, scaffold N50 length of 273.213 Mb, and includes 5,294 scaffolds. The chromosome-level assembly is composed of 6 macrochromosomes and 11 microchromosomes. A total of 20,764 genes were annotated in the assembly. GC content and gene density are higher for microchromosomes than macrochromosomes, while repeat element distributions show the opposite trend. Pathway analyses provide preliminary evidence that microchromosome and macrochromosome gene content are functionally distinct. Synteny analysis indicates that large microchromosome blocks are conserved among closely related species, whereas macrochromosomes show evidence of frequent fusion and fission events among reptiles, even between closely related species.

CONCLUSIONS

Our results demonstrate dynamic karyotypic evolution across Reptilia, with frequent inferred splits, fusions, and rearrangements that have resulted in shuffling of chromosomal blocks between macrochromosomes and microchromosomes. Our analyses also provide new evidence for distinct gene content and chromosomal structure between microchromosomes and macrochromosomes within reptiles.