The transcriptome of the veiled chameleon (Chamaeleo calyptratus): A resource for studying the evolution and development of vertebrates

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.

A lizard is never late: Squamate genomics as a recent catalyst for understanding sex chromosome and microchromosome evolution

In 2011, the first high-quality genome assembly of a squamate reptile (lizard or snake) was published for the green anole. Dozens of genome assemblies were subsequently published over the next decade, yet these assemblies were largely inadequate for answering fundamental questions regarding genome evolution in squamates due to their lack of contiguity or annotation. As the "genomics age" was beginning to hit its stride in many organismal study systems, progress in squamates was largely stagnant following the publication of the green anole genome. In fact, zero high-quality (chromosome-level) squamate genomes were published between the years 2012 and 2017. However, since 2018, an exponential increase in high-quality genome assemblies has materialized with 24 additional high-quality genomes published for species across the squamate tree of life. As the field of squamate genomics is rapidly evolving, we provide a systematic review from an evolutionary genomics perspective. We collated a near-complete list of publicly available squamate genome assemblies from more than half-a-dozen international and third-party repositories and systematically evaluated them with regard to their overall quality, phylogenetic breadth, and usefulness for continuing to provide accurate and efficient insights into genome evolution across squamate reptiles. This review both highlights and catalogs the currently available genomic resources in squamates and their ability to address broader questions in vertebrates, specifically sex chromosome and microchromosome evolution, while addressing why squamates may have received less historical focus and has caused their progress in genomics to lag behind peer taxa.

A lizard is never late: squamate genomics as a recent catalyst for understanding sex chromosome and microchromosome evolution

In 2011, the first high-quality genome assembly of a squamate reptile (lizard or snake) was published for the green anole. Dozens of genome assemblies were subsequently published over the next decade, yet these assemblies were largely inadequate for answering fundamental questions regarding genome evolution in squamates due to their lack of contiguity or annotation. As the "genomics age" was beginning to hit its stride in many organismal study systems, progress in squamates was largely stagnant following the publication of the green anole genome. In fact, zero high-quality (chromosome-level) squamate genomes were published between the years 2012-2017. However, since 2018, an exponential increase in high-quality genome assemblies has materialized with 24 additional high-quality genomes published for species across the squamate tree of life. As the field of squamate genomics is rapidly evolving, we provide a systematic review from an evolutionary genomics perspective. We collated a near-complete list of publicly available squamate genome assemblies from more than half-a-dozen international and third-party repositories and systematically evaluated them with regard to their overall quality, phylogenetic breadth, and usefulness for continuing to provide accurate and efficient insights into genome evolution across squamate reptiles. This review both highlights and catalogs the currently available genomic resources in squamates and their ability to address broader questions in vertebrates, specifically sex chromosome and microchromosome evolution, while addressing why squamates may have received less historical focus and has caused their progress in genomics to lag behind peer taxa.

Morphological changes and two Nodal paralogs drive left-right asymmetry in the squamate veiled chameleon (C. calyptratus)

The ancestral mode of left-right (L-R) patterning involves cilia in the L-R organizer. However, the mechanisms regulating L-R patterning in non-avian reptiles remains an enigma, since most squamate embryos are undergoing organogenesis at oviposition. In contrast, veiled chameleon (Chamaeleo calyptratus) embryos are pre-gastrula at oviposition, making them an excellent organism for studying L-R patterning evolution. Here we show that veiled chameleon embryos lack motile cilia at the time of L-R asymmetry establishment. Thus, the loss of motile cilia in the L-R organizers is a synapomorphy of all reptiles. Furthermore, in contrast to avians, geckos and turtles, which have one Nodal gene, veiled chameleon exhibits expression of two paralogs of Nodal in the left lateral plate mesoderm, albeit in non-identical patterns. Using live imaging, we observed asymmetric morphological changes that precede, and likely trigger, asymmetric expression of the Nodal cascade. Thus, veiled chameleons are a new and unique model for studying the evolution of L-R patterning.

Identification of Iguania Ancestral Syntenic Blocks and Putative Sex Chromosomes in the Veiled Chameleon (Chamaeleo calyptratus, Chamaeleonidae, Iguania)

The veiled chameleon (Chamaeleo calyptratus) is a typical member of the family Chamaeleonidae and a promising object for comparative cytogenetics and genomics. The karyotype of C. calyptratus differs from the putative ancestral chameleon karyotype (2n = 36) due to a smaller chromosome number (2n = 24) resulting from multiple chromosome fusions. The homomorphic sex chromosomes of an XX/XY system were described recently using male-specific RADseq markers. However, the chromosomal pair carrying these markers was not identified. Here we obtained chromosome-specific DNA libraries of C. calyptratus by chromosome flow sorting that were assigned by FISH and sequenced. Sequence comparison with three squamate reptiles reference genomes revealed the ancestral syntenic regions in the C. calyptratus chromosomes. We demonstrated that reducing the chromosome number in the C. calyptratus karyotype occurred through two fusions between microchromosomes and four fusions between micro-and macrochromosomes. PCR-assisted mapping of a previously described Y-specific marker indicates that chromosome 5 may be the sex chromosome pair. One of the chromosome 5 conserved synteny blocks shares homology with the ancestral pleurodont X chromosome, assuming parallelism in the evolution of sex chromosomes from two basal Iguania clades (pleurodonts and acrodonts). The comparative chromosome map produced here can serve as the foundation for future genome assembly of chameleons and vertebrate-wide comparative genomic studies.

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.

Two novel bornaviruses identified in colubrid and viperid snakes

We present the complete genome sequences of Caribbean watersnake bornavirus (CWBV) and Mexican black-tailed rattlesnake bornavirus (MRBV), which we identified in archived raw transcriptomic read data of a Caribbean watersnake (Tretanorhinus variabilis) and a Mexican black-tailed rattlesnake (Crotalus molossus nigrescens), respectively. The genomes of CWBV and MRBV have a length of about 8,900 nucleotides and comprise the complete coding regions and the untranslated regions. The overall genomic makeup and predicted gene content is typical for members of the genus Orthobornavirus within the family Bornaviridae. Alternative splicing was detected for the L and M genes. Based on a phylogenetic analysis of all viral proteins, we consider both viruses to be members of a single novel species within the genus Orthobornavirus. Both viruses form a distinct outgroup to all currently known orthobornaviruses. Based on the novel virus genomes, we furthermore identified closely related endogenous bornavirus-like nucleoprotein sequences in transcriptomic data of veiled chameleons (Chamaeleo calyptratus) and a common lancehead (Bothrops atrox).

A chromosome-level genome assembly of the parasitoid wasp, Cotesia glomerata (Hymenoptera: Braconidae)

Hymenopterans make up about 20% of all animal species, but most are poorly known and lack high-quality genomic resources. One group of important, yet under-studied hymenopterans, are parasitoid wasps in the family Braconidae. Among this under-studied group are braconid wasps in the genus Cotesia; a clade of ~1,000 species routinely used in studies of physiology, ecology, biological control, and genetics. However, our ability to understand these organisms has been hindered by a lack of genomic resources. We helped bridge this gap by generating a high-quality genome assembly for the parasitoid wasp, Cotesia glomerata (Braconidae; Microgastrinae). We generated this assembly using multiple sequencing technologies, including Oxford Nanopore, whole-genome shotgun sequencing, and 3-D chromatin contact information (Hi-C). Our assembly is one of the most contiguous, complete, and publicly available hymenopteran genomes, represented by 3,355 scaffolds with a scaffold N50 of ~28Mb and a BUSCO score of ~99%. Given the genome sizes found in closely related species, our genome assembly was ~50% larger than expected, which was apparently induced by runaway amplification of three types of repetitive elements: simple repeats, Long Terminal Repeats (LTRs), and Long Interspersed Nuclear Elements (LINEs). This assembly is another step forward for genomics across this hyper-diverse, yet understudied, order of insects. The assembled genomic data and metadata files are publicly available via Figshare (https://doi.org/10.6084/m9.figshare.13010549).

The development, patterning and evolution of neural crest cell differentiation into cartilage and bone

Neural crest cells are a vertebrate-specific migratory, multipotent cell population that give rise to a diverse array of cells and tissues during development. Cranial neural crest cells, in particular, generate cartilage, bone, tendons and connective tissue in the head and face as well as neurons, glia and melanocytes. In this review, we focus on the chondrogenic and osteogenic potential of cranial neural crest cells and discuss the roles of Sox9, Runx2 and Msx1/2 transcription factors and WNT, FGF and TGFβ signaling pathways in regulating neural crest cell differentiation into cartilage and bone. We also describe cranioskeletal defects and disorders arising from gain or loss-of-function of genes that are required for patterning and differentiation of cranial neural crest cells. Finally, we discuss the evolution of skeletogenic potential in neural crest cells and their function as a conduit for intraspecies and interspecies variation, and the evolution of craniofacial novelties.

Polarized Sonic Hedgehog Protein Localization and a Shift in the Expression of Region-Specific Molecules Is Associated With the Secondary Palate Development in the Veiled Chameleon

Secondary palate development is characterized by the formation of two palatal shelves on the maxillary prominences, which fuse in the midline in mammalian embryos. However, in reptilian species, such as turtles, crocodilians, and lizards, the palatal shelves of the secondary palate develop to a variable extent and morphology. While in most Squamates, the palate is widely open, crocodilians develop a fully closed secondary palate. Here, we analyzed developmental processes that underlie secondary palate formation in chameleons, where large palatal shelves extend horizontally toward the midline. The growth of the palatal shelves continued during post-hatching stages and closure of the secondary palate can be observed in several adult animals. The massive proliferation of a multilayered oral epithelium and mesenchymal cells in the dorsal part of the palatal shelves underlined the initiation of their horizontal outgrowth, and was decreased later in development. The polarized cellular localization of primary cilia and Sonic hedgehog protein was associated with horizontal growth of the palatal shelves. Moreover, the development of large palatal shelves, supported by the pterygoid and palatine bones, was coupled with the shift in Meox2, Msx1, and Pax9 gene expression along the rostro-caudal axis. In conclusion, our results revealed distinctive developmental processes that contribute to the expansion and closure of the secondary palate in chameleons and highlighted divergences in palate formation across amniote species.

Transcriptomic data support a nocturnal bottleneck in the ancestor of gecko lizards

Gecko lizards are a species-rich clade of primarily-nocturnal squamate reptiles. In geckos, adaptations to nocturnality have dramatically reshaped the eye. Perhaps the most notable change is the loss of rod cells in the retina and subsequent "transmutation" of cones into a rod-like morphology and physiology. While many studies have noted the absence of some rod-specific genes, such as the visual pigment Rhodopsin (RH1), these studies have focused on just a handful of species that are nested deep in the gecko phylogeny. Thus, it is not clear whether these changes arose through convergence, are homologous and ubiquitous across geckos, or restricted to a subset of species. Here, we used de novo eye transcriptomes from five gecko species, and genomes from two additional gecko species, representing the breadth of extant gecko diversity (i.e. 4 of the 7 gecko families, spanning the deepest divergence of crown Gekkota), to show that geckos lost expression of almost the entire suite of necessary rod-cell phototransduction genes in the eye, distinct from all other squamate reptiles. Geckos are the first vertebrate group to have lost their complete rod-cell expression pathway, not just the visual pigment. In addition, all sampled species have also lost expression of the cone-opsin SWS2 visual pigment. These results strongly suggest a single loss of rod cells and subsequent cone-to-rod transmutation that occurred prior to the diversification of extant geckos.

Primary myogenesis in the sand lizard (Lacerta agilis) limb bud

Our studies conducted on reptilian limb muscle development revealed, for the first time, early forelimb muscle differentiation at the morphological and molecular level. Sand lizard skeletal muscle differentiation in the early forelimb bud was investigated by light, confocal, and transmission electron microscopy as well as western blot. The early forelimb bud, filled with mesenchymal cells, is surrounded by monolayer epithelium cells. The immunocytochemical analysis revealed the presence of Pax3- and Lbx-positive cells in the vicinity of the ventro-lateral lip (VLL) of the dermomyotome, suggesting that VLL is the source of limb muscle progenitor cells. Furthermore, Pax3- and Lbx-positive cells were observed in the dorsal and ventral myogenic pools of the forelimb bud. Skeletal muscle development in the early limb bud is asynchronous, which is manifested by the presence of myogenic cells in different stages of differentiation: multinucleated myotubes with well-developed contractile apparatus, myoblasts, and mitotically active premyoblasts. The western blot analysis revealed the presence of MyoD and Myf5 proteins in all investigated developmental stages. The MyoD western blot analysis showed two bands corresponding to monomeric (mMyoD) and dimeric (dMyoD) fractions. Two separate bands were also detected in the case of Myf5. The observed bands were related to non-phosphorylated (Myf5) and phosphorylated (pMyf5) fractions of Myf5. Our investigations on sand lizard forelimb myogenesis showed that the pattern of muscle differentiation in the early forelimb bud shares many features with rodents and chicks.