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     

The Role of MIP in Lens Fiber Cell Membrane Transport

MIP has been hypothesized to be a gap junction protein, a membrane ion channel, a membrane water channel and a facilitator of glycerol transport and metabolism. These possible roles have been indirectly suggested by the localization of MIP in lens gap junctional plaques and the properties of MIP when reconstituted into artificial membranes or exogenously expressed in oocytes. We have examined lens fiber cells to see if these functions are present and whether they are affected by a mutation of MIP found in Cat ^Fr mouse lens. Of these five hypothesized functions, only one, the role of water channel, appears to be true of fiber cells in situ. Based on the rate of volume change of vesicles placed in a hypertonic solution, fiber cell membrane lipids have a low water permeability (p _H2O ) on the order of 1 μm/sec whereas normal fiber cell membrane p _H2O was 17 μm/sec frog, 32 μm/sec rabbit and 43 μm/sec mouse. Cat ^Fr mouse lens fiber cell p _H2O was reduced by 13 μm/sec for heterozygous and 30 μm/sec for homozygous mutants when compared to wild type. Lastly, when expressed in oocytes, the p _H2O conferred by MIP is not sensitive to Hg²⁺ whereas that of CHIP28 (AQP1) is blocked by Hg²⁺. The fiber cell membrane p _H2O was also not sensitive to Hg²⁺ whereas lens epithelial cell p _H2O (136 μm/sec in rabbit) was blocked by Hg²⁺. With regard to the other hypothesized roles, fiber cell membrane or lipid vesicles had a glycerol permeability on the order of 1 nm/sec, an order of magnitude less than that conferred by MIP when expressed in oocytes. Impedance studies were employed to determine gap junctional coupling and fiber cell membrane conductance in wild-type and heterozygous Cat ^Fr mouse lenses. There was no detectable difference in either coupling or conductance between the wild-type and the mutant lenses. Received: 17 February 1999/Revised: 16 April 1999     

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).     

Families of Proteins Forming Transmembrane Channels

Channel-forming proteins/peptides fall into over 100 currently recognized families, most of which are restricted to prokaryotes or eukaryotes, but a few of which are ubiquitous. These proteins fall into three major currently recognized classes: (i) α-helix-type channels present in bacterial, archaeal and eukaryotic cytoplasmic and organellar membranes, (ii) β-barrel-type porins present in the outer membranes of Gram-negative bacterial cells, mitochondria and chloroplasts, and (iii) protein/peptide toxins targeted to the cytoplasmic membranes of cells other than those that synthesize the toxins. High-resolution 3-dimensional structural data are available for representative proteins/peptides of all three of these channel-forming types. Each type exhibits distinctive features that distinguish them from the other channel protein types and from carriers. Structural, functional, and evolutionary aspects of transmembrane channel-formers are discussed. Received: 10 September 1999/Revised: 11 February 2000     

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     

Identification of Major Phylogenetic Branches of Inhibitory Ligand-Gated Channel Receptors

The gene superfamily of ligand-gated ion channel (LGIC) receptors is composed of members of excitatory LGIC receptors (ELGIC) and inhibitory LGIC receptors (ILGIC), all using amino acids as ligands. The ILGICs, including GABA_A, Gly, and GluCl receptors, conduct Cl⁻ when the ligand is bound. To evaluate the phylogenetic relationships among ILGIC members, 90 protein sequences were analyzed by both maximum-parsimony and distance matrix-based methods. The strength of the resulting phylogenetic trees was evaluated by means of bootstrap. Four major phylogenetic branches are recognized. Branch I, called BZ, for the majority of the members are known to be related to benzodiazepine binding, is subdivided into IA, composed of all GABA_A receptor α subunits, and IB, composed of the γ and ε subunits, which are shown to be tightly linked. Branch II, named NB for non–benzodiazepine binding, and consisting of GABA_A receptor β, δ, π, and ρ subunits, is further subdivided into IIA, containing β subunits; IIB, containing δ, and π subunits; and IIC, containing ρ subunits. Branch IIIA, composed of vertebrate Gly receptors, is loosely clustered with Branch IIIB, composed of invertebrate GluCl receptors, to form Branch III, which is designated NA for being non–GABA responsive. Branch IV is called UD for being undefined in specificity. The existence of primitive forms of GABA_A receptor non-β subunits in invertebrates is first suggested by the present analysis, and the identities of sequences p25123 from Drosophila melanogaster, s34469 from Lymnaea stagnalis, and u14635 and p41849 from C. aenorhabditis elegans are determined to be different from their previously given annotations. The proposed branching classification of ILGICs provides a phylogenetic map, based on protein sequences, for tracing the evolutionary pathways of ILGIC receptor subunits and determining the identities of newly discovered subunits on the basis of their protein sequences. Received: 15 April 1997 / Accepted: 11 March 1998     

Isolation and Characterization of Five Actin cDNAs from the Cestode Diphyllobothrium dendriticum: A Phylogenetic Study of the Multigene Family

Five cDNAs (pDidact2–pDidact6), representing different actin genes, were isolated from a Diphyllobothrium dendriticum cDNA library, and the DNA as well as the putative amino acid sequences were determined. The corresponding Didact2 and Didact4 genes code for peptides 376 amino acids long, with molecular weights 41,772 and 41,744 Da, respectively, while the deduced Didact3 protein is 377 amino acids long and weighs 41,912 Da. The pDidact5 and -6 cDNAs lack nucleotides corresponding to three to six amino acids at the amino-terminus. Two of the five cDNAs contain the conventional AATAAA as the putative polyadenylation signal, one has the common variant ATTAAA, whereas the hexanucleotide AATAGA is found 15 and 18 nucleotides, respectively, upstream of the poly(A) site in two of the cDNAs. Phylogenetic studies including 102 actin protein sequences revealed that there are at least four different types of cestode actins. In this study three of these types were found to be expressed in the adult D. dendriticum tapeworm. Structurally the cestode actin groupings differ from each other to an extent seen only among the metazoan actins between the vertebrate muscle and cytoplasmic isoforms. In the phylogenetic trees constructed, cestode actins were seen to map to two different regions, one on the border of the metazoan actins and the other within this group. It is, however, difficult to say whether the cestode actins branched off early in the metazoan evolution or if this position in the phylogenetic tree only reflects upon differences in evolutionary rate. Received: 19 June 1996 / Accepted: 20 August 1996     

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,Structural and Functional Characteristics of the Na-K-Cl Cotransporter Family

Summary Bumetanide-sensitive Na-K-Cl cotransporters and thiazide-sensitive Na-Cl cotransporters comprise a family of integral membrane transport proteins, the Na-K-Cl cotransporter (NKCC) family. Each of the members of this family is over 1,000 amino acids in length. We have multiply aligned the ten currently sequenced members of this family from human, rabbit, rodent, shark, flounder, moth, worm and yeast sources. Phylogenetic analyses suggest the presence of at least six isoforms of these full length proteins in eukaryotes. Average hydropathy and average similarity plots have been derived revealing that each of these proteins possesses a central, well conserved, hydrophobic domain of almost invariant length, possibly consisting of twelve transmembrane α-helical spanners, an N-terminal, poorly conserved, hydrophilic domain of variable length, and a C-terminal, moderately conserved, hydrophilic domain of moderately constant length. A functionally uncharacterized homologue of this family occurs in the cyanobacterium Synechococcus sp. Limited sequence similarity of these proteins with members of a family of basic amino acid transporters suggests that the NKCC family may be distantly related to the previously characterized, ubiquitous, amino acid-polyamine-choline (APC) family of facilitators. These observations suggest that the NKCC family is an old family that has its roots in the prokaryotic kingdom. Received: 27 July 1995/Revised: 8 November 1995     

Near Homogeneity of PR2-Bias Fingerprints in the Human Genome and Their Implications in Phylogenetic Analyses

Genes of a multicellular organism are heterogeneous in the G+C content, which is particularly true in the third codon position. The extent of deviation from intra-strand equality rule of A = T and G = C (Parity Rule 2, or PR2) is specific for individual amino acids and has been expressed as the PR2-bias fingerprint. Previous results suggested that the PR2-bias fingerprints tend to be similar among the genes of an organism, and the fingerprint of the organism is specific for different taxa, reflecting phylogenetic relationships of organisms. In this study, using coding sequences of a large number of human genes, we examined the intragenomic heterogeneity of their PR2-bias fingerprints in relation to the G+C content of the third codon position (P ₃ ). Result shows that the PR2-bias fingerprint is similar in the wide range of the G+C content at the third codon position (0.30–0.80). This range covers approximately 89% of the genes, and further analysis of the high G+C range (0.80–1.00), where genes with normal PR2-bias fingerprints and those with anomalous fingerprints are mixed, shows that the total of 95% of genes have the similar finger prints. The result indicates that the PR2-bias fingerprint is a unique property of an organism and represents the overall characteristics of the genome. Combined with the previous results that the evolutionary change of the PR2-bias fingerprint is a slow process, PR2-bias fingerprints may be used for the phylogenetic analyses to supplement and augment the conventional methods that use the differences of the sequences of orthologous proteins and nucleic acids. Potential advantages and disadvantages of the PR2-bias fingerprint analysis are discussed. Received: 21 December 2000 / Accepted: 16 February 2001     

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     

A Phylogenetic Assessment of the Eukaryotic Light-Harvesting Antenna Proteins, with Implications for Plastid Evolution

The light-harvesting complexes (LHCs) are a superfamily of chlorophyll-binding proteins present in all photosynthetic eukaryotes. The Lhc genes are nuclear-encoded, yet the pigment–protein complexes are localized to the thylakoid membrane and provide a marker to follow the evolutionary paths of plastids with different pigmentation. The LHCs are divided into the chlorophyll a/b-binding proteins of the green algae, euglenoids, and higher plants and the chlorophyll a/c-binding proteins of various algal taxa. This work examines the phylogenetic position of the LHCs from three additional taxa: the rhodophytes, the cryptophytes, and the chlorarachniophytes. Phylogenetic analysis of the LHC sequences provides strong statistical support for the clustering of the rhodophyte and cryptomonad LHC sequences within the chlorophyll a/c-binding protein lineage, which includes the fucoxanthin–chlorophyll proteins (FCP) of the heterokonts and the intrinsic peridinin–chlorophyll proteins (iPCP) of the dinoflagellates. These associations suggest that plastids from the heterokonts, haptophytes, cryptomonads, and the dinoflagellate, Amphidinium, evolved from a red algal-like ancestor. The Chlorarachnion LHC is part of the chlorophyll a/b-binding protein assemblage, consistent with pigmentation, providing further evidence that its plastid evolved from a green algal secondary endosymbiosis. The Chlorarachnion LHC sequences cluster with the green algal LHCs that are predominantly associated with photosystem II (LHCII). This suggests that the green algal endosymbiont that evolved into the Chlorarachnion plastid was acquired following the emergence of distinct LHCI and LHCII complexes. Received: 25 February 1998 / Accepted: 13 May 1998     

Associated Proteins and Renal Epithelial Na+ Channel Function

The hypothesis that amiloride-sensitive Na⁺ channel complexes immunopurified from bovine renal papillary collecting tubules contain, as their core conduction component, an ENaC subunit, was tested by functional and immunological criteria. Disulfide bond reduction with dithiothreitol (DTT) of renal Na⁺ channels incorporated into planar lipid bilayers caused a reduction of single channel conductance from 40 pS to 13 pS, and uncoupled PKA regulation of this channel. The cation permeability sequence, as assessed from bi-ionic reversal potential measurements, and apparent amiloride equilibrium dissociation constant (K ^amil _i ) of the Na⁺ channels were unaltered by DTT treatment. Like ENaC, the DTT treated renal channel became mechanosensitive, and displayed a substantial decrease in K ^amil _i following stretch (0.44 ± 0.12 μm versus 6.9 ± 1.0 μm). Moreover, stretch activation induced a loss in the channel's ability to discriminate between monovalent cations, and even allowed Ca²⁺ to permeate. Polyclonal antibodies generated against a fusion protein of αbENaC recognized a 70 kDa polypeptide component of the renal Na⁺ channel complex. These data suggest that ENaC is present in the immunopurified renal Na⁺ channel protein complex, and that PKA sensitivity is conferred by other associated proteins. Received: 5 June 1995/Revised: 29 September 1995     

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 the Friedreich Ataxia GAA Trinucleotide Repeat

Friedreich ataxia is an autosomal recessive neurodegenerative disorder associated with a GAA repeat expansion in the first intron of the gene (FRDA) encoding a novel, highly conserved, 210 amino acid protein known as frataxin. Normal variation in repeat size was determined by analysis of more than 600 DNA samples from seven human populations. This analysis showed that the most frequent allele had nine GAA repeats, and no alleles with fewer than five GAA repeats were found. The European and Syrian populations had the highest percentage of alleles with 10 or more GAA repeats, while the Papua New Guinea population did not have any alleles carrying more than 10 GAA repeats. The distributions of repeat sizes in the European, Syrian, and African American populations were significantly different from those in the Asian and Papua New Guinea populations (p < 0.001). The GAA repeat size was also determined in five nonhuman primates. Samples from 10 chimpanzees, 3 orangutans, 1 gorilla, 1 rhesus macaque, 1 mangabey, and 1 tamarin were analyzed. Among those primates belonging to the Pongidae family, the chimpanzees were found to carry three or four GAA repeats, the orangutans had four or five GAA repeats, and the gorilla carried three GAA repeats. In primates belonging to the Cercopithecidae family, three GAA repeats were found in the mangabey and two in the rhesus macaque. However, an AluY subfamily member inserted in the poly(A) tract preceding the GAA repeat region in the rhesus macaque, making the amplified sequence approximately 300 bp longer. The GAA repeat was also found in the tamarin, suggesting that it arose at least 40 million years ago and remained relatively small throughout the majority of primate evolution, with a punctuated expansion in the human genome. Received: 18 August 2000 / Accepted: 10 November 2000     

Phylogenetic Analysis of Three Lipocalin-Like Proteins Present in the Milk of Trichosurus vulpecula (Phalangeridae, Marsupialia)

Three proteins have been identified in the milk of the common brush tail possum, Trichosurus vulpecula that from sequence analysis are members of the lipocalin family. They include β-lactoglobulin, which appears to have two forms; a homologue to the late-lactation protein found in tammar, Macropus eugenii; milk; and a novel protein termed trichosurin. Whereas β-lactoglobulin and trichosurin are both expressed throughout lactation, the late-lactation protein is not detected in samples taken before days 100–110 of lactation. The cDNAs encoding each of these proteins have been isolated from cDNA libraries prepared using possum mammary mRNA and sequenced. Phylogenetic analysis showed that the T. vulpeculaβ-lactoglobulin, along with two other macropod β-lactoglobulins, forms a subclass of β-lactoglobulins distinct from those for eutherian mammals; both marsupial late-lactation proteins appear to have similarities to a family of odorant-binding proteins, whereas trichosurin has similarities to the major urinary proteins of rodents. Received: 28 October 1996 / Accepted: 19 May 1997     

大麦黄花叶病毒属成员UTR系统进化树分析及编码蛋白跨膜结构推测

对大麦黄花叶病毒属成员的同源性和系统进化树分析表明,同一病毒RNA1和2基因组5′-UTR区域在进化上的相关性高于不同病毒同一RNA组分间的关系,而3′-UTR则相反.蛋白质二级结构分析显示RNA1编码的P3和14K蛋白存在跨膜结构.BaMMV ALS1分离物P3蛋白跨膜结构在大肠杆菌中原核表达时有致死作用.类比Potyviridae科其它属成员功能,此2个蛋白可能与膜附着功能相关.RNA2编码的P2蛋白存在两个结构保守的跨膜区域,一些摩擦接种保存的大麦和性花叶病毒(BaMMV)分离物和燕麦花叶病毒(OMV) P2蛋白此两个结构(或一个)缺失后,不能由介体禾谷多粘菌(P. graminis)传播,揭示P2蛋白跨膜结构与禾谷多粘菌传播能力相关.     

The Family of Major Royal Jelly Proteins and Its Evolution

A cDNA encoding a new member of the gene family of major royal jelly proteins (MRJPs) from the honeybee, Apis mellifera, was isolated and sequenced. Royal jelly (RJ) is a secretion of the cephalic glands of nurse bees. The origin and biological function of the protein component (12.5%, w/w) of RJ is unknown. We show that the MRJP gene family encodes a group of closely related proteins that share a common evolutionary origin with the yellow protein of Drosophila melanogaster. Yellow protein functions in cuticle pigmentation in D. melanogaster. The MRJPs appear to have evolved a novel nutritional function in the honeybee. Received: 26 September 1998 / Accepted: 28 February 1999     

Phylogenetic Relationships of Fungi, Plantae, and Animalia Inferred from Homologous Comparison of Ribosomal Proteins

The complete set of available ribosomal proteins was utilized, at both the peptidic and the nucleotidic level, to establish that plants and metazoans form two sister clades relative to fungi. Different phylogenetic inference methods are applied to the sequence data, using archeans as the outgroup. The evolutionary length of the internal branch within the eukaryotic crown trichotomy is demonstrated to be, at most, one-tenth of the evolutionary length of the branch leading to the cenancester of these three kingdoms. Received: 1 November 1997 / Accepted: 7 January 1998