Molecular Characterization and Evolution of the Amylase Multigene Family of Drosophila ananassae

Drosophila ananassae is known to produce numerous alpha-amylase variants. We have cloned seven different Amy genes in an African strain homozygous for the AMY1,2,3,4 electrophoretic pattern. These genes are organized as two main clusters: the first one contains three intronless copies on the 2L chromosome arm, two of which are tandemly arranged. The other cluster, on the 3L arm, contains two intron-bearing copies. The amylase variants AMY1 and AMY2 have been assigned to the intronless cluster, and AMY3 and AMY4 to the second one. The divergence of coding sequences between clusters is moderate (6.1% in amino acids), but the flanking regions are very different, which could explain their differential regulation. Within each cluster, coding and noncoding regions are conserved. Two very divergent genes were also cloned, both on chromosome 3L, but very distant from each other and from the other genes. One is the Amyrel homologous (41% divergent), the second one, Amyc1 (21.6% divergent) is unknown outside the D. ananassae subgroup. These two genes have unknown functions. Received: 30 May 2000 / Accepted: 17 July 2000     

Phylogenetic Analysis of Vertebrate Lactate Dehydrogenase (LDH) Multigene Families

In this paper we analyzed 49 lactate dehydrogenase (LDH) sequences, mostly from vertebrates. The amino acid sequence differences were found to be larger for a human–killifish pair than a human–lamprey pair. This indicates that some protein sequence convergence may occur and reduce the sequence differences in distantly related species. We also examined transitions and transversions separately for several species pairs and found that the transitions tend to be saturated in the distantly related species pair, while transversions are increasing. We conclude that transversions maintain a conservative rate through the evolutionary time. Kimura's two-parameter model for multiple-hit correction on transversions only was used to derive a distance measure and then construct a neighbor-joining (NJ) tree. Three findings were revealed from the NJ tree: (i) the branching order of the tree is consistent with the common branch pattern of major vertebrates; (ii) Ldh-A and Ldh-B genes were duplicated near the origin of vertebrates; and (iii) Ldh-C and Ldh-A in mammals were produced by an independent gene duplication in early mammalian history. Furthermore, a relative rate test showed that mammalian Ldh-C evolved more rapidly than mammalian Ldh-A. Under a two-rate model, this duplication event was calibrated to be approximately 247 million years ago (mya), dating back to the Triassic period. Other gene duplication events were also discovered in Xenopus, the first duplication occurring approximately 60–70 mya in both Ldh-A and Ldh-B, followed by another recent gene duplication event, approximately 20 mya, in Ldh-B. Received: 5 October 2001 / Accepted: 24 October 2001     

Gene Conversions May Obscure Actin Gene Family Relationships

Phylogenetic hypotheses of muscle actin evolution are significantly different when a sea urchin is used as a representative echinoderm than when a sea star is used. While sea urchin muscle actins support an echinoderm–chordate sister relationship, sea star sequences suggest that echinoderm muscle actins are convergent with chordate muscle actins. Our results suggest that gene conversion in the sea star muscle actin may be responsible for these discordant results. Received: 19 July 1999 / Accepted: 1 October 1999     

Evolution of Chordate Actin Genes: Evidence from Genomic Organization and Amino Acid Sequences

The origin and evolutionary relationship of actin isoforms was investigated in chordates by isolating and characterizing two new ascidian cytoplasmic and muscle actin genes. The exon–intron organization and sequences of these genes were compared with those of other invertebrate and vertebrate actin genes. The gene HrCA1 encodes a cytoplasmic (nonmuscle)-type actin, whereas the MocuMA2 gene encodes an adult muscle-type actin. Our analysis of these genes showed that intron positions are conserved among the deuterostome actin genes. This suggests that actin gene families evolved from a single actin gene in the ancestral deuterostome. Sequence comparisons and molecular phylogenetic analyses also suggested a close relationship between the ascidian and vertebrate actin isoforms. It was also found that there are two distinct lineages of muscle actin isoforms in ascidians: the larval muscle and adult body-wall isoforms. The four muscle isoforms in vertebrates show a closer relationship to each other than to the ascidian muscle isoforms. Similarly, the two cytoplasmic isoforms in vertebrates show a closer relationship to each other than to the ascidian and echinoderm cytoplasmic isoforms. In contrast, the two types of ascidian muscle actin diverge from each other. The close relationship between the ascidian larval muscle actin and the vertebrate muscle isoforms was supported by both neighbor-joining and maximum parsimony analyses. These results suggest that the chordate ancestor had at least two muscle actin isoforms and that the vertebrate actin isoforms evolved after the separation of the vertebrates and urochordates. Received: 20 June 1996 / Accepted: 16 October 1996     

Phylogenetic Characterization of the MIP Family of Transmembrane Channel Proteins

The ubiquitous major intrinsic protein (MIP) family includes several transmembrane channel proteins known to exhibit specificity for water and/or neutral solutes. We have identified 84 fully or partially sequenced members of this family, have multiply aligned over 50 representative, divergent, fully sequenced members, have used the resultant multiple alignment to derive current MIP family-specific signature sequences, and have constructed a phylogenetic tree. The tree reveals novel features relevant to the evolutionary history of this protein family. These features plus an evaluation of functional studies lead to the postulates: (i) that all current MIP family proteins derived from two divergent bacterial paralogues, one a glycerol facilitator, the other an aquaporin, and (ii) that most or all current members of the family have retained these or closely related physiological functions. Received: 19 April 1996/Revised: 3 June 1996     

Sequence Analyses and Phylogenetic Characterization of the ZIP Family of Metal Ion Transport Proteins

Several novel but similar heavy metal ion transporters, Zrt1, Zrt2, Zip1-4 and Irt1, have recently been characterized. Zrt1, Zrt2 and Zip1-4 are probably zinc transporters in Saccharomyces cerevisiae and Arabidopsis thaliana whereas Irt1 appears to play a role in iron uptake in A. thaliana. The family of proteins including these functionally characterized transporters has been designated the Zrt- and Irt-related protein (ZIP) family. In this report, ZIP family proteins in the current databases were identified and multiply aligned, and a phylogenetic tree for the family was constructed. A family specific signature sequence was derived, and the available sequences were analyzed for residues of potential functional significance. A fully conserved intramembranous histidyl residue, present within a putative amphipathic, α-helical, transmembrane spanning segment, was identified which may serve as a part of an intrachannel heavy metal ion binding site. The occurrence of a proposed extramembranal metal binding motif (H X H X H) was examined in order to evaluate its potential functional significance for various members of the family. The computational analyses reported in this topical review should serve as a guide to future researchers interested in the structure-function relationships of ZIP family proteins. Received: 31 March 1997/Revised: 14 May 1998     

Molecular Phylogenetic Characterization of Streptomyces Protease Inhibitor Family

We previously found that proteinaceous protease inhibitors homologous to Streptomyces subtilisin inhibitor (SSI) are widely produced by various Streptomyces species, and we designated them ``SSI-like proteins' (Taguchi S, Kikuchi H, Suzuki M, Kojima S, Terabe M, Miura K, Nakase T, Momose H [1993] Appl Environ Microbiol 59:4338–4341). In this study, SSI-like proteins from five strains of the genus Streptoverticillium were purified and sequenced, and molecular phylogenetic trees were constructed on the basis of the determined amino acid sequences together with those determined previously for Streptomyces species. The phylogenetic trees showed that SSI-like proteins from Streptoverticillium species are phylogenetically included in Streptomyces SSI-like proteins but form a monophyletic group as a distinct lineage within the Streptomyces proteins. This provides an alternative phylogenetic framework to the previous one based on partial small ribosomal RNA sequences, and it may indicate that the phylogenetic affiliation of the genus Streptoverticillium should be revised. The phylogenetic trees also suggested that SSI-like proteins possessing arginine or methionine at the P1 site, the major reactive center site toward target proteases, arose multiple times on independent lineages from ancestral proteins possessing lysine at the P1 site. Most of the codon changes at the P1 site inferred to have occurred during the evolution of SSI-like proteins are consistent with those inferred from the extremely high G + C content of Streptomyces genomes. The inferred minimum number of amino acid replacements at the P1 site was nearly equal to the average number for all the variable sites. It thus appears that positive Darwinian selection, which has been postulated to account for accelerated rates of amino acid replacement at the major reaction center site of mammalian protease inhibitors, may not have dictated the evolution of the bacterial SSI-like proteins. Received: 23 August 1996 / Accepted: 20 November 1996     

A Phylogenetic Framework for the Aquaporin Family in Eukaryotes

A comprehensive evolutionary analysis of aquaporins, a family of intrinsic membrane proteins that function as water channels, was conducted to establish groups of homology (i.e., to identify orthologues and paralogues) within the family and to gain insights into the functional constraints acting on the structure of the aquaporin molecule structure. Aquaporins are present in all living organisms, and therefore, they provide an excellent opportunity to further our understanding of the broader biological significance of molecular evolution by gene duplication followed by functional and structural specialization. Based on the resulting phylogeny, the 153 channel proteins analyzed were classified into six major paralogous groups: (1) GLPs, or glycerol-transporting channel proteins, which include mammalian AQP3, AQP7, and AQP9, several nematode paralogues, a yeast paralogue, and Escherichia coli GLP; (2) AQPs, or aquaporins, which include metazoan AQP0, AQP1, AQP2, AQP4, AQP5, and AQP6; (3) PIPs, or plasma membrane intrinsic proteins of plants, which include PIP1 and PIP2; (4) TIPs, or tonoplast intrinsic proteins of plants, which include alphaTIP, gammaTIP, and deltaTIP; (5) NODs, or nodulins of plants; and (6) AQP8s, or metazoan aquaporin 8 proteins. Of these groups, AQPs, PIPs, and TIPs cluster together. According to the results, the capacity to transport glycerol shown by several members of the family was acquired only early in the history of the family. The new phylogeny reveals that several water channel proteins are misclassified and require reassignment, whereas several previously undetermined ones can now be classified with confidence. The deduced phylogenetic framework was used to characterize the molecular features of water channel proteins. Three motifs are common to all family members: AEF (Ala-Glu-Phe), which is located in the N-terminal domain; and two NPA (Asp-Pro-Ala) boxes, which are located in the center and C-terminal domains, respectively. Other residues are found to be conserved within the major groups but not among them. Overall, the PIP subfamily showed the least variation. In general, no radical amino acid replacements affecting tertiary structure were identified, with the exception of Ala-->Ser in the TIP subfamily. Constancy of rates of evolution was demonstrated within the different paralogues but rejected among several of them (GLP and NOD).     

A Phylogenetic Study of the Origin of the Domestic Pig Estimated from the Near-Complete mtDNA Genome

The near-complete pig mtDNA genome sequence (15,997 bp) was determined from two domestic pigs (one Chinese Meishan and one Swedish Landrace) and two European wild boars. The sequences were analyzed together with a previously published sequence representing a Swedish domestic pig. The sequences formed three distinct clades, denoted A, E1, and E2, with considerable sequence divergence between them (0.8–1.2%). The results confirm our previous study (based on the sequence of the cytochrome B gene and the control region only) and provide compelling evidence that domestication of pigs must have occurred from both an Asian and a European subspecies of the wild boar. We estimated the time since the divergence of clade A (found in Chinese Meishan pigs) and E1 (found in European domestic pigs) at about 900,000 years before present, long before domestication about 9000 years ago. The pattern of nucleotide substitutions among the sequences was in good agreement with previous interspecific comparisons of mammalian mtDNA; the lowest substitution rates were observed at nonsynonymous sites in protein-coding genes, in the tRNA and rRNA genes, while the highest rates were observed at synonymous sites and in the control region. The presence of Asian clade A in some major European breeds (Large White and Landrace) most likely reflects the documented introgression of Asian germplasm into European stocks during the 18th and 19th centuries. The coexistence of such divergent mtDNA haplotypes for 100+ generations is expected to lead to the presence of recombinant haplotypes if paternal transmission and recombination occur at a low frequency. We found no evidence of such recombination events in the limited sample studied so far. Received: 19 April 2000; Accepted: 15 November 2000     

Phylogenetic Analysis of a Retroposon Family in African Great Apes

The SINE-R retroposon family has been identified by its relationship with the long terminal repeats (LTRs) of human endogenous retrovirus class K (HERV-K) as a mobile element that has evolved recently in the human genome. Here we examined the recent evolutionary history of this class of elements by a PCR approach to genomic DNA from the African great apes and by phylogenetic analysis including comparison with the HERV K10 parent sequence. With primers derived from a cDNA sequence from human brain, we identified 27 sequences from the chimpanzee and 16 from the gorilla. Phylogenetic comparisons with previously recognized sequences from the human and from the orangutan and gibbon revealed wide overlap of elements across species, suggesting multiple origins in the course of hominoid evolution. Two human elements SINE-R.C2 and HS307 were the furthest removed from the HERV-K10 sequence but these two elements were closely related to three elements from the chimpanzee and four elements from the gorilla. This group of elements (our clusters 14 and 15) appears to have transposed late in hominoid evolution. One element (Ch-M16) showed 99.1% sequence identity with the SINE-R.C2 element, which is human-specific. Thus the SINE-R family appears to have continued to be active in transposition throughout the course of primate evolution. Received: 12 March 1999 / Accepted: 25 May 1999     

Phylogenetic Analysis of Components of the Eukaryotic Vesicle Transport System Reveals a Common Origin of Adaptor Protein Complexes 1, 2, and 3 and the F Subcomplex of the Coatomer COPI

Eukaryotic vesicular transport requires the recognition of membranes through specific protein complexes. The heterotetrameric adaptor protein complexes 1, 2, and 3 (AP1/2/3) are composed of two large, one small, and one medium adaptin subunit. We isolated and characterized the cDNA for Arabidopsisγ-adaptin and performed a phylogenetic analysis of all adaptin subunits (proteins) in the context of all known homologous proteins. This analysis revealed (i) that the large subunits of AP1/2/3 are homologous and (ii) two subunits of the heptameric coatomer I (COPI) complex belong to this gene family. In addition, all small subunits and the aminoterminal domain of the medium subunits of the heterotetramers are homologous to each other; this also holds for two corresponding subunits of the COPI complex. AP1/2/3 and a substructure (heterotetrameric, F-COPI subcomplex) of the heptameric COPI had a common ancestral complex (called pre-F-COPI). Since all large and all small/medium subunits share sequence similarity, the ancestor of this complex is inferred to have been a heterodimer composed of one large and one small subunit. The situation encountered today is the result of successive rounds of coordinated gene duplications of both the large and the small/medium subunits, with F-COPI being the first that separated from the ancestral pre-F-COPI. Received: 1 October 1998 / Accepted: 4 January 1999     

Phylogenetic Position of the Hexactinellida Within the Phylum Porifera Based on the Amino Acid Sequence of the Protein Kinase C from Rhabdocalyptus dawsoni

Recent analyses of genes encoding proteins typical for multicellularity, especially adhesion molecules and receptors, favor the conclusion that all metazoan phyla, including the phylum Porifera (sponges), are of monophyletic origin. However, none of these data includes cDNA encoding a protein from the sponge class Hexactinellida. We have now isolated and characterized the cDNA encoding a protein kinase C, belonging to the C subfamily (cPKC), from the hexactinellid sponge Rhabdocalyptus dawsoni. The two conserved regions, the regulatory part with the pseudosubstrate site, the two zinc fingers, and the C2 domain, as well as the catalytic domain were used for phylogenetic analyses. Sequence alignment and construction of a phylogenetic tree from the catalytic domains revealed that the yeast Saccharomyces cerevisiae and the protozoan Trypanosoma brucei are at the base of the tree, while the hexactinellid R. dawsoni branches off first among the metazoan sequences; the other two classes of the Porifera, the Calcarea (the sequence from Sycon raphanus was used) and the Demospongiae (sequences from Geodia cydonium and Suberites domuncula were used), branch off later. The statistically robust tree also shows that the two cPKC sequences from the higher invertebrates Drosophila melanogaster and Lytechinus pictus are most closely related to the calcareous sponge. This finding was also confirmed by comparing the regulatory part of the kinase gene. We suggest, that (i) within the phylum Porifera, the class Hexactinellida diverged first from a common ancestor to the Calcarea and the Demospongiae, which both appeared later, and (ii) the higher invertebrates are more closely related to the calcareous sponges. Received: 6 August 1997 / Accepted: 24 October 1997     

Classification and Phylogenetic Analysis of the cAMP-Dependent Protein Kinase Regulatory Subunit Family

The members of the PKA regulatory subunit family (PKA-R family) were analyzed by multiple sequence alignment and clustering based on phylogenetic tree construction. According to the phylogenetic trees generated from multiple sequence alignment of the complete sequences, the PKA-R family was divided into four subfamilies (types I to IV). Members of each subfamily were exclusively from animals (types I and II), fungi (type III), and alveolates (type IV). Application of the same methodology to the cAMP-binding domains, and subsequently to the region delimited by β-strands 6 and 7 of the crystal structures of bovine RIα and rat RIIβ (the phosphate-binding cassette; PBC), proved that this highly conserved region was enough to classify unequivocally the members of the PKA-R family. A single signature sequence, F–G–E–[LIV]–A–L–[LIMV]–x(3)–[PV]–R–[ANQV]–A, corresponding to the PBC was identified which is characteristic of the PKA-R family and is sufficient to distinguish it from other members of the cyclic nucleotide-binding protein superfamily. Specific determinants for the A and B domains of each R-subunit type were also identified. Conserved residues defining the signature motif are important for interaction with cAMP or for positioning the residues that directly interact with cAMP. Conversely, residues that define subfamilies or domain types are not conserved and are mostly located on the loop that connects α-helix B′ and β strand 7. Received: 2 November 2000/Accepted: 14 June 2001     

Isolation and Characterization of a Carbonic Anhydrase Homologue from the Zebrafish (Danio rerio)

We have isolated a 29,000-Da carbonic anhydrase (CA) protein from the zebrafish, Danio rerio, sequenced two peptide fragments, and tentatively identified it as a high-activity CA by inhibition kinetics. We have also characterized a 1,537-bp message whose deduced sequence of 260 amino acids matches that of the isolated protein. This CA is clearly an α-CA based on the similarity of its sequence to that of other members of the α-CA gene family. A phylogenetic analysis suggested CAH-Z diverged after the branching of the CA-V and CA-VII genes and prior to the duplications that generated the CA-I, CA-II, and CA-III genes of amniotes. This marks the first characterization of the mRNA and its protein product from the CA gene of a teleost. Received: 31 March 1996 / Accepted: 8 September 1996     

Positive Darwinian Selection Promotes Heterogeneity Among Members of the Antifreeze Protein Multigene Family

A variety of organisms have independently evolved proteins exhibiting antifreeze activity that allows survival at subfreezing temperatures. The antifreeze proteins (AFPs) bind ice nuclei and depress the freezing point by a noncolligative absorption–inhibition mechanism. Many organisms have a heterogeneous suite of AFPs with variation in primary sequence between paralogous loci. Here, we demonstrate that the diversification of the AFP paralogues is promoted by positive Darwinian selection in two independently evolved AFPs from fish and beetle. First, we demonstrate an elevated rate of nonsynonymous substitutions compared to synonymous substitutions in the mature protein coding region. Second, we perform phylogeny-based tests of selection to demonstrate a subset of codons is subjected to positive selection. When mapped onto the three-dimensional structure of the fish antifreeze type III antifreeze structure, these codons correspond to amino acid positions that surround but do not interrupt the putative ice-binding surface. The selective agent may be related to efficient binding to diverse ice surfaces or some other aspect of AFP function. Received: 27 February 2001 / Accepted: 12 September 2001     

DNA Polymerase C of the Thermophilic Bacterium Thermus aquaticus: Classification and Phylogenetic Analysis of the Family C DNA Polymerases

Bacterial family C DNA polymerases (DNA pol IIIs), the major chromosomal replicative enzymes, have been provisionally classified based on primary sequences and domain structures into three classes: class I (Escherichia coli DNA pol C-type), class II (Bacillus subtilis DNA pol C-type), and class III (cyanobacterial DNA pol C-type), respectively. We have sequenced the structural gene encoding the DNA pol C catalytic subunit of the thermophilic bacterium Thermus aquaticus. This gene, designated the Taq DNA pol C gene, contains a 3660-bp open reading frame which specifies a polypeptide of molecular weight of 137,388 daltons. Comparative sequence analyses revealed that Taq DNA pol C is a class I family C DNA polymerase. The Taq DNA pol C is most closely related to the Deinococcus radiodurans DNA pol C. Although a phylogenetic tree based on the class I family C DNA pols is still in the provisional stage, some important conclusion can be drawn. First, the high-G+C and the low-G+C Gram-positive bacteria are not monophyletic. Second, the low-G+C Gram-positive bacteria contain multigenes of family C DNA pols (classes I and II). Third, the cyanobacterial family C DNA pol, classified as class III because it is encoded by a split gene, forms a group with the high-G+C Gram-positive bacteria. Received: 7 October 1998 / Accepted: 12 January 1999     

Origin and Evolution of Paralogous rRNA Gene Clusters Within the Flatworm Family Dugesiidae (Platyhelminthes, Tricladida)

Analysis of the 18S rDNA sequences of five species of the family Dugesiidae (phylum Platyhelminthes, suborder Tricladida, infraorder Paludicola) and eight species belonging to families Dendrocoelidae and Planaridae and to the infraorder Maricola showed that members of the family Dugesiidae have two types of 18S rDNA genes, while the rest of the species have only one. The duplication event also affected the ITS-1, 5.8S, ITS-2 region and probably the 28S gene. The mean divergence value between the type I and the type II sequences is 9% and type II 18S rDNA genes are evolving 2.3 times more rapidly than type I. The evolutionary rates of type I and type II genes were calibrated from biogeographical data, and an approximate date for the duplication event of 80–120 million years ago was calculated. The type II gene was shown, by RT-PCR, to be transcribed in adult individuals of Schmidtea polychroa, though at very low levels. This result, together with the fact that most of the functionally important positions for small-subunit rRNA in prokaryotes have been conserved, indicates that the type II gene is probably functional. Received: 24 March 1998 / Accepted: 17 March 1999     

A Survey of the Molecular Evolutionary Dynamics of Twenty-Five Multigene Families from Four Grass Taxa

We surveyed the molecular evolutionary characteristics of 25 plant gene families, with the goal of better understanding general processes in plant gene family evolution. The survey was based on 247 GenBank sequences representing four grass species (maize, rice, wheat, and barley). For each gene family, orthology and paralogy relationships were uncertain. Recognizing this uncertainty, we characterized the molecular evolution of each gene family in four ways. First, we calculated the ratio of nonsynonymous to synonymous substitutions (d _N/d _S) both on branches of gene phylogenies and across codons. Our results indicated that the d _N/d _S ratio was statistically heterogeneous across branches in 17 of 25 (68%) gene families. The vast majority of d _N/d _S estimates were <<1.0, suggestive of selective constraint on amino acid replacements, and no estimates were >1.0, either across phylogenetic lineages or across codons. Second, we tested separately for nonsynonymous and synonymous molecular clocks. Sixty-eight percent of gene families rejected a nonsynonymous molecular clock, and 52% of gene families rejected a synonymous molecular clock. Thus, most gene families in this study deviated from clock-like evolution at either synonymous or nonsynonymous sites. Third, we calculated the effective number of codons and the proportion of G+C synonymous sites for each sequence in each gene family. One or both quantities vary significantly within 18 of 25 gene families. Finally, we tested for gene conversion, and only six gene families provided evidence of gene conversion events. Altogether, evolution for these 25 gene families is marked by selective constraint that varies among gene family members, a lack of molecular clock at both synonymous and nonsynonymous sites, and substantial variation in codon usage. Received: 25 May 2000 / Accepted: 16 October 2000     

Elucidating the Population Histories and Transmission Dynamics of Papillomaviruses Using Phylogenetic Trees

Using gene genealogies constructed from gene sequence data, we show that both the mucosal and cutaneous papillomaviruses (PV)—supergroups A and B—appear to have been transmitted through susceptible populations faster than exponentially. The data and methods involved (1) examining the PV database for phylogenetic signal in an L1 open reading frame (ORF) fragment and an E1 ORF segment, (2) demonstrating that the same two fragments have evolved in a way consistent with a molecular clock, and (3) applying methods of phylogenetic tree analysis that test different scenarios for the dynamics of viral transmission within populations. The results indicate increases in PV populations of both supergroups A and B in the recent past. This form of the increases, which fit a null model of population growth with an exponent increasing in time, is compatible with the fact that human populations have grown at a faster than exponential rate, thus increasing the numbers of susceptible hosts for HPVs. There are, however, indications that the population of supergroup A has now stopped increasing in size. Received: 4 June 1996 / Accepted: 12 August 1996     

Molecular Evolution of the Aldo-keto Reductase Gene Superfamily

The aldo-keto reductase enzymes comprise a functionally diverse gene family which catalyze the NADPH-dependant reduction of a variety of carbonyl compounds. The protein sequences of 45 members of this family were aligned and phylogenetic trees were deduced from this alignment using the neighbor-joining and Fitch algorithms. The branching order of these trees indicates that the vertebrate enzymes cluster in three groups, which have a monophyletic origin distinct from the bacterial, plant, and invertebrate enzymes. A high level of conservation was observed between the vertebrate hydroxysteroid dehydrogenase enzymes, prostaglandin F synthase, and ρ-crystallin of Xenopus laevis. We infer from the phylogenetic analysis that prostaglandin F synthase may represent a recent recruit to the eicosanoid biosynthetic pathway from the hydroxysteroid dehydrogenase pathway and furthermore that, in the context of gene recruitment, Xenopus laevisρ-crystallin may represent a shared gene. Received: 26 August 1996 / Accepted: 5 June 1997