Unresolved orthology and peculiar coding sequence properties of lamprey genes: the KCNA gene family as test case

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

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

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

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

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.

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

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

Evolutionary history of the vertebrate Piwi gene family

PIWIs are regulatory proteins that belong to the Argonaute family. Piwis are primarily expressed in gonads and protect the germline against the mobilization and propagation of transposable elements (TEs) through transcriptional gene silencing. Vertebrate genomes encode up to four Piwi genes: Piwil1, Piwil2, Piwil3 and Piwil4, but their duplication history is unresolved. We leveraged phylogenetics, synteny and expression analyses to address this void. Our phylogenetic analysis suggests Piwil1 and Piwil2 were retained in all vertebrate members. Piwil4 was the result of Piwil1 duplication in the ancestor of gnathostomes, but was independently lost in ray-finned fishes and birds. Further, Piwil3 was derived from a tandem Piwil1 duplication in the common ancestor of marsupial and placental mammals, but was secondarily lost in Atlantogenata (Xenarthra and Afrotheria) and some rodents. The evolutionary rate of Piwil3 is considerably faster than any Piwi among all lineages, but an explanation is lacking. Our expression analyses suggest Piwi expression has mostly been constrained to gonads throughout vertebrate evolution. Vertebrate evolution is marked by two early rounds of whole genome duplication and many multigene families are linked to these events. However, our analyses suggest Piwi expansion was independent of whole genome duplications.

Analysis of lamprey meis genes reveals that conserved inputs from Hox, Meis and Pbx proteins control their expression in the hindbrain and neural tube

Meis genes are known to play important roles in the hindbrain and neural crest cells of jawed vertebrates. To explore the roles of Meis genes in head development during evolution of vertebrates, we have identified four meis genes in the sea lamprey genome and characterized their patterns of expression and regulation, with a focus on the hindbrain and pharynx. Each of the lamprey meis genes displays temporally and spatially dynamic patterns of expression, some of which are coupled to rhombomeric domains in the developing hindbrain and select pharyngeal arches. Studies of Meis loci in mouse and zebrafish have identified enhancers that are bound by Hox and TALE (Meis and Pbx) proteins, implicating these factors in the direct regulation of Meis expression. We examined the lamprey meis loci and identified a series of cis-elements conserved between lamprey and jawed vertebrate meis genes. In transgenic reporter assays we demonstrated that these elements act as neural enhancers in lamprey embryos, directing reporter expression in appropriate domains when compared to expression of their associated endogenous meis gene. Sequence alignments reveal that these conserved elements are in similar relative positions of the meis loci and contain a series of consensus binding motifs for Hox and TALE proteins. This suggests that ancient Hox and TALE-responsive enhancers regulated expression of ancestral vertebrate meis genes in segmental domains in the hindbrain and have been retained in the meis loci during vertebrate evolution. The presence of conserved Meis, Pbx and Hox binding sites in these lamprey enhancers links Hox and TALE factors to regulation of lamprey meis genes in the developing hindbrain, indicating a deep ancestry for these regulatory interactions prior to the divergence of jawed and jawless vertebrates.

Chromosome karyotype and stability of new synthetic hexaploid wheat

Synthetic hexaploid wheat offers breeders ready access to potentially novel genetic variation in wild ancestral species. In this study, we crossed MY3478 (2n = 4x = 28, AABB) as the maternal parent with the stripe rust-resistant SY41 (2n = 2x = 14, DD) as the paternal parent to construct the new hexaploid wheat line NA0928 through natural allopolyploidization. Agronomic traits and the cytology of the S8-S9 generations of NA0928 were analyzed. Abundant variation in agronomic traits was observed among each strain of NA0928 in the S8 generation. Agronomic traits were superior in strains resistant to stripe rust compared with those of highly susceptible strains. The rank order of the coefficients of variation were tiller number (55.3%) > spike length (15.3%) > number of spikelets (13.9%) > plant height (8.7). Number of tillers and spike length are important traits in wheat breeding to improve yield. Cytological observation and fluorescence in situ hybridization showed that the chromosome number and configuration showed rich variation among NA0928 strains in the S9 generation. Chromosome number ranged from 36 to 44. Variation in chromosome karyotype was detected in the A and B subgenomes. Meiotic chromosome behavior in pollen mother cells and multicolor genomic in situ hybridization revealed that two new synthetic hexaploid wheat strains showed genetic stability; one strain was resistant to stripe rust and developed multiple tillers, and the other strain was susceptible to stripe rust, but both showed improved thousand-kernel weight (TKW) weight and produced multiple tillers. The two strains will be valuable germplasm resources for use in wheat breeding.

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.

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.

Shark and ray genomics for disentangling their morphological diversity and vertebrate evolution

Developmental studies of sharks and rays (elasmobranchs) have provided much insight into the process of morphological evolution of vertebrates. Although those studies are supposedly fueled by large-scale molecular sequencing information, whole-genome sequences of sharks and rays were made available only recently. One compelling difficulty of elasmobranch developmental biology is the low accessibility to embryonic study materials and their slow development. Another limiting factor is the relatively large size of their genomes. Moreover, their large body sizes restrict sustainable captive breeding, while their high body fluid osmolarity prevents reproducible cell culturing for in vitro experimentation, which has also limited our knowledge of their chromosomal organization for validation of genome sequencing products. This article focuses on egg-laying elasmobranch species used in developmental biology and provides an overview of the characteristics of the shark and ray genomes revealed to date. Developmental studies performed on a gene-by-gene basis are also reviewed from a whole-genome perspective. Among the popular regulatory genes studied in developmental biology, I scrutinize shark homologs of Wnt genes that highlight vanishing repertoires in many other vertebrate lineages, as well as Hox genes that underwent an unexpected modification unique to the elasmobranch lineage. These topics are discussed together with insights into the reconstruction of developmental programs in the common ancestor of vertebrates and its subsequent evolutionary trajectories that mark the features that are unique to, and those characterizing the diversity among, cartilaginous fishes.

Evolutionary history of histamine receptors: Early vertebrate origin and expansion of the H3-H4 subtypes

Histamine receptors belonging to the superfamily of G protein-coupled receptors (GPCRs) mediate the diverse biological effects of biogenic histamine. They are classified into four phylogenetically distinct subtypes H₁-H₄, each with a different binding affinity for histamine and divergent downstream signaling pathways. Here we present the evolutionary history of the histamine receptors using a phylogenetic approach complemented with comparative genomics analyses of the sequences, gene structures, and synteny of gene neighborhoods. The data indicate the earliest emergence of histamine-mediated GPCR signaling by a H₂ in a prebilaterian ancestor. The analyses support a revised classification of the vertebrate H₃-H₄ receptor subtypes. We demonstrate the presence of the H₄ across vertebrates, contradicting the currently held notion that H₄ is restricted to mammals. These non-mammalian vertebrate H₄ orthologs have been mistaken for H₃. We also identify the presence of a new H₃ subtype (H₃B), distinct from the canonical H₃ (H₃A), and propose that the H₃A, H₃B, and H₄ likely emerged from a H₃ progenitor through the 1R/2R whole genome duplications in an ancestor of the vertebrates. It is apparent that the ability of the H₁, H₂, and H_3-4 to bind histamine was acquired convergently. We identified genomic signatures suggesting that the H₁ and H₃-H₄ shared a last common ancestor with the muscarinic receptor in a bilaterian predecessor whereas, the H₂ and the α-adrenoreceptor shared a progenitor in a prebilaterian ancestor. Furthermore, site-specific analysis of the vertebrate subtypes revealed potential residues that may account for the functional divergence between them.

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

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

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

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

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

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

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

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

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.

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

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

The p53 gene family in vertebrates: Evolutionary considerations

The origin of the p53 gene family predates multicellular life since TP53 members of this gene family have been found in unicellular eukaryotes. In invertebrates one or two genes attributable to a TP53-like or TP63/73-like gene are present. The radiation into three genes, TP53, TP63, and TP73, has been reported as a vertebrate invention. TP53 is considered the "guardian of the genome" given its role in protecting cells against the DNA damage and cellular stressors. TP63 and TP73 play a role in epithelial development and neurogenesis, respectively. The evolution of the p53 gene family has been the subject of considerable analyses even if several questions remain still open. In this study we addressed the evolutionary history of the p53 gene family in vertebrates performing an extended microsyntenic investigation coupled with a phylogenetic analysis, together with protein domain organization and structure assessment. On the basis of our results we discussed a possible evolutionary scenario according to which a TP53/63/73 ancestor form gave rise to the current TP53 and a TP63/73 form, which in turn independently duplicated into two genes in agnathe and gnathostome lineages.

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.

The Reprimo Gene Family: A Novel Gene Lineage in Gastric Cancer with Tumor Suppressive Properties

The reprimo (RPRM) gene family is a group of single exon genes present exclusively within the vertebrate lineage. Two out of three members of this family are present in humans: RPRM and RPRM-Like (RPRML). RPRM induces cell cycle arrest at G2/M in response to p53 expression. Loss-of-expression of RPRM is related to increased cell proliferation and growth in gastric cancer. This evidence suggests that RPRM has tumor suppressive properties. However, the molecular mechanisms and signaling partners by which RPRM exerts its functions remain unknown. Moreover, scarce studies have attempted to characterize RPRML, and its functionality is unclear. Herein, we highlight the role of the RPRM gene family in gastric carcinogenesis, as well as its potential applications in clinical settings. In addition, we summarize the current knowledge on the phylogeny and expression patterns of this family of genes in embryonic zebrafish and adult humans. Strikingly, in both species, RPRM is expressed primarily in the digestive tract, blood vessels and central nervous system, supporting the use of zebrafish for further functional characterization of RPRM. Finally, drawing on embryonic and adult expression patterns, we address the potential relevance of RPRM and RPRML in cancer. Active investigation or analytical research in the coming years should contribute to novel translational applications of this poorly understood gene family as potential biomarkers and development of novel cancer therapies.

Molecular and evolutionary perspectives of the renin-angiotensin system from lamprey

The recent advance and revision on the renin-angiotensin system in lamprey were summarized and we emphasized that presence of two types of angiotensins (Angs) in lamprey. Due to the parasitic nature on fish blood, teleost-type Angs were produced in their buccal gland and secreted into the lamphredin to evade the host immunorejection. A native lamprey angiotensinogen (AGT) was identified in genome and it retains serine-protease inhibitor activity for thrombin that regulates the blood coagulation pathway. The native lamprey angiotensin II (Lp-Ang II) is hypotensive instead of hypertensive, suggesting a functional divergence on cardiovascular regulation from the main vertebrate groups. The renin gene was absent from the lamprey genome so far, and the mutation on the renin-recognition site on lamprey AGT suggested that other proteases may have replaced the role of renin. Lp-Ang II was shown to bind to AT1 receptor and internalized, but the downstream signaling was still unknown. Molecular and phylogenetic evidence on invertebrate ACE-like proteins indicated that they were not homologous to those in vertebrates and could be acting on other native peptides. Although it was generally believed that the RAS was a well-conserved hormone system in vertebrates and invertebrates, revision by molecular data indicated that invertebrates lack homologous RAS components while lamprey possess an almost complete RAS. This suggests that the hormone cascade system was first evolved around cyclostome emergence and invertebrates could have taken up the RAS components from vertebrates through horizontal gene transfer.

The sea lamprey germline genome provides insights into programmed genome rearrangement and vertebrate evolution

The sea lamprey (Petromyzon marinus) serves as a comparative model for reconstructing vertebrate evolution. To enable more informed analyses, we developed a new assembly of the lamprey germline genome that integrates several complementary datasets. Analysis of this highly contiguous (chromosome-scale) assembly reveals that both chromosomal and whole-genome duplications have played significant roles in the evolution of ancestral vertebrate and lamprey genomes, including chromosomes that carry the six lamprey HOX clusters. The assembly also contains several hundred genes that are reproducibly eliminated from somatic cells during early development in lamprey. Comparative analyses show that gnathostome (mouse) homologs of these genes are frequently marked by Polycomb Repressive Complexes (PRCs) in embryonic stem cells, suggesting overlaps in the regulatory logic of somatic DNA elimination and repressive/bivalent states that are regulated by early embryonic PRCs. This new assembly will enhance diverse studies that are informed by lampreys’ unique biology and evolutionary/comparative perspective.

Voltage-gated sodium channel gene repertoire of lampreys: gene duplications, tissue-specific expression and discovery of a long-lost gene

Studies of the voltage-gated sodium (Nav) channels of extant gnathostomes have made it possible to deduce that ancestral gnathostomes possessed four voltage-gated sodium channel genes derived from a single ancestral chordate gene following two rounds of genome duplication early in vertebrates. We investigated the Nav gene family in two species of lampreys (the Japanese lamprey Lethenteron japonicum and sea lamprey Petromyzon marinus) (jawless vertebrates-agnatha) and compared them with those of basal vertebrates to better understand the origin of Nav genes in vertebrates. We noted six Nav genes in both lamprey species, but orthology with gnathostome (jawed vertebrate) channels was inconclusive. Surprisingly, the Nav2 gene, ubiquitously found in invertebrates and believed to have been lost in vertebrates, is present in lampreys, elephant shark (Callorhinchus milii) and coelacanth (Latimeria chalumnae). Despite repeated duplication of the Nav1 family in vertebrates, Nav2 is only in single copy in those vertebrates in which it is retained, and was independently lost in ray-finned fishes and tetrapods. Of the other five Nav channel genes, most were expressed in brain, one in brain and heart, and one exclusively in skeletal muscle. Invertebrates do not express Nav channel genes in muscle. Thus, early in the vertebrate lineage Nav channels began to diversify and different genes began to express in heart and muscle.

Evolution of the α2-adrenoreceptors in vertebrates: ADRA2D is absent in mammals and crocodiles

Evolutionary studies of genes that have been functionally characterized and whose variation has been associated with pathological conditions represent an opportunity to understand the genetic basis of pathologies. α₂-Adrenoreceptors (ADRA2) are a class of G protein-coupled receptors that regulate several physiological processes including blood pressure, platelet aggregation, insulin secretion, lipolysis, and neurotransmitter release. This gene family has been extensively studied from a molecular/physiological perspective, yet much less is known about its evolutionary history. Accordingly, the goal of this study was to investigate the evolutionary history of α₂-adrenoreceptors (ADRA2) in vertebrates. Our results show that in addition to the three well-recognized α₂-adrenoreceptor genes (ADRA2A, ADRA2B and ADRA2C), we recovered a clade that corresponds to the fourth member of the α₂-adrenoreceptor gene family (ADRA2D). We also recovered a clade that possesses two ADRA2 sequences found in two lamprey species. Furthermore, our results show that mammals and crocodiles are characterized by possessing three α₂-adrenoreceptor genes, whereas all other vertebrate groups possess the full repertoire of α₂-adrenoreceptor genes. Among vertebrates ADRA2D seems to be a dispensable gene, as it was lost two independent times during the evolutionary history of the group. Additionally, we found that most examined species possess the most common alleles described for humans; however, there are cases in which non-human mammals possess the alternative variant. Finally, transcript abundance profiles revealed that during the early evolutionary history of gnathostomes, the expression of ADRA2D in different taxonomic groups became specialized to different tissues, but in the ancestor of sarcopterygians this specialization would have been lost.

Lampreys, the jawless vertebrates, contain only two ParaHox gene clusters

ParaHox genes (Gsx, Pdx, and Cdx) are an ancient family of developmental genes closely related to the Hox genes. They play critical roles in the patterning of brain and gut. The basal chordate, amphioxus, contains a single ParaHox cluster comprising one member of each family, whereas nonteleost jawed vertebrates contain four ParaHox genomic loci with six or seven ParaHox genes. Teleosts, which have experienced an additional whole-genome duplication, contain six ParaHox genomic loci with six ParaHox genes. Jawless vertebrates, represented by lampreys and hagfish, are the most ancient group of vertebrates and are crucial for understanding the origin and evolution of vertebrate gene families. We have previously shown that lampreys contain six Hox gene loci. Here we report that lampreys contain only two ParaHox gene clusters (designated as α- and β-clusters) bearing five ParaHox genes (Gsxα, Pdxα, Cdxα, Gsxβ, and Cdxβ). The order and orientation of the three genes in the α-cluster are identical to that of the single cluster in amphioxus. However, the orientation of Gsxβ in the β-cluster is inverted. Interestingly, Gsxβ is expressed in the eye, unlike its homologs in jawed vertebrates, which are expressed mainly in the brain. The lamprey Pdxα is expressed in the pancreas similar to jawed vertebrate Pdx genes, indicating that the pancreatic expression of Pdx was acquired before the divergence of jawless and jawed vertebrate lineages. It is likely that the lamprey Pdxα plays a crucial role in pancreas specification and insulin production similar to the Pdx of jawed vertebrates.

Evolution of the β-adrenoreceptors in vertebrates

The study of the evolutionary history of genes related to human disease lies at the interface of evolution and medicine. These studies provide the evolutionary context on which medical researchers should work, and are also useful in providing information to suggest further genetic experiments, especially in model species where genetic manipulations can be made. Here we studied the evolution of the β-adrenoreceptor gene family in vertebrates with the aim of adding an evolutionary framework to the already abundant physiological information. Our results show that in addition to the three already described vertebrate β-adrenoreceptor genes there is an additional group containing cyclostome sequences. We suggest that β-adrenoreceptors diversified as a product of the two whole genome duplications that occurred in the ancestor of vertebrates. Gene expression patterns are in general consistent across species, suggesting that expression dynamics were established early in the evolutionary history of vertebrates, and have been maintained since then. Finally, amino acid polymorphisms that are associated to pathological conditions in humans appear to be common in non-human mammals, suggesting that the phenotypic effects of these mutations depend on epistatic interaction with other positions. The evolutionary analysis of the β-adrenoreceptors delivers new insights about the diversity of these receptors in vertebrates, the evolution of the expression patterns and a comparative perspective regarding the polymorphisms that in humans are linked to pathological conditions.

Evolution of the RH gene family in vertebrates revealed by brown hagfish (Eptatretus atami) genome sequences

In vertebrates, there are four major genes in the RH (Rhesus) gene family, RH, RHAG, RHBG, and RHCG. These genes are thought to have been formed by the two rounds of whole-genome duplication (2R-WGD) in the common ancestor of all vertebrates. In our previous work, where we analyzed details of the gene duplications process of this gene family, three nucleotide sequences belonging to this family were identified in Far Eastern brook lamprey (Lethenteron reissneri), and the phylogenetic positions of the genes were determined. Lampreys, along with hagfishes, are cyclostomata (jawless fishes), which is a sister group of gnathostomata (jawed vertebrates). Although those results suggested that one gene was orthologous to the gnathostome RHCG genes, we did not identify clear orthologues for other genes. In this study, therefore, we identified three novel cDNA sequences that belong to the RH gene family using de novo transcriptome analysis of another cyclostome: the brown hagfish (Eptatretus atami). We also determined the nucleotide sequences for the RHBG and RHCG genes in a red stingray (Dasyatis akajei), which belongs to the cartilaginous fishes. The phylogenetic tree showed that two brown hagfish genes, which were probably duplicated in the cyclostome lineage, formed a cluster with the gnathostome RHAG genes, whereas another brown hagfish gene formed a cluster with the gnathostome RHCG genes. We estimated that the RH genes had a higher evolutionary rate than the RHAG, RHBG, and RHCG genes. Interestingly, in the RHBG genes, only the bird lineage showed a higher rate of nonsynonymous substitutions. It is likely that this higher rate was caused by a state of relaxed functional constraints rather than positive selection nor by pseudogenization.

Identification of three somatostatin genes in lampreys

Somatostatins (SSs) are a structurally diverse family of neuropeptides that play important roles in the regulation of growth, development and metabolism in vertebrates. It has been recently proposed that the common ancestor of gnathostomes possessed three SS genes, namely SS1, SS2 and SS5. SS1 and SS2 are still present in most extant gnathostome species investigated so far while SS5 primarily occurs in chondrichthyes, actinopterygians and actinistia but not in tetrapods. Very little is known about the repertoire of SSs in cyclostomes, which are extant jawless vertebrates. In the present study, we report the cloning of the cDNAs encoding three distinct lamprey SS variants that we call SSa, SSb and SSc. SSa and SSb correspond to the two SS variants previously characterized in lamprey, while SSc appears to be a totally novel one. SSa exhibits the same sequence as gnathostome SS1. SSb differs from SSa by only one substitution (Thr¹²→Ser). SSc exhibits a totally unique structure (ANCRMFYWKTMAAC) that shares only 50% identity with SSa and SSb. SSa, SSb and SSc precursors do not exhibit any appreciable sequence similarity outside the C-terminal region containing the SS sequence. Phylogenetic analyses failed to clearly assign orthology relationships between lamprey and gnathostome SS genes. Synteny analysis suggests that the SSc gene arose before the split of the three gnathostome genes SS1, SS2 and SS5.

Evolutionary history of the reprimo tumor suppressor gene family in vertebrates with a description of a new reprimo gene lineage

Genes related to human diseases should be natural targets for evolutionary studies, since they could provide clues regarding the genetic bases of pathologies and potential treatments. Here we studied the evolution of the reprimo gene family, a group of tumor-suppressor genes that are implicated in p53-mediated cell cycle arrest. These genes, especially the reprimo duplicate located on human chromosome 2, have been associated with epigenetic modifications correlated with transcriptional silencing and cancer progression. We demonstrate the presence of a third reprimo lineage that, together with the reprimo and reprimo-like genes, appears to have been differentially retained during the evolutionary history of vertebrates. We present evidence that these reprimo lineages originated early in vertebrate evolution and expanded as a result of the two rounds of whole genome duplications that occurred in the last common ancestor of vertebrates. The reprimo gene has been lost in birds, and the third reprimo gene lineage has been retained in only a few distantly related species, such as coelacanth and gar. Expression analyses revealed that the reprimo paralogs are mainly expressed in the nervous system. Different vertebrate lineages have retained different reprimo paralogs, and even in species that have retained multiple copies, only one of them is heavily expressed.

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.

Early Evolution of Vertebrate Mybs: An Integrative Perspective Combining Synteny, Phylogenetic, and Gene Expression Analyses

The genes in the Myb superfamily encode for three related transcription factors in most vertebrates, A-, B-, and c-Myb, with functionally distinct roles, whereas most invertebrates have a single Myb. B-Myb plays an essential role in cell division and cell cycle progression, c-Myb is involved in hematopoiesis, and A-Myb is involved in spermatogenesis and regulating expression of pachytene PIWI interacting RNAs, a class of small RNAs involved in posttranscriptional gene regulation and the maintenance of reproductive tissues. Comparisons between teleost fish and tetrapods suggest that the emergence and functional divergence of the Myb genes were linked to the two rounds of whole-genome duplication early in vertebrate evolution. We combined phylogenetic, synteny, structural, and gene expression analyses of the Myb paralogs from elephant shark and lampreys with data from 12 bony vertebrates to reconstruct the early evolution of vertebrate Mybs. Phylogenetic and synteny analyses suggest that the elephant shark and Japanese lamprey have copies of the A-, B-, and c-Myb genes, implying their origin could be traced back to the common ancestor of lampreys and gnathostomes. However, structural and gene expression analyses suggest that their functional roles diverged between gnathostomes and cyclostomes. In particular, we did not detect A-Myb expression in testis suggesting that the involvement of A-Myb in the pachytene PIWI interacting RNA pathway is probably a gnathostome-specific innovation. We speculate that the secondary loss of a central domain in lamprey A-Myb underlies the functional differences between the cyclostome and gnathostome A-Myb proteins.

Bioinformatic Characterization of Genes and Proteins Involved in Blood Clotting in Lampreys

Lampreys and hagfish are the earliest diverging of extant vertebrates and are obvious targets for investigating the origins of complex biochemical systems found in mammals. Currently, the simplest approach for such inquiries is to search for the presence of relevant genes in whole genome sequence (WGS) assemblies. Unhappily, in the past a high-quality complete genome sequence has not been available for either lampreys or hagfish, precluding the possibility of proving gene absence. Recently, improved but still incomplete genome assemblies for two species of lamprey have been posted, and, taken together with an extensive collection of short sequences in the NCBI trace archive, they have made it possible to make reliable counts for specific gene families. Particularly, a multi-source tactic has been used to study the lamprey blood clotting system with regard to the presence and absence of genes known to occur in higher vertebrates. As was suggested in earlier studies, lampreys lack genes for coagulation factors VIII and IX, both of which are critical for the "intrinsic" clotting system and responsible for hemophilia in humans. On the other hand, they have three each of genes for factors VII and X, participants in the "extrinsic" clotting system. The strategy of using raw trace sequence "reads" together with partial WGS assemblies for lampreys can be used in studies on the early evolution of other biochemical systems in vertebrates.

The sea lamprey meiotic map improves resolution of ancient vertebrate genome duplications

It is generally accepted that many genes present in vertebrate genomes owe their origin to two whole-genome duplications that occurred deep in the ancestry of the vertebrate lineage. However, details regarding the timing and outcome of these duplications are not well resolved. We present high-density meiotic and comparative genomic maps for the sea lamprey (Petromyzon marinus), a representative of an ancient lineage that diverged from all other vertebrates ∼550 million years ago. Linkage analyses yielded a total of 95 linkage groups, similar to the estimated number of germline chromosomes (1n ∼ 99), spanning a total of 5570.25 cM. Comparative mapping data yield strong support for the hypothesis that a single whole-genome duplication occurred in the basal vertebrate lineage, but do not strongly support a hypothetical second event. Rather, these comparative maps reveal several evolutionarily independent segmental duplications occurring over the last 600+ million years of chordate evolution. This refined history of vertebrate genome duplication should permit more precise investigations of vertebrate evolution.

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.

Ancient Duplications and Expression Divergence in the Globin Gene Superfamily of Vertebrates: Insights from the Elephant Shark Genome and Transcriptome

Comparative analyses of vertebrate genomes continue to uncover a surprising diversity of genes in the globin gene superfamily, some of which have very restricted phyletic distributions despite their antiquity. Genomic analysis of the globin gene repertoire of cartilaginous fish (Chondrichthyes) should be especially informative about the duplicative origins and ancestral functions of vertebrate globins, as divergence between Chondrichthyes and bony vertebrates represents the most basal split within the jawed vertebrates. Here, we report a comparative genomic analysis of the vertebrate globin gene family that includes the complete globin gene repertoire of the elephant shark (Callorhinchus milii). Using genomic sequence data from representatives of all major vertebrate classes, integrated analyses of conserved synteny and phylogenetic relationships revealed that the last common ancestor of vertebrates possessed a repertoire of at least seven globin genes: single copies of androglobin and neuroglobin, four paralogous copies of globin X, and the single-copy progenitor of the entire set of vertebrate-specific globins. Combined with expression data, the genomic inventory of elephant shark globins yielded four especially surprising findings: 1) there is no trace of the neuroglobin gene (a highly conserved gene that is present in all other jawed vertebrates that have been examined to date), 2) myoglobin is highly expressed in heart, but not in skeletal muscle (reflecting a possible ancestral condition in vertebrates with single-circuit circulatory systems), 3) elephant shark possesses two highly divergent globin X paralogs, one of which is preferentially expressed in gonads, and 4) elephant shark possesses two structurally distinct α-globin paralogs, one of which is preferentially expressed in the brain. Expression profiles of elephant shark globin genes reveal distinct specializations of function relative to orthologs in bony vertebrates and suggest hypotheses about ancestral functions of vertebrate globins.

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.

Evolution of the vertebrate Pax4/6 class of genes with focus on its novel member, the Pax10 gene

The members of the paired box (Pax) family regulate key developmental pathways in many metazoans as tissue-specific transcription factors. Vertebrate genomes typically possess nine Pax genes (Pax1-9), which are derived from four proto-Pax genes in the vertebrate ancestor that were later expanded through the so-called two-round (2R) whole-genome duplication. A recent study proposed that pax6a genes of a subset of teleost fishes (namely, acanthopterygians) are remnants of a paralog generated in the 2R genome duplication, to be renamed pax6.3, and reported one more group of vertebrate Pax genes (Pax6.2), most closely related to the Pax4/6 class. We propose to designate this new member Pax10 instead and reconstruct the evolutionary history of the Pax4/6/10 class with solid phylogenetic evidence. Our synteny analysis showed that Pax4, -6, and -10 originated in the 2R genome duplications early in vertebrate evolution. The phylogenetic analyses of relationships between teleost pax6a and other Pax4, -6, and -10 genes, however, do not support the proposed hypothesis of an ancient origin of the acanthopterygian pax6a genes in the 2R genome duplication. Instead, we confirmed the traditional scenario that the acanthopterygian pax6a is derived from the more recent teleost-specific genome duplication. Notably, Pax6 is present in all vertebrates surveyed to date, whereas Pax4 and -10 were lost multiple times in independent vertebrate lineages, likely because of their restricted expression patterns: Among Pax6-positive domains, Pax10 has retained expression in the adult retina alone, which we documented through in situ hybridization and quantitative reverse transcription polymerase chain reaction experiments on zebrafish, Xenopus, and anole lizard.

Origin and loss of nested LRRTM/α-catenin genes during vertebrate evolution

Leucine-rich repeat transmembrane neuronal proteins (LRRTMs) form in mammals a family of four postsynaptic adhesion proteins, which have been shown to bind neurexins and heparan sulphate proteoglycan (HSPG) glypican on the presynaptic side. Mutations in the genes encoding LRRTMs and neurexins are implicated in human cognitive disorders such as schizophrenia and autism. Our analysis shows that in most jawed vertebrates, lrrtm1, lrrtm2, and lrrtm3 genes are nested on opposite strands of large conserved intron of α-catenin genes ctnna2, ctnna1, and ctnna3, respectively. No lrrtm genes could be found in tunicates or lancelets, while two lrrtm genes are found in the lamprey genome, one of which is adjacent to a single ctnna homolog. Based on similar highly positive net charge of lamprey LRRTMs and the HSPG-binding LRRTM3 and LRRTM4 proteins, we speculate that the ancestral LRRTM might have bound HSPG before acquiring neurexins as binding partners. Our model suggests that lrrtm gene translocated into the large ctnna intron in early vertebrates, and that subsequent duplications resulted in three lrrtm/ctnna gene pairs present in most jawed vertebrates. However, we detected three prominent exceptions: (1) the lrrtm3/ctnna3 gene structure is absent in the ray-finned fish genomes, (2) the genomes of clawed frogs contain ctnna1 but lack the corresponding nested (lrrtm2) gene, and (3) contain lrrtm3 gene in the syntenic position but lack the corresponding host (ctnna3) gene. We identified several other protein-coding nested gene structures of which either the host or the nested gene has presumably been lost in the frog or chicken lineages. Interestingly, majority of these nested genes comprise LRR domains.

Insight from the lamprey genome: glimpsing early vertebrate development via neuroendocrine-associated genes and shared synteny of gonadotropin-releasing hormone (GnRH)

Study of the ancient lineage of jawless vertebrates is key to understanding the origins of vertebrate biology. The establishment of the neuroendocrine system with the hypothalamic-pituitary axis at its crux is of particular interest. Key neuroendocrine hormones in this system include the pivotal gonadotropin-releasing hormones (GnRHs) responsible for controlling reproduction via the pituitary. Previous data incorporating several lines of evidence showed all known vertebrate GnRHs were grouped into four paralogous lineages: GnRH1, 2, 3 and 4; with proposed evolutionary paths. Using the currently available lamprey genome assembly, we searched genes of the neuroendocrine system and summarize here the details representing the state of the current lamprey genome assembly. Additionally, we have analyzed in greater detail the evolutionary history of the GnRHs based on the information of the genomic neighborhood of the paralogs in lamprey as compared to other gnathostomes. Significantly, the current evidence suggests that two genome duplication events (both 1R and 2R) that generated the different fish and tetrapod paralogs took place before the divergence of the ancestral agnathans and gnathostome lineages. Syntenic analysis supports this evidence in that the previously-classified type IV GnRHs in lamprey (lGnRH-I and -III) share a common ancestry with GnRH2 and 3, and thus are no longer considered type IV GnRHs. Given the single amino acid difference between lGnRH-II and GnRH2 we propose that a GnRH2-like gene existed before the lamprey/gnathostome split giving rise to lGnRH-II and GnRH2. Furthermore, paralogous type 3 genes (lGnRH-I/III and GnRH3) evolved divergent structure/function in lamprey and gnathostome lineages.

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

Cyclostomes, comprising jawless vertebrates such as lampreys and hagfishes, are the sister group of living jawed vertebrates (gnathostomes) and hence an important group for understanding the origin and diversity of vertebrates. In vertebrates and other metazoans, Hox genes determine cell fate along the anteroposterior axis of embryos and are implicated in driving morphological diversity. Invertebrates contain a single Hox cluster (either intact or fragmented), whereas elephant shark, coelacanth, and tetrapods contain four Hox clusters owing to two rounds of whole-genome duplication ("1R" and "2R") during early vertebrate evolution. By contrast, most teleost fishes contain up to eight Hox clusters because of an additional "teleost-specific" genome duplication event. By sequencing bacterial artificial chromosome (BAC) clones and the whole genome, here we provide evidence for at least six Hox clusters in the Japanese lamprey (Lethenteron japonicum). This suggests that the lamprey lineage has experienced an additional genome duplication after 1R and 2R. The relative age of lamprey and human paralogs supports this hypothesis. Compared with gnathostome Hox clusters, lamprey Hox clusters are unusually large. Several conserved noncoding elements (CNEs) were predicted in the Hox clusters of lamprey, elephant shark, and human. Transgenic zebrafish assay indicated the potential of CNEs to function as enhancers. Interestingly, CNEs in individual lamprey Hox clusters are frequently conserved in multiple Hox clusters in elephant shark and human, implying a many-to-many orthology relationship between lamprey and gnathostome Hox clusters. Such a relationship suggests that the first two rounds of genome duplication may have occurred independently in the lamprey and gnathostome lineages.

Vertebrate paralogous conserved noncoding sequences may be related to gene expressions in brain

Vertebrate genomes include gene regulatory elements in protein-noncoding regions. A part of gene regulatory elements are expected to be conserved according to their functional importance, so that evolutionarily conserved noncoding sequences (CNSs) might be good candidates for those elements. In addition, paralogous CNSs, which are highly conserved among both orthologous loci and paralogous loci, have the possibility of controlling overlapping expression patterns of their adjacent paralogous protein-coding genes. The two-round whole-genome duplications (2R WGDs), which most probably occurred in the vertebrate common ancestors, generated large numbers of paralogous protein-coding genes and their regulatory elements. These events could contribute to the emergence of vertebrate features. However, the evolutionary history and influences of the 2R WGDs are still unclear, especially in noncoding regions. To address this issue, we identified paralogous CNSs. Region-focused Basic Local Alignment Search Tool (BLAST) search of each synteny block revealed 7,924 orthologous CNSs and 309 paralogous CNSs conserved among eight high-quality vertebrate genomes. Paralogous CNSs we found contained 115 previously reported ones and newly detected 194 ones. Through comparisons with VISTA Enhancer Browser and available ChIP-seq data, one-third (103) of paralogous CNSs detected in this study showed gene regulatory activity in the brain at several developmental stages. Their genomic locations are highly enriched near the transcription factor-coding regions, which are expressed in brain and neural systems. These results suggest that paralogous CNSs are conserved mainly because of maintaining gene expression in the vertebrate brain.

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.