Orderly order in protein intrinsic disorder distribution: disorder in 3500 proteomes from viruses and the three domains of life

Intrinsically disordered proteins and intrinsically disordered protein regions are highly abundant in nature. However, the quantitative and qualitative measures of protein intrinsic disorder in species with known genomes are still not available. Furthermore, although the correlation between high fraction of disordered residues and advanced species has been reported, the details of this correlation and the connection between the disorder content and proteome complexity have not been reported as of yet. To fill this gap, we analysed entire proteomes of 3484 species from three domains of life (archaea, bacteria and eukaryotes) and from viruses. Our analysis revealed that the evolution process is characterized by distinctive patterns of changes in the protein intrinsic disorder content. We are showing here that viruses are characterized by the widest spread of the proteome disorder content (the percentage of disordered residues ranges from 7.3% in human coronavirus NL63 to 77.3% in Avian carcinoma virus). For several organisms, a clear correlation is seen between their disorder contents and habitats. In multicellular eukaryotes, there is a weak correlation between the complexity of an organism (evaluated as a number of different cell types) and its overall disorder content. For both the prokaryotes and eukaryotes, the disorder content is generally independent of the proteome size. However, disorder shows a sharp increase associated with the transition from prokaryotic to eukaryotic cells. This suggests that the increased disorder content in eukaryotic proteomes might be used by nature to deal with the increased cell complexity due to the appearance of the various cellular compartments.     

Dehydrogenases from all three domains of life cleave RNA

Specific interactions of glyceraldehyde-3-phosphate dehydrogenase (GAPDH) with RNA have been reported both in vitro and in vivo. We show that eukaryotic and bacterial GAPDH and two proteins from the hyperthermophilic archaeon Sulfolobus solfataricus, which are annotated as dehydrogenases, cleave RNA producing similar degradation patterns. RNA cleavage is most efficient at 60 degrees C, at MgCl(2) concentrations up to 5 mm, and takes place between pyrimidine and adenosine. The RNase active center of the putative aspartate semialdehyde dehydrogenase from S. solfataricus is located within the N-terminal 73 amino acids, which comprise the first mononucleotide-binding site of the predicted Rossmann fold. Thus, RNA cleavage has to be taken into account in the ongoing discussion of the possible biological function of RNA binding by dehydrogenases.     

Zinc through the three domains of life

Zinc is one of the metal ions essential for life, as it is required for the proper functioning of a large number of proteins. Despite its importance, the annotation of zinc-binding proteins in gene banks or protein domain databases still has significant room for improvement. In the present work, we compiled a list of known zinc-binding protein domains and of known zinc-binding sequence motifs (zinc-binding patterns), and then used them jointly to analyze the proteome of 57 different organisms to obtain an overview of zinc usage by archaeal, bacterial, and eukaryotic organisms. Zinc-binding proteins are an abundant fraction of these proteomes, ranging between 4% and 10%. The number of zinc-binding proteins correlates linearly with the total number of proteins encoded by the genome of an organism, but the proportionality constant of Eukaryota (8.8%) is significantly higher than that observed in Bacteria and Archaea (from 5% to 6%). Most of this enrichment is due to the larger portfolio of regulatory proteins in Eukaryota.     

Non-heme iron through the three domains of life

Metalloproteins are proteins capable of binding one or more metal ions, which are often required for their biological function or for regulation of their activities or for structural purposes. In high-throughput genome-level protein investigation efforts, such as Structural Genomics, the systematic experimental characterization of metal-binding properties (i.e. the investigation of the metalloproteome) is not always pursued, and remains far from trivial. In the present work we have applied a bioinformatic approach to investigate the occurrence of (putative) non-heme iron-binding proteins in 57 different organisms spanning the entire tree of life. It is found that the non-heme iron-proteome constitutes between 1% and 10% of the entire proteome of an organism. However, the iron-proteome constitutes a higher fraction of the proteome in archaea (on average 7.1% +/- 2.1%) than in bacteria (3.9% +/- 1.6%) and in eukaryota (1.1% +/- 0.4%). The analysis of the function of each putative iron-protein identified suggests that extant organisms have inherited the large majority of their iron-proteome from the last common ancestor.     

Noncellulosomal cohesin- and dockerin-like modules in the three domains of life

The high-affinity cohesin–dockerin interaction was originally discovered as modular components, which mediate the assembly of the various subunits of the multienzyme cellulosome complex that characterizes some cellulolytic bacteria. Until recently, the presence of cohesins and dockerins within a bacterial proteome was considered a definitive signature of a cellulosome-producing bacterium. Widespread genome sequencing has since revealed a wealth of putative cohesin- and dockerin-containing proteins in Bacteria, Archaea, and in primitive eukaryotes. The newly identified modules appear to serve diverse functions that are clearly distinct from the classical cellulosome archetype, and the vast majority of parent proteins are not predicted glycoside hydrolases. In most cases, only a few such genes have been identified in a given microorganism, which encode proteins containing but a single cohesin and/or dockerin. In some cases, one or the other module appears to be missing from a given species, and in other cases both modules occur within the same protein. This review provides a bioinformatics-based survey of the current status of cohesin- and dockerin-like sequences in species from the Bacteria, Archaea, and Eukarya. Surprisingly, many identified modules and their parent proteins are clearly unrelated to cellulosomes. The cellulosome paradigm may thus be the exception rather than the rule for bacterial, archaeal, and eukaryotic employment of cohesin and dockerin modules.     

Origins of DNA replication in the three domains of life

Replication of DNA is essential for the propagation of life. It is somewhat surprising then that, despite the vital nature of this process, cellular organisms show a great deal of variety in the mechanisms that they employ to ensure appropriate genome duplication. This diversity is manifested along classical evolutionary lines, with distinct combinations of replicon architecture and replication proteins being found in the three domains of life: the Bacteria, the Eukarya and the Archaea. Furthermore, although there are mechanistic parallels, even within a given domain of life, the way origins of replication are defined shows remarkable variation.     

Estimation of phylogenetic inconsistencies in the three domains of life

Discrepancies in phylogenetic trees of bacteria and archaea are often explained as lateral gene transfer events. However, such discrepancies may also be due to phylogenetic artifacts or orthology assignment problems. A first step that may help to resolve this dilemma is to estimate the extent of phylogenetic inconsistencies in trees of prokaryotes in comparison with those of higher eukaryotes, where no lateral gene transfer is expected. To test this, we used 21 proteomes each of eukaryotes (mainly opisthokonts), proteobacteria, and archaea that spanned equivalent levels of genetic divergence. In each domain of life, we defined a set of putative orthologous sequences using a phylogenetic-based orthology protocol and, as a reference topology, we used a tree constructed with concatenated genes of each domain. Our results show, for most of the tests performed, that the magnitude of topological inconsistencies with respect to the reference tree was very similar in the trees of proteobacteria and eukaryotes. When clade support was taken into account, prokaryotes showed some more inconsistencies, but then all values were very low. Discrepancies were only consistently higher in archaea but, as shown by simulation analysis, this is likely due to the particular tree of the archaeal species used here being more difficult to reconstruct, whereas the trees of proteobacteria and eukaryotes were of similar difficulty. Although these results are based on a relatively small number of genes, it seems that phylogenetic reconstruction problems, including orthology assignment problems, have a similar overall effect over prokaryotic and eukaryotic trees based on single genes. Consequently, lateral gene transfer between distant prokaryotic species may have been more rare than previously thought, which opens the way to obtain the tree of life of bacterial and archaeal species using genomic data and the concatenation of adequate genes, in the same way as it is usually done in eukaryotes.     

Transitional forms between the three domains of life and evolutionary implications

The question as to the origin and relationship between the three domains of life is lodged in a phylogenetic impasse. The dominant paradigm is to see the three domains as separated. However, the recently characterized bacterial species have suggested continuity between the three domains. Here, we review the evidence in support of this hypothesis and evaluate the implications for and against the models of the origin of the three domains of life. The existence of intermediate steps between the three domains discards the need for fusion to explain eukaryogenesis and suggests that the last universal common ancestor was complex. We propose a scenario in which the ancestor of the current bacterial Planctomycetes, Verrucomicrobiae and Chlamydiae superphylum was related to the last archaeal and eukaryotic common ancestor, thus providing a way out of the phylogenetic impasse.     

Spatial and temporal variability in a stratified hypersaline microbial mat community

Hypersaline microbial mat communities have recently been shown to be more diverse than once thought. The variability in community composition of hypersaline mats, both in terms of spatial and temporal dimensions, is still poorly understood. Because this information is essential to understanding the complex biotic and abiotic interactions within these communities, terminal restriction fragment analysis and 16S rRNA gene sequencing were used to characterize the near-surface community of a hypersaline microbial mat in Guerrero Negro, Mexico. Core samples were analyzed to assay community variability over large regional scales (centimeter to kilometer) and to track depth-related changes in population distribution at 250-μm intervals over a diel period. Significant changes in total species diversity were observed at increasing distances across the mat surface; however, key species (e.g. Microcoleus sp.) were identified throughout the mat. The vertical position and abundance of >50% of the 60 peaks detected varied dramatically over a diel cycle, including Beggiatoa sp., cyanobacteria, Chloroflexus sp., Halochromatium sp., Bacteroidetes sp. and several as-yet-identified bacteria. Many of these migrations correlated strongly with diel changes in redox conditions within the mat, contributing to strong day–night community structure differences.     

A network-based approach to classify the three domains of life

Background

Chromosomal orthologs can reveal the shared ancestral gene set and their evolutionary trends. Additionally, physico-chemical properties of encoded proteins could provide information about functional adaptation and ecological niche requirements.

Results

We analyzed 7080 genes (five groups of 1416 orthologs each) from Rhizobiales species (S. meliloti, R. etli, and M. loti, plant symbionts; A. tumefaciens, a plant pathogen; and B. melitensis, an animal pathogen). We evaluated their phylogenetic relationships and observed three main topologies. The first, with closer association of R. etli to A. tumefaciens; the second with R. etli closer to S. meliloti; and the third with A. tumefaciens and S. meliloti as the closest pair. This was not unusual, given the close relatedness of these three species. We calculated the synonymous (dS) and nonsynonymous (dN) substitution rates of these orthologs, and found that informational and metabolic functions showed relatively low dN rates; in contrast, genes from hypothetical functions and cellular processes showed high dN rates. An alternative measure of sequence variability, percentage of changes by species, was used to evaluate the most specific proportion of amino acid residues from alignments. When dN was compared with that measure a high correlation was obtained, revealing that much of evolutive information was extracted with the percentage of changes by species at the amino acid level. By analyzing the sequence variability of orthologs with a set of five properties (polarity, electrostatic charge, formation of secondary structures, molecular volume, and amino acid composition), we found that physico-chemical characteristics of proteins correlated with specific functional roles, and association of species did not follow their typical phylogeny, probably reflecting more adaptation to their life styles and niche preferences. In addition, orthologs with low dN rates had residues with more positive values of polarity, volume and electrostatic charge.

Conclusions

These findings revealed that even when orthologs perform the same function in each genomic background, their sequences reveal important evolutionary tendencies and differences related to adaptation. This article was reviewed by: Dr. Purificación López-García, Prof. Jeffrey Townsend (nominated by Dr. J. Peter Gogarten), and Ms. Olga Kamneva.     

Cladogenesis, coalescence and the evolution of the three domains of life

In this article, we explore the large-scale structure of the tree of life by using a simple model with a constant number of species and rates of speciation that equal the rates of extinction. In addition, we discuss the consequences of horizontal gene transfer for the concept of a most recent common ancestor of all living organisms (cenancestor). A simple null hypothesis based on coalescence theory explains some features of the observed topologies of the tree of life. Simulations of genes and organismal lineages suggest that there was no single common ancestor that contained all the genes ancestral to those shared among the three domains of life. Each contemporary molecule has its own history that traces back to an individual molecular cenancestor. However, these molecular ancestors were likely to be present in different organisms and at different times.