Evidence for at least six Hox clusters in the Japanese lamprey (Lethenteron japonicum)

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.

Evolution of the expression and regulation of the nuclear hormone receptor ERR gene family in the chordate lineage

The Estrogen Related Receptor (ERR) nuclear hormone receptor genes have a wide diversity of roles in vertebrate development. In embryos, ERR genes are expressed in several tissues, including the central and peripheral nervous systems. Here we seek to establish the evolutionary history of chordate ERR genes, their expression and their regulation. We examine ERR expression in mollusc, amphioxus and sea squirt embryos, finding the single ERR orthologue is expressed in the nervous system in all three, with muscle expression also found in the two chordates. We show that most jawed vertebrates and lampreys have four ERR paralogues, and that vertebrate ERR genes were ancestrally linked to Estrogen Receptor genes. One of the lamprey paralogues shares conserved expression domains with jawed vertebrate ERRγ in the embryonic vestibuloacoustic ganglion, eye, brain and spinal cord. Hypothesising that conserved expression derives from conserved regulation, we identify a suite of pan-vertebrate conserved non-coding sequences in ERR introns. We use transgenesis in lamprey and chicken embryos to show that these sequences are regulatory and drive reporter gene expression in the nervous system. Our data suggest an ancient association between ERR and the nervous system, including expression in cells associated with photosensation and mechanosensation. This includes the origin in the vertebrate common ancestor of a suite of regulatory elements in the 3' introns that drove nervous system expression and have been conserved from this point onwards.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

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.

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

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

A Revised Perspective on the Evolution of Troponin I and Troponin T Gene Families in Vertebrates

The troponin (Tn) complex, responsible for the Ca2+ activation of striated muscle, is composed of three interacting protein subunits: TnC, TnI, and TnT, encoded by TNNC, TNNI, and TNNT genes. TNNI and TNNT are sister gene families, and in mammals the three TNNI paralogs (TNNI1, TNNI2, TNNI3), which encode proteins with tissue-specific expression, are each in close genomic proximity with one of the three TNNT paralogs (TNNT2, TNNT3, TNNT1, respectively). It has been widely presumed that all vertebrates broadly possess genes of these same three classes, although earlier work has overlooked jawless fishes (cyclostomes) and cartilaginous fishes (chimeras, rays, and sharks), which are distantly related to other jawed vertebrates. With a new phylogenetic and synteny analysis of a diverse array of vertebrates including these taxonomic groups, we define five distinct TNNI classes (TNNI1-5), with TNNI4 and TNNI5 being only present in non-amniote vertebrates and typically found in tandem, and four classes of TNNT (TNNT1-4). These genes are located in four genomic loci that were generated by the 2R whole-genome duplications. TNNI3, encoding "cardiac TnI" in tetrapods, was independently lost in cartilaginous and ray-finned fishes. Instead, ray-finned fishes predominantly express TNNI1 in the heart. TNNI5 is highly expressed in shark hearts and contains a N-terminal extension similar to that of TNNI3 found in tetrapod hearts. Given that TNNI3 and TNNI5 are distantly related, this supports the hypothesis that the N-terminal extension may be an ancestral feature of vertebrate TNNI and not an innovation unique to TNNI3, as has been commonly believed.

Seeking Sense in the Hox Gene Cluster

The Hox gene cluster, responsible for patterning of the head-tail axis, is an ancestral feature of all bilaterally symmetrical animals (the Bilateria) that remains intact in a wide range of species. We can say that the Hox cluster evolved successfully only once since it is commonly the same in all groups, with labial-like genes at one end of the cluster expressed in the anterior embryo, and Abd-B-like genes at the other end of the cluster expressed posteriorly. This review attempts to make sense of the Hox gene cluster and to address the following questions. How did the Hox cluster form in the protostome-deuterostome last common ancestor, and why was this with a particular head-tail polarity? Why is gene clustering usually maintained? Why is there collinearity between the order of genes along the cluster and the positions of their expressions along the embryo? Why do the Hox gene expression domains overlap along the embryo? Why have vertebrates duplicated the Hox cluster? Why do Hox gene knockouts typically result in anterior homeotic transformations? How do animals adapt their Hox clusters to evolve new structural patterns along the head-tail axis?

Structure/Function Studies on the Activation Motif of Two Non-Mammalian Mrap1 Orthologs, and Observations on the Phylogeny of Mrap1, Including a Novel Characterization of an Mrap1 from the Chondrostean Fish, Polyodon spathula

In derived bony vertebrates, activation of the melanocortin-2 receptor (Mc2r) by its ACTH ligand requires chaperoning by the Mc2r accessory protein (Mrap1). The N-terminal domain of the non-mammalian tetrapod MRAP1 from chicken (c; Gallus gallus) has the putative activation motif, W18D19Y20I21, and the N-terminal domain in the neopterygian ray-finned fish Mrap1 from bowfin (bf; Amia calva) has the putative activation motif, Y18D19Y20I21. The current study used an alanine-substitution paradigm to test the hypothesis that only the Y20 position in the Mrap1 ortholog of these non-mammalian vertebrates is required for activation of the respective Mc2r ortholog. Instead, we found that for cMRAP1, single alanine-substitution resulted in a gradient of inhibition of activation (Y20 >> D19 = W18 > I21). For bfMrap1, single alanine-substitution also resulted in a gradient of inhibition of activation (Y20 >> D19 > I21 > Y18). This study also included an analysis of Mc2r activation in an older lineage of ray-finned fish, the paddlefish (p), Polyodon spathula (subclass Chondronstei). Currently no mrap1 gene has been found in the paddlefish genome. When pmc2r was expressed alone in our CHO cell/cAMP reporter gene assay, no activation was observed following stimulation with ACTH. However, when pmc2r was co-expressed with either cmrap1 or bfmrap1 robust dose response curves were generated. These results indicate that the formation of an Mc2r/Mrap1 heterodimer emerged early in the radiation of the bony vertebrates.

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.

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.

Hmx gene conservation identifies the origin of vertebrate cranial ganglia

The evolutionary origin of vertebrates included innovations in sensory processing associated with the acquisition of a predatory lifestyle¹. Vertebrates perceive external stimuli through sensory systems serviced by cranial sensory ganglia, whose neurons arise predominantly from cranial placodes; however, the understanding of the evolutionary origin of placodes and cranial sensory ganglia is hampered by the anatomical differences between living lineages and the difficulty in assigning homology between cell types and structures. Here we show that the homeobox transcription factor Hmx is a constitutive component of vertebrate sensory ganglion development and that in the tunicate Ciona intestinalis, Hmx is necessary and sufficient to drive the differentiation programme of bipolar tail neurons, cells previously thought to be homologues of neural crest^2,3. Using Ciona and lamprey transgenesis, we demonstrate that a unique, tandemly duplicated enhancer pair regulated Hmx expression in the stem-vertebrate lineage. We also show notably robust vertebrate Hmx enhancer function in Ciona, demonstrating that deep conservation of the upstream regulatory network spans the evolutionary origin of vertebrates. These experiments demonstrate regulatory and functional conservation between Ciona and vertebrate Hmx, and point to bipolar tail neurons as homologues of cranial sensory ganglia.

Origin and Evolution of the Multifaceted Adherens Junction Component Plekha7

Plekha7 is a key adherens junction component involved in numerous functions in mammalian cells. Plekha7 is the most studied member of the PLEKHA protein family, which includes eight members with diverse functions. However, the evolutionary history of Plekha7 remains unexplored. Here, we outline the phylogeny and identify the origins of this gene and its paralogs. We show that Plekha7, together with Plekha4, Plekha5, and Plekha6, belong to a subfamily that we name PLEKHA4/5/6/7. This subfamily is distinct from the other Plekha proteins, which form two additional separate subfamilies, namely PLEKHA1/2 and PLEKHA3/8. Sequence, phylogenetic, exon-intron organization, and syntenic analyses reveal that the PLEKHA4/5/6/7 subfamily is represented by a single gene in invertebrates, which remained single in the last common ancestor of all chordates and underwent gene duplications distinctly in jawless and jawed vertebrates. In the latter species, a first round of gene duplications gave rise to the Plekha4/7 and Plekha5/6 pairs and a second round to the four extant members of the subfamily. These observations are consistent with the 1R/2R hypothesis of vertebrate genome evolution. Plekha7 and Plekha5 also exist in two copies in ray-finned fishes, due to the Teleostei-specific whole genome duplication. Similarities between the vertebrate Plekha4/5/6/7 members and non-chordate sequences are restricted to their N-terminal PH domains, whereas similarities across the remaining protein molecule are only sporadically found among few invertebrate species and are limited to the coiled-coil and extreme C-terminal ends. The vertebrate Plekha4/5/6/7 proteins contain extensive intrinsically disordered domains, which are topologically and structurally conserved in all chordates, but not in non-chordate invertebrates. In summary, our study sheds light on the origins and evolution of Plekha7 and the PLEKHA4/5/6/7 subfamily and unveils new critical information suitable for future functional studies of this still understudied group of proteins.

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.

Generation of knock-in lampreys by CRISPR-Cas9-mediated genome engineering

The lamprey represents the oldest group of living vertebrates and has been a key organism in various research fields such as evolutionary developmental biology and neuroscience. However, no knock-in technique for this animal has been established yet, preventing application of advanced genetic techniques. Here, we report efficient generation of F₀ knock-in lampreys by CRISPR-Cas9-mediated genome editing. A donor plasmid containing a heat-shock promoter was co-injected with a short guide RNA (sgRNA) for genome digestion, a sgRNA for donor plasmid digestion, and Cas9 mRNA. Targeting different genetic loci, we succeeded in generating knock-in lampreys expressing photoconvertible protein Dendra2 as well as those expressing EGFP. With its simplicity, design flexibility, and high efficiency, we propose that the present method has great versatility for various experimental uses in lamprey research and that it can also be applied to other “non-model” organisms.

Whole-Genome Duplications and the Diversification of the Globin-X Genes of Vertebrates

Globin-X (GbX) is an enigmatic member of the vertebrate globin gene family with a wide phyletic distribution that spans protostomes and deuterostomes. Unlike canonical globins such as hemoglobins and myoglobins, functional data suggest that GbX does not have a primary respiratory function. Instead, evidence suggests that the monomeric, membrane-bound GbX may play a role in cellular signaling or protection against the oxidation of membrane lipids. Recently released genomes from key vertebrates provide an excellent opportunity to address questions about the early stages of the evolution of GbX in vertebrates. We integrate bioinformatics, synteny, and phylogenetic analyses to characterize the diversity of GbX genes in nonteleost ray-finned fishes, resolve relationships between the GbX genes of cartilaginous fish and bony vertebrates, and demonstrate that the GbX genes of cyclostomes and gnathostomes derive from independent duplications. Our study highlights the role that whole-genome duplications (WGDs) have played in expanding the repertoire of genes in vertebrate genomes. Our results indicate that GbX paralogs have a remarkably high rate of retention following WGDs relative to other globin genes and provide an evolutionary framework for interpreting results of experiments that examine functional properties of GbX and patterns of tissue-specific expression. By identifying GbX paralogs that are products of different WGDs, our results can guide the design of experimental work to explore whether gene duplicates that originate via WGDs have evolved novel functional properties or expression profiles relative to singleton or tandemly duplicated copies of GbX.

Why Do Some Vertebrates Have Microchromosomes?

With more than 70,000 living species, vertebrates have a huge impact on the field of biology and research, including karyotype evolution. One prominent aspect of many vertebrate karyotypes is the enigmatic occurrence of tiny and often cytogenetically indistinguishable microchromosomes, which possess distinctive features compared to macrochromosomes. Why certain vertebrate species carry these microchromosomes in some lineages while others do not, and how they evolve remain open questions. New studies have shown that microchromosomes exhibit certain unique characteristics of genome structure and organization, such as high gene densities, low heterochromatin levels, and high rates of recombination. Our review focuses on recent concepts to expand current knowledge on the dynamic nature of karyotype evolution in vertebrates, raising important questions regarding the evolutionary origins and ramifications of microchromosomes. We introduce the basic karyotypic features to clarify the size, shape, and morphology of macro- and microchromosomes and report their distribution across different lineages. Finally, we characterize the mechanisms of different evolutionary forces underlying the origin and evolution of microchromosomes.

Protein kinase C family evolution in jawed vertebrates

Protein kinase C (PKC) was one of the first kinases identified in human cells. It is now known to constitute a family of kinases that respond to diacylglycerol, phosphatidylserine and for some family members, Ca²⁺. They have a plethora of different functions, such as cell cycle regulation, immune response and memory formation. In mammals, 12 PKC family members have been described, usually divided into 4 different subfamilies. We present here a comprehensive evolutionary analysis of the PKC genes in jawed vertebrates with special focus on the impact of the two tetraploidizations (1R and 2R) before the radiation of jawed vertebrates and the teleost tetraploidization (3R), as illuminated by synteny and paralogon analysis including many neighboring gene families. We conclude that the vertebrate predecessor had five PKC genes, as tunicates and lancelets still do, and that the PKC family should therefore ideally be organized into five subfamilies. The 1R and 2R events led to a total of 12 genes distributed among these five subfamilies. All 12 genes are still present in some of the major lineages of jawed vertebrates, including mammals, whereas birds and cartilaginous fishes have lost one member. The 3R event added another nine genes in teleosts, bringing the total to 21 genes. The zebrafish, a common experimental model animal, has retained 19. We have found no independent gene duplications. Thus, the genome doublings completely account for the complexity of this gene family in jawed vertebrates and have thereby had a huge impact on their evolution.

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.

The Japanese lamprey (Lethenteron camtschaticum) expresses functional lysophosphatidic acid receptors

Lysophosphatidic acid (LPA) signaling plays diverse roles in the development of various vertebrates such as mammals and fish. The lamprey is a fish that retains ancestral features of vertebrates, but information regarding lamprey LPA receptor genes is limited. Here, using information from the lamprey genome database, we cloned two LPA receptor genes, Lpar1 and Lpar5, from the Japanese lamprey (Lethenteron camtschaticum). Lamprey Lpar1 had a high amino acid identity to mouse and medaka fish Lpar1, whereas Lpar5 amino acid sequences were more diverse between species. Our functional analyses using a heterologous expression system demonstrated that Lpar1 and Lpar5 responded to LPA treatment with G_12/13-associated cellular responses, which are indicative of cytoskeletal actions. The existence of functional LPA receptors in the Japanese lamprey suggests that LPA receptor-dependent signals contribute to lamprey growth and development.

Evolutionary history and functional characterization of Lj-TICAM-a and Lj-TICAM-b formed via lineage-specific tandem duplication in lamprey (Lampetra japonica)

Toll/interleukin-1 receptor domain-containing adaptor molecule (TICAM) genes respond to infections. We identified TICAM-a and TICAM-b in Lampetra japonica and investigated their evolutionary history and potential function via comparative genomics and molecular evolution analyses. They are arranged in tandem and evolved from a multi-exon to a single-exon structure. Lj-TICAM-a and Lj-TICAM-b might be the ancestral gene of the vertebrate TICAM genes. Lj-TICAM-b arose via a lamprey-specific tandem duplication event. Both genes are expressed in many tissues during an immune response, and exhibit different responses to peptidoglycan, indicating their functional divergence. Simultaneous overexpression of both proteins activated nuclear factor κB expression and co-immunoprecipitation assays indicated that they might form a complex for signal transduction. However, unlike in mammals, the TICAM-dependent signaling pathway in lamprey might rely on TRAF3 rather than on TRAF6. These results suggest that both Lj-TICAM-a and Lj-TICAM-b play a role in host defenses.

Functional Domains and Evolutionary History of the PMEL and GPNMB Family Proteins

The ancient paralogs premelanosome protein (PMEL) and glycoprotein nonmetastatic melanoma protein B (GPNMB) have independently emerged as intriguing disease loci in recent years. Both proteins possess common functional domains and variants that cause a shared spectrum of overlapping phenotypes and disease associations: melanin-based pigmentation, cancer, neurodegenerative disease and glaucoma. Surprisingly, these proteins have yet to be shown to physically or genetically interact within the same cellular pathway. This juxtaposition inspired us to compare and contrast this family across a breadth of species to better understand the divergent evolutionary trajectories of two related, but distinct, genes. In this study, we investigated the evolutionary history of PMEL and GPNMB in clade-representative species and identified TMEM130 as the most ancient paralog of the family. By curating the functional domains in each paralog, we identified many commonalities dating back to the emergence of the gene family in basal metazoans. PMEL and GPNMB have gained functional domains since their divergence from TMEM130, including the core amyloid fragment (CAF) that is critical for the amyloid potential of PMEL. Additionally, the PMEL gene has acquired the enigmatic repeat domain (RPT), composed of a variable number of imperfect tandem repeats; this domain acts in an accessory role to control amyloid formation. Our analyses revealed the vast variability in sequence, length and repeat number in homologous RPT domains between craniates, even within the same taxonomic class. We hope that these analyses inspire further investigation into a gene family that is remarkable from the evolutionary, pathological and cell biology perspectives.

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.

3D genomics across the tree of life reveals condensin II as a determinant of architecture type

We investigated genome folding across the eukaryotic tree of life. We find two types of three-dimensional (3D) genome architectures at the chromosome scale. Each type appears and disappears repeatedly during eukaryotic evolution. The type of genome architecture that an organism exhibits correlates with the absence of condensin II subunits. Moreover, condensin II depletion converts the architecture of the human genome to a state resembling that seen in organisms such as fungi or mosquitoes. In this state, centromeres cluster together at nucleoli, and heterochromatin domains merge. We propose a physical model in which lengthwise compaction of chromosomes by condensin II during mitosis determines chromosome-scale genome architecture, with effects that are retained during the subsequent interphase. This mechanism likely has been conserved since the last common ancestor of all eukaryotes.

Functional and Evolutionary Diversification of Otx2 and Crx in Vertebrate Retinal Photoreceptor and Bipolar Cell Development

Otx family homeoproteins Otx2 and Crx are expressed in photoreceptor precursor cells and bind to the common DNA-binding consensus sequence, but these two proteins have distinct functions in retinal development. To examine the functional substitutability of Otx2 and Crx, we generate knockin mouse lines in which Crx is replaced by Otx2 and vice versa. We find that Otx2 and Crx cannot be substituted in photoreceptor development. Subsequently, we investigate the function of Otx2 in photoreceptor and bipolar cell development. High Otx2 levels induce photoreceptor cell fate but not bipolar cell fate, whereas reduced Otx2 expression impairs bipolar cell maturation and survival. Furthermore, we identify Otx2 and Crx in the lamprey genome by using synteny analysis, suggesting that the last common ancestor of vertebrates possesses both Otx2 and Crx. We find that the retinal Otx2 expression pattern is different between lampreys and mice, suggesting that neofunctionalization of Otx2 occurred in the jawed vertebrate lineage.

Differential expression of somatostatin genes in the central nervous system of the sea lamprey

The identification of three somatostatin (SST) genes (SSTa, SSTb, and SSTc) in lampreys (Tostivint et al. Gen Comp Endocrinol 237:89-97 https://doi.org/10.1016/j.ygcen.2016.08.006 , 2016) prompted us to study their expression in the brain and spinal cord of the sea lamprey by in situ hybridization. These three genes were only expressed in equivalent neuronal populations in the hypothalamus. In other regions, SST transcripts showed clear differential expression. In the telencephalon, SSTc-positive cells were observed in the medial pallium, ventral part of the lateral pallium, striatum, subhippocampal lobe, and preoptic region. In the diencephalon, SSTa-positive cells were observed in the thalamus and SSTc-positive cells in the prethalamus, posterior tubercle, pretectal area, and nucleus of the medial longitudinal fascicle. In the midbrain, SSTc-positive cells were observed in the torus semicircularis, lateral reticular area, and perioculomotor tegmentum. Different SSTa- and SSTc-positive populations were observed in the isthmus. SSTc neurons were also observed in the rostral octavolateralis area and caudal rhombencephalon. In the spinal cord, SSTa was expressed in cerebrospinal-fluid-contacting (CSF-c) neurons and SSTc in non-CSF-c interneurons. Comparison with previous immunohistochemical studies using anti-SST-14 antibodies strongly suggests that SST-14-like neurons correspond with the SSTa populations. Thus, the SSTc populations were not reported previously in immunohistochemical studies. Cluster-based analyses and alignments of mature peptides suggested that SSTa is an ortholog of SST1 and that SSTb is closely related to SST2 and SST6. These results provide important new insights into the evolution of the somatostatinergic system in vertebrates.

Connexins during 500 Million Years-From Cyclostomes to Mammals

It was previously shown that the connexin gene family had relatively similar subfamily structures in several vertebrate groups. Still, many details were left unclear. There are essentially no data between tunicates, which have connexins that cannot be divided into the classic subfamilies, and teleosts, where the subfamilies are easily recognized. There are also relatively few data for the groups that diverged between the teleosts and mammals. As many of the previously analyzed genomes have been improved, and many more genomes are available, we reanalyzed the connexin gene family and included species from all major vertebrate groups. The major results can be summarized as follows: (i) The same connexin subfamily structures are found in all Gnathostomata (jawed vertebrates), with some variations due to genome duplications, gene duplications and gene losses. (ii) In contrast to previous findings, birds do not have a lower number of connexins than other tetrapods. (iii) The cyclostomes (lampreys and hagfishes) possess genes in the alpha, beta, gamma and delta subfamilies, but only some of the genes show a phylogenetic affinity to specific genes in jawed vertebrates. Thus, two major evolutionary transformations have occurred in this gene family, from tunicates to cyclostomes and from cyclostomes to jawed vertebrates.

Comparative Genomics within and across Bilaterians Illuminates the Evolutionary History of ALK and LTK Proto-Oncogene Origination and Diversification

Comparative genomic analyses have enormous potential for identifying key genes central to human health phenotypes, including those that promote cancers. In particular, the successful development of novel therapeutics using model species requires phylogenetic analyses to determine molecular homology. Accordingly, we investigate the evolutionary histories of anaplastic lymphoma kinase (ALK)—which can underlie tumorigenesis in neuroblastoma, nonsmall cell lung cancer, and anaplastic large-cell lymphoma—its close relative leukocyte tyrosine kinase (LTK) and their candidate ligands. Homology of ligands identified in model organisms to those functioning in humans remains unclear. Therefore, we searched for homologs of the human genes across metazoan genomes, finding that the candidate ligands Jeb and Hen-1 were restricted to nonvertebrate species. In contrast, the ligand augmentor (AUG) was only identified in vertebrates. We found two ALK-like and four AUG-like protein-coding genes in lamprey. Of these six genes, only one ALK-like and two AUG-like genes exhibited early embryonic expression that parallels model mammal systems. Two copies of AUG are present in nearly all jawed vertebrates. Our phylogenetic analysis strongly supports the presence of previously unrecognized functional convergences of ALK and LTK between actinopterygians and sarcopterygians—despite contemporaneous, highly conserved synteny of ALK and LTK. These findings provide critical guidance regarding the propriety of fish and mammal models with regard to model organism-based investigation of these medically important genes. In sum, our results provide the phylogenetic context necessary for effective investigations of the functional roles and biology of these critically important receptors.

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.

VLRs expression were significantly affected by complement C3 knockdown morphants in Lampetra morii

The complement component 3 of the lamprey, a jawless vertebrate, functions as an opsonin during the phagocytosis of rabbit red cells. Furthermore, lamprey C3 may be activated and cleaved into C3b, which is attached to the surface of target cells in the cytolytic process. However, the mechanism mediating the biological function of C3 in the lamprey is unknown. To our knowledge, this study is the first to show that variable lymphocyte receptors (VLRs) expression were significantly affected by complement C3 knockdown morphants in Lampetra morii. We identified the C3 gene in the lamprey genome based on its orthologs, conserved synteny, functional domains, phylogenetic tree, and conserved motifs. Additionally, we determined the optimal infection concentration of Aeromonas hydrophila to perform immune stimulation experiments in the lamprey larvae. The quantitative real-time polymerase chain reaction and immunofluorescence analyses revealed that the expression of Lampetra morii C3 (lmC3) was significantly upregulated in the larvae infected with 10⁷ CFU/mL of A. hydrophila. The lmC3 morphants (lmC3 MO) of lamprey larvae were generated by morpholino-mediated knockdown. The lmC3 MO larvae were highly susceptible to A. hydrophila infection, which indicated that lmC3 is critical in lamprey immune response. The expression of a selected panel of orthologous genes was comparatively analyzed in the infected wild type, infected lmC3 MO, infected control MO, uninfected wild type and uninfected lmC3 MO one-month-old ammocoete larvae. The knockdown of lmC3 strongly affected the expression of VLRA⁺/VLRB⁺/VLRC⁺-associated genes, which was also confirmed by immunohistochemical analysis. Thus, VLR expression were significantly affected by complement C3 knockdown morphants in Lampetra morii.

Complement component C1q plays a critical role in VLRA/VLRC-mediated immune response

In jawless vertebrates, the lamprey complement component C1q (LC1q) acts as a lectin and activates lamprey complement component C3 (LC3) in association with mannose-binding lectin (MBL)-associated serine protease (MASP) via the lectin pathway. Furthermore, LC1q may interact with variable lymphocyte receptor B (VLRB) in a complex with antigens and mediate the activation of LC3, leading to cytolysis. In the present study, we found, for the first time, that LC1q plays a critical role in VLRA/VLRC-mediated immune response. Escherichia coli, Shigella flexneri, Aeromonas hydrophila, Pseudomonas plecoglossicida, Aeromonas allosaccharophila, P. luteola, Brevundimonas diminuta, and Bacillus cereus were isolated from infected Lampetra morii in our laboratory and identified using the 16s rRNA method. A. hydrophila was confirmed as a rapidly spreading lethal pathogen in the larvae of L. morii and was used in subsequent immune stimulation experiments. The results of real-time quantitative polymerase chain reaction (Q-PCR) and immunofluorescence analyses indicated that the RNA and protein expression levels of LC1q were upregulated following exposure to 10⁷ cfu/mL of A. hydrophila, compared to the levels of the naïve group. We obtained LC1q morphants (LC1q MO) of lamprey larvae by morpholino-mediated knockdowns. We found that LC1q played key roles in the embryonic development of lamprey. The median lethal time (LT50) of LC1q MO larvae was 2 d after being exposed to the pathogens, whereas the LT50 of control MO was 5 d. The drastic decrease in LT50 values after LC1q knockdown implies that LC1q plays a critical role in lamprey immune response. Gene expression profiles of LC1q-deficient A. hydrophila, control MO A. hydrophila, wild type A. hydrophila, and naive 1-month-old ammocoetes larvae were compared by examining the expression levels of a selected panel of orthologous genes. It is worth mentioning that LC1q MO affected the VLRA+/VLRC + population genes but did not affect the VLRB + populations. Immunohistochemical data indicated that LC1q deficiency also affected VLRA and VLRC but not VLRB. Thus, LC1q plays a critical role in VLRA/VLRC-mediated immune response in lamprey.

Lamprey VLRB participates in pathogen detection, VLRB/L-BLNK/L-NF-κB (B-like cells) signal transduction, and development

In jawed vertebrates, B cell receptors (BCR) are primary pathogen detectors that activate downstream signaling pathways to express adaptive immune effectors. In jawless vertebrates, the variable lymphocyte receptors (VLR) B positive lymphocytes can express and secrete specific VLRB molecules in an analogous manner to that of immunoglobulins by B cells in jawed vertebrates. Our study is the first to demonstrate the possibility of incubation of fertilized eggs and artificial breeding of Lampetra morii larvae throughout their life cycle under laboratory condition. We also found that VLRB, lamprey B-cell linker (L-BLNK), and lamprey nuclear factor-kappa B (L-NF-κB) play key roles in early larval development. Aeromonas hydrophila was found to be a lethal pathogen of L. morii larvae causing rapid infection at a concentration of 10⁷ cfu/mL qRT-PCR results revealed that gene expression levels of VLRB, L-BLNK, and L-NF-κB were up-regulated significantly. Ten-day infection trials showed that VLRB, L-BLNK, and L-NF-κB are crucial for lamprey immune response. Furthermore, the expression levels of L-BLNK and L-NF-κB were down-regulated drastically both at mRNA and protein levels after bacterial infection than in the naive group of VLRB morphants. A similar expression pattern of VLRB and L-BLNK was found in L-NF-κB morphants post bacterial infection. The results were strikingly different in the other two morphants. The VLRB and L-NF-κB expression levels were found to be down-regulated at mRNA and protein levels by less than 30% and 45%, respectively, in L-BLNK morphants compared to those in the naive group. These results indicate that L-BLNK and L-NF-κB might participate in VLRB-mediated immune response. Additionally, in VLRB morphants, the mRNA expression levels of some genes, especially the ones expressed in VLRB⁺ lymphocytes but not in VLRA⁺ lymphocytes, were found to be affected. Therefore, these findings of B-like lymphocytes in lamprey offer key evidence with regard to the evolution of adaptive immunity.

Discovery of four Noggin genes in lampreys suggests two rounds of ancient genome duplication

The secreted protein Noggin1 was the first discovered natural embryonic inducer produced by cells of the Spemann organizer. Thereafter, it was shown that vertebrates have a whole family of Noggin genes with different expression patterns and functional properties. For example, Noggin1 and Noggin2 inhibit the activity of BMP, Nodal/Activin and Wnt-beta-catenin signalling, while Noggin4 cannot suppress BMP but specifically modulates Wnt signalling. In this work, we described and investigated phylogeny and expression patterns of four Noggin genes in lampreys, which represent the most basally divergent group of extant vertebrates, the cyclostomes, belonging to the superclass Agnatha. Assuming that lampreys have Noggin homologues in all representatives of another superclass of vertebrates, the Gnathostomata, we propose a model for Noggin family evolution in vertebrates. This model is in agreement with the hypotheses suggesting two rounds of genome duplication in the ancestor of vertebrates before the divergence of Agnatha and Gnathostomata.

Ermakova et al. report four Noggin genes in lampreys and using phylogenetics, gene synteny analysis, and in situ hybridization, suggest that the Noggin gene underwent two rounds of duplication and evolved specific functions before the divergence of vertebrate and lamprey lineages. These findings offer insight into early vertebrate genome and developmental evolution.

Novel developmental bases for the evolution of hypobranchial muscles in vertebrates

BACKGROUND

Vertebrates are characterized by possession of hypobranchial muscles (HBMs). Cyclostomes, or modern jawless vertebrates, possess a rudimentary and superficial HBM lateral to the pharynx, whereas the HBM in jawed vertebrates is internalized and anteroposteriorly specified. Precursor cells of the HBM, marked by expression of Lbx1, originate from somites and undergo extensive migration before becoming innervated by the hypoglossal nerve. How the complex form of HBM arose in evolution is relevant to the establishment of the vertebrate body plan, but despite having long been assumed to be similar to that of limb muscles, modification of developmental mechanisms of HBM remains enigmatic.

RESULTS

Here we characterize the expression of Lbx genes in lamprey and hagfish (cyclostomes) and catshark (gnathostome; jawed vertebrates). We show that the expression patterns of the single cyclostome Lbx homologue, Lbx-A, do not resemble the somitic expression of mammalian Lbx1. Disruption of Lbx-A revealed that LjLbx-A is required for the formation of both HBM and body wall muscles, likely due to the insufficient extension of precursor cells rather than to hindered muscle differentiation. Both homologues of Lbx in the catshark were expressed in the somitic muscle primordia, unlike in amniotes. During catshark embryogenesis, Lbx2 is expressed in the caudal HBM as well as in the abdominal rectus muscle, similar to lamprey Lbx-A, whereas Lbx1 marks the rostral HBM and pectoral fin muscle.

CONCLUSIONS

We conclude that the vertebrate HBM primarily emerged as a specialized somatic muscle to cover the pharynx, and the anterior internalized HBM of the gnathostomes is likely a novelty added rostral to the cyclostome-like HBM, for which duplication and functionalization of Lbx genes would have been a prerequisite.

Two highly related odorant receptors specifically detect α-bile acid pheromones in sea lamprey (Petromyzon marinus)

Pheromones play critical roles in habitat identification and reproductive behavior synchronization in the sea lamprey (Petromyzon marinus). The bile acid 3-keto petromyzonol sulfate (3kPZS) is a major component of the sex pheromone mixture from male sea lamprey that induces specific olfactory and behavioral responses in conspecific individuals. Olfactory receptors interact directly with pheromones, which is the first step in their detection, but identifying the cognate receptors of specific pheromones is often challenging. Here, we deorphanized two highly related odorant receptors (ORs), OR320a and OR320b, of P. marinus that respond to 3kPZS. In a heterologous expression system coupled to a cAMP-responsive CRE-luciferase, OR320a and OR320b specifically responded to C24 5α-bile acids, and both receptors were activated by the same set of 3kPZS analogs. OR320a displayed larger responses to all 3kPZS analogs than did OR320b. This difference appeared to be largely determined by a single amino acid residue, Cys-792.56, the C-terminal sixth residue relative to the most conserved residue in the second transmembrane domain (2.56) of OR320a. This region of TM2 residues 2.56-2.60 apparently is critical for the detection of steroid compounds by odorant receptors in lamprey, zebrafish, and humans. Finally, we identified OR320 orthologs in Japanese lamprey (Lethenteron camtschaticum), suggesting that the OR320 family may be widely present in lamprey species and that OR320 may be under purifying selection. Our results provide a system to examine the origin of olfactory steroid detection in vertebrates and to define a highly conserved molecular mechanism for steroid-ligand detection by G protein-coupled receptors.

Corticotropin-Releasing Hormone (CRH) Gene Family Duplications in Lampreys Correlate With Two Early Vertebrate Genome Doublings

The ancestor of gnathostomes (jawed vertebrates) is generally considered to have undergone two rounds of whole genome duplication (WGD). The timing of these WGD events relative to the divergence of the closest relatives of the gnathostomes, the cyclostomes, has remained contentious. Lampreys and hagfishes are extant cyclostomes whose gene families can shed light on the relationship between the WGDs and the cyclostome-gnathostome divergence. Previously, we have characterized in detail the evolution of the gnathostome corticotropin-releasing hormone (CRH) family and found that its five members arose from two ancestral genes that existed before the WGDs. The two WGDs resulted, after secondary losses, in one triplet consisting of CRH1, CRH2, and UCN1, and one pair consisting of UCN2 and UCN3. All five genes exist in representatives for cartilaginous fishes, ray-finned fishes, and lobe-finned fishes. Differential losses have occurred in some lineages. We present here analyses of CRH-family members in lamprey and hagfish by comparing sequences and gene synteny with gnathostomes. We found five CRH-family genes in each of two lamprey species (Petromyzon marinus and Lethenteron camtschaticum) and two genes in a hagfish (Eptatretus burgeri). Synteny analyses show that all five lamprey CRH-family genes have similar chromosomal neighbors as the gnathostome genes. The most parsimonious explanation is that the lamprey CRH-family genes are orthologs of the five gnathostome genes and thus arose in the same chromosome duplications. This suggests that lampreys and gnathostomes share the same two WGD events and that these took place before the lamprey-gnathostome divergence.

Reconstruction of the Carbohydrate 6-O Sulfotransferase Gene Family Evolution in Vertebrates Reveals Novel Member, CHST16, Lost in Amniotes

Glycosaminoglycans are sulfated polysaccharide molecules, essential for many biological processes. The 6-O sulfation of glycosaminoglycans is carried out by carbohydrate 6-O sulfotransferases (C6OSTs), previously named Gal/GalNAc/GlcNAc 6-O sulfotransferases. Here, for the first time, we present a detailed phylogenetic reconstruction, analysis of gene synteny conservation and propose an evolutionary scenario for the C6OST family in major vertebrate groups, including mammals, birds, nonavian reptiles, amphibians, lobe-finned fishes, ray-finned fishes, cartilaginous fishes, and jawless vertebrates. The C6OST gene expansion likely started early in the chordate lineage, giving rise to four ancestral genes after the divergence of tunicates and before the emergence of extant vertebrates. The two rounds of whole-genome duplication in early vertebrate evolution (1R/2R) only contributed two additional C6OST subtype genes, increasing the vertebrate repertoire from four genes to six, divided into two branches. The first branch includes CHST1 and CHST3 as well as a previously unrecognized subtype, CHST16 that was lost in amniotes. The second branch includes CHST2, CHST7, and CHST5. Subsequently, local duplications of CHST5 gave rise to CHST4 in the ancestor of tetrapods, and to CHST6 in the ancestor of primates. The teleost-specific gene duplicates were identified for CHST1, CHST2, and CHST3 and are result of whole-genome duplication (3R) in the teleost lineage. We could also detect multiple, more recent lineage-specific duplicates. Thus, the vertebrate repertoire of C6OST genes has been shaped by gene duplications and gene losses at several stages of vertebrate evolution, with implications for the evolution of skeleton, nervous system, and cell-cell interactions.

The non-visual opsins expressed in deep brain neurons projecting to the retina in lampreys

In lower vertebrates, brain photoreceptor cells express vertebrate-specific non-visual opsins. We previously revealed that a pineal-related organ-specific opsin, parapinopsin, is UV-sensitive and allows pineal wavelength discrimination in lampreys and teleost. The Australian pouched lamprey was recently reported as having two parapinopsin-related genes. We demonstrate that a parapinopsin-like opsin from the Japanese river lamprey exhibits different molecular properties and distribution than parapinopsin. This opsin activates Gi-type G protein in a mammalian cell culture assay in a light-dependent manner. Heterologous action spectroscopy revealed that the opsin forms a violet to blue-sensitive pigment. Interestingly, the opsin is co-localised with green-sensitive P-opsin in the cells of the M5 nucleus of Schober (M5NS) in the mesencephalon of the river and brook lamprey. Some opsins-containing cells of the river lamprey have cilia and others an axon projecting to the retina. The opsins of the brook lamprey are co-localised in the cilia of centrifugal neurons projecting to the retina, suggesting that cells expressing the parapinopsin-like opsin and P-opsin are sensitive to violet to green light. Moreover, we found neural connections between M5NS cells expressing the opsins and the retina. These findings suggest that the retinal activity might be modulated by brain photoreception.

Deeply conserved synteny resolves early events in vertebrate evolution

Although it is widely believed that early vertebrate evolution was shaped by ancient whole-genome duplications, the number, timing and mechanism of these events remain elusive. Here, we infer the history of vertebrates through genomic comparisons with a new chromosome-scale sequence of the invertebrate chordate amphioxus. We show how the karyotypes of amphioxus and diverse vertebrates are derived from 17 ancestral chordate linkage groups (and 19 ancestral bilaterian groups) by fusion, rearrangement and duplication. We resolve two distinct ancient duplications based on patterns of chromosomal conserved synteny. All extant vertebrates share the first duplication, which occurred in the mid/late Cambrian by autotetraploidization (that is, direct genome doubling). In contrast, the second duplication is found only in jawed vertebrates and occurred in the mid-late Ordovician by allotetraploidization (that is, genome duplication following interspecific hybridization) from two now-extinct progenitors. This complex genomic history parallels the diversification of vertebrate lineages in the fossil record.

Constitutive activity of a spermine receptor is maintained by a single site in the C-terminal

Olfactory receptors are G-protein coupled receptors (GPCRs) that enable olfactory epithelia to detect odorants. These GPCRs may also show constitutive activity, which play important roles in the development and responses of odorant receptor neurons. However, little is known about the molecular characteristics that support the constitutive activities in olfactory receptors. Here, we characterize a pair of olfactory receptor orthologs that show similar ligand-dependent activity but different levels of constitutive activity, and elucidate the molecular characteristics that maintain the constitutive activity. Previously, PmTAAR348, a sea lamprey (Petromyzon marinus) olfactory receptor that is activated by the male sex pheromone spermine has been reported. In this study, we identified LmTAAR348 of Northeast Chinese lamprey (Lethenteron morii) as an ortholog of PmTAAR348. When expressed in HEK293T cell lines, both receptors showed similar levels of activation when exposed to spermine. However, the constitutive activity of LmTAAR348 was 100-fold higher than that of PmTAAR348. Using site-directed mutagenesis, we screened all 13 amino acid residues (aa) that differed between the two orthologs and found that a switch in position 340 reversed the constitutive activity levels between LmTAAR348 and PmTAAR348. Mutating the remaining 12 aa did not affect the ligand-dependent or constitutive activation. Moreover, both the ligand-dependent and constitutive activation of TAAR348 are G_olf (G protein ⍺ subunit olfactory type) independent. We conclude that a single aa in the C-terminal maintains the constitutive activity in a spermine receptor.

Origin and evolution of the Rax homeobox gene by comprehensive evolutionary analysis

Rax is one of the key transcription factors crucial for vertebrate eye development. In this study, we conducted comprehensive evolutionary analysis of Rax. We found that Bilateria and Cnidaria possess Rax, but Placozoa, Porifera, and Ctenophora do not, implying that the origin of the Rax gene dates back to the common ancestor of Cnidaria and Bilateria. The results of molecular phylogenetic and synteny analyses on Rax loci between jawed and jawless vertebrates indicate that segmental duplication of the Rax locus occurred in an early common ancestor of jawed vertebrates, resulting in two Rax paralogs in jawed vertebrates, Rax and Rax2. By analyzing 86 mammalian genomes from all four major groups of mammals, we found that at least five independent Rax2 gene loss events occurred in mammals. This study may provide novel insights into the evolution of the eye.

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.

Lampreys, the jawless vertebrates, contain three Pax6 genes with distinct expression in eye, brain and pancreas

The transcription factor Pax6 is crucial for the development of the central nervous system, eye, olfactory system and pancreas, and is implicated in human disease. While a single Pax6 gene exists in human and chicken, Pax6 occurs as a gene family in other vertebrates, with two members in elephant shark, Xenopus tropicalis and Anolis lizard and three members in teleost fish such as stickleback and medaka. However, the complement of Pax6 genes in jawless vertebrates (cyclostomes), the sister group of jawed vertebrates (gnathostomes), is unknown. Using a combination of BAC sequencing and genome analysis, we discovered three Pax6 genes in lampreys. Unlike the paired-less Pax6 present in some gnathostomes, all three lamprey Pax6 have a highly conserved full-length paired domain. All three Pax6 genes are expressed in the eye and brain, with variable expression in other tissues. Notably, lamprey Pax6α transcripts are found in the pancreas, a vertebrate-specific organ, indicating the involvement of Pax6 in development of the pancreas in the vertebrate ancestor. Multi-species sequence comparisons revealed only a single conserved non-coding element, in the lamprey Pax6β locus, with similarity to the PAX6 neuroretina enhancer. Using a transgenic zebrafish enhancer assay we demonstrate functional conservation of this element over 500 million years of vertebrate evolution.

Functional characterization and osmoregulatory role of the Na+-K+-2Cl- cotransporter in the gill of sea lamprey (Petromyzon marinus), a basal vertebrate

The present study provides molecular and functional characterization of Na⁺-K⁺-2Cl^- cotransporter (NKCC1/Slc12a2) in the gills of sea lamprey (Petromyzon marinus), the most basal extant vertebrate with an osmoregulatory strategy. We report the full-length peptide sequence for the lamprey Na-K-Cl cotransporter 1 (NKCC1), which we show groups strongly with and occupies a basal position among other vertebrate NKCC1 sequences. In postmetamorphic juvenile lamprey, nkcc1 mRNA was present in many tissues but was fivefold higher in the gill than any other examined tissue, and NKCC1 protein was only detected in the gill. Gill mRNA and protein abundances of NKCC1 and Na⁺-K⁺-ATPase (NKA/Atp1a1) were significantly upregulated (20- to 200-fold) during late metamorphosis in fresh water, coinciding with the development of salinity tolerance, and were upregulated an additional twofold after acclimation to seawater (SW). Immunohistochemistry revealed that NKCC1 in the gill is found in filamental ionocytes coexpressing NKA, which develop during metamorphosis in preparation for SW entry. Lamprey treated with bumetanide, a widely used pharmacological inhibitor of NKCC1, exhibited higher plasma Cl^- and osmolality as well as reduced muscle water content after 24 h in SW; there were no effects of bumetanide in freshwater-acclimated lamprey. This work provides the first functional characterization of NKCC1 as a mechanism for branchial salt secretion in lampreys, providing evidence that this mode of Cl^- secretion has been present among vertebrates for ~550 million years.

Characterization and Evolution of Germ1, an Element that Undergoes Diminution in Lampreys (Cyclostomata: Petromyzontidae)

The sea lamprey (Petromyzon marinus) undergoes substantial genomic alterations during embryogenesis in which specific sequences are deleted from the genome of somatic cells yet retained in cells of the germ line. One element that undergoes diminution in P. marinus is Germ1, which consists of a somatically rare (SR) region and a fragment of 28S rDNA. Although the SR-region has been used as a marker for genomic alterations in lampreys, the evolutionary significance of its diminution is unknown. We examined the Germ1 element in five additional species of lamprey to better understand its evolutionary significance. Each representative species contained sequences similar enough to the Germ1 element of P. marinus to be detected via PCR and Southern hybridizations, although the SR-regions of Lampetra aepyptera and Lethenteron appendix are quite divergent from the homologous sequences of Petromyzon and three species of Ichthyomyzon. Lamprey Germ1 sequences have a number of features characteristic of the R2 retrotransposon, a mobile element that specifically targets 28S rDNA. Phylogenetic analyses of the SR-regions revealed patterns generally consistent with relationships among the species included in our study, although the 28S-fragments of each species/genus were most closely related to its own functional rDNA, suggesting that the two components of Germ1 were assembled independently in each lineage. Southern hybridizations showed evidence of genomic alterations involving Germ1 in each species. Our results suggest that Germ1 is a R2 retroelement that occurs in the genome of P. marinus and other petromyzontid lampreys, and that its diminution is incidental to the reduction in rDNA copies during embryogenesis.

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.