Comparative analyses of fundamental differences in membrane transport capabilities in prokaryotes and eukaryotes

Whole-genome transporter analyses have been conducted on 141 organisms whose complete genome sequences are available. For each organism, the complete set of membrane transport systems was identified with predicted functions, and classified into protein families based on the transporter classification system. Organisms with larger genome sizes generally possessed a relatively greater number of transport systems. In prokaryotes and unicellular eukaryotes, the significant factor in the increase in transporter content with genome size was a greater diversity of transporter types. In contrast, in multicellular eukaryotes, greater number of paralogs in specific transporter families was the more important factor in the increase in transporter content with genome size. Both eukaryotic and prokaryotic intracellular pathogens and endosymbionts exhibited markedly limited transport capabilities. Hierarchical clustering of phylogenetic profiles of transporter families, derived from the presence or absence of a certain transporter family, showed that clustering patterns of organisms were correlated to both their evolutionary history and their overall physiology and lifestyles.     

Voltage-gated H+ channels associated with human phagocyte superoxide-generating NADPH oxidases: sequence comparisons,structural predictions,and phylogenetic analyses

The N-terminal domain of the human phagocyte flavocytochrome b558 NADPH oxidase, gp91phox, is believed to be a heme-containing voltage-gated H+ channel. The authors have conducted structural, sequence and phylogenetic analyses of the putative transmembrane channel/heme-binding domains of all homologous proteins in the NCBI GenBank database as of May 2001, as well as of the full-length proteins. Fifty-six homologues were identified, including 26 from animals, 19 from plants, seven from yeast, one from a slime mould and three from bacteria. Six well-defined sub-families were revealed by phylogenetic tree construction, two consisting of animal proteins, two of plant proteins, and one each of yeast and bacterial homologues, with the slime mould protein clustering loosely with one of the animal clusters. Signature sequences for the entire family as well as for the sub-families were determined. Most proteins have six putative TMSs, four of which may comprise the heme-binding H+ channel. The hydrophobic and amphipathic characteristics of each of the putative alpha-helical transmembrane segments were defined, and conserved residues that may be involved in heme binding, channel formation, and/or conformational changes were identified. The analyses lead to the suggestion that the oxidase domain became associated with the channel/heme-binding domain to form a single polypeptide chain early in evolutionary history, before eukaryotes diverged from prokaryotes, and that genetic transmission to present day organisms occurred primarily by vertical descent.     

Voltage-gated H + channels associated with human phagocyte superoxide-generating NADPH oxidases: Sequence comparisons,structural predictions,and phylogenetic analyses

The N-terminal domain of the human phagocyte flavocytochrome b 558 NADPH oxidase, gp91 phox, is believed to be a heme-containing voltage-gated H + channel. The authors have conducted structural, sequence and phylogenetic analyses of the putative transmembrane channel/heme-binding domains of all homologous proteins in the NCBI GenBank database as of May 2001, as well as of the full-length proteins. Fifty-six homologues were identified, including 26 from animals, 19 from plants, seven from yeast, one from a slime mould and three from bacteria. Six well-defined sub-families were revealed by phylogenetic tree construction, two consisting of animal proteins, two of plant proteins, and one each of yeast and bacterial homologues, with the slime mould protein clustering loosely with one of the animal clusters. Signature sequences for the entire family as well as for the sub-families were determined. Most proteins have six putative TMSs, four of which may comprise the heme-binding H + channel. The hydrophobic and amphipathic characteristics of each of the putative &#102 -helical transmembrane segments were defined, and conserved residues that may be involved in heme binding, channel formation, and/or conformational changes were identified. The analyses lead to the suggestion that the oxidase domain became associated with the channel/heme-binding domain to form a single polypeptide chain early in evolutionary history, before eukaryotes diverged from prokaryotes, and that genetic transmission to present day organisms occurred primarily by vertical descent.     

Phylogenetic analyses of the constituents of Type III protein secretion systems

Multicomponent Type III protein secretion systems transfer gram-negative bacterial virulence factors directly from the bacterial cytoplasm to the cytoplasm of a host eukaryotic cell in a process that may involve a single energy-coupled step. Extensive evidence supports the conclusion that the genetic apparatuses that encode these systems have been acquired independently by different gram-negative bacteria, presumably by lateral transfer. In this paper we conduct phylogenetic analyses of currently sequenced constituents of these systems and their homologues. The results reveal the relative relatedness of these systems and show that they evolved with little or no exchange of constituents between systems. This fact suggests that horizontal transmission of the genes encoding these systems always occurred as a unit without the formation of hybrid gene clusters. Moreover, homologous flagellar proteins show phylogenetic clustering that suggests that the flagellar systems and Type III protein secretory systems diverged from each other following very early duplication of a gene cluster sharing many (but not all) genes. Phylogenies of most or all of the flagellar proteins follow those of the source organisms with little or no lateral gene transfer suggesting that homologous flagellar proteins are true orthologues. We suggest that the flagellar apparatus was the evolutionary precursor of Type III protein secretion systems.     

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.     

Identification of rice TUBBY-like genes and their evolution

The identification of TUBBY-like genes in organisms ranging from single-celled to multicellular eukaryotes has allowed the phylogenetic history of this gene family to be traced back to the early evolutionary stages of eukaryote development. Rice TUBBY-like genes were located on chromosomes 1, 2, 3, 4, 5, 7, 8, 11 and 12 without any obvious clustering. On a genomic scale, it was revealed that the rice TUBBY-like gene family probably evolved mainly through segmental duplication produced by polyploidy. The altered selective constraints (or site-specific rate changes), related to functional divergence during protein evolution between plant and animal TUBBY-like genes, were statistically significant. Based on posterior probability analysis, five amino acid sites (103, 312, 315, 317 and 319) are thought to be responsible for functional divergence.     

Proteome-wide analysis of protein function composition reveals the clustering and phylogenetic properties of organisms

A 17-dimensional vector named the proteome vector is defined to represent an organism. The components of the vector reflect the relative contents of protein-encoding genes of the 17 cluster of orthologous groups of proteins (COGs) classes in the whole genome of the relevant organism. Based on the definition of this proteome vector, the fuzzy clustering of 36 completely sequenced organisms (8 archaea, 24 bacteria, and 4 eukarya) was performed and a proteome tree was constructed. Our results show that (1) the 36 organisms can be 100% correctly classified into three clusters corresponding to the three primary kingdoms, (2) our proteome tree is remarkably similar to that derived from 16S rRNA, and (3) the chromosomes and/or plasmids belonging to the same organism have very similar gene composition. Based on these results, we argue that the 17-dimensional proteome vector could be a good criterion for clustering approaches and to a large extent reveals the phylogenetic properties of organisms; the Three Primary Kingdoms Hypothesis is trustworthy although the existence of lateral gene transfer (LGT) brings controversy to the construction of the "universal tree of life."     

Evolution and cellular function of monothiol glutaredoxins: involvement in iron-sulphur cluster assembly

A number of bacterial species, mostly proteobacteria, possess monothiol glutaredoxins homologous to the Saccharomyces cerevisiae mitochondrial protein Grx5, which is involved in iron-sulphur cluster synthesis. Phylogenetic profiling is used to predict that bacterial monothiol glutaredoxins also participate in the iron-sulphur cluster (ISC) assembly machinery, because their phylogenetic profiles are similar to the profiles of the bacterial homologues of yeast ISC proteins. High evolutionary co-occurrence is observed between the Grx5 homologues and the homologues of the Yah1 ferredoxin, the scaffold proteins Isa1 and Isa2, the frataxin protein Yfh1 and the Nfu1 protein. This suggests that a specific functional interaction exists between these ISC machinery proteins. Physical interaction analyses using low-definition protein docking predict the formation of strong and specific complexes between Grx5 and several components of the yeast ISC machinery. Two-hybrid analysis has confirmed the in vivo interaction between Grx5 and Isa1. Sequence comparison techniques and cladistics indicate that the other two monothiol glutaredoxins of S. cerevisiae, Grx3 and Grx4, have evolved from the fusion of a thioredoxin gene with a monothiol glutaredoxin gene early in the eukaryotic lineage, leading to differential functional specialization. While bacteria do not contain these chimaeric glutaredoxins, in many eukaryotic species Grx5 and Grx3/4-type monothiol glutaredoxins coexist in the cell.     

Constructing phylogenetic trees using interacting pathways

Phylogenetic trees are used to represent evolutionary relationships among biological species or organisms. The construction of phylogenetic trees is based on the similarities or differences of their physical or genetic features. Traditional approaches of constructing phylogenetic trees mainly focus on physical features. The recent advancement of high-throughput technologies has led to accumulation of huge amounts of biological data, which in turn changed the way of biological studies in various aspects. In this paper, we report our approach of building phylogenetic trees using the information of interacting pathways. We have applied hierarchical clustering on two domains of organisms—eukaryotes and prokaryotes. Our preliminary results have shown the effectiveness of using the interacting pathways in revealing evolutionary relationships.     

The evolutionary origin of peroxisomes: an ER-peroxisome connection

The peroxisome is an essential eukaryotic organelle, crucial for lipid metabolism and free radical detoxification, development, differentiation, and morphogenesis from yeasts to humans. Loss of peroxisomes invariably leads to fatal peroxisome biogenesis disorders in man. The evolutionary origin of peroxisomes remains unsolved; proposals for either a symbiogenetic or cellular membrane invagination event are unconclusive. To address this question, we have probed with a peroxisomal proteome, an "ensemble" of 19 representative eukaryotic complete genomes. Molecular phylogenetic and sequence comparison tools allowed us to identify four proteins as peroxisomal markers for unequivocal in silico peroxisome detection. We have then detected the Apicomplexa phylum as the first group of organisms devoid of peroxisomes, in the presence of mitochondria. Finally, we deliver evidence against a prokaryotic ancestor of peroxisomes: (1) the peroxisomal membrane is composed of purely eukaryotic bricks and is thus useful to trace the eukaryotes in their evolutionary paths and (2) the peroxisomal matrix protein import system shares mechanistic similarities with the endoplasmic reticulum/proteasome degradation process, indicating a common evolutionary history.     

Reversing transmembrane electron flow: the DsbD and DsbB protein families

DsbD and DsbB are two proteins that in Escherichia coli catalyze transmembrane electron flow in opposite directions, thereby allowing reversible oxidoreduction of periplasmic dithiol/disulfide-containing proteins. We have identified all recognizable homologues of these two proteins in the databases and have conducted structural and phylogenetic analyses of the two families. The larger DsbD family is more diverse in sequence, topology, function and organismal distribution than the smaller DsbB family. DsbB homologues are rarely found outside of the proteobacteria, although DsbD homologues are found in many bacterial kingdoms as well as archaea and plant chloroplasts. Few organisms with a fully sequenced genome and a DsbB homologue lack a DsbD homologue, and most of these DsbD homologues fall within two clusters in the DsbD tree, exhibiting phylogenetic relationships that are the same as those observed for the DsbB proteins. These observations suggest that a subset of the DsbD homologues evolved in parallel with the DsbB family to perform a single unified function involving reversible extracytoplasmic protein dithiol-disulfide bond interchange. DsbD family proteins are shown to have arisen by an internal gene duplication event, and this observation leads to prediction of the pathway taken for the evolutionary appearance of the different protein topological types found within this family.     

Hallmarks of mycolic acid biosynthesis: a comparative genomics study

Mycolic acids, which render unique qualities to mycobacteria, are known to be important for mycobacterial growth, survival, and pathogenicity. It is of interest to understand the evolutionary origins of the mycolic acid pathway (MAP), as well as the common minimum principles critical for generating the capability of mycolic acid biosynthesis. The recent curation of a comprehensive model of the MAP in Mycobacterium tuberculosis and the availability of a large number of genome sequences make it feasible to carry out detailed sequence and phylogenetic analyses, to address these questions. A comprehensive phylogenetic pathway profile analysis was carried out for 318 fully sequenced bacterial genomes, for each of the proteins present in the MAP. The organisms were clustered on the basis of co-occurrence of the MAP proteins in their proteome, while the proteins were clustered on the basis of their phylogenetic profiles. The MAP proteins were also searched against the nonredundant sequence database, to identify similar proteins from other phyla. The pathway profiles indicate that four proteins and certain protein domains stand out as more characteristic to mycolate producing organisms. Further analysis leads to the identification of the desaturases DesA1 and DesA2 and certain domains of Fas and Pks13 as hallmarks of the pathway. The roles of these proteins in some other organisms, as well as the distribution of these proteins across all known genome sequences are also briefly discussed. The clustering of organisms, carried out to group organisms with similar profiles, provides a means of obtaining finer classification as compared to the standard taxonomic method. The results indicate that the MAP and hence the capacity of mycolic acid production in mycobacteria is an example of an emergent property that has come about by recruiting enzymes from unrelated pathways in plants, presumably through lateral gene transfer. The understanding of the hallmarks of mycolic acid biosynthesis will also find application in evaluating drug targets.     

Comparison of algorithms for prediction of related proteins using the method of phylogenetic profiles

Computational interactomics deals with prediction of functionally related proteins. One approach for solving this problem using comparative genomics consists in analysis of similarities between phylogenetic profiles of proteins. In contrast to most methods, which predict only pairwise interactions between proteins, in the present work we have applied the cluster analysis techniques in order to find modules of functionally related proteins. We have performed the cluster analysis of phylogenetic profiles of E. coli proteins using several clustering techniques and various modes for estimation of distances between profiles. We report here, that the best correspondence in the composition of resultant clusters to known metabolic pathways is achieved using Ward’s clustering together with Hamming’s distance. The proposed technique of assessing predictions of the modules of functionally related proteins can be used for comparative analysis of different algorithms for computational interactomics.     

Protein-translocating outer membrane porins of Gram-negative bacteria

Five families of outer membrane porins that function in protein secretion in Gram-negative bacteria are currently recognized. In this report, these five porin families are analyzed from structural and phylogenetic standpoints. They are the fimbrial usher protein (FUP), outer membrane factor (OMF), autotransporter (AT), two-partner secretion (TPS) and outer membrane secretin (Secretin) families. All members of these families in the current databases were identified, and all full-length homologues were multiply aligned for structural and phylogenetic analyses. The organismal distribution of homologues in each family proved to be unique with some families being restricted to proteobacteria and others being widespread in other bacterial kingdoms as well as eukaryotes. The compositions of and size differences between subfamilies provide evidence for specific orthologous relationships, which agree with available functional information and intra-subfamily phylogeny. The results reveal that horizontal transfer of genes encoding these proteins between phylogenetically distant organisms has been exceptionally rare although transfer within select bacterial kingdoms may have occurred. The resultant in silico analyses are correlated with available experimental evidence to formulate models relevant to the structures and evolutionary origins of these proteins.     

Yeasts illustrate the molecular mechanisms of eukaryotic genome evolution

Hemiascomycetous yeasts have the greatest number of sequenced species for a single phylum, and are at the forefront of evolutionary genomics of eukaryotes. Yeast genomes show the dynamic interplay between the formation and loss of genes and help to characterize the mechanisms involved and their functional and evolutionary consequences. These mechanisms have equivalents in the genomes of multicellular organisms. Yeast genomes show extensive loss of introns and a reduced role of transposable elements, and so probably have a more limited potential to form novel genes and functions than multicellular organisms, possibly explaining their conserved biological and morphological properties despite their considerable evolutionary range.     

Structure-Based Phylogenetic Analysis of the Lipocalin Superfamily

Lipocalins constitute a superfamily of extracellular proteins that are found in all three kingdoms of life. Although very divergent in their sequences and functions, they show remarkable similarity in 3-D structures. Lipocalins bind and transport small hydrophobic molecules. Earlier sequence-based phylogenetic studies of lipocalins highlighted that they have a long evolutionary history. However the molecular and structural basis of their functional diversity is not completely understood. The main objective of the present study is to understand functional diversity of the lipocalins using a structure-based phylogenetic approach. The present study with 39 protein domains from the lipocalin superfamily suggests that the clusters of lipocalins obtained by structure-based phylogeny correspond well with the functional diversity. The detailed analysis on each of the clusters and sub-clusters reveals that the 39 lipocalin domains cluster based on their mode of ligand binding though the clustering was performed on the basis of gross domain structure. The outliers in the phylogenetic tree are often from single member families. Also structure-based phylogenetic approach has provided pointers to assign putative function for the domains of unknown function in lipocalin family. The approach employed in the present study can be used in the future for the functional identification of new lipocalin proteins and may be extended to other protein families where members show poor sequence similarity but high structural similarity.