Genome analysis of Minibacterium massiliensis highlights the convergent evolution of water-living bacteria

Unraveling the genetic basis of microbial metal resistance: Shift from mendelian to systems biology

Microbial metal resistance, a trait that enables microorganisms to withstand high levels of toxic metals, has been studied for over a century. The significance of uncovering these mechanisms goes beyond basic science as they have implications for human health through their connection to microbial pathogenesis, metal bioremediation, and biomining. Recent advances in analytical chemistry and molecular biology have accelerated the discovery and understanding of genetic mechanisms underlying microbial metal resistance, identifying specific metal resistance genes and their operons. The emergence of omics tools has further propelled research towards a comprehensive understanding of how cells respond to metal stress at the systemic level, revealing the complex regulatory networks and evolutionary dynamics that drive microbial adaptation to metal-rich environments. In this article, we present a historical overview of the evolving understanding of the genetic determinants of metal resistance in microbes. Through multiple narrative threads, we illustrate how our knowledge of microbial metal resistance and genetics has interacted with genetic tools and concept development. This review also discusses how our understanding of microbial metal resistance has progressed from the Mendelian perspective to the current systems biology viewpoint, particularly as omics approaches have considerably enhanced our understanding. This system-level understanding has opened new possibilities for genetically engineered microorganisms to regulate metal homeostasis.

Novel heavy metal resistance gene clusters are present in the genome of Cupriavidus neocaledonicus STM 6070, a new species of Mimosa pudica microsymbiont isolated from heavy-metal-rich mining site soil

BACKGROUND

Cupriavidus strain STM 6070 was isolated from nickel-rich soil collected near Koniambo massif, New Caledonia, using the invasive legume trap host Mimosa pudica. STM 6070 is a heavy metal-tolerant strain that is highly effective at fixing nitrogen with M. pudica. Here we have provided an updated taxonomy for STM 6070 and described salient features of the annotated genome, focusing on heavy metal resistance (HMR) loci and heavy metal efflux (HME) systems.

RESULTS

The 6,771,773 bp high-quality-draft genome consists of 107 scaffolds containing 6118 protein-coding genes. ANI values show that STM 6070 is a new species of Cupriavidus. The STM 6070 symbiotic region was syntenic with that of the M. pudica-nodulating Cupriavidus taiwanensis LMG 19424T. In contrast to the nickel and zinc sensitivity of C. taiwanensis strains, STM 6070 grew at high Ni2+ and Zn2+ concentrations. The STM 6070 genome contains 55 genes, located in 12 clusters, that encode HMR structural proteins belonging to the RND, MFS, CHR, ARC3, CDF and P-ATPase protein superfamilies. These HMR molecular determinants are putatively involved in arsenic (ars), chromium (chr), cobalt-zinc-cadmium (czc), copper (cop, cup), nickel (nie and nre), and silver and/or copper (sil) resistance. Seven of these HMR clusters were common to symbiotic and non-symbiotic Cupriavidus species, while four clusters were specific to STM 6070, with three of these being associated with insertion sequences. Within the specific STM 6070 HMR clusters, three novel HME-RND systems (nieIC cep nieBA, czcC2B2A2, and hmxB zneAC zneR hmxS) were identified, which constitute new candidate genes for nickel and zinc resistance.

CONCLUSIONS

STM 6070 belongs to a new Cupriavidus species, for which we have proposed the name Cupriavidus neocaledonicus sp. nov.. STM6070 harbours a pSym with a high degree of gene conservation to the pSyms of M. pudica-nodulating C. taiwanensis strains, probably as a result of recent horizontal transfer. The presence of specific HMR clusters, associated with transposase genes, suggests that the selection pressure of the New Caledonian ultramafic soils has driven the specific adaptation of STM 6070 to heavy-metal-rich soils via horizontal gene transfer.

Genome-wide analysis of horizontally acquired genes in the genus Mycobacterium

Horizontal gene transfer (HGT) was attributed as a major driving force for the innovation and evolution of prokaryotic genomes. Previously, multiple research endeavors were undertaken to decipher HGT in different bacterial lineages. The genus Mycobacterium houses some of the most deadly human pathogens; however, the impact of HGT in Mycobacterium has never been addressed in a systematic way. Previous initiatives to explore the genomic imprints of HGTs in Mycobacterium were focused on few selected species, specifically among the members of Mycobacterium tuberculosis complex. Considering the recent availability of a large number of genomes, the current study was initiated to decipher the probable events of HGTs among 109 completely sequenced Mycobacterium species. Our comprehensive phylogenetic analysis with more than 9,000 families of Mycobacterium proteins allowed us to list several instances of gene transfers spread across the Mycobacterium phylogeny. Moreover, by examining the topology of gene phylogenies here, we identified the species most likely to donate and receive these genes and provided a detailed overview of the putative functions these genes may be involved in. Our study suggested that horizontally acquired foreign genes had played an enduring role in the evolution of Mycobacterium genomes and have contributed to their metabolic versatility and pathogenicity.

Strategies for high-altitude adaptation revealed from high-quality draft genome of non-violacein producing Janthinobacterium lividum ERGS5:01

A light pink coloured bacterial strain ERGS5:01 isolated from glacial stream water of Sikkim Himalaya was affiliated to Janthinobacterium lividum based on 16S rRNA gene sequence identity and phylogenetic clustering. Whole genome sequencing was performed for the strain to confirm its taxonomy as it lacked the typical violet pigmentation of the genus and also to decipher its survival strategy at the aquatic ecosystem of high elevation. The PacBio RSII sequencing generated genome of 5,168,928 bp with 4575 protein-coding genes and 118 RNA genes. Whole genome-based multilocus sequence analysis clustering, in silico DDH similarity value of 95.1% and, the ANI value of 99.25% established the identity of the strain ERGS5:01 (MCC 2953) as a non-violacein producing J. lividum. The genome comparisons across genus Janthinobacterium revealed an open pan-genome with the scope of the addition of new orthologous cluster to complete the genomic inventory. The genomic insight provided the genetic basis of freezing and frequent freeze-thaw cycle tolerance and, for industrially important enzymes. Extended insight into the genome provided clues of crucial genes associated with adaptation in the harsh aquatic ecosystem of high altitude.

Physiological and genomic insights into the lifestyle of arsenite-oxidizing Herminiimonas arsenitoxidans

Arsenic, a representative toxic metalloid, is responsible for serious global health problems. Most organisms possess arsenic resistance strategies to mitigate this toxicity. Here, we reported a microorganism, strain AS8, from heavy metal/metalloid-contaminated soil that is able to oxidize arsenite, and investigated its physiological and genomic traits. Its cells were rod-shaped and Gram-negative, and formed small beige-pigmented colonies. 16S rRNA-based phylogenetic analysis indicated that the strain belongs to the genus Herminiimonas and is closely related to Herminiimonas glaciei UMB49^T (98.7% of 16S rRNA gene sequence similarity), Herminiimonas arsenicoxydans ULPAs1^T (98.4%), and Herminiimonas saxobsidens NS11^T (98.4%). Under chemolithoheterotrophic conditions, the strain utilized some organic acids and amino acids as carbon and/or nitrogen sources but not electron sources. Further, the strain grew as a sulfur oxidizer in a complex medium (trypticase soy agar). Unexpectedly, most carbohydrates failed to support its growth as sole carbon sources. Genome sequencing supported these observations, and very few ABC transporters capable of oligo/monosaccharide uptake were identified in the AS8 genome. The genome harbored genes required for the colonization, flagella biosynthesis, urea degradation, and heavy metal and antibiotic resistance. Based on these polyphasic and genomic analyses, we propose that the strain AS8 be named Herminiimonas arsenitoxidans.

Saudi Moumouvirus, the First Group B Mimivirus Isolated from Asia

The number of novel giant viruses identified and characterized from the recently proposed order Megavirales has increased in recent years and new questions have been raised regarding viral diversity and evolution. Here, we describe the isolation and characterization of Saudi moumouvirus (SDMV), a new giant virus belonging to Mimivirus lineage B, isolated from a sewage sample taken from the King Abdulaziz University hospital in Jeddah, Saudi Arabia. SDMV presented 500 nm icosahedral particles with a 1,046,087 bp genome, which is larger than moumouvirus-like genomes which have been described in the past. The SDMV genome was predicted to encode 868 ORFs, ranging in size from 54 to 2,914 amino acids, with a mean size of 349 aa. Furthermore, this genome was predicted to encode 40 new genes (ORFans) without similarity with other sequences (ORFan L850 transcript was detected by qPCR in infected amoeba), in addition to 42 hypothetical proteins (pseudo-ORFs) with less than 100 aa, which matched other sequences in the NCBI nr database. Phylogenetic analysis showed that SDMV clustered together with mimiviruses from lineage B, including moumouvirus-like strains. It is, therefore, the third Mimivirus to be isolated in Asia and the first of group B.

Molecular Keys to the Janthinobacterium and Duganella spp. Interaction with the Plant Pathogen Fusarium graminearum

Janthinobacterium and Duganella are well-known for their antifungal effects. Surprisingly, almost nothing is known on molecular aspects involved in the close bacterium-fungus interaction. To better understand this interaction, we established the genomes of 11 Janthinobacterium and Duganella isolates in combination with phylogenetic and functional analyses of all publicly available genomes. Thereby, we identified a core and pan genome of 1058 and 23,628 genes. All strains encoded secondary metabolite gene clusters and chitinases, both possibly involved in fungal growth suppression. All but one strain carried a single gene cluster involved in the biosynthesis of alpha-hydroxyketone-like autoinducer molecules, designated JAI-1. Genome-wide RNA-seq studies employing the background of two isolates and the corresponding JAI-1 deficient strains identified a set of 45 QS-regulated genes in both isolates. Most regulated genes are characterized by a conserved sequence motif within the promoter region. Among the most strongly regulated genes were secondary metabolite and type VI secretion system gene clusters. Most intriguing, co-incubation studies of J. sp. HH102 or its corresponding JAI-1 synthase deletion mutant with the plant pathogen Fusarium graminearum provided first evidence of a QS-dependent interaction with this pathogen.

Host-Associated Genomic Features of the Novel Uncultured Intracellular Pathogen Ca. Ichthyocystis Revealed by Direct Sequencing of Epitheliocysts

Advances in single-cell and mini-metagenome sequencing have enabled important investigations into uncultured bacteria. In this study, we applied the mini-metagenome sequencing method to assemble genome drafts of the uncultured causative agents of epitheliocystis, an emerging infectious disease in the Mediterranean aquaculture species gilthead seabream. We sequenced multiple cyst samples and constructed 11 genome drafts from a novel beta-proteobacterial lineage, Candidatus Ichthyocystis. The draft genomes demonstrate features typical of pathogenic bacteria with an obligate intracellular lifestyle: a reduced genome of up to 2.6 Mb, reduced G + C content, and reduced metabolic capacity. Reconstruction of metabolic pathways reveals that Ca Ichthyocystis genomes lack all amino acid synthesis pathways, compelling them to scavenge from the fish host. All genomes encode type II, III, and IV secretion systems, a large repertoire of predicted effectors, and a type IV pilus. These are all considered to be virulence factors, required for adherence, invasion, and host manipulation. However, no evidence of lipopolysaccharide synthesis could be found. Beyond the core functions shared within the genus, alignments showed distinction into different species, characterized by alternative large gene families. These comprise up to a third of each genome, appear to have arisen through duplication and diversification, encode many effector proteins, and are seemingly critical for virulence. Thus, Ca Ichthyocystis represents a novel obligatory intracellular pathogenic beta-proteobacterial lineage. The methods used: mini-metagenome analysis and manual annotation, have generated important insights into the lifestyle and evolution of the novel, uncultured pathogens, elucidating many putative virulence factors including an unprecedented array of novel gene families.

Evolution of a Heavy Metal Homeostasis/Resistance Island Reflects Increasing Copper Stress in Enterobacteria

Copper homeostasis in bacteria is challenged by periodic elevation of copper levels in the environment, arising from both natural sources and human inputs. Several mechanisms have evolved to efflux copper from bacterial cells, including thecus(copper sensing copper efflux system), andpco(plasmid-borne copper resistance system) systems. The genes belonging to these two systems can be physically clustered in a Copper Homeostasis and Silver Resistance Island (CHASRI) on both plasmids and chromosomes in Enterobacteria. Increasing use of copper in agricultural and industrial applications raises questions about the role of human activity in the evolution of novel copper resistance mechanisms. Here we present evidence that CHASRI emerged and diversified in response to copper deposition across aerobic and anaerobic environments. An analysis of diversification rates and a molecular clock model suggest that CHASRI experienced repeated episodes of elevated diversification that could correspond to peaks in human copper production. Phylogenetic analyses suggest that CHASRI originated in a relative ofEnterobacter cloacaeas the ultimate product of sequential assembly of several pre-existing two-gene modules. Once assembled, CHASRI dispersed via horizontal gene transfer within Enterobacteriaceae and also to certain members of Shewanellaceae, where the originalpcomodule was replaced by a divergentpcohomolog. Analyses of copper stress mitigation suggest that CHASRI confers increased resistance aerobically, anaerobically, and during shifts between aerobic and anaerobic environments, which could explain its persistence in facultative anaerobes and emergent enteric pathogens.

Draft genome of Gemmata massiliana sp. nov, a water-borne Planctomycetes species exhibiting two variants

Gemmata massiliana is a new Planctomycetes bacterium isolated from a hospital water network in France, using a new culture medium. It is an aerobic microorganism with optimal growth at pH 8, at 30 °C and salinity ≤ 1.25 % NaCl. G. massiliana is resistant to β-lactam antibiotics, due to lack of peptidoglycan in its cell wall.G. massiliana shares a 97 % 16S rRNA gene sequence similarity with the nearest species, Gemmata obscuriglobus; and 99 % similarity with unnamed soil isolates. Its 9,249,437-bp genome consists in one chromosome and no detectable plasmid and has a 64.07 % G + C content, 32.94 % of genes encoding for hypothetical proteins. The genome contains an incomplete 19.6-kb phage sequence, 26 CRISPRs, 3 CAS and 15 clusters of secondary metabolites. G. massiliana genome increases knowledge of a poorly known world of bacteria.

Pan-genomic analysis to redefine species and subspecies based on quantum discontinuous variation: the Klebsiella paradigm

Background

Various methods are currently used to define species and are based on the phylogenetic marker 16S ribosomal RNA gene sequence, DNA-DNA hybridization and DNA GC content. However, these are restricted genetic tools and showed significant limitations.

Results

In this work, we describe an alternative method to build taxonomy by analyzing the pan-genome composition of different species of the Klebsiella genus. Klebsiella species are Gram-negative bacilli belonging to the large Enterobacteriaceae family. Interestingly, when comparing the core/pan-genome ratio; we found a clear discontinuous variation that can define a new species.

Conclusions

Using this pan-genomic approach, we showed that Klebsiella pneumoniae subsp. ozaenae and Klebsiella pneumoniae subsp. rhinoscleromatis are species of the Klebsiella genus, rather than subspecies of Klebsiella pneumoniae. This pan-genomic analysis, helped to develop a new tool for defining species introducing a quantic perspective for taxonomy.

Reviewers

This article was reviewed by William Martin, Pierre Pontarotti and Pere Puigbo (nominated by Dr Yuri Wolf).

Electronic supplementary material

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

Genome Sequence of the Soil Bacterium Janthinobacterium sp. KBS0711

We present a draft genome of Janthinobacterium sp. KBS0711 that was isolated from agricultural soil. The genome provides insight into the ecological strategies of this bacterium in free-living and host-associated environments.

Genome Sequence of Mushroom Soft-Rot Pathogen Janthinobacterium agaricidamnosum

Janthinobacterium agaricidamnosum causes soft-rot disease of the cultured button mushroom Agaricus bisporus and is thus responsible for agricultural losses. Here, we present the genome sequence of J. agaricidamnosum DSM 9628. The 5.9-Mb genome harbors several secondary metabolite biosynthesis gene clusters, which renders this neglected bacterium a promising source for genome mining approaches.

Looking at protists as a source of pathogenic viruses

In the environment, protozoa are predators of bacteria and feed on them. The possibility that some protozoa could be a source of human pathogens is consistent with the discovery that free-living amoebae were the reservoir of Legionella pneumophila, the agent of Legionnaires' disease. Later, while searching for Legionella in the environment using amoeba co-culture, the first giant virus, Acanthamoeba polyphaga mimivirus, was discovered. Since then, many other giant viruses have been isolated, including Marseilleviridae, Pithovirus sibericum, Cafeteria roenbergensis virus and Pandoravirus spp. The methods used to isolate all of these viruses are herein reviewed. By analogy to Legionella, it was originally suspected that these viruses could be human pathogens. After showing by indirect evidence, such as sero-epidemiologic studies, that it was possible for these viruses to be human pathogens, the recent isolation of some of these viruses (belonging to the Mimiviridae and Marseilleviridae families) in humans in the context of pathologic conditions shows that they are opportunistic human pathogens in some instances.

The gastrointestinal tract of the white-throated Woodrat (Neotoma albigula) harbors distinct consortia of oxalate-degrading bacteria

The microbiota inhabiting the mammalian gut is a functional organ that provides a number of services for the host. One factor that may regulate the composition and function of gut microbial communities is dietary toxins. Oxalate is a toxic plant secondary compound (PSC) produced in all major taxa of vascular plants and is consumed by a variety of animals. The mammalian herbivore Neotoma albigula is capable of consuming and degrading large quantities of dietary oxalate. We isolated and characterized oxalate-degrading bacteria from the gut contents of wild-caught animals and used high-throughput sequencing to determine the distribution of potential oxalate-degrading taxa along the gastrointestinal tract. Isolates spanned three genera: Lactobacillus, Clostridium, and Enterococcus. Over half of the isolates exhibited significant oxalate degradation in vitro, and all Lactobacillus isolates contained the oxc gene, one of the genes responsible for oxalate degradation. Although diverse potential oxalate-degrading genera were distributed throughout the gastrointestinal tract, they were most concentrated in the foregut, where dietary oxalate first enters the gastrointestinal tract. We hypothesize that unique environmental conditions present in each gut region provide diverse niches that select for particular functional taxa and communities.

A Markovian analysis of bacterial genome sequence constraints

The arrangement of nucleotides within a bacterial chromosome is influenced by numerous factors. The degeneracy of the third codon within each reading frame allows some flexibility of nucleotide selection; however, the third nucleotide in the triplet of each codon is at least partly determined by the preceding two. This is most evident in organisms with a strong G + C bias, as the degenerate codon must contribute disproportionately to maintaining that bias. Therefore, a correlation exists between the first two nucleotides and the third in all open reading frames. If the arrangement of nucleotides in a bacterial chromosome is represented as a Markov process, we would expect that the correlation would be completely captured by a second-order Markov model and an increase in the order of the model (e.g., third-, fourth-…order) would not capture any additional uncertainty in the process. In this manuscript, we present the results of a comprehensive study of the Markov property that exists in the DNA sequences of 906 bacterial chromosomes. All of the 906 bacterial chromosomes studied exhibit a statistically significant Markov property that extends beyond second-order, and therefore cannot be fully explained by codon usage. An unrooted tree containing all 906 bacterial chromosomes based on their transition probability matrices of third-order shares ∼25% similarity to a tree based on sequence homologies of 16S rRNA sequences. This congruence to the 16S rRNA tree is greater than for trees based on lower-order models (e.g., second-order), and higher-order models result in diminishing improvements in congruence. A nucleotide correlation most likely exists within every bacterial chromosome that extends past three nucleotides. This correlation places significant limits on the number of nucleotide sequences that can represent probable bacterial chromosomes. Transition matrix usage is largely conserved by taxa, indicating that this property is likely inherited, however some important exceptions exist that may indicate the convergent evolution of some bacteria.

Genome expression profile analysis reveals important transcripts in maize roots responding to the stress of heavy metal Pb

Lead (Pb) has become one of the most abundant heavy metal pollutants of the environment. With its large biomass, maize could be an important object for studying the phytoremediation of Pb-contaminated soil. In our previous research, we screened 19 inbred lines of maize for Pb concentration, and line 178 was identified to be a hyperaccumulator for Pb in both the roots and aboveground parts. To identify important genes and metabolic pathways related to Pb accumulation and tolerance, line 178 was underwent genome expression profile under Pb stress and a control (CK). A total of approximately 11 million cDNA tags were sequenced and 4 665 539 and 4 936 038 clean tags were obtained from the libraries of the test and CK, respectively. In comparison to CK, 2379 and 1832 genes were identified up- or downregulated, respectively, more than fivefolds under Pb stress. Interestingly, all the genes were related to cellular processes and signaling, information storage and processing or metabolism functions. Particularly, the genes involved in posttranslational modification, protein turnover and chaperones; signal transduction, carbohydrate transport and metabolism; and lipid transport and metabolism significantly changed under the treatment. In addition, seven pathways including ribosome, photosynthesis, and carbon fixation were affected significantly, with 118, 12, 34, 21, 18, 72 and 43 differentially expressed genes involved. The significant upregulation of the ribosome pathway may reveal an important secret for Pb tolerance of line 178. And the sharp increase of laccase transcripts and metal ion transporters were suggested to account in part for Pb hyperaccumulation in the line.

The Janthinobacterium sp. HH01 genome encodes a homologue of the V. cholerae CqsA and L. pneumophila LqsA autoinducer synthases

Janthinobacteria commonly form biofilms on eukaryotic hosts and are known to synthesize antibacterial and antifungal compounds. Janthinobacterium sp. HH01 was recently isolated from an aquatic environment and its genome sequence was established. The genome consists of a single chromosome and reveals a size of 7.10 Mb, being the largest janthinobacterial genome so far known. Approximately 80% of the 5,980 coding sequences (CDSs) present in the HH01 genome could be assigned putative functions. The genome encodes a wealth of secretory functions and several large clusters for polyketide biosynthesis. HH01 also encodes a remarkable number of proteins involved in resistance to drugs or heavy metals. Interestingly, the genome of HH01 apparently lacks the N-acylhomoserine lactone (AHL)-dependent signaling system and the AI-2-dependent quorum sensing regulatory circuit. Instead it encodes a homologue of the Legionella- and Vibrio-like autoinducer (lqsA/cqsA) synthase gene which we designated jqsA. The jqsA gene is linked to a cognate sensor kinase (jqsS) which is flanked by the response regulator jqsR. Here we show that a jqsA deletion has strong impact on the violacein biosynthesis in Janthinobacterium sp. HH01 and that a jqsA deletion mutant can be functionally complemented with the V. cholerae cqsA and the L. pneumophila lqsA genes.

Parallel changes in the taxonomical structure of bacterial communities exposed to a similar environmental disturbance

Bacterial communities play a central role in ecosystems, by regulating biogeochemical fluxes. Therefore, understanding how multiple functional interactions between species face environmental perturbations is a major concern in conservation biology. Because bacteria can use several strategies, including horizontal gene transfers (HGT), to cope with rapidly changing environmental conditions, potential decoupling between function and taxonomy makes the use of a given species as a general bioindicator problematic. The present work is a first step to characterize the impact of a recent polymetallic gradient over the taxonomical networks of five lacustrine bacterial communities. Given that evolutionary convergence represents one of the best illustration of natural selection, we focused on a system composed of two pairs of impacted and clean lakes in order to test whether similar perturbation exerts a comparable impact on the taxonomical networks of independent bacterial communities. First, we showed that similar environmental stress drove parallel structural changes at the taxonomic level on two independent bacterial communities. Second, we showed that a long-term exposure to contaminant gradients drove significant taxonomic structure changes within three interconnected bacterial communities. Thus, this model lake system is relevant to characterize the strategies, namely acclimation and/or adaptation, of bacterial communities facing environmental perturbations, such as metal contamination.

Comparative genomics of bacteria from the genus Collimonas: linking (dis)similarities in gene content to phenotypic variation and conservation

Collimonas is a genus of soil bacteria comprising three recognized species: C. fungivorans, C. pratensis and C. arenae. Collimonads share the ability to degrade chitin (chitinolysis), feed on living fungal hyphae (mycophagy), and dissolve minerals (weathering), but vary in their inhibition of fungi (fungistasis). To better understand this phenotypic variability, we analysed the genomic content of four strains representing three Collimonas species (Ter14, Ter6, Ter91 and Ter10) by hybridization to a microarray based on reference strain C. fungivorans Ter331. The analysis revealed genes unique to strain Ter331 (e.g. those on the extrachromosomal element pTer331) and genes present in some but not all of the tested strains. Among the latter were several candidates that may contribute to fungistasis, including genes for the production and secretion of antifungals. We hypothesize that differential possession of these genes underlies the specialization of Collimonas strains towards different fungal hosts. We identified a set of 136 genes that were common in all tested Collimonas strains, but absent from the genomes of three other members of the family Oxalobacteraceae. Predicted products of these 'Collimonas core' genes include lytic, secreted enzymes such as chitinases, peptidases, nucleases and phosphatases with a putative role in mycophagy and weathering.

Comparative metagenomics of two microbial mats at Cuatro Ciénegas Basin II: community structure and composition in oligotrophic environments

Microbial mats are self-sustained, functionally complex ecosystems that make good models for the understanding of past and present microbial ecosystems as well as putative extraterrestrial ecosystems. Ecological theory suggests that the composition of these communities might be affected by nutrient availability and disturbance frequency. We characterized two microbial mats from two contrasting environments in the oligotrophic Cuatro Ciénegas Basin: a permanent green pool and a red desiccation pond. We analyzed their taxonomic structure and composition by means of 16S rRNA clone libraries and metagenomics and inferred their metabolic role by the analysis of functional traits in the most abundant organisms. Both mats showed a high diversity with metabolically diverse members and strongly differed in structure and composition. The green mat had a higher species richness and evenness than the red mat, which was dominated by a lineage of Pseudomonas. Autotrophs were abundant in the green mat, and heterotrophs were abundant in the red mat. When comparing with other mats and stromatolites, we found that taxonomic composition was not shared at species level but at order level, which suggests environmental filtering for phylogenetically conserved functional traits with random selection of particular organisms. The highest diversity and composition similarity was observed among systems from stable environments, which suggests that disturbance regimes might affect diversity more strongly than nutrient availability, since oligotrophy does not appear to prevent the establishment of complex and diverse microbial mat communities. These results are discussed in light of the search for extraterrestrial life.

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.

Chromosome-encoded extended-spectrum class A β-lactamase MIN-1 from Minibacterium massiliensis

Minibacterium massiliensis strain CIP107820 is a recently discovered waterborne Gram-negative rod isolated from hospital water samples. It harbors a chromosomally located gene encoding an Ambler class A extended-spectrum β-lactamase termed MIN-1, sharing 56%, 54%, and 51% amino acid identities with β-lactamases LUT-1, KPC-2, and CTX-M-2, respectively. β-Lactamase MIN-1 hydrolyzes penicillins, narrow-spectrum cephalosporins, cefotaxime, and, less efficiently, cefepime, while ceftazidime and carbapenems are very poor substrates, and cephamycins and aztreonam are not hydrolyzed.

The genealogic tree of mycobacteria reveals a long-standing sympatric life into free-living protozoa

Free-living protozoa allow horizontal gene transfer with and between the microorganisms that they host. They host mycobacteria for which the sources of transferred genes remain unknown. Using BLASTp, we searched within the genomes of 15 mycobacteria for homologous genes with 34 amoeba-resistant bacteria and the free-living protozoa Dictyostelium discoideum. Subsequent phylogenetic analysis of these sequences revealed that eight mycobacterial open-reading frames (ORFs) were probably acquired via horizontal transfer from beta- and gamma-Proteobacteria and from Firmicutes, but the transfer histories could not be reliably established in details. One further ORF encoding a pyridine nucleotide disulfide oxidoreductase (pyr-redox) placed non-tuberculous mycobacteria in a clade with Legionella spp., Francisella spp., Coxiella burnetii, the ciliate Tetrahymena thermophila and D. discoideum with a high reliability. Co-culturing Mycobacterium avium and Legionella pneumophila with the amoeba Acanthamoeba polyphaga demonstrated that these two bacteria could live together in amoebae for five days, indicating the biological relevance of intra-amoebal transfer of the pyr-redox gene. In conclusion, the results of this study support the hypothesis that protists can serve as a source and a place for gene transfer in mycobacteria.

The influence of rickettsiologists on post-modern microbiology

Many of the definitions in microbiology are currently false. We have reviewed the great denominations of microbiology and attempted to free microorganisms from the theories of the twentieth century. The presence of compartmentation and a nucleoid in Planctomycetes clearly calls into question the accuracy of the definitions of eukaryotes and prokaryotes. Archaea are viewed as prokaryotes resembling bacteria. However, the name archaea, suggesting an archaic origin of lifestyle, is inconsistent with the lifestyle of this family. Viruses are defined as small, filterable infectious agents, but giant viruses challenge the size criteria used for the definition of a virus. Pathogenicity does not require the acquisition of virulence factors (except for toxins), and in many cases, gene loss is significantly inked to the emergence of virulence. Species classification based on 16S rRNA is useless for taxonomic purposes of human pathogens, as a 2% divergence would classify all Rickettsiae within the same species and would not identify bacteria specialized for mammal infection. The use of metagenomics helps us to understand evolution and physiology by elucidating the structure, function, and interactions of the major microbial communities, but it neglects the minority populations. Finally, Darwin’s descent with modification theory, as represented by the tree of life, no longer matches our current genomic knowledge because genomics has revealed the occurrence of de novo-created genes and the mosaic structure of genomes, the Rhizome of life is therefore more appropriate.

Genomes of the most dangerous epidemic bacteria have a virulence repertoire characterized by fewer genes but more toxin-antitoxin modules

BACKGROUND

We conducted a comparative genomic study based on a neutral approach to identify genome specificities associated with the virulence capacity of pathogenic bacteria. We also determined whether virulence is dictated by rules, or if it is the result of individual evolutionary histories. We systematically compared the genomes of the 12 most dangerous pandemic bacteria for humans ("bad bugs") to their closest non-epidemic related species ("controls").

METHODOLOGY/PRINCIPAL FINDINGS

We found several significantly different features in the "bad bugs", one of which was a smaller genome that likely resulted from a degraded recombination and repair system. The 10 Cluster of Orthologous Group (COG) functional categories revealed a significantly smaller number of genes in the "bad bugs", which lacked mostly transcription, signal transduction mechanisms, cell motility, energy production and conversion, and metabolic and regulatory functions. A few genes were identified as virulence factors, including secretion system proteins. Five "bad bugs" showed a greater number of poly (A) tails compared to the controls, whereas an elevated number of poly (A) tails was found to be strongly correlated to a low GC% content. The "bad bugs" had fewer tandem repeat sequences compared to controls. Moreover, the results obtained from a principal component analysis (PCA) showed that the "bad bugs" had surprisingly more toxin-antitoxin modules than did the controls.

CONCLUSIONS/SIGNIFICANCE

We conclude that pathogenic capacity is not the result of "virulence factors" but is the outcome of a virulent gene repertoire resulting from reduced genome repertoires. Toxin-antitoxin systems could participate in the virulence repertoire, but they may have developed independently of selfish evolution.

Defining pathogenic bacterial species in the genomic era

Actual definitions of bacterial species are limited due to the current criteria of definition and the use of restrictive genetic tools. The 16S ribosomal RNA sequence, for example, has been widely used as a marker for phylogenetic analyses; however, its use often leads to misleading species definitions. According to the first genetic studies, removing a certain number of genes from pathogenic bacteria removes their capacity to infect hosts. However, more recent studies have demonstrated that the specialization of bacteria in eukaryotic cells is associated with massive gene loss, especially for allopatric endosymbionts that have been isolated for a long time in an intracellular niche. Indeed, sympatric free-living bacteria often have bigger genomes and exhibit greater resistance and plasticity and constitute species complexes rather than true species. Specialists, such as pathogenic bacteria, escape these bacterial complexes and colonize a niche, thereby gaining a species name. Their specialization allows them to become allopatric, and their gene losses eventually favor reductive genome evolution. A pathogenic species is characterized by a gene repertoire that is defined not only by genes that are present but also by those that are lacking. It is likely that current bacterial pathogens will disappear soon and be replaced by new ones that will emerge from bacterial complexes that are already in contact with humans.

Functional convergence in reduced genomes of bacterial symbionts spanning 200 My of evolution

The main genomic changes in the evolution of host-restricted microbial symbionts are ongoing inactivation and loss of genes combined with rapid sequence evolution and extreme structural stability; these changes reflect high levels of genetic drift due to small population sizes and strict clonality. This genomic erosion includes irreversible loss of genes in many functional categories and can include genes that underlie the nutritional contributions to hosts that are the basis of the symbiotic association. Candidatus Sulcia muelleri is an ancient symbiont of sap-feeding insects and is typically coresident with another bacterial symbiont that varies among host subclades. Previously sequenced Sulcia genomes retain pathways for the same eight essential amino acids, whereas coresident symbionts synthesize the remaining two. Here, we describe a dual symbiotic system consisting of Sulcia and a novel species of Betaproteobacteria, Candidatus Zinderia insecticola, both living in the spittlebug Clastoptera arizonana. This Sulcia has completely lost the pathway for the biosynthesis of tryptophan and, therefore, retains the ability to make only 7 of the 10 essential amino acids. Zinderia has a tiny genome (208 kb) and the most extreme nucleotide base composition (13.5% G + C) reported to date, yet retains the ability to make the remaining three essential amino acids, perfectly complementing capabilities of the coresident Sulcia. Combined with the results from related symbiotic systems with complete genomes, these data demonstrate the critical role that bacterial symbionts play in the host insect’s biology and reveal one outcome following the loss of a critical metabolic activity through genome reduction.

Mycobacterium tuberculosis Rv1090 and Rv1987 encode functional β-glucan-targeting proteins

Mycobacterium tuberculosis is a facultative intracellular pathogen, and the ability of this bacterium to survive and to grow inside macrophages is central to its virulence. Multiple strategies are employed by M. tuberculosis to ensure survival in macrophages, including secretion of several proteins, which are good candidates to be virulence factors, drug targets for disease intervention, and vaccine antigens. However, some M. tuberculosis secreted proteins do not appear to play any role in the growth or survival of the bacterium in its mammalian host. Among these proteins are three putative cellulose-targeting proteins encoded by the genes Rv0062, Rv1090, and Rv1987. It has been previously shown that Rv0062 encodes an active cellulase. Here we report that Rv1090 and Rv1987 also encode functional proteins. Rv1090 is able to hydrolyze barley β-glucan while Rv1987 displays cellulose-binding activity on filter paper and on microcrystalline cellulose (Avicel). Collectively, these observations point toward a unique unknown relationship between M. tuberculosis and a cellulose-containing host. We hypothesize that amoeba could be such hosts.

The complete genome sequence of Cupriavidus metallidurans strain CH34, a master survivalist in harsh and anthropogenic environments

Many bacteria in the environment have adapted to the presence of toxic heavy metals. Over the last 30 years, this heavy metal tolerance was the subject of extensive research. The bacterium Cupriavidus metallidurans strain CH34, originally isolated by us in 1976 from a metal processing factory, is considered a major model organism in this field because it withstands milli-molar range concentrations of over 20 different heavy metal ions. This tolerance is mostly achieved by rapid ion efflux but also by metal-complexation and -reduction. We present here the full genome sequence of strain CH34 and the manual annotation of all its genes. The genome of C. metallidurans CH34 is composed of two large circular chromosomes CHR1 and CHR2 of, respectively, 3,928,089 bp and 2,580,084 bp, and two megaplasmids pMOL28 and pMOL30 of, respectively, 171,459 bp and 233,720 bp in size. At least 25 loci for heavy-metal resistance (HMR) are distributed over the four replicons. Approximately 67% of the 6,717 coding sequences (CDSs) present in the CH34 genome could be assigned a putative function, and 9.1% (611 genes) appear to be unique to this strain. One out of five proteins is associated with either transport or transcription while the relay of environmental stimuli is governed by more than 600 signal transduction systems. The CH34 genome is most similar to the genomes of other Cupriavidus strains by correspondence between the respective CHR1 replicons but also displays similarity to the genomes of more distantly related species as a result of gene transfer and through the presence of large genomic islands. The presence of at least 57 IS elements and 19 transposons and the ability to take in and express foreign genes indicates a very dynamic and complex genome shaped by evolutionary forces. The genome data show that C. metallidurans CH34 is particularly well equipped to live in extreme conditions and anthropogenic environments that are rich in metals.

Giant Marseillevirus highlights the role of amoebae as a melting pot in emergence of chimeric microorganisms

Giant viruses such as Mimivirus isolated from amoeba found in aquatic habitats show biological sophistication comparable to that of simple cellular life forms and seem to evolve by similar mechanisms, including extensive gene duplication and horizontal gene transfer (HGT), possibly in part through a viral parasite, the virophage. We report here the isolation of "Marseille" virus, a previously uncharacterized giant virus of amoeba. The virions of Marseillevirus encompass a 368-kb genome, a minimum of 49 proteins, and some messenger RNAs. Phylogenetic analysis of core genes indicates that Marseillevirus is the prototype of a family of nucleocytoplasmic large DNA viruses (NCLDV) of eukaryotes. The genome repertoire of the virus is composed of typical NCLDV core genes and genes apparently obtained from eukaryotic hosts and their parasites or symbionts, both bacterial and viral. We propose that amoebae are "melting pots" of microbial evolution where diverse forms emerge, including giant viruses with complex gene repertoires of various origins.

Herminiimonas glaciei sp. nov., a novel ultramicrobacterium from 3042 m deep Greenland glacial ice

A Gram-negative ultramicrobacterium (designated strain UMB49(T)) was isolated from a 120,000-year-old, 3,042 m deep Greenland glacier ice core using a 0.2 mum filtration enrichment procedure. Phylogenetic analysis of the 16S rRNA gene sequence indicated that this strain belonged to the genus Herminiimonas of the family Oxalobacteraceae of the class Betaproteobacteria. Strain UMB49(T) was most closely related to Herminiimonas saxobsidens (99.6 % sequence similarity), Herminiimonas arsenicoxydans (98.4 %), Herminiimonas aquatilis (97.6 %) and Herminiimonas fonticola (97.9 %). Genomic DNA-DNA hybridization showed low levels of relatedness (below 57 %) to H. saxobsidens and H. arsenicoxydans. Cells of strain UMB49(T) were small thin rods with a mean volume of 0.043 mum(3) and possessed 1 or 2 polar and/or 1-3 lateral very long flagella. The original colony pigmentation was brown-purple but after recultivation the colonies were translucent white to tan coloured. Strain UMB49(T) grew aerobically and under microaerophilic conditions. The strain produced catalase and oxidase, but did not reduce nitrate. Sole carbon sources included citrate, succinate, malate, lactate and alanine. The strain produced acid from l-arabinose, d-arabinose, l-xylose, d-xylose and d-ribose. The DNA G+C content was 59.0 mol%. Based on differential characteristics of strain UMB49(T) and recognized Herminiimonas species, it was concluded that strain UMB49(T) represents a novel species of the genus Herminiimonas, for which the name Herminiimonas glaciei sp. nov. is proposed. The type strain is UMB49(T) (=ATCC BAA-1623(T)=DSM 21140(T)).

Stenotrophomonas maltophilia genomes: a start-up comparison

The whole DNA sequences of 2 Stenotrophomonas maltophilia strains isolated from the blood of a cancer patient (K279a) and the poplar Populus trichocarpa (R551-3) have been compared. The 2 chromosomes exhibit extensive synteny, but each is punctuated by about 40 genomic islands (GEIs), which vary in size from 3 to 70kb, and may encode up to about 50 proteins. A large set of smaller DNA sequences, encoding strain-specific 'solo' orfs, contributes to genetic heterogeneity in a significant manner. S. maltophilia GEIs potentially encode several proteins mediating interactions with the environment such as transmembrane proteins, haemagglutinins, components of type I and IV secretion systems, and efflux proteins having a role in metal and/or drug resistance. The presence of specific GEIs in the S. maltophilia population was monitored by PCR and slot-blot analyses. Data suggest that some islands are present at sites different from those identified in K279a and that alternative islands may be integrated at mapped sites.

Massive comparative genomic analysis reveals convergent evolution of specialized bacteria

BACKGROUND

Genome size and gene content in bacteria are associated with their lifestyles. Obligate intracellular bacteria (i.e., mutualists and parasites) have small genomes that derived from larger free-living bacterial ancestors; however, the different steps of bacterial specialization from free-living to intracellular lifestyle have not been studied comprehensively. The growing number of available sequenced genomes makes it possible to perform a statistical comparative analysis of 317 genomes from bacteria with different lifestyles.

RESULTS

Compared to free-living bacteria, host-dependent bacteria exhibit fewer rRNA genes, more split rRNA operons and fewer transcriptional regulators, linked to slower growth rates. We found a function-dependent and non-random loss of the same 100 orthologous genes in all obligate intracellular bacteria. Thus, we showed that obligate intracellular bacteria from different phyla are converging according to their lifestyle. Their specialization is an irreversible phenomenon characterized by translation modification and massive gene loss, including the loss of transcriptional regulators. Although both mutualists and parasites converge by genome reduction, these obligate intracellular bacteria have lost distinct sets of genes in the context of their specific host associations: mutualists have significantly more genes that enable nutrient provisioning whereas parasites have genes that encode Types II, IV, and VI secretion pathways.

CONCLUSION

Our findings suggest that gene loss, rather than acquisition of virulence factors, has been a driving force in the adaptation of parasites to eukaryotic cells. This comparative genomic analysis helps to explore the strategies by which obligate intracellular genomes specialize to particular host-associations and contributes to advance our knowledge about the mechanisms of bacterial evolution.

The genome of Borrelia recurrentis, the agent of deadly louse-borne relapsing fever, is a degraded subset of tick-borne Borrelia duttonii

In an effort to understand how a tick-borne pathogen adapts to the body louse, we sequenced and compared the genomes of the recurrent fever agents Borrelia recurrentis and B. duttonii. The 1,242,163–1,574,910-bp fragmented genomes of B. recurrentis and B. duttonii contain a unique 23-kb linear plasmid. This linear plasmid exhibits a large polyT track within the promoter region of an intact variable large protein gene and a telomere resolvase that is unique to Borrelia. The genome content is characterized by several repeat families, including antigenic lipoproteins. B. recurrentis exhibited a 20.4% genome size reduction and appeared to be a strain of B. duttonii, with a decaying genome, possibly due to the accumulation of genomic errors induced by the loss of recA and mutS. Accompanying this were increases in the number of impaired genes and a reduction in coding capacity, including surface-exposed lipoproteins and putative virulence factors. Analysis of the reconstructed ancestral sequence compared to B. duttonii and B. recurrentis was consistent with the accelerated evolution observed in B. recurrentis. Vector specialization of louse-borne pathogens responsible for major epidemics was associated with rapid genome reduction. The correlation between gene loss and increased virulence of B. recurrentis parallels that of Rickettsia prowazekii, with both species being genomic subsets of less-virulent strains.

Borreliae are vector-borne spirochetes that are responsible for Lyme disease and recurrent fevers. We completed the genome sequences of the tick-borne Borrelia duttonii and the louse-borne B. recurrentis. The former of these is responsible for emerging infections that mimic malaria in Africa and in travellers, and the latter is responsible for severe recurrent fever in poor African populations. Diagnostic tools for these pathogens remain poor with regard to sensitivity and specificity due, in part, to the lack of genomic sequences. In this study, we show that the genomic content of B. recurrentis is a subset of that of B. duttonii, the genes of which are undergoing a decay process. These phenomena are common to all louse-borne pathogens compared to their tick-borne counterparts. In B. recurrentis, this process may be due to the inactivation of genes encoding DNA repair mechanisms, implying the accumulation of errors in the genome. The increased virulence of B. recurrentis could not be traced back to specific virulence factors, illustrating the lack of correlation between the virulence of a pathogen and so-called virulence genes. Knowledge of these genomes will allow for the development of new molecular tools that provide a more-accurate, sensitive, and specific diagnosis of these emerging infections.

Identification and characterization of genes underlying chitinolysis in Collimonas fungivorans Ter331

Through a combinatorial approach of plasposon mutagenesis, genome mining, and heterologous expression, we identified genes contributing to the chitinolytic phenotype of bacterium Collimonas fungivorans Ter331. One of five mutants with abolished ability to hydrolyze colloidal chitin carried its plasposon in the chiI gene coding for an extracellular endochitinase. Two mutants were affected in the promoter of chiP-II coding for an outer-membrane transporter of chitooligosaccharides. The remaining two mutations were linked to chitobiose/N-acetylglucosamine uptake. Thus, our model for the Collimonas chitinolytic system assumes a positive feedback regulation of chitinase activity by chitin degradation products. A second chitinase gene, chiII, coded for an exochitinase that preferentially released chitobiose from chitin analogs. Genes hexI and hexII showed coding resemblance to N-acetylglucosaminidases, and the activity of purified HexI protein towards chitin analogs suggested its role in converting chitobiose to N-acetylglucosamine. The hexI gene clustered with chiI, chiII, and chiP-II in one locus, while chitobiose/N-acetylglucosamine uptake genes colocalized in another. Both loci contained genes for conversion of N-acetylglucosamine to fructose-6-phosphate, confirming that C. fungivorans Ter331 features a complete chitin pathway. No link could be established between chitinolysis and antifungal activity of C. fungivorans Ter331, suggesting that the bacterium's reported antagonism towards fungi relies on other mechanisms.

The genome of Bacillus coahuilensis reveals adaptations essential for survival in the relic of an ancient marine environment

The Cuatro Ciénegas Basin (CCB) in the central part of the Chihuahan desert (Coahuila, Mexico) hosts a wide diversity of microorganisms contained within springs thought to be geomorphological relics of an ancient sea. A major question remaining to be answered is whether bacteria from CCB are ancient marine bacteria that adapted to an oligotrophic system poor in NaCl, rich in sulfates, and with extremely low phosphorus levels (<0.3 microM). Here, we report the complete genome sequence of Bacillus coahuilensis, a sporulating bacterium isolated from the water column of a desiccation lagoon in CCB. At 3.35 Megabases this is the smallest genome sequenced to date of a Bacillus species and provides insights into the origin, evolution, and adaptation of B. coahuilensis to the CCB environment. We propose that the size and complexity of the B. coahuilensis genome reflects the adaptation of an ancient marine bacterium to a novel environment, providing support to a "marine isolation origin hypothesis" that is consistent with the geology of CCB. This genomic adaptation includes the acquisition through horizontal gene transfer of genes involved in phosphorous utilization efficiency and adaptation to high-light environments. The B. coahuilensis genome sequence also revealed important ecological features of the bacterial community in CCB and offers opportunities for a unique glimpse of a microbe-dominated world last seen in the Precambrian.

De novo bacterial genome sequencing: millions of very short reads assembled on a desktop computer

Novel high-throughput DNA sequencing technologies allow researchers to characterize a bacterial genome during a single experiment and at a moderate cost. However, the increase in sequencing throughput that is allowed by using such platforms is obtained at the expense of individual sequence read length, which must be assembled into longer contigs to be exploitable. This study focuses on the Illumina sequencing platform that produces millions of very short sequences that are 35 bases in length. We propose a de novo assembler software that is dedicated to process such data. Based on a classical overlap graph representation and on the detection of potentially spurious reads, our software generates a set of accurate contigs of several kilobases that cover most of the bacterial genome. The assembly results were validated by comparing data sets that were obtained experimentally for Staphylococcus aureus strain MW2 and Helicobacter acinonychis strain Sheeba with that of their published genomes acquired by conventional sequencing of 1.5- to 3.0-kb fragments. We also provide indications that the broad coverage achieved by high-throughput sequencing might allow for the detection of clonal polymorphisms in the set of DNA molecules being sequenced.

Redefining viruses: lessons from Mimivirus

Viruses are the most abundant living entities and probably had a major role in the evolution of life, but are still defined using negative criteria. Here, we propose to divide biological entities into two groups of organisms: ribosome-encoding organisms, which include eukaryotic, archaeal and bacterial organisms, and capsid-encoding organisms, which include viruses. Other replicons (for example, plasmids and viroids) can be termed 'orphan replicons'. Based on this suggested classification system, we propose a new definition for a virus--a capsid-encoding organism that is composed of proteins and nucleic acids, self-assembles in a nucleocapsid and uses a ribosome-encoding organism for the completion of its life cycle.