Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome

Comparative genomics approach to detecting split-coding regions in a low-coverage genome: lessons from the chimaera Callorhinchus milii (Holocephali, Chondrichthyes)

Recent development of deep sequencing technologies has facilitated de novo genome sequencing projects, now conducted even by individual laboratories. However, this will yield more and more genome sequences that are not well assembled, and will hinder thorough annotation when no closely related reference genome is available. One of the challenging issues is the identification of protein-coding sequences split into multiple unassembled genomic segments, which can confound orthology assignment and various laboratory experiments requiring the identification of individual genes. In this study, using the genome of a cartilaginous fish, Callorhinchus milii, as test case, we performed gene prediction using a model specifically trained for this genome. We implemented an algorithm, designated ESPRIT, to identify possible linkages between multiple protein-coding portions derived from a single genomic locus split into multiple unassembled genomic segments. We developed a validation framework based on an artificially fragmented human genome, improvements between early and recent mouse genome assemblies, comparison with experimentally validated sequences from GenBank, and phylogenetic analyses. Our strategy provided insights into practical solutions for efficient annotation of only partially sequenced (low-coverage) genomes. To our knowledge, our study is the first formulation of a method to link unassembled genomic segments based on proteomes of relatively distantly related species as references.

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

BACKGROUND

In understanding the evolutionary process of vertebrates, cyclostomes (hagfishes and lamprey) occupy crucial positions. Resolving molecular phylogenetic relationships of cyclostome genes with gnathostomes (jawed vertebrates) genes is indispensable in deciphering both the species tree and gene trees. However, molecular phylogenetic analyses, especially those including lamprey genes, have produced highly discordant results between gene families. To efficiently scrutinize this problem using partial genome assemblies of early vertebrates, we focused on the potassium voltage-gated channel, shaker-related (KCNA) family, whose members are mostly single-exon.

RESULTS

Seven sea lamprey KCNA genes as well as six elephant shark genes were identified, and their orthologies to bony vertebrate subgroups were assessed. In contrast to robustly supported orthology of the elephant shark genes to gnathostome subgroups, clear orthology of any sea lamprey gene could not be established. Notably, sea lamprey KCNA sequences displayed unique codon usage pattern and amino acid composition, probably associated with exceptionally high GC-content in their coding regions. This lamprey-specific property of coding sequences was also observed generally for genes outside this gene family.

CONCLUSIONS

Our results suggest that secondary modifications of sequence properties unique to the lamprey lineage may be one of the factors preventing robust orthology assessments of lamprey genes, which deserves further genome-wide validation. The lamprey lineage-specific alteration of protein-coding sequence properties needs to be taken into consideration in tackling the key questions about early vertebrate evolution.

Hox clusters of the bichir (Actinopterygii, Polypterus senegalus) highlight unique patterns of sequence evolution in gnathostome phylogeny

Teleost fishes have extra Hox gene clusters owing to shared or lineage-specific genome duplication events in rayfinned fish (actinopterygian) phylogeny. Hence, extrapolating between genome function of teleosts and human or even between different fish species is difficult. We have sequenced and analyzed Hox gene clusters of the Senegal bichir (Polypterus senegalus), an extant representative of the most basal actinopterygian lineage. Bichir possesses four Hox gene clusters (A, B, C, D); phylogenetic analysis supports their orthology to the four Hox gene clusters of the gnathostome ancestor. We have generated a comprehensive database of conserved Hox noncoding sequences that include cartilaginous, lobe-finned, and ray-finned fishes (bichir and teleosts). Our analysis identified putative and known Hox cis-regulatory sequences with differing depths of conservation in Gnathostoma. We found that although bichir possesses four Hox gene clusters, its pattern of conservation of noncoding sequences is mosaic between outgroups, such as human, coelacanth, and shark, with four Hox gene clusters and teleosts, such as zebrafish and pufferfish, with seven or eight Hox gene clusters. Notably, bichir Hox gene clusters have been invaded by DNA transposons and this trend is further exemplified in teleosts, suggesting an as yet unrecognized mechanism of genome evolution that may explain Hox cluster plasticity in actinopterygians. Taken together, our results suggest that actinopterygian Hox gene clusters experienced a reduction in selective constraints that surprisingly predates the teleost-specific genome duplication.

Occurrence of two distinct urotensin II-related peptides in zebrafish provides new insight into the evolutionary history of the urotensin II gene family

The urotensin II (UII) family is currently known to consist of two paralogous peptides, namely UII and UII-related peptide (URP). In contrast to UII, which has been identified in all vertebrate classes so far, URP has only been characterized in tetrapods. We report here the occurrence of two distinct URP genes in teleosts, which we have named URP1 and URP2. Synteny analysis revealed that teleost URP1 and URP2 genes and tetrapod URP genes represent three distinct paralog genes that, together with the UII gene, probably arose from the two rounds of tetraploidization, which took place early in vertebrate evolution. The absence of URP in fish indicates that the corresponding gene has been lost in the teleost lineage, whereas it is likely that both the URP1 and URP2 genes have been lost in the tetrapod lineage. Quantitative RT-PCR analysis revealed that the URP2 gene is mainly expressed in the spinal cord and the brain in adult zebrafish. In situ hybridization experiments showed that in zebrafish embryos, URP2 mRNA-containing cells are located in the floor plate of the neural tube. In adult, URP2-expressing cells occur in close contact with the ventral side of the ependymal canal along the whole spinal cord, whereas in the brain, they are located below the fourth ventricle. These URP-expressing cells may correspond to cerebrospinal fluid-contacting neurons. In conclusion, our study reveals the occurrence of four distinct UII paralogous systems in vertebrates that may exert distinct functions, both in tetrapods and teleosts.

Identification of members of the gonadotropin-releasing hormone (GnRH), corticotropin-releasing factor (CRF) families in the genome of the holocephalan, Callorhinchus milii (elephant shark)

The gonadotropin-releasing hormone (GnRH) and corticotropin-releasing family (CRF) are two neuropeptides families that are strongly conserved throughout evolution. Recently, the genome of the holocephalan, Callorhinchus milii (elephant shark) has been sequenced. The phylogenetic position of C. milii, along with the relatively slow evolution of the cartilaginous fish suggests that neuropeptides in this species may resemble the earliest gnathostome forms. The genome of the elephant shark was screened, in silico, using the various conserved motifs of both the vertebrate CRF paralogs and the insect diuretic hormone sequences to identify the structure of the C. milii CRF/DH-like peptides. A similar approach was taken to identify the GnRH peptides using conserved motifs in both vertebrate and invertebrate forms. Two CRF peptides, a urotensin-1 peptide and a urocortin 3 peptide were found in the genome. There was only about 50% sequence identity between the two CRF peptides suggesting an early divergence. In addition, the urocortin 2 peptide seems to have been lost and was identified as a pseudogene in C. milii. In contrast to the number of CRF family peptides, only a GnRH-II preprohormone with the conserved mature decapeptide was found. This confirms early studies about the identity of GnRH in the Holocephali, and suggests that the Holocephali and Elasmobranchii differ with respect to GnRH structure and function.

Ion regulation in fish gills: recent progress in the cellular and molecular mechanisms

Fish encounter harsh ionic/osmotic gradients on their aquatic environments, and the mechanisms through which they maintain internal homeostasis are more challenging compared with those of terrestrial vertebrates. Gills are one of the major organs conducting the internal ionic and acid-base regulation, with specialized ionocytes as the major cells carrying out active transport of ions. Exploring the iono/osmoregulatory mechanisms in fish gills, extensive literature proposed several models, with many conflicting or unsolved issues. Recent studies emerged, shedding light on these issues with new opened windows on other aspects, on account of available advanced molecular/cellular physiological approaches and animal models. Respective types of ionocytes and ion transporters, and the relevant regulators for the mechanisms of NaCl secretion, Na(+) uptake/acid secretion/NH(4)(+) excretion, Ca(2+) uptake, and Cl(-) uptake/base secretion, were identified and functionally characterized. These new ideas broadened our understanding of the molecular/cellular mechanisms behind the functional modification/regulation of fish gill ion transport during acute and long-term acclimation to environmental challenges. Moreover, a model for the systematic and local carbohydrate energy supply to gill ionocytes during these acclimation processes was also proposed. These provide powerful platforms to precisely study transport pathways and functional regulation of specific ions, transporters, and ionocytes; however, very few model species were established so far, whereas more efforts are needed in other species.

Primordial linkage of β2-microglobulin to the MHC

β2-Microglobulin (β2M) is believed to have arisen in a basal jawed vertebrate (gnathostome) and is the essential L chain that associates with most MHC class I molecules. It contains a distinctive molecular structure called a constant-1 Ig superfamily domain, which is shared with other adaptive immune molecules including MHC class I and class II. Despite its structural similarity to class I and class II and its conserved function, β2M is encoded outside the MHC in all examined species from bony fish to mammals, but it is assumed to have translocated from its original location within the MHC early in gnathostome evolution. We screened a nurse shark bacterial artificial chromosome library and isolated clones containing β2M genes. A gene present in the MHC of all other vertebrates (ring3) was found in the bacterial artificial chromosome clone, and the close linkage of ring3 and β2M to MHC class I and class II genes was determined by single-strand conformational polymorphism and allele-specific PCR. This study satisfies the long-held conjecture that β2M was linked to the primordial MHC (Ur MHC); furthermore, the apparent stability of the shark genome may yield other genes predicted to have had a primordial association with the MHC specifically and with immunity in general.

Evolution of the gonadotropin-releasing hormone (GnRH) gene family in relation to vertebrate tetraploidizations

The neuropeptide gonadotropin-releasing hormone (GnRH) plays an important role in the control of reproductive functions. Vertebrates possess multiple GnRH isoforms that are classified into three main groups, namely GnRH1, GnRH2 and GnRH3. In the present study, we show that the chromosomal organization of the three GnRH loci is very well conserved among gnathostome species. We analyzed genes belonging to several other multigenic families that are present in the vicinity of GnRH genes. Five of them were seen to occur in four chromosomal regions that clearly form a paralogon. Moreover, we show that the homologous regions in the amphioxus genome are present on a single locus. Taken together, these observations indicate that GnRH1, GnRH2 and GnRH3 genes represent three paralogous genes that resulted from the two rounds of tetraploidization that took place early in vertebrate evolution. They confirm that the GnRH3 gene which is currently known only in teleost has most likely been lost in the tetrapod lineage. Finally, they suggest the existence of a fourth GnRH gene, named GnRH4. Whether the GnRH4 gene still exists in extant vertebrates is currently unknown. A search for this putative gene would be particularly useful in basal groups such as agnathans and cartilaginous fish.

Expressed 2-5A synthetase genes and pseudogenes in the marine sponge Geodia barretti

The 2',5'-oligoadenylate synthetases (2-5A synthetases, OAS) form a family of proteins presented in many branches of Metazoa. The phylum Porifera (sponges) contains OAS proteins which are different from those in vertebrates and form a distinct OAS subfamily. In turn, OAS proteins from different genera of Demospongia show rather low similarities in their primary structures. To ascertain divergence of the OAS genes within a particular sponge genus, we identified the OAS genes from the marine sponge Geodia barretti and compared them with those from another member of the genus Geodia, Geodia cydonium. The identity and similarity of the OAS sequences found in G. barretti with those from G. cydonium were considerably higher than identities and similarities compared with those from other sponges, 75% and 85% versus 27-30% and 42-47%, respectively. We also established the presence of a transcriptionally active polymorphic OAS pseudogene in the genome of G. barretti. The transcripts of the OAS pseudogene(s) lack several internal exons encoding necessary motifs for OAS enzymatic activity. The maintenance and further diversification of OAS gene(s) and pseudogene(s) suggest the prevalence of gene duplication events over the loss of gene duplicates in Geodia genomes during the evolution.

Holocephalan embryos provide evidence for gill arch appendage reduction and opercular evolution in cartilaginous fishes

Chondrichthyans possess endoskeletal appendages called branchial rays that extend laterally from their hyoid and gill-bearing (branchial) arches. Branchial ray outgrowth, like tetrapod limb outgrowth, is maintained by Sonic hedgehog (Shh) signaling. In limbs, distal endoskeletal elements fail to form in the absence of normal Shh signaling, whereas shortened duration of Shh expression correlates with distal endoskeletal reduction in naturally variable populations. Chondrichthyans also exhibit natural variation with respect to branchial ray distribution--elasmobranchs (sharks and batoids) possess a series of ray-supported septa on their hyoid and gill arches, whereas holocephalans (chimaeras) possess a single hyoid arch ray-supported operculum. Here we show that the elongate hyoid rays of the holocephalan Callorhinchus milii grow in association with sustained Shh expression within an opercular epithelial fold, whereas Shh is only transiently expressed in the gill arches. Coincident with this transient Shh expression, branchial ray outgrowth is initiated in C. milii but is not maintained, yielding previously unrecognized vestigial gill arch branchial rays. This is in contrast to the condition seen in sharks, where sustained Shh expression corresponds to the presence of fully formed branchial rays on the hyoid and gill arches. Considered in light of current hypotheses of chondrichthyan phylogeny, our data suggest that the holocephalan operculum evolved in concert with gill arch appendage reduction by attenuation of Shh-mediated branchial ray outgrowth, and that chondrichthyan branchial rays and tetrapod limbs exhibit parallel developmental mechanisms of evolutionary reduction.

Ancestral and derived attributes of the dlx gene repertoire, cluster structure and expression patterns in an African cichlid fish

Background

Cichlid fishes have undergone rapid, expansive evolutionary radiations that are manifested in the diversification of their trophic morphologies, tooth patterning and coloration. Understanding the molecular mechanisms that underlie the cichlids' unique patterns of evolution requires a thorough examination of genes that pattern the neural crest, from which these diverse phenotypes are derived. Among those genes, the homeobox-containing Dlx gene family is of particular interest since it is involved in the patterning of the brain, jaws and teeth.

Results

In this study, we characterized the dlx genes of an African cichlid fish, Astatotilapia burtoni, to provide a baseline to later allow cross-species comparison within Cichlidae. We identified seven dlx paralogs (dlx1a, -2a, -4a, -3b, -4b, -5a and -6a), whose orthologies were validated with molecular phylogenetic trees. The intergenic regions of three dlx gene clusters (dlx1a-2a, dlx3b-4b, and dlx5a-6a) were amplified with long PCR. Intensive cross-species comparison revealed a number of conserved non-coding elements (CNEs) that are shared with other percomorph fishes. This analysis highlighted additional lineage-specific gains/losses of CNEs in different teleost fish lineages and a novel CNE that had previously not been identified. Our gene expression analyses revealed overlapping but distinct expression of dlx orthologs in the developing brain and pharyngeal arches. Notably, four of the seven A. burtoni dlx genes, dlx2a, dlx3b, dlx4a and dlx5a, were expressed in the developing pharyngeal teeth.

Conclusion

This comparative study of the dlx genes of A. burtoni has deepened our knowledge of the diversity of the Dlx gene family, in terms of gene repertoire, expression patterns and non-coding elements. We have identified possible cichlid lineage-specific changes, including losses of a subset of dlx expression domains in the pharyngeal teeth, which will be the targets of future functional studies.

The early origin of melanocortin receptors, agouti-related peptide, agouti signalling peptide, and melanocortin receptor-accessory proteins, with emphasis on pufferfishes, elephant shark, lampreys, and amphioxus

There are conflicting theories about the evolution of melanocortin MC receptors while only few studies have addressed the evolution of agouti-related peptide (AgRP) and agouti signalling peptide (ASIP), which are antagonists at the melanocortin receptors (MCRs), or the melanocortin MC(2) receptor accessory proteins (MRAP1 and MRAP2). Previously we have cloned melanocortin MC receptors (MC(a) and MC(b)) genes in river lamprey and here we identify orthologues to these melanocortin MC receptor sequences in the sea lamprey. We investigate the putative presence of the melanocortin MC receptor genes in lancelet (amphioxus; Branchiostoma floridae) but we find it unlikely that such gene exists, due to a sharp drop in sequence similarity beyond sequence clusters of known receptors. We show the presence of AgRP and ASIP in elephant shark, a cartilaginous fish belonging to the subclass of Elasmobranchii. However, we do not find any of these genes in lamprey or lancelet after detailed analysis of both targeted and whole proteome regular expression scans. We found MRAP2, but not MRAP1, to be present in elephant shark and sea lamprey while Fugu (T. rubripes) has both genes. This study shows that the most ancient presence of these melanocortin-related sequences is found in elephant shark and lampreys considering the current available sequence data.

Evolution of the parathyroid hormone family and skeletal formation pathways

Bone is considered to be a feature of higher vertebrates and one of the features that was required for the movement from water onto land. But there are a number of evolutionarily important species that have cartilaginous skeletons, including sharks. Both bony and cartilaginous fish are believed to have a common ancestor who had a bony skeleton. A number of factors and pathways have been shown to be involved in the development and maintenance of bony skeleton including the Wnt pathway and the parathyroid hormone gene family. The study of these pathways and factors in cartilaginous animals may shed light on the evolution of the vertebrate skeleton.

The vertebrate makorin ubiquitin ligase gene family has been shaped by large-scale duplication and retroposition from an ancestral gonad-specific, maternal-effect gene

BACKGROUND

Members of the makorin (mkrn) gene family encode RING/C3H zinc finger proteins with U3 ubiquitin ligase activity. Although these proteins have been described in a variety of eukaryotes such as plants, fungi, invertebrates and vertebrates including human, almost nothing is known about their structural and functional evolution.

RESULTS

Via partial sequencing of a testis cDNA library from the poeciliid fish Xiphophorus maculatus, we have identified a new member of the makorin gene family, that we called mkrn4. In addition to the already described mkrn1 and mkrn2, mkrn4 is the third example of a makorin gene present in both tetrapods and ray-finned fish. However, this gene was not detected in mouse and rat, suggesting its loss in the lineage leading to rodent murids. Mkrn2 and mkrn4 are located in large ancient duplicated regions in tetrapod and fish genomes, suggesting the possible involvement of ancestral vertebrate-specific genome duplication in the formation of these genes. Intriguingly, many mkrn1 and mkrn2 intronless retrocopies have been detected in mammals but not in other vertebrates, most of them corresponding to pseudogenes. The nature and number of zinc fingers were found to be conserved in Mkrn1 and Mkrn2 but much more variable in Mkrn4, with lineage-specific differences. RT-qPCR analysis demonstrated a highly gonad-biased expression pattern for makorin genes in medaka and zebrafish (ray-finned fishes) and amphibians, but a strong relaxation of this specificity in birds and mammals. All three mkrn genes were maternally expressed before zygotic genome activation in both medaka and zebrafish early embryos.

CONCLUSION

Our analysis demonstrates that the makorin gene family has evolved through large-scale duplication and subsequent lineage-specific retroposition-mediated duplications in vertebrates. From the three major vertebrate mkrn genes, mkrn4 shows the highest evolutionary dynamics, with lineage-specific loss of zinc fingers and even complete gene elimination from certain groups of vertebrates. Comparative expression analysis strongly suggests that the ancestral E3 ubiquitin ligase function of the single copy mkrn gene before duplication in vertebrates was gonad-specific, with maternal expression in early embryos.

Parathyroid hormone gene family in a cartilaginous fish, the elephant shark (Callorhinchus milii)

The development of bone was a major step in the evolution of vertebrates. A bony skeleton provided structural support and a calcium reservoir essential for the movement from an aquatic to a terrestrial environment. Cartilaginous fishes are the oldest living group of jawed vertebrates. In this study we have identified three members of the parathyroid hormone (Pth) gene family in a cartilaginous fish, the elephant shark (Callorhinchus milii). The three genes include two Pth genes, designated as Pth1 and Pth2, and a Pthrp gene. Phylogenetic analysis suggested that elephant shark Pth2 is an ancient gene whose orthologue is lost in bony vertebrates. The Pth1 and Pth2 genes have the same structure as the Pth gene in bony vertebrates, whereas the structure of the Pthrp gene is more complex in tetrapods compared with elephant shark. The three elephant shark genes showed distinct patterns of expression, with Pth2 being expressed only in the brain and spleen. This contrasts with localization of the corresponding proteins, which showed considerable overlap in their distribution. There were conserved sites of localization for Pthrp between elephant shark and mammals, including tissues such as kidney, skin, skeletal and cardiac muscle, pancreas, and cartilage. The elephant shark Pth1(1-34) and Pthrp(1-34) peptides were able to stimulate cAMP accumulation in mammalian UMR106.01 cells. However, Pth2(1-34) peptide did not show such PTH-like biologic activity. The presence of Pth and Pthrp genes in the elephant shark indicates that these genes played fundamental roles before their recruitment to bone development in bony jawed vertebrates.

Differential expression of Na+-Cl- cotransporter and Na+-K+-Cl- cotransporter 2 in the distal nephrons of euryhaline and seawater pufferfishes

The process of NaCl reabsorption in the distal nephron allows freshwater fishes to excrete hypotonic urine and seawater fishes to excrete urine containing high concentrations of divalent ions; the relevant transporters, however, have not yet been identified. In the mammalian distal nephron, NaCl absorption is mediated by Na(+)-K(+)-Cl(-) cotransporter 2 (NKCC2, Slc12a1) in the thick ascending limb, Na(+)-Cl(-) cotransporter (NCC, Slc12a3) in the distal convoluted tubule, and epithelial sodium channel (ENaC) in the collecting duct. In this study, we compared the expression profiles of these proteins in the kidneys of euryhaline and seawater pufferfishes. Mining the fugu genome identified one NKCC2 gene and one NCC gene, but no ENaC gene. RT-PCR and in situ hybridization analyses demonstrated that NKCC2 was highly expressed in the distal tubules and NCC was highly expressed in the collecting ducts of euryhaline pufferfish (mefugu, Takifugu obscurus). On the other hand, the kidney of seawater pufferfish (torafugu, Takifugu rubripes), which lacked distal tubules, expressed very low levels of NCC, and, in the collecting ducts, high levels of NKCC2. Acclimation of mefugu to seawater resulted in a 2.7× decrease in NCC expression, whereas NKCC2 expression was not markedly affected. Additionally, internalization of NCC from the apical surface of the collecting ducts was observed. These results suggest that NaCl reabsorption in the distal nephron of the fish kidney is mediated by NCC and NKCC2 in freshwater and by NKCC2 in seawater.

Ancient vertebrate conserved noncoding elements have been evolving rapidly in teleost fishes

Vertebrate genomes contain thousands of conserved noncoding elements (CNEs) that often function as tissue-specific enhancers. In this study, we have identified CNEs in human, dog, chicken, Xenopus, and four teleost fishes (zebrafish, stickleback, medaka, and fugu) using elephant shark, a cartilaginous vertebrate, as the base genome and investigated the evolution of these ancient vertebrate CNEs (aCNEs) in bony vertebrate lineages. Our analysis shows that aCNEs have been evolving at different rates in different bony vertebrate lineages. Although 78-83% of CNEs have diverged beyond recognition ("lost") in different teleost fishes, only 24% and 40% have been lost in the chicken and mammalian lineages, respectively. Relative rate tests of substitution rates in CNEs revealed that the teleost fish CNEs have been evolving at a significantly higher rate than those in other bony vertebrates. In the ray-finned fish lineage, 68% of aCNEs were lost before the divergence of the four teleosts. This implicates the "fish-specific" whole-genome duplication in the accelerated evolution and the loss of a large number of both copies of duplicated CNEs in teleost fishes. The aCNEs are rich in tissue-specific enhancers and thus many of them are likely to be evolutionarily constrained cis-regulatory elements. The rapid evolution of aCNEs might have affected the expression patterns driven by them. Transgenic zebrafish assay of some human CNE enhancers that have been lost in teleosts has indicated instances of conservation or changes in trans-acting factors between mammals and fishes.

Human brain evolution: harnessing the genomics (r)evolution to link genes, cognition, and behavior

The evolution of the human brain has resulted in numerous specialized features including higher cognitive processes such as language. Knowledge of whole-genome sequence and structural variation via high-throughput sequencing technology provides an unprecedented opportunity to view human evolution at high resolution. However, phenotype discovery is a critical component of these endeavors and the use of nontraditional model organisms will also be critical for piecing together a complete picture. Ultimately, the union of developmental studies of the brain with studies of unique phenotypes in a myriad of species will result in a more thorough model of the groundwork the human brain was built upon. Furthermore, these integrative approaches should provide important insights into human diseases.

Parallel retention of Pdx2 genes in cartilaginous fish and coelacanths

The Pdx1 or Ipf1 gene encodes an important homeodomain-containing protein with key roles in pancreas development and function. Mutations in human PDX1 are implicated in developmental defects and disease of the pancreas. Extensive research, including genome sequencing, has indicated that Pdx1 is the only member of its gene family in mammals, birds, amphibians, and ray-finned fish, and with the exception of teleost fish, this gene forms part of the ParaHox gene cluster along with Gsx1 and Cdx2. The ParaHox cluster, however, is a remnant of a 4-fold genome duplication; the three other ParaHox paralogues lack a Pdx-like gene in all vertebrate genomes examined to date. We have used bacterial artificial chromosome cloning and synteny analysis to show that the ancestor of living jawed vertebrates in fact had more ParaHox genes, including two Pdx genes (Pdx1 and Pdx2). Surprisingly, the two Pdx genes have been retained in parallel in two quite distantly related lineages, the cartilaginous fish (sharks, skates, and chimeras) and the Indonesian coelacanth, Latimeria menadoensis. The Pdx2 gene has been lost independently in ray-finned fish and in tetrapods.

Cloning and functional characterization of Chondrichthyes, cloudy catshark, Scyliorhinus torazame and whale shark, Rhincodon typus estrogen receptors

Sex-steroid hormones are essential for normal reproductive activity in both sexes in all vertebrates. Estrogens are required for ovarian differentiation during a critical developmental stage and promote the growth and differentiation of the female reproductive system following puberty. Recent studies have shown that environmental estrogens influence the developing reproductive system as well as gametogenesis, especially in males. To understand the molecular mechanisms of estrogen actions and to evaluate estrogen receptor-ligand interactions in Elasmobranchii, we cloned a single estrogen receptor (ESR) from two shark species, the cloudy catshark (Scyliorhinus torazame) and whale shark (Rhincodon typus) and used an ERE-luciferase reporter assay system to characterize the interaction of these receptors with steroidal and other environmental estrogens. In the transient transfection ERE-luciferase reporter assay system, both shark ESR proteins displayed estrogen-dependent activation of transcription, and shark ESRs were more sensitive to 17beta-estradiol compared with other natural and synthetic estrogens. Further, the environmental chemicals, bisphenol A, nonylphenol, octylphenol and DDT could activate both shark ESRs. The assay system provides a tool for future studies examining the receptor-ligand interactions and estrogen disrupting mechanisms in Elasmobranchii.

Phylogenetic analysis of ferlin genes reveals ancient eukaryotic origins

BACKGROUND

The ferlin gene family possesses a rare and identifying feature consisting of multiple tandem C2 domains and a C-terminal transmembrane domain. Much currently remains unknown about the fundamental function of this gene family, however, mutations in its two most well-characterised members, dysferlin and otoferlin, have been implicated in human disease. The availability of genome sequences from a wide range of species makes it possible to explore the evolution of the ferlin family, providing contextual insight into characteristic features that define the ferlin gene family in its present form in humans.

RESULTS

Ferlin genes were detected from all species of representative phyla, with two ferlin subgroups partitioned within the ferlin phylogenetic tree based on the presence or absence of a DysF domain. Invertebrates generally possessed two ferlin genes (one with DysF and one without), with six ferlin genes in most vertebrates (three DysF, three non-DysF). Expansion of the ferlin gene family is evident between the divergence of lamprey (jawless vertebrates) and shark (cartilaginous fish). Common to almost all ferlins is an N-terminal C2-FerI-C2 sandwich, a FerB motif, and two C-terminal C2 domains (C2E and C2F) adjacent to the transmembrane domain. Preservation of these structural elements throughout eukaryotic evolution suggests a fundamental role of these motifs for ferlin function. In contrast, DysF, C2DE, and FerA are optional, giving rise to subtle differences in domain topologies of ferlin genes. Despite conservation of multiple C2 domains in all ferlins, the C-terminal C2 domains (C2E and C2F) displayed higher sequence conservation and greater conservation of putative calcium binding residues across paralogs and orthologs. Interestingly, the two most studied non-mammalian ferlins (Fer-1 and Misfire) in model organisms C. elegans and D. melanogaster, present as outgroups in the phylogenetic analysis, with results suggesting reproduction-related divergence and specialization of species-specific functions within their genus.

CONCLUSIONS

Our phylogenetic studies provide evolutionary insight into the ferlin gene family. We highlight the existence of ferlin-like proteins throughout eukaryotic evolution, from unicellular phytoplankton and apicomplexan parasites, through to humans. We characterise the preservation of ferlin structural motifs, not only of C2 domains, but also the more poorly characterised ferlin-specific motifs representing the DysF, FerA and FerB domains. Our data suggest an ancient role of ferlin proteins, with lessons from vertebrate biology and human disease suggesting a role relating to vesicle fusion and plasma membrane specialization.

Widespread presence of human BOULE homologs among animals and conservation of their ancient reproductive function

Sex-specific traits that lead to the production of dimorphic gametes, sperm in males and eggs in females, are fundamental for sexual reproduction and accordingly widespread among animals. Yet the sex-biased genes that underlie these sex-specific traits are under strong selective pressure, and as a result of adaptive evolution they often become divergent. Indeed out of hundreds of male or female fertility genes identified in diverse organisms, only a very small number of them are implicated specifically in reproduction in more than one lineage. Few genes have exhibited a sex-biased, reproductive-specific requirement beyond a given phylum, raising the question of whether any sex-specific gametogenesis factors could be conserved and whether gametogenesis might have evolved multiple times. Here we describe a metazoan origin of a conserved human reproductive protein, BOULE, and its prevalence from primitive basal metazoans to chordates. We found that BOULE homologs are present in the genomes of representative species of each of the major lineages of metazoans and exhibit reproductive-specific expression in all species examined, with a preponderance of male-biased expression. Examination of Boule evolution within insect and mammalian lineages revealed little evidence for accelerated evolution, unlike most reproductive genes. Instead, purifying selection was the major force behind Boule evolution. Furthermore, loss of function of mammalian Boule resulted in male-specific infertility and a global arrest of sperm development remarkably similar to the phenotype in an insect boule mutation. This work demonstrates the conservation of a reproductive protein throughout eumetazoa, its predominant testis-biased expression in diverse bilaterian species, and conservation of a male gametogenic requirement in mice. This shows an ancient gametogenesis requirement for Boule among Bilateria and supports a model of a common origin of spermatogenesis.

Loss of fish actinotrichia proteins and the fin-to-limb transition

The early development of teleost paired fins is strikingly similar to that of tetrapod limb buds and is controlled by similar mechanisms. One early morphological divergence between pectoral fins and limbs is in the fate of the apical ectodermal ridge (AER), the distal epidermis that rims the bud. Whereas the AER of tetrapods regresses after specification of the skeletal progenitors, the AER of teleost fishes forms a fold that elongates. Formation of the fin fold is accompanied by the synthesis of two rows of rigid, unmineralized fibrils called actinotrichia, which keep the fold straight and guide the migration of mesenchymal cells within the fold. The actinotrichia are made of elastoidin, the components of which, apart from collagen, are unknown. Here we show that two zebrafish proteins, which we name actinodin 1 and 2 (And1 and And2), are essential structural components of elastoidin. The presence of actinodin sequences in several teleost fishes and in the elephant shark (Callorhinchus milii, which occupies a basal phylogenetic position), but not in tetrapods, suggests that these genes have been lost during tetrapod species evolution. Double gene knockdown of and1 and and2 in zebrafish embryos results in the absence of actinotrichia and impaired fin folds. Gene expression profiles in embryos lacking and1 and and2 function are consistent with pectoral fin truncation and may offer a potential explanation for the polydactyly observed in early tetrapod fossils. We propose that the loss of both actinodins and actinotrichia during evolution may have led to the loss of lepidotrichia and may have contributed to the fin-to-limb transition.

Evolutionarily conserved TCR binding sites, identification of T cells in primary lymphoid tissues, and surprising trans-rearrangements in nurse shark

Cartilaginous fish are the oldest animals that generate RAG-based Ag receptor diversity. We have analyzed the genes and expressed transcripts of the four TCR chains for the first time in a cartilaginous fish, the nurse shark (Ginglymostoma cirratum). Northern blotting found TCR mRNA expression predominantly in lymphoid and mucosal tissues. Southern blotting suggested translocon-type loci encoding all four chains. Based on diversity of V and J segments, the expressed combinatorial diversity for gamma is similar to that of human, alpha and beta may be slightly lower, and delta diversity is the highest of any organism studied to date. Nurse shark TCRdelta have long CDR3 loops compared with the other three chains, creating binding site topologies comparable to those of mammalian TCR in basic paratope structure; additionally, nurse shark TCRdelta CDR3 are more similar to IgH CDR3 in length and heterogeneity than to other TCR chains. Most interestingly, several cDNAs were isolated that contained IgM or IgW V segments rearranged to other gene segments of TCRdelta and alpha. Finally, in situ hybridization experiments demonstrate a conservation of both alpha/beta and gamma/delta T cell localization in the thymus across 450 million years of vertebrate evolution, with gamma/delta TCR expression especially high in the subcapsular region. Collectively, these data make the first cellular identification of TCR-expressing lymphocytes in a cartilaginous fish.

Evolutionary origin and phylogeny of the modern holocephalans (Chondrichthyes: Chimaeriformes): a mitogenomic perspective

With our increasing ability for generating whole-genome sequences, comparative analysis of whole genomes has become a powerful tool for understanding the structure, function, and evolutionary history of human and other vertebrate genomes. By virtue of their position basal to bony vertebrates, cartilaginous fishes (class Chondrichthyes) are a valuable outgroup in comparative studies of vertebrates. Recently, a holocephalan cartilaginous fish, the elephant shark, Callorhinchus milii (Subclass Holocephali: Order Chimaeriformes), has been proposed as a model genome, and low-coverage sequence of its genome has been generated. Despite such an increasing interest, the evolutionary history of the modern holocephalans-a previously successful and diverse group but represented by only 39 extant species-and their relationship with elasmobranchs and other jawed vertebrates has been poorly documented largely owing to a lack of well-preserved fossil materials after the end-Permian about 250 Ma. In this study, we assembled the whole mitogenome sequences for eight representatives from all the three families of the modern holocephalans and investigated their phylogenetic relationships and evolutionary history. Unambiguously aligned sequences from these holocephalans together with 17 other vertebrates (9,409 nt positions excluding entire third codon positions) were subjected to partitioned maximum likelihood analysis. The resulting tree strongly supported a single origin of the modern holocephalans and their sister-group relationship with elasmobranchs. The mitogenomic tree recovered the most basal callorhinchids within the chimaeriforms, which is sister to a clade comprising the remaining two families (rhinochimaerids and chimaerids). The timetree derived from a relaxed molecular clock Bayesian method suggests that the holocephalans originated in the Silurian about 420 Ma, having survived from the end-Permian (250 Ma) mass extinction and undergoing familial diversifications during the late Jurassic to early Cretaceous (170-120 Ma). This postulated evolutionary scenario agrees well with that based on the paleontological observations.

Functional conservation of a forebrain enhancer from the elephant shark (Callorhinchus milii ) in zebrafish and mice

Background

The phylogenetic position of the elephant shark (Callorhinchus milii ) is particularly relevant to study the evolution of genes and gene regulation in vertebrates. Here we examine the evolution of Dlx homeobox gene regulation during vertebrate embryonic development with a particular focus on the forebrain. We first identified the elephant shark sequence orthologous to the URE2 cis -regulatory element of the mouse Dlx1/Dlx2 locus (herein named CmURE2). We then conducted a comparative study of the sequence and enhancer activity of CmURE2 with that of orthologous regulatory sequences from zebrafish and mouse.

Results

The CmURE2 sequence shows a high percentage of identity with its mouse and zebrafish counterparts but is overall more similar to mouse URE2 (MmURE2) than to zebrafish URE2 (DrURE2). In transgenic zebrafish and mouse embryos, CmURE2 displayed enhancer activity in the forebrain that overlapped with that of DrURE2 and MmURE2. However, we detected notable differences in the activity of the three sequences in the diencephalon. Outside of the forebrain, CmURE2 shows enhancer activity in areas such as the pharyngeal arches and dorsal root ganglia where its' counterparts are also active.

Conclusions

Our transgenic assays show that part of the URE2 enhancer activity is conserved throughout jawed vertebrates but also that new characteristics have evolved in the different groups. Our study demonstrates that the elephant shark is a useful outgroup to study the evolution of regulatory mechanisms in vertebrates and to address how changes in the sequence of cis -regulatory elements translate into changes in their regulatory activity.

Leu85 in the beta1-beta2 linker of ASIC1 slows activation and decreases the apparent proton affinity by stabilizing a closed conformation

Acid-sensing ion channels (ASICs) are proton-activated channels expressed in neurons of the central and peripheral nervous systems where they modulate neuronal activity in response to external increases in proton concentration. The size of ASIC1 currents evoked by a given local acidification is determined by the number of channels in the plasma membrane and by the apparent proton affinities for activation and steady-state desensitization of the channel. Thus, the magnitude of the pH drop and the value of the baseline pH both are functionally important. Recent characterization of ASIC1s from an increasing number of species has made evident that proton affinities of these channels vary across vertebrates. We found that in species with high baseline plasma pH, e.g. frog, shark, and fish, ASIC1 has high proton affinity compared with the mammalian channel. The beta1-beta2 linker in the extracellular domain, specifically by the substitution M85L, determines the interspecies differences in proton affinities and also the time course of ASIC1 macroscopic currents. The mechanism underlying these observations is a delay in channel opening after application of protons, most likely by stabilizing a closed conformation that decreases the apparent affinity to protons and also slows the rise and decay phases of the current. Together, the results suggest evolutionary adaptation of ASIC1 to match the value of the species-specific plasma pH. At the molecular level, adaptation is achieved by substitutions of nonionizable residues rather than by modification of the channel proton sensor.

Phylogenetic analysis of the vertebrate excitatory/neutral amino acid transporter (SLC1/EAAT) family reveals lineage specific subfamilies

Background

The composition and expression of vertebrate gene families is shaped by species specific gene loss in combination with a number of gene and genome duplication events (R1, R2 in all vertebrates, R3 in teleosts) and depends on the ecological and evolutionary context. In this study we analyzed the evolutionary history of the solute carrier 1 (SLC1) gene family. These genes are supposed to be under strong selective pressure (purifying selection) due to their important role in the timely removal of glutamate at the synapse.

Results

In a genomic survey where we manually annotated and analyzing sequences from more than 300 SLC1 genes (from more than 40 vertebrate species), we found evidence for an interesting evolutionary history of this gene family. While human and mouse genomes contain 7 SLC1 genes, in prototheria, sauropsida, and amphibia genomes up to 9 and in actinopterygii up to 13 SLC1 genes are present. While some of the additional slc1 genes in ray-finned fishes originated from R3, the increased number of SLC1 genes in prototheria, sauropsida, and amphibia genomes originates from specific genes retained in these lineages.

Phylogenetic comparison and microsynteny analyses of the SLC1 genes indicate, that theria genomes evidently lost several SLC1 genes still present in the other lineage. The genes lost in theria group into two new subfamilies of the slc1 gene family which we named slc1a8/eaat6 and slc1a9/eaat7.

Conclusions

The phylogeny of the SLC1/EAAT gene family demonstrates how multiple genome reorganization and duplication events can influence the number of active genes. Inactivation and preservation of specific SLC1 genes led to the complete loss of two subfamilies in extant theria, while other vertebrates have retained at least one member of two newly identified SLC1 subfamilies.

The evolution of thrombospondins and their ligand-binding activities

The extracellular matrix (ECM) is a complex, multiprotein network that has essential roles in tissue integrity and intercellular signaling in the metazoa. Thrombospondins (TSPs) are extracellular, calcium-binding glycoproteins that have biologically important roles in mammals in angiogenesis, vascular biology, connective tissues, immune response, and synaptogenesis. The evolution of these complex functional properties is poorly understood. We report here on the evolution of TSPs and their ligand-binding capacities, from comparative genomics of species representing the major phyla of metazoa and experimental analyses of the oligomerization properties of noncanonical TSPs of basal deuterostomes. Monomeric, dimeric, trimeric, and pentameric TSPs have arisen through separate evolutionary events involving gain, loss, or modification of a coiled-coil domain or distinct domains at the amino-terminus. The relative transience of monomeric forms under evolution implicates a biological importance for multivalency of the C-terminal region of TSPs. Most protostomes have a single TSP gene encoding a pentameric TSP. The pentameric form is also present in deuterostomes, and gene duplications at the origin of deuterostomes and gene loss and further gene duplication events in the vertebrate lineage gave rise to distinct forms and novel domain architectures. Parallel analysis of the major ligands of mammalian TSPs revealed that many binding activities are neofunctions representing either coevolutionary innovations in the deuterostome lineage or neofunctions of ancient molecules such as CD36. Contrasting widely conserved capacities include binding to heparan glycosaminoglycans, fibrillar collagen, or RGD-dependent integrins. These findings identify TSPs as fundamental components of the extracellular interaction systems of metazoa and thus impact understanding of the evolution of ECM networks. The widely conserved activities of TSPs in binding to ECM components or PS2 clade integrins will be relevant to use of TSPs in synthetic extracellular matrices or tissue engineering. In contrast, the neofunctions of vertebrate TSPs likely include interactions suitable for therapeutic targeting without general disruption of ECM.

Stage-specific gene expression during spermatogenesis in the dogfish (Scyliorhinus canicula)

In the dogfish testis, the cystic arrangement and polarization of germ cell stages make it possible to observe all stages of spermatogenesis in a single transverse section. By taking advantage of the zonation of this organ, we have used suppressive subtractive libraries construction, real-time PCR, and in situ hybridization to identify 32 dogfish genes showing differential expressions during spermatogenesis. These include homologs of genes already known to be expressed in the vertebrate testis, but found here to be specifically expressed either in pre-meiotic and/or meiotic zones (ribosomal protein S8, high-mobility group box 3, ubiquitin carboxyl-terminal esterase L3, 20beta-hydroxysteroid dehydrogenase, or cyclophilin B) or in post-meiotic zone (speriolin, Soggy, zinc finger protein 474, calreticulin, or phospholipase c-zeta). We also report, for the first time, testis-specific expression patterns for dogfish genes coding for A-kinase anchor protein 5, ring finger protein 152, or F-box only protein 7. Finally, the study highlights the differential expression of new sequences whose identity remains to be assessed. This study provides the first molecular characterization of spermatogenesis in a chondrichthyan, a key species to gain insight into the evolution of this process in gnathostomes.

Elasmobranch qPCR reference genes: a case study of hypoxia preconditioned epaulette sharks

Background

Elasmobranch fishes are an ancient group of vertebrates which have high potential as model species for research into evolutionary physiology and genomics. However, no comparative studies have established suitable reference genes for quantitative PCR (qPCR) in elasmobranchs for any physiological conditions. Oxygen availability has been a major force shaping the physiological evolution of vertebrates, especially fishes. Here we examined the suitability of 9 reference candidates from various functional categories after a single hypoxic insult or after hypoxia preconditioning in epaulette shark (Hemiscyllium ocellatum).

Results

Epaulette sharks were caught and exposed to hypoxia. Tissues were collected from 10 controls, 10 individuals with single hypoxic insult and 10 individuals with hypoxia preconditioning (8 hypoxic insults, 12 hours apart). We produced sequence information for reference gene candidates and monitored mRNA expression levels in four tissues: cerebellum, heart, gill and eye. The stability of the genes was examined with analysis of variance, geNorm and NormFinder. The best ranking genes in our study were eukaryotic translation elongation factor 1 beta (eef1b), ubiquitin (ubq) and polymerase (RNA) II (DNA directed) polypeptide F (polr2f). The performance of the ribosomal protein L6 (rpl6) was tissue-dependent. Notably, in one tissue the analysis of variance indicated statistically significant differences between treatments for genes that were ranked as the most stable candidates by reference gene software.

Conclusions

Our results indicate that eef1b and ubq are generally the most suitable reference genes for the conditions and tissues in the present epaulette shark studies. These genes could also be potential reference gene candidates for other physiological studies examining stress in elasmobranchs. The results emphasise the importance of inter-group variation in reference gene evaluation.

Exon-primed intron-crossing (EPIC) markers for non-model teleost fishes

Background

Exon-primed intron-crossing (EPIC) markers have three advantages over anonymous genomic sequences in studying evolution of natural populations. First, the universal primers designed in exon regions can be applied across a broad taxonomic range. Second, the homology of EPIC-amplified sequences can be easily determined by comparing either their exon or intron portion depending on the genetic distance between the taxa. Third, having both the exon and intron fragments could help in examining genetic variation at the intraspecific and interspecific level simultaneously, particularly helpful when studying species complex. However, the paucity of EPIC markers has hindered multilocus studies using nuclear gene sequences, particularly in teleost fishes.

Results

We introduce a bioinformatics pipeline for developing EPIC markers by comparing the whole genome sequences between two or more species. By applying this approach on five teleost fishes whose genomes were available in the Ensembl database http://www.ensembl.org, we identified 210 EPIC markers that have single-copy and conserved exon regions with identity greater than 85% among the five teleost fishes. We tested 12 randomly chosen EPIC markers in nine teleost species having a wide phylogenetic range. The success rate of amplifying and sequencing those markers varied from 44% to 100% in different species. We analyzed the exon sequences of the 12 EPIC markers from 13 teleosts. The resulting phylogeny contains many traditionally well-supported clades, indicating the usefulness of the exon portion of EPIC markers in reconstructing species phylogeny, in addition to the value of the intron portion of EPIC markers in interrogating the population history.

Conclusions

This study illustrated an effective approach to develop EPIC markers in a taxonomic group, where two or more genome sequences are available. The markers identified could be amplified across a broad taxonomic range of teleost fishes. The phylogenetic utility of individual markers varied according to intron size and amplifiability. The bioinformatics pipelines developed are readily adapted to other taxonomic groups.

The evolution of gnathostome development: Insight from chondrichthyan embryology

Chondrichthyans (cartilaginous fishes) represent one of the two lineages of gnathostomes, the other being the osteicthyans (bony fishes). Classical studies on chondrichthyan embryology have strongly impacted our views of vertebrate body plan evolution, while recent studies highlight oviparous chondrichthyans as emerging vertebrate model systems that are amenable to experimental embryological manipulation. Here, we review three particular areas of interest in the field of chondrichthyan developmental biology-gastrulation, neural development, and appendage patterning-and we discuss recent findings within a broader chondrichthyan-osteichthyan comparative framework. In some cases, comparative studies of chondrichthyan and osteichthyan development reveal conserved patterns of gene expression in common developmental contexts. Studies of chondrichthyan gastrulation reveal conserved patterns of developmental gene expression, despite highly divergent modes of mesendoderm internalization, while molecular characterization of chondrichthyan neurogenic placodes indicates a conservation of placode transcription factor expression across gnathostome phylogeny. In other cases, comparative studies of chondrichthyan and osteichthyan development yield evidence of shared patterning mechanisms functioning in different developmental contexts. This is exemplified by studies on the development of chondrichthyan appendages-paired fins, median fins, and gill rays. These have demonstrated that a retinoic acid-responsive Shh-expressing signaling center functions to pattern the endoskeleton of gnathostome paired fins and chondrichthyan gill rays, while expression patterns of Tbx18 and HoxD family members are shared by gnathostome paired fins and chondrichthyan median fins. These findings fuel novel hypotheses of developmental genetic homology, and demonstrate how comparative studies of gnathostome development can provide insight into the evolutionary processes that underlie morphological diversity.

Evolutionary dynamics of olfactory receptor genes in chordates: interaction between environments and genomic contents

Olfaction is essential for the survival of animals. Versatile odour molecules in the environment are received by olfactory receptors (ORs), which form the largest multigene family in vertebrates. Identification of the entire repertories of OR genes using bioinformatics methods from the whole-genome sequences of diverse organisms revealed that the numbers of OR genes vary enormously, ranging from ~1,200 in rats and ~400 in humans to ~150 in zebrafish and ~15 in pufferfish. Most species have a considerable fraction of pseudogenes. Extensive phylogenetic analyses have suggested that the numbers of gene gains and losses are extremely large in the OR gene family, which is a striking example of the birth-and-death evolution. It appears that OR gene repertoires change dynamically, depending on each organism's living environment. For example, higher primates equipped with a well-developed vision system have lost a large number of OR genes. Moreover, two groups of OR genes for detecting airborne odorants greatly expanded after the time of terrestrial adaption in the tetrapod lineage, whereas fishes retain diverse repertoires of genes that were present in aquatic ancestral species. The origin of vertebrate OR genes can be traced back to the common ancestor of all chordate species, but insects, nematodes and echinoderms utilise distinctive families of chemoreceptors, suggesting that chemoreceptor genes have evolved many times independently in animal evolution.

Natural occurrence of 2',5'-linked heteronucleotides in marine sponges

2′,5′-oligoadenylate synthetases (OAS) as a component of mammalian interferon-induced antiviral enzymatic system catalyze the oligomerization of cellular ATP into 2′,5′-linked oligoadenylates (2-5A). Though vertebrate OASs have been characterized as 2′-nucleotidyl transferases under in vitro conditions, the natural occurrence of 2′,5′-oligonucleotides other than 2-5A has never been demonstrated. Here we have demonstrated that OASs from the marine sponges Thenea muricata and Chondrilla nucula are able to catalyze in vivo synthesis of 2-5A as well as the synthesis of a series 2′,5′-linked heteronucleotides which accompanied high levels of 2′,5′-diadenylates. In dephosphorylated perchloric acid extracts of the sponges, these heteronucleotides were identified as A2′p5′G, A2′ p5′U, A2′p5′C, G2′p5′A and G2′ p5′U. The natural occurrence of 2′-adenylated NAD⁺ was also detected. In vitro assays demonstrated that besides ATP, GTP was a good substrate for the sponge OAS, especially for OAS from C. nucula. Pyrimidine nucleotides UTP and CTP were also used as substrates for oligomerization, giving 2′,5′-linked homo-oligomers. These data refer to the substrate specificity of sponge OASs that is remarkably different from that of vertebrate OASs. Further studies of OASs from sponges may help to elucidate evolutionary and functional aspects of OASs as proteins of the nucleotidyltransferase family.

Evolution of developmental regulation in the vertebrate FgfD subfamily

Fibroblast growth factors (Fgfs) encode small signaling proteins that help regulate embryo patterning. Fgfs fall into seven families, including FgfD. Nonvertebrate chordates have a single FgfD gene; mammals have three (Fgf8, Fgf17, and Fgf18); and teleosts have six (fgf8a, fgf8b, fgf17, fgf18a, fgf18b, and fgf24). What are the evolutionary processes that led to the structural duplication and functional diversification of FgfD genes during vertebrate phylogeny? To study this question, we investigated conserved syntenies, patterns of gene expression, and the distribution of conserved noncoding elements (CNEs) in FgfD genes of stickleback and zebrafish, and compared them with data from cephalochordates, urochordates, and mammals. Genomic analysis suggests that Fgf8, Fgf17, Fgf18, and Fgf24 arose in two rounds of whole genome duplication at the base of the vertebrate radiation; that fgf8 and fgf18 duplications occurred at the base of the teleost radiation; and that Fgf24 is an ohnolog that was lost in the mammalian lineage. Expression analysis suggests that ancestral subfunctions partitioned between gene duplicates and points to the evolution of novel expression domains. Analysis of CNEs, at least some of which are candidate regulatory elements, suggests that ancestral CNEs partitioned between gene duplicates. These results help explain the evolutionary pathways by which the developmentally important family of FgfD molecules arose and the deduced principles that guided FgfD evolution are likely applicable to the evolution of developmental regulation in many vertebrate multigene families.

Evolution of vertebrate rod and cone phototransduction genes

Vertebrate cones and rods in several cases use separate but related components for their signal transduction (opsins, G-proteins, ion channels, etc.). Some of these proteins are also used differentially in other cell types in the retina. Because cones, rods and other retinal cell types originated in early vertebrate evolution, it is of interest to see if their specific genes arose in the extensive gene duplications that took place in the ancestor of the jawed vertebrates (gnathostomes) by two tetraploidizations (genome doublings). The ancestor of teleost fishes subsequently underwent a third tetraploidization. Our previously reported analyses showed that several gene families in the vertebrate visual phototransduction cascade received new members in the basal tetraploidizations. We here expand these data with studies of additional gene families and vertebrate species. We conclude that no less than 10 of the 13 studied phototransduction gene families received additional members in the two basal vertebrate tetraploidizations. Also the remaining three families seem to have undergone duplications during the same time period but it is unclear if this happened as a result of the tetraploidizations. The implications of the many early vertebrate gene duplications for functional specialization of specific retinal cell types, particularly cones and rods, are discussed.

Evolutionary history and functional characterization of androgen receptor genes in jawed vertebrates

Vertebrates show diverse sexual characters in sexually attractive and reproductive organs, which are regulated by steroid hormones, particularly androgens. However, the evolutionary history of androgen receptor (AR) gene remains largely unknown on the basis of phylogenic and functional analyses. To elucidate the evolutionary history and functional diversification of AR genes in vertebrates, we cloned the AR cDNAs from a shark, basal ray-finned fishes (Actinopterygii), namely bichir and sturgeon (Acipenseriformes), and teleosts including a basal teleost, arowana (Osteoglossiformes). Molecular phylogenetic analysis revealed that the gene duplication event that gave rise to two different teleost ARs (alpha and beta) likely occurred in the actinopterygian lineage leading to teleosts after the divergence of Acipenseriformes but before the split of Osteoglossiformes, which is compatible with the phylogenetic timing of teleost-specific genome duplication. Searching for AR genes in the medaka genome indicated that the teleost AR gene duplication has been associated with the duplication between chromosomes 10 and 14. Our functional analysis revealed that the shark AR activates the target gene via androgen response element by classical androgens. The teleost ARalpha showed the unique intracellular localization with a significantly higher transactivating capacity than that by teleost ARbeta. These findings indicate that the most ancient type of AR, as activated by the classical androgens as ligands, emerged before the Chondrichthyes-Osteichthyes split, and the AR gene was duplicated during the teleost-specific genome duplication event. We report here for the first time the accurate evolutionary history of AR gene and functional characterization of AR duplicates in teleost lineage.

Evolution of the 2'-5'-oligoadenylate synthetase family in eukaryotes and bacteria

The 2'-5'-oligoadenylate synthetase (OAS) belongs to a nucleotidyl transferase family that includes poly(A) polymerases and CCA-adding enzymes. In mammals and birds, the OAS functions in the interferon system but it is also present in an active form in sponges, which are devoid of the interferon system. In view of these observations, we have pursued the idea that OAS genes could be present in other metazoans and in unicellular organisms as well. We have identified a number of OAS1 genes in annelids, mollusks, a cnidarian, chordates, and unicellular eukaryotes and also found a family of proteins in bacteria that contains the five OAS-specific motifs. This indicates a specific relationship to OAS. The wide distribution of the OAS genes has made it possible to suggest how the OAS1 gene could have evolved from a common ancestor to choanoflagellates and metazoans. Furthermore, we suggest that the OASL may have evolved from an ancestor of cartilaginous fishes, and that the OAS2 and the OAS3 genes evolved from a mammalian ancestor. OAS proteins function in the interferon system in mammals. This system is only found in jawed vertebrates. We therefore suggest that the original function of OAS may differ from its function in the interferon system, and that this original function of OAS is preserved even in OAS genes that code for proteins, which do not have 2'-5'-oligoadenylate synthetase activity.

Multiple urea transporter proteins in the kidney of holocephalan elephant fish (Callorhinchus milii)

Reabsorption of filtered urea by the kidney is essential for retaining high levels of urea in marine cartilaginous fish. Our previous studies on the shark facilitative urea transporter (UT) suggest that additional UT(s) comprising the urea reabsorption system could exist in the cartilaginous fish kidney. Here, we isolated three cDNAs encoding UTs from the kidney of elephant fish, Callorhinchus milii, and termed them efUT-1, efUT-2 and efUT-3. efUT-1 is orthologous to known elasmobranch UTs, while efUT-2 and efUT-3 are novel UTs in cartilaginous fish. Two variants were found for efUT-1 and efUT-2, in which the NH(2)-terminal intracellular domain was distinct between the variants. Differences in potential phosphorylation sites were found in the variant-specific NH(2)-terminal domains. When expressed in Xenopus oocytes, all five UT transcripts including the efUT-1 and efUT-2 variants induced more than a 10-fold increase in [(14)C] urea uptake. Phloretin inhibited dose-dependently the increase of urea uptake, suggesting that the identified UTs are facilitative UTs. Molecular phylogenetic analysis revealed that efUT-1 and efUT-2 had diverged in the cartilaginous fish lineage, while efUT-3 is distinct from efUT-1 and efUT-2. The present finding of multiple UTs in elephant fish provides a key to understanding the molecular mechanisms of urea reabsorption system in the cartilaginous fish kidney.