Mythical origins of the actin cytoskeleton

The origin of the eukaryotic cell is one of the greatest mysteries in modern biology. Eukaryotic-wide specific biological processes arose in the lost ancestors of eukaryotes. These distinctive features, such as the actin cytoskeleton, define what it is to be a eukaryote. Recent sequencing, characterization, and isolation of Asgard archaea have opened an intriguing window into the pre-eukaryotic cell. Firstly, sequencing of anaerobic sediments identified a group of uncultured organisms, Asgard archaea, which contain genes with homology to eukaryotic signature genes. Secondly, characterization of the products of these genes at the protein level demonstrated that Asgard archaea have related biological processes to eukaryotes. Finally, the isolation of an Asgard archaeon has produced a model organism in which the morphological consequences of the eukaryotic-like processes can be studied. Here, we consider the consequences for the Asgard actin cytoskeleton and for the evolution of a regulated actin system in the archaea-to-eukaryotic transition.     

Phylogenomic evidence for a myxococcal contribution to the mitochondrial fatty acid beta-oxidation

Background

The origin of eukaryotes remains a fundamental question in evolutionary biology. Although it is clear that eukaryotic genomes are a chimeric combination of genes of eubacterial and archaebacterial ancestry, the specific ancestry of most eubacterial genes is still unknown. The growing availability of microbial genomes offers the possibility of analyzing the ancestry of eukaryotic genomes and testing previous hypotheses on their origins.

Methodology/Principal Findings

Here, we have applied a phylogenomic analysis to investigate a possible contribution of the Myxococcales to the first eukaryotes. We conducted a conservative pipeline with homologous sequence searches against a genomic sampling of 40 eukaryotic and 357 prokaryotic genomes. The phylogenetic reconstruction showed that several eukaryotic proteins traced to Myxococcales. Most of these proteins were associated with mitochondrial lipid intermediate pathways, particularly enzymes generating reducing equivalents with pivotal roles in fatty acid β-oxidation metabolism. Our data suggest that myxococcal species with the ability to oxidize fatty acids transferred several genes to eubacteria that eventually gave rise to the mitochondrial ancestor. Later, the eukaryotic nucleocytoplasmic lineage acquired those metabolic genes through endosymbiotic gene transfer.

Conclusions/Significance

Our results support a prokaryotic origin, different from α-proteobacteria, for several mitochondrial genes. Our data reinforce a fluid prokaryotic chromosome model in which the mitochondrion appears to be an important entry point for myxococcal genes to enter eukaryotes.     

Stabilization against hyperthermal denaturation through increased CG content can explain the discrepancy between whole genome and 16S rRNA analyses

Based largely upon analysis of ribosomal RNA, a third domain of life, called archaea, had been proposed in addition to bacteria and eukaryotes. However, quantitative analysis of 73 whole genomes shows only a two-domain division of life: into eukaryotes and prokaryotes. Thousands of orthologous genes in archaea and bacteria show an essentially unimodal distribution of sequence identities. Thus, whole genome analyses indicate that archaea are a phylum of bacteria rather than a separate domain of life. In contrast, archaeal rRNA and that of hyperthermophilic bacteria differ from the rRNA of mesophilic bacteria. Thus, there is a bimodal distribution of rRNA sequence identities which differ by 12%. This discrepancy in rRNA and gene content based analyses of whole genomes is likely due to a 15% elevated C:G content of the rRNA of archaea and hyperthermophilic bacteria. The elevated C:G content is consistent with stabilization against thermal denaturation caused by additional hydrogen bonding (3 bonds) in C:G pairs compared to A:U pairs (2 bonds). Based upon this premise, there is no reliable way to correct rRNA for such differences in base composition and it is not possible to quantitatively compare hyperthermophiles with mesophiles by the rRNA method. Furthermore, quantitative study of whole genomes shows that the extent of change in both bacterial and archaeal genes, including rRNA, has reached a limit. Thus, direct sequence comparisons work with closely related genomes, but it is not possible to differentiate the most divergent prokaryotic species, which are currently designated as separate phyla. We believe that the differences in characteristics of archaeal species is based primarily upon selection of genes and pathways compatible with the extreme environmental lifestyle, i.e., hyperthermophily.     

Phylogenomics of prokaryotic ribosomal proteins

Archaeal and bacterial ribosomes contain more than 50 proteins, including 34 that are universally conserved in the three domains of cellular life (bacteria, archaea, and eukaryotes). Despite the high sequence conservation, annotation of ribosomal (r-) protein genes is often difficult because of their short lengths and biased sequence composition. We developed an automated computational pipeline for identification of r-protein genes and applied it to 995 completely sequenced bacterial and 87 archaeal genomes available in the RefSeq database. The pipeline employs curated seed alignments of r-proteins to run position-specific scoring matrix (PSSM)-based BLAST searches against six-frame genome translations, mitigating possible gene annotation errors. As a result of this analysis, we performed a census of prokaryotic r-protein complements, enumerated missing and paralogous r-proteins, and analyzed the distributions of ribosomal protein genes among chromosomal partitions. Phyletic patterns of bacterial and archaeal r-protein genes were mapped to phylogenetic trees reconstructed from concatenated alignments of r-proteins to reveal the history of likely multiple independent gains and losses. These alignments, available for download, can be used as search profiles to improve genome annotation of r-proteins and for further comparative genomics studies.     

Spirochete contributions to the eukaryotic genome

Expanded genome/proteome databases and effective use of sequence alignment tools make it possible to trace the phylogeny of individual eukaryotic proteins and ultimately to identify the prokaryotes that contributed to the last eukaryotic common ancestor (LECA). I developed an application of reciprocal BLASTp that identifies (1) the prokaryotic lineages that have contributed to the nuclear genome and (2) the specific proteins acquired from prokaryotic ancestors. Eight complete eubacterial proteomes were analyzed: two free-living spirochetes, two clostridia, two actinobacteria, and two proteobacteria (one alpha and one gamma). The data reveal a spirochete genetic contribution to the eukaryotic genome including essential proteins involved in DNA binding and repair, cyclic nucleotide metabolism, acyltransferase, and signal transduction. My results, consistent with the sulfur syntrophy hypothesis that posits LECA evolved from a merger of spirochetes (eubacteria) with sulfidogenic eocytes (archaebacteria), confirm the contribution of mitochondrial genes from alpha-proteobacteria. A contribution from clostridia to eukaryote genomes was also detected whereas none was seen from either actinobacterium or Escherichia coli. The complete spirochete and clostridial genetic contributions to eukaryotes and those of other eu-and archaebacteria can be identified by this method.     

Evidence of a large novel gene pool associated with prokaryotic genomic islands

Microbial genes that are “novel” (no detectable homologs in other species) have become of increasing interest as environmental sampling suggests that there are many more such novel genes in yet-to-be-cultured microorganisms. By analyzing known microbial genomic islands and prophages, we developed criteria for systematic identification of putative genomic islands (clusters of genes of probable horizontal origin in a prokaryotic genome) in 63 prokaryotic genomes, and then characterized the distribution of novel genes and other features. All but a few of the genomes examined contained significantly higher proportions of novel genes in their predicted genomic islands compared with the rest of their genome (Paired t test = 4.43E-14 to 1.27E-18, depending on method). Moreover, the reverse observation (i.e., higher proportions of novel genes outside of islands) never reached statistical significance in any organism examined. We show that this higher proportion of novel genes in predicted genomic islands is not due to less accurate gene prediction in genomic island regions, but likely reflects a genuine increase in novel genes in these regions for both bacteria and archaea. This represents the first comprehensive analysis of novel genes in prokaryotic genomic islands and provides clues regarding the origin of novel genes. Our collective results imply that there are different gene pools associated with recently horizontally transmitted genomic regions versus regions that are primarily vertically inherited. Moreover, there are more novel genes within the gene pool associated with genomic islands. Since genomic islands are frequently associated with a particular microbial adaptation, such as antibiotic resistance, pathogen virulence, or metal resistance, this suggests that microbes may have access to a larger “arsenal” of novel genes for adaptation than previously thought.     

A Novel Prokaryotic Promoter Identified in the Genome of Some Monopartite Begomoviruses

Geminiviruses are known to exhibit both prokaryotic and eukaryotic features in their genomes, with the ability to express their genes and even replicate in bacterial cells. We have demonstrated previously the existence of unit-length single-stranded circular DNAs of Ageratum yellow vein virus (AYVV, a species in the genus Begomovirus, family Geminiviridae) in Escherichia coli cells, which prompted our search for unknown prokaryotic functions in the begomovirus genomes. By using a promoter trapping strategy, we identified a novel prokaryotic promoter, designated AV3 promoter, in nts 762-831 of the AYVV genome. Activity assays revealed that the AV3 promoter is strong, unidirectional, and constitutive, with an endogenous downstream ribosome binding site and a translatable short open reading frame of eight amino acids. Sequence analyses suggested that the AV3 promoter might be a remnant of prokaryotic ancestors that could be related to certain promoters of bacteria from marine or freshwater environments. The discovery of the prokaryotic AV3 promoter provided further evidence for the prokaryotic origin in the evolutionary history of geminiviruses.     

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.     

On the origin of eukaryotic cytoskeleton.

The origin of eukaryote-specific cytoskeletal proteins is an issue which is closely related to the origin of the domain Eukarya. As nearly all of these proteins are not found in prokaryotes, the prokaryotic origin of eukaryotic cytoskeletal network suggested by most models is questionable. Eukaryotic cytoskeletal proteins might descend from subpopulations of pre-cells co-existing with Bacteria and Archaea prior to the origin of eukaryotes. The pre-karyote (the host for a-proteobacterial ancestors of mitochondria) might have already possessed eukaryotic-like cytoskeleton. A possible role for viruses in the origin of eukaryotic cytoskeletal proteins is discussed. Viruses parasitizing on pre-cells and/or on the pre-karyote might have themselves used several eukaryotic-like cytoskeletal proteins for segregation and packing of their genomes.     

The impact of genomics on research in diversity and evolution of archaea

Since the definition of archaea as a separate domain of life along with bacteria and eukaryotes, they have become one of the most interesting objects of modern microbiology, molecular biology, and biochemistry. Sequencing and analysis of archaeal genomes were especially important for studies on archaea because of a limited availability of genetic tools for the majority of these microorganisms and problems associated with their cultivation. Fifteen years since the publication of the first genome of an archaeon, more than one hundred complete genome sequences of representatives of different phylogenetic groups have been determined. Analysis of these genomes has expanded our knowledge of biology of archaea, their diversity and evolution, and allowed identification and characterization of new deep phylogenetic lineages of archaea. The development of genome technologies has allowed sequencing the genomes of uncultivated archaea directly from enrichment cultures, metagenomic samples, and even from single cells. Insights have been gained into the evolution of key biochemical processes in archaea, such as cell division and DNA replication, the role of horizontal gene transfer in the evolution of archaea, and new relationships between archaea and eukaryotes have been revealed.     

A universal compositional correlation among codon positions.

We have investigated the compositional distributions of third codon positions of genes from the 16 prokaryotes and seven eukaryotes for which the largest numbers of coding sequences are available in data banks. In prokaryotes, both narrow and broad distributions were found. In eukaryotes, distributions were very broad (except for Saccharomyces cerevisiae) and remarkably different for different genomes. In low-GC genomes, third codon positions were lower in GC than first + second codon positions and trailed towards high GC; the opposite situation was found for high-GC genomes. In all genomes, first codon positions were higher in GC than second codon positions. We then investigated the compositional correlations between third and first + second codon positions in prokaryotic genomes (the 16 mentioned above plus 87 additional ones) and in genome compartments of eukaryotes. A general, common relationship was found, which also holds within the same (heterogeneous) genomes. This universal correlation is due to the fact that the relative effects of compositional constraints on different codon positions are the same, on the average, whatever the genome under consideration.     

Genome compaction and stability in microsporidian intracellular parasites

Microsporidian genomes are extraordinary among eukaryotes for their extreme reduction: although they are similar in form to other eukaryotic genomes, they are typically smaller than many prokaryotic genomes. At the same time, their rates of sequence evolution are among the highest for eukaryotic organisms. To explore the effects of compaction on nuclear genome evolution, we sequenced 685,000 bp of the Antonospora locustae genome (formerly Nosema locustae) and compared its organization with the recently completed genome of the human parasite Encephalitozoon cuniculi. Despite being very distantly related, the genomes of these two microsporidian species have retained an unexpected degree of synteny: 13% of genes are in the same context, and 30% of the genes were separated by a small number of short rearrangements. Microsporidian genomes are, therefore, paradoxically composed of rapidly evolving sequences harbored within a slowly evolving genome, although these two processes are sometimes considered to be coupled. Microsporidian genomes show that eukaryotic genomes (like genes) do not evolve in a clock-like fashion, and genome stability may result from compaction in addition to a lack of recombination, as has been traditionally thought to occur in bacterial and organelle genomes.     

A unique ATG triplet downstream of gene start in archaea: implications for translation initiation and evolution

Searching for unique features of archaeal genome may shed light on the mechanism of gene regulation in primitive life forms. Statistical analysis of ATG frequency on the complete genome sequences of 16 archaea, 20 bacteria and 2 eukaryotes revealed that most of the archaeal genomes have a remarkably high ATG frequency at the position of nine nucleotide (nt) downstream of the translation initiation site (the first nucleotide of the translation initiation codon is designated as 0). To understand the role of this unique ATG in archaea, we further analyzed the ATG-initiated genes and non-ATG-initiated genes separately, and the results indicated that only the non-ATG-initiated genes contribute to the high ATG frequency at position +9. This led us to speculate that the in-frame ATG at +9 may serve as a remedial initiation site for archaea in case of initiation failure at the regular site. In addition, it seems that this phenomenon does not result from the harsh environment that archaea are usually viable according to the fact that no considerably high ATG frequency at +9 was observed in all the four thermophilic bacteria that also live in harsh environment. We proposed that the high ATG frequency at position +9 might reflect the decreased efficiency of the translation initiation machinery in archaea. Since archaea evolve very slowly, this unique characteristic of high ATG frequency at position +9 may present the primitive state of the Universal Ancestor.     

Prokaryotic structural maintenance of chromosomes (SMC) proteins: distribution, phylogeny, and comparison with MukBs and additional prokaryotic and eukaryotic coiled-coil proteins.

Structural maintenance of chromosomes (SMC) proteins are known to be essential for chromosome segregation in some prokaryotes and in eukaryotes. A systematic search for the distribution of SMC proteins in prokaryotes with fully or partially sequenced genomes showed that they form a larger family than previously anticipated and raised the number of known prokaryotic homologs to 54. Secondary structure predictions revealed that the length of the globular N-terminal and C-terminal domains is extremely well conserved in contrast to the hinge domain and coiled-coil domains which are considerably shorter in several bacterial species. SMC proteins are present in all gram-positive bacteria and in nearly all archaea while they were found in less than half of the gram-negative bacteria. Phylogenetic analyses indicate that the SMC tree roughly resembles the 16S rRNA tree, but that cyanobacteria and Aquifex aeolicus obtained smc genes by lateral transfer from archaea. Fourteen out of 22 smc genes located in fully sequenced genomes seem to be co-transcribed with a second gene out of six different gene families, indicating that the deduced gene products might be involved in similar functions. The SMC proteins were compared with other prokaryotic proteins with long coiled-coil domains. The lengths of different protein domains and signature sequences allowed to differentiate SMCs, MukBs, which were found to be confined to gamma proteobacteria, and two subfamilies of COG 0419 including the SbcC nuclease from E. coli. A phylogenetic analysis was performed including the prokaryotic coiled-coil proteins as well as SMCs and Rad18 proteins from selected eukaryotes.     

Microbial ecology in the age of genomics and metagenomics: concepts, tools, and recent advances

Microbial ecology examines the diversity and activity of micro-organisms in Earth's biosphere. In the last 20 years, the application of genomics tools have revolutionized microbial ecological studies and drastically expanded our view on the previously underappreciated microbial world. This review first introduces the basic concepts in microbial ecology and the main genomics methods that have been used to examine natural microbial populations and communities. In the ensuing three specific sections, the applications of the genomics in microbial ecological research are highlighted. The first describes the widespread application of multilocus sequence typing and representational difference analysis in studying genetic variation within microbial species. Such investigations have identified that migration, horizontal gene transfer and recombination are common in natural microbial populations and that microbial strains can be highly variable in genome size and gene content. The second section highlights and summarizes the use of four specific genomics methods (phylogenetic analysis of ribosomal RNA, DNA-DNA re-association kinetics, metagenomics, and micro-arrays) in analysing the diversity and potential activity of microbial populations and communities from a variety of terrestrial and aquatic environments. Such analyses have identified many unexpected phylogenetic lineages in viruses, bacteria, archaea, and microbial eukaryotes. Functional analyses of environmental DNA also revealed highly prevalent, but previously unknown, metabolic processes in natural microbial communities. In the third section, the ecological implications of sequenced microbial genomes are briefly discussed. Comparative analyses of prokaryotic genomic sequences suggest the importance of ecology in determining microbial genome size and gene content. The significant variability in genome size and gene content among strains and species of prokaryotes indicate the highly fluid nature of prokaryotic genomes, a result consistent with those from multilocus sequence typing and representational difference analyses. The integration of various levels of ecological analyses coupled to the application and further development of high throughput technologies are accelerating the pace of discovery in microbial ecology.     

Ab initio gene identification in metagenomic sequences

We describe an algorithm for gene identification in DNA sequences derived from shotgun sequencing of microbial communities. Accurate ab initio gene prediction in a short nucleotide sequence of anonymous origin is hampered by uncertainty in model parameters. While several machine learning approaches could be proposed to bypass this difficulty, one effective method is to estimate parameters from dependencies, formed in evolution, between frequencies of oligonucleotides in protein-coding regions and genome nucleotide composition. Original version of the method was proposed in 1999 and has been used since for (i) reconstructing codon frequency vector needed for gene finding in viral genomes and (ii) initializing parameters of self-training gene finding algorithms. With advent of new prokaryotic genomes en masse it became possible to enhance the original approach by using direct polynomial and logistic approximations of oligonucleotide frequencies, as well as by separating models for bacteria and archaea. These advances have increased the accuracy of model reconstruction and, subsequently, gene prediction. We describe the refined method and assess its accuracy on known prokaryotic genomes split into short sequences. Also, we show that as a result of application of the new method, several thousands of new genes could be added to existing annotations of several human and mouse gut metagenomes.     

Clusters of orthologous genes for 41 archaeal genomes and implications for evolutionary genomics of archaea

Background  

An evolutionary classification of genes from sequenced genomes that distinguishes between orthologs and paralogs is indispensable for genome annotation and evolutionary reconstruction. Shortly after multiple genome sequences of bacteria, archaea, and unicellular eukaryotes became available, an attempt on such a classification was implemented in Clusters of Orthologous Groups of proteins (COGs). Rapid accumulation of genome sequences creates opportunities for refining COGs but also represents a challenge because of error amplification. One of the practical strategies involves construction of refined COGs for phylogenetically compact subsets of genomes.     

Argonaute Proteins and Mechanisms of RNA Interference in Eukaryotes and Prokaryotes

Noncoding RNAs play essential roles in genetic regulation in all organisms. In eukaryotic cells, many small non-coding RNAs act in complex with Argonaute proteins and regulate gene expression by recognizing complementary RNA targets. The complexes of Argonaute proteins with small RNAs also play a key role in silencing of mobile genetic elements and, in some cases, viruses. These processes are collectively called RNA interference. RNA interference is a powerful tool for specific gene silencing in both basic research and therapeutic applications. Argonaute proteins are also found in prokaryotic organisms. Recent studies have shown that prokaryotic Argonautes can also cleave their target nucleic acids, in particular DNA. This activity of prokaryotic Argonautes might potentially be used to edit eukaryotic genomes. However, the molecular mechanisms of small nucleic acid biogenesis and the functions of Argonaute proteins, in particular in bacteria and archaea, remain largely unknown. Here we briefly review available data on the RNA interference processes and Argonaute proteins in eukaryotes and prokaryotes.     

Analysis of evolution of exon-intron structure of eukaryotic genes

The availability of multiple, complete eukaryotic genome sequences allows one to address many fundamental evolutionary questions on genome scale. One such important, long-standing problem is evolution of exon-intron structure of eukaryotic genes. Analysis of orthologous genes from completely sequenced genomes revealed numerous shared intron positions in orthologous genes from animals and plants and even between animals, plants and protists. The data on shared and lineage-specific intron positions were used as the starting point for evolutionary reconstruction with parsimony and maximum-likelihood approaches. Parsimony methods produce reconstructions with intron-rich ancestors but also infer lineage-specific, in many cases, high levels of intron loss and gain. Different probabilistic models gave opposite results, apparently depending on model parameters and assumptions, from domination of intron loss, with extremely intron-rich ancestors, to dramatic excess of gains, to the point of denying any true conservation of intron positions among deep eukaryotic lineages. Development of models with adequate, realistic parameters and assumptions seems to be crucial for obtaining more definitive estimates of intron gain and loss in different eukaryotic lineages. Many shared intron positions were detected in ancestral eukaryotic paralogues which evolved by duplication prior to the divergence of extant eukaryotic lineages. These findings indicate that numerous introns were present in eukaryotic genes already at the earliest stages of evolution of eukaryotes and are compatible with the hypothesis that the original, catastrophic intron invasion accompanied the emergence of the eukaryotic cells. Comparison of various features of old and younger introns starts shedding light on probable mechanisms of intron insertion, indicating that propagation of old introns is unlikely to be a major mechanism for origin of new ones. The existence and structure of ancestral protosplice sites were addressed by examining the context of introns inserted within codons that encode amino acids conserved in all eukaryotes and, accordingly, are not subject to selection for splicing efficiency. It was shown that introns indeed predominantly insert into or are fixed in specific protosplice sites which have the consensus sequence (A/C)AG|Gt.     

The slow road to the eukaryotic genome

The eukaryotic genome is a mosaic of eubacterial and archaeal genes in addition to those unique to itself. The mosaic may have arisen as the result of two prokaryotes merging their genomes, or from genes acquired from an endosymbiont of eubacterial origin. A third possibility is that the eukaryotic genome arose from successive events of lateral gene transfer over long periods of time. This theory does not exclude the endosymbiont, but questions whether it is necessary to explain the peculiar set of eukaryotic genes. We use phylogenetic studies and reconstructions of ancestral first appearances of genes on the prokaryotic phylogeny to assess evidence for the lateral gene transfer scenario. We find that phylogenies advanced to support fusion can also arise from a succession of lateral gene transfer events. Our reconstructions of ancestral first appearances of genes reveal that the various genes that make up the eukaryotic mosaic arose at different times and in diverse lineages on the prokaryotic tree, and were not available in a single lineage. Successive events of lateral gene transfer can explain the unusual mosaic structure of the eukaryotic genome, with its content linked to the immediate adaptive value of the genes its acquired. Progress in understanding eukaryotes may come from identifying ancestral features such as the eukaryotic splicesome that could explain why this lineage invaded, or created, the eukaryotic niche.