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.

Microbial culturomics: paradigm shift in the human gut microbiome study.

Comprehensive determination of the microbial composition of the gut microbiota and the relationships with health and disease are major challenges in the 21st century. Metagenomic analysis of the human gut microbiota detects mostly uncultured bacteria. We studied stools from two lean Africans and one obese European, using 212 different culture conditions (microbial culturomics), and tested the colonies by using mass spectrometry and 16S rRNA amplification and sequencing. In parallel, we analysed the same three samples by pyrosequencing 16S rRNA amplicons targeting the V6 region. The 32 500 colonies obtained by culturomics have yielded 340 species of bacteria from seven phyla and 117 genera, including two species from rare phyla (Deinococcus-Thermus and Synergistetes, five fungi, and a giant virus (Senegalvirus). The microbiome identified by culturomics included 174 species never described previously in the human gut, including 31 new species and genera for which the genomes were sequenced, generating c. 10 000 new unknown genes (ORFans), which will help in future molecular studies. Among these, the new species Microvirga massiliensis has the largest bacterial genome so far obtained from a human, and Senegalvirus is the largest virus reported in the human gut. Concurrent metagenomic analysis of the same samples produced 698 phylotypes, including 282 known species, 51 of which overlapped with the microbiome identified by culturomics. Thus, culturomics complements metagenomics by overcoming the depth bias inherent in metagenomic approaches.

The phylogenetic forest and the quest for the elusive tree of life.

Extensive horizontal gene transfer (HGT) among prokaryotes seems to undermine the tree of life (TOL) concept. However, the possibility remains that the TOL can be salvaged as a statistical central trend in the phylogenetic "forest of life" (FOL). A comprehensive comparative analysis of 6901 phylogenetic trees for prokaryotic genes revealed a signal of vertical inheritance that was particularly strong among the 102 nearly universal trees (NUTs), despite the high topological inconsistency among the trees in the FOL, most likely, caused by HGT. The topologies of the NUTs are similar to the topologies of numerous other trees in the FOL; although the NUTs cannot represent the FOL completely, they reflect a significant central trend. Thus, the original TOL concept becomes obsolete but the idea of a "weak" TOL as the dominant trend in the FOL merits further investigation. The totality of gene trees comprising the FOL appears to be a natural representation of the history of life given the inherent tree-like character of the replication process.

Hax1 lacks BH modules and is peripherally associated to heavy membranes: implications for Omi/HtrA2 and PARL activity in the regulation of mitochondrial stress and apoptosis.

Hax1 has an important role in immunodeficiency syndromes and apoptosis. A recent report (Chao et al., Nature, 2008) proposed that the Bcl-2-family-related protein, Hax1, suppresses apoptosis in lymphocytes and neurons through a mechanism that involves its association to the inner mitochondrial membrane rhomboid protease PARL, to proteolytically activate the serine protease Omi/HtrA2 and eliminate active Bax. This model implies that the control of cell-type sensitivity to pro-apoptotic stimuli is governed by the PARL/Hax1 complex in the mitochondria intermembrane space and, more generally, that Bcl-2-family-related proteins can control mitochondrial outer-membrane permeabilization from inside the mitochondrion. Further, it defines a novel, anti-apoptotic Opa1-independent pathway for PARL. In this study, we present evidence that, in vivo, the activity of Hax1 cannot be mechanistically coupled to PARL because the two proteins are confined in distinct cellular compartments and their interaction in vitro is an artifact. We also show by sequence analysis and secondary structure prediction that Hax1 is extremely unlikely to be a Bcl-2-family-related protein because it lacks Bcl-2 homology modules. These results indicate a different function and mechanism of Hax1 in apoptosis and re-opens the question of whether mammalian PARL, in addition to apoptosis, regulates mitochondrial stress response through Omi/HtrA2 processing.

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.

A bipolar DNA helicase gene, herA, clusters with rad50, mre11 and nurA genes in thermophilic archaea.

We showed previously that rad50 and mre11 genes of thermophilic archaea are organized in an operon-like structure with a third gene (nurA) encoding a 5' to 3' exonuclease. Here, we show that the rad50, mre11 and nurA genes from the hyperthermophilic archaeon Sulfolobus acidocaldarius are co-transcribed with a fourth gene encoding a DNA helicase. This enzyme (HerA) is the prototype of a new class of DNA helicases able to utilize either 3' or 5' single-stranded DNA extensions for loading and subsequent DNA duplex unwinding. To our knowledge, DNA helicases capable of translocating along the DNA in both directions have not been identified previously. Sequence analysis of HerA shows that it is a member of the TrwB, FtsK and VirB4/VirD4 families of the PilT class NTPases. HerA homologs are found in all thermophilic archaeal species and, in all cases except one, the rad50, mre11, nurA and herA genes are grouped together. These results suggest that the archaeal Rad50-Mre11 complex might act in association with a 5' to 3' exonuclease (NurA) and a bipolar DNA helicase (HerA) indicating a probable involvement in the initiation step of homologous recombination.

Scores of RINGS but no PHDs in ubiquitin signaling.

Recently, it has been reported that PHD fingers of MEKK1 kinase and a family of viral and cellular membrane proteins have E3 ubiquitin ligase activity. Here we describe unique sequence and structural signatures that distinguish PHD fingers from RING fingers, which function primarily as E3 ubiquitin ligases, and demonstrate that the Zn-binding modules of the above proteins are distinct versions of the RING domain rather than PHD fingers. Thus, currently available data reveal extreme versatility of RINGs and their derivatives that function as E3 ubiquitin ligases but provide no evidence of this activity among PHD fingers whose principal function appears to involve specific protein-protein and possibly protein-DNA interactions in chromatin.

Origin and evolution of eukaryotic apoptosis: the bacterial connection.

The availability of numerous complete genome sequences of prokaryotes and several eukaryotic genome sequences provides for new insights into the origin of unique functional systems of the eukaryotes. Several key enzymes of the apoptotic machinery, including the paracaspase and metacaspase families of the caspase-like protease superfamily, apoptotic ATPases and NACHT family NTPases, and mitochondrial HtrA-like proteases, have diverse homologs in bacteria, but not in archaea. Phylogenetic analysis strongly suggests a mitochondrial origin for metacaspases and the HtrA-like proteases, whereas acquisition from Actinomycetes appears to be the most likely scenario for AP-ATPases. The homologs of apoptotic proteins are particularly abundant and diverse in bacteria that undergo complex development, such as Actinomycetes, Cyanobacteria and alpha-proteobacteria, the latter being progenitors of the mitochondria. In these bacteria, the apoptosis-related domains typically form multidomain proteins, which are known or inferred to participate in signal transduction and regulation of gene expression. Some of these bacterial multidomain proteins contain fusions between apoptosis-related domains, such as AP-ATPase fused with a metacaspase or a TIR domain. Thus, bacterial homologs of eukaryotic apoptotic machinery components might functionally and physically interact with each other as parts of signaling pathways that remain to be investigated. An emerging scenario of the origin of the eukaryotic apoptotic system involves acquisition of several central apoptotic effectors as a consequence of mitochondrial endosymbiosis and probably also as a result of subsequent, additional horizontal gene transfer events, which was followed by recruitment of newly emerging eukaryotic domains as adaptors.

Origin of alternative splicing by tandem exon duplication.

Genes with new functions often evolve by gene duplication. Alternative splicing is another means of evolutionary innovation in eukaryotes, which allows a single gene to encode functionally diverse proteins. We investigate a connection between these two evolutionary phenomena. For approximately 10% of the described cases of substitution alternative splicing, such that either one or another amino acid sequence is included into the protein, evidence of origin by tandem exon duplication was found. This is a conservative estimate because alternative exons are typically short and, on many occasions, duplicates may have diverged beyond recognition. Dating exon duplications through a combination of the available experimental data on alternative splicing in orthologous genes from different species and computational analysis indicates that most of the duplications antedate at least the radiation of mammalian orders or even the radiation of vertebrate classes. At present, tandem exon duplication is the only mechanism of evolution of substitution alternative splicing that can be specifically demonstrated. Along with gene duplication, this could be a major route for generating functional diversity during evolution of multicellular eukaryotes.

Genome trees constructed using five different approaches suggest new major bacterial clades.

BACKGROUND: The availability of multiple complete genome sequences from diverse taxa prompts the development of new phylogenetic approaches, which attempt to incorporate information derived from comparative analysis of complete gene sets or large subsets thereof. Such attempts are particularly relevant because of the major role of horizontal gene transfer and lineage-specific gene loss, at least in the evolution of prokaryotes. RESULTS: Five largely independent approaches were employed to construct trees for completely sequenced bacterial and archaeal genomes: i) presence-absence of genomes in clusters of orthologous genes; ii) conservation of local gene order (gene pairs) among prokaryotic genomes; iii) parameters of identity distribution for probable orthologs; iv) analysis of concatenated alignments of ribosomal proteins; v) comparison of trees constructed for multiple protein families. All constructed trees support the separation of the two primary prokaryotic domains, bacteria and archaea, as well as some terminal bifurcations within the bacterial and archaeal domains. Beyond these obvious groupings, the trees made with different methods appeared to differ substantially in terms of the relative contributions of phylogenetic relationships and similarities in gene repertoires caused by similar life styles and horizontal gene transfer to the tree topology. The trees based on presence-absence of genomes in orthologous clusters and the trees based on conserved gene pairs appear to be strongly affected by gene loss and horizontal gene transfer. The trees based on identity distributions for orthologs and particularly the tree made of concatenated ribosomal protein sequences seemed to carry a stronger phylogenetic signal. The latter tree supported three potential high-level bacterial clades,: i) Chlamydia-Spirochetes, ii) Thermotogales-Aquificales (bacterial hyperthermophiles), and ii) Actinomycetes-Deinococcales-Cyanobacteria. The latter group also appeared to join the low-GC Gram-positive bacteria at a deeper tree node. These new groupings of bacteria were supported by the analysis of alternative topologies in the concatenated ribosomal protein tree using the Kishino-Hasegawa test and by a census of the topologies of 132 individual groups of orthologous proteins. Additionally, the results of this analysis put into question the sister-group relationship between the two major archaeal groups, Euryarchaeota and Crenarchaeota, and suggest instead that Euryarchaeota might be a paraphyletic group with respect to Crenarchaeota. CONCLUSIONS: We conclude that, the extensive horizontal gene flow and lineage-specific gene loss notwithstanding, extension of phylogenetic analysis to the genome scale has the potential of uncovering deep evolutionary relationships between prokaryotic lineages.

Common origin of four diverse families of large eukaryotic DNA viruses.

Comparative analysis of the protein sequences encoded in the genomes of three families of large DNA viruses that replicate, completely or partly, in the cytoplasm of eukaryotic cells (poxviruses, asfarviruses, and iridoviruses) and phycodnaviruses that replicate in the nucleus reveals 9 genes that are shared by all of these viruses and 22 more genes that are present in at least three of the four compared viral families. Although orthologous proteins from different viral families typically show weak sequence similarity, because of which some of them have not been identified previously, at least five of the conserved genes appear to be synapomorphies (shared derived characters) that unite these four viral families, to the exclusion of all other known viruses and cellular life forms. Cladistic analysis with the genes shared by at least two viral families as evolutionary characters supports the monophyly of poxviruses, asfarviruses, iridoviruses, and phycodnaviruses. The results of genome comparison allow a tentative reconstruction of the ancestral viral genome and suggest that the common ancestor of all of these viral families was a nucleocytoplasmic virus with an icosahedral capsid, which encoded complex systems for DNA replication and transcription, a redox protein involved in disulfide bond formation in virion membrane proteins, and probably inhibitors of apoptosis. The conservation of the disulfide-oxidoreductase, a major capsid protein, and two virion membrane proteins indicates that the odd-shaped virions of poxviruses have evolved from the more common icosahedral virion seen in asfarviruses, iridoviruses, and phycodnaviruses.

Rapid evolution of a cyclin A inhibitor gene, roughex, in Drosophila.

The recent sequencing of the complete genome of the fruit fly Drosophila melanogaster has yielded about 30% of the predicted genes with no obvious counterparts in other organisms. These rapidly evolving genes remain largely unexplored. Here, we present evidence for a striking variability in an important Drosophila cell cycle regulator encoded by the gene roughex (rux) in closely related fly species. The unusual level of Rux protein variability indicates that there are very low overall constraints on amino acid substitutions. Despite the lack of sequence similarity, certain common features, including the presence of a C-terminal nuclear localization signal and a functionally important N-terminal RXL cyclin-binding motif, exist between Rux and cyclin-dependent kinase inhibitors of the Cip/Kip family. These results indicate that even some genes involved in key regulatory processes in eukaryotes evolve at extremely high rates.

Independent evolution of heavy metal-associated domains in copper chaperones and copper-transporting atpases.

Copper chaperones are small cytoplasmic proteins that bind intracellular copper (Cu) and deliver it to Cu-dependent enzymes such as cytochrome oxidase, superoxide dismutase, and amine oxidase. Copper chaperones are similar in sequence and structure to the Cu-binding heavy metal-associated (HMA) domains of Cu-transporting ATPases (Cu-ATPases), and the genes for copper chaperones and Cu-ATPases are often located in the same operon. Phylogenetic analysis shows that Cu chaperones and HMA domains of Cu-ATPases represent ancient and distinct lineages that have evolved largely independently since their initial separation. Copper chaperone-Cu-ATPase operons appear to have evolved independently in different prokaryotic lineages, probably due to a strong selective pressure for coexpression of these genes.

Presence of ATG triplets in 5' untranslated regions of eukaryotic cDNAs correlates with a 'weak' context of the start codon.

MOTIVATION: The context of the start codon (typically, AUG) and the features of the 5' Untranslated Regions (5' UTRs) are important for understanding translation regulation in eukaryotic mRNAs and for accurate prediction of the coding region in genomic and cDNA sequences. The presence of AUG triplets in 5' UTRs (upstream AUGs) might effect the initiation rate and, in the context of gene prediction, could reduce the accuracy of the identification of the authentic start. To reveal potential connections between the presence of upstream AUGs and other features of 5' UTRs, such as their length and the start codon context, we undertook a systematic analysis of the available eukaryotic 5' UTR sequences. RESULTS: We show that a large fraction of 5' UTRs in the available cDNA sequences, 15-53% depending on the organism, contain upstream ATGs. A negative correlation was observed between the information content of the translation start signal and the length of the 5' UTR. Similarly, a negative correlation exists between the 'strength' of the start context and the number of upstream ATGs. Typically, cDNAs containing long 5' UTRs with multiple upstream ATGs have a 'weak' start context, and in contrast, cDNAs containing short 5' UTRs without ATGs have 'strong' starts. These counter-intuitive results may be interpreted in terms of upstream AUGs having an important role in the regulation of translation efficiency by ensuring low basal translation level via double negative control and creating the potential for additional regulatory mechanisms. One of such mechanisms, supported by experimental studies of some mRNAs, includes removal of the AUG-containing portion of the 5' UTR by alternative splicing. AVAILABILITY: An ATG_ EVALUATOR program is available upon request or at www.itba.mi.cnr.it/webgene. CONTACT: rogozin@ncbi.nlm.nih.gov, milanesi@itba.mi.cnr.it.

Horizontal gene transfer in prokaryotes: quantification and classification.

Comparative analysis of bacterial, archaeal, and eukaryotic genomes indicates that a significant fraction of the genes in the prokaryotic genomes have been subject to horizontal transfer. In some cases, the amount and source of horizontal gene transfer can be linked to an organism's lifestyle. For example, bacterial hyperthermophiles seem to have exchanged genes with archaea to a greater extent than other bacteria, whereas transfer of certain classes of eukaryotic genes is most common in parasitic and symbiotic bacteria. Horizontal transfer events can be classified into distinct categories of acquisition of new genes, acquisition of paralogs of existing genes, and xenologous gene displacement whereby a gene is displaced by a horizontally transferred ortholog from another lineage (xenolog). Each of these types of horizontal gene transfer is common among prokaryotes, but their relative contributions differ in different lineages. The fixation and long-term persistence of horizontally transferred genes suggests that they confer a selective advantage on the recipient organism. In most cases, the nature of this advantage remains unclear, but detailed examination of several cases of acquisition of eukaryotic genes by bacteria seems to reveal the evolutionary forces involved. Examples include isoleucyl-tRNA synthetases whose acquisition from eukaryotes by several bacteria is linked to antibiotic resistance, ATP/ADP translocases acquired by intracellular parasitic bacteria, Chlamydia and Rickettsia, apparently from plants, and proteases that may be implicated in chlamydial pathogenesis.