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     

Analysis of Orthologous Retrovirus-Like Elements in the White-Footed Mouse,Peromyscus leucopus

Three loci in the genome of the white-footed mouse, Peromyscus leucopus, were examined for the presence or absence of orthologous copies of the retrovirus-like element mys using polymerase chain reaction. We examined these loci in 28 mice collected throughout the P. leucopus species range. Mys insertions were present in only one of the individuals examined at the mys-1 and mys-7 loci. Conversely, the mys-6 element was found in several individuals, but the presence of this element was limited to northern latitudes. Because the long terminal repeats (LTRs) of a given element are expected to be identical at the time of retrotransposition into the genome, and to accumulate changes over evolutionary time, within-element LTR sequence comparisons can be used to estimate the relative age of insertions. Within-element LTR differences are greater in mys-6 than in mys-1 or mys-7. The LTRs from orthologous mys-6 elements of six mice were sequenced. The alignment revealed 13 of the 22 differences between the right and left LTRs that were shared by all orthologous mys-6 sites, suggesting that relative to its time of insertion into the genome, mys-6 has only recently spread across the northern part of the species range. Received: 23 January 1996 / Accepted: 24 April 1996     

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     

Synonymous Nucleotide Divergence and Saturation: Effects of Site-Specific Variations in Codon Bias and Mutation Rates

The synonymous divergence between Escherichia coli and Salmonella typhimurium is explained in a model where there is a large variation between mutation rates at different nucleotide sites in the genome. The model is based on the experimental observation that spontaneous mutation rates can vary over several orders of magnitude at different sites in a gene. Such site-specific variation must be taken into account when studying synonymous divergence and will result in an apparent saturation below the level expected from an assumption of uniform rates. Recently, it has been suggested that codon preference in enterobacteria has a very large site-specific variation and that the synonymous divergence between different species, e.g., E. coli and Salmonella, is saturated. In the present communication it is shown that when site-specific variation in mutation rates is introduced, there is no need to invoke assumptions of saturation and a large variability in codon preference. The same rate variation will also bring average mutation rates as estimated from synonymous sequence divergence into numerical agreement with experimental values. Received: 10 July 1998 / Accepted: 20 August 1998     

Evolution of the nad3-rps12 Gene Cluster in Angiosperm Mitochondria: Comparison of Edited and Unedited Sequences

We have analyzed the nad3-rps12 locus for eight angiosperms in order to compare the utility of mitochondrial DNA and edited mRNA sequences in phylogenetic reconstruction. The two coding regions, containing from 25 to 35 editing sites in the various plants, have been concatenated in order to increase the significance of the analysis. Differing from the corresponding chloroplast sequences, unedited mitochondrial DNA sequences seem to evolve under a quasi-neutral substitution process which undifferentiates the nucleotide substitution rates for the three codon positions. By using complete gene sequences (all codon positions) we found that genomic sequences provide a classical angiosperm phylogenetic tree with a clear-cut grouping of monocotyledons and dicotyledons with Magnoliidae at the basal branch of the tree. Conversely, owing to their low nucleotide substitution rates, edited mRNA sequences were found not to be suitable for studying phylogenetic relationships among angiosperms. Received: 24 January 1996 / Accepted: 5 June 1996     

A theoretical method for evaluating the relative importance of positive selection and neutral drift from observed base changes

To evaluate the relative importance of positive selection and neutral drift from the nucleotide base changes observed in the homologous alignment of genes, a theoretical equation of base changes is formulated by including both the influence of selection and the base substitutions due to mutations. Under the assumption that the average rate of base substitutions estimated from synonymous changes is the ``true' mutation rate applicable at all positions, this method is applied to the vertebrate globin gene family, and evaluates the departures of base change rates from the ``true' mutation rate at the first and second codon positions as a consequence of preferential selection for the conservation of important function. In addition to the strong effect of selection on the amino acid residues in the internal region mostly common to myoglobin and hemoglobin chains, the distinctive directions of selective parameter values are seen at sites on the globin surface, distinguishing the subunit contact residues of hemoglobins from the polar residues on the surface of myoglobins. Moreover, this effect of selection distinguishing between the myoglobin and hemoglobin chain genes becomes weaker in cold-blooded vertebrates, especially in fish, strongly suggesting the possibility that the clear distinction between these globins is a result of selection out of the changes regarded as neutral ones in an ancestor of vertebrates. Thus, the present method may also serve to investigate the homology of many other proteins from the aspect of molecular evolution, mainly focusing on the evolution of their biological functions. Received: 2 January 1996 / Accepted: 20 February 1997     

Selection and Methionine Accumulation in the Fat Body Protein 2 Gene (FBP2), a Duplicate of the Drosophila Alcohol Dehydrogenase (ADH) Gene

The Drosophila fat body protein 2 gene (Fbp2) is an ancient duplication of the alcohol dehydrogenase gene (Adh) which encodes a protein that differs substantially from ADH in its methionine content. In D. melanogaster, there is one methionine in ADH, while there are 51 (20% of all amino acids) in FBP2. Methionine is involved in 46% of amino acid replacements when Fbp2 DNA sequences are compared between D. melanogaster and D. pseudoobscura. Methionine accumulation does not affect conserved residues of the ADH-ADH^r-FBP2 multigene family. The multigene family has evolved by replacement of mildly hydrophobic amino acids by methionine with no apparent reversion. Its short-term evolution was compared between two Drosophila species, while its long-term evolution was compared between two genera belonging respectively to acalyptrate and calyptrate Diptera, Drosophila and Sarcophaga. The pattern of nucleotide substitution was consistent with an independent accumulation of methionines at the Fbp2 locus in each lineage. Under a steady-state model, the rate of methionine accumulation was constant in the lineage leading to Drosophila, and was twice as fast as that in the calyptrate lineage. Substitution rates were consistent with a slight positive selective advantage for each methionine change in about one-half of amino acid sites in Drosophila. This shows that selection can potentially account for a large proportion of amino acid replacements in the molecular evolution of proteins. Received: 12 December 1994 / Accepted: 15 April 1996     

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     

Nucleotide Composition Bias Affects Amino Acid Content in Proteins Coded by Animal Mitochondria

We show that in animal mitochondria homologous genes that differ in guanine plus cytosine (G + C) content code for proteins differing in amino acid content in a manner that relates to the G + C content of the codons. DNA sequences were analyzed using square plots, a new method that combines graphical visualization and statistical analysis of compositional differences in both DNA and protein. Square plots divide codons into four groups based on first and second position A + T (adenine plus thymine) and G + C content and indicate differences in amino acid content when comparing sequences that differ in G + C content. When sequences are compared using these plots, the amino acid content is shown to correlate with the nucleotide bias of the genes. This amino acid effect is shown in all protein-coding genes in the mitochondrial genome, including cox I, cox II, and cyt b, mitochondrial genes which are commonly used for phylogenetic studies. Furthermore, nucleotide content differences are shown to affect the content of all amino acids with A + T- and G + C-rich codons. We speculate that phylogenetic analysis of genes so affected may tend erroneously to indicate relatedness (or lack thereof) based only on amino acid content. Received: 3 July 1996 / Accepted: 6 November 1996     

Relatedness and Phylogeny Within the Family of Periplasmic Chaperones Involved in the Assembly of Pili or Capsule-Like Structures of Gram-Negative Bacteria

The structure of a Salmonella enterica serovar typhi gene located within the fim gene cluster and encoding a putative periplasmic chaperone-like protein involved in the assembly of type 1 pili was determined. This gene, named fimC, has the ability to encode a 26-kDa polypeptide which is similar, at the sequence level, to the PapD periplasmic chaperonin mediating the assembly of P pili of Escherichia coli, as well as to other periplasmic chaperone-like proteins involved in the biogenesis of pili or capsule-like structures of various Gram-negative bacteria. A comprehensive search through the literature and sequence databases identified 31 (putative) bacterial proteins that can be included in this protein family on the basis of sequence similarity. Results of a multiple sequence comparison analysis showed that several residues, including most of those known to be critical in maintaining the three-dimensional structure of PapD, are either conserved or conservatively substituted in all these proteins, suggesting an overall similar folding for all of them. It was also evident that members of this family are clustered into different subfamilies according to structural and phyletic data. Received: 15 February 1996 / Accepted: 3 October 1996     

Evolution of the Larval Cuticle Proteins Coded by the Secondary Sex Chromosome Pair: X2 and Neo-Y of Drosophila miranda: II. Comparison at the Amino Acid Sequence Level

The larval cuticle proteins (LCPs) are encoded by a multigene family, Lcp1–4, located at the right arm of the metacentric autosome 2 (2R) in Drosophila melanogaster. Due to a chromosome fusion the Lcp locus of Drosophila miranda is situated on a pair of secondary sex chromosomes, the X2 and neo-Y chromosomes. Comparing the deduced amino acid sequences of the autosomal D. melanogaster loci with the sex-chromosomal loci of D. miranda, we were able to trace the evolution of the Lcp loci with respect to their different chromosomal inheritance. The length of the signal peptide is conserved in all four LCPs, while the size of the mature LCPs varies. Conserved protein motifs became obvious from the alignment, indicating regions of structural and functional importance. Analyzing intra- and interspecific sequence similarities of the Lcp gene families allowed us to reconstruct the phylogeny of the gene cluster. Alignment with cuticle amino acid sequences originating from divergent insect species reveals motifs already present in the primordial insect LCPs. These motifs indicate different levels of constraint acting during the evolution of the LCPs. Received: 27 December 1995 / Accepted: 30 April 1996     

Evolutionary Comparisons of RecA-Like Proteins Across All Major Kingdoms of Living Organisms

Protein sequences with similarities to Escherichia coli RecA were compared across the major kingdoms of eubacteria, archaebacteria, and eukaryotes. The archaeal sequences branch monophyletically and are most closely related to the eukaryotic paralogous Rad51 and Dmc1 groups. A multiple alignment of the sequences suggests a modular structure of RecA-like proteins consisting of distinct segments, some of which are conserved only within subgroups of sequences. The eukaryotic and archaeal sequences share an N-terminal domain which may play a role in interactions with other factors and nucleic acids. Several positions in the alignment blocks are highly conserved within the eubacteria as one group and within the eukaryotes and archaebacteria as a second group, but compared between the groups these positions display nonconservative amino acid substitutions. Conservation within the RecA-like core domain identifies possible key residues involved in ATP-induced conformational changes. We propose that RecA-like proteins derive evolutionarily from an assortment of independent domains and that the functional homologs of RecA in noneubacteria comprise an array of RecA-like proteins acting in series or cooperatively. Received: 25 October 1996 / Accepted: 31 December 1996     

Multiple origins of fungal group I introns located in the same position of nuclear SSU rRNA gene

The archiascomycetous fungus Protomyces pachydermus has two group I introns within the nuclear small subunit (nSSU) rRNA gene. One of these introns has an internal open reading frame (ORF) that encodes a predicted protein of 228 amino acid residues. On the other hand, Protomyces macrosporus has two group I introns that insert at the same positions as P. pachydermus, which have no ORF. Each alignment was constructed with Protomyces group I introns located in the same position and other introns retrieved by the BLAST Search. Each phylogenetic tree based on the alignment shows that Protomyces introns are monophyletic but the relationships among fungal introns do not reflect on the fungal phylogeny. Therefore, it is suggested that two different horizontal transfers of group I introns occurred at the early stage of Protomyces species diversification. Received: 11 June 1997 / Accepted: 2 September 1997     

Ribonucleotide reductases: divergent evolution of an ancient enzyme

Ribonucleotide reductases (RNRs) are uniquely responsible for converting nucleotides to deoxynucleotides in all dividing cells. The three known classes of RNRs operate through a free radical mechanism but differ in the way in which the protein radical is generated. Class I enzymes depend on oxygen for radical generation, class II uses adenosylcobalamin, and the anaerobic class III requires S-adenosylmethionine and an iron–sulfur cluster. Despite their metabolic prominence, the evolutionary origin and relationships between these enzymes remain elusive. This gap in RNR knowledge can, to a major extent, be attributed to the fact that different RNR classes exhibit greatly diverged polypeptide chains, rendering homology assessments inconclusive. Evolutionary studies of RNRs conducted until now have focused on comparison of the amino acid sequence of the proteins, without considering how they fold into space. The present study is an attempt to understand the evolutionary history of RNRs taking into account their three-dimensional structure. We first infer the structural alignment by superposing the equivalent stretches of the three-dimensional structures of representatives of each family. We then use the structural alignment to guide the alignment of all publicly available RNR sequences. Our results support the hypothesis that the three RNR classes diverged from a common ancestor currently represented by the anaerobic class III. Also, lateral transfer appears to have played a significant role in the evolution of this protein family.     

Effects of Density and Gating of Delayed-Rectifier Potassium Channels on Resting Membrane Potential and its Fluctuations

The aim of this study is to evaluate directly, using a reduced experimental system, the nature of interactions between voltage-gated potassium channels and the resting membrane potential. Xenopus oocytes were injected with various concentrations of cRNA coding for a delayed-rectifier potassium channel Shaker-IR. The effects of the density and kinetics of the expressed channels on resting membrane potential is explored in isolated (``inside-out') patches. The channel density is given in terms of maximal conductance (G _max), measured from the maximal slope of the I-V curve under voltage clamp conditions. The capacitance of the experimental setup is approximately 1 pF. At high channel densities (G _max > 10 pA/mV) the mean membrane potential is stabilized at approximately −60 mV. This resting membrane potential is more than 35 mV positive to the reversal potential for potassium ions under the same experimental conditions. Analyses of voltage clamp experiments indicate that at high channel densities the mean membrane potential is determined by the rates of channel activation and deactivation, but is not affected by the rates involved in the process of slow (C-type) inactivation. In contrast, at lower channel densities membrane potential is very unstable, and its mean value and amplitude of fluctuations are strongly affected by the process of slow (C-type) inactivation. Received: 21 March 1996/Revised: 6 August 1996     

Molecular Evolution of Cytochrome c Oxidase Subunit IV: Evidence for Positive Selection in Simian Primates

Cytochrome c oxidase (COX) is a multi-subunit enzyme complex that catalyzes the final step of electron transfer through the respiratory chain on the mitochondrial inner membrane. Up to 13 subunits encoded by both the mitochondrial (subunits I, II, and III) and nuclear genomes occur in eukaryotic organisms ranging from yeast to human. Previously, we observed a high number of amino acid replacements in the human COX IV subunit compared to mouse, rat, and cow orthologues. Here we examined COX IV evolution in the two groups of anthropoid primates, the catarrhines (hominoids, cercopithecoids) and platyrrhines (ceboids), as well as one prosimian primate (lorisiform), by sequencing PCR-amplified portions of functional COX4 genes from genomic DNAs. Phylogenetic analysis of the COX4 sequence data revealed that accelerated nonsynonymous substitution rates were evident in the early evolution of both catarrhines and, to a lesser extent, platyrrhines. These accelerated rates were followed later by decelerated rates, suggesting that positive selection for adaptive amino acid replacement became purifying selection, preserving replacements that had occurred. The evidence for positive selection was especially pronounced along the catarrhine lineage to hominoids in which the nonsynonymous rate was first faster than the synonymous rate, then later much slower. The rates of three types of ``neutral DNA' nucleotide substitutions (synonymous substitutions, pseudogene nucleotide substitutions, and intron nucleotide substitutions) are similar and are consistent with previous observations of a slower rate of such substitutions in the nuclear genomes of hominoids than in the nuclear genomes of other primate and mammalian lineages. Received: 22 May 1996 / Accepted: 24 November 1996     

Duplication and divergence of the genes of the α-esterase cluster ofDrosophila melanogaster

The α-esterase cluster of D. melanogaster contains 11 esterase genes dispersed over 60 kb. Embedded in the cluster are two unrelated open reading frames that have sequence similarity with genes encoding ubiquitin-conjugating enzyme and tropomyosin. The esterase amino acid sequences show 37–66% identity with one another and all but one have all the motifs characteristic of functional members of the carboxyl/cholinesterase multigene family. The exception has several frameshift mutations and appears to be a pseudogene. Patterns of amino acid differences among cluster members in relation to generic models of carboxyl/cholinesterase protein structure are broadly similar to those among other carboxyl/cholinesterases sequenced to date. However the α-esterases differ from most other members of the family in: their lack of a signal peptide; the lack of conservation in cysteines involved in disulfide bridges; and in four indels, two of which occur in or adjacent to regions that align with proposed substrate-binding sites of other carboxyl/cholinesterases. Phylogenetic analyses clearly identify three simple gene duplication events within the cluster. The most recent event involved the pseudogene which is located in an intron of another esterase gene. However, relative rate tests suggest that the pseudogene remained functional after the duplication event and has become inactive relatively recently. The distribution of indels also suggests a deeper node in the gene phylogeny that separates six genes at the two ends of the cluster from a block of five in the middle. Received: 18 January 1996 / Accepted: 12 March 1996     

Intron Distribution in Ancient Paralogs Supports Random Insertion and Not Random Loss

The intron positions of ten different protein families were examined to determine (the statistical likelihood of) whether spliceosomal introns are the result of random insertion events into previously intronless genes, on the one hand, or the result of random loss from common ancestral introns, on the other. The number of expected matches for the alternative scenarios was calculated for a binomial distribution by considering currently observed introns relative to all possible locations for insertion or loss. Introns occurring at approximately the same location (hereafter called a ``match') were tallied for each of the paired proteins. Matches were identified by their positions in the multiple alignment and were defined as any two introns occurring within a window of 11 possible nucleotide positions, thereby allowing for possible alignment errors and ``intron sliding.' Matches were tallied from the raw data and compared with the expected number of matches for the two different scenarios. The results suggest that the distribution of introns in genes encoding proteins is due to random insertion and not random loss. Received: 8 September 1996 / Accepted: 24 January 1997     

Accumulation of Species-Specific Amino Acid Replacements That Cause Loss of Particular Protein Functions in Buchnera, an Endocellular Bacterial Symbiont

Endosymbiotic bacteria live in animal cells and are transmitted vertically at the time of the host's reproduction. In view of their small and asexual populations with infrequent chances of recombination, these endocellular bacteria are expected to accumulate mildly deleterious mutations. Previous studies showed that the DNA sequences of these bacteria evolved faster than those of free-living bacteria. In this study, we compared all the ORFs of Buchnera, an endocellular bacterial symbiont of aphids, with those of 34 other prokaryotic organisms and estimated the effect of the accelerated evolution of Buchnera on the functions of its proteins. It was revealed that Buchnera proteins contain many mutations at the sites where sequences are conserved in their orthologues in many other organisms. In addition, amino acid replacements at the conserved sites are mostly changes to physicochemically different amino acids. These results suggest that functions and conformations of Buchnera proteins have been seriously impaired or strongly modified. Indeed, extensive loss of functional motifs was observed in some Buchnera proteins. In many Buchnera proteins mutations were not detected evenly throughout each molecule but tended to accumulate in some functional units, possibly leading to loss of specific functions. As Buchnera has an unusual and limited gene repertory, it is conceivable that the manner of interactions among its proteins has been changed, and thus, functional constraints over their amino acid residues have also been changed during evolution. This may account for the loss of some functional units only in the Buchnera proteins. We obtained evidence that amino acid replacements in Buchnera were not always deleterious, but neutral or, in some cases, even positively selected. Received: 14 December 2000 / Accepted: 12 March 2001