Thermodynamic selection: mechanisms and scenarios.

Thermodynamic selection is an indirect competition between agents feeding on the same energy resource and obeying the laws of thermodynamics. We examine scenarios of this selection, where the agent is modeled as a heat-engine coupled to two thermal baths and extracting work from the high-temperature bath. The agents can apply different work-extracting, game-theoretical strategies, e.g. the maximum power or the maximum efficiency. They can also have a fixed structure or be adaptive. Depending on whether the resource (i.e. the high-temperature bath) is infinite or finite, the fitness of the agent relates to the work-power or the total extracted work. These two selection scenarios lead to increasing or decreasing efficiencies of the work-extraction, respectively. The scenarios are illustrated via plant competition for sunlight, and the competition between different ATP production pathways. We also show that certain general concepts of game-theory and ecology-the prisoner's dilemma and the maximal power principle-emerge from the thermodynamics of competing agents. We emphasize the role of adaptation in developing efficient work-extraction mechanisms.

Ancient Systems of Sodium/Potassium Homeostasis as Predecessors of Membrane Bioenergetics.

Cell cytoplasm of archaea, bacteria, and eukaryotes contains substantially more potassium than sodium, and potassium cations are specifically required for many key cellular processes, including protein synthesis. This distinct ionic composition and requirements have been attributed to the emergence of the first cells in potassium-rich habitats. Different, albeit complementary, scenarios have been proposed for the primordial potassium-rich environments based on experimental data and theoretical considerations. Specifically, building on the observation that potassium prevails over sodium in the vapor of inland geothermal systems, we have argued that the first cells could emerge in the pools and puddles at the periphery of primordial anoxic geothermal fields, where the elementary composition of the condensed vapor would resemble the internal milieu of modern cells. Marine and freshwater environments generally contain more sodium than potassium. Therefore, to invade such environments, while maintaining excess of potassium over sodium in the cytoplasm, primordial cells needed means to extrude sodium ions. The foray into new, sodium-rich habitats was the likely driving force behind the evolution of diverse redox-, light-, chemically-, or osmotically-dependent sodium export pumps and the increase of membrane tightness. Here we present a scenario that details how the interplay between several, initially independent sodium pumps might have triggered the evolution of sodium-dependent membrane bioenergetics, followed by the separate emergence of the proton-dependent bioenergetics in archaea and bacteria. We also discuss the development of systems that utilize the sodium/potassium gradient across the cell membranes.

Predicting functions from protein sequences--where are the bottlenecks?

The exponential growth of sequence data does not necessarily lead to an increase in knowledge about the functions of genes and their products. Prediction of function using comparative sequence analysis is extremely powerful but, if not performed appropriately, may also lead to the creation and propagation of assignment errors. While current homology detection methods can cope with the data flow, the identification, verification and annotation of functional features need to be drastically improved.

Molecular characterization of the plant virus genus Ourmiavirus and evidence of inter-kingdom reassortment of viral genome segments as its possible route of origin.

Ourmia melon virus (OuMV), Epirus cherry virus (EpCV) and Cassava virus C (CsVC) are three species placed in the genus Ourmiavirus. We cloned and sequenced their RNA genomes. The sizes of the three genomic RNAs of OuMV, the type member of the genus, were 2814, 1064 and 974 nt and each had one open reading frame. RNA1 potentially encoded a 97.5 kDa protein carrying the GDD motif typical of RNA-dependent RNA polymerases (RdRps). The putative RdRps of ourmiaviruses are distantly related to known viral RdRps, with the closest similarity and phylogenetic affinity observed with fungal viruses of the genus Narnaviridae. RNA2 encoded a 31.6 kDa protein which, expressed in bacteria as a His-tag fusion protein and in plants through agroinfiltration, reacted specifically with antibodies made against tubular structures found in the cytoplasm. The ORF2 product is significantly similar to movement proteins of the genus Tombusviridae, and phylogenetic analysis supported this evolutionary relationship. The product of OuMV ORF3 is a 23.8 kDa protein. This protein was also expressed in bacteria and plants, and reacted specifically with antisera against the OuMV coat protein. The sequence of the ORF3 protein showed limited but significant similarity to capsid proteins of several plant and animal viruses, although phylogenetic analysis failed to reveal its most likely origin. Taken together, these results indicate that ourmiaviruses comprise a unique group of plant viruses that might have evolved by reassortment of genomic segments of RNA viruses infecting hosts belonging to different eukaryotic kingdoms, in particular, fungi and plants.

Apoptotic molecular machinery: vastly increased complexity in vertebrates revealed by genome comparisons.

A comparison of the proteins encoded in the recently (nearly) completed human genome to those from the fly and nematode genomes reveals a major increase in the complexity of the apoptotic molecular machinery in vertebrates, in terms of both the number of proteins involved and their domain architecture. Several components of the apoptotic system are shared by humans and flies, to the exclusion of nematodes, which seems to support the existence of a coelomate clade in animal evolution. A considerable repertoire of apoptotic protein domains was detected in Actinomycetes and Cyanobacteria, which suggests a major contribution of horizontal gene transfer to the early evolution of apoptosis.

A genomic perspective on protein families.

In order to extract the maximum amount of information from the rapidly accumulating genome sequences, all conserved genes need to be classified according to their homologous relationships. Comparison of proteins encoded in seven complete genomes from five major phylogenetic lineages and elucidation of consistent patterns of sequence similarities allowed the delineation of 720 clusters of orthologous groups (COGs). Each COG consists of individual orthologous proteins or orthologous sets of paralogs from at least three lineages. Orthologs typically have the same function, allowing transfer of functional information from one member to an entire COG. This relation automatically yields a number of functional predictions for poorly characterized genomes. The COGs comprise a framework for functional and evolutionary genome analysis.

Genome sequence of a human tumorigenic poxvirus: prediction of specific host response-evasion genes.

Molluscum contagiosum virus (MCV) commonly causes asymptomatic cutaneous neoplasms in children and sexually active adults as well as persistent opportunistic acquired immunodeficiency syndrome (AIDS)-associated disease. Sequencing the 190-kilobase pair genome of MCV has now revealed that the virus potentially encodes 163 proteins, of which 103 have homologs in the smallpox virus. MCV lacks counterparts to 83 genes of the smallpox virus, including those important in suppression of host responses to infection, nucleotide biosynthesis, and cell proliferation. MCV possesses 59 genes that are predicted to encode previously uncharacterized proteins, including major histocompatibility complex class I, chemokine, and glutathione peroxidase homologs, which suggests that there are MCV-specific strategies for coexistence with the human host.

Pineal serotonin N-acetyltransferase: expression cloning and molecular analysis.

Pineal serotonin N-acetyltransferase (arylalkylamine N-acetyltransferase, or AA-NAT) generates the large circadian rhythm in melatonin, the hormone that coordinates daily and seasonal physiology in some mammals. Complementary DNA encoding ovine AA-NAT was cloned. The abundance of AA-NAT messenger RNA (mRNA) during the day was high in the ovine pineal gland and somewhat lower in retina. AA-NAT mRNA was found unexpectedly in the pituitary gland and in some brain regions. The night-to-day ratio of ovine pineal AA-NAT mRNA is less than 2. In contrast, the ratio exceeds 150 in rats. AA-NAT represents a family within a large superfamily of acetyltransferases.

Comparison of the complete protein sets of worm and yeast: orthology and divergence.

Comparative analysis of predicted protein sequences encoded by the genomes of Caenorhabditis elegans and Saccharomyces cerevisiae suggests that most of the core biological functions are carried out by orthologous proteins (proteins of different species that can be traced back to a common ancestor) that occur in comparable numbers. The specialized processes of signal transduction and regulatory control that are unique to the multicellular worm appear to use novel proteins, many of which re-use conserved domains. Major expansion of the number of some of these domains seen in the worm may have contributed to the advent of multicellularity. The proteins conserved in yeast and worm are likely to have orthologs throughout eukaryotes; in contrast, the proteins unique to the worm may well define metazoans.

Genome sequence of an obligate intracellular pathogen of humans: Chlamydia trachomatis.

Analysis of the 1,042,519-base pair Chlamydia trachomatis genome revealed unexpected features related to the complex biology of chlamydiae. Although chlamydiae lack many biosynthetic capabilities, they retain functions for performing key steps and interconversions of metabolites obtained from their mammalian host cells. Numerous potential virulence-associated proteins also were characterized. Several eukaryotic chromatin-associated domain proteins were identified, suggesting a eukaryotic-like mechanism for chlamydial nucleoid condensation and decondensation. The phylogenetic mosaic of chlamydial genes, including a large number of genes with phylogenetic origins from eukaryotes, implies a complex evolution for adaptation to obligate intracellular parasitism.

Bacterial rhodopsin: evidence for a new type of phototrophy in the sea.

Extremely halophilic archaea contain retinal-binding integral membrane proteins called bacteriorhodopsins that function as light-driven proton pumps. So far, bacteriorhodopsins capable of generating a chemiosmotic membrane potential in response to light have been demonstrated only in halophilic archaea. We describe here a type of rhodopsin derived from bacteria that was discovered through genomic analyses of naturally occuring marine bacterioplankton. The bacterial rhodopsin was encoded in the genome of an uncultivated gamma-proteobacterium and shared highest amino acid sequence similarity with archaeal rhodopsins. The protein was functionally expressed in Escherichia coli and bound retinal to form an active, light-driven proton pump. The new rhodopsin exhibited a photochemical reaction cycle with intermediates and kinetics characteristic of archaeal proton-pumping rhodopsins. Our results demonstrate that archaeal-like rhodopsins are broadly distributed among different taxa, including members of the domain Bacteria. Our data also indicate that a previously unsuspected mode of bacterially mediated light-driven energy generation may commonly occur in oceanic surface waters worldwide.

Chromosome 2 sequence of the human malaria parasite Plasmodium falciparum.

Chromosome 2 of Plasmodium falciparum was sequenced; this sequence contains 947,103 base pairs and encodes 210 predicted genes. In comparison with the Saccharomyces cerevisiae genome, chromosome 2 has a lower gene density, introns are more frequent, and proteins are markedly enriched in nonglobular domains. A family of surface proteins, rifins, that may play a role in antigenic variation was identified. The complete sequencing of chromosome 2 has shown that sequencing of the A+T-rich P. falciparum genome is technically feasible.

THUMP--a predicted RNA-binding domain shared by 4-thiouridine, pseudouridine synthases and RNA methylases.

Sequence profile searches were used to identify an ancient domain in ThiI-like thiouridine synthases, conserved RNA methylases, archaeal pseudouridine synthases and several uncharacterized proteins. We predict that this domain is an RNA-binding domain that adopts an alpha/beta fold similar to that found in the C-terminal domain of translation initiation factor 3 and ribosomal protein S8.

The domains of death: evolution of the apoptosis machinery.

Recent progress in research into programmed cell death has resulted in the identification of the principal protein domains involved in this process. The evolution of many of these domains can be traced back in evolution to unicellular eukaryotes or even bacteria, where the domains appear to be involved in other regulatory functions. Cell-death systems in animals and plants share several conserved domains, in particular the family of apoptotic ATPases; this allows us to suggest a plausible, even if still incomplete, scenario for the evolution of apoptosis.

Genome of the extremely radiation-resistant bacterium Deinococcus radiodurans viewed from the perspective of comparative genomics.

The bacterium Deinococcus radiodurans shows remarkable resistance to a range of damage caused by ionizing radiation, desiccation, UV radiation, oxidizing agents, and electrophilic mutagens. D. radiodurans is best known for its extreme resistance to ionizing radiation; not only can it grow continuously in the presence of chronic radiation (6 kilorads/h), but also it can survive acute exposures to gamma radiation exceeding 1,500 kilorads without dying or undergoing induced mutation. These characteristics were the impetus for sequencing the genome of D. radiodurans and the ongoing development of its use for bioremediation of radioactive wastes. Although it is known that these multiple resistance phenotypes stem from efficient DNA repair processes, the mechanisms underlying these extraordinary repair capabilities remain poorly understood. In this work we present an extensive comparative sequence analysis of the Deinococcus genome. Deinococcus is the first representative with a completely sequenced genome from a distinct bacterial lineage of extremophiles, the Thermus-Deinococcus group. Phylogenetic tree analysis, combined with the identification of several synapomorphies between Thermus and Deinococcus, supports the hypothesis that it is an ancient group with no clear affinities to any of the other known bacterial lineages. Distinctive features of the Deinococcus genome as well as features shared with other free-living bacteria were revealed by comparison of its proteome to the collection of clusters of orthologous groups of proteins. Analysis of paralogs in Deinococcus has revealed several unique protein families. In addition, specific expansions of several other families including phosphatases, proteases, acyltransferases, and Nudix family pyrophosphohydrolases were detected. Genes that potentially affect DNA repair and recombination and stress responses were investigated in detail. Some proteins appear to have been horizontally transferred from eukaryotes and are not present in other bacteria. For example, three proteins homologous to plant desiccation resistance proteins were identified, and these are particularly interesting because of the correlation between desiccation and radiation resistance. Compared to other bacteria, the D. radiodurans genome is enriched in repetitive sequences, namely, IS-like transposons and small intergenic repeats. In combination, these observations suggest that several different biological mechanisms contribute to the multiple DNA repair-dependent phenotypes of this organism.

Prokaryotic genomes: the emerging paradigm of genome-based microbiology.

Comparative analysis of the complete sequences of seven bacterial and three archaeal genomes leads to the first generalizations of emerging genome-based microbiology. Protein sequences are, generally, highly conserved, with -70% of the gene products in bacteria and archaea containing ancient conserved regions. In contrast, there is little conservation of genome organization, except for a few essential operons. The most striking conclusions derived by comparison of multiple genomes from phylogenetically distant species are that the number of universally conserved gene families is very small and that multiple events of horizontal gene transfer and genome fusion are major forces in evolution.

Complete genome sequences of cellular life forms: glimpses of theoretical evolutionary genomics.

The availability of complete genome sequences of cellular life forms creates the opportunity to explore the functional content of the genomes and evolutionary relationships between them at a new qualitative level. With the advent of these sequences, the construction of a minimal gene set sufficient for sustaining cellular life and reconstruction of the genome of the last common ancestor of bacteria, eukaryotes, and archaea become realistic, albeit challenging, research projects. A version of the minimal gene set for modern-type cellular life derived by comparative analysis of two bacterial genomes, those of Haemophilus influenzae and Mycoplasma genitalium, consists of approximately 250 genes. A comparison of the protein sequences encoded in these genes with those of the proteins encoded in the complete yeast genome suggests that the last common ancestor of all extant life might have had an RNA genome.

Computational analysis of human disease-associated genes and their protein products.

The complete genome sequences for human, Drosophila melanogaster and Arabidopsis thaliana have been reported recently. With the availability of complete sequences for many bacteria and archaea, and five eukaryotes, comparative genomics and sequence analysis are enabling us to identify counterparts of many human disease genes in model organisms, which in turn should accelerate the pace of research and drug development to combat human diseases. Continuous improvement of specialized protein databases, together with sensitive computational tools, have enhanced the power and reliability of computational prediction of protein function.