Sequence similarity analysis of Escherichia coli proteins: functional and evolutionary implications.

A computer analysis of 2328 protein sequences comprising about 60% of the Escherichia coli gene products was performed using methods for database screening with individual sequences and alignment blocks. A high fraction of E. coli proteins--86%--shows significant sequence similarity to other proteins in current databases; about 70% show conservation at least at the level of distantly related bacteria, and about 40% contain ancient conserved regions (ACRs) shared with eukaryotic or Archaeal proteins. For > 90% of the E. coli proteins, either functional information or sequence similarity, or both, are available. Forty-six percent of the E. coli proteins belong to 299 clusters of paralogs (intraspecies homologs) defined on the basis of pairwise similarity. Another 10% could be included in 70 superclusters using motif detection methods. The majority of the clusters contain only two to four members. In contrast, nearly 25% of all E. coli proteins belong to the four largest superclusters--namely, permeases, ATPases and GTPases with the conserved "Walker-type" motif, helix-turn-helix regulatory proteins, and NAD(FAD)-binding proteins. We conclude that bacterial protein sequences generally are highly conserved in evolution, with about 50% of all ACR-containing protein families represented among the E. coli gene products. With the current sequence databases and methods of their screening, computer analysis yields useful information on the functions and evolutionary relationships of the vast majority of genes in a bacterial genome. Sequence similarity with E. coli proteins allows the prediction of functions for a number of important eukaryotic genes, including several whose products are implicated in human diseases.

Systematic derivation of partition functions for ligand binding to two-dimensional lattices

The Ising problem consists in finding the analytical solution of the partition function of a lattice once the interaction geometry among its elements is specified. No general analytical solution is available for this problem, except for the one-dimensional case. Using site-specific thermodynamics, it is shown that the partition function for ligand binding to a two-dimensional lattice can be obtained from those of one-dimensional lattices with known solution. The complexity of the lattice is reduced recursively by application of a contact transformation that involves a relatively small number of steps. The transformation implemented in a computer code solves the partition function of the lattice by operating on the connectivity matrix of the graph associated with it. This provides a powerful new approach to the Ising problem, and enables a systematic analysis of two-dimensional lattices that model many biologically relevant phenomena. Application of this approach to finite two-dimensional lattices with positive cooperativity indicates that the binding capacity per site diverges as Na (N = number of sites in the lattice) and experiences a phase-transition-like discontinuity in the thermodynamic limit N → ∞. The zeroes of the partition function tend to distribute on a slightly distorted unit circle in complex plane and approach the positive real axis already for a 5×5 square lattice. When the lattice has negative cooperativity, its properties mimic those of a system composed of two classes of independent sites with the apparent population of low-affinity binding sites increasing with the size of the lattice, thereby accounting for a phenomenon encountered in many ligand-receptor interactions.

Evolutionary analyses of hedgehog and Hoxd-10 genes in fish species closely related to the zebrafish

The study of development has relied primarily on the isolation of mutations in genes with specific functions in development and on the comparison of their expression patterns in normal and mutant phenotypes. Comparative evolutionary analyses can complement these approaches. Phylogenetic analyses of Sonic hedgehog (Shh) and Hoxd-10 genes from 18 cyprinid fish species closely related to the zebrafish provide novel insights into the functional constraints acting on Shh. Our results confirm and extend those gained from expression and crystalline structure analyses of this gene. Unexpectedly, exon 1 of Shh is found to be almost invariant even in third codon positions among these morphologically divergent species suggesting that this exon encodes for a functionally important domain of the hedgehog protein. This is surprising because the main functional domain of Shh had been thought to be that encoded by exon 2. Comparisons of Shh and Hoxd-10 gene sequences and of resulting gene trees document higher evolutionary constraints on the former than on the latter. This might be indicative of more general evolutionary patterns in networks of developmental regulatory genes interacting in a hierarchical fashion. The presence of four members of the hedgehog gene family in cyprinid fishes was documented and their homologies to known hedgehog genes in other vertebrates were established.

A mammalian mitochondrial drug receptor functions as a bacterial “oxygen” sensor

The rat mitochondrial outer membrane-localized benzodiazepine receptor (MBR) was expressed in wild-type and TspO− (tryptophan-rich sensory protein) strains of the facultative photoheterotroph, Rhodobacter sphaeroides 2.4.1, and was shown to retain its structure within the bacterial outer membrane as assayed by its binding properties with a variety of MBR ligands. Functionally, it was able to substitute for TspO by negatively regulating the expression of photosynthesis genes in response to oxygen. This effect was reversed pharmacologically with the MBR ligand PK11195. These results suggest a close evolutionary and functional relationship between the bacterial TspO and the MBR. This relationship provides further support for the origin of the mammalian mitochondrion from a “photosynthetic” precursor. Finally, these findings provide novel insights into the physiological role that has been obscure for the MBR in situ.

Permanence or change? The meaning of genetic variation

Selected aspects of the evolutionary process and more specifically of the genetic variation are considered, with an emphasis in studies performed by my group. One key aspect of evolution seems to be the concomitant occurrence of dichotomic, contradictory (dialect) processes. Genetic variation is structured, and the dynamics of change at one level is not necessarily paralleled by that in another. The pathogenesis-related protein superfamily can be cited as an example in which permanence (the maintenance of certain key genetic features) coexists with change (modifications that led to different functions in different classes of organisms). Relationships between structure and function are exemplified by studies with hemoglobin Porto Alegre. The genetic structure of tribal populations may differ in important aspects from that of industrialized societies. Evolutionary histories also may differ when considered through the investigation of patrilineal or matrilineal lineages. Global evaluations taking into consideration all of these aspects are needed if we really want to understand the meaning of genetic variation.

Computation, prediction, and experimental tests of fitness for bacteriophage T7 mutants with permuted genomes

We created a simulation based on experimental data from bacteriophage T7 that computes the developmental cycle of the wild-type phage and also of mutants that have an altered genome order. We used the simulation to compute the fitness of more than 105 mutants. We tested these computations by constructing and experimentally characterizing T7 mutants in which we repositioned gene 1, coding for T7 RNA polymerase. Computed protein synthesis rates for ectopic gene 1 strains were in moderate agreement with observed rates. Computed phage-doubling rates were close to observations for two of four strains, but significantly overestimated those of the other two. Computations indicate that the genome organization of wild-type T7 is nearly optimal for growth: only 2.8% of random genome permutations were computed to grow faster, the highest 31% faster, than wild type. Specific discrepancies between computations and observations suggest that a better understanding of the translation efficiency of individual mRNAs and the functions of qualitatively “nonessential” genes will be needed to improve the T7 simulation. In silico representations of biological systems can serve to assess and advance our understanding of the underlying biology. Iteration between computation, prediction, and observation should increase the rate at which biological hypotheses are formulated and tested.

Evolutionary conservation of apoptosis mechanisms: Lepidopteran and baculoviral inhibitor of apoptosis proteins are inhibitors of mammalian caspase-9

We cloned a new inhibitor of apoptosis protein (IAP) homolog, SfIAP, from Spodoptera frugiperda Sf-21 cells, a host of insect baculoviruses. SfIAP contains two baculovirus IAP repeat domains followed by a RING domain. SfIAP has striking amino acid sequence similarity with baculoviral IAPs, CpIAP and OpIAP, suggesting that baculoviral IAPs may be host-derived genes. SfIAP and baculoviral CpIAP inhibit Bax but not Fas-induced apoptosis in human cells. Their apoptosis-suppressing activity in mammalian cells requires both baculovirus IAP repeat and RING domains. Further biochemical data suggest that SfIAP and CpIAP are specific inhibitors of mammalian caspase-9, the pinnacle caspase in the mitochondria/cytochrome c pathway for apoptosis, but are not inhibitors of downstream caspase-3 and caspase-7. Thus the mechanisms by which insect and baculoviral IAPs suppress apoptosis may involve inhibition of an insect caspase-9 homologue. Peptides representing the IAP-binding domain of the Drosophila cell death protein Grim abrogated human caspase suppression by SfIAP and CpIAP, implying evolutionary conservation of the functions of IAPs and their inhibitors.

Cryptocyanin, a crustacean molting protein: Evolutionary link with arthropod hemocyanins and insect hexamerins

Cryptocyanin, a copper-free hexameric protein in crab (Cancer magister) hemolymph, has been characterized and the amino acid sequence has been deduced from its cDNA. It is markedly similar in sequence, size, and structure to hemocyanin, the copper-containing oxygen-transport protein found in many arthropods. Cryptocyanin does not bind oxygen, however, and lacks three of the six highly conserved copper-binding histidine residues of hemocyanin. Cryptocyanin has no phenoloxidase activity, although a phenoloxidase is present in the hemolymph. The concentration of cryptocyanin in the hemolymph is closely coordinated with the molt cycle and reaches levels higher than hemocyanin during premolt. Cryptocyanin resembles insect hexamerins in the lack of copper, molt cycle patterns of biosynthesis, and potential contributions to the new exoskeleton. Phylogenetic analysis of sequence similarities between cryptocyanin and other members of the hemocyanin gene family shows that cryptocyanin is closely associated with crustacean hemocyanins and suggests that cryptocyanin arose as a result of a hemocyanin gene duplication. The presence of both hemocyanin and cryptocyanin in one animal provides an example of how insect hexamerins might have evolved from hemocyanin. Our results suggest that multiple members of the hemocyanin gene family—hemocyanin, cryptocyanin, phenoloxidase, and hexamerins—may participate in two vital functions of molting animals, oxygen binding and molting. Cryptocyanin may provide important molecular data to further investigate evolutionary relationships among all molting animals.

Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships

Pairwise sequence comparison methods have been assessed using proteins whose relationships are known reliably from their structures and functions, as described in the scop database [Murzin, A. G., Brenner, S. E., Hubbard, T. & Chothia C. (1995) J. Mol. Biol. 247, 536–540]. The evaluation tested the programs blast [Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. 215, 403–410], wu-blast2 [Altschul, S. F. & Gish, W. (1996) Methods Enzymol. 266, 460–480], fasta [Pearson, W. R. & Lipman, D. J. (1988) Proc. Natl. Acad. Sci. USA 85, 2444–2448], and ssearch [Smith, T. F. & Waterman, M. S. (1981) J. Mol. Biol. 147, 195–197] and their scoring schemes. The error rate of all algorithms is greatly reduced by using statistical scores to evaluate matches rather than percentage identity or raw scores. The E-value statistical scores of ssearch and fasta are reliable: the number of false positives found in our tests agrees well with the scores reported. However, the P-values reported by blast and wu-blast2 exaggerate significance by orders of magnitude. ssearch, fasta ktup = 1, and wu-blast2 perform best, and they are capable of detecting almost all relationships between proteins whose sequence identities are >30%. For more distantly related proteins, they do much less well; only one-half of the relationships between proteins with 20–30% identity are found. Because many homologs have low sequence similarity, most distant relationships cannot be detected by any pairwise comparison method; however, those which are identified may be used with confidence.

Genomic evidence for two functionally distinct gene classes

Analyses of complete genomes indicate that a massive prokaryotic gene transfer (or transfers) preceded the formation of the eukaryotic cell. In comparisons of the entire set of Methanococcus jannaschii genes with their orthologs from Escherichia coli, Synechocystis 6803, and the yeast Saccharomyces cerevisiae, it is shown that prokaryotic genomes consist of two different groups of genes. The deeper, diverging informational lineage codes for genes which function in translation, transcription, and replication, and also includes GTPases, vacuolar ATPase homologs, and most tRNA synthetases. The more recently diverging operational lineage codes for amino acid synthesis, the biosynthesis of cofactors, the cell envelope, energy metabolism, intermediary metabolism, fatty acid and phospholipid biosynthesis, nucleotide biosynthesis, and regulatory functions. In eukaryotes, the informational genes are most closely related to those of Methanococcus, whereas the majority of operational genes are most closely related to those of Escherichia, but some are closest to Methanococcus or to Synechocystis.

Unusual horizontal transfer of a long interspersed nuclear element between distant vertebrate classes

We have shown previously by Southern blot analysis that Bov-B long interspersed nuclear elements (LINEs) are present in different Viperidae snake species. To address the question as to whether Bov-B LINEs really have been transmitted horizontally between vertebrate classes, the analysis has been extended to a larger number of vertebrate, invertebrate, and plant species. In this paper, the evolutionary origin of Bov-B LINEs is shown unequivocally to be in Squamata. The previously proposed horizontal transfer of Bov-B LINEs in vertebrates has been confirmed by their discontinuous phylogenetic distribution in Squamata (Serpentes and two lizard infra-orders) as well as in Ruminantia, by the high level of nucleotide identity, and by their phylogenetic relationships. The horizontal transfer of Bov-B LINEs from Squamata to the ancestor of Ruminantia is evident from the genetic distances and discontinuous phylogenetic distribution. The ancestor of Colubroidea snakes is a possible donor of Bov-B LINEs to Ruminantia. The timing of horizontal transfer has been estimated from the distribution of Bov-B LINEs in Ruminantia and the fossil data of Ruminantia to be 40–50 My ago. The phylogenetic relationships of Bov-B LINEs from the various Squamata species agrees with that of the species phylogeny, suggesting that Bov-B LINEs have been maintained stably by vertical transmission since the origin of Squamata in the Mesozoic era.

Regulation of number and size of digits by posterior Hox genes: A dose-dependent mechanism with potential evolutionary implications

The proper development of digits, in tetrapods, requires the activity of several genes of the HoxA and HoxD homeobox gene complexes. By using a variety of loss-of-function alleles involving the five Hox genes that have been described to affect digit patterning, we report here that the group 11, 12, and 13 genes control both the size and number of murine digits in a dose-dependent fashion, rather than through a Hox code involving differential qualitative functions. A similar dose–response is observed in the morphogenesis of the penian bone, the baculum, which further suggests that digits and external genitalia share this genetic control mechanism. A progressive reduction in the dose of Hox gene products led first to ectrodactyly, then to olygodactyly and adactyly. Interestingly, this transition between the pentadactyl to the adactyl formula went through a step of polydactyly. We propose that in the distal appendage of polydactylous short-digited ancestral tetrapods, such as Acanthostega, the HoxA complex was predominantly active. Subsequent recruitment of the HoxD complex contributed to both reductions in digit number and increase in digit length. Thus, transition through a polydactylous limb before reaching and stabilizing the pentadactyl pattern may have relied, at least in part, on asynchronous and independent changes in the regulation of HoxA and HoxD gene complexes.

Requirement of the Lec35 Gene for All Known Classes of Monosaccharide-P-Dolichol-dependent Glycosyltransferase Reactions in Mammals

The Lec35 gene product (Lec35p) is required for utilization of the mannose donor mannose-P-dolichol (MPD) in synthesis of both lipid-linked oligosaccharides (LLOs) and glycosylphosphatidylinositols, which are important for functions such as protein folding and membrane anchoring, respectively. The hamster Lec35 gene is shown to encode the previously identified cDNA SL15, which corrects the Lec35 mutant phenotype and predicts a novel endoplasmic reticulum membrane protein. The mutant hamster alleles Lec35.1 and Lec35.2 are characterized, and the human Lec35 gene (mannose-P-dolichol utilization defect 1) was mapped to 17p12-13. To determine whether Lec35p was required only for MPD-dependent mannosylation of LLO and glycosylphosphatidylinositol intermediates, two additional lipid-mediated reactions were investigated: MPD-dependent C-mannosylation of tryptophanyl residues, and glucose-P-dolichol (GPD)-dependent glucosylation of LLO. Both were found to require Lec35p. In addition, the SL15-encoded protein was selective for MPD compared with GPD, suggesting that an additional GPD-selective Lec35 gene product remains to be identified. The predicted amino acid sequence of Lec35p does not suggest an obvious function or mechanism. By testing the water-soluble MPD analog mannose-β-1-P-citronellol in an in vitro system in which the MPD utilization defect was preserved by permeabilization with streptolysin-O, it was determined that Lec35p is not directly required for the enzymatic transfer of mannose from the donor to the acceptor substrate. These results show that Lec35p has an essential role for all known classes of monosaccharide-P-dolichol-dependent reactions in mammals. The in vitro data suggest that Lec35p controls an aspect of MPD orientation in the endoplasmic reticulum membrane that is crucial for its activity as a donor substrate.

Evolutionary relationships among Rel domains indicate functional diversification by recombination

The recent sequencing of several complete genomes has made it possible to track the evolution of large gene families by their genomic structure. Following the large-scale association of exons encoding domains with well defined functions in invertebrates could be useful in predicting the function of complex multidomain proteins in mammals produced by accretion of domains. With this objective, we have determined the genomic structure of the 14 genes in invertebrates and vertebrates that contain rel domains. The sequence encoding the rel domain is defined by intronic boundaries and has been recombined with at least three structurally and functionally distinct genomic sequences to generate coding sequences for: (i) the rel/Dorsal/NFκB proteins that are retained in the cytoplasm by IkB-like proteins; (ii) the NFATc proteins that sense calcium signals and undergo cytoplasmic-to-nuclear translocation in response to dephosphorylation by calcineurin; and (iii) the TonEBP tonicity-responsive proteins. Remarkably, a single exon in each NFATc family member encodes the entire Ca2+/calcineurin sensing region, including nuclear import/export, calcineurin-binding, and substrate regions. The Rel/Dorsal proteins and the TonEBP proteins are present in Drosophila but not Caenorhabditis elegans. On the other hand, the calcium-responsive NFATc proteins are present only in vertebrates, suggesting that the NFATc family is dedicated to functions specific to vertebrates such as a recombinational immune response, cardiovascular development, and vertebrate-specific aspects of the development and function of the nervous system.

From symmetry to asymmetry: Phylogenetic patterns of asymmetry variation in animals and their evolutionary significance

Phylogenetic analyses of asymmetry variation offer a powerful tool for exploring the interplay between ontogeny and evolution because (i) conspicuous asymmetries exist in many higher metazoans with widely varying modes of development, (ii) patterns of bilateral variation within species may identify genetically and environmentally triggered asymmetries, and (iii) asymmetries arising at different times during development may be more sensitive to internal cytoplasmic inhomogeneities compared to external environmental stimuli. Using four broadly comparable asymmetry states (symmetry, antisymmetry, dextral, and sinistral), and two stages at which asymmetry appears developmentally (larval and postlarval), I evaluated relations between ontogenetic and phylogenetic patterns of asymmetry variation. Among 140 inferred phylogenetic transitions between asymmetry states, recorded from 11 classes in five phyla, directional asymmetry (dextral or sinistral) evolved directly from symmetrical ancestors proportionally more frequently among larval asymmetries. In contrast, antisymmetry, either as an end state or as a transitional stage preceding directional asymmetry, was confined primarily to postlarval asymmetries. The ontogenetic origin of asymmetry thus significantly influences its subsequent evolution. Furthermore, because antisymmetry typically signals an environmentally triggered asymmetry, the phylogenetic transition from antisymmetry to directional asymmetry suggests that many cases of laterally fixed asymmetries evolved via genetic assimilation.