Genomics of bacteria and archaea: the emerging dynamic view of the prokaryotic world

Investigations into the Biosynthesis, Regulation, and Self-Resistance of Toxoflavin in Pseudomonas protegens Pf-5

Pseudomonas spp. are prolific producers of natural products from many structural classes. Here we show that the soil bacterium Pseudomonas protegens Pf-5 is capable of producing trace levels of the triazine natural product toxoflavin (1) under microaerobic conditions. We evaluated toxoflavin production by derivatives of Pf-5 with deletions in specific biosynthesis genes, which led us to propose a revised biosynthetic pathway for toxoflavin that shares the first two steps with riboflavin biosynthesis. We also report that toxM, which is not present in the well-characterized cluster of Burkholderia glumae, encodes a monooxygenase that degrades toxoflavin. The toxoflavin degradation product of ToxM is identical to that of TflA, the toxoflavin lyase from Paenibacillus polymyxa. Toxoflavin production by P. protegens causes inhibition of several plant-pathogenic bacteria, and introduction of toxM into the toxoflavin-sensitive strain Pseudomonas syringae DC3000 results in resistance to toxoflavin.

Evolutionary Advantage Conferred by an Eukaryote-to-Eukaryote Gene Transfer Event in Wine Yeasts

Although an increasing number of horizontal gene transfers have been reported in eukaryotes, experimental evidence for their adaptive value is lacking. Here, we report the recent transfer of a 158-kb genomic region between Torulaspora microellipsoides and Saccharomyces cerevisiae wine yeasts or closely related strains. This genomic region has undergone several rearrangements in S. cerevisiae strains, including gene loss and gene conversion between two tandemly duplicated FOT genes encoding oligopeptide transporters. We show that FOT genes confer a strong competitive advantage during grape must fermentation by increasing the number and diversity of oligopeptides that yeast can utilize as a source of nitrogen, thereby improving biomass formation, fermentation efficiency, and cell viability. Thus, the acquisition of FOT genes has favored yeast adaptation to the nitrogen-limited wine fermentation environment. This finding indicates that anthropic environments offer substantial ecological opportunity for evolutionary diversification through gene exchange between distant yeast species.

Development of New Modular Genetic Tools for Engineering the Halophilic Archaeon Halobacterium salinarum

Our ability to genetically manipulate living organisms is usually constrained by the efficiency of the genetic tools available for the system of interest. In this report, we present the design, construction and characterization of a set of four new modular vectors, the pHsal series, for engineering Halobacterium salinarum, a model halophilic archaeon widely used in systems biology studies. The pHsal shuttle vectors are organized in four modules: (i) the E. coli’s specific part, containing a ColE1 origin of replication and an ampicillin resistance marker, (ii) the resistance marker and (iii) the replication origin, which are specific to H. salinarum and (iv) the cargo, which will carry a sequence of interest cloned in a multiple cloning site, flanked by universal M13 primers. Each module was constructed using only minimal functional elements that were sequence edited to eliminate redundant restriction sites useful for cloning. This optimization process allowed the construction of vectors with reduced sizes compared to currently available platforms and expanded multiple cloning sites. Additionally, the strong constitutive promoter of the fer2 gene was sequence optimized and incorporated into the platform to allow high-level expression of heterologous genes in H. salinarum. The system also includes a new minimal suicide vector for the generation of knockouts and/or the incorporation of chromosomal tags, as well as a vector for promoter probing using a GFP gene as reporter. This new set of optimized vectors should strongly facilitate the engineering of H. salinarum and similar strategies could be implemented for other archaea.

Horizontal gene transfer and genome evolution in Methanosarcina

Background

Genomes of Methanosarcina spp. are among the largest archaeal genomes. One suggested reason for that is massive horizontal gene transfer (HGT) from bacteria. Genes of bacterial origin may be involved in the central metabolism and solute transport, in particular sugar synthesis, sulfur metabolism, phosphate metabolism, DNA repair, transport of small molecules etc. Horizontally transferred (HT) genes are considered to play the key role in the ability of Methanosarcina spp. to inhabit diverse environments. At the moment, genomes of three Methanosarcina spp. have been sequenced, and while these genomes vary in length and number of protein-coding genes, they all have been shown to accumulate HT genes. However, previous estimates had been made when fewer archaeal genomes were known. Moreover, several Methanosarcinaceae genomes from other genera have been sequenced recently. Here, we revise the census of genes of bacterial origin in Methanosarcinaceae.

Results

About 5 % of Methanosarcina genes have been shown to be horizontally transferred from various bacterial groups to the last common ancestor either of Methanosarcinaceae, or Methanosarcina, or later in the evolution. Simulation of the composition of the NCBI protein non-redundant database for different years demonstrates that the estimates of the HGT rate have decreased drastically since 2002, the year of publication of the first Methanosarcina genome.

The phylogenetic distribution of HT gene donors is non-uniform. Most HT genes were transferred from Firmicutes and Proteobacteria, while no HGT events from Actinobacteria to the common ancestor of Methanosarcinaceae were found. About 50 % of HT genes are involved in metabolism. Horizontal transfer of transcription factors is not common, while 46 % of horizontally transferred genes have demonstrated differential expression in a variety of conditions. HGT of complete operons is relatively infrequent and half of HT genes do not belong to operons.

Conclusions

While genes of bacterial origin are still more frequent in Methanosarcinaceae than in other Archaea, most HGT events described earlier as Methanosarcina-specific seem to have occurred before the divergence of Methanosarcinaceae. Genes horizontally transferred from bacteria to archaea neither tend to be transferred with their regulators, nor in long operons.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0393-2) contains supplementary material, which is available to authorized users.

Frequent, independent transfers of a catabolic gene from bacteria to contrasted filamentous eukaryotes

Even genetically distant prokaryotes can exchange genes between them, and these horizontal gene transfer events play a central role in adaptation and evolution. While this was long thought to be restricted to prokaryotes, certain eukaryotes have acquired genes of bacterial origin. However, gene acquisitions in eukaryotes are thought to be much less important in magnitude than in prokaryotes. Here, we describe the complex evolutionary history of a bacterial catabolic gene that has been transferred repeatedly from different bacterial phyla to stramenopiles and fungi. Indeed, phylogenomic analysis pointed to multiple acquisitions of the gene in these filamentous eukaryotes-as many as 15 different events for 65 microeukaryotes. Furthermore, once transferred, this gene acquired introns and was found expressed in mRNA databases for most recipients. Our results show that effective inter-domain transfers and subsequent adaptation of a prokaryotic gene in eukaryotic cells can happen at an unprecedented magnitude.

Strategies and approaches in plasmidome studies-uncovering plasmid diversity disregarding of linear elements?

The term plasmid was originally coined for circular, extrachromosomal genetic elements. Today, plasmids are widely recognized not only as important factors facilitating genome restructuring but also as vehicles for the dissemination of beneficial characters within bacterial communities. Plasmid diversity has been uncovered by means of culture-dependent or -independent approaches, such as endogenous or exogenous plasmid isolation as well as PCR-based detection or transposon-aided capture, respectively. High-throughput-sequencing made possible to cover total plasmid populations in a given environment, i.e., the plasmidome, and allowed to address the quality and significance of self-replicating genetic elements. Since such efforts were and still are rather restricted to circular molecules, here we put equal emphasis on the linear plasmids which—despite their frequent occurrence in a large number of bacteria—are largely neglected in prevalent plasmidome conceptions.

Gene expression homeostasis and chromosome architecture

In rapidly growing populations of bacterial cells, including those of the model organism Escherichia coli, genes essential for growth--such as those involved in protein synthesis--are expressed at high levels; this is in contrast to many horizontally-acquired genes, which are maintained at low transcriptional levels. (1) This balance in gene expression states between 2 distinct classes of genes is established by a galaxy of transcriptional regulators, including the so-called nucleoid associated proteins (NAP) that contribute to shaping the chromosome. (2) Besides these active players in gene regulation, it is not too far-fetched to anticipate that genome organization in terms of how genes are arranged on the chromosome, (3) which is the result of long-drawn transactions among genome rearrangement processes and selection, and the manner in which it is structured inside the cell, plays a role in establishing this balance. A recent study from our group has contributed to the literature investigating the interplay between global transcriptional regulators and genome organization in establishing gene expression homeostasis. (4) In particular, we address a triangle of functional interactions among genome organization, gene expression homeostasis and horizontal gene transfer.

Comparative genomics reveals diversified CRISPR-Cas systems of globally distributed Microcystis aeruginosa, a freshwater bloom-forming cyanobacterium

Microcystis aeruginosa is one of the most common and dominant bloom-forming cyanobacteria in freshwater lakes around the world. Microcystis cells can produce toxic secondary metabolites, such as microcystins, which are harmful to human health. Two M. aeruginosa strains were isolated from two highly eutrophic lakes in China and their genomes were sequenced. Comparative genomic analysis was performed with the 12 other available M. aeruginosa genomes and closely related unicellular cyanobacterium. Each genome of M. aeruginosa containing at least one clustered regularly interspaced short palindromic repeat (CRISPR) locus and total 71 loci were identified, suggesting it is ubiquitous in M. aeruginosa genomes. In addition to the previously reported subtype I-D cas gene sets, three CAS subtypes I-A, III-A and III-B were identified and characterized in this study. Seven types of CRISPR direct repeat have close association with CAS subtype, confirming that different and specific secondary structures of CRISPR repeats are important for the recognition, binding and process of corresponding cas gene sets. Homology search of the CRISPR spacer sequences provides a history of not only resistance to bacteriophages and plasmids known to be associated with M. aeruginosa, but also the ability to target much more exogenous genetic material in the natural environment. These adaptive and heritable defense mechanisms play a vital role in keeping genomic stability and self-maintenance by restriction of horizontal gene transfer. Maintaining genomic stability and modulating genomic plasticity are both important evolutionary strategies for M. aeruginosa in adaptation and survival in various habitats.

Production level of tetrodotoxin in Aeromonas is associated with the copy number of a plasmid

Tetrodotoxin (TTX) has been identified from taxonomically diverse organisms. Artificial synthesis of TTX has been reported, but the biosynthetic pathway of TTX remains elusive. In this study, we found TTX producing ability was associated with the copy number of plasmid pNe-1 in Aeromonas strain Ne-1 during fermentation, suggesting that at least one gene encoding a TTX-synthesis enzyme is located on this plasmid. Compared with bacterial genomes, plasmids are small and easier to screen for genes associated with TTX biosynthesis. The approximately 100 kb genome of pNe-1 was sequenced. The plasmid contains 60 complete open reading frames (orfs) of which 32 (53.3%) encode hypothetical proteins. Seven genes are related to the type IV secretion system (T4SS) and 2 genes are related to transposons, indicating that the TTX-producing bacterium Aeromonas might have the ability to transfer the TTX biosynthesis gene via the conjugation and contagion of plasmid pNe-1. In addition, we unexpectedly found that Aeromonas Ne-1 contains unknown TTX-degrading materials, indicating there is a homeostatic mechanism to maintain a stable amount of TTX in the bacterium. These results will help us to better understand TTX biosynthesis, the bacterial origin of TTX, and TTX degradation.

The animal gut as a melting pot for horizontal gene transfer

In this minireview, we examine horizontal gene transfer (HGT) events in the mammalian gastrointestinal tract and their role in the evolutionary adaptation of microorganisms to the gut environment. We explore the notion of the mammalian gut as a melting pot of genetic exchange, resulting in the large extent of HGT occurrence.

The Turbulent Network Dynamics of Microbial Evolution and the Statistical Tree of Life

The wide spread and high rate of gene exchange and loss in the prokaryotic world translate into “network genomics”. The rates of gene gain and loss are comparable with the rate of point mutations but are substantially greater than the duplication rate. Thus, evolution of prokaryotes is primarily shaped by gene gain and loss. These processes are essential to prevent mutational meltdown of microbial populations by stopping Muller’s ratchet and appear to trigger emergence of major novel clades by opening up new ecological niches. At least some bacteria and archaea seem to have evolved dedicated devices for gene transfer. Despite the dominance of gene gain and loss, evolution of genes is intrinsically tree-like. The significant coherence between the topologies of numerous gene trees, particularly those for (nearly) universal genes, is compatible with the concept of a statistical tree of life, which forms the framework for reconstruction of the evolutionary processes in the prokaryotic world.

Classification of prokaryotic genetic replicators: between selfishness and altruism

Prokaryotes harbor a variety of genetic replicators, including plasmids, viruses, and chromosomes, each having different effects on the phenotype of the hosting cell. Here, we propose a classification for replicators of bacteria and archaea on the basis of their horizontal-transfer potential and the type of relationships (mutualistic, symbiotic, commensal, or parasitic) that they have with the host cell vehicle. Horizontal movement of replicators can be either active or passive, reflecting whether or not the replicator encodes the means to mediate its own transfer from one cell to another. Some replicators also have an infectious extracellular state, thus separating viruses from other mobile elements. From the perspective of the cell vehicle, the different types of replicators form a continuum from genuinely mutualistic to completely parasitic replicators. This classification provides a general framework for dissecting prokaryotic systems into evolutionarily meaningful components.

Babela massiliensis, a representative of a widespread bacterial phylum with unusual adaptations to parasitism in amoebae

Background

Only a small fraction of bacteria and archaea that are identifiable by metagenomics can be grown on standard media. Recent efforts on deep metagenomics sequencing, single-cell genomics and the use of specialized culture conditions (culturomics) increasingly yield novel microbes some of which represent previously uncharacterized phyla and possess unusual biological traits.

Results

We report isolation and genome analysis of Babela massiliensis, an obligate intracellular parasite of Acanthamoeba castellanii. B. massiliensis shows an unusual, fission mode of cell multiplication whereby large, polymorphic bodies accumulate in the cytoplasm of infected amoeba and then split into mature bacterial cells. This unique mechanism of cell division is associated with a deep degradation of the cell division machinery and delayed expression of the ftsZ gene. The genome of B. massiliensis consists of a circular chromosome approximately 1.12 megabase in size that encodes, 981 predicted proteins, 38 tRNAs and one typical rRNA operon. Phylogenetic analysis shows that B. massiliensis belongs to the putative bacterial phylum TM6 that so far was represented by the draft genome of the JCVI TM6SC1 bacterium obtained by single cell genomics and numerous environmental sequences.

Conclusions

Currently, B. massiliensis is the only cultivated member of the putative TM6 phylum. Phylogenomic analysis shows diverse taxonomic affinities for B. massiliensis genes, suggestive of multiple gene acquisitions via horizontal transfer from other bacteria and eukaryotes. Horizontal gene transfer is likely to be facilitated by the cohabitation of diverse parasites and symbionts inside amoeba. B. massiliensis encompasses many genes encoding proteins implicated in parasite-host interaction including the greatest number of ankyrin repeats among sequenced bacteria and diverse proteins related to the ubiquitin system. Characterization of B. massiliensis, a representative of a distinct bacterial phylum, thanks to its ability to grow in amoeba, reaffirms the critical role of diverse culture approaches in microbiology.

Reviewers

This article was reviewed by Dr. Igor Zhulin, Dr. Jeremy Selengut, and Pr Martijn Huynen.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0043-z) contains supplementary material, which is available to authorized users.

Distinguishing between productive and abortive promoters using a random forest classifier in Mycoplasma pneumoniae

Distinguishing between promoter-like sequences in bacteria that belong to true or abortive promoters, or to those that do not initiate transcription at all, is one of the important challenges in transcriptomics. To address this problem, we have studied the genome-reduced bacterium Mycoplasma pneumoniae, for which the RNAs associated with transcriptional start sites have been recently experimentally identified. We determined the contribution to transcription events of different genomic features: the –10, extended –10 and –35 boxes, the UP element, the bases surrounding the –10 box and the nearest-neighbor free energy of the promoter region. Using a random forest classifier and the aforementioned features transformed into scores, we could distinguish between true, abortive promoters and non-promoters with good –10 box sequences. The methods used in this characterization of promoters can be extended to other bacteria and have important applications for promoter design in bacterial genome engineering.

Archaeal Clusters of Orthologous Genes (arCOGs): An Update and Application for Analysis of Shared Features between Thermococcales, Methanococcales, and Methanobacteriales

With the continuously accelerating genome sequencing from diverse groups of archaea and bacteria, accurate identification of gene orthology and availability of readily expandable clusters of orthologous genes are essential for the functional annotation of new genomes. We report an update of the collection of archaeal Clusters of Orthologous Genes (arCOGs) to cover, on average, 91% of the protein-coding genes in 168 archaeal genomes. The new arCOGs were constructed using refined algorithms for orthology identification combined with extensive manual curation, including incorporation of the results of several completed and ongoing research projects in archaeal genomics. A new level of classification is introduced, superclusters that unit two or more arCOGs and more completely reflect gene family evolution than individual, disconnected arCOGs. Assessment of the current archaeal genome annotation in public databases indicates that consistent use of arCOGs can significantly improve the annotation quality. In addition to their utility for genome annotation, arCOGs also are a platform for phylogenomic analysis. We explore this aspect of arCOGs by performing a phylogenomic study of the Thermococci that are traditionally viewed as the basal branch of the Euryarchaeota. The results of phylogenomic analysis that involved both comparison of multiple phylogenetic trees and a search for putative derived shared characters by using phyletic patterns extracted from the arCOGs reveal a likely evolutionary relationship between the Thermococci, Methanococci, and Methanobacteria. The arCOGs are expected to be instrumental for a comprehensive phylogenomic study of the archaea.

No evidence of inhibition of horizontal gene transfer by CRISPR-Cas on evolutionary timescales

The CRISPR (clustered, regularly, interspaced, short, palindromic repeats)–Cas (CRISPR-associated genes) systems of archaea and bacteria provide adaptive immunity against viruses and other selfish elements and are believed to curtail horizontal gene transfer (HGT). Limiting acquisition of new genetic material could be one of the sources of the fitness cost of CRISPR–Cas maintenance and one of the causes of the patchy distribution of CRISPR–Cas among bacteria, and across environments. We sought to test the hypothesis that the activity of CRISPR–Cas in microbes is negatively correlated with the extent of recent HGT. Using three independent measures of HGT, we found no significant dependence between the length of CRISPR arrays, which reflects the activity of the immune system, and the estimated number of recent HGT events. In contrast, we observed a significant negative dependence between the estimated extent of HGT and growth temperature of microbes, which could be explained by the lower genetic diversity in hotter environments. We hypothesize that the relevant events in the evolution of resistance to mobile elements and proclivity for HGT, to which CRISPR–Cas systems seem to substantially contribute, occur on the population scale rather than on the timescale of species evolution.

Deconstruction of archaeal genome depict strategic consensus in core pathways coding sequence assembly

A comprehensive in silico analysis of 71 species representing the different taxonomic classes and physiological genre of the domain Archaea was performed. These organisms differed in their physiological attributes, particularly oxygen tolerance and energy metabolism. We explored the diversity and similarity in the codon usage pattern in the genes and genomes of these organisms, emphasizing on their core cellular pathways. Our thrust was to figure out whether there is any underlying similarity in the design of core pathways within these organisms. Analyses of codon utilization pattern, construction of hierarchical linear models of codon usage, expression pattern and codon pair preference pointed to the fact that, in the archaea there is a trend towards biased use of synonymous codons in the core cellular pathways and the Nc-plots appeared to display the physiological variations present within the different species. Our analyses revealed that aerobic species of archaea possessed a larger degree of freedom in regulating expression levels than could be accounted for by codon usage bias alone. This feature might be a consequence of their enhanced metabolic activities as a result of their adaptation to the relatively O²-rich environment. Species of archaea, which are related from the taxonomical viewpoint, were found to have striking similarities in their ORF structuring pattern. In the anaerobic species of archaea, codon bias was found to be a major determinant of gene expression. We have also detected a significant difference in the codon pair usage pattern between the whole genome and the genes related to vital cellular pathways, and it was not only species-specific but pathway specific too. This hints towards the structuring of ORFs with better decoding accuracy during translation. Finally, a codon-pathway interaction in shaping the codon design of pathways was observed where the transcription pathway exhibited a significantly different coding frequency signature.

Horizontal gene transfers in insects

Horizontal gene transfer is the transfer of genetic material across species boundaries. Although horizontal gene transfers are relatively rare in animals, the recent rapid accumulation of genomic data has identified increasing amounts of exogenous DNA inserts in insect genomes. Most of the horizontally acquired sequences appear to be non-functional; however, there is growing evidence that some genes are truly expressed and confer novel functions on the recipient insects. These include previously unavailable metabolic properties including digesting food, degrading toxins, providing resistance to pathogens, and facilitating an obligate mutualistic relationship with intracellular bacteria. A recent analysis revealed that an aphid gene of bacterial origin encodes a protein that is transported into the obligate symbiont, paralleling the evolution of endosymbiotic organelles.

Identification of horizontally transferred genes in the genus Colletotrichum reveals a steady tempo of bacterial to fungal gene transfer

BACKGROUND

Horizontal gene transfer (HGT) is the stable transmission of genetic material between organisms by means other than vertical inheritance. HGT has an important role in the evolution of prokaryotes but is relatively rare in eukaryotes. HGT has been shown to contribute to virulence in eukaryotic pathogens. We studied the importance of HGT in plant pathogenic fungi by identifying horizontally transferred genes in the genomes of three members of the genus Colletotrichum.

RESULTS

We identified eleven HGT events from bacteria into members of the genus Colletotrichum or their ancestors. The HGT events include genes involved in amino acid, lipid and sugar metabolism as well as lytic enzymes. Additionally, the putative minimal dates of transference were calculated using a time calibrated phylogenetic tree. This analysis reveals a constant flux of genes from bacteria to fungi throughout the evolution of subphylum Pezizomycotina.

CONCLUSIONS

Genes that are typically transferred by HGT are those that are constantly subject to gene duplication and gene loss. The functions of some of these genes suggest roles in niche adaptation and virulence. We found no evidence of a burst of HGT events coinciding with major geological events. In contrast, HGT appears to be a constant, albeit rare phenomenon in the Pezizomycotina, occurring at a steady rate during their evolution.

Robust identification of orthologues and paralogues for microbial pan-genomics using GET_HOMOLOGUES: a case study of pIncA/C plasmids

GET_HOMOLOGUES is an open-source software package written in Perl and R to define robust core- and pan-genomes by computing consensus clusters of orthologous gene families from whole-genome sequences using the bidirectional best-hit, COGtriangles, and OrthoMCL clustering algorithms. The granularity of the clusters can be fine-tuned by a user-configurable filtering strategy based on a combination of blastp pairwise alignment parameters, hmmscan-based scanning of Pfam domain composition of the proteins in each cluster, and a partial synteny criterion. We present detailed protocols to fit exponential and binomial mixture models to estimate core- and pan-genome sizes, compute pan-genome trees from the pan-genome matrix using a parsimony criterion, analyze and graphically represent the pan-genome structure, and identify lineage-specific gene families for the 12 complete pIncA/C plasmids currently available in NCBI's RefSeq. The software package, license, and detailed user manual can be downloaded for free for academic use from two mirrors: http://www.eead.csic.es/compbio/soft/gethoms.php and http://maya.ccg.unam.mx/soft/gethoms.php.

Why There Are No Essential Genes on Plasmids

Mobile genetic elements such as plasmids are important for the evolution of prokaryotes. It has been suggested that there are differences between functions coded for by mobile genes and those in the "core" genome and that these differences can be seen between plasmids and chromosomes. In particular, it has been suggested that essential genes, such as those involved in the formation of structural proteins or in basic metabolic functions, are rarely located on plasmids. We model competition between genotypically varying bacteria within a single population to investigate whether selection favors a chromosomal location for essential genes. We find that in general, chromosomal locations for essential genes are indeed favored. This is because the inheritance of chromosomes is more stable than that for plasmids. We define the "degradation" rate as the rate at which chance genetic processes, for example, mutation, deletion, or translocation, render essential genes nonfunctioning. The only way in which plasmids can be a location for functioning essential genes is if chromosomal genes degrade faster than plasmid genes. If the two degradation rates are equal, or if plasmid genes degrade faster than chromosomal genes, functioning essential genes will be found only on chromosomes.

The fundamental tradeoff in genomes and proteomes of prokaryotes established by the genetic code, codon entropy, and physics of nucleic acids and proteins

Background

Mutations in nucleotide sequences provide a foundation for genetic variability, and selection is the driving force of the evolution and molecular adaptation. Despite considerable progress in the understanding of selective forces and their compositional determinants, the very nature of underlying mutational biases remains unclear.

Results

We explore here a fundamental tradeoff, which analytically describes mutual adjustment of the nucleotide and amino acid compositions and its possible effect on the mutational biases. The tradeoff is determined by the interplay between the genetic code, optimization of the codon entropy, and demands on the structure and stability of nucleic acids and proteins.

Conclusion

The tradeoff is the unifying property of all prokaryotes regardless of the differences in their phylogenies, life styles, and extreme environments. It underlies mutational biases characteristic for genomes with different nucleotide and amino acid compositions, providing foundation for evolution and adaptation.

Reviewers

This article was reviewed by Eugene Koonin, Michael Gromiha, and Alexander Schleiffer.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-014-0029-2) contains supplementary material, which is available to authorized users.

Systematic genomic analysis reveals the complementary aerobic and anaerobic respiration capacities of the human gut microbiota

Because of the specific anatomical and physiological properties of the human intestine, a specific oxygen gradient builds up within this organ that influences the intestinal microbiota. The intestinal microbiome has been intensively studied in recent years, and certain respiratory substrates used by gut inhabiting microbes have been shown to play a crucial role in human health. Unfortunately, a systematic analysis has not been previously performed to determine the respiratory capabilities of human gut microbes (HGM). Here, we analyzed the distribution of aerobic and anaerobic respiratory reductases in 254 HGM genomes. In addition to the annotation of known enzymes, we also predicted a novel microaerobic reductase and novel thiosulfate reductase. Based on this comprehensive assessment of respiratory reductases in the HGM, we proposed a number of exchange pathways among different bacteria involved in the reduction of various nitrogen oxides. The results significantly expanded our knowledge of HGM metabolism and interactions in bacterial communities.

Integron associated mobile genes: Just a collection of plug in apps or essential components of cell network hardware?

Lateral gene transfer (LGT) impacts on the evolution of prokaryotes in both the short and long-term. The short-term impacts of mobilized genes are a concern to humans since LGT explains the global rise of multi drug resistant pathogens seen in the past 70 years. However, LGT has been a feature of prokaryotes from the earliest days of their existence and the concept of a bifurcating tree of life is not entirely applicable to prokaryotes since most genes in extant prokaryotic genomes have probably been acquired from other lineages. Successful transfer and maintenance of a gene in a new host is understandable if it acts independently of cell networks and confers an advantage. Antibiotic resistance provides an example of this whereby a gene can be advantageous in virtually any cell across broad species backgrounds. In a longer evolutionary context however laterally transferred genes can be assimilated into even essential cell networks. How this happens is not well understood and we discuss recent work that identifies a mobile gene, unique to a cell lineage, which is detrimental to the cell when lost. We also present some additional data and believe our emerging model will be helpful in understanding how mobile genes integrate into cell networks.

Large-scale dynamics of horizontal transfers

The widespread exchange of genes between bacteria must have consequences on the global architecture of their genomes, which are being found in the abundant genomic data available today. Most of the expansion of bacterial protein families can be attributed to transfer events, which are positively biased for smaller evolutionary distances between genomes, and more frequent for classes that are larger, when summed over all known bacteria. Moreover, “innovation” events where horizontal transfers carry exogenous evolutionary families appear to be less frequent for larger genomes. This dynamic expansion of evolutionary families is interconnected with the acquisition of new biological functions and thus with the size and distribution of the genes’ functional categories found on a genome. This commentary presents our recent contributions to this line of work and possible future directions.

Noncanonical cell-to-cell DNA transfer in Thermus spp. is insensitive to argonaute-mediated interference

Horizontal gene transfer drives the rapid evolution of bacterial populations. Classical processes that promote the lateral flow of genetic information are conserved throughout the prokaryotic world. However, some species have nonconserved transfer mechanisms that are not well known. This is the case for the ancient extreme thermophile Thermus thermophilus. In this work, we show that T. thermophilus strains are capable of exchanging large DNA fragments by a novel mechanism that requires cell-to-cell contacts and employs components of the natural transformation machinery. This process facilitates the bidirectional transfer of virtually any DNA locus but favors by 10-fold loci found in the megaplasmid over those in the chromosome. In contrast to naked DNA acquisition by transformation, the system does not activate the recently described DNA-DNA interference mechanism mediated by the prokaryotic Argonaute protein, thus allowing the organism to distinguish between DNA transferred from a mate and exogenous DNA acquired from unknown hosts. This Argonaute-mediated discrimination may be tentatively viewed as a strategy for safe sharing of potentially "useful" traits by the components of a given population of Thermus spp. without increasing the genome sizes of its individuals.

A meta-approach for improving the prediction and the functional annotation of ortholog groups

BACKGROUND

In comparative genomics, orthologs are used to transfer annotation from genes already characterized to newly sequenced genomes. Many methods have been developed for finding orthologs in sets of genomes. However, the application of different methods on the same proteome set can lead to distinct orthology predictions.

METHODS

We developed a method based on a meta-approach that is able to combine the results of several methods for orthologous group prediction. The purpose of this method is to produce better quality results by using the overlapping results obtained from several individual orthologous gene prediction procedures. Our method proceeds in two steps. The first aims to construct seeds for groups of orthologous genes; these seeds correspond to the exact overlaps between the results of all or several methods. In the second step, these seed groups are expanded by using HMM profiles.

RESULTS

We evaluated our method on two standard reference benchmarks, OrthoBench and Orthology Benchmark Service. Our method presents a higher level of accurately predicted groups than the individual input methods of orthologous group prediction. Moreover, our method increases the number of annotated orthologous pairs without decreasing the annotation quality compared to twelve state-of-the-art methods.

CONCLUSIONS

The meta-approach based method appears to be a reliable procedure for predicting orthologous groups. Since a large number of methods for predicting groups of orthologous genes exist, it is quite conceivable to apply this meta-approach to several combinations of different methods.

Estimation of prokaryotic supergenome size and composition from gene frequency distributions

BACKGROUND

Because prokaryotic genomes experience a rapid flux of genes, selection may act at a higher level than an individual genome. We explore a quantitative model of the distributed genome whereby groups of genomes evolve by acquiring genes from a fixed reservoir which we denote as supergenome. Previous attempts to understand the nature of the supergenome treated genomes as random, independent collections of genes and assumed that the supergenome consists of a small number of homogeneous sub-reservoirs. Here we explore the consequences of relaxing both assumptions.

RESULTS

We surveyed several methods for estimating the size and composition of the supergenome. The methods assumed that genomes were either random, independent samples of the supergenome or that they evolved from a common ancestor along a known tree via stochastic sampling from the reservoir. The reservoir was assumed to be either a collection of homogeneous sub-reservoirs or alternatively composed of genes with Gamma distributed gain probabilities. Empirical gene frequencies were used to either compute the likelihood of the data directly or first to reconstruct the history of gene gains and then compute the likelihood of the reconstructed numbers of gains.

CONCLUSIONS

Supergenome size estimates using the empirical gene frequencies directly are not robust with respect to the choice of the model. By contrast, using the gene frequencies and the phylogenetic tree to reconstruct multiple gene gains produces reliable estimates of the supergenome size and indicates that a homogeneous supergenome is more consistent with the data than a supergenome with Gamma distributed gain probabilities.

Processing of abasic site damaged lesions by APE1 enzyme on DNA adsorbed over normal and organomodified clay

The efficiency of the apurinic/apyrimidinic endonuclease (APE1) DNA repair enzyme in the processing of abasic site DNA damage lesions at precise location in DNA oligomer duplexes that are adsorbed on clay surfaces was evaluated. Three different forms of clay namely montmorillonite, quaternary ammonium salt modified montmorillonite and its boiled counterpart i.e. partially devoid of organic moiety were used for a comparative study of adsorption, desorption and DNA repair efficiency on their surfaces. The interaction between the DNA and the clay was analysed by X-ray diffraction, Atomic force microscopy, UV-Vis spectroscopy and Infrared spectroscopy. The abasic site cleavage efficiency of APE1 enzyme was quantitatively evaluated by polyacrylamide gel electrophoresis. Apart from the difference in the DNA adsorption or desorption capacity of the various forms of clay, substantial variation in the repair efficiency of abasic sites initiated by the APE1 enzyme on the clay surfaces was observed. The incision efficiency of APE1 enzyme at abasic sites was found to be greatly diminished, when the DNA was adsorbed over organomodified montmorillonite. The reduced repair activity indicates an important role of the pendant surfactant groups on the clay surfaces in directing APE1 mediated cleavage of abasic site DNA damage lesions.

Conserved patterns in bacterial genomes: a conundrum physically tailored by evolutionary tinkering

The proper functioning of bacteria is encoded in their genome at multiple levels or scales, each of which is constrained by specific physical forces. At the smallest spatial scales, interatomic forces dictate the folding and function of proteins and nucleic acids. On longer length scales, stochastic forces emerging from the thermal jiggling of proteins and RNAs impose strong constraints on the organization of genes along chromosomes, more particularly in the context of the building of nucleoprotein complexes and the operational mode of regulatory agents. At the cellular level, transcription, replication and cell division activities generate forces that act on both the internal structure and cellular location of chromosomes. The overall result is a complex multi-scale organization of genomes that reflects the evolutionary tinkering of bacteria. The goal of this review is to highlight avenues for deciphering this complexity by focusing on patterns that are conserved among evolutionarily distant bacteria. To this end, I discuss three different organizational scales: the protein structures, the chromosomal organization of genes and the global structure of chromosomes.

Genomes in turmoil: quantification of genome dynamics in prokaryote supergenomes

BACKGROUND

Genomes of bacteria and archaea (collectively, prokaryotes) appear to exist in incessant flux, expanding via horizontal gene transfer and gene duplication, and contracting via gene loss. However, the actual rates of genome dynamics and relative contributions of different types of event across the diversity of prokaryotes are largely unknown, as are the sizes of microbial supergenomes, i.e. pools of genes that are accessible to the given microbial species.

RESULTS

We performed a comprehensive analysis of the genome dynamics in 35 groups (34 bacterial and one archaeal) of closely related microbial genomes using a phylogenetic birth-and-death maximum likelihood model to quantify the rates of gene family gain and loss, as well as expansion and reduction. The results show that loss of gene families dominates the evolution of prokaryotes, occurring at approximately three times the rate of gain. The rates of gene family expansion and reduction are typically seven and twenty times less than the gain and loss rates, respectively. Thus, the prevailing mode of evolution in bacteria and archaea is genome contraction, which is partially compensated by the gain of new gene families via horizontal gene transfer. However, the rates of gene family gain, loss, expansion and reduction vary within wide ranges, with the most stable genomes showing rates about 25 times lower than the most dynamic genomes. For many groups, the supergenome estimated from the fraction of repetitive gene family gains includes about tenfold more gene families than the typical genome in the group although some groups appear to have vast, 'open' supergenomes.

CONCLUSIONS

Reconstruction of evolution for groups of closely related bacteria and archaea reveals an extremely rapid and highly variable flux of genes in evolving microbial genomes, demonstrates that extensive gene loss and horizontal gene transfer leading to innovation are the two dominant evolutionary processes, and yields robust estimates of the supergenome size.

Horizontal transfer of short and degraded DNA has evolutionary implications for microbes and eukaryotic sexual reproduction

Horizontal gene transfer in the form of long DNA fragments has changed our view of bacterial evolution. Recently, we discovered that such processes may also occur with the massive amounts of short and damaged DNA in the environment, and even with truly ancient DNA. Although it presently remains unclear how often it takes place in nature, horizontal gene transfer of short and damaged DNA opens up the possibility for genetic exchange across distinct species in both time and space. In this essay, we speculate on the potential evolutionary consequences of this phenomenon. We argue that it may challenge basic assumptions in evolutionary theory; that it may have distant origins in life's history; and that horizontal gene transfer should be viewed as an evolutionary strategy not only preceding but causally underpinning the evolution of sexual reproduction.

Linkage, mobility, and selfishness in the MazF family of bacterial toxins: a snapshot of bacterial evolution

Prokaryotic MazF family toxins cooccur with cognate antitoxins having divergent DNA-binding folds and can be of chromosomal or plasmid origin. Sequence similarity search was carried out to identify the Toxin–Antitoxin (TA) operons of MazF family followed by sequence analysis and phylogenetic studies. The genomic DNA upstream of the TA operons was searched for the presence of regulatory motifs. The MazF family toxins showed a conserved hydrophobic pocket in a multibinding site and are present in pathogenic bacteria. The toxins of the MazF family are associated with four main types of cognate antitoxin partners and cluster as a subfamily on the branches of the phylogenetic tree. This indicates that transmission of the entire operon is the dominant mode of inheritance. The plasmid borne TA modules were interspersed between the chromosomal TA modules of the same subfamily, compatible with a frequent interchange of TA genes between the chromosome and the plasmid akin to that observed for antibiotic resistance gens. The split network of the MazF family toxins showed the AbrB-linked toxins as a hub of horizontal gene transfer. Distinct motifs are present in the upstream region of each subfamily. The presence of MazF family TA modules in pathogenic bacteria and identification of a conserved binding pocket are significant for the development of novel antibacterials to disrupt the TA interaction. However, the role of TAs in stress resistance needs to be established. Phylogenetic studies provide insight into the evolution of MazF family TAs and effect on the bacterial genome.

Dark matter in archaeal genomes: a rich source of novel mobile elements, defense systems and secretory complexes

Microbial genomes encompass a sizable fraction of poorly characterized, narrowly spread fast-evolving genes. Using sensitive methods for sequences comparison and protein structure prediction, we performed a detailed comparative analysis of clusters of such genes, which we denote "dark matter islands", in archaeal genomes. The dark matter islands comprise up to 20% of archaeal genomes and show remarkable heterogeneity and diversity. Nevertheless, three classes of entities are common in these genomic loci: (a) integrated viral genomes and other mobile elements; (b) defense systems, and (c) secretory and other membrane-associated systems. The dark matter islands in the genome of thermophiles and mesophiles show similar general trends of gene content, but thermophiles are substantially enriched in predicted membrane proteins whereas mesophiles have a greater proportion of recognizable mobile elements. Based on this analysis, we predict the existence of several novel groups of viruses and mobile elements, previously unnoticed variants of CRISPR-Cas immune systems, and new secretory systems that might be involved in stress response, intermicrobial conflicts and biogenesis of novel, uncharacterized membrane structures.

Origin of giant viruses from smaller DNA viruses not from a fourth domain of cellular life

The numerous and diverse eukaryotic viruses with large double-stranded DNA genomes that at least partially reproduce in the cytoplasm of infected cells apparently evolved from a single virus ancestor. This major group of viruses is known as Nucleocytoplasmic Large DNA Viruses (NCLDV) or the proposed order Megavirales. Among the "Megavirales", there are three groups of giant viruses with genomes exceeding 500kb, namely Mimiviruses, Pithoviruses, and Pandoraviruses that hold the current record of viral genome size, about 2.5Mb. Phylogenetic analysis of conserved, ancestral NLCDV genes clearly shows that these three groups of giant viruses have three distinct origins within the "Megavirales". The Mimiviruses constitute a distinct family that is distantly related to Phycodnaviridae, Pandoraviruses originate from a common ancestor with Coccolithoviruses within the Phycodnaviridae family, and Pithoviruses are related to Iridoviridae and Marseilleviridae. Maximum likelihood reconstruction of gene gain and loss events during the evolution of the "Megavirales" indicates that each group of giant viruses evolved from viruses with substantially smaller and simpler gene repertoires. Initial phylogenetic analysis of universal genes, such as translation system components, encoded by some giant viruses, in particular Mimiviruses, has led to the hypothesis that giant viruses descend from a fourth, probably extinct domain of cellular life. The results of our comprehensive phylogenomic analysis of giant viruses refute the fourth domain hypothesis and instead indicate that the universal genes have been independently acquired by different giant viruses from their eukaryotic hosts.

The origins of cellular life

All life on earth can be naturally classified into cellular life forms and virus-like selfish elements, the latter being fully dependent on the former for their reproduction. Cells are reproducers that not only replicate their genome but also reproduce the cellular organization that depends on semipermeable, energy-transforming membranes and cannot be recovered from the genome alone, under the famous dictum of Rudolf Virchow, Omnis cellula e cellula. In contrast, simple selfish elements are replicators that can complete their life cycles within the host cell starting from genomic RNA or DNA alone. The origin of the cellular organization is the central and perhaps the hardest problem of evolutionary biology. I argue that the origin of cells can be understood only in conjunction with the origin and evolution of selfish genetic elements. A scenario of precellular evolution is presented that involves cohesion of the genomes of the emerging cellular life forms from primordial pools of small genetic elements that eventually segregated into hosts and parasites. I further present a model of the coevolution of primordial membranes and membrane proteins, discuss protocellular and non-cellular models of early evolution, and examine the habitats on the primordial earth that could have been conducive to precellular evolution and the origin of cells.

Comparative genomics of the bacterial genus Streptococcus illuminates evolutionary implications of species groups

Members of the genus Streptococcus within the phylum Firmicutes are among the most diverse and significant zoonotic pathogens. This genus has gone through considerable taxonomic revision due to increasing improvements of chemotaxonomic approaches, DNA hybridization and 16S rRNA gene sequencing. It is proposed to place the majority of streptococci into “species groups”. However, the evolutionary implications of species groups are not clear presently. We use comparative genomic approaches to yield a better understanding of the evolution of Streptococcus through genome dynamics, population structure, phylogenies and virulence factor distribution of species groups. Genome dynamics analyses indicate that the pan-genome size increases with the addition of newly sequenced strains, while the core genome size decreases with sequential addition at the genus level and species group level. Population structure analysis reveals two distinct lineages, one including Pyogenic, Bovis, Mutans and Salivarius groups, and the other including Mitis, Anginosus and Unknown groups. Phylogenetic dendrograms show that species within the same species group cluster together, and infer two main clades in accordance with population structure analysis. Distribution of streptococcal virulence factors has no obvious patterns among the species groups; however, the evolution of some common virulence factors is congruous with the evolution of species groups, according to phylogenetic inference. We suggest that the proposed streptococcal species groups are reasonable from the viewpoints of comparative genomics; evolution of the genus is congruent with the individual evolutionary trajectories of different species groups.

An operon of three transcriptional regulators controls horizontal gene transfer of the integrative and conjugative element ICEclc in Pseudomonas knackmussii B13

The integrative and conjugative element ICEclc is a mobile genetic element in Pseudomonas knackmussii B13, and an experimental model for a widely distributed group of elements in Proteobacteria. ICEclc is transferred from specialized transfer competent cells, which arise at a frequency of 3-5% in a population at stationary phase. Very little is known about the different factors that control the transfer frequency of this ICE family. Here we report the discovery of a three-gene operon encoded by ICEclc, which exerts global control on transfer initiation. The operon consists of three consecutive regulatory genes, encoding a TetR-type repressor MfsR, a MarR-type regulator and a LysR-type activator TciR. We show that MfsR autoregulates expression of the operon, whereas TciR is a global activator of ICEclc gene expression, but no clear role was yet found for MarR. Deletion of mfsR increases expression of tciR and marR, causing the proportion of transfer competent cells to reach almost 100% and transfer frequencies to approach 1 per donor. mfsR deletion also caused a two orders of magnitude loss in population viability, individual cell growth arrest and loss of ICEclc. This indicates that autoregulation is an important feature maintaining ICE transfer but avoiding fitness loss. Bioinformatic analysis showed that the mfsR-marR-tciR operon is unique for ICEclc and a few highly related ICE, whereas tciR orthologues occur more widely in a large variety of suspected ICE among Proteobacteria.

Integrative and conjugative elements (ICEs) are a relatively newly recognized class of mobile elements in bacteria, which integrate at one or more positions in a host chromosome, can be excised, circularized, and transfer by conjugation to a new recipient cell. Genome sequencing indicated that ICEs often carry genes with potential adaptive functions for the host. Various ICE-types have been described and ICEclc is a useful model for a wide class of elements found in Beta- and Gammaproteobacteria. Because ICEs normally remain “silent” in the host chromosome and often lack selectable markers, their lifestyle is difficult to study. One of the characteristics of ICEclc is that transfer is initiated in only 3-5% of donor cells in a population during stationary phase. Here, we describe an operon of three regulatory genes, two of which control the transfer initiation of ICEclc. Our findings suggest that the low transfer rate results from the repression of an activator and that this is essential to minimize the deleterious effect of hyper-activation of transfer initiation. While the individual regulatory genes are quite common on ICEs, they rarely occur in this configuration.

Cross-species gene-family fluctuations reveal the dynamics of horizontal transfers

Prokaryotes vary their protein repertoire mainly through horizontal transfer and gene loss. To elucidate the links between these processes and the cross-species gene-family statistics, we perform a large-scale data analysis of the cross-species variability of gene-family abundance (the number of members of the family found on a given genome). We find that abundance fluctuations are related to the rate of horizontal transfers. This is rationalized by a minimal theoretical model, which predicts this link. The families that are not captured by the model show abundance profiles that are markedly peaked around a mean value, possibly because of specific abundance selection. Based on these results, we define an abundance variability index that captures a family's evolutionary behavior (and thus some of its relevant functional properties) purely based on its cross-species abundance fluctuations. Analysis and model, combined, show a quantitative link between cross-species family abundance statistics and horizontal transfer dynamics, which can be used to analyze genome ‘flux’. Groups of families with different values of the abundance variability index correspond to genome sub-parts having different plasticity in terms of the level of horizontal exchange allowed by natural selection.

Predicting the minimal translation apparatus: lessons from the reductive evolution of mollicutes

Mollicutes is a class of parasitic bacteria that have evolved from a common Firmicutes ancestor mostly by massive genome reduction. With genomes under 1 Mbp in size, most Mollicutes species retain the capacity to replicate and grow autonomously. The major goal of this work was to identify the minimal set of proteins that can sustain ribosome biogenesis and translation of the genetic code in these bacteria. Using the experimentally validated genes from the model bacteria Escherichia coli and Bacillus subtilis as input, genes encoding proteins of the core translation machinery were predicted in 39 distinct Mollicutes species, 33 of which are culturable. The set of 260 input genes encodes proteins involved in ribosome biogenesis, tRNA maturation and aminoacylation, as well as proteins cofactors required for mRNA translation and RNA decay. A core set of 104 of these proteins is found in all species analyzed. Genes encoding proteins involved in post-translational modifications of ribosomal proteins and translation cofactors, post-transcriptional modifications of t+rRNA, in ribosome assembly and RNA degradation are the most frequently lost. As expected, genes coding for aminoacyl-tRNA synthetases, ribosomal proteins and initiation, elongation and termination factors are the most persistent (i.e. conserved in a majority of genomes). Enzymes introducing nucleotides modifications in the anticodon loop of tRNA, in helix 44 of 16S rRNA and in helices 69 and 80 of 23S rRNA, all essential for decoding and facilitating peptidyl transfer, are maintained in all species. Reconstruction of genome evolution in Mollicutes revealed that, beside many gene losses, occasional gains by horizontal gene transfer also occurred. This analysis not only showed that slightly different solutions for preserving a functional, albeit minimal, protein synthetizing machinery have emerged in these successive rounds of reductive evolution but also has broad implications in guiding the reconstruction of a minimal cell by synthetic biology approaches.

Implications of streamlining theory for microbial ecology

Whether a small cell, a small genome or a minimal set of chemical reactions with self-replicating properties, simplicity is beguiling. As Leonardo da Vinci reportedly said, 'simplicity is the ultimate sophistication'. Two diverging views of simplicity have emerged in accounts of symbiotic and commensal bacteria and cosmopolitan free-living bacteria with small genomes. The small genomes of obligate insect endosymbionts have been attributed to genetic drift caused by small effective population sizes (Ne). In contrast, streamlining theory attributes small cells and genomes to selection for efficient use of nutrients in populations where Ne is large and nutrients limit growth. Regardless of the cause of genome reduction, lost coding potential eventually dictates loss of function. Consequences of reductive evolution in streamlined organisms include atypical patterns of prototrophy and the absence of common regulatory systems, which have been linked to difficulty in culturing these cells. Recent evidence from metagenomics suggests that streamlining is commonplace, may broadly explain the phenomenon of the uncultured microbial majority, and might also explain the highly interdependent (connected) behavior of many microbial ecosystems. Streamlining theory is belied by the observation that many successful bacteria are large cells with complex genomes. To fully appreciate streamlining, we must look to the life histories and adaptive strategies of cells, which impose minimum requirements for complexity that vary with niche.

The precarious prokaryotic chromosome

Evolutionary selection for optimal genome preservation, replication, and expression should yield similar chromosome organizations in any type of cells. And yet, the chromosome organization is surprisingly different between eukaryotes and prokaryotes. The nuclear versus cytoplasmic accommodation of genetic material accounts for the distinct eukaryotic and prokaryotic modes of genome evolution, but it falls short of explaining the differences in the chromosome organization. I propose that the two distinct ways to organize chromosomes are driven by the differences between the global-consecutive chromosome cycle of eukaryotes and the local-concurrent chromosome cycle of prokaryotes. Specifically, progressive chromosome segregation in prokaryotes demands a single duplicon per chromosome, while other "precarious" features of the prokaryotic chromosomes can be viewed as compensations for this severe restriction.

Light-dependent expression of four cryptic archaeal circadian gene homologs

Circadian rhythms are important biological signals that have been found in almost all major groups of life from bacteria to man, yet it remains unclear if any members of the second major prokaryotic domain of life, the Archaea, also possess a biological clock. As an initial investigation of this question, we examined the regulation of four cyanobacterial-like circadian gene homologs present in the genome of the haloarchaeon Haloferax volcanii. These genes, designated cirA, cirB, cirC, and cirD, display similarity to the KaiC-family of cyanobacterial clock proteins, which act to regulate rhythmic gene expression and to control the timing of cell division. Quantitative RT-PCR analysis was used to examine the expression of each of the four cir genes in response to 12 h light/12 h dark cycles (LD 12:12) in H. volcanii during balanced growth. Our data reveal that there is an approximately two to sixteen-fold increase in cir gene expression when cells are shifted from light to constant darkness, and this pattern of gene expression oscillates with the light conditions in a rhythmic manner. Targeted single- and double-gene knockouts in the H. volcanii cir genes result in disruption of light-dependent, rhythmic gene expression, although it does not lead to any significant effect on growth under these conditions. Restoration of light-dependent, rhythmic gene expression was demonstrated by introducing, in trans, a wild-type copy of individual cir genes into knockout strains. These results are noteworthy as this is the first attempt to characterize the transcriptional expression and regulation of the ubiquitous kaiC homologs found among archaeal genomes.

Inferring gene function from evolutionary change in signatures of translation efficiency

Background

The genetic code is redundant, meaning that most amino acids can be encoded by more than one codon. Highly expressed genes tend to use optimal codons to increase the accuracy and speed of translation. Thus, codon usage biases provide a signature of the relative expression levels of genes, which can, uniquely, be quantified across the domains of life.

Results

Here we describe a general statistical framework to exploit this phenomenon and to systematically associate genes with environments and phenotypic traits through changes in codon adaptation. By inferring evolutionary signatures of translation efficiency in 911 bacterial and archaeal genomes while controlling for confounding effects of phylogeny and inter-correlated phenotypes, we linked 187 gene families to 24 diverse phenotypic traits. A series of experiments in Escherichia coli revealed that 13 of 15, 19 of 23, and 3 of 6 gene families with changes in codon adaptation in aerotolerant, thermophilic, or halophilic microbes. Respectively, confer specific resistance to, respectively, hydrogen peroxide, heat, and high salinity. Further, we demonstrate experimentally that changes in codon optimality alone are sufficient to enhance stress resistance. Finally, we present evidence that multiple genes with altered codon optimality in aerobes confer oxidative stress resistance by controlling the levels of iron and NAD(P)H.

Conclusions

Taken together, these results provide experimental evidence for a widespread connection between changes in translation efficiency and phenotypic adaptation. As the number of sequenced genomes increases, this novel genomic context method for linking genes to phenotypes based on sequence alone will become increasingly useful.

Bacterial genome instability

Bacterial genomes are remarkably stable from one generation to the next but are plastic on an evolutionary time scale, substantially shaped by horizontal gene transfer, genome rearrangement, and the activities of mobile DNA elements. This implies the existence of a delicate balance between the maintenance of genome stability and the tolerance of genome instability. In this review, we describe the specialized genetic elements and the endogenous processes that contribute to genome instability. We then discuss the consequences of genome instability at the physiological level, where cells have harnessed instability to mediate phase and antigenic variation, and at the evolutionary level, where horizontal gene transfer has played an important role. Indeed, this ability to share DNA sequences has played a major part in the evolution of life on Earth. The evolutionary plasticity of bacterial genomes, coupled with the vast numbers of bacteria on the planet, substantially limits our ability to control disease.