Timing of genome duplications relative to the origin of the vertebrates: did cyclostomes diverge before or after?

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.

Divergence time estimates for the hypoxia-inducible factor-1 alpha (HIF1α) reveal an ancient emergence of animals in low-oxygen environments

Unveiling the tempo and mode of animal evolution is necessary to understand the links between environmental changes and biological innovation. Although the earliest unambiguous metazoan fossils date to the late Ediacaran period, molecular clock estimates agree that the last common ancestor (LCA) of all extant animals emerged ~850 Ma, in the Tonian period, before the oldest evidence for widespread ocean oxygenation at ~635-560 Ma in the Ediacaran period. Metazoans are aerobic organisms, that is, they are dependent on oxygen to survive. In low-oxygen conditions, most animals have an evolutionarily conserved pathway for maintaining oxygen homeostasis that triggers physiological changes in gene expression via the hypoxia-inducible factor (HIFa). However, here we confirm the absence of the characteristic HIFa protein domain responsible for the oxygen sensing of HIFa in sponges and ctenophores, indicating the LCA of metazoans lacked the functional protein domain as well, and so could have maintained their transcription levels unaltered under the very low-oxygen concentrations of their environments. Using Bayesian relaxed molecular clock dating, we inferred that the ancestral gene lineage responsible for HIFa arose in the Mesoproterozoic Era, ~1273 Ma (Credibility Interval 957-1621 Ma), consistent with the idea that important genetic machinery associated with animals evolved much earlier than the LCA of animals. Our data suggest at least two duplication events in the evolutionary history of HIFa, which generated three vertebrate paralogs, products of the two successive whole-genome duplications that occurred in the vertebrate LCA. Overall, our results support the hypothesis of a pre-Tonian emergence of metazoans under low-oxygen conditions, and an increase in oxygen response elements during animal evolution.

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.

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.

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.

Evolutionary Overview of Aquaporin Superfamily

Aquaporins (AQPs) are present not only in three domains of life, bacteria, eukaryotes, and archaea, but also in viruses. With the accumulating arrays of AQP superfamily, the evolutional relationship has attracted much attention with multiple publications on "the genome-wide identification and phylogenetic analysis" of AQP superfamily. A pair of NPA boxes forming a pore is highly conserved throughout the evolution and renders key residues for the classification of AQP superfamily into four groups: AQP1-like, AQP3-like, AQP8-like, and AQP11-like. The complexity of AQP family has mostly been achieved in nematodes and subsequent evolution has been directed toward increasing the number of AQPs through whole-genome duplications (WGDs) to extend the tissue specific expression and regulation. The discovery of the intracellular AQP (iAQP: AQP8-like and AQP11-like) and substrate transports by the plasma membrane AQP (pAQP: AQP1-like and AQP3-like) have accelerated the AQP research much more toward the transport of substrates with complex profiles. This evolutionary overview based on a simple classification of AQPs into four subfamilies will provide putative structural, functional, and localization information and insights into the role of AQP as well as clues to understand the complex diversity of AQP superfamily.

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.

Evolution of calcitonin/calcitonin gene-related peptide family in chordates: Identification of CT/CGRP family peptides in cartilaginous fish genome

The calcitonin (CT)/CT gene-related peptide (CGRP) family is a peptide gene family that is widely found in bilaterians. CT, CGRP, adrenomedullin (AM), amylin (AMY), and CT receptor-stimulating peptide (CRSP) are members of the CT/CGRP family. In mammals, CT is involved in calcium homeostasis, while CGRP and AM primarily function in vasodilation. AMY and CRSP are associated with anorectic effects. Diversification of the molecular features and physiological functions of the CT/CGRP family in vertebrate lineages have been extensively reported. However, the origin and diversification mechanisms of the vertebrate CT/CGRP family of peptides remain unclear. In this review, the molecular characteristics of CT/CGRP family peptides and their receptors, along with their major physiological functions in mammals and teleosts, are introduced. Furthermore, novel candidates of the CT/CGRP family in cartilaginous fish are presented based on genomic information. The CT/CGRP family peptides and receptors in urochordates and cephalochordates, which are closely related to vertebrates, are also described. Finally, a putative evolutionary scenario of the CT/CGRP family peptides and receptors in chordates is discussed.

Genome-Wide Identification and Characterization of the BRD Family in Nile Tilapia (Oreochromis niloticus)

Simple Summary

Nile tilapia is a good model for genome-wide identification and examination of the expression and role of gene families. In this study, we identified 54 bromodomain genes (BRDs) divided into eight subfamilies in Nile tilapia. Phylogenetic analysis revealed a high conservation of the BRDs family in vertebrates, with BRDs expansion due to fish-specific duplications. Most of the BRDs displayed sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and 9 ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testes of treated fish.

Abstract

The bromodomain (BRD) proteins specifically recognize the N-acetyllysine motifs, which is a key event in the reading process of epigenetic marks. BRDs are evolutionarily highly conserved. Over recent years, BRDs attracted great interest because of their important roles in biological processes. However, the genome-wide identification of this family was not carried out in many animal groups, in particular, in teleosts. Moreover, the expression patterns were not reported for any of the members in this family, and the role of the BRD family was not extensively studied in fish reproduction. In this study, we identified 16 to 120 BRD genes in 24 representative species. BRDs expanded significantly in vertebrates. Phylogenetic analysis showed that the BRD family was divided into eight subfamilies (I–VIII). Transcriptome analysis showed that BRDs in Nile tilapia (Oreochromis niloticus) exhibited different expression patterns in different tissues, suggesting that these genes may play different roles in growth and development. Gonadal transcriptome analysis showed that most of the BRDs display sexually dimorphic expression in the gonads at 90 and 180 dah (days after hatching), including 21 testis-dominated genes (brdt, brd4a and brd2b, etc.), and nine ovary-dominated genes (brd3b, brd2a and kat2a, etc.). Consistent with transcriptomic data, the results of qRT-PCR and fluorescence in situ hybridization showed that brdt expression was higher in the testis than in the ovary, suggesting its critical role in the spermatogenesis of the tilapia. Male fish treated with JQ1 (BET subfamily inhibitor) displayed abnormal spermatogenesis. The numbers of germ cells were reduced, and the expression of steroidogenic enzyme genes was downregulated, while the expression of apoptosis-promoting genes was elevated in the testis tissue of treated fish. Our data provide insights into the evolution and expression of BRD genes, which is helpful for understanding their critical roles in sex differentiation and gonadal development in teleosts.

An EvoDevo Study of Salmonid Visual Opsin Dynamics and Photopigment Spectral Sensitivity

Salmonids are ideal models as many species follow a distinct developmental program from demersal eggs and a large yolk sac to hatching at an advanced developmental stage. Further, these economically important teleosts inhabit both marine- and freshwaters and experience diverse light environments during their life histories. At a genome level, salmonids have undergone a salmonid-specific fourth whole genome duplication event (Ss4R) compared to other teleosts that are already more genetically diverse compared to many non-teleost vertebrates. Thus, salmonids display phenotypically plastic visual systems that appear to be closely related to their anadromous migration patterns. This is most likely due to a complex interplay between their larger, more gene-rich genomes and broad spectrally enriched habitats; however, the molecular basis and functional consequences for such diversity is not fully understood. This study used advances in genome sequencing to identify the repertoire and genome organization of visual opsin genes (those primarily expressed in retinal photoreceptors) from six different salmonids [Atlantic salmon (Salmo salar), brown trout (Salmo trutta), Chinook salmon (Oncorhynchus tshawytcha), coho salmon (Oncorhynchus kisutch), rainbow trout (Oncorhynchus mykiss), and sockeye salmon (Oncorhynchus nerka)] compared to the northern pike (Esox lucius), a closely related non-salmonid species. Results identified multiple orthologues for all five visual opsin classes, except for presence of a single short-wavelength-sensitive-2 opsin gene. Several visual opsin genes were not retained after the Ss4R duplication event, which is consistent with the concept of salmonid rediploidization. Developmentally, transcriptomic analyzes of Atlantic salmon revealed differential expression within each opsin class, with two of the long-wavelength-sensitive opsins not being expressed before first feeding. Also, early opsin expression in the retina was located centrally, expanding dorsally and ventrally as eye development progressed, with rod opsin being the dominant visual opsin post-hatching. Modeling by spectral tuning analysis and atomistic molecular simulation, predicted the greatest variation in the spectral peak of absorbance to be within the Rh2 class, with a ∼40 nm difference in λ max values between the four medium-wavelength-sensitive photopigments. Overall, it appears that opsin duplication and expression, and their respective spectral tuning profiles, evolved to maximize specialist color vision throughout an anadromous lifecycle, with some visual opsin genes being lost to tailor marine-based vision.

G protein-coupled receptors (GPCRs) in rotifers and cladocerans: Potential applications in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology

The G protein-coupled receptor (GPCR) superfamily plays a fundamental role in both sensory functions and the regulation of homeostasis, and is highly conserved across the eukaryote taxa. Its functional diversity is related to a conserved seven-transmembrane core and invariant set of intracellular signaling mechanisms. The interplay between these properties is key to the evolutionary success of GPCR. As this superfamily originated from a common ancestor, GPCR genes have evolved via lineage-specific duplications through the process of adaptation. Here we summarized information on GPCR gene families in rotifers and cladocerans based on their evolutionary position in aquatic invertebrates and their potential application in ecotoxicology, ecophysiology, comparative endocrinology, and pharmacology. Phylogenetic analyses were conducted to examine the evolutionary significance of GPCR gene families and to provide structural insight on their role in aquatic invertebrates. In particular, most GPCR gene families have undergone sporadic evolutionary processes, but some GPCRs are highly conserved across species despite the dynamics of GPCR evolution. Overall, this review provides a better understanding of GPCR evolution in aquatic invertebrates and expand our knowledge of the potential application of these receptors in various fields.

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.

Comprehensive expression analysis for the core cell cycle regulators in the chicken embryo reveals novel tissue-specific synexpression groups and similarities and differences with expression in mouse, frog and zebrafish

The core cell cycle machinery is conserved from yeast to humans, and hence it is assumed that all vertebrates share the same set of players. Yet during vertebrate evolution, the genome was duplicated twice, followed by a further genome duplication in teleost fish. Thereafter, distinct genes were retained in different vertebrate lineages; some individual gene duplications also occurred. To which extent these diversifying tendencies were compensated by retaining the same expression patterns across homologous genes is not known. This study for the first time undertook a comprehensive expression analysis for the core cell cycle regulators in the chicken, focusing in on early neurula and pharyngula stages of development, with the latter representing the vertebrate phylotypic stage. We also compared our data with published data for the mouse, Xenopus and zebrafish, the other established vertebrate models. Our work shows that, while many genes are expressed widely, some are upregulated or specifically expressed in defined tissues of the chicken embryo, forming novel synexpression groups with markers for distinct developmental pathways. Moreover, we found that in the neural tube and in the somite, mRNAs of some of the genes investigated accumulate in a specific subcellular localisation, pointing at a novel link between the site of mRNA translation, cell cycle control and interkinetic nuclear movements. Finally, we show that expression patterns of orthologous genes may differ in the four vertebrate models. Thus, for any study investigating cell proliferation, cell differentiation, tissue regeneration, stem cell behaviour and cancer/cancer therapy, it has to be carefully examined which of the observed effects are due to the specific model organism used, and which can be generalised.

Transcriptomic Analyses Reveal 2 and 4 Family Members of Cytochromes P450 (CYP) Involved in LPS Inflammatory Response in Pharynx of Ciona robusta

Cytochromes P450 (CYP) are enzymes responsible for the biotransformation of most endogenous and exogenous agents. The expression of each CYP is influenced by a unique combination of mechanisms and factors including genetic polymorphisms, induction by xenobiotics, and regulation by cytokines and hormones. In recent years, Ciona robusta, one of the closest living relatives of vertebrates, has become a model in various fields of biology, in particular for studying inflammatory response. Using an in vivo LPS exposure strategy, next-generation sequencing (NGS) and qRT-PCR combined with bioinformatics and in silico analyses, compared whole pharynx transcripts from naïve and LPS-exposed C. robusta, and we provide the first view of cytochrome genes expression and miRNA regulation in the inflammatory response induced by LPS in a hematopoietic organ. In C. robusta, cytochromes belonging to 2B,2C, 2J, 2U, 4B and 4F subfamilies were deregulated and miRNA network interactions suggest that different conserved and species-specific miRNAs are involved in post-transcriptional regulation of cytochrome genes and that there could be an interplay between specific miRNAs regulating both inflammation and cytochrome molecules in the inflammatory response in C. robusta.

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.

High degree of conservation of the enzymes synthesizing the laminin-binding glycoepitope of α-dystroglycan

The dystroglycan (DG) complex plays a pivotal role for the stabilization of muscles in Metazoa. It is formed by two subunits, extracellular α-DG and transmembrane β-DG, originating from a unique precursor via a complex post-translational maturation process. The α-DG subunit is extensively glycosylated in sequential steps by several specific enzymes and employs such glycan scaffold to tightly bind basement membrane molecules. Mutations of several of these enzymes cause an alteration of the carbohydrate structure of α-DG, resulting in severe neuromuscular disorders collectively named dystroglycanopathies. Given the fundamental role played by DG in muscle stability, it is biochemically and clinically relevant to investigate these post-translational modifying enzymes from an evolutionary perspective. A first phylogenetic history of the thirteen enzymes involved in the fabrication of the so-called 'M3 core' laminin-binding epitope has been traced by an overall sequence comparison approach, and interesting details on the primordial enzyme set have emerged, as well as substantial conservation in Metazoa. The optimization along with the evolution of a well-conserved enzymatic set responsible for the glycosylation of α-DG indicate the importance of the glycosylation shell in modulating the connection between sarcolemma and surrounding basement membranes to increase skeletal muscle stability, and eventually support movement and locomotion.

The Natterin Proteins Diversity: A Review on Phylogeny, Structure, and Immune Function

Since the first record of the five founder members of the group of Natterin proteins in the venom of the medically significant fish Thalassophryne nattereri, new sequences have been identified in other species. In this work, we performed a detailed screening using available genome databases across a wide range of species to identify sequence members of the Natterin group, sequence similarities, conserved domains, and evolutionary relationships. The high-throughput tools have enabled us to dramatically expand the number of members within this group of proteins, which has a remote origin (around 400 million years ago) and is spread across Eukarya organisms, even in plants and primitive Agnathans jawless fish. Overall, the survey resulted in 331 species presenting Natterin-like proteins, mainly fish, and 859 putative genes. Besides fish, the groups with more species included in our analysis were insects and birds. The number and variety of annotations increased the knowledge of the obtained sequences in detail, such as the conserved motif AGIP in the pore-forming loop involved in the transmembrane barrel insertion, allowing us to classify them as important constituents of the innate immune defense system as effector molecules activating immune cells by interacting with conserved intracellular signaling mechanisms in the hosts.

Molecular Phylogenetic Analysis of the AIG Family in Vertebrates

Androgen-inducible genes (AIGs), which can be regulated by androgen level, constitute a group of genes characterized by the presence of the AIG/FAR-17a domain in its protein sequence. Previous studies on AIGs demonstrated that one member of the gene family, AIG1, is involved in many biological processes in cancer cell lines and that ADTRP is associated with cardiovascular diseases. It has been shown that the numbers of AIG paralogs in humans, mice, and zebrafish are 2, 2, and 3, respectively, indicating possible gene duplication events during vertebrate evolution. Therefore, classifying subgroups of AIGs and identifying the homologs of each AIG member are important to characterize this novel gene family further. In this study, vertebrate AIGs were phylogenetically grouped into three major clades, ADTRP, AIG1, and AIG-L, with AIG-L also evident in an outgroup consisting of invertebrsate species. In this case, AIG-L, as the ancestral AIG, gave rise to ADTRP and AIG1 after two rounds of whole-genome duplications during vertebrate evolution. Then, the AIG family, which was exposed to purifying forces during evolution, lost or gained some of its members in some species. For example, in eutherians, Neognathae, and Percomorphaceae, AIG-L was lost; in contrast, Salmonidae and Cyprinidae acquired additional AIG copies. In conclusion, this study provides a comprehensive molecular phylogenetic analysis of vertebrate AIGs, which can be employed for future functional characterization of AIGs.

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

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

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.

Bird aquaporins: Molecular machinery for urine concentration

Water conservation is one of the most challenging processes for terrestrial vertebrates and is necessary for their survival. Birds are the only vertebrate animals other than mammals that have the ability to concentrate their urine. Previously, we identified and characterized aquaporins (AQP)1-4 responsible for urine concentration in Japanese quail kidneys. Today, a total of 13 orthologs for these genes have been reported in birds. Bird AQPs can be classified into four subfamilies: 1) Classical AQPs (AQP0-5 and novel member, AQP4-like) that conserve the selectivity filter; 2) aquaglyceroporins (AQP3, 7, 9 and 10) that retain an aspartic acid residue in the second NPA box and expand the pore to accept larger molecules; 3) unorthodox AQPs (AQP11-12) which structurally resemble their mammalian counterparts; 4) AQP8-type, a subfamily that differs from mammalian AQP8. Interestingly, over the course of time, birds lost their mammalian counterpart AQP6 but obtained a novel AQP4-like aquaporin member. In quail and/or chicken kidneys, at least six AQPs are expressed. Quail AQP1 (qAQP1) is expressed in both cortical and medullary proximal tubules but is absent in the descending limb (DL) and the thick ascending limb (TAL), supporting our previous finding that the DL and TAL are water impermeable. AQP2, an arginine vasotocin (AVT)-sensitive water channel, is exclusively expressed in the principal cells of the collecting duct (CD). AQP4 is unlikely to participate in free water resorption from the collecting duct (CD), and only AQP3 may represent an exit pathway for water reabsorbed apically via AQP2. While AQP9 is not expressed in mammalian kidneys, AQP9 was recently found in chicken kidneys. This review summarizes the current knowledge of the structure, function and expression of bird AQPs.

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.

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.

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.

Lamprey lecticans link new vertebrate genes to the origin and elaboration of vertebrate tissues

The evolution of vertebrates from an invertebrate chordate ancestor involved the evolution of new organs, tissues, and cell types. It was also marked by the origin and duplication of new gene families. If, and how, these morphological and genetic innovations are related is an unresolved question in vertebrate evolution. Hyaluronan is an extracellular matrix (ECM) polysaccharide important for water homeostasis and tissue structure. Vertebrates possess a novel family of hyaluronan binding proteins called Lecticans, and studies in jawed vertebrates (gnathostomes) have shown they function in many of the cells and tissues that are unique to vertebrates. This raises the possibility that the origin and/or expansion of this gene family helped drive the evolution of these vertebrate novelties. In order to better understand the evolution of the lectican gene family, and its role in the evolution of vertebrate morphological novelties, we investigated the phylogeny, genomic arrangement, and expression patterns of all lecticans in the sea lamprey (Petromyzon marinus), a jawless vertebrate. Though both P. marinus and gnathostomes each have four lecticans, our phylogenetic and syntenic analyses are most consistent with the independent duplication of one of more lecticans in the lamprey lineage. Despite the likely independent expansion of the lamprey and gnathostome lectican families, we find highly conserved expression of lecticans in vertebrate-specific and mesenchyme-derived tissues. We also find that, unlike gnathostomes, lamprey expresses its lectican paralogs in distinct subpopulations of head skeleton precursors, potentially reflecting an ancestral diversity of skeletal tissue types. Together, these observations suggest that the ancestral pre-duplication lectican had a complex expression pattern, functioned to support mesenchymal histology, and likely played a role in the evolution of vertebrate-specific cell and tissue types.

Evolution of thymopoietic microenvironments

In vertebrates, the development of lymphocytes from undifferentiated haematopoietic precursors takes place in so-called primary lymphoid organs, such as the thymus. Therein, lymphocytes undergo a complex differentiation and selection process that culminates in the generation of a pool of mature T cells that collectively express a self-tolerant repertoire of somatically diversified antigen receptors. Throughout this entire process, the microenvironment of the thymus in large parts dictates the sequence and outcome of the lymphopoietic activity. In vertebrates, direct genetic evidence in some species and circumstantial evidence in others suggest that the formation of a functional thymic microenvironment is controlled by members of the Foxn1/4 family of transcription factors. In teleost fishes, both Foxn1 and Foxn4 contribute to thymopoietic activity, whereas Foxn1 is both necessary and sufficient in the mammalian thymus. The evolutionary history of Foxn1/4 genes suggests that an ancient Foxn4 gene lineage gave rise to the Foxn1 genes in early vertebrates, raising the question of the thymopoietic capacity of the ancestor common to all vertebrates. Recent attempts to reconstruct the early events in the evolution of thymopoietic tissues by replacement of the mouse Foxn1 gene by Foxn1-like genes isolated from various chordate species suggest a plausible scenario. It appears that the primordial thymus was a bi-potent lymphoid organ, supporting both B cell and T cell development; however, during the course of vertebrate, evolution B cell development was gradually diminished converting the thymus into a site specialized in T cell development.

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.

The Evolution of Oxytocin and Vasotocin Receptor Genes in Jawed Vertebrates: A Clear Case for Gene Duplications Through Ancestral Whole-Genome Duplications

The neuronal and neuroendocrine peptides oxytocin (OT) and vasotocin (VT), including vasopressins, have six cognate receptors encoded by six receptor subtype genes in jawed vertebrates. The peptides elicit a broad range of responses that are specifically mediated by the receptor subtypes including neuronal functions regulating behavior and hormonal actions on reproduction and water/electrolyte balance. Previously, we have demonstrated that these six receptor subtype genes, which we designated VTR1A, VTR1B, OTR, VTR2A, VTR2B and VTR2C, arose from a syntenic ancestral gene pair, one VTR1/OTR ancestor and one VTR2 ancestor, through the early vertebrate whole-genome duplications (WGD) called 1R and 2R. This was supported by both phylogenetic and chromosomal conserved synteny data. More recently, other studies have focused on confounding factors, such as the OTR/VTR orthologs in cyclostomes, to question this scenario for the origin of the OTR/VTR gene family; proposing instead less parsimonious interpretations involving only one WGD followed by complex series of chromosomal or segmental duplications. Here, we have updated the phylogeny of the OTR/VTR gene family, including a larger number of vertebrate species, and revisited seven representative neighboring gene families from our previous conserved synteny analyses, adding chromosomal information from newer high-coverage genome assemblies from species that occupy key phylogenetic positions: the polypteriform fish reedfish (Erpetoichthys calabaricus), the cartilaginous fish thorny skate (Amblyraja radiata) and a more recent high-quality assembly of the Western clawed frog (Xenopus tropicalis) genome. Our analyses once again add strong support for four-fold symmetry, i.e., chromosome quadruplication in the same time window as the WGD events early in vertebrate evolution, prior to the jawed vertebrate radiation. Thus, the evolution of the OTR/VTR gene family can be most parsimoniously explained by two WGD events giving rise to the six ancestral genes, followed by differential gene losses of VTR2 genes in different lineages. We also argue for more coherence and clarity in the nomenclature of OT/VT receptors, based on the most parsimonious scenario.

Visual opsin expression and morphological characterization of retinal photoreceptors in the pouched lamprey (Geotria australis, Gray)

Lampreys are extant members of the agnathan (jawless) vertebrates that diverged ~500 million years ago, during a critical stage of vertebrate evolution when image-forming eyes first emerged. Among lamprey species assessed thus far, the retina of the southern hemisphere pouched lamprey, Geotria australis, is unique, in that it possesses morphologically distinct photoreceptors and expresses five visual photopigments. This study focused on determining the number of different photoreceptors present in the retina of G. australis and whether each cell type expresses a single opsin class. Five photoreceptor subtypes were identified based on ultrastructure and differential expression of one of each of the five different visual opsin classes (lws, sws1, sws2, rh1, and rh2) known to be expressed in the retina. This suggests, therefore, that the retina of G. australis possesses five spectrally and morphologically distinct photoreceptors, with the potential for complex color vision. Each photoreceptor subtype was shown to have a specific spatial distribution in the retina, which is potentially associated with changes in spectral radiance across different lines of sight. These results suggest that there have been strong selection pressures for G. australis to maintain broad spectral sensitivity for the brightly lit surface waters that this species inhabits during its marine phase. These findings provide important insights into the functional anatomy of the early vertebrate retina and the selection pressures that may have led to the evolution of complex color vision.

Evolution and expansion of the RUNX2 QA repeat corresponds with the emergence of vertebrate complexity

Runt-related transcription factor 2 (RUNX2) is critical for the development of the vertebrate bony skeleton. Unlike other RUNX family members, RUNX2 possesses a variable poly-glutamine, poly-alanine (QA) repeat domain. Natural variation within this repeat is able to alter the transactivation potential of RUNX2, acting as an evolutionary 'tuning knob' suggested to influence mammalian skull shape. However, the broader role of the RUNX2 QA repeat throughout vertebrate evolution is unknown. In this perspective, we examine the role of the RUNX2 QA repeat during skeletal development and discuss how its emergence and expansion may have facilitated the evolution of morphological novelty in vertebrates.

The immune system of jawless vertebrates: insights into the prototype of the adaptive immune system

Jawless vertebrates diverged from an ancestor of jawed vertebrates approximately 550 million years ago. They mount adaptive immune responses to repetitive antigenic challenges, despite lacking major histocompatibility complex molecules, immunoglobulins, T cell receptors, and recombination-activating genes. Instead of B cell and T cell receptors, agnathan lymphocytes express unique antigen receptors named variable lymphocyte receptors (VLRs), which generate diversity through a gene conversion-like mechanism. Although gnathostome antigen receptors and VLRs are structurally unrelated, jawed and jawless vertebrates share essential features of lymphocyte-based adaptive immunity, including the expression of a single type of receptor on each lymphocyte, clonal expansion of antigen-stimulated lymphocytes, and the dichotomy of cellular and humoral immunity, indicating that the backbone of the adaptive immune system was established in a common ancestor of all vertebrates. Furthermore, recent evidence indicates that, unlike previously thought, agnathans have a unique classical pathway of complement activation where VLRB molecules act as antibodies instead of immunoglobulins. It seems likely that the last common ancestor of all vertebrates had an adaptive immune system resembling that of jawless vertebrates, suggesting that, as opposed to jawed vertebrates, agnathans have retained the prototype of vertebrate adaptive immunity.

Tumor Necrosis Factor Superfamily: Ancestral Functions and Remodeling in Early Vertebrate Evolution

The evolution of the tumor necrosis factor superfamily (TNFSF) in early vertebrates is inferred by comparing the TNFSF genes found in humans and nine fishes: three agnathans, two chondrichthyans, three actinopterygians, and the sarcopterygian Latimeria chalumnae. By combining phylogenetic and synteny analyses, the TNFSF sequences detected are classified into five clusters of genes and 24 orthology groups. A model for their evolution since the origin of vertebrates is proposed. Fifteen TNFSF genes emerged from just three progenitors due to the whole-genome duplications (WGDs) that occurred before the agnathan/gnathostome split. Later, gnathostomes not only kept most of the genes emerged in the WGDs but soon added several tandem duplicates. More recently, complex, lineage-specific patterns of duplications and losses occurred in different gnathostome lineages. In agnathan species only seven to eight TNFSF genes are detected, because this lineage soon lost six of the genes emerged in the ancestral WGDs and additional losses in both hagfishes and lampreys later occurred. The orthologs of many of these lost genes are, in mammals, ligands of death-domain-containing TNFSF receptors, indicating that the extrinsic apoptotic pathway became simplified in the agnathan lineage. From the patterns of emergence of these genes, it is deduced that both the regulation of apoptosis and the control of the NF-κB pathway that depends in modern mammals on TNFSF members emerged before the ancestral vertebrate WGDs.

Polyploidization and pseudogenization in allotetraploid frog Xenopus laevis promote the evolution of aquaporin family in higher vertebrates

BACKGROUND

Aquaporins (AQPs), as members of the major intrinsic protein (MIP) superfamily, facilitated the permeation of water and other solutes and are involved in multiple biological processes. AQP family exists in almost all living organisms and is highly diversified in vertebrates in both classification and function due to genome wide duplication. While some AQP orthologs have been lost in higher vertebrates through evolution.

RESULT

Genome-wide comparative analyses of the AQP family between allotetraploid frog Xenopus laevis (Xla) and diploid frog Xenopus tropicalis (Xtr), based on the genome assemblies, revealed that the number of AQPs in Xla genome nearly doubled that in Xtr (32 vs. 19). Synteny analysis indicated that the distribution of the retained AQPs in Xla subgenomes (17 in Xla. L, the longer homeolog of Xla genome and 15 in Xla. S, the shorter homeolog of Xla genome) were highly symmetrical when compared with that in Xtr genome. Remarkably, two members in Xla. L and four members in Xla. S were lost through evolution. Blast analysis revealed that the lost AQPs in Xla are pseudogenized via either the deletion of some exons or some single nucleotide insertions or deletions that lead the reading frame shift. Additionally, comparative genomic analyses suggested that the orthologs of AQPs that with one copy absence in Xla are also prone to be lost in higher vertebrates.

CONCLUSION

This study revealed that polyploidization and subsequent pseudogenization and deletion in Xla genome promote the evolution of AQP family in higher vertebrates. Besides, our results would also contribute to understanding the evolution of AQP family.

Mechanistic insights into the evolution of the differential expression of tandemly arrayed cone opsin genes in zebrafish

The genome of many organisms contains several loci consisting of duplicated genes that are arrayed in tandem. The daughter genes produced by duplication typically exhibit differential expression patterns with each other or otherwise experience pseudogenization. Remarkably, opsin genes in fish are preserved after many duplications in different lineages. This fact indicates that fish opsin genes are characterized by a regulatory mechanism that could intrinsically facilitate the differentiation of the expression patterns. However, little is known about the mechanisms that underlie the differential expression patterns or how they were established during evolution. The loci of green (RH2)- and red (LWS)-sensitive cone opsin genes in zebrafish have been used as model systems to study the differential regulation of tandemly arrayed opsin genes. Over a decade of studies have uncovered several mechanistic features that might have assisted the differentiation and preservation of duplicated genes. Furthermore, recent progress in the understanding of the transcriptional process in general has added essential insights. In this article, the current understanding of the transcriptional regulation of differentially expressed tandemly arrayed cone opsin genes in zebrafish is summarized and a possible evolutionary scenario that could achieve this differentiation is discussed.

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.

Evolutionary and Molecular Characterization of liver-enriched gene 1

Liver-enriched gene 1 (Leg1) is a newly identified gene with little available functional information. To evolutionarily and molecularly characterize Leg1 genes, a phylogenetic study was first conducted, which indicated that Leg1 is a conserved gene that exists from bacteria to mammals. During the evolution of mammals, Leg1s underwent tandem duplications, which gave rise to Leg1a, Leg1b, and Leg1c clades. Analysis of the pig genome showed the presence of all three paralogs of pig Leg1 genes (pLeg1s), whereas only Leg1a could be found in the human (hLeg1a) or mouse (mLeg1a) genomes. Purifying force acts on the evolution of Leg1 genes, likely subjecting them to functional constraint. Molecularly, pLeg1a and its coded protein, pig LEG1a (pLEG1a), displayed high similarities to its human and mouse homologs in terms of gene organization, expression patterns, and structures. Hence, pLeg1a, hLeg1a, and mLeg1a might preserve similar functions. Additionally, expression analysis of the three Leg1as suggested that eutherian Leg1as might have different functions from those of zebrafish and platypus due to subfunctionalization. Therefore, pLeg1a might provide essential information about eutherian Leg1a. Moreover, a preliminary functional study using RNA-seq suggested that pLeg1a is involved in the lipid homeostasis. In conclusion, our study provides some basic information on the aspects of evolution and molecular function, which could be applied for further validation of Leg1 using pig models.

A Comparison of the Innate and Adaptive Immune Systems in Cartilaginous Fish, Ray-Finned Fish, and Lobe-Finned Fish

The immune system is composed of two subsystems-the innate immune system and the adaptive immune system. The innate immune system is the first to respond to pathogens and does not retain memory of previous responses. Innate immune responses are evolutionarily older than adaptive responses and elements of innate immunity can be found in all multicellular organisms. If a pathogen persists, the adaptive immune system will engage the pathogen with specificity and memory. Several components of the adaptive system including immunoglobulins (Igs), T cell receptors (TCR), and major histocompatibility complex (MHC), are assumed to have arisen in the first jawed vertebrates-the Gnathostomata. This review will discuss and compare components of both the innate and adaptive immune systems in Gnathostomes, particularly in Chondrichthyes (cartilaginous fish) and in Osteichthyes [bony fish: the Actinopterygii (ray-finned fish) and the Sarcopterygii (lobe-finned fish)]. While many elements of both the innate and adaptive immune systems are conserved within these species and with higher level vertebrates, some elements have marked differences. Components of the innate immune system covered here include physical barriers, such as the skin and gastrointestinal tract, cellular components, such as pattern recognition receptors and immune cells including macrophages and neutrophils, and humoral components, such as the complement system. Components of the adaptive system covered include the fundamental cells and molecules of adaptive immunity: B lymphocytes (B cells), T lymphocytes (T cells), immunoglobulins (Igs), and major histocompatibility complex (MHC). Comparative studies in fish such as those discussed here are essential for developing a comprehensive understanding of the evolution of the immune system.

Phylogenetic Reclassification of Vertebrate Melatonin Receptors To Include Mel1d

The circadian and seasonal actions of melatonin are mediated by high affinity G-protein coupled receptors (melatonin receptors, MTRs), classified into phylogenetically distinct subtypes based on sequence divergence and pharmacological characteristics. Three vertebrate MTR subtypes are currently described: MT1 (MTNR1A), MT2 (MTNR1B), and Mel1c (MTNR1C / GPR50), which exhibit distinct affinities, tissue distributions and signaling properties. We present phylogenetic and comparative genomic analyses supporting a revised classification of the vertebrate MTR family. We demonstrate four ancestral vertebrate MTRs, including a novel molecule hereafter named Mel1d. We reconstructed the evolution of each vertebrate MTR, detailing genetic losses in addition to gains resulting from whole genome duplication events in teleost fishes. We show that Mel1d was lost separately in mammals and birds and has been previously mistaken for an MT1 paralogue. The genetic and functional diversity of vertebrate MTRs is more complex than appreciated, with implications for our understanding of melatonin actions in different taxa. The significance of our findings, including the existence of Mel1d, are discussed in an evolutionary and functional context accommodating a robust phylogenetic assignment of MTR gene family structure.

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.

Inference of the ancestral vertebrate phenotype through vestiges of the whole-genome duplications

Inferring the phenotype of the last common ancestor of living vertebrates is a challenging problem because of several unresolvable factors. They include the lack of reliable out-groups of living vertebrates, poor information about less fossilizable organs and specialized traits of phylogenetically important species, such as lampreys and hagfishes (e.g. secondary loss of vertebrae in adult hagfishes). These factors undermine the reliability of ancestral reconstruction by traditional character mapping approaches based on maximum parsimony. In this article, we formulate an approach to hypothesizing ancestral vertebrate phenotypes using information from the phylogenetic and functional properties of genes duplicated by genome expansions in early vertebrate evolution. We named the conjecture as ‘chronological reconstruction of ohnolog functions (CHROF)’. This CHROF conjecture raises the possibility that the last common ancestor of living vertebrates may have had more complex traits than currently thought.

Cartography of rhodopsin-like G protein-coupled receptors across vertebrate genomes

We conduct a cartography of rhodopsin-like non-olfactory G protein-coupled receptors in the Ensembl database. The most recent genomic data (releases 90–92, 90 vertebrate genomes) are analyzed through the online interface and receptors mapped on phylogenetic guide trees that were constructed based on a set of ~14.000 amino acid sequences. This snapshot of genomic data suggest vertebrate genomes to harbour 142 clades of GPCRs without human orthologues. Among those, 69 have not to our knowledge been mentioned or studied previously in the literature, of which 28 are distant from existing receptors and likely new orphans. These newly identified receptors are candidates for more focused evolutionary studies such as chromosomal mapping as well for in-depth pharmacological characterization. Interestingly, we also show that 37 of the 72 human orphan (or recently deorphanized) receptors included in this study cluster into nineteen closely related groups, which implies that there are less ligands to be identified than previously anticipated. Altogether, this work has significant implications when discussing nomenclature issues for GPCRs.

Asymmetric paralog evolution between the "cryptic" gene Bmp16 and its well-studied sister genes Bmp2 and Bmp4

The vertebrate gene repertoire is characterized by “cryptic” genes whose identification has been hampered by their absence from the genomes of well-studied species. One example is the Bmp16 gene, a paralog of the developmental key genes Bmp2 and -4. We focus on the Bmp2/4/16 group of genes to study the evolutionary dynamics following gen(om)e duplications with special emphasis on the poorly studied Bmp16 gene. We reveal the presence of Bmp16 in chondrichthyans in addition to previously reported teleost fishes and reptiles. Using comprehensive, vertebrate-wide gene sampling, our phylogenetic analysis complemented with synteny analyses suggests that Bmp2, -4 and -16 are remnants of a gene quartet that originated during the two rounds of whole-genome duplication (2R-WGD) early in vertebrate evolution. We confirm that Bmp16 genes were lost independently in at least three lineages (mammals, archelosaurs and amphibians) and report that they have elevated rates of sequence evolution. This finding agrees with their more “flexible” deployment during development; while Bmp16 has limited embryonic expression domains in the cloudy catshark, it is broadly expressed in the green anole lizard. Our study illustrates the dynamics of gene family evolution by integrating insights from sequence diversification, gene repertoire changes, and shuffling of expression domains.

Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates

Modern cartilaginous fishes are divided into elasmobranchs (sharks, rays and skates) and chimaeras, and the lack of established whole-genome sequences for the former has prevented our understanding of early vertebrate evolution and the unique phenotypes of elasmobranchs. Here we present de novo whole-genome assemblies of brownbanded bamboo shark and cloudy catshark and an improved assembly of the whale shark genome. These relatively large genomes (3.8-6.7 Gbp) contain sparse distributions of coding genes and regulatory elements and exhibit reduced molecular evolutionary rates. Our thorough genome annotation revealed Hox C genes previously hypothesized to have been lost, as well as distinct gene repertories of opsins and olfactory receptors that would be associated with adaptation to unique underwater niches. We also show the early establishment of the genetic machinery governing mammalian homoeostasis and reproduction at the jawed vertebrate ancestor. This study, supported by genomic, transcriptomic and epigenomic resources, provides a foundation for the comprehensive, molecular exploration of phenotypes unique to sharks and insights into the evolutionary origins of vertebrates.