Hagfish and lamprey Hox genes reveal conservation of temporal colinearity in vertebrates

Global level of methylation in the sea lamprey (jawless vertebrate) genome is intermediate between invertebrate and jawed vertebrate genomes

In eukaryotes, cytosine methylation is a primary heritable epigenetic modification of the genome that regulates many cellular processes. In invertebrate, methylated cytosine generally located on specific genomic elements (e.g., gene bodies and silenced repetitive elements) to show a "mosaic" pattern. While in jawed vertebrate (teleost and tetrapod), highly methylated cytosine located genome-wide but only absence at regulatory regions (e.g., promoter and enhancer). Many studies imply that the evolution of DNA methylation reprogramming may have helped the transition from invertebrates to jawed vertebrates, but the detail remains largely elusive. In this study, we used the whole-genome bisulfite-sequencing technology to investigate the genome-wide methylation in three tissues (heart, muscle, and sperm) from the sea lamprey, an extant agnathan (jawless) vertebrate. Strikingly, we found that the methylation level of the sea lamprey is very similar to that in sea urchin (a deuterostome) and sea squirt (a chordate) invertebrates. In sum, the global pattern in sea lamprey is intermediate methylation level (around 30%), that is higher than methylation level in the genomes of pre-bilaterians and protostomes (1%-10%), but lower than methylation level appeared in jawed vertebrates (around 70%, teleost and tetrapod). We anticipate that, in addition to genetic dynamics such as genome duplications, epigenetic dynamics such as global methylation reprograming was also orchestrated toward the emergence and evolution of vertebrates.

Colonization of the ocean floor by jawless vertebrates across three mass extinctions

BACKGROUND

The deep (> 200 m) ocean floor is often considered to be a refugium of biodiversity; many benthic marine animals appear to share ancient common ancestry with nearshore and terrestrial relatives. Whether this pattern holds for vertebrates is obscured by a poor understanding of the evolutionary history of the oldest marine vertebrate clades. Hagfishes are jawless vertebrates that are either the living sister to all vertebrates or form a clade with lampreys, the only other surviving jawless fishes.

RESULTS

We use the hagfish fossil record and molecular data for all recognized genera to construct a novel hypothesis for hagfish relationships and diversification. We find that crown hagfishes persisted through three mass extinctions after appearing in the Permian ~ 275 Ma, making them one of the oldest living vertebrate lineages. In contrast to most other deep marine vertebrates, we consistently infer a deep origin of continental slope occupation by hagfishes that dates to the Paleozoic. Yet, we show that hagfishes have experienced marked body size diversification over the last hundred million years, contrasting with a view of this clade as morphologically stagnant.

CONCLUSION

Our results establish hagfishes as ancient members of demersal continental slope faunas and suggest a prolonged accumulation of deep sea jawless vertebrate biodiversity.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

Comparative analysis of Hmx expression and the distribution of neuronal somata in the trigeminal ganglion in lamprey and shark: insights into the homology of the trigeminal nerve branches and the evolutionary origin of the vertebrate jaw

The evolutionary origin of the jaw remains one of the most enigmatic events in vertebrate evolution. The trigeminal nerve is a key component for understanding jaw evolution, as it plays a crucial role as a sensorimotor interface for the effective manipulation of the jaw. This nerve is also found in the lamprey, an extant jawless vertebrate. The trigeminal nerve has three major branches in both the lamprey and jawed vertebrates. Although each of these branches was classically thought to be homologous between these two taxa, this homology is now in doubt. In the present study, we compared expression patterns of Hmx, a candidate genetic marker of the mandibular nerve (rV3, the third branch of the trigeminal nerve in jawed vertebrates), and the distribution of neuronal somata of trigeminal nerve branches in the trigeminal ganglion in lamprey and shark. We first confirmed the conserved expression pattern of Hmx1 in the shark rV3 neuronal somata, which are distributed in the caudal part of the trigeminal ganglion. By contrast, lamprey Hmx genes showed peculiar expression patterns, with expression in the ventrocaudal part of the trigeminal ganglion similar to Hmx1 expression in jawed vertebrates, which labeled the neuronal somata of the second branch. Based on these results, we propose two alternative hypotheses regarding the homology of the trigeminal nerve branches, providing new insights into the evolutionary origin of the vertebrate jaw.

Bioinformatic and fine-scale chromosomal mapping reveal the nature and evolution of eliminated chromosomes in the Japanese hagfish, Eptatretus burgeri, through analysis of repetitive DNA families

In the Japanese hagfish, Eptatretus burgeri, approximately 21% of the genomic DNA in germ cells (2n = 52) consists of 16 chromosomes (eliminated [E]-chromosomes) that are eliminated from presumptive somatic cells (2n = 36). To uncover the eliminated genome (E-genome), we have identified 16 eliminated repetitive DNA families from eight hagfish species, with 11 of these repeats being selectively amplified in the germline genome of E. burgeri. Furthermore, we have demonstrated that six of these sequences, namely EEEb1-6, are exclusively localized on all 16 E-chromosomes. This has led to the hypothesis that the eight pairs of E-chromosomes are derived from one pair of ancestral chromosomes via multiple duplication events over a prolonged evolutionary period. NGS analysis has recently facilitated the re-assembly of two distinct draft genomes of E. burgeri, derived from the testis and liver. This advancement allows for the prediction of not only nonrepetitive eliminated sequences but also over 100 repetitive and eliminated sequences, accomplished through K-mer-based analysis. In this study, we report four novel eliminated repetitive DNA sequences (designated as EEEb7-10) and confirm the relative chromosomal localization of all eliminated repeats (EEEb1-10) by fluorescence in situ hybridization (FISH). With the exception of EEEb10, all sequences were exclusively detected on EEEb1-positive chromosomes. Surprisingly, EEEb10 was detected as an intense signal on EEEb1-positive chromosomes and as a scattered signal on other chromosomes in germ cells. The study further divided the eight pairs of E-chromosomes into six groups based on the signal distribution of each DNA family, and fiber-FISH experiments showed that the EEEb2-10 family was dispersed in the EEEb1-positive extended chromatin fiber. These findings provide new insights into the mechanisms underlying chromosome elimination and the evolution of E-chromosomes, supporting our previous hypothesis.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

Comparative expression analysis of the Atoh7 gene regulatory network in the mouse and chicken auditory hindbrain

The mammalian and avian auditory brainstem likely arose by independent evolution. To compare the underlying molecular mechanisms, we focused on Atoh7, as its expression pattern in the mammalian hindbrain is restricted to bushy cells in the ventral cochlear nucleus. We thereby took advantage of an Atoh7 centered gene regulatory network (GRN) in the retina including upstream regulators, Hes1 and Pax6, and downstream targets, Ebf3 and Eya2. In situ hybridization demonstrated for the latter four genes broad expression in all three murine cochlear nuclei at postnatal days (P) 4 and P30, contrasting the restricted expression of Atoh7. In chicken, all five transcription factors were expressed in all auditory hindbrain nuclei at embryonic day (E) 13 and P14. Notably, all five genes showed graded expression in the embryonic nucleus magnocellularis (NM). Atoh7 was highly expressed in caudally located neurons, whereas the other four transcription factors were highly expressed in rostrally located neurons. Thus, Atoh7 shows a strikingly different expression between the mammalian and avian auditory hindbrain. This together with the consistent absence of graded expression of GRN components in developing mammalian nuclei provide the first molecular support to the current view of convergent evolution as a major mechanism in the amniote auditory hindbrain. The graded expression of five transcription factors specifically in the developing NM confirms this nucleus as a central organizer of tonotopic features in birds. Finally, the expression of all five retinal GRN components in the auditory system suggests co-options of genes for development of sensory systems of distinct modalities.

Gonadal Transcriptome Analysis and Sequence Characterization of Sex-Related Genes in Cranoglanis bouderius

In China, the Cranoglanis bouderius is classified as a national class II-protected animal. The development of C. bouderius populations has been affected by a variety of factors over the past few decades, with severe declines occurring. Considering the likelihood of continued population declines of the C. bouderius in the future, it is critical to investigate the currently unknown characteristics of gonadal differentiation and sex-related genes for C. bouderius conservation. In this study, the Illumina sequencing platform was used to sequence the gonadal transcriptome of the C. bouderius to identify the pathways and genes related to gonadal development and analyze the expression differences in the gonads. A total of 12,002 DEGs were identified, with 7220 being significantly expressed in the ovary and 4782 being significantly expressed in the testis. According to the functional enrichment results, the cell cycle, RNA transport, apoptosis, Wnt signaling pathway, p53 signaling pathway, and prolactin signaling pathway play important roles in sex development in the C. bouderius. Furthermore, the sequence characterization and evolutionary analysis revealed that AMH, DAX1, NANOS1, and AR of the C. bouderius are highly conserved. Specifically, the qRT-PCR results from various tissues showed significant differences in AMH, DAX1, NANOS1, and AR expression levels in the gonads of both sexes of C. bouderius. These analyses indicated that AMH, DAX1, NANOS1, and AR may play important roles in the differentiation and development of C. bouderius gonads. To our best knowledge, this study is the first to analyze the C. bouderius gonadal transcriptome and identify the structures of sex-related genes, laying the foundation for future research.

Novel selectively amplified DNA sequences in the germline genome of the Japanese hagfish, Eptatretus burgeri

In the Japanese hagfish Eptatretus burgeri, 16 chromosomes (eliminated [E]-chromosomes) have been lost in somatic cells (2n = 36), which is equivalent to approx. 21% of the genomic DNA in germ cells (2n = 52). At least seven of the 12 eliminated repetitive DNA families isolated in eight hagfish species were selectively amplified in the germline genome of this species. One of them, EEEb1 (eliminated element of E. burgeri 1) is exclusively localized on all E-chromosomes. Herein, we identified four novel eliminated repetitive DNA families (named EEEb3-6) through PCR amplification and suppressive subtractive hybridization (SSH) combined with Southern-blot hybridization. EEEb3 was mosaic for 5S rDNA and SINE elements. EEEb4 was GC-rich repeats and has one pair of direct and inverted repeats, whereas EEEb5 and EEEb6 were AT-rich repeats with one pair and two pairs of sub-repeats, respectively. Interestingly, all repeat classes except EEEb3 were transcribed in the testes, although no open reading frames (ORF) were identified. We conducted fluorescence in situ hybridization (FISH) to examine the chromosomal localizations of EEEb3-6 and EEEb2, which was previously isolated from the germline genome of E. burgeri. All sequences were only found on all EEEb1-positive E-chromosomes. Copy number estimation of the repeated elements by slot-blot hybridization revealed that (i) the EEEb1-6 family members occupied 39.9% of the total eliminated DNA, and (ii) a small number of repeats were retained in somatic cells, suggesting that there is incomplete elimination of the repeated elements. These results provide new insights into the mechanisms involved in the chromosome elimination and the evolution of E-chromosomes.

Discovery of prolactin-like in lamprey: Role in osmoregulation and new insight into the evolution of the growth hormone/prolactin family

We used a representative of one of the oldest extant vertebrate lineages (jawless fish or agnathans) to investigate the early evolution and function of the growth hormone (GH)/prolactin (PRL) family. We identified a second member of the GH/PRL family in an agnathan, the sea lamprey (Petromyzon marinus). Structural, phylogenetic, and synteny analyses supported the identification of this hormone as prolactin-like (PRL-L), which has led to added insight into the evolution of the GH/PRL family. At least two ancestral genes were present in early vertebrates, which gave rise to distinct GH and PRL-L genes in lamprey. A series of gene duplications, gene losses, and chromosomal rearrangements account for the diversity of GH/PRL-family members in jawed vertebrates. Lamprey PRL-L is produced in the proximal pars distalis of the pituitary and is preferentially bound by the lamprey PRL receptor, whereas lamprey GH is preferentially bound by the lamprey GH receptor. Pituitary PRL-L messenger RNA (mRNA) levels were low in larvae, then increased significantly in mid-metamorphic transformers (stage 3); thereafter, levels subsided in final-stage transformers and metamorphosed juveniles. The abundance of PRL-L mRNA and immunoreactive protein increased in the pituitary of juveniles under hypoosmotic conditions, and treatment with PRL-L blocked seawater-associated inhibition of freshwater ion transporters. These findings clarify the origin and divergence of GH/PRL family genes in early vertebrates and reveal a function of PRL-L in osmoregulation of sea lamprey, comparable to a role of PRLs that is conserved in jawed vertebrates.

Echiuran Hox genes provide new insights into the correspondence between Hox subcluster organization and collinearity pattern

In many bilaterians, Hox genes are generally clustered along the chromosomes and expressed in spatial and temporal order. In vertebrates, the expression of Hox genes follows a whole-cluster spatio-temporal collinearity (WSTC) pattern, whereas in some invertebrates the expression of Hox genes exhibits a subcluster-level spatio-temporal collinearity pattern. In bilaterians, the diversity of collinearity patterns and the cause of collinearity differences in Hox gene expression remain poorly understood. Here, we investigate genomic organization and expression pattern of Hox genes in the echiuran worm Urechis unicinctus (Annelida, Echiura). Urechis unicinctus has a split cluster with four subclusters divided by non-Hox genes: first subcluster (Hox1 and Hox2), second subcluster (Hox3), third subcluster (Hox4, Hox5, Lox5, Antp and Lox4), fourth subcluster (Lox2 and Post2). The expression of U. unicinctus Hox genes shows a subcluster-based whole-cluster spatio-temporal collinearity (S-WSTC) pattern: the anterior-most genes in each subcluster are activated in a spatially and temporally colinear manner (reminiscent of WSTC), with the subsequent genes in each subcluster then being very similar to their respective anterior-most subcluster gene. Combining genomic organization and expression profiles of Hox genes in different invertebrate lineages, we propose that the spatio-temporal collinearity of invertebrate Hox is subcluster-based.

Evolution of a chordate-specific mechanism for myoblast fusion

Vertebrate myoblast fusion allows for multinucleated muscle fibers to compound the size and strength of mononucleated cells, but the evolution of this important process is unknown. We investigated the evolutionary origins and function of membrane-coalescing agents Myomaker and Myomixer in various groups of chordates. Here, we report that Myomaker likely arose through gene duplication in the last common ancestor of tunicates and vertebrates, while Myomixer appears to have evolved de novo in early vertebrates. Functional tests revealed a complex evolutionary history of myoblast fusion. A prevertebrate phase of muscle multinucleation driven by Myomaker was followed by the later emergence of Myomixer that enables the highly efficient fusion system of vertebrates. Evolutionary comparisons between vertebrate and nonvertebrate Myomaker revealed key structural and mechanistic insights into myoblast fusion. Thus, our findings suggest an evolutionary model of chordate fusogens and illustrate how new genes shape the emergence of novel morphogenetic traits and mechanisms.

Thyroid and endostyle development in cyclostomes provides new insights into the evolutionary history of vertebrates

Background

The endostyle is an epithelial exocrine gland found in non-vertebrate chordates (amphioxi and tunicates) and the larvae of modern lampreys. It is generally considered to be an evolutionary precursor of the thyroid gland of vertebrates. Transformation of the endostyle into the thyroid gland during the metamorphosis of lampreys is thus deemed to be a recapitulation of a past event in vertebrate evolution. In 1906, Stockard reported that the thyroid gland in hagfish, the sister cyclostome group of lampreys, develops through an endostyle-like primordium, strongly supporting the plesiomorphy of the lamprey endostyle. However, the findings in hagfish thyroid development were solely based on this single study, and these have not been confirmed by modern molecular, genetic, and morphological data pertaining to hagfish thyroid development over the last century.

Results

Here, we showed that the thyroid gland of hagfish undergoes direct development from the ventrorostral pharyngeal endoderm, where the previously described endostyle-like primordium was not found. The developmental pattern of the hagfish thyroid, including histological features and regulatory gene expression profiles, closely resembles that found in modern jawed vertebrates (gnathostomes). Meanwhile, as opposed to gnathostomes but similar to non-vertebrate chordates, lamprey and hagfish share a broad expression domain of Nkx2-1/2-4, a key regulatory gene, in the pharyngeal epithelium during early developmental stages.

Conclusions

Based on the direct development of the thyroid gland both in hagfish and gnathostomes, and the shared expression profile of thyroid-related transcription factors in the cyclostomes, we challenge the plesiomorphic status of the lamprey endostyle and propose an alternative hypothesis where the lamprey endostyle could be obtained secondarily in crown lampreys.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01282-7.

Forebrain Architecture and Development in Cyclostomes, with Reference to the Early Morphology and Evolution of the Vertebrate Head

The vertebrate head and brain are characterized by highly complex morphological patterns. The forebrain, the most anterior division of the brain, is subdivided into the diencephalon, hypothalamus, and telencephalon from the neuromeric subdivision into prosomeres. Importantly, the telencephalon contains the cerebral cortex, which plays a key role in higher order cognitive functions in humans. To elucidate the evolution of the forebrain regionalization, comparative analyses of the brain development between extant jawed and jawless vertebrates are crucial. Cyclostomes - lampreys and hagfishes - are the only extant jawless vertebrates, and diverged from jawed vertebrates (gnathostomes) over 500 million years ago. Previous developmental studies on the cyclostome brain were conducted mainly in lampreys because hagfish embryos were rarely available. Although still scarce, the recent availability of hagfish embryos has propelled comparative studies of brain development and gene expression. By integrating findings with those of cyclostomes and fossil jawless vertebrates, we can depict the morphology, developmental mechanism, and even the evolutionary path of the brain of the last common ancestor of vertebrates. In this review, we summarize the development of the forebrain in cyclostomes and suggest what evolutionary changes each cyclostome lineage underwent during brain evolution. In addition, together with recent advances in the head morphology in fossil vertebrates revealed by CT scanning technology, we discuss how the evolution of craniofacial morphology and the changes of the developmental mechanism of the forebrain towards crown gnathostomes are causally related.

Genetic Mechanism for the Cyclostome Cerebellar Neurons Reveals Early Evolution of the Vertebrate Cerebellum

The vertebrate cerebellum arises at the dorsal part of rhombomere 1, induced by signals from the isthmic organizer. Two major cerebellar neuronal subtypes, granule cells (excitatory) and Purkinje cells (inhibitory), are generated from the anterior rhombic lip and the ventricular zone, respectively. This regionalization and the way it develops are shared in all extant jawed vertebrates (gnathostomes). However, very little is known about early evolution of the cerebellum. The lamprey, an extant jawless vertebrate lineage or cyclostome, possesses an undifferentiated, plate-like cerebellum, whereas the hagfish, another cyclostome lineage, is thought to lack a cerebellum proper. In this study, we found that hagfish Atoh1 and Wnt1 genes are co-expressed in the rhombic lip, and Ptf1a is expressed ventrally to them, confirming the existence of r1's rhombic lip and the ventricular zone in cyclostomes. In later stages, lamprey Atoh1 is downregulated in the posterior r1, in which the NeuroD increases, similar to the differentiation process of cerebellar granule cells in gnathostomes. Also, a continuous Atoh1-positive domain in the rostral r1 is reminiscent of the primordium of valvula cerebelli of ray-finned fishes. Lastly, we detected a GAD-positive domain adjacent to the Ptf1a-positive ventricular zone in lampreys, suggesting that the Ptf1a-positive cells differentiate into some GABAergic inhibitory neurons such as Purkinje and other inhibitory neurons like in gnathostomes. Altogether, we conclude that the ancestral genetic programs for the formation of a distinct cerebellum were established in the last common ancestor of vertebrates.

Reconstruction of proto-vertebrate, proto-cyclostome and proto-gnathostome genomes provides new insights into early vertebrate evolution

Ancient polyploidization events have had a lasting impact on vertebrate genome structure, organization and function. Some key questions regarding the number of ancient polyploidization events and their timing in relation to the cyclostome-gnathostome divergence have remained contentious. Here we generate de novo long-read-based chromosome-scale genome assemblies for the Japanese lamprey and elephant shark. Using these and other representative genomes and developing algorithms for the probabilistic macrosynteny model, we reconstruct high-resolution proto-vertebrate, proto-cyclostome and proto-gnathostome genomes. Our reconstructions resolve key questions regarding the early evolutionary history of vertebrates. First, cyclostomes diverged from the lineage leading to gnathostomes after a shared tetraploidization (1R) but before a gnathostome-specific tetraploidization (2R). Second, the cyclostome lineage experienced an additional hexaploidization. Third, 2R in the gnathostome lineage was an allotetraploidization event, and biased gene loss from one of the subgenomes shaped the gnathostome genome by giving rise to remarkably conserved microchromosomes. Thus, our reconstructions reveal the major evolutionary events and offer new insights into the origin and evolution of vertebrate genomes.

Improved Understanding of the Role of Gene and Genome Duplications in Chordate Evolution With New Genome and Transcriptome Sequences

Comparative approaches to understanding chordate genomes have uncovered a significant role for gene duplications, including whole genome duplications (WGDs), giving rise to and expanding gene families. In developmental biology, gene families created and expanded by both tandem and WGDs are paramount. These genes, often involved in transcription and signalling, are candidates for underpinning major evolutionary transitions because they are particularly prone to retention and subfunctionalisation, neofunctionalisation, or specialisation following duplication. Under the subfunctionalisation model, duplication lays the foundation for the diversification of paralogues, especially in the context of gene regulation. Tandemly duplicated paralogues reside in the same regulatory environment, which may constrain them and result in a gene cluster with closely linked but subtly different expression patterns and functions. Ohnologues (WGD paralogues) often diversify by partitioning their expression domains between retained paralogues, amidst the many changes in the genome during rediploidisation, including chromosomal rearrangements and extensive gene losses. The patterns of these retentions and losses are still not fully understood, nor is the full extent of the impact of gene duplication on chordate evolution. The growing number of sequencing projects, genomic resources, transcriptomics, and improvements to genome assemblies for diverse chordates from non-model and under-sampled lineages like the coelacanth, as well as key lineages, such as amphioxus and lamprey, has allowed more informative comparisons within developmental gene families as well as revealing the extent of conserved synteny across whole genomes. This influx of data provides the tools necessary for phylogenetically informed comparative genomics, which will bring us closer to understanding the evolution of chordate body plan diversity and the changes underpinning the origin and diversification of vertebrates.

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.

Evolutionary divergence of a Hoxa2b hindbrain enhancer in syngnathids mimics results of functional assays

Hoxa2 genes provide critical patterning signals during development, and their regulation and function have been extensively studied. We report a previously uncharacterized significant sequence divergence of a highly conserved hindbrain hoxa2b enhancer element in the family syngnathidae (pipefishes, seahorses, pipehorses, seadragons). We compared the hox cis-regulatory element variation in the Gulf pipefish and two species of seahorse against eight other species of fish, as well as human and mouse. We annotated the hoxa2b enhancer element binding sites across three species of seahorse, four species of pipefish, and one species of ghost pipefish. Finally, we performed in situ hybridization analysis of hoxa2b expression in Gulf pipefish embryos. We found that all syngnathid fish examined share a modified rhombomere 4 hoxa2b enhancer element, despite the fact that this element has been found to be highly conserved across all vertebrates examined previously. Binding element sequence motifs and spacing between binding elements have been modified for the hoxa2b enhancer in several species of pipefish and seahorse, and that the loss of the Prep/Meis binding site and further space shortening happened after ghost pipefish split from the rest of the syngnathid clade. We showed that expression of this gene in rhombomere 4 is lower relative to the surrounding rhombomeres in developing Gulf pipefish embryos, reflecting previously published functional tests for this enhancer. Our findings highlight the benefits of studying highly derived, diverse taxa for understanding of gene regulatory evolution and support the hypothesis that natural mutations can occur in deeply conserved pathways in ways potentially related to phenotypic diversity.

Genome-Wide Identification, Comparison, and Expression Analysis of Transcription Factors in Ascidian Styela clava

Tunicates include diverse species, as they are model animals for evolutionary developmental biology study. The embryonic development of tunicates is known to be extensively regulated by transcription factors (TFs). Styela clava, the globally distributed invasive tunicate, exhibits a strong capacity for environmental adaptation. However, the TFs were not systematically identified and analyzed. In this study, we reported 553 TFs categorized into 60 families from S. clava, based on the whole genome data. Comparison of TFs analysis among the tunicate species revealed that the gene number in the zinc finger superfamily displayed the most significant discrepancy, indicating this family was under the highly evolutionary selection and might be related to species differentiation and environmental adaptation. The greatest number of TFs was discovered in the Cys2His2-type zinc finger protein (zf-C2H2) family in S. clava. From the point of temporal view, more than half the TFs were expressed at the early embryonic stage. The expression correlation analysis revealed the existence of a transition for TFs expression from early embryogenesis to the later larval development in S. clava. Eight Hox genes were identified to be located on one chromosome, exhibiting different arrangement and expression patterns, compared to Ciona robusta (C. intestinalis type A). In addition, a total of 23 forkhead box (fox) genes were identified in S. clava, and their expression profiles referred to their potential roles in neurodevelopment and sensory organ development. Our data, thus, provides crucial clues to the potential functions of TFs in development and environmental adaptation in the leathery sea squirt.

Cytidine deaminase 2 is required for VLRB antibody gene assembly in lampreys

The antibodies of jawless vertebrates consist of leucine-rich repeat arrays encoded by somatically assembled VLRB genes. It is unknown how the incomplete germline VLRB loci are converted into functional antibody genes during B lymphocyte development in lampreys. In Lampetra planeri larvae lacking the cytidine deaminase CDA2 gene, VLRB assembly fails, whereas the T lineage-associated VLRA and VLRC antigen receptor gene assemblies occur normally. Thus, CDA2 acts in a B cell lineage-specific fashion to support the somatic diversification of VLRB antibody genes. CDA2 is closely related to activation-induced cytidine deaminase (AID), which is essential for the elaboration of immunoglobulin gene repertoires in jawed vertebrates. Our results thus identify a convergent mechanism of antigen receptor gene assembly and diversification that independently evolved in the two sister branches of vertebrates.

Evo-devo studies of cyclostomes and the origin and evolution of jawed vertebrates

Modern vertebrates consist of two sister groups: cyclostomes and gnathostomes. Cyclostomes are a monophyletic jawless group that can be further divided into hagfishes and lampreys, which show conspicuously different developmental and morphological patterns. However, during early pharyngula development, there appears to be a stage when the embryos of hagfishes and lampreys resemble each other by showing an "ancestral" craniofacial pattern; this pattern enables morphological comparison of hagfish and lamprey craniofacial development at late stages. This cyclostome developmental pattern, or more accurately, this developmental pattern of the jawless grade of vertebrates in early pharyngula was very likely shared by the gnathostome stem before the division of the nasohypophyseal placode led to the jaw and paired nostrils. The craniofacial pattern of the modern jawed vertebrates seems to have begun in fossil ostracoderms (including galeaspids), and was completed by the early placoderm lineages. The transition from jawless to jawed vertebrates was thus driven by heterotopy of development, mainly caused by separation and shift of ectodermal placodes and resultant ectomesenchymal distribution, and shifts of the epithelial-mesenchymal interactions that underlie craniofacial differentiation. Thus, the evolution of the jaw was not a simple modification of the mandibular arch, but a heterotopic shift of the developmental interactions involving not only the mandibular arch, but also the premandibular region rostral to the mandibular arch.

Chromosome-level genome assembly of Lethenteron reissneri provides insights into lamprey evolution

The reissner lamprey Lethenteron reissneri, belonging to the class Cyclostomata, serves as a bridge between invertebrates and jawed vertebrates, and is considered the sister group of jawed vertebrates. However, despite this evolutionary significance, the genetic mechanisms underlying the adaptive evolution of the lamprey lineage remain unclear. Here, we assembled a 1.06 Gb chromosome-level draft genome of L. reissneri, with 72 chromosomes (ranging in length from 4.5 Mb to 25.9 Mb) and a scaffold N50 length of 13.23 Mb. Genome quality comparisons revealed that the reissner lamprey genome has higher completeness and contiguity than the previously published sea lamprey and Japanese lamprey genomes. Moreover, reissner lamprey, sea lamprey, and Japanese lamprey species share similar transposable element profiles and Hox gene cluster compositions, suggesting that a burst of transposable element activity and whole genome duplication occurred before their divergence. Additionally, the Lip gene copy numbers, which have been studied for their functions in the host defence system, were found to be expanded uniquely in lamprey lineages, suggesting key roles for these genes in lamprey evolution and adaptation. We also identified two neural-related genes, Nrn1 and Unc13a, with copy number expansions in jawed vertebrates, which may be functionally relevant to the origin of lamprey brains. Hence, this study not only provides the first chromosome-level reference genome for Cyclostomata, but also highlights features of the unique biology and adaptive evolution of the lamprey lineage.

Programmed DNA Elimination in Vertebrates

Over the last few decades, an increasing number of vertebrate taxa have been identified that undergo programmed genome rearrangement, or programmed DNA loss, during development. In these organisms, the genome of germ cells is often reproducibly different from the genome of all other cells within the body. Although we clearly have not identified all vertebrate taxa that undergo programmed genome loss, the list of species known to undergo loss now represents ∼10% of vertebrate species, including several basally diverging lineages. Recent studies have shed new light on the targets and mechanisms of DNA loss and their association with canonical modes of DNA silencing. Ultimately, expansion of these studies into a larger collection of taxa will aid in reconstructing patterns of shared/independent ancestry of programmed DNA loss in the vertebrate lineage, as well as more recent evolutionary events that have shaped the structure and content of eliminated DNA.

Hox gene expression during development of the phoronid Phoronopsis harmeri

Background

Phoronida is a small group of marine worm-like suspension feeders, which together with brachiopods and bryozoans form the clade Lophophorata. Although their development is well studied on the morphological level, data regarding gene expression during this process are scarce and restricted to the analysis of relatively few transcription factors. Here, we present a description of the expression patterns of Hox genes during the embryonic and larval development of the phoronid Phoronopsis harmeri.

Results

We identified sequences of eight Hox genes in the transcriptome of Ph. harmeri and determined their expression pattern during embryonic and larval development using whole mount in situ hybridization. We found that none of the Hox genes is expressed during embryonic development. Instead their expression is initiated in the later developmental stages, when the larval body is already formed. In the investigated initial larval stages the Hox genes are expressed in the non-collinear manner in the posterior body of the larvae: in the telotroch and the structures that represent rudiments of the adult worm. Additionally, we found that certain head-specific transcription factors are expressed in the oral hood, apical organ, preoral coelom, digestive system and developing larval tentacles, anterior to the Hox-expressing territories.

Conclusions

The lack of Hox gene expression during early development of Ph. harmeri indicates that the larval body develops without positional information from the Hox patterning system. Such phenomenon might be a consequence of the evolutionary intercalation of the larval form into an ancestral life cycle of phoronids. The observed Hox gene expression can also be a consequence of the actinotrocha representing a “head larva”, which is composed of the most anterior body region that is devoid of Hox gene expression. Such interpretation is further supported by the expression of head-specific transcription factors. This implies that the Hox patterning system is used for the positional information of the trunk rudiments and is, therefore, delayed to the later larval stages. We propose that a new body form was intercalated to the phoronid life cycle by precocious development of the anterior structures or by delayed development of the trunk rudiment in the ancestral phoronid larva.

Divergent genes encoding the putative receptors for growth hormone and prolactin in sea lamprey display distinct patterns of expression

Growth hormone receptor (GHR) and prolactin receptor (PRLR) in jawed vertebrates were thought to arise after the divergence of gnathostomes from a basal vertebrate. In this study we have identified two genes encoding putative GHR and PRLR in sea lamprey (Petromyzon marinus) and Arctic lamprey (Lethenteron camtschaticum), extant members of one of the oldest vertebrate groups, agnathans. Phylogenetic analysis revealed that lamprey GHR and PRLR cluster at the base of gnathostome GHR and PRLR clades, respectively. This indicates that distinct GHR and PRLR arose prior to the emergence of the lamprey branch of agnathans. In the sea lamprey, GHR and PRLR displayed a differential but overlapping pattern of expression; GHR had high expression in liver and heart tissues, whereas PRLR was expressed highly in the brain and moderately in osmoregulatory tissues. Branchial PRLR mRNA levels were significantly elevated by stage 5 of metamorphosis and remained elevated through stage 7, whereas levels of GHR mRNA were only elevated in the final stage (7). Branchial expression of GHR increased following seawater (SW) exposure of juveniles, but expression of PRLR was not significantly altered. The results indicate that GHR and PRLR may both participate in metamorphosis and that GHR may mediate SW acclimation.

Current Principles of Motor Control, with Special Reference to Vertebrate Locomotion

The vertebrate control of locomotion involves all levels of the nervous system from cortex to the spinal cord. Here, we aim to cover all main aspects of this complex behavior, from the operation of the microcircuits in the spinal cord to the systems and behavioral levels and extend from mammalian locomotion to the basic undulatory movements of lamprey and fish. The cellular basis of propulsion represents the core of the control system, and it involves the spinal central pattern generator networks (CPGs) controlling the timing of different muscles, the sensory compensation for perturbations, and the brain stem command systems controlling the level of activity of the CPGs and the speed of locomotion. The forebrain and in particular the basal ganglia are involved in determining which motor programs should be recruited at a given point of time and can both initiate and stop locomotor activity. The propulsive control system needs to be integrated with the postural control system to maintain body orientation. Moreover, the locomotor movements need to be steered so that the subject approaches the goal of the locomotor episode, or avoids colliding with elements in the environment or simply escapes at high speed. These different aspects will all be covered in the review.

Microsyntenic Clusters Reveal Conservation of lncRNAs in Chordates Despite Absence of Sequence Conservation

Homologous long non-coding RNAs (lncRNAs) are elusive to identify by sequence similarity due to their fast-evolutionary rate. Here we develop LincOFinder, a pipeline that finds conserved intergenic lncRNAs (lincRNAs) between distant related species by means of microsynteny analyses. Using this tool, we have identified 16 bona fide homologous lincRNAs between the amphioxus and human genomes. We characterized and compared in amphioxus and Xenopus the expression domain of one of them, Hotairm1, located in the anterior part of the Hox cluster. In addition, we analyzed the function of this lincRNA in Xenopus, showing that its disruption produces a severe headless phenotype, most probably by interfering with the regulation of the Hox cluster. Our results strongly suggest that this lincRNA has probably been regulating the Hox cluster since the early origin of chordates. Our work pioneers the use of syntenic searches to identify non-coding genes over long evolutionary distances and helps to further understand lncRNA evolution.

Evolution of vertebrate spinal cord patterning

The vertebrate spinal cord is organized across three developmental axes, anterior-posterior (AP), dorsal-ventral (DV), and medial-lateral (ML). Patterning of these axes is regulated by canonical intercellular signaling pathways: the AP axis by Wnt, fibroblast growth factor, and retinoic acid (RA), the DV axis by Hedgehog, Tgfβ, and Wnt, and the ML axis where proliferation is controlled by Notch. Developmental time plays an important role in which signal does what and when. Patterning across the three axes is not independent, but linked by interactions between signaling pathway components and their transcriptional targets. Combined this builds a sophisticated organ with many different types of cell in specific AP, DV, and ML positions. Two living lineages share phylum Chordata with vertebrates, amphioxus, and tunicates, while the jawless fish such as lampreys, survive as the most basally divergent vertebrate lineage. Genes and mechanisms shared between lampreys and other vertebrates tell us what predated vertebrates, while those also shared with other chordates tell us what evolved early in chordate evolution. Between these lie vertebrate innovations: genetic and developmental changes linked to evolution of new morphology. These include gene duplications, differences in how signals are received, and new regulatory connections between signaling pathways and their target genes.

The constrained architecture of mammalian Hox gene clusters

Vertebrate Hox genes are clustered. This organization has a functional relevance, as the transcription of each gene in time and space depends upon its relative position within the gene cluster. Hox clusters display a high organization, and all genes are transcribed from the same DNA strand. Here, we investigate the importance of this uniform transcriptional polarity by engineering alleles where one or several transcription units are inverted, with or without a CTCF site. We observe that inversions are likely detrimental to the proper implementation of this genetic system. We propose that the enhanced organization of Hox clusters in vertebrates evolved in conjunction with the emergence of global gene regulation to optimize a coordinated response of selected subsets of target genes.

In many animal species with a bilateral symmetry, Hox genes are clustered either at one or at several genomic loci. This organization has a functional relevance, as the transcriptional control applied to each gene depends upon its relative position within the gene cluster. It was previously noted that vertebrate Hox clusters display a much higher level of genomic organization than their invertebrate counterparts. The former are always more compact than the latter, they are generally devoid of repeats and of interspersed genes, and all genes are transcribed by the same DNA strand, suggesting that particular factors constrained these clusters toward a tighter structure during the evolution of the vertebrate lineage. Here, we investigate the importance of uniform transcriptional orientation by engineering several alleles within the HoxD cluster, such as to invert one or several transcription units, with or without a neighboring CTCF site. We observe that the association between the tight structure of mammalian Hox clusters and their regulation makes inversions likely detrimental to the proper implementation of this complex genetic system. We propose that the consolidation of Hox clusters in vertebrates, including transcriptional polarity, evolved in conjunction with the emergence of global gene regulation via the flanking regulatory landscapes, to optimize a coordinated response of selected subsets of target genes in cis.

An atlas of anterior hox gene expression in the embryonic sea lamprey head: Hox-code evolution in vertebrates

In the hindbrain and the adjacent cranial neural crest (NC) cells of jawed vertebrates (gnathostomes), nested and segmentally-restricted domains of Hox gene expression provide a combinatorial Hox-code for specifying regional properties during head development. Extant jawless vertebrates, such as the sea lamprey (Petromyzon marinus), can provide insights into the evolution and diversification of this Hox-code in vertebrates. There is evidence for gnathostome-like spatial patterns of Hox expression in lamprey; however, the expression domains of the majority of lamprey hox genes from paralogy groups (PG) 1-4 are yet to be characterized, so it is unknown whether they are coupled to hindbrain segments (rhombomeres) and NC. In this study, we systematically describe the spatiotemporal expression of all 14 sea lamprey hox genes from PG1-PG4 in the developing hindbrain and pharynx to investigate the extent to which their expression conforms to the archetypal gnathostome hindbrain and pharyngeal hox-codes. We find many similarities in Hox expression between lamprey and gnathostome species, particularly in rhombomeric domains during hindbrain segmentation and in the cranial neural crest, enabling inference of aspects of Hox expression in the ancestral vertebrate embryonic head. These data are consistent with the idea that a Hox regulatory network underlying hindbrain segmentation is a pan vertebrate trait. We also reveal differences in hindbrain domains at later stages, as well as expression in the endostyle and in pharyngeal arch (PA) 1 mesoderm. Our analysis suggests that many Hox expression domains that are observed in extant gnathostomes were present in ancestral vertebrates but have been partitioned differently across Hox clusters in gnathostome and cyclostome lineages after duplication.

A Hox-TALE regulatory circuit for neural crest patterning is conserved across vertebrates

In jawed vertebrates (gnathostomes), Hox genes play an important role in patterning head and jaw formation, but mechanisms coupling Hox genes to neural crest (NC) are unknown. Here we use cross-species regulatory comparisons between gnathostomes and lamprey, a jawless extant vertebrate, to investigate conserved ancestral mechanisms regulating Hox2 genes in NC. Gnathostome Hoxa2 and Hoxb2 NC enhancers mediate equivalent NC expression in lamprey and gnathostomes, revealing ancient conservation of Hox upstream regulatory components in NC. In characterizing a lamprey hoxα2 NC/hindbrain enhancer, we identify essential Meis, Pbx, and Hox binding sites that are functionally conserved within Hoxa2/Hoxb2 NC enhancers. This suggests that the lamprey hoxα2 enhancer retains ancestral activity and that Hoxa2/Hoxb2 NC enhancers are ancient paralogues, which diverged in hindbrain and NC activities. This identifies an ancestral mechanism for Hox2 NC regulation involving a Hox-TALE regulatory circuit, potentiated by inputs from Meis and Pbx proteins and Hox auto-/cross-regulatory interactions.

De novo transcription of multiple Hox cluster genes takes place simultaneously in early Xenopus tropicalis embryos

hox genes are found as clusters in the genome in most bilaterians. The order of genes in the cluster is supposed to be correlated with the site of expression along the anterior-posterior body axis and the timing of expression during development, and these correlations are called spatial and temporal collinearity, respectively. Here we studied the expression dynamics of all hox genes of the diploid species Xenopus tropicalis in four Hox clusters (A–D) by analyzing high-temporal-resolution RNA-seq databases and the results showed that temporal collinearity is not supported, which is consistent with our previous data from allotetraploid Xenopuslaevis. Because the temporal collinearity hypothesis implicitly assumes the collinear order of gene activation, not mRNA accumulation, we determined for the first time the timing of when new transcripts of hox genes are produced, by detecting pre-spliced RNA in whole embryos with reverse transcription and quantitative PCR (RT-qPCR) for all hoxa genes as well as several selected hoxb, hoxc and hoxd genes. Our analyses showed that, coinciding with the RNA-seq results, hoxa genes started to be transcribed in a non-sequential order, and found that multiple genes start expression almost simultaneously or more posterior genes could be expressed earlier than anterior ones. This tendency was also found in hoxb and hoxc genes. These results suggest that temporal collinearity of hox genes is not held during early development of Xenopus.

Summary: qPCR analysis for de novo transcription of hox genes suggest that temporal collinearity is not held for all hox genes during early development of Xenopus tropicalis.

Vertebrate hox temporal collinearity: does it exist and what is it's function?

Hox temporal collinearity (TC) is a mysterious feature of embryogenesis. This article is opportune because of a recent challenge to TC’s existence This challenge is examined and the evidence that TC does exist is presented. Its function is discussed. Temporal collinearity is thought to be important because it lays the basis for Hox spatial collinearity and the vertebrate A-P axial pattern. The time-space translation mechanism whereby this occurs is examined.

The phylum Vertebrata: a case for zoological recognition

The group Vertebrata is currently placed as a subphylum in the phylum Chordata, together with two other subphyla, Cephalochordata (lancelets) and Urochordata (ascidians). The past three decades, have seen extraordinary advances in zoological taxonomy and the time is now ripe for reassessing whether the subphylum position is truly appropriate for vertebrates, particularly in light of recent advances in molecular phylogeny, comparative genomics, and evolutionary developmental biology. Four lines of current research are discussed here. First, molecular phylogeny has demonstrated that Deuterostomia comprises Ambulacraria (Echinodermata and Hemichordata) and Chordata (Cephalochordata, Urochordata, and Vertebrata), each clade being recognized as a mutually comparable phylum. Second, comparative genomic studies show that vertebrates alone have experienced two rounds of whole-genome duplication, which makes the composition of their gene family unique. Third, comparative gene-expression profiling of vertebrate embryos favors an hourglass pattern of development, the most conserved stage of which is recognized as a phylotypic period characterized by the establishment of a body plan definitively associated with a phylum. This mid-embryonic conservation is supported robustly in vertebrates, but only weakly in chordates. Fourth, certain complex patterns of body plan formation (especially of the head, pharynx, and somites) are recognized throughout the vertebrates, but not in any other animal groups. For these reasons, we suggest that it is more appropriate to recognize vertebrates as an independent phylum, not as a subphylum of the phylum Chordata.

A new look at an old question: when did the second whole genome duplication occur in vertebrate evolution?

A recent study used 61 extant animal genomes to reconstruct the chromosomes of the hypothetical amniote ancestor. Comparison of this karyotype to the 17 chordate linkage groups previously inferred in the ancestral chordate indicated that two whole genome duplications probably occurred in the lineage preceding the ancestral vertebrate.

Stepwise participation of HGF/MET signaling in the development of migratory muscle precursors during vertebrate evolution

BACKGROUND

The skeletal musculature of gnathostomes, which is derived from embryonic somites, consists of epaxial and hypaxial portions. Some hypaxial muscles, such as tongue and limb muscles, undergo de-epithelialization and migration during development. Delamination and migration of these myoblasts, or migratory muscle precursors (MMPs), is generally thought to be regulated by hepatocyte growth factor (HGF) and receptor tyrosine kinase (MET) signaling. However, the prevalence of this mechanism and the expression patterns of the genes involved in MMP development across different vertebrate species remain elusive.

RESULTS

We performed a comparative analysis of Hgf and Met gene expression in several vertebrates, including mouse, chicken, dogfish (Scyliorhinus torazame), and lamprey (Lethenteron camtschaticum). While both Hgf and Met were expressed during development in the mouse tongue muscle, and in limb muscle formation in the mouse and chicken, we found no clear evidence for the involvement of HGF/MET signaling in MMP development in shark or lamprey embryos.

CONCLUSIONS

Our results indicate that the expressions and functions of both Hgf and Met genes do not represent shared features of vertebrate MMPs, suggesting a stepwise participation of HGF/MET signaling in MMP development during vertebrate evolution.