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

Comparative genomics illuminates karyotype and sex chromosome evolution of sharks

Chondrichthyes is an important lineage to reconstruct the evolutionary history of vertebrates. Here, we analyzed genome synteny for six chondrichthyan chromosome-level genomes. Our comparative analysis reveals a slow evolutionary rate of chromosomal changes, with infrequent but independent fusions observed in sharks, skates, and chimaeras. The chondrichthyan common ancestor had a proto-vertebrate-like karyotype, including the presence of 18 microchromosome pairs. The X chromosome is a conversed microchromosome shared by all sharks, suggesting a likely common origin of the sex chromosome at least 181 million years ago. We characterized the Y chromosomes of two sharks that are highly differentiated from the X except for a small young evolutionary stratum and a small pseudoautosomal region. We found that shark sex chromosomes lack global dosage compensation but that dosage-sensitive genes are locally compensated. Our study on shark chromosome evolution enhances our understanding of shark sex chromosomes and vertebrate chromosome evolution.

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.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

Whole-genome Sequencing Reveals Autooctoploidy in Chinese Sturgeon and Its Evolutionary Trajectories

The order Acipenseriformes, which includes sturgeons and paddlefishes, represents "living fossils" with complex genomes that are good models for understanding whole-genome duplication (WGD) and ploidy evolution in fishes. Here, we sequenced and assembled the first high-quality chromosome-level genome for the complex octoploid Acipenser sinensis (Chinese sturgeon), a critically endangered species that also represents a poorly understood ploidy group in Acipenseriformes. Our results show that A. sinensis is a complex autooctoploid species containing four kinds of octovalents (8n), a hexavalent (6n), two tetravalents (4n), and a divalent (2n). An analysis taking into account delayed rediploidization reveals that the octoploid genome composition of Chinese sturgeon results from two rounds of homologous WGDs, and further provides insights into the timing of its ploidy evolution. This study provides the first octoploid genome resource of Acipenseriformes for understanding ploidy compositions and evolutionary trajectories of polyploid fishes.

Shark genome size evolution and its relationship with cellular, life-history, ecological, and diversity traits

Among vertebrates, sharks exhibit both large and heterogeneous genome sizes ranging from 2.86 to 17.05 pg. Aiming for a better understanding of the patterns and causalities of shark genome size evolution, we applied phylogenetic comparative methods to published genome-size estimates for 71 species representing the main phylogenetic lineages, life-histories and ecological traits. The sixfold range of genome size variation was strongly traceable throughout the phylogeny, with a major expansion preceding shark diversification during the late Paleozoic and an ancestral state (6.33 pg) close to the present-day average (6.72 pg). Subsequent deviations from this average occurred at higher rates in squalomorph than in galeomorph sharks and were unconnected to evolutionary changes in the karyotype architecture, which were dominated by descending disploidy events. Genome size was positively correlated with cell and nucleus sizes and negatively with metabolic rate. The metabolic constraints on increasing genome size also manifested at higher phenotypic scales, with large genomes associated with slow lifestyles and purely marine waters. Moreover, large genome sizes were also linked to non-placental reproductive modes, which may entail metabolically less demanding embryological developments. Contrary to ray-finned fishes, large genome size was associated neither with the taxonomic diversity of affected clades nor with low genetic diversity.

Extensive DNA methylome rearrangement during early lamprey embryogenesis

DNA methylation (5mC) is a repressive gene regulatory mark widespread in vertebrate genomes, yet the developmental dynamics in which 5mC patterns are established vary across species. While mammals undergo two rounds of global 5mC erasure, teleosts, for example, exhibit localized maternal-to-paternal 5mC remodeling. Here, we studied 5mC dynamics during the embryonic development of sea lamprey, a jawless vertebrate which occupies a critical phylogenetic position as the sister group of the jawed vertebrates. We employed 5mC quantification in lamprey embryos and tissues, and discovered large-scale maternal-to-paternal epigenome remodeling that affects ~30% of the embryonic genome and is predominantly associated with partially methylated domains. We further demonstrate that sequences eliminated during programmed genome rearrangement (PGR), are hypermethylated in sperm prior to the onset of PGR. Our study thus unveils important insights into the evolutionary origins of vertebrate 5mC reprogramming, and how this process might participate in diverse developmental strategies.

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.

Chimeric and mutant CARD9 constructs enable analyses of conserved and diverged autoinhibition mechanisms in the CARD-CC protein family

Caspase recruitment domain-containing protein (CARD)9, CARD10, CARD11, and CARD14 all belong to the CARD-coiled coil (CC) protein family and originated from a single common ancestral protein early in vertebrate evolution. All four proteins form CARD-CC/BCL10/MALT1 (CBM) complexes leading to nuclear factor-kappa-B (NF-κB) activation after upstream phosphorylation by various protein kinase C (PKC) isoforms. CBM complex signaling is critical for innate and adaptive immunity, but aberrant activation can cause autoimmune or autoinflammatory diseases, or be oncogenic. CARD9 shows a superior auto-inhibition compared with other CARD-CC family proteins, with very low spontaneous activity when overexpressed in HEK293T cells. In contrast, the poor auto-inhibition of other CARD-CC family proteins, especially CARD10 (CARMA3) and CARD14 (CARMA2), is hampering characterization of upstream activators or activating mutations in overexpression studies. We grafted different domains from CARD10, 11, and 14 on CARD9 to generate chimeric CARD9 backbones for functional characterization of activating mutants using NF-κB reporter gene activation in HEK293T cells as readout. CARD11 (CARMA1) activity was not further reduced by grafting on CARD9 backbones. The chimeric CARD9 approach was subsequently validated by using several known disease-associated mutations in CARD10 and CARD14, and additional screening allowed us to identify several previously unknown activating natural variants in human CARD9 and CARD10. Using Genebass as a resource of exome-based disease association statistics, we found that activated alleles of CARD9 correlate with irritable bowel syndrome (IBS), constipation, osteoarthritis, fibromyalgia, insomnia, anxiety, and depression, which can occur as comorbidities.

Gene expression in the liver of the hagfish (Eptatretus burgeri) belonging to the Cyclostomata is ancestral to that of mammals

Although the liver of the hagfish, an earliest diverged lineage among vertebrates, has a histological architecture similar to that of mammals, its gene expression has not been explored yet. The present study was undertaken to comparatively characterize gene expression in the liver of the hagfish with that of the mouse, using in situ hybridization technique. Expression of alb (albumin) was detectable in all hepatocytes of the hagfish liver, but was negative in intrahepatic bile ducts. Their expression in abundant periportal ductules was weak. The expression pattern basically resembled that in mammalian livers, indicating that the differential expression of hepatocyte markers in hepatocytes and biliary cells may have been acquired in ancestral vertebrates. alb expression was almost homogeneous in the hagfish liver, whereas that in the mouse liver lobule was zonal. The glul (glutamate-ammonia ligase) expression was also homogeneously detectable in hepatocytes without zonation, and weakly so in biliary cells of the hagfish, which contrasted with its restricted pericentral expression in mouse livers. These findings indicated that the hagfish liver did not have mammalian-type zonation. Whereas tetrapods had Hnf (hepatocyte nuclear factor) 1a and Hnf1b genes encoding the transcription factors, the hagfish had a single gene of their orthologue hnf1. Although HNF1α and HNF1β were immunohistochemically detected in hepatocytes and biliary cells of the mouse, respectively, hnf1 was expressed in both hepatocytes and biliary cells of the hagfish. These data indicate that gene expression of hnf1 in the hagfish liver may be ancestral with that of alb and glul during vertebrate evolution.

Models for the retention of duplicate genes and their biological underpinnings

Gene content in genomes changes through several different processes, with gene duplication being an important contributor to such changes. Gene duplication occurs over a range of scales from individual genes to whole genomes, and the dynamics of this process can be context dependent. Still, there are rules by which genes are retained or lost from genomes after duplication, and probabilistic modeling has enabled characterization of these rules, including their context-dependence. Here, we describe the biology and corresponding mathematical models that are used to understand duplicate gene retention and its contribution to the set of biochemical functions encoded in a genome.

Hox genes control homocercal caudal fin development and evolution

Ancient bony fishes had heterocercal tails, like modern sharks and sturgeons, with asymmetric caudal fins and a vertebral column extending into an elongated upper lobe. Teleost fishes, in contrast, developed a homocercal tail characterized by two separate equal-sized fin lobes and the body axis not extending into the caudal fin. A similar heterocercal-to-homocercal transition occurs during teleost ontogeny, although the underlying genetic and developmental mechanisms for either transition remain unresolved. Here, we investigated the role of hox13 genes in caudal fin formation as these genes control posterior identity in animals. Analysis of expression profiles of zebrafish hox13 paralogs and phenotypes of CRISPR/Cas9-induced mutants showed that double hoxb13a and hoxc13a mutants fail to form a caudal fin. Furthermore, single mutants display heterocercal-like morphologies not seen since Mesozoic fossil teleosteomorphs. Relaxation of functional constraints after the teleost genome duplication may have allowed hox13 duplicates to neo- or subfunctionalize, ultimately contributing to the evolution of a homocercal tail in teleost fishes.

Three foxg1 paralogues in lampreys and gnathostomes-brothers or cousins?

Foxg1 is a key regulator of the early development of the vertebrate forebrain and sensory organs. In this study, we describe for the first time three foxg1 paralogues in lamprey, representative of one of two basally diverged lineages of vertebrates-the agnathans. We also first describe three foxg1 genes in sterlet-representative of one of the evolutionarily ancient clades of gnathostomes. According to the analysis of local genomic synteny, three foxg1 genes of agnathans and gnathostomes have a common origin as a result of two rounds of genomic duplications in the early evolution of vertebrates. At the same time, it is difficult to reliably establish pairwise orthology between foxg1 genes of agnathans and gnathostomes based on the analysis of phylogeny and local genomic synteny, as well as our studies of the spatiotemporal expression of foxg1 genes in the river lamprey Lampetra fluviatilis and the sterlet Acipenser ruthenus. Thus, the appearance of three foxg1 paralogues in agnathans and gnathostomes could have occurred either as a result of two rounds of duplication of the vertebrate common ancestor genome (2R hypothesis) or as a result of the first common round followed by subsequent independent polyploidizations in two evolutionary lineages (1R hypothesis).

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Evolution of lysine-specific demethylase 1 and REST corepressor gene families and their molecular interaction

Lysine-specific demethylase 1A (LSD1) binds to the REST corepressor (RCOR) protein family of corepressors to erase transcriptionally active marks on histones. Functional diversity in these complexes depends on the type of RCOR included, which modulates the catalytic activity of the complex. Here, we studied the duplicative history of the RCOR and LSD gene families and analyzed the evolution of their interaction. We found that RCOR genes are the product of the two rounds of whole-genome duplications that occurred early in vertebrate evolution. In contrast, the origin of the LSD genes traces back before to the divergence of animals and plants. Using bioinformatics tools, we show that the RCOR and LSD1 interaction precedes the RCOR repertoire expansion that occurred in the last common ancestor of jawed vertebrates. Overall, we trace LSD1-RCOR complex evolution and propose that animal non-model species offer advantages in addressing questions about the molecular biology of this epigenetic complex.

Genomic view of the origins of cell-mediated immunity

NKp30 is an activating natural killer cell receptor (NKR) with a single-exon variable (VJ)-type immunoglobulin superfamily (IgSF) domain. Such VJ-IgSF domains predate the emergence of the antigen receptors (immunoglobulin and T cell receptor), which possess the same domain but undergo gene rearrangement. NCR3, the gene encoding NKp30, is present in jawed vertebrates from sharks to mammals; thus, unlike most NKR that are highly divergent among vertebrate taxa, NKp30 is uniquely conserved. We previously hypothesized that an ancestral NCR3 gene was encoded in the proto-major histocompatibility complex (MHC), the region where many immune-related genes have accumulated. Herein, we searched in silico databases to identify NCR3 paralogues and examined their genomic locations. We found a paralogue, NCR3H, in many vertebrates but was lost in mammals. Additionally, we identified a set of voltage-gated sodium channel beta (SCNB) genes as NCR3-distantly-related genes. Like NCR3, both NCR3H and SCNB proteins contain a single VJ-IgSF domain followed by a transmembrane region. These genes map to MHC paralogous regions, originally described in an invertebrate, along with genes encoding cell adhesion molecules involved in NK cell recognition networks. Other genes having no obvious relationship to immunity also map to these paralogous regions. These gene complexes were traced to several invertebrates, suggesting that the foundation of these cellular networks emerged before the genome-wide duplications in early gnathostome history. Here, we propose that this ancestral region was involved in cell-mediated immunity prior to the emergence of adaptive immunity and that NCR3 piggybacked onto this primordial complex, heralding the emergence of vertebrate NK cell/T cells.

Evolution and phylogenetic distribution of endo-α-mannosidase

While glycans underlie many biological processes, such as protein folding, cell adhesion, and cell-cell recognition, deep evolution of glycosylation machinery remains an understudied topic. N-linked glycosylation is a conserved process in which mannosidases are key trimming enzymes. One of them is the glycoprotein endo-α-1,2-mannosidase which participates in the initial trimming of mannose moieties from an N-linked glycan inside the cis-Golgi. It is unique as the only endo-acting mannosidase found in this organelle. Relatively little is known about its origins and evolutionary history; so far it was reported to occur only in vertebrates. In this work, a taxon-rich bioinformatic survey to unravel the evolutionary history of this enzyme, including all major eukaryotic clades and a wide representation of animals, is presented. The endomannosidase was found to be more widely distributed in animals and other eukaryotes. The protein motif changes in context of the canonical animal enzyme were tracked. Additionally, the data show the two canonical vertebrate endomannosidase genes, MANEA and MANEAL, arose at the second round of the two vertebrate genome duplications and one more vertebrate paralog, CMANEAL, is uncovered. Finally, a framework where N-glycosylation co-evolved with complex multicellularity is described. A better understanding of the evolution of core glycosylation pathways is pivotal to understanding biology of eukaryotes in general, and the Golgi apparatus in particular. This systematic analysis of the endomannosidase evolution is one step toward this goal.

Ohnologs and SSD Paralogs Differ in Genomic and Expression Features Related to Dosage Constraints

Gene duplication is recognized as a critical process in genome evolution; however, many questions about this process remain unanswered. Although gene duplicability has been observed to differ by duplication mechanism and evolutionary rate, there is so far no broad characterization of its determinants. Many features correlate with this difference in duplicability; however, our ability to exploit these observations to advance our understanding of the role of duplication in evolution is hampered by limitations within existing work. In particular, the existence of methodological differences across studies impedes meaningful comparison. Here, we use consistent definitions of duplicability in the human lineage to explore these associations, allow resolution of the impact of confounding factors, and define the overall relevance of individual features. Using a classifier approach and controlling for the confounding effect of duplicate longevity, we find a subset of gene features important in differentiating genes duplicable by small-scale duplication from those duplicable by whole-genome duplication, revealing critical roles for gene dosage and expression costs in duplicability. We further delve into patterns of functional enrichment and find a lack of constraint on duplicate retention in any context for genes duplicable by small-scale duplication.

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.

Evolution and divergence of teleost adrenergic receptors: why sometimes 'the drugs don't work' in fish

Adrenaline and noradrenaline, released as hormones and/or neurotransmitters, exert diverse physiological functions in vertebrates, and teleost fishes are widely used as model organisms to study adrenergic regulation; however, such investigations often rely on receptor subtype-specific pharmacological agents (agonists and antagonists; see Glossary) developed and validated in mammals. Meanwhile, evolutionary (phylogenetic and comparative genomic) studies have begun to unravel the diversification of adrenergic receptors (ARs) and reveal that whole-genome duplications and pseudogenization events in fishes results in notable distinctions from mammals in their genomic repertoire of ARs, while lineage-specific gene losses within teleosts have generated significant interspecific variability. In this Review, we visit the evolutionary history of ARs (including α1-, α2- and β-ARs) to highlight the prominent interspecific differences in teleosts, as well as between teleosts and other vertebrates. We also show that structural modelling of teleost ARs predicts differences in ligand binding affinity compared with mammalian orthologs. To emphasize the difficulty of studying the roles of different AR subtypes in fish, we collate examples from the literature of fish ARs behaving atypically compared with standard mammalian pharmacology. Thereafter, we focus on specific case studies of the liver, heart and red blood cells, where our understanding of AR expression has benefited from combining pharmacological approaches with molecular genetics. Finally, we briefly discuss the ongoing advances in 'omics' technologies that, alongside classical pharmacology, will provide abundant opportunities to further explore adrenergic signalling in teleosts.

Elasmobranch genome sequencing reveals evolutionary trends of vertebrate karyotype organization

Genomic studies of vertebrate chromosome evolution have long been hindered by the scarcity of chromosome-scale DNA sequences of some key taxa. One of those limiting taxa has been the elasmobranchs (sharks and rays), which harbor species often with numerous chromosomes and enlarged genomes. Here, we report the chromosome-scale genome assembly for the zebra shark Stegostoma tigrinum, an endangered species that has a relatively small genome among sharks (3.71 Gb), as well as for the whale shark Rhincodon typus Our analysis using a male-female comparison identified an X Chromosome, the first genomically characterized shark sex chromosome. The X Chromosome harbors the Hox C cluster whose intact linkage has not been shown for an elasmobranch fish. The sequenced shark genomes show a gradualism of chromosome length with remarkable length-dependent characteristics-shorter chromosomes tend to have higher GC content, gene density, synonymous substitution rate, and simple tandem repeat content as well as smaller gene length and lower interspersed repeat content. We challenge the traditional binary classification of karyotypes as with and without so-called microchromosomes. Even without microchromosomes, the length-dependent characteristics persist widely in nonmammalian vertebrates. Our investigation of elasmobranch karyotypes underpins their unique characteristics and provides clues for understanding how vertebrate karyotypes accommodate intragenomic heterogeneity to realize a complex readout. It also paves the way to dissecting more genomes with variable sizes to be sequenced at high quality.

Reduced steroid activation of elephant shark GR and MR after inserting four amino acids from the DNA-binding domain of lamprey corticoid receptor-1

Atlantic sea lamprey contains two corticoid receptors (CRs), CR1 and CR2, that have identical amino acid sequences, except for a four amino acid insert (Thr-Arg-Gln-Gly) in the CR1 DNA-binding domain (DBD). Steroids are stronger transcriptional activators of CR2 than of CR1 suggesting that the insert reduces the transcriptional response of lamprey CR1 to steroids. The DBD in elephant shark mineralocorticoid receptor (MR) and glucocorticoid receptor (GR), which are descended from a CR, lack these four amino acids, suggesting that a CR2 is their common ancestor. To determine if, similar to lamprey CR1, the presence of this insert in elephant shark MR and GR decreases transcriptional activation by corticosteroids, we inserted these four CR1-specific residues into the DBD of elephant shark MR and GR. Compared to steroid activation of wild-type elephant shark MR and GR, cortisol, corticosterone, aldosterone, 11-deoxycorticosterone and 11-deoxycortisol had lower transcriptional activation of these mutant MR and GR receptors, indicating that the absence of this four-residue segment in the DBD in wild-type elephant shark MR and GR increases transcriptional activation by corticosteroids.

Expansion and collapse of VEGF diversity in major clades of the animal kingdom

Together with the platelet-derived growth factors (PDGFs), the vascular endothelial growth factors (VEGFs) form the PDGF/VEGF subgroup among cystine knot growth factors. The evolutionary relationships within this subgroup have not been examined thoroughly to date. Here, we comprehensively analyze the PDGF/VEGF growth factors throughout all animal phyla and propose a phylogenetic tree. Vertebrate whole-genome duplications play a role in expanding PDGF/VEGF diversity, but several limited duplications are necessary to account for the temporal pattern of emergence. The phylogenetically oldest PDGF/VEGF-like growth factor likely featured a C-terminus with a BR3P signature, a hallmark of the modern-day lymphangiogenic growth factors VEGF-C and VEGF-D. Some younger VEGF genes, such as VEGFB and PGF, appeared completely absent in important vertebrate clades such as birds and amphibia, respectively. In contrast, individual PDGF/VEGF gene duplications frequently occurred in fish on top of the known fish-specific whole-genome duplications. The lack of precise counterparts for human genes poses limitations but also offers opportunities for research using organisms that diverge considerably from humans. Sources for the graphical abstract: 326 MYA and older [1]; 72-240 MYA [2]; 235-65 MYA [3].

Emergence of the Synucleins

This study establishes the origin and evolutionary history of the synuclein genes. A combination of phylogenetic analyses of the synucleins from twenty-two model species, characterization of local synteny similarities among humans, sharks and lampreys, and statistical comparisons among lamprey and human chromosomes, provides conclusive evidence for the current diversity of synuclein genes arising from the whole-genome duplications (WGDs) that occurred in vertebrates. An ancestral synuclein gene was duplicated in a first WGD, predating the diversification of all living vertebrates. The two resulting genes are still present in agnathan vertebrates. The second WGD, specific to the gnathostome lineage, led to the emergence of the three classical synuclein genes, SNCA, SNCB and SNCG, which are present in all jawed vertebrate lineages. Additional WGDs have added new genes in both agnathans and gnathostomes, while some gene losses have occurred in particular species. The emergence of synucleins through WGDs prevented these genes from experiencing dosage effects, thus avoiding the potential detrimental effects associated with individual duplications of genes that encode proteins prone to aggregation. Additional insights into the structural and functional features of synucleins are gained through the analysis of the highly divergent synuclein proteins present in chondrichthyans and agnathans.

Nitric Oxide Function and Nitric Oxide Synthase Evolution in Aquatic Chordates

Nitric oxide (NO) is a key signaling molecule in almost all organisms and is active in a variety of physiological and pathological processes. Our understanding of the peculiarities and functions of this simple gas has increased considerably by extending studies to non-mammal vertebrates and invertebrates. In this review, we report the nitric oxide synthase (Nos) genes so far characterized in chordates and provide an extensive, detailed, and comparative analysis of the function of NO in the aquatic chordates tunicates, cephalochordates, teleost fishes, and amphibians. This comprehensive set of data adds new elements to our understanding of Nos evolution, from the single gene commonly found in invertebrates to the three genes present in vertebrates.

A brief history of galectin evolution

Galectins are a family of carbohydrate-binding proteins found in vertebrates in great abundance and diversity in terms of both structure and ligand-binding properties as well as physiological function. Proteins with clear relationships to vertebrate galectins are already found in primitive Bilateria. The increasing amount of accessible well-annotated bilaterian genomes has allowed us to reveal, through synteny analyses, a new hypothesis about the phylogenetic history of the galectin family in this animal group. Thus, we can trace the genomic localization of the putative ancestral Bilateria galectin back to the scallops as a still very primitive slow-evolving bilaterian lineage. Intriguingly, our analyses show that the primordial galectin of the Deuterostomata most likely exhibited galectin-8-like characteristics. This basal standing galectin is characterized by a tandem-repeat type with two carbohydrate recognition domains as well as by a sialic acid binding property of the N-terminal domain, which is typical for galectin-8. With the help of synteny, the amplification of this potential primordial galectin to the broad galectin cosmos of modern jawed vertebrates can be reconstructed. Therefore, it is possible to distinguish between the paralogs resulting from small-scale duplication and the ohnologues generated by whole-genome duplication. Our findings support a substantially new hypothesis about the origin of the various members of the galectin family in vertebrates. This allows us to reveal new theories on the kinship relationships of the galectins of Gnatostomata. In addition, we focus for the first time on the galectines of the Cyclostomata, which as a sister group of jawed vertebrates providing important insights into the evolutionary history of the entire subphylum. Our studies also highlight a previously neglected member of the galectin family, galectin-related protein 2. This protein appears to be a widespread ohnologue of the original tandem-repeat ancestor within Gnathostomata that has not been the focus of galectin research due to its nonclassical galactose binding sequence motif and the fact that it was lost during mammalian evolution.

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

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

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

Phylogenetic and functional properties of hagfish neurohypophysial hormone receptors distinct from their jawed vertebrate counterparts

Vertebrate neurohypophysial hormones, i.e., vasopressin- and oxytocin-family peptides, exert versatile physiological actions via distinct G protein-coupled receptors. The neurohypophysial hormone receptor (NHR) family was classically categorized into four subtypes (V1aR, V1bR, V2R and OTR), while recent studies have identified seven subtypes (V1aR, V1bR, V2aR, V2bR, V2cR, V2dR and OTR; V2aR corresponds to the conventional V2R). The vertebrate NHR family were diversified via multiple gene duplication events at different scales. Despite intensive research effort in non-osteichthyes vertebrates such as cartilaginous fish and lamprey, the molecular phylogeny of the NHR family has not been fully understood. In the present study, we focused on the inshore hagfish (Eptatretus burgeri), another group of cyclostomes, and Arctic lamprey (Lethenteron camtschaticum) for comparison. Two putative NHR homologs, which were previously identified only in silico, were cloned from the hagfish and designated as ebV1R and ebV2R. In vitro, ebV1R, as well as two out of five Arctic lamprey NHRs, increased intracellular Ca2+ in response to exogenous neurohypophysial hormones. None of the examined cyclostome NHRs altered intracellular cAMP levels. Transcripts of ebV1R were detected in multiple tissues including the brain and gill, with intense hybridization signals in the hypothalamus and adenohypophysis, while ebV2R was predominantly expressed in the systemic heart. Similarly, Arctic lamprey NHRs showed distinct expression patterns, underscoring the multifunctionality of VT in the cyclostomes as in the gnathostomes. These results and exhaustive gene synteny comparisons provide new insights into the molecular and functional evolution of the neurohypophysial hormone system in vertebrates.

Molecular evolution and gene expression of ferritin family involved in immune defense of lampreys

Ferritin, one of the key regulators of iron homeostasis, is widely present throughout almost all species. The vertebrate ferritin family, which originates from a single gene in the ancestral invertebrates, contains the widest variety of ferritin subtypes among all animal species. However, the evolutionary history of the vertebrate ferritin family remains to be further clarified. In this study, genome-wide identification of the ferritin homologs is performed in lampreys, which are the extant representatives of jawless vertebrates that diverged from the future jawed vertebrates more than 500 million years ago. Molecular evolutionary analyses show that four members of the lamprey ferritin family, L-FT1-4, are derived from a common ancestor with jawed vertebrate ferritins prior to the divergence of the jawed vertebrate ferritin subtypes. The lamprey ferritin family shares evolutionarily conserved characteristics of the ferritin H subunit with higher vertebrates, but certain members such as L-FT1 additionally accumulate some features of the M or L subunits. Expression profiling reveals that lamprey ferritins are highly expressed in the liver. The transcription of L-FT1 is significantly induced in the liver and heart during lipopolysaccharide stimulation, indicating that L-FTs may play a role in the innate immune response to bacterial infection in lampreys. Furthermore, the transcriptional expression of L-FT1 in quiescent and LPS-activated leukocytes is up- and down-regulated by the lamprey TGF-β2, an essential regulator of the inflammatory response, respectively. Our results provide new insights into the origin and evolution of the vertebrate ferritin family and reveal that lamprey ferritins may be involved in immune regulation as target genes of the TGF-β signaling pathway.

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.

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

An improved germline genome assembly for the sea lamprey Petromyzon marinus illuminates the evolution of germline-specific chromosomes

Programmed DNA loss is a gene silencing mechanism that is employed by several vertebrate and nonvertebrate lineages, including all living jawless vertebrates and songbirds. Reconstructing the evolution of somatically eliminated (germline-specific) sequences in these species has proven challenging due to a high content of repeats and gene duplications in eliminated sequences and a corresponding lack of highly accurate and contiguous assemblies for these regions. Here, we present an improved assembly of the sea lamprey (Petromyzon marinus) genome that was generated using recently standardized methods that increase the contiguity and accuracy of vertebrate genome assemblies. This assembly resolves highly contiguous, somatically retained chromosomes and at least one germline-specific chromosome, permitting new analyses that reconstruct the timing, mode, and repercussions of recruitment of genes to the germline-specific fraction. These analyses reveal major roles of interchromosomal segmental duplication, intrachromosomal duplication, and positive selection for germline functions in the long-term evolution of germline-specific chromosomes.

Na+/Cl- cotransporter 2 is not fish-specific and is widely found in amphibians, non-avian reptiles, and select mammals

Solute carrier 12 (Slc12) is a family of electroneutral cation-coupled chloride (Cl-) cotransporters. Na+/K+/2Cl- (Nkcc) and Na+/Cl- cotransporters (Ncc) belong to the Nkcc/Ncc subfamily. Human and mouse possess one gene for the Na+/Cl- cotransporter (ncc gene: slc12a3), whereas teleost fishes possess multiple ncc genes, slc12a3 (ncc1) and slc12a10 (ncc2), in addition to their species-specific paralogs. Amphibians and squamates have two ncc genes: slc12a3 (ncc1) and ncc3. However, the evolutionary relationship between slc12a10 and ncc3 remains unresolved, and the presence of slc12a10 (ncc2) in mammals has not been clarified. Synteny and phylogenetic analyses of vertebrate genome databases showed that ncc3 is the ortholog of slc12a10, and slc12a10 is present in most ray-finned fishes, coelacanths, amphibians, reptiles, and a few mammals (e.g., platypus and horse) but pseudogenized or deleted in birds, most mammals, and some ray-finned fishes (pufferfishes). This shows that slc12a10 is widely present among bony vertebrates and pseudogenized or deleted independently in multiple lineages. Notably, as compared with some fish that show varied slc12a10 tissue expression profile, spotted gar, African clawed frog, red-eared slider turtle, and horse express slc12a10 in the ovaries or premature gonads. In horse tissues, an unexpectedly large number of splicing variants for Slc12a10 have been cloned, many of which encode truncated forms of Slc12a10, suggesting that the functional constraints of horse slc12a10 are weakened, which may be in the process of becoming a pseudogene. Our results elaborate on the evolution of Nkcc/Ncc subfamily of Slc12 in vertebrates.NEW & NOTEWORTHY slc12a10 is not a fish-specific gene and is present in a few mammals (e.g., platypus and horse), non-avian reptiles, amphibians, but was pseudogenized or deleted in most mammals (e.g., human, mouse, cat, cow, and rhinoceros), birds, and some ray-finned fishes (pufferfishes).

Evolutionary analysis of a complete chicken genome

Microchromosomes are prevalent in nonmammalian vertebrates [P. D. Waters et al., Proc. Natl. Acad. Sci. U.S.A. 118 (2021)], but a few of them are missing in bird genome assemblies. Here, we present a new chicken reference genome containing all autosomes, a Z and a W chromosome, with all gaps closed except for the W. We identified ten small microchromosomes (termed dot chromosomes) with distinct sequence and epigenetic features, among which six were newly assembled. Those dot chromosomes exhibit extremely high GC content and a high level of DNA methylation and are enriched for housekeeping genes. The pericentromeric heterochromatin of dot chromosomes is disproportionately large and continues to expand with the proliferation of satellite DNA and testis-expressed genes. Our analyses revealed that the 41-bp CNM repeat frequently forms higher-order repeats (HORs) at the centromeres of acrocentric chromosomes. The centromere core regions where the kinetochore attaches often encompass telomeric sequence (TTAGGG)n, and in a one of the dot chromosomes, the centromere core recruits an endogenous retrovirus (ERV). We further demonstrate that the W chromosome shares some common features with dot chromosomes, having large arrays of hypermethylated tandem repeats. Finally, using the complete chicken chromosome models, we reconstructed a fine picture of chordate karyotype evolution, revealing frequent chromosomal fusions before and after vertebrate whole-genome duplications. Our sequence and epigenetic characterization of chicken chromosomes shed insights into the understanding of vertebrate genome evolution and chromosome biology.

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.

Single-cell transcriptomics of the goldfish retina reveals genetic divergence in the asymmetrically evolved subgenomes after allotetraploidization

The recent whole-genome duplication (WGD) in goldfish (Carassius auratus) approximately 14 million years ago makes it a valuable model for studying gene evolution during the early stages after WGD. We analyzed the transcriptome of the goldfish retina at the level of single-cell (scRNA-seq) and open chromatin regions (scATAC-seq). We identified a group of genes that have undergone dosage selection, accounting for 5% of the total 11,444 ohnolog pairs. We also identified 306 putative sub/neo-functionalized ohnolog pairs that are likely to be under cell-type-specific genetic variation at single-cell resolution. Diversification in the expression patterns of several ohnolog pairs was observed in the retinal cell subpopulations. The single-cell level transcriptome analysis in this study uncovered the early stages of evolution in retinal cell of goldfish after WGD. Our results provide clues for understanding the relationship between the early stages of gene evolution after WGD and the evolution of diverse vertebrate retinal functions.

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.

Sharks Provide Evidence for a Highly Complex TNFSF Repertoire in the Jawed Vertebrate Ancestor

Cytokines of the TNF superfamily (TNFSF) control many immunological processes and are implicated in the etiology of many immune disorders and diseases. Despite their obvious biological importance, the TNFSF repertoires of many species remain poorly characterized. In this study, we perform detailed bioinformatic, phylogenetic, and syntenic analyses of five cartilaginous fish genomes to identify their TNFSF repertoires. Strikingly, we find that shark genomes harbor ∼30 TNFSF genes, more than any other vertebrate examined to date and substantially more than humans. This is due to better retention of the ancestral jawed vertebrate TNFSF repertoire than any other jawed vertebrate lineage, combined with lineage-specific gene family expansions. All human TNFSFs appear in shark genomes, except for lymphotoxin-α (LTA; TNFSF1) and TNF (TNFSF2), and CD70 (TNFSF7) and 4-1BBL (TNFSF9), which diverged by tandem duplications early in tetrapod and mammalian evolution, respectively. Although lacking one-to-one LTA and TNF orthologs, sharks have evolved lineage-specific clusters of LTA/TNF co-orthologs. Other key findings include the presence of two BAFF (TNFSF13B) genes along with orthologs of APRIL (TNFSF13) and BALM (TNFSF13C) in sharks, and that all cartilaginous fish genomes harbor an ∼400-million-year-old cluster of multiple FASLG (TNFSF6) orthologs. Finally, sharks have retained seven ancestral jawed vertebrate TNFSF genes lost in humans. Taken together, our data indicate that the jawed vertebrate ancestor possessed a much larger and diverse TNFSF repertoire than previously hypothesized and oppose the idea that the cartilaginous fish immune system is “primitive” compared with that of mammals.

Photoreceptor physiology and evolution: cellular and molecular basis of rod and cone phototransduction

The detection of light in the vertebrate retina utilizes a duplex system of closely related rod and cone photoreceptors: cones respond extremely rapidly, and operate at 'photopic' levels of illumination, from moonlight upwards; rods respond much more slowly, thereby obtaining greater sensitivity, and function effectively only at 'scotopic' levels of moonlight and lower. Rods and cones employ distinct isoforms of many of the proteins in the phototransduction cascade, and they thereby represent a unique evolutionary system, whereby the same process (the detection of light) uses a distinct set of genes in two classes of cell. The molecular mechanisms of phototransduction activation are described, and the classical quantitative predictions for the onset phase of the electrical response to light are developed. Recent work predicting the recovery phase of the rod's response to intense flashes is then presented, that provides an accurate account of the time that the response spends in saturation. Importantly, this also provides a new estimate for the rate at which a single rhodopsin activates molecules of the G-protein, transducin, that is substantially higher than other estimates in the literature. Finally, the evolutionary origin of the phototransduction proteins in rods and cones is examined, and it is shown that most of the rod/cone differences were established at the first of the two rounds of whole-genome duplication more than 500 million years ago.

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.

Cytokine Receptor Diversity in the Lamprey Predicts the Minimal Essential Cytokine Networks of Vertebrates

The vertebrate adaptive immune systems (Agnatha and Gnathostomata) use sets of T and B lymphocyte lineages that somatically generate highly diverse repertoires of Ag-specific receptors and Abs. In Gnathostomata, cytokine networks regulate the activation of lymphoid and myeloid cells, whereas little is known about these components in Agnathans. Most gnathostome cytokines are four-helix bundle cytokines with poorly conserved primary sequences. In contrast, sequence conservation across bilaterians has been observed for cognate cytokine receptor chains, allowing their structural classification into two classes, and for downstream JAK/STAT signaling mediators. With conserved numbers among Gnathostomata, human cytokine receptor chains (comprising 34 class I and 12 class II) are able to interact with 28 class I helical cytokines (including most ILs) and 16 class II cytokines (including all IFNs), respectively. Hypothesizing that the arsenal of cytokine receptors and transducers may reflect homologous cytokine networks, we analyzed the lamprey genome and transcriptome to identify genes and transcripts for 23 class I and five class II cytokine receptors alongside one JAK signal mediator and four STAT transcription factors. On the basis of deduction of their respective orthologs, we predict that these receptors may interact with 16 class I and 3 class II helical cytokines (including IL-4, IL-6, IL-7, IL-12, IL-10, IFN-γ, and thymic stromal lymphoprotein homologs). On the basis of their respective activities in mammals, this analysis suggests the existence of lamprey cytokine networks that may regulate myeloid and lymphoid cell differentiation, including potential Th1/Th2 polarization. The predicted networks thus appear remarkably homologous to those of Gnathostomata, albeit reduced to essential functions.

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.

Evolution and developmental expression of the sodium-iodide symporter (NIS, slc5a5) gene family: Implications for perchlorate toxicology

The vertebrate sodium-iodide symporter (NIS or SLC5A5) transports iodide into the thyroid follicular cells that synthesize thyroid hormone. The SLC5A protein family includes transporters of vitamins, minerals, and nutrients. Disruption of SLC5A5 function by perchlorate, a pervasive environmental contaminant, leads to human pathologies, especially hypothyroidism. Perchlorate also disrupts the sexual development of model animals, including threespine stickleback (Gasterosteus aculeatus) and zebrafish (Danio rerio), but the mechanism of action is unknown. To test the hypothesis that SLC5A5 paralogs are expressed in tissues necessary for the development of reproductive organs, and therefore are plausible candidates to mediate the effects of perchlorate on sexual development, we first investigated the evolutionary history of Slc5a paralogs to better understand potential functional trajectories of the gene family. We identified two clades of slc5a paralogs with respect to an outgroup of sodium/choline cotransporters (slc5a7); these clades are the NIS clade of sodium/iodide and lactate cotransporters (slc5a5, slc5a6, slc5a8, slc5a8, and slc5a12) and the SGLT clade of sodium/glucose cotransporters (slc5a1, slc5a2, slc5a3, slc5a4, slc5a10, and slc5a11). We also characterized expression patterns of slc5a genes during development. Stickleback embryos and early larvae expressed NIS clade genes in connective tissue, cartilage, teeth, and thyroid. Stickleback males and females expressed slc5a5 and its paralogs in gonads. Single-cell transcriptomics (scRNA-seq) on zebrafish sex-genotyped gonads revealed that NIS clade-expressing cells included germ cells (slc5a5, slc5a6a, and slc5a6b) and gonadal soma cells (slc5a8l). These results are consistent with the hypothesis that perchlorate exerts its effects on sexual development by interacting with slc5a5 or its paralogs in reproductive tissues. These findings show novel expression domains of slc5 genes in stickleback and zebrafish, which suggest similar functions across vertebrates including humans, and provide candidates to mediate the effects of perchlorate on sexual development.

A glimpse into the past: phylogenesis and protein domain analysis of the group XIV of C-type lectins in vertebrates

Background

The group XIV of C-type lectin domain-containing proteins (CTLDcps) is one of the seventeen groups of CTLDcps discovered in mammals and composed by four members: CD93, Clec14A, CD248 and Thrombomodulin, which have shown to be important players in cancer and vascular biology. Although these proteins belong to the same family, their phylogenetic relationship has never been dissected. To resolve their evolution and characterize their protein domain composition we investigated CTLDcp genes in gnathostomes and cyclostomes and, by means of phylogenetic approaches as well as synteny analyses, we inferred an evolutionary scheme that attempts to unravel their evolution in modern vertebrates.

Results

Here, we evidenced the paralogy of the group XIV of CTLDcps in gnathostomes and discovered that a gene loss of CD248 and Clec14A occurred in different vertebrate groups, with CD248 being lost due to chromosome disruption in birds, while Clec14A loss in monotremes and marsupials did not involve chromosome rearrangements. Moreover, employing genome annotations of different lampreys as well as one hagfish species, we investigated the origin and evolution of modern group XIV of CTLDcps. Furthermore, we carefully retrieved and annotated gnathostome CTLDcp domains, pointed out important differences in domain composition between gnathostome classes, and assessed codon substitution rate of each domain by analyzing nonsynonymous (Ka) over synonymous (Ks) substitutions using one representative species per gnathostome order.

Conclusions

CTLDcps appeared with the advent of early vertebrates after a whole genome duplication followed by a sporadic tandem duplication. These duplication events gave rise to three CTLDcps in the ancestral vertebrate that underwent further duplications caused by the independent polyploidizations that characterized the evolution of cyclostomes and gnathostomes. Importantly, our analyses of CTLDcps in gnathostomes revealed critical inter-class differences in both extracellular and intracellular domains, which might help the interpretation of experimental results and the understanding of differences between animal models.

Pervasive male-biased expression throughout the germline-specific regions of the sea lamprey genome supports key roles in sex differentiation and spermatogenesis

Sea lamprey undergo programmed genome rearrangement (PGR) in which ∼20% of the genome is jettisoned from somatic cells during embryogenesis. Although the role of PGR in embryonic development has been studied, the role of the germline-specific region (GSR) in gonad development is unknown. We analysed RNA-sequence data from 28 sea lamprey gonads sampled across life-history stages, generated a genome-guided de novo superTranscriptome with annotations, and identified germline-specific genes (GSGs). Overall, we identified 638 GSGs that are enriched for reproductive processes and exhibit 36x greater odds of being expressed in testes than ovaries. Next, while 55% of the GSGs have putative somatic paralogs, the somatic paralogs are not differentially expressed between sexes. Further, putative orthologs of some the male-biased GSGs have known functions in sex determination or differentiation in other vertebrates. We conclude that the GSR of sea lamprey plays an important role in testicular differentiation and potentially sex determination.

RNA-sequencing of sea lamprey gonads at different life-history stage identifies germline-specific genes which are highly expressed in males during spermatogenesis. This suggests a link between male-biased germline expression and sex differentiation in the sea lamprey.

Comparison of Reptilian Genomes Reveals Deletions Associated with the Natural Loss of γδ T Cells in Squamates

T lymphocytes or T cells are key components of the vertebrate response to pathogens and cancer. There are two T cell classes based on their TCRs, αβ T cells and γδ T cells, and each plays a critical role in immune responses. The squamate reptiles may be unique among the vertebrate lineages by lacking an entire class of T cells, the γδ T cells. In this study, we investigated the basis of the loss of the γδ T cells in squamates. The genome and transcriptome of a sleepy lizard, the skink Tiliqua rugosa, were compared with those of tuatara, Sphenodon punctatus, the last living member of the Rhynchocephalian reptiles. We demonstrate that the lack of TCRγ and TCRδ transcripts in the skink are due to large deletions in the T. rugosa genome. We also show that tuataras are on a growing list of species, including sharks, frogs, birds, alligators, and platypus, that can use an atypical TCRδ that appears to be a chimera of a TCR chain with an Ab-like Ag-binding domain. Tuatara represents the nearest living relative to squamates that retain γδ T cells. The loss of γδTCR in the skink is due to genomic deletions that appear to be conserved in other squamates. The genes encoding the αβTCR chains in the skink do not appear to have increased in complexity to compensate for the loss of γδ T cells.

Evolution of Somite Compartmentalization: A View From Xenopus

Somites are transitory metameric structures at the basis of the axial organization of vertebrate musculoskeletal system. During evolution, somites appear in the chordate phylum and compartmentalize mainly into the dermomyotome, the myotome, and the sclerotome in vertebrates. In this review, we summarized the existing literature about somite compartmentalization in Xenopus and compared it with other anamniote and amniote vertebrates. We also present and discuss a model that describes the evolutionary history of somite compartmentalization from ancestral chordates to amniote vertebrates. We propose that the ancestral organization of chordate somite, subdivided into a lateral compartment of multipotent somitic cells (MSCs) and a medial primitive myotome, evolves through two major transitions. From ancestral chordates to vertebrates, the cell potency of MSCs may have evolved and gave rise to all new vertebrate compartments, i.e., the dermomyome, its hypaxial region, and the sclerotome. From anamniote to amniote vertebrates, the lateral MSC territory may expand to the whole somite at the expense of primitive myotome and may probably facilitate sclerotome formation. We propose that successive modifications of the cell potency of some type of embryonic progenitors could be one of major processes of the vertebrate evolution.