The nature and dynamics of bacterial genomes

Horizontally Acquired Biosynthesis Genes Boost Coxiella burnetii's Physiology

Coxiella burnetii, the etiologic agent of acute Q fever and chronic endocarditis, has a unique biphasic life cycle, which includes a metabolically active intracellular form that occupies a large lysosome-derived acidic vacuole. C. burnetii is the only bacterium known to thrive within such an hostile intracellular niche, and this ability is fundamental to its pathogenicity; however, very little is known about genes that facilitate Coxiella's intracellular growth. Recent studies indicate that C. burnetii evolved from a tick-associated ancestor and that the metabolic capabilities of C. burnetii are different from that of Coxiella-like bacteria found in ticks. Horizontally acquired genes that allow C. burnetii to infect and grow within mammalian cells likely facilitated the host shift; however, because of its obligate intracellular replication, C. burnetii would have lost most genes that have been rendered redundant due to the availability of metabolites within the host cell. Based on these observations, we reasoned that horizontally derived biosynthetic genes that have been retained in the reduced genome of C. burnetii are ideal candidates to begin to uncover its intracellular metabolic requirements. Our analyses identified a large number of putative foreign-origin genes in C. burnetii, including tRNA^Glu2 that is potentially required for heme biosynthesis, and genes involved in the production of lipopolysaccharide—a virulence factor, and of critical metabolites such as fatty acids and biotin. In comparison to wild-type C. burnetii, a strain that lacks tRNA^Glu2 exhibited reduced growth, indicating its importance to Coxiella's physiology. Additionally, by using chemical agents that block heme and biotin biosyntheses, we show that these pathways are promising targets for the development of new anti-Coxiella therapies.

Vector Integration Sites Identification for Gene-Trap Screening in Mammalian Haploid Cells

Forward genetic screens using retroviral (or transposon) gene-trap vectors in a haploid genome revolutionized the investigation of molecular networks in mammals. However, the sequencing data generated by Phenotypic interrogation followed by Tag sequencing (PhiT-seq) were not well characterized. The analysis of human and mouse haploid screens allowed us to describe PhiT-seq data and to define quality control steps. Moreover, we identified several blind spots in both haploid genomes where gene-trap vectors can hardly integrate. Integration of transcriptomic data improved the performance of candidate gene identification. Furthermore, we experimented with various statistical tests to account for biological replicates in PhiT-seq and investigated the effect of normalization methods and other parameters on the performance. Finally, we developed: VISITs, a dedicated pipeline for analyzing PhiT-seq data (https://sourceforge.net/projects/visits/).

Defying Muller's Ratchet: Ancient Heritable Endobacteria Escape Extinction through Retention of Recombination and Genome Plasticity

Heritable endobacteria, which are transmitted from one host generation to the next, are subjected to evolutionary forces that are different from those experienced by free-living bacteria. In particular, they suffer consequences of Muller’s ratchet, a mechanism that leads to extinction of small asexual populations due to fixation of slightly deleterious mutations combined with the random loss of the most-fit genotypes, which cannot be recreated without recombination. Mycoplasma-related endobacteria (MRE) are heritable symbionts of fungi from two ancient lineages, Glomeromycota (arbuscular mycorrhizal fungi) and Mucoromycotina. Previous studies revealed that MRE maintain unusually diverse populations inside their hosts and may have been associated with fungi already in the early Paleozoic. Here we show that MRE are vulnerable to genomic degeneration and propose that they defy Muller’s ratchet thanks to retention of recombination and genome plasticity. We suggest that other endobacteria may be capable of raising similar defenses against Muller’s ratchet.

Pyomelanin-producing Pseudomonas aeruginosa selected during chronic infections have a large chromosomal deletion which confers resistance to pyocins

When bacterial lineages make the transition from free-living to permanent association with hosts, they can undergo massive gene losses, for which the selective forces within host tissues are unknown. We identified here melanogenic clinical isolates of Pseudomonas aeruginosa with large chromosomal deletions (66 to 270 kbp) and characterized them to investigate how they were selected. When compared with their wild-type parents, melanogenic mutants (i) exhibited a lower fitness in growth conditions found in human tissues, such as hyperosmolarity and presence of aminoglycoside antibiotics, (ii) narrowed their metabolic spectrum with a growth disadvantage with particular carbon sources, including aromatic amino acids and acyclic terpenes, suggesting a reduction of metabolic flexibility. Despite an impaired fitness in rich media, melanogenic mutants can inhibit their wild-type parents and compete with them in coculture. Surprisingly, melanogenic mutants became highly resistant to two intraspecific toxins, the S-pyocins AP41 and S1. Our results suggest that pyocins produced within a population of infecting P. aeruginosa may have selected for bacterial mutants that underwent massive gene losses and that were adapted to the life in diverse bacterial communities in the human host. Intraspecific interactions may therefore be an important factor driving the continuing evolution of pathogens during host infections.

Comparative genomic analysis of six new-found integrative conjugative elements (ICEs) in Vibrio alginolyticus

Background

Vibrio alginolyticus is ubiquitous in marine and estuarine environments. In 2012–2013, SXT/R391-like integrative conjugative elements (ICEs) in environmental V. alginolyticus strains were discovered and found to occur in 8.9 % of 192 V. alginolyticus strains, which suggests that V. alginolyticus may be a natural pool possessing resourceful ICEs. However, complete ICE sequences originating from this bacterium have not been reported, which represents a significant barrier to characterizing the ICEs of this bacterium and exploring their relationships with other ICEs. In the present study, we acquired six ICE sequences from five V. alginolyticus strains and performed a comparative analysis of these ICE genomes.

Results

A sequence analysis showed that there were only 14 variable bases dispersed between ICEValE0601 and ICEValHN492. ICEValE0601 and ICEValHN492 were treated as the same ICE. ICEValA056-1, ICEValE0601 and ICEValHN492 integrate into the 5′ end of the host’s prfC gene, and their Int and Xis share at least 97 % identity with their counterparts from SXT. ICEValE0601 or ICEValHN492 contain 50 of 52 conserved core genes in the SXT/R391 ICEs (not s025 or s026). ICEValA056-2, ICEValHN396 and ICEValHN437 have a different tRNA-ser integration site and a distinct int/xis module; however, the remaining backbone genes are highly similar to their counterparts in SXT/R391 ICEs. DNA sequences inserted into hotspot and variable regions of the ICEs are of various sizes. The variable genes of six ICEs encode a large array of functions to bestow various adaptive abilities upon their hosts, and only ICEValA056-1 contains drug-resistant genes. Many variable genes have orthologous and functionally related genes to those found in SXT/R391 ICEs, such as genes coding for a toxin-antitoxin system, a restriction-modification system, helicases and endonucleases. Six ICEs also contain a large number of unique genes or gene clusters that were not found in other ICEs. Six ICEs harbor more abundant transposase genes compared with other parts of their host genomes. A phylogenetic analysis indicated that transposase genes in these ICEs are highly diverse.

Conclusions

ICEValA056-1, ICEValE0601 and ICEValHN492 are typical members of the SXT/R391 family. ICEValA056-2, ICEValHN396 and ICEValHN437 form a new atypical group belonging to the SXT/R391 family. In addition to the many genes found to be present in other ICEs, six ICEs contain a large number of unique genes or gene clusters that were not found in other ICEs. ICEs may serve as a carrier for transposable genetic elements (TEs) and largely facilitate the dissemination of TEs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-016-0692-9) contains supplementary material, which is available to authorized users.

The genomic signatures of Shigella evolution, adaptation and geographical spread

Shigella spp. are some of the key pathogens responsible for the global burden of diarrhoeal disease. These facultative intracellular bacteria belong to the family Enterobacteriaceae, together with other intestinal pathogens, such as Escherichia coli and Salmonella spp. The genus Shigella comprises four different species, each consisting of several serogroups, all of which show phenotypic similarity, including invasive pathogenicity. DNA sequencing suggests that this similarity results from the convergent evolution of different Shigella spp. founders. Here, we review the evolutionary relationships between Shigella spp. and E . coli, and we highlight how the genomic plasticity of these bacteria and their acquisition of a distinctive virulence plasmid have enabled the development of such highly specialized pathogens. Furthermore, we discuss the insights that genotyping and whole-genome sequencing have provided into the phylogenetics and intercontinental spread of Shigella spp.

Microevolution Analysis of Bacillus coahuilensis Unveils Differences in Phosphorus Acquisition Strategies and Their Regulation

Bacterial genomes undergo numerous events of gene losses and gains that generate genome variability among strains of the same species (microevolution). Our aim was to compare the genomes and relevant phenotypes of three Bacillus coahuilensis strains from two oligotrophic hydrological systems in the Cuatro Ciénegas Basin (México), to unveil the environmental challenges that this species cope with, and the microevolutionary differences in these genotypes. Since the strains were isolated from a low P environment, we placed emphasis on the search of different phosphorus acquisition strategies. The three B. coahuilensis strains exhibited similar numbers of coding DNA sequences, of which 82% (2,893) constituted the core genome, and 18% corresponded to accessory genes. Most of the genes in this last group were associated with mobile genetic elements (MGEs) or were annotated as hypothetical proteins. Ten percent of the pangenome consisted of strain-specific genes. Alignment of the three B. coahuilensis genomes indicated a high level of synteny and revealed the presence of several genomic islands. Unexpectedly, one of these islands contained genes that encode the 2-keto-3-deoxymannooctulosonic acid (Kdo) biosynthesis enzymes, a feature associated to cell walls of Gram-negative bacteria. Some microevolutionary changes were clearly associated with MGEs. Our analysis revealed inconsistencies between phenotype and genotype, which we suggest result from the impossibility to map regulatory features to genome analysis. Experimental results revealed variability in the types and numbers of auxotrophies between the strains that could not consistently be explained by in silico metabolic models. Several intraspecific differences in preferences for carbohydrate and phosphorus utilization were observed. Regarding phosphorus recycling, scavenging, and storage, variations were found between the three genomes. The three strains exhibited differences regarding alkaline phosphatase that revealed that in addition to gene gain and loss, regulation adjustment of gene expression also has contributed to the intraspecific diversity of B. coahuilensis.

The Phylogeny of Bacillus cereus sensu lato

The three main species of the Bacillus cereus sensu lato , B. cereus , B. thuringiensis , and B. anthracis , were recognized and established by the early 1900s because they each exhibited distinct phenotypic traits. B. thuringiensis isolates and their parasporal crystal proteins have long been established as a natural pesticide and insect pathogen. B. anthracis , the etiological agent for anthrax, was used by Robert Koch in the 19th century as a model to develop the germ theory of disease, and B. cereus , a common soil organism, is also an occasional opportunistic pathogen of humans. In addition to these three historical species designations, are three less-recognized and -understood species: B. mycoides , B. weihenstephanensis , and B. pseudomycoides . All of these “species” combined comprise the Bacillus cereus sensu lato group. Despite these apparently clear phenotypic definitions, early molecular approaches to separate the first three by various DNA hybridization and 16S/23S ribosomal sequence analyses led to some “confusion” because there were limited differences to differentiate between these species. These and other results have led to frequent suggestions that a taxonomic change was warranted to reclassify this group to a single species. But the pathogenic properties of B. anthracis and the biopesticide applications of B. thuringiensis appear to “have outweighed pure taxonomic considerations” and the separate species categories are still being maintained. B. cereus sensu lato represents a classic example of a now common bacterial species taxonomic quandary.

Comparative Genomics between Two Xenorhabdus bovienii Strains Highlights Differential Evolutionary Scenarios within an Entomopathogenic Bacterial Species

Bacteria of the genus Xenorhabdus are symbionts of soil entomopathogenic nematodes of the genus Steinernema. This symbiotic association constitutes an insecticidal complex active against a wide range of insect pests. Within Xenorhabdus bovienii species, the X. bovienii CS03 strain (Xb CS03) is nonvirulent when directly injected into lepidopteran insects, and displays a low virulence when associated with its Steinernema symbiont. The genome of Xb CS03 was sequenced and compared with the genome of a virulent strain, X. bovienii SS-2004 (Xb SS-2004). The genome size and content widely differed between the two strains. Indeed, Xb CS03 had a large genome containing several specific loci involved in the inhibition of competitors, including a few NRPS-PKS loci (nonribosomal peptide synthetases and polyketide synthases) producing antimicrobial molecules. Consistently, Xb CS03 had a greater antimicrobial activity than Xb SS-2004. The Xb CS03 strain contained more pseudogenes than Xb SS-2004. Decay of genes involved in the host invasion and exploitation (toxins, invasins, or extracellular enzymes) was particularly important in Xb CS03. This may provide an explanation for the nonvirulence of the strain when injected into an insect host. We suggest that Xb CS03 and Xb SS-2004 followed divergent evolutionary scenarios to cope with their peculiar life cycle. The fitness strategy of Xb CS03 would involve competitor inhibition, whereas Xb SS-2004 would quickly and efficiently kill the insect host. Hence, Xenorhabdus strains would have widely divergent host exploitation strategies, which impact their genome structure.

Evidence for the recent origin of a bacterial protein-coding, overlapping orphan gene by evolutionary overprinting

BACKGROUND

Gene duplication is believed to be the classical way to form novel genes, but overprinting may be an important alternative. Overprinting allows entirely novel proteins to evolve de novo, i.e., formerly non-coding open reading frames within functional genes become expressed. Only three cases have been described for Escherichia coli. Here, a fourth example is presented.

RESULTS

RNA sequencing revealed an open reading frame weakly transcribed in cow dung, coding for 101 residues and embedded completely in the -2 reading frame of citC in enterohemorrhagic E. coli. This gene is designated novel overlapping gene, nog1. The promoter region fused to gfp exhibits specific activities and 5' rapid amplification of cDNA ends indicated the transcriptional start 40-bp upstream of the start codon. nog1 was strand-specifically arrested in translation by a nonsense mutation silent in citC. This Nog1-mutant showed a phenotype in competitive growth against wild type in the presence of MgCl2. Small differences in metabolite concentrations were also found. Bioinformatic analyses propose Nog1 to be inner membrane-bound and to possess at least one membrane-spanning domain. A phylogenetic analysis suggests that the orphan gene nog1 arose by overprinting after Escherichia/Shigella separated from the other γ-proteobacteria.

CONCLUSIONS

Since nog1 is of recent origin, non-essential, short, weakly expressed and only marginally involved in E. coli's central metabolism, we propose that this gene is in an initial stage of evolution. While we present specific experimental evidence for the existence of a fourth overlapping gene in enterohemorrhagic E. coli, we believe that this may be an initial finding only and overlapping genes in bacteria may be more common than is currently assumed by microbiologists.

The Phylogeny of Rickettsia Using Different Evolutionary Signatures: How Tree-Like is Bacterial Evolution?

Rickettsia is a genus of intracellular bacteria whose hosts and transmission strategies are both impressively diverse, and this is reflected in a highly dynamic genome. Some previous studies have described the evolutionary history of Rickettsia as non-tree-like, due to incongruity between phylogenetic reconstructions using different portions of the genome. Here, we reconstruct the Rickettsia phylogeny using whole-genome data, including two new genomes from previously unsampled host groups. We find that a single topology, which is supported by multiple sources of phylogenetic signal, well describes the evolutionary history of the core genome. We do observe extensive incongruence between individual gene trees, but analyses of simulations over a single topology and interspersed partitions of sites show that this is more plausibly attributed to systematic error than to horizontal gene transfer. Some conflicting placements also result from phylogenetic analyses of accessory genome content (i.e., gene presence/absence), but we argue that these are also due to systematic error, stemming from convergent genome reduction, which cannot be accommodated by existing phylogenetic methods. Our results show that, even within a single genus, tests for gene exchange based on phylogenetic incongruence may be susceptible to false positives.

Evolutionary Divergence of Aggregatibacter actinomycetemcomitans

Gram-negative facultative Aggregatibacter actinomycetemcomitans is an oral pathogen associated with periodontitis. The genetic heterogeneity among A. actinomycetemcomitans strains has been long recognized. This study provides a comprehensive genomic analysis of A. actinomycetemcomitans and the closely related nonpathogenic Aggregatibacter aphrophilus. Whole genome sequencing by Illumina MiSeq platform was performed for 31 A. actinomycetemcomitans and 2 A. aphrophilus strains. Sequence similarity analysis shows a total of 3,220 unique genes across the 2 species, where 1,550 are core genes present in all genomes and 1,670 are variable genes (accessory genes) missing in at least 1 genome. Phylogenetic analysis based on 397 concatenated core genes distinguished A. aphrophilus and A. actinomycetemcomitans. The latter was in turn divided into 5 clades: clade b (serotype b), clade c (serotype c), clade e/f (serotypes e and f), clade a/d (serotypes a and d), and clade e' (serotype e strains). Accessory genes accounted for 14.1% to 23.2% of the A. actinomycetemcomitans genomes, with a majority belonging to the category of poorly characterized by Cluster of Orthologous Groups classification. These accessory genes were often organized into genomic islands (n = 387) with base composition biases, suggesting their acquisitions via horizontal gene transfer. There was a greater degree of similarity in gene content and genomic islands among strains within clades than between clades. Strains of clade e' isolated from human were found to be missing the genomic island that carries genes encoding cytolethal distending toxins. Taken together, the results suggest a pattern of sequential divergence, starting from the separation of A. aphrophilus and A. actinomycetemcomitans through gain and loss of genes and ending with the divergence of the latter species into distinct clades and serotypes. With differing constellations of genes, the A. actinomycetemcomitans clades may have evolved distinct adaptation strategies to the human oral cavity.

Genomic and transcriptomic insights into the efficient entomopathogenicity of Bacillus thuringiensis

Bacillus thuringiensis has been globally used as a microbial pesticide for over 70 years. However, information regarding its various adaptions and virulence factors and their roles in the entomopathogenic process remains limited. In this work, we present the complete genomes of two industrially patented Bacillus thuringiensis strains (HD-1 and YBT-1520). A comparative genomic analysis showed a larger and more complicated genome constitution that included novel insecticidal toxicity-related genes (ITRGs). All of the putative ITRGs were summarized according to the steps of infection. A comparative genomic analysis showed that highly toxic strains contained significantly more ITRGs, thereby providing additional strategies for infection, immune evasion, and cadaver utilization. Furthermore, a comparative transcriptomic analysis suggested that a high expression of these ITRGs was a key factor in efficient entomopathogenicity. We identified an active extra urease synthesis system in the highly toxic strains that may aid B. thuringiensis survival in insects (similar to previous results with well-known pathogens). Taken together, these results explain the efficient entomopathogenicity of B. thuringiensis. It provides novel insights into the strategies used by B. thuringiensis to resist and overcome host immune defenses and helps identify novel toxicity factors.

Found and Lost: The Fates of Horizontally Acquired Genes in Arthropod-Symbiotic Spiroplasma

Horizontal gene transfer (HGT) is an important mechanism that contributed to biological diversity, particularly in bacteria. Through acquisition of novel genes, the recipient cell may change its ecological preference and the process could promote speciation. In this study, we determined the complete genome sequence of two Spiroplasma species for comparative analyses and inferred the putative gene gains and losses. Although most Spiroplasma species are symbionts of terrestrial insects, Spiroplasma eriocheiris has evolved to be a lethal pathogen of freshwater crustaceans. We found that approximately 7% of the genes in this genome may have originated from HGT and these genes expanded the metabolic capacity of this organism. Through comparison with the closely related Spiroplasma atrichopogonis, as well as other more divergent lineages, our results indicated that these HGT events could be traced back to the most recent common ancestor of these two species. However, most of these horizontally acquired genes have been pseudogenized in S. atrichopogonis, suggesting that they did not contribute to the fitness of this lineage that maintained the association with terrestrial insects. Thus, accumulation of small deletions that disrupted these foreign genes was not countered by natural selection. On the other hand, the long-term survival of these horizontally acquired genes in the S. eriocheiris genome hinted that they might play a role in the ecological shift of this species. Finally, the implications of these findings and the conflicts among gene content, 16S rRNA gene sequencing, and serological typing, are discussed in light of defining bacterial species.

Trichodesmium genome maintains abundant, widespread noncoding DNA in situ, despite oligotrophic lifestyle

Understanding the evolution of the free-living, cyanobacterial, diazotroph Trichodesmium is of great importance because of its critical role in oceanic biogeochemistry and primary production. Unlike the other >150 available genomes of free-living cyanobacteria, only 63.8% of the Trichodesmium erythraeum (strain IMS101) genome is predicted to encode protein, which is 20-25% less than the average for other cyanobacteria and nonpathogenic, free-living bacteria. We use distinctive isolates and metagenomic data to show that low coding density observed in IMS101 is a common feature of the Trichodesmium genus, both in culture and in situ. Transcriptome analysis indicates that 86% of the noncoding space is expressed, although the function of these transcripts is unclear. The density of noncoding, possible regulatory elements predicted in Trichodesmium, when normalized per intergenic kilobase, was comparable and twofold higher than that found in the gene-dense genomes of the sympatric cyanobacterial genera Synechococcus and Prochlorococcus, respectively. Conserved Trichodesmium noncoding RNA secondary structures were predicted between most culture and metagenomic sequences, lending support to the structural conservation. Conservation of these intergenic regions in spatiotemporally separated Trichodesmium populations suggests possible genus-wide selection for their maintenance. These large intergenic spacers may have developed during intervals of strong genetic drift caused by periodic blooms of a subset of genotypes, which may have reduced effective population size. Our data suggest that transposition of selfish DNA, low effective population size, and high-fidelity replication allowed the unusual "inflation" of noncoding sequence observed in Trichodesmium despite its oligotrophic lifestyle.

Making novel proteins from pseudogenes

MOTIVATION

Recently, we made synthetic proteins from non-coding DNA of Escherichia coli. Encouraged by this, we asked: can we artificially express pseudogenes into novel and functional proteins? What kind of structures would be generated? Would these proteins be stable? How would the organism respond to the artificial reactivation of pseudogenes?

RESULTS

To answer these questions, we studied 16 full-length protein equivalents of pseudogenes. The sequence-based predictions indicated interesting molecular and cellular functional roles for pseudogene-derived proteins. Most of the proteins were predicted to be involved in the amino acid biosynthesis, energy metabolism, purines and pyrimidine biosynthesis, central intermediary metabolism, transport and binding. Interestingly, many of the pseudogene-derived proteins were predicted to be enzymes. Furthermore, proteins showed strong evidence of stable tertiary structures. The prediction scores for structure, function and stability were found to be favorable in most of the cases.

IMPACT

To our best knowledge, this is the first such report that predicts the possibility of making functional and stable proteins from pseudogenes. In future, it would be interesting to experimentally synthesize and validate these predictions.

Taking the pseudo out of pseudogenes

Pseudogenes are defined as fragments of once-functional genes that have been silenced by one or more nonsense, frameshift or missense mutations. Despite continuing increases in the speed of sequencing and annotating bacterial genomes, the identification and categorisation of pseudogenes remains problematic. Even when identified, pseudogenes are considered to be rare and tend to be ignored. On the contrary, pseudogenes are surprisingly prevalent and can persist for long evolutionary time periods, representing a record of once-functional genetic characteristics. Most importantly, pseudogenes provide an insight into prokaryotic evolutionary history as a record of phenotypic traits that have been lost. Focusing on the intracellular and symbiotic bacteria in which pseudogenes predominate, this review discusses the importance of identifying pseudogenes to fully understand the abilities of bacteria, and to understand prokaryotes within their evolutionary context.

The interactions of bacteria with fungi in soil: emerging concepts

In this chapter, we review the existing literature on bacterial-fungal interactions in soil, exploring the role fungi may play for soil bacteria as providers of hospitable niches. A focus is placed on the mycosphere, i.e., the narrow zone of influence of fungal hyphae on the external soil milieu, in which hypha-associated bacterial cells dwell. Evidence is brought forward for the contention that the hyphae of both mycorrhizal and saprotrophic fungi serve as providers of ecological opportunities in a grossly carbon-limited soil, as a result of their release of carbonaceous compounds next to the provision of a colonizable surface. Soil bacteria of particular nature are postulated to have adapted to such selection pressures, evolving to the extent that they acquired capabilities that allow them to thrive in the novel habitat created by the emerging fungal hyphae. The mechanisms involved in the interactions and the modes of genetic adaptation of the mycosphere dwellers are discussed, with an emphasis on one key mycosphere-adapted bacterium, Burkholderia terrae BS001. In this discussion, we interrogate the positive interactions between soil fungi and bacteria, and refrain from considering negative interactions.

Comparative analyses of nonpathogenic, opportunistic, and totally pathogenic mycobacteria reveal genomic and biochemical variabilities and highlight the survival attributes of Mycobacterium tuberculosis

Mycobacterial evolution involves various processes, such as genome reduction, gene cooption, and critical gene acquisition. Our comparative genome size analysis of 44 mycobacterial genomes revealed that the nonpathogenic (NP) genomes were bigger than those of opportunistic (OP) or totally pathogenic (TP) mycobacteria, with the TP genomes being smaller yet variable in size—their genomic plasticity reflected their ability to evolve and survive under various environmental conditions. From the 44 mycobacterial species, 13 species, representing TP, OP, and NP, were selected for genomic-relatedness analyses. Analysis of homologous protein-coding genes shared between Mycobacterium indicus pranii (NP), Mycobacterium intracellulare ATCC 13950 (OP), and Mycobacterium tuberculosis H37Rv (TP) revealed that 4,995 (i.e., ~95%) M. indicaus pranii proteins have homology with M. intracellulare, whereas the homologies among M. indicus pranii, M. intracellulare ATCC 13950, and M. tuberculosis H37Rv were significantly lower. A total of 4,153 (~79%) M. indicus pranii proteins and 4,093 (~79%) M. intracellulare ATCC 13950 proteins exhibited homology with the M. tuberculosis H37Rv proteome, while 3,301 (~82%) and 3,295 (~82%) M. tuberculosis H37Rv proteins showed homology with M. indicus pranii and M. intracellulare ATCC 13950 proteomes, respectively. Comparative metabolic pathway analyses of TP/OP/NP mycobacteria showed enzymatic plasticity between M. indicus pranii (NP) and M. intracellulare ATCC 13950 (OP), Mycobacterium avium 104 (OP), and M. tuberculosis H37Rv (TP). Mycobacterium tuberculosis seems to have acquired novel alternate pathways with possible roles in metabolism, host-pathogen interactions, virulence, and intracellular survival, and by implication some of these could be potential drug targets.

The complete sequence analysis of Mycobacterium indicus pranii, a novel species of Mycobacterium shown earlier to have strong immunomodulatory properties and currently in use for the treatment of leprosy, places it evolutionarily at the point of transition to pathogenicity. With the purpose of establishing the importance of M. indicus pranii in providing insight into the virulence mechanism of tuberculous and nontuberculous mycobacteria, we carried out comparative genomic and proteomic analyses of 44 mycobacterial species representing nonpathogenic (NP), opportunistic (OP), and totally pathogenic (TP) mycobacteria. Our results clearly placed M. indicus pranii as an ancestor of the M. avium complex. Analyses of comparative metabolic pathways between M. indicus pranii (NP), M. tuberculosis (TP), and M. intracellulare (OP) pointed to the presence of novel alternative pathways in M. tuberculosis with implications for pathogenesis and survival in the human host and identification of new drug targets.

Molecular evolution of the substrate utilization strategies and putative virulence factors in mosquito-associated Spiroplasma species

Comparative genomics provides a powerful tool to characterize the genetic differences among species that may be linked to their phenotypic variations. In the case of mosquito-associated Spiroplasma species, such approach is useful for the investigation of their differentiations in substrate utilization strategies and putative virulence factors. Among the four species that have been assessed for pathogenicity by artificial infection experiments, Spiroplasma culicicola and S. taiwanense were found to be pathogenic, whereas S. diminutum and S. sabaudiense were not. Intriguingly, based on the species phylogeny, the association with mosquito hosts and the gain or loss of pathogenicity in these species appears to have evolved independently. Through comparison of their complete genome sequences, we identified the genes and pathways that are shared by all or specific to one of these four species. Notably, we found that a glycerol-3-phosphate oxidase gene (glpO) is present in S. culicicola and S. taiwanense but not in S. diminutum or S. sabaudiense. Because this gene is involved in the production of reactive oxygen species and has been demonstrated as a major virulence factor in Mycoplasma, this distribution pattern suggests that it may be linked to the observed differences in pathogenicity among these species as well. Moreover, through comparative analysis with other Spiroplasma, Mycoplasma, and Mesoplasma species, we found that the absence of glpO in S. diminutum and S. sabaudiense is best explained by independent losses. Finally, our phylogenetic analyses revealed possible recombination of glpO between distantly related lineages and local rearrangements of adjacent genes.

Proteome-wide subcellular topologies of E. coli polypeptides database (STEPdb)

Cell compartmentalization serves both the isolation and the specialization of cell functions. After synthesis in the cytoplasm, over a third of all proteins are targeted to other subcellular compartments. Knowing how proteins are distributed within the cell and how they interact is a prerequisite for understanding it as a whole. Surface and secreted proteins are important pathogenicity determinants. Here we present the STEP database (STEPdb) that contains a comprehensive characterization of subcellular localization and topology of the complete proteome of Escherichia coli. Two widely used E. coli proteomes (K-12 and BL21) are presented organized into thirteen subcellular classes. STEPdb exploits the wealth of genetic, proteomic, biochemical, and functional information on protein localization, secretion, and targeting in E. coli, one of the best understood model organisms. Subcellular annotations were derived from a combination of bioinformatics prediction, proteomic, biochemical, functional, topological data and extensive literature re-examination that were refined through manual curation. Strong experimental support for the location of 1553 out of 4303 proteins was based on 426 articles and some experimental indications for another 526. Annotations were provided for another 320 proteins based on firm bioinformatic predictions. STEPdb is the first database that contains an extensive set of peripheral IM proteins (PIM proteins) and includes their graphical visualization into complexes, cellular functions, and interactions. It also summarizes all currently known protein export machineries of E. coli K-12 and pairs them, where available, with the secretory proteins that use them. It catalogs the Sec- and TAT-utilizing secretomes and summarizes their topological features such as signal peptides and transmembrane regions, transmembrane topologies and orientations. It also catalogs physicochemical and structural features that influence topology such as abundance, solubility, disorder, heat resistance, and structural domain families. Finally, STEPdb incorporates prediction tools for topology (TMHMM, SignalP, and Phobius) and disorder (IUPred) and implements the BLAST2STEP that performs protein homology searches against the STEPdb.

Whole genome sequencing and analysis of plant growth promoting bacteria isolated from the rhizosphere of plantation crops coconut, cocoa and arecanut

Coconut, cocoa and arecanut are commercial plantation crops that play a vital role in the Indian economy while sustaining the livelihood of more than 10 million Indians. According to 2012 Food and Agricultural organization's report, India is the third largest producer of coconut and it dominates the production of arecanut worldwide. In this study, three Plant Growth Promoting Rhizobacteria (PGPR) from coconut (CPCRI-1), cocoa (CPCRI-2) and arecanut (CPCRI-3) characterized for the PGP activities have been sequenced. The draft genome sizes were 4.7 Mb (56% GC), 5.9 Mb (63.6% GC) and 5.1 Mb (54.8% GB) for CPCRI-1, CPCRI-2, CPCRI-3, respectively. These genomes encoded 4056 (CPCRI-1), 4637 (CPCRI-2) and 4286 (CPCRI-3) protein-coding genes. Phylogenetic analysis revealed that both CPCRI-1 and CPCRI-3 belonged to Enterobacteriaceae family, while, CPCRI-2 was a Pseudomonadaceae family member. Functional annotation of the genes predicted that all three bacteria encoded genes needed for mineral phosphate solubilization, siderophores, acetoin, butanediol, 1-aminocyclopropane-1-carboxylate (ACC) deaminase, chitinase, phenazine, 4-hydroxybenzoate, trehalose and quorum sensing molecules supportive of the plant growth promoting traits observed in the course of their isolation and characterization. Additionally, in all the three CPCRI PGPRs, we identified genes involved in synthesis of hydrogen sulfide (H2S), which recently has been proposed to aid plant growth. The PGPRs also carried genes for central carbohydrate metabolism indicating that the bacteria can efficiently utilize the root exudates and other organic materials as energy source. Genes for production of peroxidases, catalases and superoxide dismutases that confer resistance to oxidative stresses in plants were identified. Besides these, genes for heat shock tolerance, cold shock tolerance and glycine-betaine production that enable bacteria to survive abiotic stress were also identified.

On the limits of computational functional genomics for bacterial lifestyle prediction

We review the level of genomic specificity regarding actinobacterial pathogenicity. As they occupy various niches in diverse habitats, one may assume the existence of lifestyle-specific genomic features. We include 240 actinobacteria classified into four pathogenicity classes: human pathogens (HPs), broad-spectrum pathogens (BPs), opportunistic pathogens (OPs) and non-pathogenic (NP). We hypothesize: (H1) Pathogens (HPs and BPs) possess specific pathogenicity signature genes. (H2) The same holds for OPs. (H3) Broad-spectrum and exclusively HPs cannot be distinguished from each other because of an observation bias, i.e. many HPs might yet be unclassified BPs. (H4) There is no intrinsic genomic characteristic of OPs compared with pathogens, as small mutations are likely to play a more dominant role to survive the immune system. To study these hypotheses, we implemented a bioinformatics pipeline that combines evolutionary sequence analysis with statistical learning methods (Random Forest with feature selection, model tuning and robustness analysis). Essentially, we present orthologous gene sets that computationally distinguish pathogens from NPs (H1). We further show a clear limit in differentiating OPs from both NPs (H2) and pathogens (H4). HPs may also not be distinguished from bacteria annotated as BPs based only on a small set of orthologous genes (H3), as many HPs might as well target a broad range of mammals but have not been annotated accordingly. In conclusion, we illustrate that even in the post-genome era and despite next-generation sequencing technology, our ability to efficiently deduce real-world conclusions, such as pathogenicity classification, remains quite limited.

Recent molecular insights into rickettsial pathogenesis and immunity

Human infections with arthropod-borne Rickettsia species remain a major global health issue, causing significant morbidity and mortality. Epidemic typhus due to Rickettsia prowazekii has an established reputation as the 'scourge of armies', and as a major determinant of significant 'historical turning points'. No suitable vaccines for human use are currently available to prevent rickettsial diseases. The unique lifestyle features of rickettsiae include obligate intracellular parasitism, intracytoplasmic niche within the host cell, predilection for infection of microvascular endothelium in mammalian hosts, association with arthropods and the tendency for genomic reduction. The fundamental research in the field of Rickettsiology has witnessed significant recent progress in the areas of pathogen adhesion/invasion and host immune responses, as well as the genomics, proteomics, metabolomics, phylogenetics, motility and molecular manipulation of important rickettsial pathogens. The focus of this review article is to capture a snapshot of the latest developments pertaining to the mechanisms of rickettsial pathogenesis and immunity.

SearchDOGS bacteria, software that provides automated identification of potentially missed genes in annotated bacterial genomes

We report the development of SearchDOGS Bacteria, software to automatically detect missing genes in annotated bacterial genomes by combining BLAST searches with comparative genomics. Having successfully applied the approach to yeast genomes, we redeveloped SearchDOGS to function as a standalone, downloadable package, requiring only a set of GenBank annotation files as input. The software automatically generates a homology structure using reciprocal BLAST and a synteny-based method; this is followed by a scan of the entire genome of each species for unannotated genes. Results are provided in a HTML interface, providing coordinates, BLAST results, syntenic location, omega values (Ka/Ks, where Ks is the number of synonymous substitutions per synonymous site and Ka is the number of nonsynonymous substitutions per nonsynonymous site) for protein conservation estimates, and other information for each candidate gene. Using SearchDOGS Bacteria, we identified 155 gene candidates in the Shigella boydii sb227 genome, including 56 candidates of length < 60 codons. SearchDOGS Bacteria has two major advantages over currently available annotation software. First, it outperforms current methods in terms of sensitivity and is highly effective at identifying small or highly diverged genes. Second, as a freely downloadable package, it can be used with unpublished or confidential data.

Engineering reduced evolutionary potential for synthetic biology

The field of synthetic biology seeks to engineer reliable and predictable behaviors in organisms from collections of standardized genetic parts. However, unlike other types of machines, genetically encoded biological systems are prone to changes in their designed sequences due to mutations in their DNA sequences after these devices are constructed and deployed. Thus, biological engineering efforts can be confounded by undesired evolution that rapidly breaks the functions of parts and systems, particularly when they are costly to the host cell to maintain. Here, we explain the fundamental properties that determine the evolvability of biological systems. Then, we use this framework to review current efforts to engineer the DNA sequences that encode synthetic biology devices and the genomes of their microbial hosts to reduce their ability to evolve and therefore increase their genetic reliability so that they maintain their intended functions over longer timescales.

Experimental evolution of pseudogenization and gene loss in a plant RNA virus

Viruses have evolved highly streamlined genomes and a variety of mechanisms to compress them, suggesting that genome size is under strong selection. Horizontal gene transfer has, on the other hand, played an important role in virus evolution. However, evolution cannot integrate initially nonfunctional sequences into the viral genome if they are rapidly purged by selection. Here we report on the experimental evolution of pseudogenization in virus genomes using a plant RNA virus expressing a heterologous gene. When long 9-week passages were performed, the added gene was lost in all lineages, whereas viruses with large genomic deletions were fixed in only two out of ten 3-week lineages and none in 1-week lineages. Illumina next-generation sequencing revealed considerable convergent evolution in the 9- and 3-week lineages with genomic deletions. Genome size was correlated to within-host competitive fitness, although there was no correlation with virus accumulation or virulence. Within-host competitive fitness of the 3-week virus lineages without genomic deletions was higher than for the 1-week lineages. Our results show that the strength of selection for a reduced genome size and the rate of pseudogenization depend on demographic conditions. Moreover, for the 3-week passage condition, we observed increases in within-host fitness, whereas selection was not strong enough to quickly remove the nonfunctional heterologous gene. These results suggest a demographically determined "sweet spot" might exist, where heterologous insertions are not immediately lost while evolution can act to integrate them into the viral genome.

Complete genomes of two dipteran-associated spiroplasmas provided insights into the origin, dynamics, and impacts of viral invasion in spiroplasma

Spiroplasma is a genus of wall-less, low-GC, Gram-positive bacteria with helical morphology. As commensals or pathogens of plants, insects, ticks, or crustaceans, they are closely related with mycoplasmas and form a monophyletic group (Spiroplasma–Entomoplasmataceae–Mycoides) with Mycoplasma mycoides and its relatives. In this study, we report the complete genome sequences of Spiroplasma chrysopicola and S. syrphidicola from the Chrysopicola clade. These species form the sister group to the Citri clade, which includes several well-known pathogenic spiroplasmas. Surprisingly, these two newly available genomes from the Chrysopicola clade contain no plectroviral genes, which were found to be highly repetitive in the previously sequenced genomes from the Citri clade. Based on the genome alignment and patterns of GC-skew, these two Chrysopicola genomes appear to be relatively stable, rather than being highly rearranged as those from the Citri clade. Phylogenetic analyses suggest that the susceptibility to plectroviral invasion probably originated in the common ancestor of the Citri clade or one of its subclades. This susceptibility may be attributed to the absence of antiviral systems found in the Chrysopicola clade. Using the virus-free genomes of the Chrysopicola clade as references, we inferred the putative viral integration sites in the Citri genomes. Comparisons of syntenic regions suggest that the extensive viral invasion in the Citri clade promoted genome rearrangements and expansions. More importantly, the viral invasion may have facilitated horizontal gene transfers that contributed to adaptation in the Citri clade.

Comparison of metabolic capacities and inference of gene content evolution in mosquito-associated Spiroplasma diminutum and S. taiwanense

Mosquitoes are hosts of several Spiroplasma species that belong to different serogroups. To investigate the genetic mechanisms that may be involved in the utilization of similar hosts in these phylogenetically distinct bacteria, we determined the complete genome sequences of Spiroplasma diminutum and S. taiwanense for comparative analysis. The genome alignment indicates that their chromosomal organization is highly conserved, which is in sharp contrast to the elevated genome instabilities observed in other Spiroplasma lineages. Examination of the substrate utilization strategies revealed that S. diminutum can use a wide range of carbohydrates, suggesting that it is well suited to living in the gut (and possibly the circulatory system) of its mosquito hosts. In comparison, S. taiwanense has lost several carbohydrate utilization genes and acquired additional sets of oligopeptide transporter genes through tandem duplications, suggesting that proteins from digested blood meal or lysed host cells may be an important nutrient source. Moreover, one glycerol-3-phosphate oxidase gene (glpO) was found in S. taiwanense but not S. diminutum. This gene is linked to the production of reactive oxygen species and has been shown to be a major virulence factor in Mycoplasma mycoides. This finding may explain the pathogenicity of S. taiwanense observed in previous artificial infection experiments, while no apparent effect was found for S. diminutum. To infer the gene content evolution at deeper divergence levels, we incorporated other Mollicutes genomes for comparative analyses. The results suggest that the losses of biosynthetic pathways are a recurrent theme in these host-associated bacteria.

Predicting genetic and ecological characteristics of bacterial species by comparing the trajectories of dN/dS and dI/dS in bacterial genomes

Indel (insertion/deletion) causes gene disruption and is considered to be deleterious like non-synonymous mutation during the evolution of bacterial genomes. The trajectory of dN/dS (the ratio of non-synonymous to synonymous mutation) has been found to decrease exponentially over time, but the trajectory of dI/dS (the ratio of indel to synonymous mutation) has not been thoroughly explored yet. Here we compared the patterns of dN/dS and dI/dS for several bacterial species. The majority of them showed a much steeper dI/dS trajectory than the dN/dS trajectory, suggesting that indel was more deleterious than non-synonymous mutation and therefore was more difficult to fix in genomes. However, the naturally competent bacteria, or those with a much lower genetic barrier for DNA exchange, presented a reverse relationship between dI/dS and dN/dS, indicative of an exceptional ability to tolerate horizontal genetic transfer. The result suggests that plotting of dN/dS and dI/dS trajectories can help to predict the genetic and ecological characteristics of bacterial species.

Comparative genomic analysis indicates that niche adaptation of terrestrial Flavobacteria is strongly linked to plant glycan metabolism

Flavobacteria are important members of aquatic and terrestrial bacterial communities, displaying extreme variations in lifestyle, geographical distribution and genome size. They are ubiquitous in soil, but are often strongly enriched in the rhizosphere and phyllosphere of plants. In this study, we compared the genome of a root-associated Flavobacterium that we recently isolated, physiologically characterized and sequenced, to 14 additional Flavobacterium genomes, in order to pinpoint characteristics associated with its high abundance in the rhizosphere. Interestingly, flavobacterial genomes vary in size by approximately two-fold, with terrestrial isolates having predominantly larger genomes than those from aquatic environments. Comparative functional gene analysis revealed that terrestrial and aquatic Flavobacteria generally segregated into two distinct clades. Members of the aquatic clade had a higher ratio of peptide and protein utilization genes, whereas members of the terrestrial clade were characterized by a significantly higher abundance and diversity of genes involved in metabolism of carbohydrates such as xylose, arabinose and pectin. Interestingly, genes encoding glycoside hydrolase (GH) families GH78 and GH106, responsible for rhamnogalacturonan utilization (exclusively associated with terrestrial plant hemicelluloses), were only present in terrestrial clade genomes, suggesting adaptation of the terrestrial strains to plant-related carbohydrate metabolism. The Peptidase/GH ratio of aquatic clade Flavobacteria was significantly higher than that of terrestrial strains (1.7±0.7 and 9.7±4.7, respectively), supporting the concept that this relation can be used to infer Flavobacterium lifestyles. Collectively, our research suggests that terrestrial Flavobacteria are highly adapted to plant carbohydrate metabolism, which appears to be a key to their profusion in plant environments.

Horizontal gene transfer from diverse bacteria to an insect genome enables a tripartite nested mealybug symbiosis

The smallest reported bacterial genome belongs to Tremblaya princeps, a symbiont of Planococcus citri mealybugs (PCIT). Tremblaya PCIT not only has a 139 kb genome, but possesses its own bacterial endosymbiont, Moranella endobia. Genome and transcriptome sequencing, including genome sequencing from a Tremblaya lineage lacking intracellular bacteria, reveals that the extreme genomic degeneracy of Tremblaya PCIT likely resulted from acquiring Moranella as an endosymbiont. In addition, at least 22 expressed horizontally transferred genes from multiple diverse bacteria to the mealybug genome likely complement missing symbiont genes. However, none of these horizontally transferred genes are from Tremblaya, showing that genome reduction in this symbiont has not been enabled by gene transfer to the host nucleus. Our results thus indicate that the functioning of this three-way symbiosis is dependent on genes from at least six lineages of organisms and reveal a path to intimate endosymbiosis distinct from that followed by organelles.

Whole genome comparison of six Crocosphaera watsonii strains with differing phenotypes

Crocosphaera watsonii, a unicellular nitrogen-fixing cyanobacterium found in oligotrophic oceans, is important in marine carbon and nitrogen cycles. Isolates of C. watsonii can be separated into at least two phenotypes with environmentally important differences, indicating possibly distinct ecological roles and niches. To better understand the evolutionary history and variation in metabolic capabilities among strains and phenotypes, this study compared the genomes of six C. watsonii strains, three from each phenotypic group, which had been isolated over several decades from multiple ocean basins. While a substantial portion of each genome was nearly identical to sequences in the other strains, a few regions were identified as specific to each strain and phenotype, some of which help explain observed phenotypic features. Overall, the small-cell type strains had smaller genomes and a relative loss of genetic capabilities, while the large-cell type strains were characterized by larger genomes, some genetic redundancy, and potentially increased adaptations to iron and phosphorus limitation. As such, strains with shared phenotypes were evolutionarily more closely related than those with the opposite phenotype, regardless of isolation location or date. Unexpectedly, the genome of the type-strain for the species, C. watsonii WH8501, was quite unusual even among strains with a shared phenotype, indicating it may not be an ideal representative of the species. The genome sequences and analyses reported in this study will be important for future investigations of the proposed differences in adaptation of the two phenotypes to nutrient limitation, and to identify phenotype-specific distributions in natural Crocosphaera populations.

Speedy speciation in a bacterial microcosm: new species can arise as frequently as adaptations within a species

Microbiologists are challenged to explain the origins of enormous numbers of bacterial species worldwide. Contributing to this extreme diversity may be a simpler process of speciation in bacteria than in animals and plants, requiring neither sexual nor geographical isolation between nascent species. Here, we propose and test a novel hypothesis for the extreme diversity of bacterial species-that splitting of one population into multiple ecologically distinct populations (cladogenesis) may be as frequent as adaptive improvements within a single population's lineage (anagenesis). We employed a set of experimental microcosms to address the relative rates of adaptive cladogenesis and anagenesis among the descendants of a Bacillus subtilis clone, in the absence of competing species. Analysis of the evolutionary trajectories of genetic markers indicated that in at least 7 of 10 replicate microcosm communities, the original population founded one or more new, ecologically distinct populations (ecotypes) before a single anagenetic event occurred within the original population. We were able to support this inference by identifying putative ecotypes formed in these communities through differences in genetic marker association, colony morphology and microhabitat association; we then confirmed the ecological distinctness of these putative ecotypes in competition experiments. Adaptive mutations leading to new ecotypes appeared to be about as common as those improving fitness within an existing ecotype. These results suggest near parity of anagenesis and cladogenesis rates in natural populations that are depauperate of bacterial diversity.

Applying Shannon's information theory to bacterial and phage genomes and metagenomes

All sequence data contain inherent information that can be measured by Shannon's uncertainty theory. Such measurement is valuable in evaluating large data sets, such as metagenomic libraries, to prioritize their analysis and annotation, thus saving computational resources. Here, Shannon's index of complete phage and bacterial genomes was examined. The information content of a genome was found to be highly dependent on the genome length, GC content, and sequence word size. In metagenomic sequences, the amount of information correlated with the number of matches found by comparison to sequence databases. A sequence with more information (higher uncertainty) has a higher probability of being significantly similar to other sequences in the database. Measuring uncertainty may be used for rapid screening for sequences with matches in available database, prioritizing computational resources, and indicating which sequences with no known similarities are likely to be important for more detailed analysis.

Natural genome diversity of AI-2 quorum sensing in Escherichia coli: conserved signal production but labile signal reception

Quorum sensing (QS) regulates the onset of bacterial social responses in function to cell density having an important impact in virulence. Autoinducer-2 (AI-2) is a signal that has the peculiarity of mediating both intra- and interspecies bacterial QS. We analyzed the diversity of all components of AI-2 QS across 44 complete genomes of Escherichia coli and Shigella strains. We used phylogenetic tools to study its evolution and determined the phenotypes of single-deletion mutants to predict phenotypes of natural strains. Our analysis revealed many likely adaptive polymorphisms both in gene content and in nucleotide sequence. We show that all natural strains possess the signal emitter (the luxS gene), but many lack a functional signal receptor (complete lsr operon) and the ability to regulate extracellular signal concentrations. This result is in striking contrast with the canonical species-specific QS systems where one often finds orphan receptors, without a cognate synthase, but not orphan emitters. Our analysis indicates that selection actively maintains a balanced polymorphism for the presence/absence of a functional lsr operon suggesting diversifying selection on the regulation of signal accumulation and recognition. These results can be explained either by niche-specific adaptation or by selection for a coercive behavior where signal-blind emitters benefit from forcing other individuals in the population to haste in cooperative behaviors.

Factors behind junk DNA in bacteria

Although bacterial genomes have been traditionally viewed as being very compact, with relatively low amounts of repetitive and non-coding DNA, this view has dramatically changed in recent years. The increase of available complete bacterial genomes has revealed that many species present abundant repetitive DNA (i.e., insertion sequences, prophages or paralogous genes) and that many of these sequences are not functional but can have evolutionary consequences as concerns the adaptation to specialized host-related ecological niches. Comparative genomics analyses with close relatives that live in non-specialized environments reveal the nature and fate of this bacterial junk DNA. In addition, the number of insertion sequences and pseudogenes, as well as the size of the intergenic regions, can be used as markers of the evolutionary stage of a genome.

OASIS: an automated program for global investigation of bacterial and archaeal insertion sequences

Insertion sequences (ISs) are simple transposable elements present in most bacterial and archaeal genomes and play an important role in genomic evolution. The recent expansion of sequenced genomes offers the opportunity to study ISs comprehensively, but this requires efficient and accurate tools for IS annotation. We have developed an open-source program called OASIS, or Optimized Annotation System for Insertion Sequences, which automatically annotates ISs within sequenced genomes. OASIS annotations of 1737 bacterial and archaeal genomes offered an unprecedented opportunity to examine IS evolution. At a broad scale, we found that most IS families are quite widespread; however, they are not present randomly across taxa. This may indicate differential loss, barriers to exchange and/or insufficient time to equilibrate across clades. The number of ISs increases with genome length, but there is both tremendous variation and no increase in IS density for genomes >2 Mb. At the finer scale of recently diverged genomes, the proportion of shared IS content falls sharply, suggesting loss and/or emergence of barriers to successful cross-infection occurs rapidly. Surprisingly, even after controlling for 16S rRNA sequence divergence, the same ISs were more likely to be shared between genomes labeled as the same species rather than as different species.

Metagenomics and its connection to microbial community organization

Microbes dominate most global biogeochemical cycles, and microbial metagenomics (studying the collective microbial genomes) provides invaluable new insights into microbial systems, independent of cultivation. Metagenomic approaches targeting specific genes, e.g. small subunit (ssu) ribosomal RNA (rRNA), can be used to investigate microbial community organization by efficiently showing which taxa of organisms are present, while shotgun approaches show all genes and can indicate what functions the organisms are capable of. But collecting and organizing comprehensive shotgun data is extremely challenging and costly, and, in theory, predicting functionalities from microbial identities alone would save immense effort. However, we don’t yet know to what extent such predictions are applicable.

Next-generation and whole-genome sequencing in the diagnostic clinical microbiology laboratory

The identification and/or prediction of the antimicrobial resistance of microorganisms in clinical materials solely by molecular means in the diagnostic microbiology laboratory is not novel. However, the ability to sequence multitudes of bacterial genomes and deliver and interpret the resultant sequence information in near "real-time" is the basis of next-generation sequencing (NGS) technologies. There have been numerous applications and successes of NGS applications in the clinical and public health domain. However, none have, as yet, delivered perhaps the most sought after application, i.e., the generation of microbial sequence data for "real-time" patient management. In this review, we discuss the use of NGS and whole-genome sequencing (WGS) of microorganisms as a logical next step for the routine diagnosis of infection and the prediction of antimicrobial susceptibility in the clinical microbiology laboratory.

The major facilitator superfamily (MFS) revisited

The major facilitator superfamily (MFS) is the largest known superfamily of secondary carriers found in the biosphere. It is ubiquitously distributed throughout virtually all currently recognized organismal phyla. This superfamily currently (2012) consists of 74 families, each of which is usually concerned with the transport of a certain type of substrate. Many of these families, defined phylogenetically, do not include even a single member that is functionally characterized. In this article, we probe the evolutionary origins of these transporters, providing evidence that they arose from a single 2-transmembrane segment (TMS) hairpin structure that triplicated to give a 6-TMS unit that duplicated to a 12-TMS protein, the most frequent topological type of these permeases. We globally examine MFS protein topologies, focusing on exceptional proteins that deviate from the norm. Nine distantly related families appear to have members with 14 TMSs in which the extra two are usually centrally localized between the two 6-TMS repeat units. They probably have arisen by intragenic duplication of an adjacent hairpin. This alternative topology probably arose multiple times during MFS evolution. Convincing evidence for MFS permeases with fewer than 12 TMSs was not forthcoming, leading to the suggestion that all 12 TMSs are required for optimal function. Some homologs appear to have 13, 14, 15 or 16 TMSs, and the probable locations of the extra TMSs were identified. A few MFS permeases are fused to other functional domains or are fully duplicated to give 24-TMS proteins with dual functions. Finally, the MFS families with no known function were subjected to genomic context analyses leading to functional predictions.

How microbiology helps define the rhizome of life

In contrast to the tree of life (TOF) theory, species are mosaics of gene sequences with different origins. Observations of the extensive lateral sequence transfers in all organisms have demonstrated that the genomes of all life forms are collections of genes with different evolutionary histories that cannot be represented by a single TOF. Moreover, genes themselves commonly have several origins due to recombination. The human genome is not free from recombination events, so it is a mosaic like other organisms' genomes. Recent studies have demonstrated evidence for the integration of parasitic DNA into the human genome. Lateral transfer events have been accepted as major contributors of genome evolution in free-living bacteria. Furthermore, the accumulation of genomic sequence data provides evidence for extended genetic exchanges in intracellular bacteria and suggests that such events constitute an agent that promotes and maintains all bacterial species. Archaea and viruses also form chimeras containing primarily bacterial but also eukaryotic sequences. In addition to lateral transfers, orphan genes are indicative of the fact that gene creation is a permanent and unsettled phenomenon. Currently, a rhizome may more adequately represent the multiplicity and de novo creation of a genome. We wanted to confirm that the term “rhizome” in evolutionary biology applies to the entire cellular life history. This view of evolution should resemble a clump of roots representing the multiple origins of the repertoires of the genes of each species.