Molecular Evolutionary Analysis of the Thiamine-Diphosphate-Dependent Enzyme,Transketolase

Members of the transketolase group of thiamine-diphosphate-dependent enzymes from 17 different organisms including mammals, yeast, bacteria, and plants have been used for phylogenetic reconstruction. Alignment of the amino acid and DNA sequences for 21 transketolase enzymes and one putative transketolase reveals a number of highly conserved regions and invariant residues that are of predicted importance for enzyme activity, based on the crystal structure of yeast transketolase. One particular sequence of 36 residues has some similarities to the nucleotide-binding motif and we designate it as the transketolase motif. We report further evidence that the recP protein from Streptococcus pneumoniae might be a transketolase and we list a number of invariant residues which might be involved in substrate binding. Phylogenies derived from the nucleotide and the amino acid sequences by various methods show a conventional clustering for mammalian, plant, and gram-negative bacterial transketolases. The branching order of the gram-positive bacteria could not be inferred reliably. The formaldehyde transketolase (sometimes known as dihydroxyacetone synthase) of the yeast Hansenula polymorpha appears to be orthologous to the mammalian enzymes but paralogous to the other yeast transketolases. The occurrence of more than one transketolase gene in some organisms is consistent with several gene duplications. The high degree of similarity in functionally important residues and the fact that the same kinetic mechanism is applicable to all characterized transketolase enzymes is consistent with the proposition that they are all derived from one common ancestral gene. Transketolase appears to be an ancient enzyme that has evolved slowly and might serve as a model for a molecular clock, at least within the mammalian clade. Received: 13 September 1995 / Accepted: 14 November 1996     

Gene structure of a chlorophyll a/c-binding protein from a brown alga: Presence of an intron and phylogenetic implications

A Laminaria saccharina genomic library in the phage EMBL 4 was used to isolate and sequence a full-length gene encoding a fucoxanthin-chlorophyll a/c-binding protein. Contrary to diatom homologues, the coding sequence is interrupted by an intron of about 900 bp which is located in the middle of the transit peptide. The deduced amino acid sequence of the mature protein is very similar to those of related proteins from Macrocystis pyrifera (Laminariales) and, to a lesser extent, to those from diatoms and Chrysophyceae. Seven of the eight putative chlorophyll-binding amino acids determined in green plants are also present. Alignments of different sequences related to the light-harvesting proteins (LHC) demonstrate a structural similarity among the three transmembrane helices and suggest a unique ancestral helix preceded by two β-turns. The β-turns are conserved in front of the second helices of the chlorophyll a/c proteins more so than in chlorophyll a/b proteins. Phylogenetic trees generated from sequence data indicate that fucoxanthin-chlorophyll-binding proteins diverged prior to the separation of photosystem I and photosystem II LHC genes of green plants. Among the fucoxanthin-containing algae, LHC I or II families could not be distinguished at this time. Received: 14 February 1996 / Accepted: 4 April 1996     

Apolipophorin II/I, Apolipoprotein B, Vitellogenin, and Microsomal Triglyceride Transfer Protein Genes Are Derived from a Common Ancestor

Large lipid transfer proteins (LLTP) are nonexchangeable apolipoproteins and intracellular lipid-exchange proteins involved in the assembly, secretion, and metabolism of lipoproteins. We have identified contiguous conserved sequence motifs in alignments of insect apolipophorin II/I precursor (apoLp-II/I), human apolipoprotein B (apoB), invertebrate and vertebrate vitellogenins (VTG), and the large subunit of mammalian microsomal triglyceride transfer protein (MTP). Conserved motifs present in the N-terminal part of nonexchangeable apolipoproteins encompass almost completely the large subunit of MTP, suggesting a derivation from a common ancestral functional unit, termed large lipid transfer (LLT) module. Divergence of LLTP from a common ancestor is supported by (1) the statistical significance of the combined match scores obtained after motif-based database searches, (2) the presence of several identical amino acid residues in all LLTP sequences currently available, (3) the conservation of hydrophobic clusters in an α-helical domain, (4) the phylogenetic analysis of the conserved sequences related to the von Willebrand factor D (VWD) module identified in nonexchangeable apolipoproteins, and (5) the presence of four and one ancestral exon boundaries in the LLT and VWD modules, respectively. Our data indicate that the genes coding for apoLp-II/I, apoB, VTG, and the MTP large subunit are members of the same multigene superfamily. LLTP have emerged from an ancestral molecule designed to ensure a pivotal event in the intracellular and extracellular transfer of lipids and liposoluble substances. Received: 8 June 1998 / Accepted: 15 February 1999     

Two isoforms of a member of the arthropod defensin family from the soft tick, Ornithodoros moubata (Acari: Argasidae)

We previously purified and determined the partial amino acid sequence of a 4 kDa peptide having high homology with scorpion defensin from the hemolymph of adult fed female soft ticks, Ornithodoros moubata. In this study, the full length sequences of two defensin isoforms were obtained. Deduced amino acid sequences reveal a precursor protein of 73 amino acid residues with a mature portion consisting of 37 amino acid residues. This mature peptide contains six cysteine residues conserved in the same location as other invertebrate defensins. Phylogenetic analysis reveals that Ornithodoros defensin is most closely related to scorpion defensin and other more ancient arthropods. Ornithodoros defensin mRNA is constitutively expressed and up-regulated by blood-feeding and bacterial injection. Ornithodoros defensin gene expression occurs mainly in the midgut. This is the first report of the cloning and gene expression of an antibacterial peptide from the Acari.     

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     

The Molecular Evolution of Visual Pigments of Freshwater Crayfishes (Decapoda: Cambaridae)

This study examines the diverse maximum wavelength absorption (λ_max) found in crayfishes (Decapoda: Cambaridae and Parastacidae) and the associated genetic variation in their opsin locus. We measured the wavelength absorption in the photoreceptors of six species that inhabit environments of different light intensities (i.e., burrows, streams, standing waters, and subterranean waters). Our results indicate that there is relatively little variation in λ_max (522–530 nm) among species from different genera and families. The existing variation did not correlate with the habitat differences of the crayfishes studied. We simultaneously sequenced the rhodopsin gene to identify the amino acid replacements that affect shifts in maximum wavelength absorption. We then related these to changes that correlated with shifts in λ_max by reconstructing ancestral character states using a maximum-likelihood approach. Using amino acid sequences obtained from five species (all were 301 amino acids in length), we identified a number of candidates for producing shifts of 4 to 8 nm in λ_max. These amino acid replacements occurred in similar regions to those involved in spectral shifts in vertebrates. Received: 12 March 1997 / Accepted: 3 June 1997     

Characterization of β-defensin prepropeptide mRNA from chicken and turkey bone marrow

Four avian β-defensin prepropeptide cDNA sequences [gallinacins: Gal 1 (synonym CHP 1, chicken heterophil peptide 1), and Gal 2; turkey heterophil peptides: THP 1 and THP 2] were amplified from chicken or turkey bone marrow mRNA samples, respectively. Partial chicken β-defensin cDNA sequences were obtained using degenerate primers based on chicken peptide sequences (Gal 1/CHP 1 and Gal 2). The complete cDNA sequences of the chicken β-defensins were then determined by designing specific intrapeptidal primers, from the newly acquired sequence, and pairing one primer with a specific poly A primer tail sequence (3' end) and the other primer with an adapter primer in a 5' rapid amplification of cDNA ends (RACE) reaction. The two, turkey β-defensins were amplified from turkey marrow using primers designed from chicken β-defensin preproregions. The complete amino acid sequences for the prepropeptides were deduced for all four avian β-defensins. Previously, only partial mature peptide sequences for the turkey β-defensins and complete mature peptide sequences for the chicken β-defensins were known. All sequences obtained translated accurately to complete and partial amino acid sequences reported for β-defensins purified from chicken and turkey heterophil granules except for one additional amino acid for Gal 1/CHP 1. The four deduced β-defensin proregions lack the long, negatively charged propiece reported in classical defensin proregions. These regions are thought to stabilize and inactivate the positively charged mature peptide and target the propeptide to the storage granule. Instead, these β-defensin proregions are shorter and similar to storage granule-free β-defensins proregions reported for bovine tracheal antimicrobial peptide (TAP) and lingual antimicrobial peptide (LAP). These are the first prepropeptide β-defensins from leukocyte granules to be completely characterized.     

The Origin of Trypsin: Evidence for Multiple Gene Duplications in Trypsins

The trypsin family of serine proteases is one of the most studied protein families, with a wealth of amino acid sequence information available in public databases. Since trypsin-like enzymes are widely distributed in living organisms in nature, likely evolutionary scenarios have been proposed. A novel methodology for Fourier transformation of biological sequences (FOTOBIS) is presented. The methodology is well suited for the identification of the size and extent of short repeats in protein sequences. In the present paper the trypsin family of enzymes is analyzed with FOTOBIS and strong evidence for tandem gene duplication is found. A likely evolutionary path for the development of present-day trypsins involved an intrinsic extensive tandem gene duplication of a small DNA fragment of 15–18 nucleotides, corresponding to five or six amino acids. This ancestral trypsin gene was subsequently duplicated, leading to the earliest version of a full-sized trypsin, from which the contemporary trypsins have developed. Received: 22 November 1997 / Accepted: 26 January 1998     

Molecular Evolution of the Domain Structures of Protein Disulfide Isomerases

Protein disulfide isomerase (PDI) is an enzyme that promotes protein folding by catalyzing disulfide bridge isomerization. PDI and its relatives form a diverse protein family whose members are characterized by thioredoxin-like (TX) domains in the primary structures. The family was classified into four classes by the number and the relative positions of the TX domains. To investigate the evolution of the domain structures, we aligned the amino acid sequences of the TX domains, and the molecular phylogeny was examined by the NJ and ML methods. We found that all of the current members of the PDI family have evolved from an ancestral enzyme, which has two TX domains in the primary structure. The diverse domain structures of the members have been generated through domain duplications and deletions.     

Molecular Evolution of the Photolyase–Blue-Light Photoreceptor Family

The photolyase–blue-light photoreceptor family is composed of cyclobutane pyrimidine dimer (CPD) photolyases, (6-4) photolyases, and blue-light photoreceptors. CPD photolyase and (6-4) photolyase are involved in photoreactivation for CPD and (6-4) photoproducts, respectively. CPD photolyase is classified into two subclasses, class I and II, based on amino acid sequence similarity. Blue-light photoreceptors are essential light detectors for the early development of plants. The amino acid sequence of the receptor is similar to those of the photolyases, although the receptor does not show the activity of photoreactivation. To investigate the functional divergence of the family, the amino acid sequences of the proteins were aligned. The alignment suggested that the recognition mechanisms of the cofactors and the substrate of class I CPD photolyases (class I photolyases) are different from those of class II CPD photolyases (class II photolyases). We reconstructed the phylogenetic trees based on the alignment by the NJ method and the ML method. The phylogenetic analysis suggested that the ancestral gene of the family had encoded CPD photolyase and that the gene duplication of the ancestral proteins had occurred at least eight times before the divergence between eubacteria and eukaryotes. Received: 23 October 1996 / Accepted: 1 April 1997     

Phylogenetic Analysis of the Homologous Proteins of the Terminal Complement Complex Supports the Emergence of C6 and C7 Followed by C8 and C9

The plasma complement system comprises several activation pathways that share a common terminal route involving the assembly of the terminal complement complex (TCC), formed by C5b–C9. The order of emergence of the homologous components of TCC (C6, C7, C8α, C8β, and C9) has been determined by phylogenetic analyses of their amino acid sequences. Using all the sequence data available for C6–C9 proteins, as well as for perforins, the results suggested that these TCC components originated from a single ancestral gene and that C6 and C7 were the earliest to emerge. Our evidence supports the notion that the ancestral gene had a complex modular composition. A series of gene duplications in combination with a tendency to lose modules resulted in successive complement proteins with decreasing modular complexity. C9 and perforin apparently are the result of different selective conditions to acquire pore-forming function. Thus C9 and perforin are examples of evolutionary parallelism. Received: 16 August 1998 / Accepted: 12 March 1999     

Evolution of Histone H4 and H3 Genes in Different Ciliate Lineages

The histones H4 are known as highly conserved proteins. However, in ciliates a high degree of variation was found compared both to other eukaryotes and between the ciliate species. To date, only H4 histones of species belonging to two distantly related classes have been investigated. In order to obtain more detailed information on histone H4 variation in ciliates we undertook a comprehensive sequence analysis of PCR-amplified internal H4 fragments from 12 species belonging to seven out of the nine currently recognized ciliate classes. In addition, we used PCR primers to amplify longer fragments of H3 and H4 genes including the intergenic region. The encoded amino acid sequences reveal a high number of differences when compared with those of other eukaryotes and the ciliate species investigated. Furthermore, in some species H4 gene variants were detected, which result in amino acid differences. The greatest number of substitutions and insertions found was in the amino terminal region of the H4 histones. However, all sequences possess a conserved region corresponding to those of all other eukaryotic H4 histones. The histone gene variations were used to reconstruct phylogenetic relationships. The tree from our data matches perfectly with the ribosomal RNA data: The heterotrichs, which were considered as a late branching lineage, diverge at the base of the ciliate tree and groups formerly thought to represent ancestral lineages now appear as highly derived ciliates. Received: 4 April 1997 / Accepted: 1 August 1997     

The Evolutionary History of Prosaposin: Two Successive Tandem-Duplication Events Gave Rise to the Four Saposin Domains in Vertebrates

Prosaposin is a multifunctional protein encoded by a single-copy gene. It contains four saposin domains (A, B, C, and D) occurring as tandem repeats connected by linker sequences. Because the saposin domains are similar to one another, it is deduced that they were created by sequential duplications of an ancestral domain. There are two types of evolutionary scenarios that may explain the creation of the four-domain gene: (1) two rounds of tandem internal gene duplication and (2) three rounds of duplications. An evolutionary and phylogenetic analysis of saposin DNA and amino acid sequences from human, mouse, rat, chicken, and zebrafish indicates that the first evolutionary scenario is the most likely. Accordingly, an ancestral saposin-unit duplication produced a two-domain gene, which, subsequently, underwent a second complete tandem duplication to give rise to the present four-domain structure of the prosaposin gene. Received: 8 February 2001 / Accepted: 29 June 2001     

Origin and Evolution of Viral Interleukin-10 and Other DNA Virus Genes with Vertebrate Homologues

Phylogenies of gene families including members in both vertebrates and DNA viruses of the poxvirus and/or herpesvirus families showed that the viral genes originated at widely different times over the history of life. Certain of these viral genes (for example, the genes encoding the large and small subunits of ribonucleoside–diphosphate reductase) originated before animals diverged from fungi, while others originated much more recently. The most striking examples of recent origin involved viral genes encoding the cytokine interleukin-10 (IL-10), which originated independently in viruses at least three times since the divergence of the orders of eutherian mammals, presumably by viral capture of host genes. In certain domains, viral IL-10 genes showed significantly higher rates of nonsynonymous substitution than their nearest mammalian homologues. Though the mutation rate in these viral genes is up to 20 times that of the corresponding mammalian genes, a high mutation rate alone did not account for these differences because they were not seen in all domains. Rather, in certain domains it appears that functional constraints present in the case of mammalian IL-10 are relaxed in the case of the viral homologues. Furthermore, a nonrandom pattern of change with respect to amino acid residue charge in the N-terminal portion of the mature protein has occurred repeatedly in independently derived viral IL-10 genes, strongly suggesting that positive selection has led to divergence of this functionally important domain in viral IL-10. Received: 11 January 2001 / Accepted: 23 May 2001     

Multiple β-defensin isoforms identified in early developmental stages of the teleost Paralichthys olivaceus

The β-defensin-like gene and its cloned isoforms (fBDI-1 to -5) were identified in an expressed sequence tag (EST) library from the early developmental stages of the olive flounder, Paralichthys olivaceus. The fBDI cDNA clones show identical amino acid sequences in 24 residues of the signal peptide and 38 residues of the mature peptide; however, the propiece region varies in sequence and length, from 5 to 15 amino acid residues. The predicted molecular weight of the mature peptide is 3.83 kDa, and its predicted isoelectric point is 4.1, showing anionic properties. The genomic organisation of the isoforms was analysed using bacterial artificial chromosome (BAC) DNA containing the fBDI gene. Southern blotting and sequence analyses of fBDI BAC DNA confirmed that the fBDI isoforms cluster at the same locus and exhibit the conserved gene organisation reported for other fish defensin genes. The fBDI mRNA was expressed constitutively in early developmental stages after hatching, and pathogen challenge induced fBDI expression in the head kidney of juvenile fish. We also produced a recombinant fBDI peptide (smfBD) using the expression plasmid pET32 and examined its bioactivity toward Escherichia coli.     

Rates of Conservative and Radical Nonsynonymous Nucleotide Substitutions in Mammalian Nuclear Genes

To understand the process and mechanism of protein evolution, it is important to know what types of amino acid substitutions are more likely to be under selection and what types are mostly neutral. An amino acid substitution can be classified as either conservative or radical, depending on whether it involves a change in a certain physicochemical property of the amino acid. Assuming Kimura's two-parameter model of nucleotide substitution, I present a method for computing the numbers of conservative and radical nonsynonymous (amino acid altering) nucleotide substitutions per site and estimate these rates for 47 nuclear genes from mammals. The results are as follows. (1) The average radical/conservative rate ratio is 0.81 for charge changes, 0.85 for polarity changes, and 0.49 when both polarity and volume changes are considered. (2) The radical/conservative rate ratio is positively correlated with the nonsynonymous/synonymous rate ratio for charge changes or when both polarity and volume changes are considered. (3) Both the conservative/synonymous rate ratio and the radical/synonymous rate ratio are lower in the rodent lineage than in the primate or artiodactyl lineage, suggesting more intense purifying selection in the rodent lineage, for both conservative and radical nonsynonymous substitutions. (4) Neglecting transition/transversion bias would cause an underestimation of both radical and conservative rates and the ratio thereof. (5) Transversions induce more dramatic genetic alternations than transitions in that transversions produce more amino acid altering changes and among which, more radical changes. Received: 6 April 1999 / Accepted: 16 August 1999     

Evolutionary Analysis of the ErbB Receptor and Ligand Families

We have compared all available deduced protein sequences of the ErbB family of receptors and their ligands. Analysis of the aligned sequences of the receptors indicates that there are some differences in the receptors that are specific to invertebrates. In addition, comparison of the vertebrate ErbB receptors suggest that a gene duplication event generated two ancestral receptors, the ErbB3/ErbB4 precursor and the ErbB1/ErbB2 precursor. Subsequent gene duplications of these precursors generated the four receptors present in mammals. Analysis of the sequences for the known ligands of the ErbB receptors suggests that the vertebrate ligands segregate into the ErbB1 ligands and the ErbB3/ErbB4 ligands, paralleling the evolution of the receptors; however, it is difficult to ascertain any correlation between the invertebrate and the vertebrate ligands. Even though ErbB3 is kinase-impaired, there is significant conservation of the kinase domain within the vertebrate lineage (human, rat, and F. rubripes), suggesting some function for this domain other than kinase activity, such as mediating protein–protein interactions that are involved in receptor dimerization and/or activation of the kinase domain of the heterodimerization partner. To date, no ligand for ErbB2 has been identified, and comparison of the extracellular domains of ErbB2 reveals two regions that are not conserved across the mammalian species. These two regions of divergence align with sequences in ErbB1 that have been shown to be proximal to the amino-terminus and to the carboxyl-terminal region, respectively, of bound EGF. Further, one of these regions contains an insertion, relative to the other members of the mammalian ErbB family, which might affect the ligand binding site and provide a structural basis for this receptor's apparent inability to bind ligand independently. Received: 8 September 1999 / Accepted: 17 January 2000