Madagascar ground gecko genome analysis characterizes asymmetric fates of duplicated genes

Identification of apolipoprotein E-derived amyloid within cholesterol granulomas of leopard geckos (Eublepharis macularius)

Apolipoprotein E (ApoE) is involved in cholesterol transport among cells and also plays an important role in amyloid formation, co-depositing with amyloid fibrils in various types of amyloidosis. Although the in vivo amyloidogenicity of ApoE has not been previously demonstrated, this study provides evidence of ApoE amyloidogenicity in leopard geckos (Eublepharis macularius), belonging to the class Reptilia. Histologically, amyloid deposits were localized within cholesterol granulomas and exhibited positive Congo red staining, with yellow to green birefringence under polarized light. On mass spectrometry-based proteomic analysis, ApoE was detected as a dominant component of amyloid; of the full length of the 274 amino acid residues, peptides derived from Leu185-Arg230 were frequently detected with non-tryptic truncations. Immunohistochemistry with anti-leopard gecko ApoE antibody showed positive reactions of amyloid deposits. These results show that ApoE is an amyloid precursor protein within the cholesterol granulomas of leopard geckos. Although further investigations are needed, the C-terminal region of ApoE involved in amyloid formation is a lipid-binding region, and there should be a relationship between amyloidogenesis and the development of cholesterol granulomas in leopard geckos. This study provides novel insights into the pathogenesis of ApoE-related diseases.

Comparative Cytogenetics of the Malagasy Ground Geckos of the Paroedura bastardi and Paroedura picta Species Groups

We present a comparative chromosome study of several taxa of the Malagasy ground geckos of the Paroedura bastardi and P. picta species groups. We employed a preliminary molecular analysis using a trait of the mitochondrial 16S rRNA gene (of about 570 bp) to assess the taxonomic status of the samples studied and a cytogenetic analysis with standard karyotyping (5% Giemsa solution), silver staining (Ag-NOR staining) and sequential C-banding (C-banding + Giemsa and + fluorochromes). Our results show that all the taxa studied of the P. bastardi group (P. ibityensis, P. rennerae and P. cf. guibeae) have a similar karyotype composed of 2n = 34 chromosomes, with two metacentric pairs (1 and 3) and all other pairs being acrocentric. Chromosome diversification in the P. bastardi group was mainly linked to the diversification of heteromorphic sex chromosome systems (ZZ/ZW) in P. ibityensis and P. rennerae, while no heteromorphic sex chromosome pair was found in P. cf. guibeae. The two taxa investigated of the P. picta species group (here named P. picta and P. cf. picta based on molecular data) showed the same chromosome number of 2n = 36, mostly acrocentric elements, but differed in the number of metacentric elements, probably as a result of an inversion at chromosome pair 2. We highlight that the genus Paroedura is characterized by the independent diversification of heterogametic sex chromosomes in different evolutionary lineages and, similarly to other phylogenetically related gecko genera, by a progressive formation of a biarmed element by means of tandem fusions and inversions of distinct pairs.

ETV2 induces endothelial, but not hematopoietic, lineage specification in birds

Cardiovascular system develops from the lateral plate mesoderm. Its three primary cell lineages (hematopoietic, endothelial, and muscular) are specified by the sequential actions of conserved transcriptional factors. ETV2, a master regulator of mammalian hemangioblast development, however, is absent in the chicken genome and acts downstream of NPAS4L in zebrafish. Here, we investigated the epistatic relationship between NPAS4L and ETV2 in avian hemangioblast development. We showed that ETV2 is deleted in all 363 avian genomes analyzed. Mouse ETV2 induced LMO2, but not NPAS4L or SCL, expression in chicken mesoderm. Squamate (lizards, geckos, and snakes) genomes contain both NPAS4L and ETV2 In Madagascar ground gecko, both genes were expressed in developing hemangioblasts. Gecko ETV2 induced only LMO2 in chicken mesoderm. We propose that both NPAS4L and ETV2 were present in ancestral amniote, with ETV2 acting downstream of NPAS4L in endothelial lineage specification. ETV2 may have acted as a pioneer factor by promoting chromatin accessibility of endothelial-specific genes and, in parallel with NPAS4L loss in ancestral mammals, has gained similar function in regulating blood-specific genes.

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.

A vertebrate-wide catalogue of T1R receptors reveals diversity in taste perception

Taste is a vital chemical sense for feeding behaviour. In mammals, the umami and sweet taste receptors comprise three members of the taste receptor type 1 (T1R/TAS1R) family: T1R1, T1R2 and T1R3. Because their functional homologues exist in teleosts, only three TAS1R genes generated by gene duplication are believed to have been inherited from the common ancestor of bony vertebrates. Here, we report five previously uncharacterized TAS1R members in vertebrates, TAS1R4, TAS1R5, TAS1R6, TAS1R7 and TAS1R8, based on genome-wide survey of diverse taxa. We show that mammalian and teleost fish TAS1R2 and TAS1R3 genes are paralogues. Our phylogenetic analysis suggests that the bony vertebrate ancestor had nine TAS1Rs resulting from multiple gene duplications. Some TAS1Rs were lost independently in descendent lineages resulting in retention of only three TAS1Rs in mammals and teleosts. Combining functional assays and expression analysis of non-teleost fishes we show that the novel T1Rs form heterodimers in taste-receptor cells and recognize a broad range of ligands such as essential amino acids, including branched-chain amino acids, which have not been previously considered as T1R ligands. This study reveals diversity of taste sensations in both modern vertebrates and their ancestors, which might have enabled vertebrates to adapt to diverse habitats on Earth.

De Novo Whole Genome Assemblies for Two Southern African Dwarf Chameleons (Bradypodion, Chamaeleonidae)

A complete and high-quality reference genome has become a fundamental tool for the study of functional, comparative, and evolutionary genomics. However, efforts to produce high-quality genomes for African taxa are lagging given the limited access to sufficient resources and technologies. The southern African dwarf chameleons (Bradypodion) are a relatively young lineage, with a large body of evidence demonstrating the highly adaptive capacity of these lizards. Bradypodion are known for their habitat specialization, with evidence of convergent phenotypes across the phylogeny. However, the underlying genetic architecture of these phenotypes remains unknown for Bradypodion, and without adequate genomic resources, many evolutionary questions cannot be answered. We present de novo assembled whole genomes for Bradypodion pumilum and Bradypodion ventrale, using Pacific Biosciences long-read sequencing data. BUSCO analysis revealed that 96.36% of single copy orthologs were present in the B. pumilum genome and 94% in B. ventrale. Moreover, these genomes boast scaffold N50 of 389.6 and 374.9 Mb, respectively. Based on a whole genome alignment of both Bradypodion genomes, B. pumilum is highly syntenic with B. ventrale. Furthermore, Bradypodion is also syntenic with Anolis lizards, despite the divergence between these lineages estimated to be nearly 170 Ma. Coalescent analysis of the genomic data also suggests that historical changes in effective population size for these species correspond to notable shifts in the southern African environment. These high-quality Bradypodion genome assemblies will support future research on the evolutionary history, diversification, and genetic underpinnings of adaptation in Bradypodion.

Functional evolutionary convergence of long noncoding RNAs involved in embryonic development

Long noncoding RNAs have been identified in most vertebrates, but the functional characterization of these molecules is challenging, mainly due to the lack of linear sequence homology between species. In this work, we aimed to find functional evolutionary convergent lncRNAs involved in development by screening of k-mer content (nonlinear similarity) and secondary structure-based approaches combining in silico, in vitro and in vivo validation analysis. From the Madagascar gecko genes, we have found a non-orthologous lncRNA with a similar k-mer content and structurally concordant with the human lncRNA EVX1AS. Analysis of function-related characteristics together with locus-specific targeting of human EVX1AS and gecko EVX1AS-like (i.e., CRISPR Display) in human neuroepithelial cells and chicken mesencephalon have confirmed that gecko EVX1AS-like lncRNA mimics human EVX1AS function and induces EVX1 expression independently of the target species. Our data shows functional convergence of non-homologous lncRNAs and presents a useful approach for the definition and manipulation of lncRNA function within different model organisms.

The impact of local genomic properties on the evolutionary fate of genes

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

A reverse genetic approach in geckos with the CRISPR/Cas9 system by oocyte microinjection

Reptiles are important model organisms in developmental and evolutionary biology, but are used less widely than other amniotes such as mouse and chicken. One of the main reasons for this is that has proven difficult to conduct CRISPR/Cas9-mediated genome editing in many reptile species despite the widespread use of this technology in other taxa. Certain features of reptile reproductive systems make it difficult to access one-cell or early-stage zygotes, which represents a key impediment to gene editing techniques. Recently, Rasys and colleagues reported a genome editing method using oocyte microinjection that allowed them to produce genome-edited Anolis lizards. This method opened a new avenue to reverse genetics studies in reptiles. In the present article, we report the development of a related method for genome editing in the Madagascar ground gecko (Paroedura picta), a well-established experimental model, and describe the generation of Tyr and Fgf10 gene-knockout geckos in the F0 generation.

Evolution of Vertebrate Hormones and Their Receptors: Insights from Non-Osteichthyan Genomes

Homeostatic control and reproductive functions of humans are regulated at the molecular levels largely by peptide hormones secreted from endocrine and/or neuroendocrine cells in the central nervous system and peripheral organs. Homologs of those hormones and their receptors function similarly in many vertebrate species distantly related to humans, but the evolutionary history of the endocrine system involving those factors has been obscured by the scarcity of genome DNA sequence information of some taxa that potentially contain their orthologs. Focusing on non-osteichthyan vertebrates, namely jawless and cartilaginous fishes, this article illustrates how investigating genome sequence information assists our understanding of the diversification of vertebrate gene repertoires in four broad themes: (a) the presence or absence of genes, (b) multiplication and maintenance of paralogs, (c) differential fates of duplicated paralogs, and (d) the evolutionary timing of gene origins.

Chromosome-Level Genome Assembly of the Cape Cliff Lizard (Hemicordylus capensis)

Squamates represent a highly diverse and species-rich vertebrate group that is remarkably understudied from a genomic perspective. A scarcity of genomic data is particularly evident for scincomorph lizards, which encompass over 10% of all living squamates, and for which high-quality genomic resources are currently lacking. To address this knowledge gap, we present the first chromosome-level reference genome for this group, generated from a male Cape cliff lizard (Hemicordylus capensis), using highly accurate PacBio HiFi long-read sequencing data, long-range Omni-C chromosomal conformation capture data and transcriptomic data for annotation. The rHemCap1.1 genome assembly spans 2.29 Gb, with a scaffold N50 of 359.65 Mb, and includes 25,300 protein-coding genes, with a BUSCO completeness score of 95.5% (sauropsida_odb10). We have generated the most contiguous and complete chromosome-level squamate reference genome assembly publicly available to date. Furthermore, we used short-read resequencing of 35 males and females and applied a differential coverage approach to infer the sex-determination system of the species, which was previously unknown. Our results suggest this species has XX/XY sex chromosomes, representing the first evidence of sex determination in the family Cordylidae. This reference genome will help to establish this species as an evolutionary model for studying variation in body armor, a key trait in cordylids and other squamate groups. Lastly, this is the first squamate reference genome from a continental African species and, as such, represents a valuable resource not only for further evolutionary research in cordylids but also in closely related groups.

Left-right asymmetric expression of the Nodal-Lefty-Pitx2 module in developing turtle forebrain

The epithalamus of zebrafish shows morphological and molecular left-right (L-R) asymmetry, but such asymmetry is not apparent in tetrapods. To provide further insight into the evolutionary diversity of brain L-R asymmetry, we have now examined the developing brains of reptile embryos for expression of Nodal, Lefty, and Pitx2. Two turtle species, the Chinese softshell turtle and the red-eared slider turtle, showed left-sided expression of these three genes in the developing forebrain, with this expression occurring after Nodal expression at the lateral plate and the L-R organizer has disappeared. Nodal activity, as revealed by the detection of phosphorylated Smad2/3, was also apparent in the neural epithelium on the left side in both turtle species. In the Chinese softshell turtle, the habenula did not show apparent asymmetry in size and the parapineal organ was absent, but the expression of Kctd12 in the habenula showed a small yet reproducible asymmetry. In contrast to the turtles, L-R asymmetric expression of Nodal, Lefty, Pitx2, or Kctd12 was not detected in the developing brain of the Madagascar ground gecko. The transcriptional enhancer (ASE) responsible for the asymmetric expression of Nodal, Lefty, and Pitx2 was conserved among reptiles, including the Chinese softshell turtle and Madagascar ground gecko. Our findings suggest that Nodal, Lefty, and Pitx2 have the potential to be asymmetrically expressed in the developing brain of vertebrates, but that their expression varies even among reptiles.

Development of male-larger sexual size dimorphism in a lizard: IGF1 peak long after sexual maturity overlaps with pronounced growth in males

Squamate reptiles have been considered to be indeterminate growers for a long time. However, recent studies demonstrate that bone prolongation is stopped in many lizards by the closure of bone growth plates. This shift in the paradigm of lizard growth has important consequences for questions concerning the proximate causes of sexual size dimorphism. The traditional model of highly plastic and indeterminate growth would correspond more to a long-term action of a sex-specific growth regulator. On the other hand, determinate growth would be more consistent with a regulator acting in a sex-specific manner on the activity of bone growth plates operating during the phase when a dimorphism in size develops. We followed the growth of males and females of the male-larger Madagascar ground gecko (Paroedura picta) and monitored the activity of bone growth plates, gonad size, levels of steroids, expression of their receptors (AR, ESR1), and expression of genes from the insulin-like growth factor network (IGF1, IGF2, IGF1R, and IGF2R) in livers. Specifically, we measured gene expression before the onset of dimorphic growth, at the time when males have more active bone growth plates and sexual size dimorphism was clearly visible, and after a period of pronounced growth in both sexes. We found a significant spike in the expression of IGF1 in males around the time when dimorphism develops. This overexpression in males comes long after an increase in circulating testosterone levels and sexual maturation in males, and it might be suppressed by ovarian hormones in females. The results suggest that sexual size dimorphism in male-larger lizards can be caused by a positive effect of high levels of IGF1 on bone growth. The peak in IGF1 resembles the situation during the pubertal growth spurt in humans, but in lizards, it seems to be sex-specific and disconnected from sexual maturation.

Premeiotic endoreplication is essential for obligate parthenogenesis in geckos

Obligate parthenogenesis evolved in reptiles convergently several times, mainly through interspecific hybridization. The obligate parthenogenetic complexes typically include both diploid and triploid lineages. Offspring of parthenogenetic hybrids are genetic copies of their mother; however, the cellular mechanism enabling the production of unreduced cells is largely unknown. Here, we show that oocytes go through meiosis in three widespread, or even strongly invasive, obligate parthenogenetic complexes of geckos, namely in diploid and triploid Lepidodactylus lugubris, and triploid Hemiphyllodactylus typus and Heteronotia binoei. In all four lineages, the majority of oocytes enter the pachytene at the original ploidy level, but their chromosomes cannot pair properly and instead form univalents, bivalents and multivalents. Unreduced eggs with clonally inherited genomes are formed from germ cells that had undergone premeiotic endoreplication, in which appropriate segregation is ensured by the formation of bivalents made from copies of identical chromosomes. We conclude that the induction of premeiotic endoreplication in reptiles was independently co-opted at least four times as an essential component of parthenogenetic reproduction and that this mechanism enables the emergence of fertile polyploid lineages within parthenogenetic complexes.

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.

Remnant of Unrelated Amniote Sex Chromosomal Linkage Sharing on the Same Chromosome in House Gecko Lizards, Providing a Better Understanding of the Ancestral Super-Sex Chromosome

Comparative chromosome maps investigating sex chromosomal linkage groups in amniotes and microsatellite repeat motifs of a male house gecko lizard (Hemidactylus frenatus, HFR) and a flat-tailed house gecko lizard (H. platyurus, HPL) of unknown sex were examined using 75 bacterial artificial chromosomes (BACs) from chicken and zebra finch genomes. No massive accumulations of microsatellite repeat motifs were found in either of the gecko lizards, but 10 out of 13 BACs mapped on HPL chromosomes were associated with other amniote sex chromosomes. Hybridization of the same BACs onto multiple different chromosome pairs suggested transitions to sex chromosomes across amniotes. No BAC hybridization signals were found on HFR chromosomes. However, HFR diverged from HPL about 30 million years ago, possibly due to intrachromosomal rearrangements occurring in the HFR lineage. By contrast, heterochromatin likely reshuffled patterns between HPL and HFR, as observed from C-positive heterochromatin distribution. Six out of ten BACs showed partial homology with squamate reptile chromosome 2 (SR2) and snake Z and/or W sex chromosomes. The gecko lizard showed shared unrelated sex chromosomal linkages-the remnants of a super-sex chromosome. A large ancestral super-sex chromosome showed a correlation between SR2 and snake W sex chromosomes.

Evolutionary adaptation of visual pigments in geckos for their photic environment

Vertebrates generally have a single type of rod for scotopic vision and multiple types of cones for photopic vision. Noteworthily, nocturnal geckos transmuted ancestral photoreceptor cells into rods containing not rhodopsin but cone pigments, and, subsequently, diurnal geckos retransmuted these rods into cones containing cone pigments. High sensitivity of scotopic vision is underlain by the rod’s low background noise, which originated from a much lower spontaneous activation rate of rhodopsin than of cone pigments. Here, we revealed that nocturnal gecko cone pigments decreased their spontaneous activation rates to mimic rhodopsin, whereas diurnal gecko cone pigments recovered high rates similar to those of typical cone pigments. We also identified amino acid residues responsible for the alterations of the spontaneous activation rates. Therefore, we concluded that the switch between diurnality and nocturnality in geckos required not only morphological transmutation of photoreceptors but also adjustment of the spontaneous activation rates of visual pigments.

Technical considerations in Hi-C scaffolding and evaluation of chromosome-scale genome assemblies

The recent development of ecological studies has been fueled by the introduction of massive information based on chromosome‐scale genome sequences, even for species for which genetic linkage is not accessible. This was enabled mainly by the application of Hi‐C, a method for genome‐wide chromosome conformation capture that was originally developed for investigating the long‐range interaction of chromatins. Performing genomic scaffolding using Hi‐C data is highly resource‐demanding and employs elaborate laboratory steps for sample preparation. It starts with building a primary genome sequence assembly as an input, which is followed by computation for genome scaffolding using Hi‐C data, requiring careful validation. This article presents technical considerations for obtaining optimal Hi‐C scaffolding results and provides a test case of its application to a reptile species, the Madagascar ground gecko (Paroedura picta). Among the metrics that are frequently used for evaluating scaffolding results, we investigate the validity of the completeness assessment of chromosome‐scale genome assemblies using single‐copy reference orthologues.

Independent pseudogenizations and losses of sox15 during amniote diversification following asymmetric ohnolog evolution

Background

Four ohnologous genes (sox1, sox2, sox3, and sox15) were generated by two rounds of whole-genome duplication in a vertebrate ancestor. In eutherian mammals, Sox1, Sox2, and Sox3 participate in central nervous system (CNS) development. Sox15 has a function in skeletal muscle regeneration and has little functional overlap with the other three ohnologs. In contrast, the frog Xenopus laevis and zebrafish orthologs of sox15 as well as sox1-3 function in CNS development. We previously reported that Sox15 is involved in mouse placental development as neofunctionalization, but is pseudogenized in the marsupial opossum. These findings suggest that sox15 might have evolved with divergent gene fates during vertebrate evolution. However, knowledge concerning sox15 in other vertebrate lineages than therian mammals, anuran amphibians, and teleost fish is scarce. Our purpose in this study was to clarify the fate and molecular evolution of sox15 during vertebrate evolution.

Results

We searched for sox15 orthologs in all vertebrate classes from agnathans to mammals by significant sequence similarity and synteny analyses using vertebrate genome databases. Interestingly, sox15 was independently pseudogenized at least twice during diversification of the marsupial mammals. Moreover, we observed independent gene loss of sox15 at least twice during reptile evolution in squamates and crocodile-bird diversification. Codon-based phylogenetic tree and selective analyses revealed an increased d_N/d_S ratio for sox15 compared to the other three ohnologs during jawed vertebrate evolution.

Conclusions

The findings revealed an asymmetric evolution of sox15 among the four ohnologs during vertebrate evolution, which was supported by the increased d_N/d_S values in cartilaginous fishes, anuran amphibians, and amniotes. The increased d_N/d_S value of sox15 may have been caused mainly by relaxed selection. Notably, independent pseudogenizations and losses of sox15 were observed during marsupial and reptile evolution, respectively. Both might have been caused by strong relaxed selection. The drastic gene fates of sox15, including neofunctionalization and pseudogenizations/losses during amniote diversification, might be caused by a release from evolutionary constraints.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01864-z.

Shark and ray genomics for disentangling their morphological diversity and vertebrate evolution

Developmental studies of sharks and rays (elasmobranchs) have provided much insight into the process of morphological evolution of vertebrates. Although those studies are supposedly fueled by large-scale molecular sequencing information, whole-genome sequences of sharks and rays were made available only recently. One compelling difficulty of elasmobranch developmental biology is the low accessibility to embryonic study materials and their slow development. Another limiting factor is the relatively large size of their genomes. Moreover, their large body sizes restrict sustainable captive breeding, while their high body fluid osmolarity prevents reproducible cell culturing for in vitro experimentation, which has also limited our knowledge of their chromosomal organization for validation of genome sequencing products. This article focuses on egg-laying elasmobranch species used in developmental biology and provides an overview of the characteristics of the shark and ray genomes revealed to date. Developmental studies performed on a gene-by-gene basis are also reviewed from a whole-genome perspective. Among the popular regulatory genes studied in developmental biology, I scrutinize shark homologs of Wnt genes that highlight vanishing repertoires in many other vertebrate lineages, as well as Hox genes that underwent an unexpected modification unique to the elasmobranch lineage. These topics are discussed together with insights into the reconstruction of developmental programs in the common ancestor of vertebrates and its subsequent evolutionary trajectories that mark the features that are unique to, and those characterizing the diversity among, cartilaginous fishes.

Ten Years of Collaborative Progress in the Quest for Orthologs

Accurate determination of the evolutionary relationships between genes is a foundational challenge in biology. Homology-evolutionary relatedness-is in many cases readily determined based on sequence similarity analysis. By contrast, whether or not two genes directly descended from a common ancestor by a speciation event (orthologs) or duplication event (paralogs) is more challenging, yet provides critical information on the history of a gene. Since 2009, this task has been the focus of the Quest for Orthologs (QFO) Consortium. The sixth QFO meeting took place in Okazaki, Japan in conjunction with the 67th National Institute for Basic Biology conference. Here, we report recent advances, applications, and oncoming challenges that were discussed during the conference. Steady progress has been made toward standardization and scalability of new and existing tools. A feature of the conference was the presentation of a panel of accessible tools for phylogenetic profiling and several developments to bring orthology beyond the gene unit-from domains to networks. This meeting brought into light several challenges to come: leveraging orthology computations to get the most of the incoming avalanche of genomic data, integrating orthology from domain to biological network levels, building better gene models, and adapting orthology approaches to the broad evolutionary and genomic diversity recognized in different forms of life and viruses.

New Ther1-derived SINE Squam3 in scaled reptiles

BACKGROUND

SINEs comprise a significant part of animal genomes and are used to study the evolution of diverse taxa. Despite significant advances in SINE studies in vertebrates and higher eukaryotes in general, their own evolution is poorly understood.

RESULTS

We have discovered and described in detail a new Squam3 SINE specific for scaled reptiles (Squamata). The subfamilies of this SINE demonstrate different distribution in the genomes of squamates, which together with the data on similar SINEs in the tuatara allowed us to propose a scenario of their evolution in the context of reptilian evolution.

CONCLUSIONS

Ancestral SINEs preserved in small numbers in most genomes can give rise to taxa-specific SINE families. Analysis of this aspect of SINEs can shed light on the history and mechanisms of SINE variation in reptilian genomes.

Individual and age-related variation of cellular brain composition in a squamate reptile

Within-species variation in the number of neurons, other brain cells and their allocation to different brain parts is poorly studied. Here, we assess these numbers in a squamate reptile, the Madagascar ground gecko (Paroedura picta). We examined adults from two captive populations and three age groups within one population. Even though reptiles exhibit extensive adult neurogenesis, intrapopulation variation in the number of neurons is similar to that in mice. However, the two populations differed significantly in most measures, highlighting the fact that using only one population can underestimate within-species variation. There is a substantial increase in the number of neurons and decrease in neuronal density in adult geckos relative to hatchlings and an increase in the number of neurons in the telencephalon in fully grown adults relative to sexually mature young adults. This finding implies that adult neurogenesis does not only replace worn out but also adds new telencephalic neurons in reptiles during adulthood. This markedly contrasts with the situation in mammals, where the number of cortical neurons declines with age.

Immunolocalization of corneous proteins including a serine-tyrosine-rich beta-protein in the adhesive pads in the tokay gecko

Adhesive pads of geckos contain many thousands of nanoscale spatulae for the adhesion and movement along vertical or inverted surfaces. Setae are composed of interlaced corneous bundles made of small cysteine-glycine-rich corneous beta proteins (CBPs, formerly indicated as beta-keratins), embedded in a matrix material composed of cytoskeletal proteins and lipids. Negatively charged intermediate filament keratins (IFKs) and positively charged CBPs likely interact within setae, aside disulphide bonds, giving rise to a flexible and resistant corneous material. Using differernt antibodies against CBPs and IFKs an updated model of the composition of setae and spatulae is presented. Immunofluorescence and ultrastructural immunogold labeling reveal that one type of neutral serine-tyrosine-rich CBP is weakly localized in the setae while it is absent from the spatula. This uncharged protein is mainly present in the thin Oberhautchen layer sustaining the setae, although with a much lower intensity with respect to the cysteine-rich CBPs. These proteins in the spatula likely originate a positively charged or neutral contact surface with the substrate but the influence of lipids and cytoskeletal proteins present in setae on the mechanism of adhesion is not known. In the spatula, protein-lipid complexes may impart the pliability for the attachment and adapt to irregular surfaces. The presence of cysteine-glycine medium rich CBPs and softer IFKs in alpha-layers sustaining the setae forms a flexible base for compliance of the setae to substrate and improved adhesion.

Nodal paralogues underlie distinct mechanisms for visceral left-right asymmetry in reptiles and mammals

Unidirectional fluid flow generated by motile cilia at the left-right organizer (LRO) breaks left-right (L-R) symmetry during early embryogenesis in mouse, frog and zebrafish. The chick embryo, however, does not require motile cilia for L-R symmetry breaking. The diversity of mechanisms for L-R symmetry breaking among vertebrates and the trigger for such symmetry breaking in non-mammalian amniotes have remained unknown. Here we examined how L-R asymmetry is established in two reptiles, Madagascar ground gecko and Chinese softshell turtle. Both of these reptiles appear to lack motile cilia at the LRO. The expression of the Nodal gene at the LRO in the reptilian embryos was found to be asymmetric, in contrast to that in vertebrates such as mouse that are dependent on cilia for L-R patterning. Two paralogues of the Nodal gene derived from an ancient gene duplication are retained and expressed differentially in cilia-dependent and cilia-independent vertebrates. The expression of these two Nodal paralogues is similarly controlled in the lateral plate mesoderm but regulated differently at the LRO. Our in-depth analysis of reptilian embryos thus suggests that mammals and non-mammalian amniotes deploy distinct strategies dependent on different Nodal paralogues for rendering Nodal activity asymmetric at the LRO.

Comparative and evolutionary analysis of the reptilian hedgehog gene family (Shh, Dhh, and Ihh)

The hedgehog signaling pathway plays a vital role in human and animal patterning and cell proliferation during the developmental process. The hedgehog gene family of vertebrate species includes three genes, Shh, Dhh, and Ihh, which possess different functions and expression patterns. Despite the importance of hedgehog genes, genomic evidence of this gene family in reptiles is lacking. In this study, the available genomes of a number of representative reptile species were explored by utilizing adaptive evolutionary analysis methods to characterize the evolutionary patterns of the hedgehog gene family. Altogether, 33 sonic hedgehog (Shh), 25 desert hedgehog (Dhh), and 20 Indian hedgehog (Ihh) genes were obtained from reptiles, and six avian and five mammalian sequences were added to the analysis. The phylogenetic maximum likelihood (ML) tree of the Shh, Dhh, and Ihh genes revealed a similar topology, which is approximately consistent with the traditional taxonomic group. No shared positive selection site was identified by the PAML site model or the three methods in the Data Monkey Server. Branch model and Clade model C analyses revealed that the Dhh and Ihh genes experienced different evolutionary forces in reptiles and other vertebrates, while the Shh gene was not significantly different in terms of selection pressure. The different evolutionary rates of the Dhh and Ihh genes suggest that these genes may be potential contributors to the discrepant sperm and body development of different clades. The different adaptive evolutionary history of the Shh, Dhh, and Ihh genes among reptiles may be due to their different functions in regulating cellular events of development from the embryonic stages to adulthood. Overall, this study has provided meaningful information regarding the evolution of the hedgehog gene family in reptiles and a theoretical foundation for further analyses on the functional and molecular mechanisms that have shaped the reptilian hedgehog genes.

Discovery of All Three Types in Cartilaginous Fishes Enables Phylogenetic Resolution of the Origins and Evolution of Interferons

Interferons orchestrate host antiviral responses in jawed vertebrates. They are categorized into three classes; IFN1 and IFN3 are the primary antiviral cytokine lineages, while IFN2 responds to a broader variety of pathogens. The evolutionary relationships within and between these three classes have proven difficult to resolve. Here, we reassess interferon evolution, considering key phylogenetic pitfalls including taxon sampling, alignment quality, model adequacy, and outgroup choice. We reveal that cartilaginous fishes, and hence the jawed vertebrate ancestor, possess(ed) orthologs of all three interferon classes. We show that IFN3 groups sister to IFN1, resolve the origins of the human IFN3 lineages, and find that intronless IFN3s emerged at least three times. IFN2 genes are highly conserved, except for IFN-γ-rel, which we confirm resulted from a teleost-specific duplication. Our analyses show that IFN1 phylogeny is highly sensitive to phylogenetic error. By accounting for this, we describe a new backbone IFN1 phylogeny that implies several IFN1 genes existed in the jawed vertebrate ancestor. One of these is represented by the intronless IFN1s of tetrapods, including mammalian-like repertoires of reptile IFN1s and a subset of amphibian IFN1s, in addition to newly-identified intron-containing shark IFN1 genes. IFN-f, previously only found in teleosts, likely represents another ancestral jawed vertebrate IFN1 family member, suggesting the current classification of fish IFN1s into two groups based on the number of cysteines may need revision. The providence of the remaining fish IFN1s and the coelacanth IFN1s proved difficult to resolve, but they may also be ancestral jawed vertebrate IFN1 lineages. Finally, a large group of amphibian-specific IFN1s falls sister to all other IFN1s and was likely also present in the jawed vertebrate ancestor. Our results verify that intronless IFN1s have evolved multiple times in amphibians and indicate that no one-to-one orthology exists between mammal and reptile IFN1s. Our data also imply that diversification of the multiple IFN1s present in the jawed vertebrate ancestor has occurred through a rapid birth-death process, consistent with functional maintenance over a 450-million-year host-pathogen arms race. In summary, this study reveals a new model of interferon evolution important to our understanding of jawed vertebrate antiviral immunity.

Evolution of nodal and nodal-related genes and the putative composition of the heterodimers that trigger the nodal pathway in vertebrates

Nodal is a signaling molecule that belongs to the transforming growth factor-β superfamily that plays key roles during the early stages of development of animals. In vertebrates Nodal forms an heterodimer with a GDF1/3 protein to activate the Nodal pathway. Vertebrates have a paralog of nodal in their genomes labeled Nodal-related, but the evolutionary history of these genes is a matter of debate, mainly because of the presence of a variable numbers of genes in the vertebrate genomes sequenced so far. Thus, the goal of this study was to investigate the evolutionary history of the Nodal and Nodal-related genes with an emphasis in tracking changes in the number of genes among vertebrates. Our results show the presence of two gene lineages (Nodal and Nodal-related) that can be traced back to the ancestor of jawed vertebrates. These lineages have undergone processes of differential retention and lineage-specific expansions. Our results imply that Nodal and Nodal-related duplicated at the latest in the ancestor of gnathostomes, and they still retain a significant level of functional redundancy. By comparing the evolution of the Nodal/Nodal-related with GDF1/3 gene family, it is possible to infer that there are several types of heterodimers that can trigger the Nodal pathway among vertebrates.

Duplications in Corneous Beta Protein Genes and the Evolution of Gecko Adhesion

Corneous proteins are an important component of the tetrapod integument. Duplication and diversification of keratins and associated proteins are linked with the origin of most novel integumentary structures like mammalian hair, avian feathers, and scutes covering turtle shells. Accordingly, the loss of integumentary structures often coincides with the loss of genes encoding keratin and associated proteins. For example, many hair keratins in dolphins and whales have become pseudogenes. The adhesive setae of geckos and anoles are composed of both intermediate filament keratins (IF-keratins, formerly known as alpha-keratins) and corneous beta-proteins (CBPs, formerly known as beta-keratins) and recent whole genome assemblies of two gecko species and an anole uncovered duplications in seta-specific CBPs in each of these lineages. While anoles evolved adhesive toepads just once, there are two competing hypotheses about the origin(s) of digital adhesion in geckos involving either a single origin or multiple origins. Using data from three published gecko genomes, I examine CBP gene evolution in geckos and find support for a hypothesis where CBP gene duplications are associated with the repeated evolution of digital adhesion. Although these results are preliminary, I discuss how additional gecko genome assemblies, combined with phylogenies of keratin and associated protein genes and gene duplication models, can provide rigorous tests of several hypotheses related to gecko CBP evolution. This includes a taxon sampling strategy for sequencing and assembly of gecko genomes that could help resolve competing hypotheses surrounding the origin(s) of digital adhesion.

Asymmetric paralog evolution between the "cryptic" gene Bmp16 and its well-studied sister genes Bmp2 and Bmp4

The vertebrate gene repertoire is characterized by “cryptic” genes whose identification has been hampered by their absence from the genomes of well-studied species. One example is the Bmp16 gene, a paralog of the developmental key genes Bmp2 and -4. We focus on the Bmp2/4/16 group of genes to study the evolutionary dynamics following gen(om)e duplications with special emphasis on the poorly studied Bmp16 gene. We reveal the presence of Bmp16 in chondrichthyans in addition to previously reported teleost fishes and reptiles. Using comprehensive, vertebrate-wide gene sampling, our phylogenetic analysis complemented with synteny analyses suggests that Bmp2, -4 and -16 are remnants of a gene quartet that originated during the two rounds of whole-genome duplication (2R-WGD) early in vertebrate evolution. We confirm that Bmp16 genes were lost independently in at least three lineages (mammals, archelosaurs and amphibians) and report that they have elevated rates of sequence evolution. This finding agrees with their more “flexible” deployment during development; while Bmp16 has limited embryonic expression domains in the cloudy catshark, it is broadly expressed in the green anole lizard. Our study illustrates the dynamics of gene family evolution by integrating insights from sequence diversification, gene repertoire changes, and shuffling of expression domains.

Evaluating Genome Assemblies and Gene Models Using gVolante

In daily practice of de novo genome assembly and gene prediction, it would be a natural urge to evaluate their products. Different programs and parameter settings give rise to variable outputs, which leaves a decision of which output to adopt for downstream analysis for addressing biological questions. Instead of superficial assessment of length-based statistics of output sequences (e.g., N50 scaffold length), completeness assessment by means of scoring the coverage of reference orthologs has been increasingly utilized.We previously launched a web service, gVolante ( https://gvolante.riken.jp /), to provide a user-friendly interface and a uniform environment for completeness assessment with the pipelines CEGMA and BUSCO. Completeness assessments performed on gVolante report scores based on not just the coverage of reference genes but also on sequence lengths, allowing quality control in multiple aspects. This chapter focuses on the procedure for such assessment and provides technical tips for higher accuracy.

Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates

Modern cartilaginous fishes are divided into elasmobranchs (sharks, rays and skates) and chimaeras, and the lack of established whole-genome sequences for the former has prevented our understanding of early vertebrate evolution and the unique phenotypes of elasmobranchs. Here we present de novo whole-genome assemblies of brownbanded bamboo shark and cloudy catshark and an improved assembly of the whale shark genome. These relatively large genomes (3.8-6.7 Gbp) contain sparse distributions of coding genes and regulatory elements and exhibit reduced molecular evolutionary rates. Our thorough genome annotation revealed Hox C genes previously hypothesized to have been lost, as well as distinct gene repertories of opsins and olfactory receptors that would be associated with adaptation to unique underwater niches. We also show the early establishment of the genetic machinery governing mammalian homoeostasis and reproduction at the jawed vertebrate ancestor. This study, supported by genomic, transcriptomic and epigenomic resources, provides a foundation for the comprehensive, molecular exploration of phenotypes unique to sharks and insights into the evolutionary origins of vertebrates.