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.

A Comparative Transcriptomic Study and Key Gene Targeting of Lamprey Gonadal Immune Response

The mammalian testis and ovary possess special immunocompetence, which is central to provide protection against pathogens. However, the innate immune responses to immune challenges in lamprey gonads are poorly understood. In this study, we extracted RNA from testis and ovary tissues of lampreys at 0 hour, 8 hours and 17 days after lipopolysaccharides (LPS) stimulation and performed transcriptome sequencing. While the transcriptome profiles of the two tissues were different for the most part, genes LIP, LECT2, LAL2, GRN, ITLN, and C1q were found to be the most significantly up-regulated genes in both. Quantitative Real-time PCR (qRT-PCR) analysis confirmed that these genes were upregulated after stimulation. Furthermore, immunohistochemical staining showed that these genes in lamprey gonads are expressed in high quantities and have a specific distribution. Taken together, our results suggest that these genes could play an essential role in response of the gonads to LPS induction. This research establishes a basis for investigating the immune mechanism of vertebrate gonads and presents a fresh concept for gaining insight into the evolutionary development of jawless vertebrates.

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.

Phylogeny and evolution of the genus Cervus (Cervidae, Mammalia) as revealed by complete mitochondrial genomes

Mitochondrial DNA (mtDNA) lineages are recognized as important components of intra- and interspecific biodiversity, and allow to reveal colonization routes and phylogeographic structure of many taxa. Among these is the genus Cervus that is widely distributed across the Holarctic. We obtained sequences of complete mitochondrial genomes from 13 Cervus taxa and included them in global phylogenetic analyses of 71 Cervinae mitogenomes. The well-resolved phylogenetic trees confirmed Cervus to be monophyletic. Molecular dating based on several fossil calibration points revealed that ca. 2.6 Mya two main mitochondrial lineages of Cervus separated in Central Asia, the Western (including C. hanglu and C. elaphus) and the Eastern (comprising C. albirostris, C. canadensis and C. nippon). We also observed convergent changes in the composition of some mitochondrial genes in C. hanglu of the Western lineage and representatives of the Eastern lineage. Several subspecies of C. nippon and C. hanglu have accumulated a large portion of deleterious substitutions in their mitochondrial protein-coding genes, probably due to drift in the wake of decreasing population size. In contrast to previous studies, we found that the relic haplogroup B of C. elaphus was sister to all other red deer lineages and that the Middle-Eastern haplogroup E shared a common ancestor with the Balkan haplogroup C. Comparison of the mtDNA phylogenetic tree with a published nuclear genome tree may imply ancient introgressions of mtDNA between different Cervus species as well as from the common ancestor of South Asian deer, Rusa timorensis and R. unicolor, to the Cervus clade.

Biliary Atresia Animal Models: Is the Needle in a Haystack?

Biliary atresia (BA) is a progressive fibro-obliterative process with a variable degree of inflammation involving the hepatobiliary system. Its consequences are incalculable for the patients, the affected families, relatives, and the healthcare system. Scientific communities have identified a rate of about 1 case per 10,000-20,000 live births, but the percentage may be higher, considering the late diagnoses. The etiology is heterogeneous. BA, which is considered in half of the causes leading to orthotopic liver transplantation, occurs in primates and non-primates. To consolidate any model, (1) more transport and cell membrane studies are needed to identify the exact mechanism of noxa-related hepatotoxicity; (2) an online platform may be key to share data from pilot projects and new techniques; and (3) the introduction of differentially expressed genes may be useful in investigating the liver metabolism to target the most intricate bilio-toxic effects of pharmaceutical drugs and toxins. As a challenge, such methodologies are still limited to very few centers, making the identification of highly functional animal models like finding a "needle in a haystack". This review compiles models from the haystack and hopes that a combinatorial search will eventually be the root for a successful pathway.

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.

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.

New view on the organization and evolution of Palaeognathae mitogenomes poses the question on the ancestral gene rearrangement in Aves

BACKGROUND

Bird mitogenomes differ from other vertebrates in gene rearrangement. The most common avian gene order, identified first in Gallus gallus, is considered ancestral for all Aves. However, other rearrangements including a duplicated control region and neighboring genes have been reported in many representatives of avian orders. The repeated regions can be easily overlooked due to inappropriate DNA amplification or genome sequencing. This raises a question about the actual prevalence of mitogenomic duplications and the validity of the current view on the avian mitogenome evolution. In this context, Palaeognathae is especially interesting because is sister to all other living birds, i.e. Neognathae. So far, a unique duplicated region has been found in one palaeognath mitogenome, that of Eudromia elegans.

RESULTS

Therefore, we applied an appropriate PCR strategy to look for omitted duplications in other palaeognaths. The analyses revealed the duplicated control regions with adjacent genes in Crypturellus, Rhea and Struthio as well as ND6 pseudogene in three moas. The copies are very similar and were subjected to concerted evolution. Mapping the presence and absence of duplication onto the Palaeognathae phylogeny indicates that the duplication was an ancestral state for this avian group. This feature was inherited by early diverged lineages and lost two times in others. Comparison of incongruent phylogenetic trees based on mitochondrial and nuclear sequences showed that two variants of mitogenomes could exist in the evolution of palaeognaths. Data collected for other avian mitogenomes revealed that the last common ancestor of all birds and early diverging lineages of Neoaves could also possess the mitogenomic duplication.

CONCLUSIONS

The duplicated control regions with adjacent genes are more common in avian mitochondrial genomes than it was previously thought. These two regions could increase effectiveness of replication and transcription as well as the number of replicating mitogenomes per organelle. In consequence, energy production by mitochondria may be also more efficient. However, further physiological and molecular analyses are necessary to assess the potential selective advantages of the mitogenome duplications.

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

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

The Zebrafish Amygdaloid Complex - Functional Ground Plan, Molecular Delineation, and Everted Topology

In mammals and other tetrapods, a multinuclear forebrain structure, called the amygdala, forms the neuroregulatory core essential for emotion, cognition, and social behavior. Currently, higher circuits of affective behavior in anamniote non-tetrapod vertebrates (“fishes”) are poorly understood, preventing a comprehensive understanding of amygdala evolution. Through molecular characterization and evolutionary-developmental considerations, we delineated the complex amygdala ground plan of zebrafish, whose everted telencephalon has made comparisons to the evaginated forebrains of tetrapods challenging. In this radical paradigm, thirteen telencephalic territories constitute the zebrafish amygdaloid complex and each territory is distinguished by conserved molecular properties and structure-functional relationships with other amygdaloid structures. Central to our paradigm, the study identifies the teleostean amygdaloid nucleus of the lateral olfactory tract (nLOT), an olfactory integrative structure that links dopaminergic telencephalic groups to the amygdala alongside redefining the putative zebrafish olfactory pallium (“Dp”). Molecular characteristics such as the distribution of substance P and the calcium-binding proteins parvalbumin (PV) and calretinin (CR) indicate, that the zebrafish extended centromedial (autonomic and reproductive) amygdala is predominantly located in the GABAergic and isl1-negative territory. Like in tetrapods, medial amygdaloid (MeA) nuclei are defined by the presence of substance P immunoreactive fibers and calretinin-positive neurons, whereas central amygdaloid (CeA) nuclei lack these characteristics. A detailed comparison of lhx5-driven and vGLut2a-driven GFP in transgenic reporter lines revealed ancestral topological relationships between the thalamic eminence (EmT), the medial amygdala (MeA), the nLOT, and the integrative olfactory pallium. Thus, the study explains how the zebrafish amygdala and the complexly everted telencephalon topologically relate to the corresponding structures in mammals indicating that an elaborate amygdala ground plan evolved early in vertebrates, in a common ancestor of teleosts and tetrapods.

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.

Deeply conserved synteny resolves early events in vertebrate evolution

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

Resolving Phylogenetic Relationships within Passeriformes Based on Mitochondrial Genes and Inferring the Evolution of Their Mitogenomes in Terms of Duplications

Mitochondrial genes are placed on one molecule, which implies that they should carry consistent phylogenetic information. Following this advantage, we present a well-supported phylogeny based on mitochondrial genomes from almost 300 representatives of Passeriformes, the most numerous and differentiated Aves order. The analyses resolved the phylogenetic position of paraphyletic Basal and Transitional Oscines. Passerida occurred divided into two groups, one containing Paroidea and Sylvioidea, whereas the other, Passeroidea and Muscicapoidea. Analyses of mitogenomes showed four types of rearrangements including a duplicated control region (CR) with adjacent genes. Mapping the presence and absence of duplications onto the phylogenetic tree revealed that the duplication was the ancestral state for passerines and was maintained in early diverged lineages. Next, the duplication could be lost and occurred independently at least four times according to the most parsimonious scenario. In some lineages, two CR copies have been inherited from an ancient duplication and highly diverged, whereas in others, the second copy became similar to the first one due to concerted evolution. The second CR copies accumulated over twice as many substitutions as the first ones. However, the second CRs were not completely eliminated and were retained for a long time, which suggests that both regions can fulfill an important role in mitogenomes. Phylogenetic analyses based on CR sequences subjected to the complex evolution can produce tree topologies inconsistent with real evolutionary relationships between species. Passerines with two CRs showed a higher metabolic rate in relation to their body mass.

Discovery of All Three Types in Cartilaginous Fishes Enables Phylogenetic Resolution of the Origins and Evolution of Interferons

Interferons orchestrate host antiviral responses in jawed vertebrates. They are categorized into three classes; IFN1 and IFN3 are the primary antiviral cytokine lineages, while IFN2 responds to a broader variety of pathogens. The evolutionary relationships within and between these three classes have proven difficult to resolve. Here, we reassess interferon evolution, considering key phylogenetic pitfalls including taxon sampling, alignment quality, model adequacy, and outgroup choice. We reveal that cartilaginous fishes, and hence the jawed vertebrate ancestor, possess(ed) orthologs of all three interferon classes. We show that IFN3 groups sister to IFN1, resolve the origins of the human IFN3 lineages, and find that intronless IFN3s emerged at least three times. IFN2 genes are highly conserved, except for IFN-γ-rel, which we confirm resulted from a teleost-specific duplication. Our analyses show that IFN1 phylogeny is highly sensitive to phylogenetic error. By accounting for this, we describe a new backbone IFN1 phylogeny that implies several IFN1 genes existed in the jawed vertebrate ancestor. One of these is represented by the intronless IFN1s of tetrapods, including mammalian-like repertoires of reptile IFN1s and a subset of amphibian IFN1s, in addition to newly-identified intron-containing shark IFN1 genes. IFN-f, previously only found in teleosts, likely represents another ancestral jawed vertebrate IFN1 family member, suggesting the current classification of fish IFN1s into two groups based on the number of cysteines may need revision. The providence of the remaining fish IFN1s and the coelacanth IFN1s proved difficult to resolve, but they may also be ancestral jawed vertebrate IFN1 lineages. Finally, a large group of amphibian-specific IFN1s falls sister to all other IFN1s and was likely also present in the jawed vertebrate ancestor. Our results verify that intronless IFN1s have evolved multiple times in amphibians and indicate that no one-to-one orthology exists between mammal and reptile IFN1s. Our data also imply that diversification of the multiple IFN1s present in the jawed vertebrate ancestor has occurred through a rapid birth-death process, consistent with functional maintenance over a 450-million-year host-pathogen arms race. In summary, this study reveals a new model of interferon evolution important to our understanding of jawed vertebrate antiviral immunity.

Evolution of nodal and nodal-related genes and the putative composition of the heterodimers that trigger the nodal pathway in vertebrates

Nodal is a signaling molecule that belongs to the transforming growth factor-β superfamily that plays key roles during the early stages of development of animals. In vertebrates Nodal forms an heterodimer with a GDF1/3 protein to activate the Nodal pathway. Vertebrates have a paralog of nodal in their genomes labeled Nodal-related, but the evolutionary history of these genes is a matter of debate, mainly because of the presence of a variable numbers of genes in the vertebrate genomes sequenced so far. Thus, the goal of this study was to investigate the evolutionary history of the Nodal and Nodal-related genes with an emphasis in tracking changes in the number of genes among vertebrates. Our results show the presence of two gene lineages (Nodal and Nodal-related) that can be traced back to the ancestor of jawed vertebrates. These lineages have undergone processes of differential retention and lineage-specific expansions. Our results imply that Nodal and Nodal-related duplicated at the latest in the ancestor of gnathostomes, and they still retain a significant level of functional redundancy. By comparing the evolution of the Nodal/Nodal-related with GDF1/3 gene family, it is possible to infer that there are several types of heterodimers that can trigger the Nodal pathway among vertebrates.

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.

Duplicated Myosin V Genes in Teleosts Show Evolutionary Rate Variations among the Motor and Cargo-Binding Domains

We analyzed evolutionary rates of conserved, duplicated myosin V (myo5) genes in nine teleost species to examine the outcomes of duplication events. Syntenic analysis and ancestral chromosome mapping suggest one tandem gene duplication event leading to the appearance of myo5a and myo5c, two rounds of whole genome duplication for vertebrates, and an additional round of whole genome duplication for teleosts account for the presence and location of the myo5 genes and their duplicates in teleosts and other vertebrates and the timing of the duplication events. Phylogenetic analyses reveal a previously unidentified myo5 clade that we refer to now as myo5bb. Analysis using dN/dS rate comparisons revealed large regions within duplicated myo5 genes that are highly conserved. Codons identified in other studies as encoding functionally important portions of the Myo5a and Myo5b proteins are shown to be highly conserved within the newly identified myo5bb clade and in other myo5 duplicates. As much as 30% of 319 codons encoding the cargo-binding domain in the myo5aa genes are conserved in all three codon positions in nine teleost species. For the myo5bb cargo-binding domain, 6.6% of 336 codons have zero substitutions in all nine teleost species. Using molecular evolution assays, we identify the myo5bb branch as being subject to evolutionary rate variation with the cargo-binding domain, having 20% of the sites under positive selection and the motor domain having 8% of its sites under positive selection. The high number of invariant codons coupled with relatively high dN/dS values in the region of the myo5 genes encoding the ATP-binding domain suggests the encoded proteins retain function and may have acquired novel functions associated with changes to the cargo-binding domain.

Teleost Fish-Specific Preferential Retention of Pigmentation Gene-Containing Families After Whole Genome Duplications in Vertebrates

Vertebrate pigmentation is a highly diverse trait mainly determined by neural crest cell derivatives. It has been suggested that two rounds (1R/2R) of whole-genome duplications (WGDs) at the basis of vertebrates allowed changes in gene regulation associated with neural crest evolution. Subsequently, the teleost fish lineage experienced other WGDs, including the teleost-specific Ts3R before teleost radiation and the more recent Ss4R at the basis of salmonids. As the teleost lineage harbors the highest number of pigment cell types and pigmentation diversity in vertebrates, WGDs might have contributed to the evolution and diversification of the pigmentation gene repertoire in teleosts. We have compared the impact of the basal vertebrate 1R/2R duplications with that of the teleost-specific Ts3R and salmonid-specific Ss4R WGDs on 181 gene families containing genes involved in pigmentation. We show that pigmentation genes (PGs) have been globally more frequently retained as duplicates than other genes after Ts3R and Ss4R but not after the early 1R/2R. This is also true for non-pigmentary paralogs of PGs, suggesting that the function in pigmentation is not the sole key driver of gene retention after WGDs. On the long-term, specific categories of PGs have been repeatedly preferentially retained after ancient 1R/2R and Ts3R WGDs, possibly linked to the molecular nature of their proteins (e.g., DNA binding transcriptional regulators) and their central position in protein-protein interaction networks. Taken together, our results support a major role of WGDs in the diversification of the pigmentation gene repertoire in the teleost lineage, with a possible link with the diversity of pigment cell lineages observed in these animals compared to other vertebrates.

The influence of molecular markers and methods on inferring the phylogenetic relationships between the representatives of the Arini (parrots, Psittaciformes), determined on the basis of their complete mitochondrial genomes

BACKGROUND

Conures are a morphologically diverse group of Neotropical parrots classified as members of the tribe Arini, which has recently been subjected to a taxonomic revision. The previously broadly defined Aratinga genus of this tribe has been split into the 'true' Aratinga and three additional genera, Eupsittula, Psittacara and Thectocercus. Popular markers used in the reconstruction of the parrots' phylogenies derive from mitochondrial DNA. However, current phylogenetic analyses seem to indicate conflicting relationships between Aratinga and other conures, and also among other Arini members. Therefore, it is not clear if the mtDNA phylogenies can reliably define the species tree. The inconsistencies may result from the variable evolution rate of the markers used or their weak phylogenetic signal. To resolve these controversies and to assess to what extent the phylogenetic relationships in the tribe Arini can be inferred from mitochondrial genomes, we compared representative Arini mitogenomes as well as examined the usefulness of the individual mitochondrial markers and the efficiency of various phylogenetic methods.

RESULTS

Single molecular markers produced inconsistent tree topologies, while different methods offered various topologies even for the same marker. A significant disagreement in these tree topologies occurred for cytb, nd2 and nd6 genes, which are commonly used in parrot phylogenies. The strongest phylogenetic signal was found in the control region and RNA genes. However, these markers cannot be used alone in inferring Arini phylogenies because they do not provide fully resolved trees. The most reliable phylogeny of the parrots under study is obtained only on the concatenated set of all mitochondrial markers. The analyses established significantly resolved relationships within the former Aratinga representatives and the main genera of the tribe Arini. Such mtDNA phylogeny can be in agreement with the species tree, owing to its match with synapomorphic features in plumage colouration.

CONCLUSIONS

Phylogenetic relationships inferred from single mitochondrial markers can be incorrect and contradictory. Therefore, such phylogenies should be considered with caution. Reliable results can be produced by concatenated sets of all or at least the majority of mitochondrial genes and the control region. The results advance a new view on the relationships among the main genera of Arini and resolve the inconsistencies between the taxa that were previously classified as the broadly defined genus Aratinga. Although gene and species trees do not always have to be consistent, the mtDNA phylogenies for Arini can reflect the species tree.

An evaluation of the evolution of the gene structure of dystroglycan

Background

Dystroglycan (DG) is an adhesion receptor complex composed of two non-covalently associated subunits, transcribed from a single gene. The extracellular α-DG is highly and heterogeneously glycosylated and binds with high affinity to laminins, and the transmembrane β-DG binds intracellular dystrophin. Multiple cellular functions have been proposed for DG, notwithstanding that its role in skeletal muscle appears central as demonstrated by both primary and secondary severe muscular dystrophic phenotypes collectively known as dystroglycanopathies. We recently analysed the molecular phylogeny of the DG core protein and identified the α/β interface, transmembrane and cytoplasmic domains of β-DG as the most conserved region. It was also identified that the IG2_MAT_NU region has been independently duplicated in multiple lineages.

Results

To understand the evolution of dystroglycan in more depth, we investigated dystroglycan gene structure in 35 species representative of the phyla in which dystroglycan has been identified (i.e., all metazoan phyla except Ctenophora). The gene structure of three exons and two introns is remarkably conserved. However, additional lineage-specific introns were identified, which interrupt the coding sequence at distinct points, were identified in multiple metazoan groups, most prominently in ecdysozoans.

Conclusions

A coding DNA sequence (CDS) intron that interrupts the encoding of the IG1 domain is universally conserved and this intron is longer in gnathostomes (jawed vertebrates) than in other metazoans. Lineage-specific gain of additional introns has occurred notably in ecdysozoans, where multiple introns interrupt the large 3′ exon. More limited intron gain has also occurred in placozoa, cnidarians, urochordates and the DG paralogues of lamprey and teleost fish.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-016-2322-x) contains supplementary material, which is available to authorized users.

Inferring Orthologs: Open Questions and Perspectives

With the increasing number of sequenced genomes and their comparisons, the detection of orthologs is crucial for reliable functional annotation and evolutionary analyses of genes and species. Yet, the dynamic remodeling of genome content through gain, loss, transfer of genes, and segmental and whole-genome duplication hinders reliable orthology detection. Moreover, the lack of direct functional evidence and the questionable quality of some available genome sequences and annotations present additional difficulties to assess orthology. This article reviews the existing computational methods and their potential accuracy in the high-throughput era of genome sequencing and anticipates open questions in terms of methodology, reliability, and computation. Appropriate taxon sampling together with combination of methods based on similarity, phylogeny, synteny, and evolutionary knowledge that may help detecting speciation events appears to be the most accurate strategy. This review also raises perspectives on the potential determination of orthology throughout the whole species phylogeny.

Comparative analysis of zebrafish bone morphogenetic proteins 2, 4 and 16: molecular and evolutionary perspectives

BMP2, BMP4 and BMP16 form a subfamily of bone morphogenetic proteins acting as pleiotropic growth factors during development and as bone inducers during osteogenesis. BMP16 is the most recent member of this subfamily and basic data regarding protein structure and function, and spatio-temporal gene expression is still scarce. In this work, insights on BMP16 were provided through the comparative analysis of structural and functional data for zebrafish BMP2a, BMP2b, BMP4 and BMP16 genes and proteins, determined from three-dimensional models, patterns of gene expression during development and in adult tissues, regulation by retinoic acid and capacity to activate BMP-signaling pathway. Structures of Bmp2a, Bmp2b, Bmp4 and Bmp16 were found to be remarkably similar; with residues involved in receptor binding being highly conserved. All proteins could activate the BMP-signaling pathway, suggesting that they share a common function. On the contrary, stage- and tissue-specific expression of bmp2, bmp4 and bmp16 suggested the genes might be differentially regulated (e.g. different transcription factors, enhancers and/or regulatory modules) but also that they are involved in distinct physiological processes, although with the same function. Retinoic acid, a morphogen known to interact with BMP-signaling during bone formation, was shown to down-regulate the expression of bmp2, bmp4 and bmp16, although to different extents. Taxonomic and phylogenetic analyses indicated that bmp16 diverged before bmp2 and bmp4, is not restricted to teleost fish lineage as previously reported, and that it probably arose from a whole genomic duplication event that occurred early in vertebrate evolution and disappeared in various tetrapod lineages through independent events.

The emergence of the vasopressin and oxytocin hormone receptor gene family lineage: Clues from the characterization of vasotocin receptors in the sea lamprey (Petromyzon marinus)

The sea lamprey (Petromyzon marinus) is a jawless vertebrate at an evolutionary nexus between invertebrates and jawed vertebrates. Lampreys are known to possess the arginine vasotocin (AVT) hormone utilized by all non-mammalian vertebrates. We postulated that the lamprey would possess AVT receptor orthologs of predecessors to the arginine vasopressin (AVP)/oxytocin (OXT) family of G protein-coupled receptors found in mammals, providing insights into the origins of the mammalian V1A, V1B, V2 and OXT receptors. Among the earliest animals to diverge from the vertebrate lineage in which these receptors are characterized is the jawed, cartilaginous elephant shark, which has genes orthologous to all four mammalian receptor types. Therefore, our work was aimed at helping resolve the critical gap concerning the outcomes of hypothesized large-scale (whole-genome) duplication events. We sequenced one partial and four full-length putative lamprey AVT receptor genes and determined their mRNA expression patterns in 15 distinct tissues. Phylogenetically, three of the full-coding genes possess structural characteristics of the V1 clade containing the V1A, V1B and OXT receptors. Another full-length coding gene and the partial sequence are part of the V2 clade and appear to be most closely related to the newly established V2B and V2C receptor subtypes. Our synteny analysis also utilizing the Japanese lamprey (Lethenteron japonicum) genome supports the recent proposal that jawless and jawed vertebrates shared one-round (1R) of WGD as the most likely scenario.

Incorporating tree-thinking and evolutionary time scale into developmental biology

Phylogenetic approaches are indispensable in any comparative molecular study involving multiple species. These approaches are in increasing demand as the amount and availability of DNA sequence information continues to increase exponentially, even for organisms that were previously not extensively studied. Without the sound application of phylogenetic concepts and knowledge, one can be misled when attempting to infer ancestral character states as well as the timing and order of evolutionary events, both of which are frequently exerted in evolutionary developmental biology. The ignorance of phylogenetic approaches can also impact non-evolutionary studies and cause misidentification of the target gene or protein to be examined in functional characterization. This review aims to promote tree-thinking in evolutionary conjecture and stress the importance of a sense of time scale in cross-species comparisons, in order to enhance the understanding of phylogenetics in all biological fields including developmental biology. To this end, molecular phylogenies of several developmental regulatory genes, including those denoted as "cryptic pan-vertebrate genes", are introduced as examples.

Neuropeptide Y family receptors Y1 and Y2 from sea lamprey, Petromyzon marinus

The vertebrate gene family for neuropeptide Y (NPY) receptors expanded by duplication of the chromosome carrying the ancestral Y1-Y2-Y5 gene triplet. After loss of some duplicates, the ancestral jawed vertebrate had seven receptor subtypes forming the Y1 (including Y1, Y4, Y6, Y8), Y2 (including Y2, Y7) and Y5 (only Y5) subfamilies. Lampreys are considered to have experienced the same chromosome duplications as gnathostomes and should also be expected to have multiple receptor genes. However, previously only a Y4-like and a Y5 receptor have been cloned and characterized. Here we report the cloning and characterization of two additional receptors from the sea lamprey Petromyzon marinus. Sequence phylogeny alone could not with certainty assign their identity, but based on synteny comparisons of P. marinus and the Arctic lamprey, Lethenteron camtschaticum, with jawed vertebrates, the two receptors most likely are Y1 and Y2. Both receptors were expressed in human HEK293 cells and inositol phosphate assays were performed to determine the response to the three native lamprey peptides NPY, PYY and PMY. The three peptides have similar potencies in the nanomolar range for Y1. No obvious response to the three peptides was detected for Y2. Synteny analysis supports identification of the previously cloned receptor as Y4. No additional NPY receptor genes could be identified in the presently available lamprey genome assemblies. Thus, four NPY-family receptors have been identified in lampreys, orthologs of the same subtypes as in humans (Y1, Y2, Y4 and Y5), whereas many other vertebrate lineages have retained additional ancestral subtypes.

Split-alignment of genomes finds orthologies more accurately

We present a new pair-wise genome alignment method, based on a simple concept of finding an optimal set of local alignments. It gains accuracy by not masking repeats, and by using a statistical model to quantify the (un)ambiguity of each alignment part. Compared to previous animal genome alignments, it aligns thousands of locations differently and with much higher similarity, strongly suggesting that the previous alignments are non-orthologous. The previous methods suffer from an overly-strong assumption of long un-rearranged blocks. The new alignments should help find interesting and unusual features, such as fast-evolving elements and micro-rearrangements, which are confounded by alignment errors.

Reevaluating Emx gene phylogeny: homopolymeric amino acid tracts as a potential factor obscuring orthology signals in cyclostome genes

BACKGROUND

Vertebrate Emx genes, retained as multiple copies, are expressed in a nested pattern in the early embryonic forebrain and required for its regionalization. This pattern seems to have originated in a vertebrate common ancestor; however, a previous analysis, reporting two lamprey Emx genes, claimed independent Emx gene duplications in both cyclostome (extant jawless fish) and gnathostome (jawed vertebrate) lineages after their divergence. This scenario is neither parsimonious nor consistent with the hypothesis that genome expansion occurred before the cyclostome-gnathostome split, which is supported by recent genome-wide analyses.

RESULTS

We isolated and sequenced cDNA of two hagfish Emx genes and performed intensive molecular phylogenetic analyses, including the hagfish and/or lamprey Emx genes. The lamprey genes tended to attract each other in inferred phylogenetic trees, an effect that tended to be relaxed on inclusion of the hagfish genes. The results of these analyses suggest that cyclostome EmxB is orthologous to gnathostome Emx2, which was also supported by conserved synteny. Homopolymeric amino acid (HPAA) tracts represent a remarkable feature of the lamprey Emx sequences, and a comparative genome-wide scan revealed that lamprey proteins exhibit a unique pattern of HPAA tract accumulation.

CONCLUSIONS

Our analysis, including hagfish Emx genes, suggests that gene duplications gave rise to Emx1, -2 and -3 before the cyclostome-gnathostome split. We propose that independent HPAA tract accumulations in multiple ancient duplicates, as identified in lamprey Emx gene products, may have led to erroneous identification of gene duplication in the lamprey lineage. Overall, our reanalysis favors the scenario that the nested Emx expression pattern in mouse and lamprey shares a common origin.

Enigmatic orthology relationships between Hox clusters of the African butterfly fish and other teleosts following ancient whole-genome duplication

Numerous ancient whole-genome duplications (WGD) have occurred during eukaryote evolution. In vertebrates, duplicated developmental genes and their functional divergence have had important consequences for morphological evolution. Although two vertebrate WGD events (1R/2R) occurred over 525 Ma, we have focused on the more recent 3R or TGD (teleost genome duplication) event which occurred approximately 350 Ma in a common ancestor of over 26,000 species of teleost fishes. Through a combination of whole genome and bacterial artificial chromosome clone sequencing we characterized all Hox gene clusters of Pantodon buchholzi, a member of the early branching teleost subdivision Osteoglossomorpha. We find 45 Hox genes organized in only five clusters indicating that Pantodon has suffered more Hox cluster loss than other known species. Despite strong evidence for homology of the five Pantodon clusters to the four canonical pre-TGD vertebrate clusters (one HoxA, two HoxB, one HoxC, and one HoxD), we were unable to confidently resolve 1:1 orthology relationships between four of the Pantodon clusters and the eight post-TGD clusters of other teleosts. Phylogenetic analysis revealed that many Pantodon genes segregate outside the conventional "a" and "b" post-TGD orthology groups, that extensive topological incongruence exists between genes physically linked on a single cluster, and that signal divergence causes ambivalence in assigning 1:1 orthology in concatenated Hox cluster analyses. Out of several possible explanations for this phenomenon we favor a model which keeps with the prevailing view of a single TGD prior to teleost radiation, but which also considers the timing of diploidization after duplication, relative to speciation events. We suggest that although the duplicated hoxa clusters diploidized prior to divergence of osteoglossomorphs, the duplicated hoxb, hoxc, and hoxd clusters concluded diploidization independently in osteoglossomorphs and other teleosts. We use the term "tetralogy" to describe the homology relationship which exists between duplicated sequences which originate through a shared WGD, but which diploidize into distinct paralogs from a common allelic pool independently in two lineages following speciation.

The globin gene repertoire of lampreys: convergent evolution of hemoglobin and myoglobin in jawed and jawless vertebrates

Agnathans (jawless vertebrates) occupy a key phylogenetic position for illuminating the evolution of vertebrate anatomy and physiology. Evaluation of the agnathan globin gene repertoire can thus aid efforts to reconstruct the origin and evolution of the globin genes of vertebrates, a superfamily that includes the well-known model proteins hemoglobin and myoglobin. Here, we report a comprehensive analysis of the genome of the sea lamprey (Petromyzon marinus) which revealed 23 intact globin genes and two hemoglobin pseudogenes. Analyses of the genome of the Arctic lamprey (Lethenteron camtschaticum) identified 18 full length and five partial globin gene sequences. The majority of the globin genes in both lamprey species correspond to the known agnathan hemoglobins. Both genomes harbor two copies of globin X, an ancient globin gene that has a broad phylogenetic distribution in the animal kingdom. Surprisingly, we found no evidence for an ortholog of neuroglobin in the lamprey genomes. Expression and phylogenetic analyses identified an ortholog of cytoglobin in the lampreys; in fact, our results indicate that cytoglobin is the only orthologous vertebrate-specific globin that has been retained in both gnathostomes and agnathans. Notably, we also found two globins that are highly expressed in the heart of P. marinus, thus representing functional myoglobins. Both genes have orthologs in L. camtschaticum. Phylogenetic analyses indicate that these heart-expressed globins are not orthologous to the myoglobins of jawed vertebrates (Gnathostomata), but originated independently within the agnathans. The agnathan myoglobin and hemoglobin proteins form a monophyletic group to the exclusion of functionally analogous myoglobins and hemoglobins of gnathostomes, indicating that specialized respiratory proteins for O2 transport in the blood and O2 storage in the striated muscles evolved independently in both lineages. This dual convergence of O2-transport and O2-storage proteins in agnathans and gnathostomes involved the convergent co-option of different precursor proteins in the ancestral globin repertoire of vertebrates.

Revealing less derived nature of cartilaginous fish genomes with their evolutionary time scale inferred with nuclear genes

Cartilaginous fishes, divided into Holocephali (chimaeras) and Elasmoblanchii (sharks, rays and skates), occupy a key phylogenetic position among extant vertebrates in reconstructing their evolutionary processes. Their accurate evolutionary time scale is indispensable for better understanding of the relationship between phenotypic and molecular evolution of cartilaginous fishes. However, our current knowledge on the time scale of cartilaginous fish evolution largely relies on estimates using mitochondrial DNA sequences. In this study, making the best use of the still partial, but large-scale sequencing data of cartilaginous fish species, we estimate the divergence times between the major cartilaginous fish lineages employing nuclear genes. By rigorous orthology assessment based on available genomic and transcriptomic sequence resources for cartilaginous fishes, we selected 20 protein-coding genes in the nuclear genome, spanning 2973 amino acid residues. Our analysis based on the Bayesian inference resulted in the mean divergence time of 421 Ma, the late Silurian, for the Holocephali-Elasmobranchii split, and 306 Ma, the late Carboniferous, for the split between sharks and rays/skates. By applying these results and other documented divergence times, we measured the relative evolutionary rate of the Hox A cluster sequences in the cartilaginous fish lineages, which resulted in a lower substitution rate with a factor of at least 2.4 in comparison to tetrapod lineages. The obtained time scale enables mapping phenotypic and molecular changes in a quantitative framework. It is of great interest to corroborate the less derived nature of cartilaginous fish at the molecular level as a genome-wide phenomenon.

aLeaves facilitates on-demand exploration of metazoan gene family trees on MAFFT sequence alignment server with enhanced interactivity

We report a new web server, aLeaves (http://aleaves.cdb.riken.jp/), for homologue collection from diverse animal genomes. In molecular comparative studies involving multiple species, orthology identification is the basis on which most subsequent biological analyses rely. It can be achieved most accurately by explicit phylogenetic inference. More and more species are subjected to large-scale sequencing, but the resultant resources are scattered in independent project-based, and multi-species, but separate, web sites. This complicates data access and is becoming a serious barrier to the comprehensiveness of molecular phylogenetic analysis. aLeaves, launched to overcome this difficulty, collects sequences similar to an input query sequence from various data sources. The collected sequences can be passed on to the MAFFT sequence alignment server (http://mafft.cbrc.jp/alignment/server/), which has been significantly improved in interactivity. This update enables to switch between (i) sequence selection using the Archaeopteryx tree viewer, (ii) multiple sequence alignment and (iii) tree inference. This can be performed as a loop until one reaches a sensible data set, which minimizes redundancy for better visibility and handling in phylogenetic inference while covering relevant taxa. The work flow achieved by the seamless link between aLeaves and MAFFT provides a convenient online platform to address various questions in zoology and evolutionary biology.

Non-parsimonious evolution of hagfish Dlx genes

Background

The number of members of the Dlx gene family increased during the two rounds of whole-genome duplication that occurred in the common ancestor of the vertebrates. Because the Dlx genes are involved in the development of the cranial skeleton, brain, and sensory organs, their expression patterns have been analysed in various organisms in the context of evolutionary developmental biology. Six Dlx genes have been isolated in the lampreys, a group of living jawless vertebrates (cyclostomes), and their expression patterns analysed. However, little is known about the Dlx genes in the hagfish, the other cyclostome group, mainly because the embryological analysis of this animal is difficult.

Results

To identify the hagfish Dlx genes and describe their expression patterns, we cloned the cDNA from embryos of the Japanese inshore hagfish Eptatretus burgeri. Our results show that the hagfish has at least six Dlx genes and one pseudogene. In a phylogenetic analysis, the hagfish Dlx genes and those of the lampreys tended to be excluded from the clade of the gnathostome Dlx genes. In several cases, the lamprey Dlx genes clustered with the clade consisting of two hagfish genes, suggesting that independent gene duplications have occurred in the hagfish lineage. Analysis of the expression of these genes showed distinctive overlapping expression patterns in the cranial mesenchymal cells and the inner ear.

Conclusions

Independent duplication, pseudogenization, and loss of the Dlx genes probably occurred in the hagfish lineage after its split from the other vertebrate lineages. This pattern is reminiscent of the non-parsimonious evolution of its morphological traits, including its inner ear and vertebrae, which indicate that this group is an early-branching lineage that diverged before those characters evolved.

The Dlx genes as clues to vertebrate genomics and craniofacial evolution

The group of Dlx genes belongs to the homeobox-containing superfamily, and its members are involved in various morphogenetic processes. In vertebrate genomes, Dlx genes exist as multiple paralogues generated by tandem duplication followed by whole genome duplications. In this review, we provide an overview of the Dlx gene phylogeny with an emphasis on the chordate lineage. Referring to the Dlx gene repertoire, we discuss the establishment and conservation of the nested expression patterns of the Dlx genes in craniofacial development. Despite the accumulating genomic sequence resources in diverse vertebrates, embryological analyses of Dlx gene expression and function remain limited in terms of species diversity. By supplementing our original analysis of shark embryos with previous data from other osteichthyans, such as mice and zebrafish, we support the previous speculation that the nested Dlx expression in the pharyngeal arch is likely a shared feature among all the extant jawed vertebrates. Here, we highlight several hitherto unaddressed issues regarding the evolution and function of Dlx genes, with special reference to the craniofacial development of vertebrates.

Hox gene clusters of early vertebrates: do they serve as reliable markers for genome evolution?

Hox genes, responsible for regional specification along the anteroposterior axis in embryogenesis, are found as clusters in most eumetazoan genomes sequenced to date. Invertebrates possess a single Hox gene cluster with some exceptions of secondary cluster breakages, while osteichthyans (bony vertebrates) have multiple Hox clusters. In tetrapods, four Hox clusters, derived from the so-called two-round whole genome duplications (2R-WGDs), are observed. Overall, the number of Hox gene clusters has been regarded as a reliable marker of ploidy levels in animal genomes. In fact, this scheme also fits the situations in teleost fishes that experienced an additional WGD. In this review, I focus on cyclostomes and cartilaginous fishes as lineages that would fill the gap between invertebrates and osteichthyans. A recent study highlighted a possible loss of the HoxC cluster in the galeomorph shark lineage, while other aspects of cartilaginous fish Hox clusters usually mark their conserved nature. In contrast, existing resources suggest that the cyclostomes exhibit a different mode of Hox cluster organization. For this group of species, whose genomes could have differently responded to the 2R-WGDs from jawed vertebrates, therefore the number of Hox clusters may not serve as a good indicator of their ploidy level.

Detection and phylogenetic assessment of conserved synteny derived from whole genome duplications

Identification of intragenomic conservation of gene compositions in multiple chromosomal segments led to evidence of whole genome (WGDs) duplications. The process by which WGDs have been maintained and decayed provides us with clues for understanding how the genome evolves. In this chapter, we summarize current understanding of phylogenetic distribution and evolutionary impact of WGDs, introduce basic procedures to detect conserved synteny, and discuss typical pitfalls, as well as biological insights.

Genome-wide detection of gene extinction in early mammalian evolution

Detecting gene losses is a novel aspect of evolutionary genomics that has been made feasible by whole-genome sequencing. However, research to date has concentrated on elucidating evolutionary patterns of genomic components shared between species, rather than identifying disparities between genomes. In this study, we searched for gene losses in the lineage leading to eutherian mammals. First, as a pilot analysis, we selected five gene families (Wnt, Fgf, Tbx, TGFβ, and Frizzled) for molecular phylogenetic analyses, and identified mammalian lineage-specific losses of Wnt11b, Tbx6L/VegT/tbx16, Nodal-related, ADMP1, ADMP2, Sizzled, and Crescent. Second, automated genome-wide phylogenetic screening was implemented based on this pilot analysis. As a result, we detected 147 chicken genes without eutherian orthologs, which resulted from 141 gene loss events. Our inventory contained a group of regulatory genes governing early embryonic axis formation, such as Noggins, and multiple members of the opsin and prolactin-releasing hormone receptor (“PRLHR”) gene families. Our findings highlight the potential of genome-wide gene phylogeny (“phylome”) analysis in detecting possible rearrangement of gene networks and the importance of identifying losses of ancestral genomic components in analyzing the molecular basis underlying phenotypic evolution.

Analysis of lamprey clustered Fox genes: insight into Fox gene evolution and expression in vertebrates

In the human genome, members of the FoxC, FoxF, FoxL1, and FoxQ1 gene families are found in two paralagous clusters. One cluster contains the genes FOXQ1, FOXF2, FOXC1 and the second consists of FOXF1, FOXC2, and FOXL1. In jawed vertebrates these genes are known to be expressed in different pharyngeal tissues and all, except FoxQ1, are involved in patterning the early embryonic mesoderm. We have previously traced the evolution of this cluster in the bony vertebrates, and the gene content is identical in the dogfish, a member of the most basally branching lineage of the jawed vertebrates. Here we extend these analyses to jawless vertebrates. Using genomic searches and molecular approaches we have identified homologues of these genes from lampreys. We identify two FoxC genes, two FoxF genes, two FoxQ1 genes and single FoxL1 gene. We examine the embryonic expression of one predominantly mesodermally expressed gene family, FoxC, and the endodermally expressed member of the cluster, FoxQ1. We identified FoxQ1 transcripts in the pharyngeal endoderm, while the two FoxC genes are differentially expressed in the pharyngeal mesenchyme and ectoderm. Furthermore we identify conserved expression of lamprey FoxC genes in the paraxial and intermediate mesoderms. We interpret our results through a chordate-wide comparison of expression patterns and discuss gene content in the context of theories on the evolution of the vertebrate genome.

Revisiting the origin of the vertebrate Hox14 by including its relict sarcopterygian members

Bilaterian Hox genes play pivotal roles in the specification of positional identities along the anteroposterior axis. Particularly in vertebrates, their regulation is tightly coordinated by tandem arrays of genes [paralogy groups (PGs)] in four gene clusters (HoxA-D). Traditionally, the uninterrupted Hox cluster (Hox1-14) of the invertebrate chordate amphioxus was regarded as an archetype of the vertebrate Hox clusters. In contrast to Hox1-13 that are globally regulated by the "Hox code" and are often phylogenetically conserved, vertebrate Hox14 members were only recently revealed to be present in an African lungfish, a coelacanth, chondrichthyans and a lamprey, and decoupled from the Hox code. In this study we performed a PCR-based search of Hox14 members from diverse vertebrates, and identified one in the Australian lungfish, Neoceratodus forsteri. Based on a molecular phylogenetic analysis, this gene was designated NfHoxA14. Our real-time RT-PCR suggested its hindgut-associated expression, previously observed also in cloudy catshark HoxD14 and lamprey Hox14α. It is likely that this altered expression scheme was established before the Hox cluster quadruplication, probably at the base of extant vertebrates. To investigate the origin of vertebrate Hox14, by including this sarcopterygian Hox14 member, we performed focused phylogenetic analyses on its relationship with other vertebrate posterior Hox PGs (Hox9-13) as well as amphioxus posterior Hox genes. Our results confirmed the hypotheses previously proposed by other studies that vertebrate Hox14 does not have any amphioxus ortholog, and that none of 1-to-1 pairs of vertebrate and amphioxus posterior Hox genes, based on their relative location in the clusters, is orthologous.

Overview of the transcriptome profiles identified in hagfish, shark, and bichir: current issues arising from some nonmodel vertebrate taxa

Because of their crucial phylogenetic positions, hagfishes, sharks, and bichirs are recognized as key taxa in our understanding of vertebrate evolution. The expression patterns of the regulatory genes involved in developmental patterning have been analyzed in the context of evolutionary developmental studies. However, in a survey of public sequence databases, we found that the large-scale sequence data for these taxa are still limited. To address this deficit, we used conventional Sanger DNA sequencing and a next-generation sequencing technology based on 454 GS FLX sequencing to obtain expressed sequence tags (ESTs) of the Japanese inshore hagfish (Eptatretus burgeri; 161,482 ESTs), cloudy catshark (Scyliorhinus torazame; 165,819 ESTs), and gray bichir (Polypterus senegalus; 34,336 ESTs). We deposited the ESTs in a newly constructed database, designated the "Vertebrate TimeCapsule." The ESTs include sequences from genes that can be effectively used in evolutionary developmental studies; for instance, several encode cartilaginous extracellular matrix proteins, which are central to an understanding of the ways in which evolutionary processes affected the skeletal elements, whereas others encode regulatory genes involved in craniofacial development and early embryogenesis. Here, we discuss how hagfishes, sharks, and bichirs contribute to our understanding of vertebrate evolution, we review the current status of the publicly available sequence data for these three taxa, and we introduce our EST projects and newly developed database.

Two lamprey Hedgehog genes share non-coding regulatory sequences and expression patterns with gnathostome Hedgehogs

Hedgehog (Hh) genes play major roles in animal development and studies of their evolution, expression and function point to major differences among chordates. Here we focused on Hh genes in lampreys in order to characterize the evolution of Hh signalling at the emergence of vertebrates. Screening of a cosmid library of the river lamprey Lampetra fluviatilis and searching the preliminary genome assembly of the sea lamprey Petromyzon marinus indicate that lampreys have two Hh genes, named Hha and Hhb. Phylogenetic analyses suggest that Hha and Hhb are lamprey-specific paralogs closely related to Sonic/Indian Hh genes. Expression analysis indicates that Hha and Hhb are expressed in a Sonic Hh-like pattern. The two transcripts are expressed in largely overlapping but not identical domains in the lamprey embryonic brain, including a newly-described expression domain in the nasohypophyseal placode. Global alignments of genomic sequences and local alignment with known gnathostome regulatory motifs show that lamprey Hhs share conserved non-coding elements (CNE) with gnathostome Hhs albeit with sequences that have significantly diverged and dispersed. Functional assays using zebrafish embryos demonstrate gnathostome-like midline enhancer activity for CNEs contained in intron2. We conclude that lamprey Hh genes are gnathostome Shh-like in terms of expression and regulation. In addition, they show some lamprey-specific features, including duplication and structural (but not functional) changes in the intronic/regulatory sequences.