Genomic and phylogenetic perspectives on the evolution of prokaryotes

Comparative genomics analysis and transposon mutagenesis provides new insights into high menaquinone-7 biosynthetic potential of Bacillus subtilis natto

This research combined Whole-Genome sequencing, intraspecific comparative genomics and transposon mutagenesis to investigate the menaquinone-7 (MK-7) synthesis potential in Bacillus subtilis natto. First, Whole-Genome sequencing showed that Bacillus subtilis natto BN-P15-11-1 contains one single circular chromosome in size of 3,982,436 bp with a GC content of 43.85 %, harboring 4,053 predicted coding genes. Next, the comparative genomics analysis among strain BN-P15-11-1 with model Bacillus subtilis 168 and four typical Bacillus subtilis natto strains proves that the closer evolutionary relationship Bacillus subtilis natto BN-P15-11-1 and Bacillus subtilis 168 both exhibit strong biosynthetic potential. To further dig for MK-7 biosynthesis latent capacity of BN-P15-11-1, we constructed a mutant library using transposons and a high throughput screening method using microplates. We obtained a YqgQ deficient high MK-7 yield strain F4 with a yield 3.02 times that of the parent strain. Experiments also showed that the high yield mutants had defects in different transcription and translation regulatory factor genes, indicating that regulatory factor defects may affect the biosynthesis and accumulation of MK-7 by altering the overall metabolic level. The findings of this study will provide more novel insights on the precise identification and rational utilization of the Bacillus subtilis subspecies for biosynthesis latent capacity.

Amino Acid Specificity of Ancestral Aminoacyl-tRNA Synthetase Prior to the Last Universal Common Ancestor Commonote commonote

Extant organisms commonly use 20 amino acids in protein synthesis. In the translation system, aminoacyl-tRNA synthetase (ARS) selectively binds an amino acid and transfers it to the cognate tRNA. It is postulated that the amino acid repertoire of ARS expanded during the development of the translation system. In this study we generated composite phylogenetic trees for seven ARSs (SerRS, ProRS, ThrRS, GlyRS-1, HisRS, AspRS, and LysRS) which are thought to have diverged by gene duplication followed by mutation, before the evolution of the last universal common ancestor. The composite phylogenetic tree shows that the AspRS/LysRS branch diverged from the other five ARSs at the deepest node, with the GlyRS/HisRS branch and the other three ARSs (ThrRS, ProRS and SerRS) diverging at the second deepest node. ThrRS diverged next, and finally ProRS and SerRS diverged from each other. Based on the phylogenetic tree, sequences of the ancestral ARSs prior to the evolution of the last universal common ancestor were predicted. The amino acid specificity of each ancestral ARS was then postulated by comparison with amino acid recognition sites of ARSs of extant organisms. Our predictions demonstrate that ancestral ARSs had substantial specificity and that the number of amino acid types amino-acylated by proteinaceous ARSs was limited before the appearance of a fuller range of proteinaceous ARS species. From an assumption that 10 amino acid species are required for folding and function, proteinaceous ARS possibly evolved in a translation system composed of preexisting ribozyme ARSs, before the evolution of the last universal common ancestor.

Characterizing the amino acid activation center of the naturally editing-deficient aminoacyl-tRNA synthetase PheRS in Mycoplasma mobile

To ensure correct amino acids are incorporated during protein synthesis, aminoacyl-tRNA synthetases (aaRSs) employ proofreading mechanisms collectively referred to as editing. Although editing is important for viability, editing-deficient aaRSs have been identified in host-dependent organisms. In Mycoplasma mobile, editing-deficient PheRS and LeuRS have been identified. We characterized the amino acid activation site of MmPheRS and identified a previously unknown hyperaccurate mutation, L287F. Additionally, we report that m-Tyr, an oxidation byproduct of Phe which is toxic to editing-deficient cells, is poorly discriminated by MmPheRS activation and is not subjected to editing. Furthermore, expressing MmPheRS and the hyperaccurate variants renders Escherichia coli susceptible to m-Tyr stress, indicating that active site discrimination is insufficient in tolerating excess m-Tyr.

Quest for Ancestors of Eukaryal Cells Based on Phylogenetic Analyses of Aminoacyl-tRNA Synthetases

The three-domain phylogenetic system of life has been challenged, particularly with regard to the position of Eukarya. The recent increase of known genome sequences has allowed phylogenetic analyses of all extant organisms using concatenated sequence alignment of universally conserved genes; these data supported the two-domain hypothesis, which place eukaryal species as ingroups of the Domain Archaea. However, the origin of Eukarya is complicated: the closest archaeal species to Eukarya differs in single-gene phylogenetic analyses depending on the genes. In this report, we performed molecular phylogenetic analyses of 23 aminoacyl-tRNA synthetases (ARS). Cytoplasmic ARSs in 12 trees showed a monophyletic Eukaryotic branch. One ARS originated from TACK superphylum. One ARS originated from Euryarchaeota and three originated from DPANN superphylum. Four ARSs originated from different bacterial species. The other 8 cytoplasmic ARSs were split into two or three groups in respective trees, which suggested that the cytoplasmic ARSs were replaced by secondary ARSs, and the original ARSs have been lost during evolution of Eukarya. In these trees, one original cytoplasmic ARS was derived from Euryarchaeota and three were derived from DPANN superphylum. Our results strongly support the two-domain hypothesis. We discovered that rampant-independent lateral gene transfers from several archaeal species of DPANN superphylum have contributed to the formation of Eukaryal cells. Based on our phylogenetic analyses, we proposed a model for the establishment of Eukarya.

Evolution of double-stranded DNA viruses of eukaryotes: from bacteriophages to transposons to giant viruses

Diverse eukaryotes including animals and protists are hosts to a broad variety of viruses with double-stranded (ds) DNA genomes, from the largest known viruses, such as pandoraviruses and mimiviruses, to tiny polyomaviruses. Recent comparative genomic analyses have revealed many evolutionary connections between dsDNA viruses of eukaryotes, bacteriophages, transposable elements, and linear DNA plasmids. These findings provide an evolutionary scenario that derives several major groups of eukaryotic dsDNA viruses, including the proposed order “Megavirales,” adenoviruses, and virophages from a group of large virus-like transposons known as Polintons (Mavericks). The Polintons have been recently shown to encode two capsid proteins, suggesting that these elements lead a dual lifestyle with both a transposon and a viral phase and should perhaps more appropriately be named polintoviruses. Here, we describe the recently identified evolutionary relationships between bacteriophages of the family Tectiviridae, polintoviruses, adenoviruses, virophages, large and giant DNA viruses of eukaryotes of the proposed order “Megavirales,” and linear mitochondrial and cytoplasmic plasmids. We outline an evolutionary scenario under which the polintoviruses were the first group of eukaryotic dsDNA viruses that evolved from bacteriophages and became the ancestors of most large DNA viruses of eukaryotes and a variety of other selfish elements. Distinct lines of origin are detectable only for herpesviruses (from a different bacteriophage root) and polyoma/papillomaviruses (from single-stranded DNA viruses and ultimately from plasmids). Phylogenomic analysis of giant viruses provides compelling evidence of their independent origins from smaller members of the putative order “Megavirales,” refuting the speculations on the evolution of these viruses from an extinct fourth domain of cellular life.

Lactococcus lactis TrxD represents a subgroup of thioredoxins prevalent in Gram-positive bacteria containing WCXDC active site motifs

Three protein disulfide reductases of the thioredoxin superfamily from the industrially important Gram-positive Lactococcus lactis (LlTrxA, LlTrxD and LlNrdH) are compared to the "classical" thioredoxin from Escherichia coli (EcTrx1). LlTrxA resembles EcTrx1 with a WCGPC active site motif and other key residues conserved. By contrast, LlTrxD is more distantly related and contains a WCGDC motif. Bioinformatics analysis suggests that LlTrxD represents a subgroup of thioredoxins from Gram-positive bacteria. LlNrdH is a glutaredoxin-like electron donor for ribonucleotide reductase class Ib. Based on protein-protein equilibria LlTrxA (E°'=-259mV) and LlNrdH (E°'=-238mV) show approximately 10mV higher standard state redox potentials than the corresponding E. coli homologues, while E°' of LlTrxD is -243mV, more similar to glutaredoxin than "classical" thioredoxin. EcTrx1 and LlTrxA have high capacity to reduce insulin disulfides and their exposed active site thiol is alkylated at a similar rate at pH 7.0. LlTrxD on the other hand, is alkylated by iodoacetamide at almost 100 fold higher rate and shows no activity towards insulin disulfides. LlTrxA, LlTrxD and LlNrdH are all efficiently reduced by NADPH dependent thioredoxin reductase (TrxR) from L. lactis and good cross-reactivity towards E. coli TrxR was observed with LlTrxD as the notable exception.

Analysis of the population structure of Anaplasma phagocytophilum using multilocus sequence typing

Anaplasma phagocytophilum is a Gram-negative obligate intracellular bacterium that replicates in neutrophils. It is transmitted via tick-bite and causes febrile disease in humans and animals. Human granulocytic anaplasmosis is regarded as an emerging infectious disease in North America, Europe and Asia. However, although increasingly detected, it is still rare in Europe. Clinically apparent A. phagocytophilum infections in animals are mainly found in horses, dogs, cats, sheep and cattle. Evidence from cross-infection experiments that A. phagocytophilum isolates of distinct host origin are not uniformly infectious for heterologous hosts has led to several approaches of molecular strain characterization. Unfortunately, the results of these studies are not always easily comparable, because different gene regions and fragment lengths were investigated. Multilocus sequence typing is a widely accepted method for molecular characterization of bacteria. We here provide for the first time a universal typing method that is easily transferable between different laboratories. We validated our approach on an unprecedented large data set of almost 400 A. phagocytophilum strains from humans and animals mostly from Europe. The typability was 74% (284/383). One major clonal complex containing 177 strains was detected. However, 54% (49/90) of the sequence types were not part of a clonal complex indicating that the population structure of A. phagocytophilum is probably semiclonal. All strains from humans, dogs and horses from Europe belonged to the same clonal complex. As canine and equine granulocytic anaplasmosis occurs frequently in Europe, human granulocytic anaplasmosis is likely to be underdiagnosed in Europe. Further, wild boars and hedgehogs may serve as reservoir hosts of the disease in humans and domestic animals in Europe, because their strains belonged to the same clonal complex. In contrast, as they were only distantly related, roe deer, voles and shrews are unlikely to harbor A. phagocytophilum strains infectious for humans, domestic or farm animals.

The anguibactin biosynthesis and transport genes are encoded in the chromosome of Vibrio harveyi: a possible evolutionary origin for the pJM1 plasmid-encoded system of Vibrio anguillarum?

Many Vibrio anguillarum serotype O1 strains carry 65-kb pJM1-type plasmids harboring genes involved in siderophore anguibactin biosynthesis and transport. The anguibactin system is an essential factor for V. anguillarum to survive under iron-limiting conditions, and as a consequence, it is a very important virulence factor of this bacterium. Our comparative analysis of genomic data identified a cluster harboring homologs of anguibactin biosynthesis and transport genes in the chromosome of Vibrio harveyi. We have purified the putative anguibactin siderophore and demonstrated that it is indeed anguibactin by mass spectrometry and specific bioassays. Furthermore, we characterized two genes, angR and fatA, in this chromosome cluster that, respectively, participate in anguibactin biosynthesis and transport as determined by mutagenesis analysis. Furthermore, we found that the V. harveyi FatA protein is located in the outer membrane fractions as previously demonstrated in V. anguillarum. Based on our data, we propose that the anguibactin biosynthesis and transport cluster in the V. anguillarum pJM1 plasmid have likely evolved from the chromosome cluster of V. harveyi or vice versa.

Progress in genomics, metabolism and biotechnology of bifidobacteria

Members of the genus Bifidobacterium were first described over a century ago and were quickly associated with a healthy intestinal tract due to their numerical dominance in breast-fed babies as compared to bottle-fed infants. Health benefits elicited by bifidobacteria to its host, as supported by clinical trials, have led to their wide application as probiotic components of health-promoting foods, especially in fermented dairy products. However, the relative paucity of genetic tools available for bifidobacteria has impeded development of a comprehensive molecular understanding of this genus. In this review we present a summary of current knowledge on bifidobacterial metabolism, classification, physiology and genetics and outline the currently available methods for genetically accessing and manipulating the genus.

Darwinian evolution in the light of genomics

Comparative genomics and systems biology offer unprecedented opportunities for testing central tenets of evolutionary biology formulated by Darwin in the Origin of Species in 1859 and expanded in the Modern Synthesis 100 years later. Evolutionary-genomic studies show that natural selection is only one of the forces that shape genome evolution and is not quantitatively dominant, whereas non-adaptive processes are much more prominent than previously suspected. Major contributions of horizontal gene transfer and diverse selfish genetic elements to genome evolution undermine the Tree of Life concept. An adequate depiction of evolution requires the more complex concept of a network or ‘forest’ of life. There is no consistent tendency of evolution towards increased genomic complexity, and when complexity increases, this appears to be a non-adaptive consequence of evolution under weak purifying selection rather than an adaptation. Several universals of genome evolution were discovered including the invariant distributions of evolutionary rates among orthologous genes from diverse genomes and of paralogous gene family sizes, and the negative correlation between gene expression level and sequence evolution rate. Simple, non-adaptive models of evolution explain some of these universals, suggesting that a new synthesis of evolutionary biology might become feasible in a not so remote future.

Origin and evolution of the mitochondrial aminoacyl-tRNA synthetases

Many theories favor a fusion of 2 prokaryotic genomes for the origin of the Eukaryotes, but there are disagreements on the origin, timing, and cellular structures of the cells involved. Equally controversial is the source of the nuclear genes for mitochondrial proteins, although the alpha-proteobacterial contribution to the mitochondrial genome is well established. Phylogenetic inferences show that the nuclearly encoded mitochondrial aminoacyl-tRNA synthetases (aaRSs) occupy a position in the tree that is not close to any of the currently sequenced alpha-proteobacterial genomes, despite cohesive and remarkably well-resolved alpha-proteobacterial clades in 12 of the 20 trees. Two or more alpha-proteobacterial clusters were observed in 8 cases, indicative of differential loss of paralogous genes or horizontal gene transfer. Replacement and retargeting events within the nuclear genomes of the Eukaryotes was indicated in 10 trees, 4 of which also show split alpha-proteobacterial groups. A majority of the mitochondrial aaRSs originate from within the bacterial domain, but none specifically from the alpha-Proteobacteria. For some aaRS, the endosymbiotic origin may have been erased by ongoing gene replacements on the bacterial as well as the eukaryotic side. For others that accurately resolve the alpha-proteobacterial divergence patterns, the lack of affiliation with mitochondria is more surprising. We hypothesize that the ancestral eukaryotic gene pool hosted primordial "bacterial-like" genes, to which a limited set of alpha-proteobacterial genes, mostly coding for components of the respiratory chain complexes, were added and selectively maintained.

Optimal gene partition into operons correlates with gene functional order

Gene arrangement into operons varies between bacterial species. Genes in a given system can be on one operon in some organisms and on several operons in other organisms. Existing theories explain why genes that work together should be on the same operon, since this allows for advantageous lateral gene transfer and accurate stoichiometry. But what causes the frequent separation into multiple operons of co-regulated genes that act together in a pathway? Here we suggest that separation is due to benefits made possible by differential regulation of each operon. We present a simple mathematical model for the optimal distribution of genes into operons based on a balance of the cost of operons and the benefit of regulation that provides 'just-when-needed' temporal order. The analysis predicts that genes are arranged such that genes on the same operon do not skip functional steps in the pathway. This prediction is supported by genomic data from 137 bacterial genomes. Our work suggests that gene arrangement is not only the result of random historical drift, genome re-arrangement and gene transfer, but has elements that are solutions of an evolutionary optimization problem. Thus gene functional order may be inferred by analyzing the operon structure across different genomes.

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.

Loss of editing activity during the evolution of mitochondrial phenylalanyl-tRNA synthetase

Accurate selection of amino acids is essential for faithful translation of the genetic code. Errors during amino acid selection are usually corrected by the editing activity of aminoacyl-tRNA synthetases such as phenylalanyl-tRNA synthetases (PheRS), which edit misactivated tyrosine. Comparison of cytosolic and mitochondrial PheRS from the yeast Saccharomyces cerevisiae suggested that the organellar protein might lack the editing activity. Yeast cytosolic PheRS was found to contain an editing site, which upon disruption abolished both cis and trans editing of Tyr-tRNA(Phe). Wild-type mitochondrial PheRS lacked cis and trans editing and could synthesize Tyr-tRNA(Phe), an activity enhanced in active site variants with improved tyrosine recognition. Possible trans editing was investigated in isolated mitochondrial extracts, but no such activity was detected. These data indicate that the mitochondrial protein synthesis machinery lacks the tyrosine proofreading activity characteristic of cytosolic translation. This difference between the mitochondria and the cytosol suggests that either organellar protein synthesis quality control is focused on another step or that translation in this compartment is inherently less accurate than in the cytosol.

Does a tree-like phylogeny only exist at the tips in the prokaryotes?

The extent to which prokaryotic evolution has been influenced by horizontal gene transfer (HGT) and therefore might be more of a network than a tree is unclear. Here we use supertree methods to ask whether a definitive prokaryotic phylogenetic tree exists and whether it can be confidently inferred using orthologous genes. We analysed an 11-taxon dataset spanning the deepest divisions of prokaryotic relationships, a 10-taxon dataset spanning the relatively recent gamma-proteobacteria and a 61-taxon dataset spanning both, using species for which complete genomes are available. Congruence among gene trees spanning deep relationships is not better than random. By contrast, a strong, almost perfect phylogenetic signal exists in gamma-proteobacterial genes. Deep-level prokaryotic relationships are difficult to infer because of signal erosion, systematic bias, hidden paralogy and/or HGT. Our results do not preclude levels of HGT that would be inconsistent with the notion of a prokaryotic phylogeny. This approach will help decide the extent to which we can say that there is a prokaryotic phylogeny and where in the phylogeny a cohesive genomic signal exists.

Paradigm change in evolutionary microbiology

Thomas Kuhn had little to say about scientific change in biological science, and biologists are ambivalent about how applicable his framework is for their disciplines. We apply Kuhn's account of paradigm change to evolutionary microbiology, where key Darwinian tenets are being challenged by two decades of findings from molecular phylogenetics. The chief culprit is lateral gene transfer, which undermines the role of vertical descent and the representation of evolutionary history as a tree of life. To assess Kuhn's relevance to this controversy, we add a social analysis of the scientists involved to the historical and philosophical debates. We conclude that while Kuhn's account may capture aspects of the pattern (or outcome) of an episode of scientific change, he has little to say about how the process of generating new understandings is occurring in evolutionary microbiology. Once Kuhn's application is limited to that of an initial investigative probe into how scientific problem-solving occurs, his disciplinary scope becomes broader.

Post-transfer editing in vitro and in vivo by the beta subunit of phenylalanyl-tRNA synthetase

Translation of the genetic code requires attachment of tRNAs to their cognate amino acids. Errors during amino-acid activation and tRNA esterification are corrected by aminoacyl-tRNA synthetase-catalyzed editing reactions, as extensively described for aliphatic amino acids. The contribution of editing to aromatic amino-acid discrimination is less well understood. We show that phenylalanyl-tRNA synthetase misactivates tyrosine and that it subsequently corrects such errors through hydrolysis of tyrosyl-adenylate and Tyr-tRNA(Phe). Structural modeling combined with an in vivo genetic screen identified the editing site in the B3/B4 domain of the beta subunit, 40 angstroms from the active site in the alpha subunit. Replacements of residues within the editing site had no effect on Phe-tRNA(Phe) synthesis, but abolished hydrolysis of Tyr-tRNA(Phe) in vitro. Expression of the corresponding mutants in Escherichia coli significantly slowed growth, and changed the activity of a recoded beta-galactosidase variant by misincorporating tyrosine in place of phenylalanine. This loss in aromatic amino-acid discrimination in vivo revealed that editing by phenylalanyl-tRNA synthetase is essential for faithful translation of the genetic code.

Gamma carbonic anhydrases in plant mitochondria

Three genes from Arabidopsis thaliana with high sequence similarity to gamma carbonic anhydrase (gammaCA), a Zn containing enzyme from Methanosarcina thermophila (CAM), were identified and characterized. Evolutionary and structural analyses predict that these genes code for active forms of gammaCA. Phylogenetic analyses reveal that these Arabidopsis gene products cluster together with CAM and related sequences from alpha and gamma proteobacteria, organisms proposed as the mitochondrial endosymbiont ancestor. Indeed, in vitro and in vivo experiments indicate that these gene products are transported into the mitochondria as occurs with several mitochondrial protein genes transferred, during evolution, from the endosymbiotic bacteria to the host genome. Moreover, putative CAM orthologous genes are detected in other plants and green algae and were predicted to be imported to mitochondria. Structural modeling and sequence analysis performed in more than a hundred homologous sequences show a high conservation of functionally important active site residues. Thus, the three histidine residues involved in Zn coordination (His 81, 117 and 122), Arg 59, Asp 61, Gin 75, and Asp 76 of CAM are conserved and properly arranged in the active site cavity of the models. Two other functionally important residues (Glu 62 and Glu 84 of CAM) are lacking, but alternative amino acids that might serve to their roles are postulated. Accordingly, we propose that photosynthetic eukaryotic organisms (green algae and plants) contain gammaCAs and that these enzymes codified by nuclear genes are imported into mitochondria to accomplish their biological function.

Conserved Ser residues, the shutter region, and speciation in serpin evolution

The suicide inhibitory mechanism of serine protease inhibitors of the serpin superfamily depends heavily on their structural flexibility, which is controlled in large part by the breach and shutter regions of the central Abeta-sheet. We examined codon usage by the highly conserved residues, Ser-53 and Ser-56, of the shutter region and found a TCN-AGY usage dichotomy for Ser-56 that remarkably is linked to the protostome-deuterostome split. Our results suggest that serpin evolution was driven by phylogenetic speciation and not pressure to fulfill new physiologic functions mitigating against coevolution with the family of serine proteases they inhibit.

Exploring prokaryotic diversity in the genomic era

Our understanding of prokaryote biology from study of pure cultures and genome sequencing has been limited by a pronounced sampling bias towards four bacterial phyla - Proteobacteria, Firmicutes, Actinobacteria and Bacteroidetes - out of 35 bacterial and 18 archaeal phylum-level lineages. This bias is beginning to be rectified by the use of phylogenetically directed isolation strategies and by directly accessing microbial genomes from environmental samples.