The Hox code responsible for the patterning of the anterior vertebrae in zebrafish.

The vertebral column is a characteristic structure of vertebrates. Genetic studies in mice have shown that Hox-mediated patterning plays a key role in specifying discrete anatomical regions of the vertebral column. Expression pattern analyses in several vertebrate embryos have provided correlative evidence that the anterior boundaries of Hox expression coincide with distinct anatomical vertebrae. However, because functional analyses have been limited to mice, it remains unclear which Hox genes actually function in vertebral patterning in other vertebrates. In this study, various zebrafish Hox mutants were generated for loss-of-function phenotypic analysis to functionally decipher the Hox code responsible for the zebrafish anterior vertebrae between the occipital and thoracic vertebrae. We found that Hox genes in HoxB- and HoxC-related clusters participate in regulating the morphology of the zebrafish anterior vertebrae. In addition, medaka hoxc6a was found to be responsible for anterior vertebral identity, as in zebrafish. Based on phenotypic similarities with Hoxc6 knockout mice, our results suggest that the Hox patterning system, including at least Hoxc6, may have been functionally established in the vertebral patterning of the common ancestor of ray-finned and lobe-finned fishes.

Teleost Hox code defines regional identities competent for the formation of dorsal and anal fins.

The dorsal and anal fins can vary widely in position and length along the anterior-posterior axis in teleost fishes. However, the molecular mechanisms underlying the diversification of these fins remain unknown. Here, we used genetic approaches in zebrafish and medaka, in which the relative positions of the dorsal and anal fins are opposite, to demonstrate the crucial role of hox genes in the patterning of the teleost posterior body, including the dorsal and anal fins. By the CRISPR-Cas9-induced frameshift mutations and positional cloning of spontaneous dorsalfinless medaka, we show that various hox mutants exhibit the absence of dorsal or anal fins, or a stepwise posterior extension of these fins, with vertebral abnormalities. Our results indicate that multiple hox genes, primarily from hoxc-related clusters, encompass the regions responsible for the dorsal and anal fin formation along the anterior-posterior axis. These results further suggest that shifts in the anterior boundaries of hox expression which vary among fish species, lead to diversification in the position and size of the dorsal and anal fins, similar to how modulations in Hox expression can alter the number of anatomically distinct vertebrae in tetrapods. Furthermore, we show that hox genes responsible for dorsal fin formation are different between zebrafish and medaka. Our results suggest that a novel mechanism has occurred during teleost evolution, in which the gene network responsible for fin formation might have switched to the regulation downstream of other hox genes, leading to the remarkable diversity in the dorsal fin position.

Ancient vertebrate dermal armor evolved from trunk neural crest.

Bone is an evolutionary novelty of vertebrates, likely to have first emerged as part of ancestral dermal armor that consisted of osteogenic and odontogenic components. Whether these early vertebrate structures arose from mesoderm or neural crest cells has been a matter of considerable debate. To examine the developmental origin of the bony part of the dermal armor, we have performed in vivo lineage tracing in the sterlet sturgeon, a representative of nonteleost ray-finned fish that has retained an extensive postcranial dermal skeleton. The results definitively show that sterlet trunk neural crest cells give rise to osteoblasts of the scutes. Transcriptional profiling further reveals neural crest gene signature in sterlet scutes as well as bichir scales. Finally, histological and microCT analyses of ray-finned fish dermal armor show that their scales and scutes are formed by bone, dentin, and hypermineralized covering tissues, in various combinations, that resemble those of the first armored vertebrates. Taken together, our results support a primitive skeletogenic role for the neural crest along the entire body axis, that was later progressively restricted to the cranial region during vertebrate evolution. Thus, the neural crest was a crucial evolutionary innovation driving the origin and diversification of dermal armor along the entire body axis.

On the evolutionary origins and regionalization of the neural crest.

The neural crest is a vertebrate-specific embryonic stem cell population that gives rise to a vast array of cell types throughout the animal body plan. These cells are first born at the edges of the central nervous system, from which they migrate extensively and differentiate into multiple cellular derivatives. Given the unique set of structures these cells comprise, the origin of the neural crest is thought to have important implications for the evolution and diversification of the vertebrate clade. In jawed vertebrates, neural crest cells exist as distinct subpopulations along the anterior-posterior axis. These subpopulations differ in terms of their respective differentiation potential and cellular derivatives. Thus, the modern neural crest is characterized as multipotent, migratory, and regionally segregated throughout the embryo. Here, we retrace the evolutionary origins of the neural crest, from the appearance of conserved regulatory circuitry in basal chordates to the emergence of neural crest subpopulations in higher vertebrates. Finally, we discuss a stepwise trajectory by which these cells may have arisen and diversified throughout vertebrate evolution.

Making a head: Neural crest and ectodermal placodes in cranial sensory development.

During development of the vertebrate sensory system, many important components like the sense organs and cranial sensory ganglia arise within the head and neck. Two progenitor populations, the neural crest, and cranial ectodermal placodes, contribute to these developing vertebrate peripheral sensory structures. The interactions and contributions of these cell populations to the development of the lens, olfactory, otic, pituitary gland, and cranial ganglia are vital for appropriate peripheral nervous system development. Here, we review the origins of both neural crest and placode cells at the neural plate border of the early vertebrate embryo and investigate the molecular and environmental signals that influence specification of different sensory regions. Finally, we discuss the underlying molecular pathways contributing to the complex vertebrate sensory system from an evolutionary perspective, from basal vertebrates to amniotes.

Threespine Stickleback: A Model System For Evolutionary Genomics.

The repeated adaptation of oceanic threespine sticklebacks to fresh water has made it a premier organism to study parallel evolution. These small fish have multiple distinct ecotypes that display a wide range of diverse phenotypic traits. Ecotypes are easily crossed in the laboratory, and families are large and develop quickly enough for quantitative trait locus analyses, positioning the threespine stickleback as a versatile model organism to address a wide range of biological questions. Extensive genomic resources, including linkage maps, a high-quality reference genome, and developmental genetics tools have led to insights into the genomic basis of adaptation and the identification of genomic changes controlling traits in vertebrates. Recently, threespine sticklebacks have been used as a model system to identify the genomic basis of highly complex traits, such as behavior and host-microbiome and host-parasite interactions. We review the latest findings and new avenues of research that have led the threespine stickleback to be considered a supermodel of evolutionary genomics.

An atlas of seven zebrafish hox cluster mutants provides insights into sub/neofunctionalization of vertebrate Hox clusters.

Vertebrate Hox clusters are comprised of multiple Hox genes that control morphology and developmental timing along multiple body axes. Although results of genetic analyses using Hox-knockout mice have been accumulating, genetic studies in other vertebrates have not been sufficient for functional comparisons of vertebrate Hox genes. In this study, we isolated all of the seven hox cluster loss-of-function alleles in zebrafish using the CRISPR-Cas9 system. Comprehensive analysis of the embryonic phenotype and X-ray micro-computed tomography scan analysis of adult fish revealed several species-specific functional contributions of homologous Hox clusters along the appendicular axis, whereas important shared general principles were also confirmed, as exemplified by serial anterior vertebral transformations along the main body axis, observed in fish for the first time. Our results provide insights into discrete sub/neofunctionalization of vertebrate Hox clusters after quadruplication of the ancient Hox cluster. This set of seven complete hox cluster loss-of-function alleles provide a formidable resource for future developmental genetic analysis of the Hox patterning system in zebrafish.

100-My history of bornavirus infections hidden in vertebrate genomes.

Although viruses have threatened our ancestors for millions of years, prehistoric epidemics of viruses are largely unknown. Endogenous bornavirus-like elements (EBLs) are ancient bornavirus sequences derived from the viral messenger RNAs that were reverse transcribed and inserted into animal genomes, most likely by retrotransposons. These elements can be used as molecular fossil records to trace past bornaviral infections. In this study, we systematically identified EBLs in vertebrate genomes and revealed the history of bornavirus infections over nearly 100 My. We confirmed that ancient bornaviral infections have occurred in diverse vertebrate lineages, especially in primate ancestors. Phylogenetic analyses indicated that primate ancestors were infected with various bornaviral lineages during evolution. EBLs in primate genomes formed clades according to their integration ages, suggesting that bornavirus lineages infected with primate ancestors had changed chronologically. However, some bornaviral lineages may have coexisted with primate ancestors and underwent repeated endogenizations for tens of millions of years. Moreover, a bornaviral lineage that coexisted with primate ancestors also endogenized in the genomes of some ancestral bats. The habitats of these bat ancestors have been reported to overlap with the migration route of primate ancestors. These results suggest that long-term virus-host coexistence expanded the geographic distributions of the bornaviral lineage along with primate migration and may have spread their infections to these bat ancestors. Our findings provide insight into the history of bornavirus infections over geological timescales that cannot be deduced from research using extant viruses alone, thus broadening our perspective on virus-host coevolution.

A new role for joint mobility in reconstructing vertebrate locomotor evolution.

Reconstructions of movement in extinct animals are critical to our understanding of major transformations in vertebrate locomotor evolution. Estimates of joint range of motion (ROM) have long been used to exclude anatomically impossible joint poses from hypothesized gait cycles. Here we demonstrate how comparative ROM data can be harnessed in a different way to better constrain locomotor reconstructions. As a case study, we measured nearly 600,000 poses from the hindlimb joints of the Helmeted Guineafowl and American alligator, which represent an extant phylogenetic bracket for the archosaurian ancestor and its pseudosuchian (crocodilian line) and ornithodiran (bird line) descendants. We then used joint mobility mapping to search for a consistent relationship between full potential joint mobility and the subset of joint poses used during locomotion. We found that walking and running poses are predictably located within full mobility, revealing additional constraints for reconstructions of extinct archosaurs. The inferential framework that we develop here can be expanded to identify ROM-based constraints for other animals and, in turn, will help to unravel the history of vertebrate locomotor evolution.

Parallel Evolution of Ameloblastic scpp Genes in Bony and Cartilaginous Vertebrates.

In bony vertebrates, skeletal mineralization relies on the secretory calcium-binding phosphoproteins (Scpp) family whose members are acidic extracellular proteins posttranslationally regulated by the Fam20°C kinase. As scpp genes are absent from the elephant shark genome, they are currently thought to be specific to bony fishes (osteichthyans). Here, we report a scpp gene present in elasmobranchs (sharks and rays) that evolved from local tandem duplication of sparc-L 5' exons and show that both genes experienced recent gene conversion in sharks. The elasmobranch scpp is remarkably similar to the osteichthyan scpp members as they share syntenic and gene structure features, code for a conserved signal peptide, tyrosine-rich and aspartate/glutamate-rich regions, and harbor putative Fam20°C phosphorylation sites. In addition, the catshark scpp is coexpressed with sparc-L and fam20°C in tooth and scale ameloblasts, similarly to some osteichthyan scpp genes. Despite these strong similarities, molecular clock and phylogenetic data demonstrate that the elasmobranch scpp gene originated independently from the osteichthyan scpp gene family. Our study reveals convergent events at the sparc-L locus in the two sister clades of jawed vertebrates, leading to parallel diversification of the skeletal biomineralization toolkit. The molecular evolution of sparc-L and its coexpression with fam20°C in catshark ameloblasts provides a unifying genetic basis that suggests that all convergent scpp duplicates inherited similar features from their sparc-L precursor. This conclusion supports a single origin for the hypermineralized outer odontode layer as produced by an ancestral developmental process performed by Sparc-L, implying the homology of the enamel and enameloid tissues in all vertebrates.

Evolution of vertebrate gill covers via shifts in an ancient Pou3f3 enhancer.

Whereas the gill chambers of jawless vertebrates open directly into the environment, jawed vertebrates evolved skeletal appendages that drive oxygenated water unidirectionally over the gills. A major anatomical difference between the two jawed vertebrate lineages is the presence of a single large gill cover in bony fishes versus separate covers for each gill chamber in cartilaginous fishes. Here, we find that these divergent patterns correlate with the pharyngeal arch expression of Pou3f3 orthologs. We identify a deeply conserved Pou3f3 arch enhancer present in humans through sharks but undetectable in jawless fish. Minor differences between the bony and cartilaginous fish enhancers account for their restricted versus pan-arch expression patterns. In zebrafish, mutation of Pou3f3 or the conserved enhancer disrupts gill cover formation, whereas ectopic pan-arch Pou3f3b expression generates ectopic skeletal elements resembling the multimeric covers of cartilaginous fishes. Emergence of this Pou3f3 arch enhancer >430 Mya and subsequent modifications may thus have contributed to the acquisition and diversification of gill covers and respiratory strategies during gnathostome evolution.

Assessing Myf5 and Lbx1 contribution to carapace development by reproducing their turtle-specific signatures in mouse embryos.

BACKGROUND: The turtle carapace is an evolutionary novelty resulting from changes in the processes that build ribs and their associated muscles in most tetrapod species. Turtle embryos have several unique features that might play a role in this process, including the carapacial ridge, a Myf5 gene with shorter coding region that generates an alternative splice variant lacking exon 2, and unusual expression patterns of Lbx1 and HGF. RESULTS: We investigated these turtle-specific expression differences using genetic approaches in mouse embryos. At mid-gestation, mouse embryos producing Myf5 transcripts lacking exon 2 replicated some early properties of turtle somites, but still developed into viable and fertile mice. Extending Lbx1 expression into the hypaxial dermomyotomal lip of trunk somites to mimic the turtle Lbx1 expression pattern, produced fusions in the distal part of the ribs. CONCLUSIONS: Turtle-like Myf5 activity might generate a plastic state in developing trunk somites under which they can either enter carapace morphogenetic routes, possibly triggered by signals from the carapacial ridge, or still engage in the development of a standard tetrapod ribcage in the absence of those signals. In addition, trunk Lbx1 expression might play a later role in the formation of the lateral border of the carapace.

Reconstructing the evolutionary history of the oxytocin and vasotocin receptor gene family: Insights on whole genome duplication scenarios.

Vertebrate genome evolution remains a hotly debated topic, specifically as regards the number and the timing of putative rounds of whole genome duplication events. In this study, I sought to shed light to this conundrum through assessing the evolutionary history of the oxytocin/vasotocin receptor family. I performed ancestral analyses of the genomic segments containing oxytocin and vasotocin receptors (OTR-VTRs) by mapping them back to the reconstructed ancestral vertebrate/chordate karyotypes reported in five independent studies (Nakatani et al., 2007; Putnam et al., 2008; Smith and Keinath, 2015; Smith et al., 2018; Simakov et al., 2020) and found that two alternative scenarios can account for their evolution: one consistent with one round of whole genome duplication in the common ancestor of lampreys and gnathostomes, followed by segmental duplications in both lineages, and another consistent with two rounds of whole genome duplication, with the first occurring in the gnathostome-lamprey ancestor and the second in the jawed vertebrate ancestor. Combining the data reported here with synteny and phylogeny data reported in our previous study (Theofanopoulou et al., 2021), I put forward that a single round of whole genome duplication scenario is more consistent with the synteny and evolution of chromosomes where OTR-VTRs are encountered, without excluding the possibility of a scenario including two rounds of whole genome duplication. Although the analysis of one gene family is not able to capture the full complexity of vertebrate genome evolution, this study can provide solid insight, since the gene family used here has been meticulously analyzed for its genes' orthologous and paralogous relationships across species using high quality genomes.

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.

Analysis of lamprey meis genes reveals that conserved inputs from Hox, Meis and Pbx proteins control their expression in the hindbrain and neural tube.

Meis genes are known to play important roles in the hindbrain and neural crest cells of jawed vertebrates. To explore the roles of Meis genes in head development during evolution of vertebrates, we have identified four meis genes in the sea lamprey genome and characterized their patterns of expression and regulation, with a focus on the hindbrain and pharynx. Each of the lamprey meis genes displays temporally and spatially dynamic patterns of expression, some of which are coupled to rhombomeric domains in the developing hindbrain and select pharyngeal arches. Studies of Meis loci in mouse and zebrafish have identified enhancers that are bound by Hox and TALE (Meis and Pbx) proteins, implicating these factors in the direct regulation of Meis expression. We examined the lamprey meis loci and identified a series of cis-elements conserved between lamprey and jawed vertebrate meis genes. In transgenic reporter assays we demonstrated that these elements act as neural enhancers in lamprey embryos, directing reporter expression in appropriate domains when compared to expression of their associated endogenous meis gene. Sequence alignments reveal that these conserved elements are in similar relative positions of the meis loci and contain a series of consensus binding motifs for Hox and TALE proteins. This suggests that ancient Hox and TALE-responsive enhancers regulated expression of ancestral vertebrate meis genes in segmental domains in the hindbrain and have been retained in the meis loci during vertebrate evolution. The presence of conserved Meis, Pbx and Hox binding sites in these lamprey enhancers links Hox and TALE factors to regulation of lamprey meis genes in the developing hindbrain, indicating a deep ancestry for these regulatory interactions prior to the divergence of jawed and jawless vertebrates.

Intron size minimisation in teleosts.

BACKGROUND: Spliceosomal introns are parts of primary transcripts that are removed by RNA splicing. Although introns apparently do not contribute to the function of the mature transcript, in vertebrates they comprise the majority of the transcribed region increasing the metabolic cost of transcription. The persistence of long introns across evolutionary time suggests functional roles that can offset this metabolic cost. The teleosts comprise one of the largest vertebrate clades. They have unusually compact and variable genome sizes and provide a suitable system for analysing intron evolution. RESULTS: We have analysed intron lengths in 172 vertebrate genomes and show that teleost intron lengths are relatively short, highly variable and bimodally distributed. Introns that were long in teleosts were also found to be long in mammals and were more likely to be found in regulatory genes and to contain conserved sequences. Our results argue that intron length has decreased in parallel in a non-random manner throughout teleost evolution and represent a deviation from the ancestral state. CONCLUSION: Our observations indicate an accelerated rate of intron size evolution in the teleosts and that teleost introns can be divided into two classes by their length. Teleost intron sizes have evolved primarily as a side-effect of genome size evolution and small genomes are dominated by short introns (<256 base pairs). However, a non-random subset of introns has resisted this process across the teleosts and these are more likely have functional roles in all vertebrate clades.

Evolution of the nitric oxide synthase family in vertebrates and novel insights in gill development.

Nitric oxide (NO) is an ancestral key signalling molecule essential for life and has enormous versatility in biological systems, including cardiovascular homeostasis, neurotransmission and immunity. Although our knowledge of NO synthases (Nos), the enzymes that synthesize NO in vivo, is substantial, the origin of a large and diversified repertoire of nos gene orthologues in fishes with respect to tetrapods remains a puzzle. The recent identification of nos3 in the ray-finned fish spotted gar, which was considered lost in this lineage, changed this perspective. This finding prompted us to explore nos gene evolution, surveying vertebrate species representing key evolutionary nodes. This study provides noteworthy findings: first, nos2 experienced several lineage-specific gene duplications and losses. Second, nos3 was found to be lost independently in two different teleost lineages, Elopomorpha and Clupeocephala. Third, the expression of at least one nos paralogue in the gills of developing shark, bichir, sturgeon, and gar, but not in lamprey, suggests that nos expression in this organ may have arisen in the last common ancestor of gnathostomes. These results provide a framework for continuing research on nos genes' roles, highlighting subfunctionalization and reciprocal loss of function that occurred in different lineages during vertebrate genome duplications.

Conserved keratin gene clusters in ancient fish: An evolutionary seed for terrestrial adaptation.

Type I and type II keratins are subgroups of intermediate filament proteins that provide toughness to the epidermis and protect it from water loss. In terrestrial vertebrates, the keratin genes form two major clusters, clusters 1 and 2, each of which is dominated by type I and II keratin genes. By contrast, such clusters are not observed in teleost fish. Although the diversification of keratins is believed to have made a substantial contribution to terrestrial adaptation, its evolutionary process has not been clarified. Here, we performed a comprehensive genomic survey of the keratin genes of a broad range of vertebrates. As a result, we found that ancient fish lineages such as elephant shark, reedfish, spotted gar, and coelacanth share both keratin gene clusters. We also discovered an expansion of keratin genes that form a novel subcluster in reedfish. Syntenic and phylogenetic analyses revealed that two pairs of krt18/krt8 keratin genes were shared among all vertebrates, thus implying that they encode ancestral type I and II keratin protein sets. We further revealed that distinct keratin gene subclusters, which show specific expressions in the epidermis of adult amphibians, stemmed from canonical keratin genes in non-terrestrial ancestors. Molecular evolutionary analyses suggested that the selective constraints were relaxed in the adult epidermal subclusters of amphibians as well as the novel subcluster of reedfish. The results of the present study represent the process of diversification of keratins through a series of gene duplications that could have facilitated the terrestrial adaptation of vertebrates.

An ancestral role for Semaphorin3F-Neuropilin signaling in patterning neural crest within the new vertebrate head.

The origin of the vertebrate head is one of the great unresolved issues in vertebrate evolutionary developmental biology. Although many of the novelties in the vertebrate head and pharynx derive from the neural crest, it is still unknown how early vertebrates patterned the neural crest within the ancestral body plan they inherited from invertebrate chordates. Here, using a basal vertebrate, the sea lamprey, we show that homologs of Semaphorin3F (Sema3F) ligand and its Neuropilin (Nrp) receptors show complementary and dynamic patterns of expression that correlate with key periods of neural crest development (migration and patterning of cranial neural crest-derived structures). Using CRISPR/Cas9-mediated mutagenesis, we demonstrate that lamprey Sema3F is essential for patterning of neural crest-derived melanocytes, cranial ganglia and the head skeleton, but is not required for neural crest migration or patterning of trunk neural crest derivatives. Based on comparisons with jawed vertebrates, our results suggest that the deployment of Nrp-Sema3F signaling, along with other intercellular guidance cues, was pivotal in allowing early vertebrates to organize and pattern cranial neural crest cells into many of the hallmark structures that define the vertebrate head.

The evolution of ADAM gene family in eukaryotes.

The ADAM (A Disintegrin And Metalloprotease) gene family encodes proteins with adhesion and proteolytic functions. ADAM proteins are associated with diseases like cancers. Twenty ADAM genes have been identified in humans. However, little is known about the evolution of the family. We analyzed the repertoire of ADAM genes in a vast number of eukaryotic genomes to clarify the main gene copy number expansions. For the first time, we provide compelling evidence that early-branching green algae (Mamiellophyceae) have ADAM genes, suggesting that they originated in the last common ancestor of eukaryotes, before the split of plants, fungi and animals. The ADAM family expanded in early metazoans, with the most significative gene expansion happening during the first steps of vertebrate evolution. We concluded that most of mammal ADAM diversity can be explained by gene duplications in early bone fish. Our data suggest that ADAM genes were lost early in green plant evolution.