The primary structure of piscine (Oncorhynchus mykiss) retinol-binding protein and a comparison with the three-dimensional structure of mammalian retinol-binding protein.

1. The primary structures of two variants of rainbow trout (Oncorhynchus mykiss) plasma retinol-binding protein (RBP) were determined and found to be approximately 60% identical with those of both human and Xenopus laevis RBPs. The comparable sequence similarities that we have found agree with the estimate of similar divergence times between bony fishes and mammals and between bony fishes and amphibians. The two piscine RBP variants differ by six amino acid substitutions at positions that are not crucial for the interaction with retinol, on the basis of the human RBP three-dimensional structure [Cowan, S. W., Newcomer, M. E. & Jones, T. A. (1990) Proteins Struct. Func. Genet. 8, 44-61]. 2. Models were developed for the three-dimensional structures of rainbow trout and X. laevis RBPs, based on that of human RBP. The overall three-dimensional structure appears to be very well preserved for RBPs isolated from vertebrate species for which the divergence time is 350-400 million years. At variance with an almost absolute conservation for the residues that participate in the formation of the retinol binding site in mammalian RBPs, several amino acid replacements are found for this part of the RBP molecule when the comparison is extended to piscine and amphibian RBPs. However, the only allowed amino acid replacements are either conservative or more than 0.4 nm distant from retinol. Besides the retinol binding site, a few regions at the protein surface appear to be rather conserved during phylogenetic development of vertebrates and, therefore, might be involved in molecular interactions.     

Evolutionary diversification of the vertebrate transferrin multi-gene family

In a phylogenetic analysis of vertebrate transferrins (TFs), six major clades (subfamilies) were identified: (a) S, the mammalian serotransferrins; (b) ICA, the mammalian inhibitor of carbonic anhydrase (ICA) homologs; (c) L, the mammalian lactoferrins; (d) O, the ovotransferrins of birds and reptiles; (e) M, the melanotransferrins of bony fishes, amphibians, reptiles, birds, and mammals; and (f) M-like, a newly identified TF subfamily found in bony fishes, amphibians, reptiles, and birds. A phylogenetic tree based on the joint alignment of N-lobes and C-lobes supported the hypothesis that three separate events of internal duplication occurred in vertebrate TFs: (a) in the common ancestor of the M subfamily, (b) in the common ancestor of the M-like subfamily, and (c) in the common ancestor of other vertebrate TFs. The S, ICA, and L subfamilies were found only in placental mammals, and the phylogenetic analysis supported the hypothesis that these three subfamilies arose by gene duplication after the divergence of placental mammals from marsupials. The M-like subfamily was unusual in several respects, including the presence of a uniquely high proportion of clade-specific conserved residues, including distinctive but conserved residues in the sites homologous to those functioning in carbonate binding of human serotransferrin. The M-like family also showed an unusually high proportion of cationic residues in the positively charged region corresponding to human lactoferrampin, suggesting a distinctive role of this region in the M-like subfamily, perhaps in antimicrobial defense.     

Evolution of trace amine associated receptor (TAAR) gene family in vertebrates: lineage-specific expansions and degradations of a second class of vertebrate chemosensory receptors expressed in the olfactory epithelium

The trace amine-associated receptors (TAARs) form a specific family of G protein-coupled receptors in vertebrates. TAARs were initially considered neurotransmitter receptors, but recent study showed that mouse TAARs function as chemosensory receptors in the olfactory epithelium. To clarify the evolutionary dynamics of the TAAR gene family in vertebrates, near-complete repertoires of TAAR genes and pseudogenes were identified from the genomic assemblies of 4 teleost fishes (zebrafish, fugu, stickleback, and medaka), western clawed frogs, chickens, 3 mammals (humans, mice, and opossum), and sea lampreys. Database searches revealed that fishes had many putatively functional TAAR genes (13-109 genes), whereas relatively small numbers of TAAR genes (3-22 genes) were identified in tetrapods. Phylogenetic analysis of these genes indicated that the TAAR gene family was subdivided into 5 subfamilies that diverged before the divergence of ray-finned fishes and tetrapods. In tetrapods, virtually all TAAR genes were located in 1 specific region of their genomes as a gene cluster; however, in fishes, TAAR genes were scattered throughout more than 2 genomic locations. This possibly reflects a whole-genome duplication that occurred in the common ancestor of ray-finned fishes. Expression analysis of zebrafish and stickleback TAAR genes revealed that many TAARs in these fishes were expressed in the olfactory organ, suggesting the relatively high importance of TAARs as chemosensory receptors in fishes. A possible evolutionary history of the vertebrate TAAR gene family was inferred from the phylogenetic and comparative genomic analyses.     

Molecular characterization of the mitochondrial cytochrome oxidase I gene of Oestridae species causing obligate myiasis

A 688-bp region of the mitochondrial cytochrome oxidase I gene was sequenced from larvae of 18 species of Oestridae causing obligate myiasis. Larvae belonged to the four Oestridae subfamilies (Cuterebrinae, Gasterophilinae, Hypodermatinae and Oestrinae), which are commonly found throughout the world. Analysis of both nucleotide and amino acid data was performed. Nucleotide sequences included 385 conserved sites and 303 variable sites; mean nucleotide variation between all species was 18.1% and variation within each subfamily ranged from 5.3% to 13.34%. Intraspecific pairwise divergences ranged from 0.14% to 1.59%, and interspecific variation ranged from 0.7% to 27%. Of the 229 amino acids, 76 were variable (60 of which were phylogenetically informative), with some highly conserved residues identified within each subfamily. Phylogenetic analysis showed a strong divergence among the four subfamilies, concordant with classical taxonomy based on morphological and biological features. This study provides the first molecular data set for myiasis-causing Oestridae species, providing an essential database for the molecular identification of these parasites and the assessment of phylogenetic relationships within family Oestridae.     

Identification of chemokines and a chemokine receptor in cichlid fish,shark, and lamprey

Chemokines are small, inducible, structurally related proteins that guide cells expressing the right chemokine receptors to sites of immune response. They have been identified and studied extensively in mammals, but little is known about their presence in other vertebrate groups. Here we describe seven new chemokines in bony fish and one in a cartilaginous fish, as well as one chemokine receptor in a jawless vertebrate. All eight chemokines belong to the SCYA (CC) subfamily characterized by four conserved cysteine residues of which the first two are adjacent. The chemokine receptor is of the CXCR4 type. Phylogenetic analysis does not reveal any clear evidence of orthology of fish and human chemokines. Although the divergence of the subfamilies began before the fish-tetrapod split, much of the divergence within the subfamilies took place separately in the two vertebrate groups. The existence of a chemokine receptor in the lamprey indicates that chemokines are apparently also present in the Agnatha.     

Modes of evolution in the protease and kringle domains of the plasminogen-prothrombin family

Phylogenetic analysis of protease domains of the vertebrate plasminogen-prothrombin family revealed two major subfamilies: (1) a subfamily containing macrophage-stimulating protein (MSP), hepatocyte growth factor (HGF), plasminogen, and apolipoprotein(a) (APOA); and (2) a subfamily containing prothrombin, HGF activator, and plasminogen activators. There was evidence that these two subfamilies diverged prior to the divergence of amphibians and amniotes. The phylogeny indicated a close relationship of APOA from the European hedgehog, rhesus monkey, and human with plasminogen. Phylogenetic analysis of repeated kringle domains supported the hypothesis that APOA evolved independently in hedgehog and primates through numerous duplications of different kringle domains of the ancestral plasminogen. Phylogenies of kringle domains revealed two modes of evolution: (1) a conservative mode, whereby duplication of kringle domains occurred prior to cladogenesis and the same kringle structure has been maintained in different lineages (exemplified by plasminogen and prothrombin); and (2) a concerted mode, whereby kringle domains have duplicated since cladogenesis and thus orthologous relationships do not exist between kringles of different lineages (exemplified by APOA).     

Molecular evolution and functional divergence of the cytochrome P450 3 (CYP3) Family in Actinopterygii (ray-finned fish)

Background

The cytochrome P450 (CYP) superfamily is a multifunctional hemethiolate enzyme that is widely distributed from Bacteria to Eukarya. The CYP3 family contains mainly the four subfamilies CYP3A, CYP3B, CYP3C and CYP3D in vertebrates; however, only the Actinopterygii (ray-finned fish) have all four subfamilies and detailed understanding of the evolutionary relationship of Actinopterygii CYP3 family members would be valuable.

Methods and Findings

Phylogenetic relationships were constructed to trace the evolutionary history of the Actinopterygii CYP3 family genes. Selection analysis, relative rate tests and functional divergence analysis were combined to interpret the relationship of the site-specific evolution and functional divergence in the Actinopterygii CYP3 family. The results showed that the four CYP3 subfamilies in Actinopterygii might be formed by gene duplication. The first gene duplication event was responsible for divergence of the CYP3B/C clusters from ancient CYP3 before the origin of the Actinopterygii, which corresponded to the fish-specific whole genome duplication (WGD). Tandem repeat duplication in each of the homologue clusters produced stable CYP3B, CYP3C, CYP3A and CYP3D subfamilies. Acceleration of asymmetric evolutionary rates and purifying selection together were the main force for the production of new subfamilies and functional divergence in the new subset after gene duplication, whereas positive selection was detected only in the retained CYP3A subfamily. Furthermore, nearly half of the functional divergence sites appear to be related to substrate recognition, which suggests that site-specific evolution is closely related with functional divergence in the Actinopterygii CYP3 family.

Conclusions

The split of fish-specific CYP3 subfamilies was related to the fish-specific WGD, and site-specific acceleration of asymmetric evolutionary rates and purifying selection was the main force for the origin of the new subfamilies and functional divergence in the new subset after gene duplication. Site-specific evolution in substrate recognition was related to functional divergence in the Actinopterygii CYP3 family.     

Modes of Evolution in the Protease and Kringle Domains of the Plasminogen–Prothrombin Family

Phylogenetic analysis of protease domains of the vertebrate plasminogen–prothrombin family revealed two major subfamilies: (1) a subfamily containing macrophage-stimulating protein (MSP), hepatocyte growth factor (HGF), plasminogen, and apolipoprotein(a) (APOA); and (2) a subfamily containing prothrombin, HGF activator, and plasminogen activators. There was evidence that these two subfamilies diverged prior to the divergence of amphibians and amniotes. The phylogeny indicated a close relationship of APOA from the European hedgehog, rhesus monkey, and human with plasminogen. Phylogenetic analysis of repeated kringle domains supported the hypothesis that APOA evolved independently in hedgehog and primates through numerous duplications of different kringle domains of the ancestral plasminogen. Phylogenies of kringle domains revealed two modes of evolution: (1) a conservative mode, whereby duplication of kringle domains occurred prior to cladogenesis and the same kringle structure has been maintained in different lineages (exemplified by plasminogen and prothrombin); and (2) a concerted mode, whereby kringle domains have duplicated since cladogenesis and thus orthologous relationships do not exist between kringles of different lineages (exemplified by APOA).     

Evolutionary selection pressure and family relationships among connexin genes

We suggest an extension of connexin orthology relationships across the major vertebrate lineages. We first show that the conserved domains of mammalian connexins (encoding the N-terminus, four transmembrane domains and two extracellular loops) are subjected to a considerably more strict selection pressure than the full-length sequences or the variable domains (the intracellular loop and C-terminal tail). Therefore, the conserved domains are more useful for the study of family relationships over larger evolutionary distances. The conserved domains of connexins were collected from chicken, Xenopus tropicalis, zebrafish, pufferfish, green spotted pufferfish, Ciona intestinalis and Halocynthia pyriformis (two tunicates). A total of 305 connexin sequences were included in this analysis. Phylogenetic trees were constructed, from which the orthologies and the presumed evolutionary relationships between the sequences were deduced. The tunicate connexins studied had the closest, but still distant, relationships to vertebrate connexin 36, 39.2, 43.4, 45 and 47. The main structure in the connexin family known from mammals pre-dates the divergence of bony fishes, but some additional losses and gains of connexin sequences have occurred in the evolutionary lineages of subsequent vertebrates. Thus, the connexin gene family probably originated in the early evolution of chordates, and underwent major restructuring with regard to gene and subfamily structures (including the number of genes in each subfamily) during early vertebrate evolution.     

Molecular phylogeny of pontobdelline leeches and their place in the descent of fish leeches (Hirudinea, Piscicolidae)

Phylogenetic relationships of all genera of the fish leech subfamily Pontobdellinae were investigated using mitochondrial (12S rDNA, COI, tRNA-Leu, ND1) and nuclear (28S rDNA) DNA sequences under maximum likelihood, Bayesian inference and parsimony. All methods resulted in trees that corroborated the monophyly of the family Piscicolidae, but recovered their subfamily Pontobdellinae as non-monophyletic. Based on the basal position of the giant Antarctic Megaliobdella szidati , it is hypothesized that the putative ancestor of fish leeches was a free-ranging, large bodied, muscular leech. The next branch contains parasites of cartilaginous fishes, Pontobdella muricata and Pontobdella macrothela . Two remaining genera of the subfamily (the Arctic Oxytonostoma and the Antarctic Moorebdellina ) showed weak affinities to other piscicolid taxa. The obtained phylogenetic hypothesis suggests a possible transition from an ancestral free-ranging life style and temporary parasitism, to parasitism on cartilaginous fishes, followed by parasitism on bony fishes.     

Molecular phylogeny of early vertebrates: monophyly of the agnathans as revealed by sequences of 35 genes

Extant vertebrates are divided into three major groups: hagfishes (Hyperotreti, myxinoids), lampreys (Hyperoartia, petromyzontids), and jawed vertebrates (Gnathostomata). The phylogenetic relationships among the groups and within the jawed vertebrates are controversial, for both morphological and molecular studies have rendered themselves to conflicting interpretations. Here, we use the sequences of 35 nuclear protein-encoding genes to provide definitive evidence for the monophyly of the Agnatha (jawless vertebrates, a group encompassing the hagfishes and lampreys). Our analyses also give a strong support for the separation of Chondrichthyes (cartilaginous fishes) before the divergence of Osteichthyes (bony fishes) from the other gnathostomes.     

Nonlinear relationships among evolutionary rates identify regions of functional divergence in heat-shock protein 70 genes

The neutral theory predicts that, in comparisons among related genes, the number of amino acid replacements per site in a given gene region should be a linear function of that in another region of the same gene, unless the genes have diverged functionally in one region. Therefore, nonlinearity of this relationship can be used to identify regions of possible functional divergence among members of a multigene family. This method of analysis was applied to members of the heat-shock protein 70 (HSP70) gene family, which encode highly conserved ATP- dependent chaperone proteins found in all organisms. A nonlinear relationship was found between the rate of amino acid replacement in the conserved IA domain of the ATPase portion of the molecule and that in other ATPase domains and the peptide-binding domain. These results suggest that genes in the HSP70 subfamily C (dnaK of bacteria and SSC1 of yeast) may have diverged functionally from other subfamilies in the ATPase domains, especially IIB, whereas SSB1 of yeast has diverged markedly in the peptide-binding domain. Functional divergence within these regions is consistent with what is known about functional differences between the HSP70 subfamilies in yeast.      

Functional analysis of a mammalian odorant receptor subfamily

Phylogenetic analysis groups mammalian odorant receptors into two broad classes and numerous subfamilies. These subfamilies are proposed to reflect functional organization. Testing this idea requires an assay allowing detailed functional characterization of odorant receptors. Here we show that a variety of Class I and Class II mouse odorant receptors can be functionally expressed in Xenopus laevis oocytes. Receptor constructs included the N-terminal 20 residues of human rhodopsin and were co-expressed with Galphaolf and the cystic fibrosis transmembrane regulator to allow electrophysiological measurement of receptor responses. For most mouse odorant receptors tested, these conditions were sufficient for functional expression. Co-expression of accessory proteins was required to allow functional surface expression of some mouse odorant receptors. We used this assay to examine the receptive ranges of all members of the mouse odorant receptor 42 (MOR42) subfamily. MOR42-1 responded to dicarboxylic acids, preferring a 10-12 carbon chain length. MOR42-2 responded to monocarboxylic acids (7-10 carbons). MOR42-3 responded to dicarboxylic acids (8-10 carbons) and monocarboxylic acids (10-12 carbons). Thus, the receptive range of each receptor was unique. However, overlap between the individual receptive ranges suggests that the members of this subfamily form one contiguous subfamily receptive range, suggesting that odorant receptor subfamilies do constitute functional units.     

Basal jawed vertebrate phylogeny inferred from multiple nuclear DNA-coded genes

Background  

Phylogenetic analyses of jawed vertebrates based on mitochondrial sequences often result in confusing inferences which are obviously inconsistent with generally accepted trees. In particular, in a hypothesis by Rasmussen and Arnason based on mitochondrial trees, cartilaginous fishes have a terminal position in a paraphyletic cluster of bony fishes. No previous analysis based on nuclear DNA-coded genes could significantly reject the mitochondrial trees of jawed vertebrates.     

Complete nucleotide sequence of the mitochondrial genome of a salamander, Mertensiella luschani

The complete nucleotide sequence (16,650 bp) of the mitochondrial genome of the salamander Mertensiella luschani (Caudata, Amphibia) was determined. This molecule conforms to the consensus vertebrate mitochondrial gene order. However, it is characterized by a long non-coding intervening sequence with two 124-bp repeats between the tRNA(Thr) and tRNA(Pro) genes. The new sequence data were used to reconstruct a phylogeny of jawed vertebrates. Phylogenetic analyses of all mitochondrial protein-coding genes at the amino acid level recovered a robust vertebrate tree in which lungfishes are the closest living relatives of tetrapods, salamanders and frogs are grouped together to the exclusion of caecilians (the Batrachia hypothesis) in a monophyletic amphibian clade, turtles show diapsid affinities and are placed as sister group of crocodiles+birds, and the marsupials are grouped together with monotremes and basal to placental mammals. The deduced phylogeny was used to characterize the molecular evolution of vertebrate mitochondrial proteins. Amino acid frequencies were analyzed across the main lineages of jawed vertebrates, and leucine and cysteine were found to be the most and least abundant amino acids in mitochondrial proteins, respectively. Patterns of amino acid replacements were conserved among vertebrates. Overall, cartilaginous fishes showed the least variation in amino acid frequencies and replacements. Constancy of rates of evolution among the main lineages of jawed vertebrates was rejected.     

Trans-species polymorphism of class IIMhc loci in danio fishes

 A characteristic feature of the major histocompatibility complex (Mhc) polymorphism in mammals is the existence of allelic lineages shared by related species. This trans-species polymorphism has thus far been documented only in primates, rodents, and artiodactyls. In this communication we provide evidence that it also exists in cyprinid (bony) fishes at the class II A and B loci coding for the α and β polypeptide chains of the class II α : β heterodimers. The study has focused on three species of the family Cyprinidae, subfamily Rasborinae: the zebrafish (Danio rerio), the giant danio (D. malabaricus), and the pearl danio (D. albolineatus). The polymerase chain reaction was used to amplify and then sequence intron 1 and exon 2 of the class II B loci and exon 2 of the class II A loci in these species. Phylogenetic analysis of the sequences revealed the existence of allelic lineages whose divergence predates the divergence of the three species at both the A and B loci. The lineages at the B locus in particular are separated by large genetic distances. The polymorphism is concentrated in the peptide-binding region sites and is apparently maintained by balancing selection. Sharing of this unique Mhc feature by both bony fishes and mammals suggests that the main function of the Mhc (presentation of peptides to T lymphocytes) has not changed during the last 400 million years of its evolution. Received: 6 December 1995 / Revised: 6 February 1996     

Evolution of duplicated growth hormone genes in autotetraploid salmonid fishes.

A defining character of the piscine family Salmonidae is autotetraploidy resulting from a genome-doubling event some 25-100 million years ago. Initially, duplicated genes may have undergone concerted evolution and tetrasomic inheritance. Homeologous chromosomes eventually diverged and the resulting reduction in recombination and gene conversion between paralogous genes allowed the re-establishment of disomic inheritance. Among extant salmonine fishes (e.g. salmon, trout, char) the growth hormone (GH) gene is generally represented by two functional paralogs, GH1 and GH2. Sequence analyses of salmonid GH genes from species of subfamilies Coregoninae (whitefish, ciscos) and Salmoninae were used to examine the evolutionary history of the duplicated GH genes. Two divergent GH gene paralogs were also identified in Coregoninae, but they were not assignable to the GH1 and GH2 categories. The average sequence divergence between the coregonine GH genes was more than twofold lower than the corresponding divergence between the salmonine GH1 and GH2. Phylogenetic analysis of the coregonine GH paralogs did not resolve their relationship to the salmonine paralogs. These findings suggest that disomic inheritance of two GH genes was established by different mechanisms in these two subfamilies.     

Molecular evidence for the monophyly of East Asian groups of Cyprinidae (Teleostei: Cypriniformes) derived from the nuclear recombination activating gene 2 sequences

The family Cyprinidae is one of the largest families of fishes in the world and a well-known component of the East Asian freshwater fish fauna. However, the phylogenetic relationships among cyprinids are still poorly understood despite much effort paid on the cyprinid molecular phylogenetics. Original nucleotide sequence data of the nuclear recombination activating gene 2 were collected from 109 cyprinid species and four non-cyprinid cypriniform outgroup taxa and used to infer the cyprinid phylogenetic relationships and to estimate node divergence times. Phylogenetic reconstructions using maximum parsimony, maximum likelihood, and Bayesian analysis retrieved the same clades, only branching order within these clades varied slightly between trees. Although the morphological diversity is remarkable, the endemic cyprinid taxa in East Asia emerged as a monophyletic clade referred to as Xenocypridini. The monophyly for the subfamilies including Cyprininae and Leuciscinae, as well as the tribes including Labeonini, Gobionini, Acheilognathini, and Leuciscini, was also well resolved with high nodal support. Analysis of the RAG2 gene supported the following cyprinid molecular phylogeny: the Danioninae is the most basal subfamily within the family Cyprinidae and the Cyprininae is the sister group of the Leuciscinae. The divergence times were estimated for the nodes corresponding to the principal clades within the Cyprinidae. The family Cyprinidae appears to have originated in the mid-Eocene in Asia, with the cladogenic event of the key basal group Danioninae occurring in the early Oligocene (about 31-30 MYA), and the origins of the two subfamilies, Cyprininae and Leuciscinae, occurring in the mid-Oligocene (around 26 MYA).