Constraints on protein evolution and the age of the eubacteria/eukaryote split

This article discusses research done by Doolittle et al. to determine the evolutionary distance between eubacteria and eukaryotes. Several equations for estimating evolutionary distances and divergence times are considered. Molecular clock calculations and their effects on these estimates are discussed. The authors conclude that their research supports the range of dates for the cenancestor of modern eubacteria, archaebacteria, and eukaryotes estimated by Doolittle et al.     

Contextual constraints on synonymous codon choice

We have studied the statistical constraints on synonymous codon choice to evaluate various proposals regarding the origin of the bias in synonymous codon usage observed by Fiers et al. (1975), Air et al. (1976), Grantham et al. (1980) and others. We have determined the statistical dependence of the degenerate third base on either of its nearest neighbors in mitochondrial, prokaryotic, and eukaryotic coding sequences. We noted an increasing dependence of the third base on its nearest neighbors in moving from mitochrondria to prokaryotes to eukaryotes.A statistical model assuming random equiprobable selection of synonymous codons was found grossly adequate for the mitochondria, but totally indequate for prokaryotes and eukaryotes. A model assuming selection of synonymous codons reflecting a genomic strategy, i.e. the genome hypothesis of Grantham et al. (1980), gave a good approximation of the mitochondrial sequences. A statistical model which exactly maintains codon frequency, but allows the position of corresponding synonymous codons to vary was only grossly adequate for prokaryotes and totally inadequate for eukaryotes. The results of these simulations are consistent with the measures on experimental sequences and suggest that a “frequency constraint” model such as that of Grantham et al. (1980) may be an adequate explanation of the codon usage in mitochondria. However, in addition to this frequency constraint, there may be constraints on synonymous codon choice in prokaryotes due to codon context. Furthermore, any proposal to explain codon usage in eukaryotes must involve a constraint on the context of a codon in the sequence.     

Horizontal gene transfer in microbial genome evolution

Horizontal gene transfer is the collective name for processes that permit the exchange of DNA among organisms of different species. Only recently has it been recognized as a significant contribution to inter-organismal gene exchange. Traditionally, it was thought that microorganisms evolved clonally, passing genes from mother to daughter cells with little or no exchange of DNA among diverse species. Studies of microbial genomes, however, have shown that genomes contain genes that are closely related to a number of different prokaryotes, sometimes to phylogenetically very distantly related ones. (Doolittle et al., 1990, J. Mol. Evol. 31, 383-388; Karlin et al., 1997, J. Bacteriol. 179, 3899-3913; Karlin et al., 1998, Annu. Rev. Genet. 32, 185-225; Lawrence and Ochman, 1998, Proc. Natl. Acad. Sci. USA 95, 9413-9417; Rivera et al., 1998, Proc. Natl. Acad. Sci. USA 95, 6239-6244; Campbell, 2000, Theor. Popul. Biol. 57 71-77; Doolittle, 2000, Sci. Am. 282, 90-95; Ochman and Jones, 2000, Embo. J. 19, 6637-6643; Boucher et al. 2001, Curr. Opin., Microbiol. 4, 285-289; Wang et al., 2001, Mol. Biol. Evol. 18, 792-800). Whereas prokaryotic and eukaryotic evolution was once reconstructed from a single 16S ribosomal RNA (rRNA) gene, the analysis of complete genomes is beginning to yield a different picture of microbial evolution, one that is wrought with the lateral movement of genes across vast phylogenetic distances. (Lane et al., 1988, Methods Enzymol. 167, 138-144; Lake and Rivera, 1996, Proc. Natl. Acad. Sci. USA 91, 2880-2881; Lake et al., 1999, Science 283, 2027-2028).     

The P450 superfamily: update on new sequences, gene mapping, and recommended nomenclature

We provide here a list of 154 P450 genes and seven putative pseudogenes that have been characterized as of October 20, 1990. These genes have been described in a total of 23 eukaryotes (including nine mammalian and one plant species) and six prokaryotes. Of 27 gene families so far described, 10 exist in all mammals. These 10 families comprise 18 subfamilies, of which 16 and 14 have been mapped in the human and mouse genomes, respectively; to date, each subfamily appears to represent a cluster of tightly linked genes. We propose here a modest revision of the initially proposed (Nebert et al., DNA 6, 1-11, 1987) and updated (Nebert et al., DNA 8, 1-13, 1989) nomenclature system based on evolution of the superfamily. For the gene we recommend that the italicized root symbol CYP for human (Cyp for mouse), representing cytochrome P450, be followed by an Arabic number denoting the family, a letter designating the subfamily (when two or more exist), and an Arabic numeral representing the individual gene within the subfamily. A hyphen should precede the final number in mouse genes. We suggest that the human nomenclature system be used for other species. This system is consistent with our earlier proposed nomenclature for P450 of all eukaryotes and prokaryotes, except that we are discouraging the future use of cumbersome Roman numerals.     

"Anticipated" nucleosome positioning pattern in prokaryotes

Linguistic (word count) analysis of prokaryotic genome sequences, by Shannon N-gram extension, reveals that the dominant hidden motifs in A+T rich genomes are T(A)(T)A and G(A)(T)C with uncertain number of repeating A and T. Since prokaryotic sequences are largely protein-coding, the motifs would correspond to amphipathic alpha-helices with alternating lysine and phenylalanine as preferential polar and non-polar residues. The motifs are also known in eukaryotes, as nucleosome positioning patterns. Their existence in prokaryotes as well may serve for binding of histone-like proteins to DNA. In this case the above patterns in prokaryotes may be considered as "anticipated" nucleosome positioning patterns which, quite likely, existed in prokaryotic genomes before the evolutionary separation between eukaryotes and prokaryotes.     

Genomic GC level, optimal growth temperature, and genome size in prokaryotes

Two years ago, we showed that positive correlations between optimal growth temperature (T(opt)) and genome GC are observed in 15 out of the 20 families of prokaryotes we analyzed, thus indicating that "T(opt) is one of the factors that influence genomic GC in prokaryotes". Our results were disputed, but these criticisms were demonstrated to be mistaken and based on misconceptions. In a recent report, Wang et al. [H.C. Wang, E. Susko, A.J. Roger, On the correlation between genomic G+C content and optimal growth temperature in prokaryotes: data quality and confounding factors, Biochem. Biophys. Res. Commun. 342 (2006) 681-684] criticize our results by stating that "all previous simple correlation analyses of GC versus temperature have ignored the fact that genomic GC content is influenced by multiple factors including both intrinsic mutational bias and extrinsic environmental factors". This statement, besides being erroneous, is surprising because it applies in fact not to ours but to the authors' article. Here, we rebut the points raised by Wang et al. and review some issues that have been a matter of debate, regarding the influence of environmental factors upon GC content in prokaryotes. Furthermore, we demonstrate that the relationship that exists between genome size and GC level is valid for aerobic, facultative, and microaerophilic species, but not for anaerobic prokaryotes.     

Symbiont genes in host genomes: fragments with a future?

While lateral transfer is the rule in the evolutionary history of bacterial and archaeal genes, events of transfer from prokaryotes to eukaryotes are rare. Germline-transmitted animal symbionts, such as Wolbachia pipientis, are well placed to participate in such transfers. In a recent issue of Science, Dunning Hotopp et al. identified instances of transfer of Wolbachia DNA to host genomes. It is unknown whether these transfers represent innovation in animal evolution.     

Usage of the three termination codons: Compilation and analysis of the known eukaryotic and prokaryotic translation termination sequences

Summary The published translation termination sequences have been compiled and analysed to aid the interpretation of experiments on termination codon usage in the Xenopus oocyte (Bienz et al. 1981). There are significant differences between prokaryotes and eukaryotes concerning the usage of the three termination codons and of tandem stops. In addition viruses show termination strategies that differ from those of their hosts. Preferred context sequences flanking termination codons are described. Contexts vary within the last codon according to the nature of the termination codon, but are uniform within the first triplet following the terminators.     

Phylogeny of Prokaryotes and Chloroplasts Revealed by a Simple Composition Approach on All Protein Sequences from Complete Genomes Without Sequence Alignment

The complete genomes of living organisms have provided much information on their phylogenetic relationships. Similarly, the complete genomes of chloroplasts have helped to resolve the evolution of this organelle in photosynthetic eukaryotes. In this paper we propose an alternative method of phylogenetic analysis using compositional statistics for all protein sequences from complete genomes. This new method is conceptually simpler than and computationally as fast as the one proposed by Qi et al. (2004b) and Chu et al. (2004). The same data sets used in Qi et al. (2004b) and Chu et al. (2004) are analyzed using the new method. Our distance-based phylogenic tree of the 109 prokaryotes and eukaryotes agrees with the biologists tree of life based on 16S rRNA comparison in a predominant majority of basic branching and most lower taxa. Our phylogenetic analysis also shows that the chloroplast genomes are separated to two major clades corresponding to chlorophytes s.l. and rhodophytes s.l. The interrelationships among the chloroplasts are largely in agreement with the current understanding on chloroplast evolution.Reviewing Editor: Dr. John Oakeshott     

Comparing function and structure between entire proteomes

More than 30 organisms have been sequenced entirely. Here, we applied a variety of simple bioinformatics tools to analyze 29 proteomes for representatives from all three kingdoms: eukaryotes, prokaryotes, and archaebacteria. We confirmed that eukaryotes have relatively more long proteins than prokaryotes and archaes, and that the overall amino acid composition is similar among the three. We predicted that approximately 15%-30% of all proteins contained transmembrane helices. We could not find a correlation between the content of membrane proteins and the complexity of the organism. In particular, we did not find significantly higher percentages of helical membrane proteins in eukaryotes than in prokaryotes or archae. However, we found more proteins with seven transmembrane helices in eukaryotes and more with six and 12 transmembrane helices in prokaryotes. We found twice as many coiled-coil proteins in eukaryotes (10%) as in prokaryotes and archaes (4%-5%), and we predicted approximately 15%-25% of all proteins to be secreted by most eukaryotes and prokaryotes. Every tenth protein had no known homolog in current databases, and 30%-40% of the proteins fell into structural families with >100 members. A classification by cellular function verified that eukaryotes have a higher proportion of proteins for communication with the environment. Finally, we found at least one homolog of experimentally known structure for approximately 20%-45% of all proteins; the regions with structural homology covered 20%-30% of all residues. These numbers may or may not suggest that there are 1200-2600 folds in the universe of protein structures. All predictions are available at http://cubic.bioc.columbia.edu/genomes.     

Comparison of dissimilarity patterns of E coli, yeast and mammalian tRNAs.

A number of experimental approaches have been developed for identification of recognition (identity) sites in tRNAs. Along with them a theoretical methodology has been proposed by McClain et al that is based on concomitant analysis of all tRNA sequences from a given species. This approach allows an evaluation of nucleotide combinations present in isoacceptor tRNAs specific for the given amino acid, and not present in equivalent positions in cloverleaf structure in other tRNAs of the same organism. These elements predicted from computer analysis of the databank could be tested experimentally for their participation in forming recognition sites. The correlation between theoretical predictions and experimental data appeared promising. The aim of the present work consisted of introducing further improvements into McClain's procedure by: i), introducing into analysis a variable region in tRNAs which had not been previously considered; to accomplish this, 'normalization' of variable nucleotides was suggested, based on primary and tertiary structures of tRNAs; ii), developing a new procedure for comparison of patterns for synonymous and non-synonymous tRNAs from different organisms; iii), analysis of 3- and 4-positional contacts between tRNAs and enzymes in addition to a formerly used 2-positional model. A systematic application of McClain's procedure to mammalian, yeast and E coli tRNAs led to the following results: i), imitancy patterns for non-synonymous tRNAs of any amino acid specificity and from any organisms analysed so far overlap by no more than 30%, providing a structural basis for discrimination with high fidelity between cognate and non-cognate tRNAs; ii), the predicted identity sites are non-randomly distributed within tRNA molecules; the dominant role is ascribed to only two regions--anticodon and amino acid stem which are located far apart from one another at extremes of all tRNA molecules; iii), the imitancy patterns for synonymous tRNAs in lower (yeast) and higher (mammalian) eukaryotes are similar but not identical; iv), distribution of predicted identity sites in the cloverleaf structure in prokaryotes and eukaryotes is essentially different: in eubacterial tRNAs the major role in recognition plays anticodon and/or amino acid acceptor stem, whereas in eukaryotic (both unicellular and multicellular) tRNAs the remaining part of the molecules is also involved in recognition; v), the imitancy patterns of synonymous tRNAs from prokaryotes and eukaryotes are dissimilar, this observation leads to the prediction that the tRNA identity sites for the same amino acid in prokaryotes and eukaryotes may differ.     

Energy metabolism among eukaryotic anaerobes in light of Proterozoic ocean chemistry

Recent years have witnessed major upheavals in views about early eukaryotic evolution. One very significant finding was that mitochondria, including hydrogenosomes and the newly discovered mitosomes, are just as ubiquitous and defining among eukaryotes as the nucleus itself. A second important advance concerns the readjustment, still in progress, about phylogenetic relationships among eukaryotic groups and the roughly six new eukaryotic supergroups that are currently at the focus of much attention. From the standpoint of energy metabolism (the biochemical means through which eukaryotes gain their ATP, thereby enabling any and all evolution of other traits), understanding of mitochondria among eukaryotic anaerobes has improved. The mainstream formulations of endosymbiotic theory did not predict the ubiquity of mitochondria among anaerobic eukaryotes, while an alternative hypothesis that specifically addressed the evolutionary origin of energy metabolism among eukaryotic anaerobes did. Those developments in biology have been paralleled by a similar upheaval in the Earth sciences regarding views about the prevalence of oxygen in the oceans during the Proterozoic (the time from ca 2.5 to 0.6 Ga ago). The new model of Proterozoic ocean chemistry indicates that the oceans were anoxic and sulphidic during most of the Proterozoic. Its proponents suggest the underlying geochemical mechanism to entail the weathering of continental sulphides by atmospheric oxygen to sulphate, which was carried into the oceans as sulphate, fueling marine sulphate reducers (anaerobic, hydrogen sulphide-producing prokaryotes) on a global scale. Taken together, these two mutually compatible developments in biology and geology underscore the evolutionary significance of oxygen-independent ATP-generating pathways in mitochondria, including those of various metazoan groups, as a watermark of the environments within which eukaryotes arose and diversified into their major lineages.     

Insulin-related materials in the nervous system of vertebrates and non-vertebrates: possible extrapancreatic production

Studies from multiple laboratories with a range of methods raised the possibility that insulin production occurs naturally at extrapancreatic sites. Part A covers the presence of insulin-related materials in organisms that do not have an endocrine pancreas, including unicellular prokaryotes and eukaryotes as well as multicellular non-vertebrate animals (insects et al.) and plants. Part B covers possible production of insulin by extrapancreatic tissues of vertebrates that are remote from a source of pancreatic insulin e.g. early chick embryos and mammalian cells in culture. Part C covers possible extrapancreatic insulin production in mammals in vivo. Each section ends with an outline summary with evidence in favor of and against the hypothesis.     

The fidelity of translation initiation: reciprocal activities of eIF1, IF3 and YciH

Eukaryotic initiation factor eIF1 and the functional C-terminal domain of prokaryotic initiation factor IF3 maintain the fidelity of initiation codon selection in eukaryotes and prokaryotes, respectively, and bind to the same regions of small ribosomal subunits, between the platform and initiator tRNA. Here we report that these nonhomologous factors can bind to the same regions of heterologous subunits and perform their functions in heterologous systems in a reciprocal manner, discriminating against the formation of initiation complexes containing codon-anticodon mismatches. We also show that like IF3, eIF1 can influence initiator tRNA selection, which occurs at the stage of ribosomal subunit joining after eIF5-induced hydrolysis of eIF2-bound GTP. The mechanisms of initiation codon and initiator tRNA selection in prokaryotes and eukaryotes are therefore unexpectedly conserved and likely involve related conformational changes induced in the small ribosomal subunit by factor binding. YciH, a prokaryotic eIF1 homologue, could perform some of IF3's functions, which justifies the possibility that YciH and eIF1 might have a common evolutionary origin as initiation factors, and that IF3 functionally replaced YciH in prokaryotes.     

Seasonal variation in Plasmodium prevalence in a population of blue tits Cyanistes caeruleus

1. Seasonal variation in environmental conditions is ubiquitous and can affect the spread of infectious diseases. Understanding seasonal patterns of disease incidence can help to identify mechanisms, such as the demography of hosts and vectors, which influence parasite transmission dynamics. 2. We examined seasonal variation in Plasmodium infection in a blue tit Cyanistes caeruleus population over 3 years using sensitive molecular diagnostic techniques, in light of Beaudoin et al.'s (1971; Journal of Wildlife Diseases, 7, 5-13) model of seasonal variation in avian malaria prevalence in temperate areas. This model predicts a within-year bimodal pattern of spring and autumn peaks with a winter absence of infection. 3. Avian malaria infections were mostly Plasmodium (24.4%) with occasional Haemoproteus infections (0.8%). Statistical nonlinear smoothing techniques applied to longitudinal presence/absence data revealed marked temporal variation in Plasmodium prevalence, which apparently showed a within-year bimodal pattern similar to Beaudoin et al.'s model. However, of the two Plasmodium morphospecies accounting for most infections, only the seasonal pattern of Plasmodium circumflexum supported Beaudoin et al.'s model. On closer examination there was also considerable age structure in infection: Beaudoin et al.'s seasonal pattern was observed only in first year and not older birds. Plasmodium relictum prevalence was less seasonally variable. 4. For these two Plasmodium morphospecies, we reject Beaudoin et al.'s model as it does not survive closer scrutiny of the complexities of seasonal variation among Plasmodium morphospecies and host age classes. Studies of host-parasite interactions should consider seasonal variation whenever possible. We discuss the ecological and evolutionary implications of seasonal variation in disease prevalence.     

The Drosophila melanogaster RPS17 gene encoding ribosomal protein S17

A human ribosomal protein S17 cDNA [Chen et al., Proc. Natl. Acad. Sci. USA 83 (1986) 6907-6911] was used as heterologous probe to isolate S17 clones from Drosophila genomic and cDNA recombinant libraries. Five S17 genomic clones were recognized; all contained overlapping regions of a single chromosomal site. Subsequently the Drosophila RPS17 gene was mapped by in situ hybridization to chromosome 3L, band 67B1-5. The locus spans approximately 1000 bp of DNA and includes four exons. It is preceded by conventional CAAT and TATA RNA polymerase II promoter motifs. The 131 amino acid protein encoded within Drosophila RPS17 is similar to ribosomal proteins from several other eukaryotes. Comparison of eukaryotic S17 proteins' primary structures as well as the number and location of their genes' intervening sequences suggest that S17 is a relatively recent addition to the ribosomal protein family, probably post-dating divergence of eukaryotes and prokaryotes.     

Old genes in new places: A taxon-rich analysis of interdomain lateral gene transfer events

Vertical inheritance is foundational to Darwinian evolution, but fails to explain major innovations such as the rapid spread of antibiotic resistance among bacteria and the origin of photosynthesis in eukaryotes. While lateral gene transfer (LGT) is recognized as an evolutionary force in prokaryotes, the role of LGT in eukaryotic evolution is less clear. With the exception of the transfer of genes from organelles to the nucleus, a process termed endosymbiotic gene transfer (EGT), the extent of interdomain transfer from prokaryotes to eukaryotes is highly debated. A common critique of studies of interdomain LGT is the reliance on the topology of single-gene trees that attempt to estimate more than one billion years of evolution. We take a more conservative approach by identifying cases in which a single clade of eukaryotes is found in an otherwise prokaryotic gene tree (i.e. exclusive presence). Starting with a taxon-rich dataset of over 13,600 gene families and passing data through several rounds of curation, we identify and categorize the function of 306 interdomain LGT events into diverse eukaryotes, including 189 putative EGTs, 52 LGTs into Opisthokonta (i.e. animals, fungi and their microbial relatives), and 42 LGTs nearly exclusive to anaerobic eukaryotes. To assess differential gene loss as an explanation for exclusive presence, we compare branch lengths within each LGT tree to a set of vertically-inherited genes subsampled to mimic gene loss (i.e. with the same taxonomic sampling) and consistently find shorter relative distance between eukaryotes and prokaryotes in LGT trees, a pattern inconsistent with gene loss. Our methods provide a framework for future studies of interdomain LGT and move the field closer to an understanding of how best to model the evolutionary history of eukaryotes.