Variable copy number, intra-genomic heterogeneities and lateral transfers of the 16S rRNA gene in Pseudomonas

Genomic Insights into Cyanide Biodegradation in the Pseudomonas Genus

Molecular studies about cyanide biodegradation have been mainly focused on the hydrolytic pathways catalyzed by the cyanide dihydratase CynD or the nitrilase NitC. In some Pseudomonas strains, the assimilation of cyanide has been linked to NitC, such as the cyanotrophic model strain Pseudomonas pseudoalcaligenes CECT 5344, which has been recently reclassified as Pseudomonas oleovorans CECT 5344. In this work, a phylogenomic approach established a more precise taxonomic position of the strain CECT 5344 within the species P. oleovorans. Furthermore, a pan-genomic analysis of P. oleovorans and other species with cyanotrophic strains, such as P. fluorescens and P. monteilii, allowed for the comparison and identification of the cioAB and mqoAB genes involved in cyanide resistance, and the nitC and cynS genes required for the assimilation of cyanide or cyanate, respectively. While cyanide resistance genes presented a high frequency among the analyzed genomes, genes responsible for cyanide or cyanate assimilation were identified in a considerably lower proportion. According to the results obtained in this work, an in silico approach based on a comparative genomic approach can be considered as an agile strategy for the bioprospection of putative cyanotrophic bacteria and for the identification of new genes putatively involved in cyanide biodegradation.

Accounting for 16S rRNA copy number prediction uncertainty and its implications in bacterial diversity analyses

16S rRNA gene copy number (16S GCN) varies among bacterial species and this variation introduces potential biases to microbial diversity analyses using 16S rRNA read counts. To correct the biases, methods have been developed to predict 16S GCN. A recent study suggests that the prediction uncertainty can be so great that copy number correction is not justified in practice. Here we develop RasperGade16S, a novel method and software to better model and capture the inherent uncertainty in 16S GCN prediction. RasperGade16S implements a maximum likelihood framework of pulsed evolution model and explicitly accounts for intraspecific GCN variation and heterogeneous GCN evolution rates among species. Using cross-validation, we show that our method provides robust confidence estimates for the GCN predictions and outperforms other methods in both precision and recall. We have predicted GCN for 592605 OTUs in the SILVA database and tested 113842 bacterial communities that represent an exhaustive and diverse list of engineered and natural environments. We found that the prediction uncertainty is small enough for 99% of the communities that 16S GCN correction should improve their compositional and functional profiles estimated using 16S rRNA reads. On the other hand, we found that GCN variation has limited impacts on beta-diversity analyses such as PCoA, NMDS, PERMANOVA and random-forest test.

Comparative genomics and informational content analysis uncovered internal regions of the core genes rpoD, pepN and gltX for an MLSA with genome-level resolving power within the genus Pseudomonas

In the field of prokaryotic taxonomy, there has been a recent transition towards phylogenomics as the gold standard approach. However, genome-based phylogenetics is still restrictive for its cost when managing large amounts of isolates. Fast, cheap, and taxonomically competent alternatives, like multilocus sequence analysis (MLSA) are thus recommendable. Nevertheless, the criteria for selecting the conserved genes for MLSA have not been explicit for different bacterial taxa, including the broadly diverse Pseudomonas genus. Here, we have carried out an unbiased and rational workflow to select internal sequence regions of Pseudomonas core genes (CG) for a MLSA with the best phylogenetic power, and with a resolution comparable to the genome-based ANI approach. A computational workflow was established to inspect 126 complete genomes of representatives from over 60 Pseudomonas species and subspecies, in order to identify the most informative CG internal regions and determine which combinations in sets of three partial CG sequences have comparable phylogenetic resolution to that of the current ANI standard. We found that the rpoD^346-1196-pepN^1711-2571-gltX^86-909 concatenated sequences were the best performing in terms of phylogenetic robustness and resulted highly sensitive and specific when contrasted with ANI. The rpoD-pepN-gltX MLSA was validated in silico and in vitro. Altogether, the results presented here supports the proposal of the rpoD-pepN-gltX MLSA as a fast, affordable, and robust phylogenetic tool for members of the Pseudomonas genus.

Pseudomonas lalucatii sp. nov. isolated from Vallgornera, a karstic cave in Mallorca, Western Mediterranean

Caves are extreme underground environments colonized by oligotrophic bacterial communities that influence mineral transformations. The identification at the species level is important and this study aims to the taxonomic characterisation of four bacterial strains previously isolated from rock surfaces and water samples from a karstic cave located on Mallorca (Spain) that were assigned to the genus Pseudomonas according to 16S rRNA nucleotide sequence analysis. Sequence analysis of the RNA polymerase sigma factor gene (rpoD) allocated these strains to the P. fluorescens lineage within the P. anguilliseptica phylogenetic group, close to the P. benzenivorans type strain. A polyphasic taxonomic approach included phenotypic characterization, fatty acid composition analysis, and whole-cell protein profiling, together with phylogenomic data. The results supported the proposal of a new species in the Pseudomonas genus. Characteristic fatty acid methyl esters of members of the Pseudomonas genus were present (C16:0, C10:0 3-OH, C12:0 2-OH and C12:0 3-OH) and the C12:1 3OH content differentiated these strains from P. benzenivorans. The genomic G + C mol% content of the four sequenced genomes was 66.9%. The average nucleotide indices based on BLAST analysis and the calculation of genome-to-genome distance with respect to their closest relative were lower than 88% and 30%, respectively. These data confirm that the four isolates, R1b-4, R1b-52A, A2bC-1 and R1b-54^T, represent a new species, for which the name Pseudomonas lalucatii is proposed, with strain R1b-54^T as the type strain (=CCUG 74754^T = CECT 30179^T). This is the first species in the P. anguilliseptica group isolated from this extreme habitat.

Division of the genus Chryseobacterium: Observation of discontinuities in amino acid identity values, a possible consequence of major extinction events, guides transfer of nine species to the genus Epilithonimonas, eleven species to the genus Kaistella, and three species to the genus Halpernia gen. nov., with description of Kaistella daneshvariae sp. nov. and Epilithonimonas vandammei sp. nov. derived from clinical specimens

The genus Chryseobacterium in the family Weeksellaceae is known to be polyphyletic. Amino acid identity (AAI) values were calculated from whole-genome sequences of species of the genus Chryseobacterium, and their distribution was found to be multi-modal. These naturally-occurring non-continuities were leveraged to standardise genus assignment of these species. We speculate that this multi-modal distribution is a consequence of loss of biodiversity during major extinction events, leading to the concept that a bacterial genus corresponds to a set of species that diversified since the Permian extinction. Transfer of nine species (Chryseobacterium arachidiradicis, Chryseobacterium bovis, Chryseobacterium caeni, Chryseobacterium hispanicum, Chryseobacterium hominis, Chryseobacterium hungaricum,, Chryseobacterium pallidum and Chryseobacterium zeae) to the genus Epilithonimonas and eleven (Chryseobacterium anthropi, Chryseobacterium antarcticum, Chryseobacterium carnis, Chryseobacterium chaponense, Chryseobacterium haifense, Chryseobacterium jeonii, Chryseobacterium montanum, Chryseobacterium palustre, Chryseobacterium solincola, Chryseobacterium treverense and Chryseobacterium yonginense) to the genus Kaistella is proposed. Two novel species are described: Kaistella daneshvariae sp. nov. and Epilithonimonas vandammei sp. nov. Evidence is presented to support the assignment of Planobacterium taklimakanense to a genus apart from Chryseobacterium, to which Planobacterium salipaludis comb nov. also belongs. The novel genus Halpernia is proposed, to contain the type species Halpernia frigidisoli comb. nov., along with Halpernia humi comb. nov., and Halpernia marina comb. nov.

Reliable Identification of Environmental Pseudomonas Isolates Using the rpoD Gene

The taxonomic affiliation of Pseudomonas isolates is currently assessed by using the 16S rRNA gene, MultiLocus Sequence Analysis (MLSA), or whole genome sequencing. Therefore, microbiologists are facing an arduous choice, either using the universal marker, knowing that these affiliations could be inaccurate, or engaging in more laborious and costly approaches. The rpoD gene, like the 16S rRNA gene, is included in most MLSA procedures and has already been suggested for the rapid identification of certain groups of Pseudomonas. However, a comprehensive overview of the rpoD-based phylogenetic relationships within the Pseudomonas genus is lacking. In this study, we present the rpoD-based phylogeny of 217 type strains of Pseudomonas and defined a cutoff value of 98% nucleotide identity to differentiate strains at the species level. To validate this approach, we sequenced the rpoD of 145 environmental isolates and complemented this analysis with whole genome sequencing. The rpoD sequence allowed us to accurately assign Pseudomonas isolates to 20 known species and represents an excellent first diagnostic tool to identify new Pseudomonas species. Finally, rpoD amplicon sequencing appears as a reliable and low-cost alternative, particularly in the case of large environmental studies with hundreds or thousands of isolates.

Pseudomonas poae-Associated Fatal Septic Transfusion Reaction, Peoria, Illinois, USA, 2017

In the United States, fatal transfusion-transmitted infections from red blood cell units are rare. Although this pattern mostly reflects how inhospitable refrigerated red blood cell units are to contaminant growth, fatalities caused by microorganisms that can grow at storage temperature (4°C), but not in standard clinical blood cultures at 37°C, are probably underestimated. We analyzed a fatal red blood cell transfusion in Peoria, Illinois, USA, that occurred in 2017. Samples from the patient's whole blood and the red blood cell unit remained culture-negative during the investigation, despite direct visualization of gram-negative bacilli within the unit immediately after transfusion. We identified the bacteria as Pseudomonas poae, a nonpathogenic pseudomonad carrying multiple cold-shock domain protein genes, and confirmed its cold tolerance and inability to grow at 37°C. Our work indicates transfusion reaction workups need to include testing for psychrophilic organisms, which could explain the cause of other apparently culture-negative transfusion reactions.

Species Interaction and Selective Carbon Addition During Antibiotic Exposure Enhances Bacterial Survival

Biofilms are multifaceted and robust microbiological systems that enable microorganisms to withstand a multitude of environmental stresses and expand their habitat range. We have shown previously that nutritional status alters antibiotic susceptibility in a mixed-species biofilm. To further elucidate the effects of nutrient addition on inter-species dynamics and whole-biofilm susceptibility to high-dose streptomycin exposures, a CO₂ Evolution Measurement System was used to monitor the metabolic activity of early steady state pure-culture and mixed-species biofilms containing Pseudomonas aeruginosa and Stenotrophomonas maltophilia, with and without added carbon. Carbon supplementation was needed for biofilm recovery from high-dose streptomycin exposures when P. aeruginosa was either the dominant community member in a mixed-species biofilm (containing predominantly P. aeruginosa and S. maltophilia) or as a pure culture. By contrast, S. maltophilia biofilms could recover from high-dose streptomycin exposures without the need for carbon addition during antibiotic exposure. Metagenomic analysis revealed that even when inocula were dominated by Pseudomonas, the relative abundance of Stenotrophomonas increased upon biofilm development to ultimately become the dominant species post-streptomycin exposure. The combined metabolic and metagenomic results demonstrated the relevance of inter-species influence on survival and that nutritional status has a strong influence on the survival of P. aeruginosa dominated biofilms.

Transcriptional regulation of fatty acid cis-trans isomerization in the solvent-tolerant soil bacterium, Pseudomonas putida F1

One key to the success of Pseudomonas spp. is their ability to reside in hostile environments. Pseudomonas spp. possess a cis-trans isomerase (Cti) an enzyme that converts the cis-unsaturated fatty acids (FAs) of the membrane lipids to their trans-isomers to rigidify the membrane and thereby resist stresses. Whereas the posttranslational Cti regulation has been previously reported, transcriptional cti regulation remains to be studied in more details. Here, we have studied cti transcriptional regulation in the solvent-tolerant strain Pseudomonas putida F1. Two cti transcriptional start sites (cti-279 and cti-77) were identified with cti-279 transcript being dominant. Expression of cti was found to increase with temperature increase, addition of the organic solvent, octanol and in the stationary growth phase. We found that cti expression was repressed by the cyclic-AMP receptor protein (Crp) and repression required the cyclic-AMP ligand of Crp. Production of trans-unsaturated FAs was found to decrease after 24 h of growth. Although this decrease was accompanied by an increase in cyclopropane FA content, this was not at the expense of trans-unsaturated FAs demonstrating the absence of competition between Cti and Cfa in FA modification.

Disentangling the Influence of Environment, Host Specificity and Thallus Differentiation on Bacterial Communities in Siphonous Green Seaweeds

Siphonous green seaweeds, such as Caulerpa, are among the most morphologically complex algae with differentiated algal structures (morphological niches). Caulerpa is also host to a rich diversity of bacterial endo- and epibionts. The degree to which these bacterial communities are species-, or even niche-specific remains largely unknown. To address this, we investigated the diversity of bacteria associated to different morphological niches of both native and invasive species of Caulerpa from different geographic locations along the Turkish coastline of the Aegean sea. Associated bacteria were identified using the 16S rDNA marker gene for three morphological niches, such as the endobiome, epibiome, and rhizobiome. Bacterial community structure was explored and deterministic factors behind bacterial variation were investigated. Of the total variation, only 21.5% could be explained. Pronounced differences in bacterial community composition were observed and variation was partly explained by a combination of host species, biogeography and nutrient levels. The majority of the explained bacterial variation within the algal holobiont was attributed to the micro-environments established by distinct morphological niches. This study further supports the hypothesis that the bacterial assembly is largely stochastic in nature and bacterial community structure is most likely linked to functional genes rather than taxonomy.

Genome-Based Analysis Reveals the Taxonomy and Diversity of the Family Idiomarinaceae

Idiomarinaceae is a family of Gram-stain negative, mesophilic euryhalophiles. To provide a robust framework for the evolutionary and taxonomic relationships of bacteria of this family, we compared herein the genomes of 36 type strains and 43 non-type strains using 16S rRNA gene sequences, core genome based 78 single-copy orthologous proteins, digital DNA-DNA hybridization and average nucleotide identity (ANI) estimation. The 79 bacteria of this family were consistently divided into taxon I, taxon II, and taxon III corresponding to the three genera Idiomarina, Pseudidiomarina, and Aliidiomarina, which contained 13 putative new genospecies in addition to 35 well-defined species represented by each type strain. Furthermore, genetic diversity of this family was evident at the genus- and species levels, and exceeded that which is defined currently by the named species. In view of multiple genotypic characteristics clearly distinct from the other two genera, we propose reinstating the genus Pseudidiomarina as a monophyletic taxon. Taken together, this is the first genome-based study of the taxonomy and diversity of bacteria within the family Idiomarinaceae, and will contribute to further insights into microbial evolution and adaptation to saline environments.

Exploring the microbiome of the "star" freshwater diatom Asterionella formosa in a laboratory context

Most of our knowledge on the mechanisms underlying diatom-bacterial interactions has been acquired through studies involving isolation of culturable partners. Here, we established a laboratory model of intermediate complexity between complex natural communities and laboratory pure culture models. We investigated the whole community formed by the freshwater diatom Asterionella formosa and its associated bacteria in a laboratory context, including both culturable and unculturable bacteria. Combining cellular and molecular approaches, we showed that in laboratory cultures, A. formosa microbiome was dynamic and comprised of numerous bacterial species (mainly Proteobacteria and Bacteroidetes). Using metagenomics, we explored several metabolic potentials present within the bacterial community. Our analyses suggested that bacteria were heterotrophic although a third of them (Alpha- and Beta-proteobacteria) could also be phototrophic. About 60% of the bacteria, phylogenetically diverse, could metabolize glycolate. The capacity to synthesize molecules such as B vitamins appeared unevenly distributed among bacteria. Altogether, our results brought insights into the bacterial diversity found in diatom-bacterial communities and hinted at metabolic interdependencies within the community that could result in diatom-bacterial and bacterial-bacterial interactions. The present work allowed us to explore the functional architecture of the bacterial community associated with A. formosa in culture and is complementary to field studies.

Pseudomonas-Specific NGS Assay Provides Insight Into Abundance and Dynamics of Pseudomonas Species Including P. aeruginosa in a Cooling Tower

Pseudomonas species are frequent inhabitants of freshwater environments and colonizers of water supply networks via bioadhesion and biofilm formation. P. aeruginosa is the species most commonly associated with human disease, causing a wide variety of infections with links to its presence in freshwater systems. Though several other Pseudomonas species are of ecological and public health importance, little knowledge exists regarding environmental abundances of these species. In the present study, an Illumina-based next-generation sequencing (NGS) approach using Pseudomonas-specific primers targeting the 16S rRNA gene was evaluated and applied to a set of freshwater samples from different environments including a cooling tower sampled monthly during 2 years. Our approach showed high in situ specificity and accuracy. NGS read counts revealed a precise quantification of P. aeruginosa and a good correlation with the absolute number of Pseudomonas genome copies in a validated genus-specific qPCR assay, demonstrating the ability of the NGS approach to determine both relative and absolute abundances of Pseudomonas species and P. aeruginosa. The characterization of Pseudomonas communities in cooling tower water allowed us to identify 43 phylotypes, with P. aeruginosa being the most abundant. A shift existed within each year from a community dominated by phylotypes belonging to P. fluorescens and P. oleovorans phylogenetic groups to a community where P. aeruginosa was highly abundant. Co-occurrence was observed between P. aeruginosa and other phylotypes of P. aeruginosa group as well as the potentially pathogenic species P. stutzeri, but not with phylotypes of the P. fluorescens group, indicating the need to further investigate the metabolic networks and ecological traits of Pseudomonas species. This study demonstrates the potential of deep sequencing as a valuable tool in environmental diagnostics and surveillance of health-related pathogens in freshwater environments.

Multiple Ribosomal RNA Operons in Bacteria; Their Concerted Evolution and Potential Consequences on the Rate of Evolution of Their 16S rRNA

Bacterial species differ greatly in the number and location of the rRNA operons which may be present in the bacterial chromosomes and plasmids. Most bacterial species contain more than one ribosomal RNA operon copy in their genomes, with some species containing up to 15 such copies. We review the number and location of the rRNA operons and discuss evolution of 16S rRNA (rrs) genes -which are considered as ultimate chronometers for phylogenetic classification- in bacteria with multiple copies of these genes. In these bacterial species, the rrs genes must evolve in concert and sequence changes generated by mutation or horizontal gene transfer must be either erased or spread to every gene copy to avoid divergence, as it occurs when they are present in different species. Analysis of polymorphic sites in intra-genomic rrs copies identifies putative conversion events and demonstrates that sequence conversion is patchy and occurs in small conversion tracts. Sequence conversion probably arises by a non-reciprocal transfer between two or more copies where one copy contributes only a small contiguous segment of DNA, whereas the other copy contributes the rest of the genome in a fairly well understood molecular process. Because concerted evolution implies that a mutation in any of the rrs copies is either eliminated or transferred to every rrs gene in the genome, this process should slow their evolution rate relative to that of single copy genes. However, available data on the rrs genes in bacterial genomes do not show a clear relationship between their evolution rates and the number of their copies in the genome.

Validation of RT-qPCR Approaches to Monitor Pseudomonas syringae Gene Expression During Infection and Exposure to Pattern-Triggered Immunity

Pseudomonas syringae pv. tomato DC3000 is an important model plant pathogen, with a fully annotated genome and multiple compatible plant hosts. Very few studies have examined the regulation of DC3000 gene expression in vivo. We developed a quantitative reverse transcription-polymerase chain reaction assay to monitor transcriptional changes in DC3000 inoculated into Arabidopsis thaliana leaves during disease and exposure to pattern-triggered immunity (PTI). In our approach, bacterial RNA concentrations in total tissue RNA are standardized using P. syringae-specific 16S ribosomal RNA primers. We validated multiple stable reference genes for normalization in calculating the relative expression of genes of interest. We used empirically derived rates of amplification efficiency to calculate relative expression of key marker genes for virulence-associated regulation. We demonstrated that exposure to PTI alters DC3000 expression of type III secretion system, coronatine synthesis genes, and flagellar marker genes.

Genomic insight into the taxonomy of Rhizobium genospecies that nodulate Phaseolus vulgaris

Due to the wide cultivation of bean (Phaseolus vulgaris L.), rhizobia associated with this plant have been isolated from many different geographical regions. In order to investigate the species diversity of bean rhizobia, comparative genome sequence analysis was performed in the present study for 69 Rhizobium strains mainly isolated from root nodules of bean and clover (Trifolium spp.). Based on genome average nucleotide identity, digital DNA:DNA hybridization, and phylogenetic analysis of 1,458 single-copy core genes, these strains were classified into 28 clusters, consistent with their species definition based on multilocus sequence analysis (MLSA) of atpD, glnII, and recA. The bean rhizobia were found in 16 defined species and nine putative novel species; in addition, 35 strains previously described as Rhizobium etli, Rhizobium phaseoli, Rhizobium vallis, Rhizobium gallicum, Rhizobium leguminosarum and Rhizobium spp. should be renamed. The phylogenetic patterns of symbiotic genes nodC and nifH were highly host-specific and inconsistent with the genomic phylogeny. Multiple symbiovars (sv.) within the Rhizobium species were found as a common feature: sv. phaseoli, sv. trifolii and sv. viciae in Rhizobium anhuiense; sv. phaseoli and sv. mimosae in Rhizobium sophoriradicis/R. etli/Rhizobium sp. III; sv. phaseoli and sv. trifolii in Rhizobium hidalgonense/Rhizobium acidisoli; sv. phaseoli and sv. viciae in R. leguminosarum/Rhizobium sp. IX; sv. trifolii and sv. viciae in Rhizobium laguerreae. Thus, genomic comparison revealed great species diversity in bean rhizobia, corrected the species definition of some previously misnamed strains, and demonstrated the MLSA a valuable and simple method for defining Rhizobium species.

Highly divergent 16S rRNA sequences in ribosomal operons of Scytonema hyalinum (Cyanobacteria)

A highly divergent 16S rRNA gene was found in one of the five ribosomal operons present in a species complex currently circumscribed as Scytonema hyalinum (Nostocales, Cyanobacteria) using clone libraries. If 16S rRNA sequence macroheterogeneity among ribosomal operons due to insertions, deletions or truncation is excluded, the sequence heterogeneity observed in S. hyalinum was the highest observed in any prokaryotic species thus far (7.3-9.0%). The secondary structure of the 16S rRNA molecules encoded by the two divergent operons was nearly identical, indicating possible functionality. The 23S rRNA gene was examined for a few strains in this complex, and it was also found to be highly divergent from the gene in Type 2 operons (8.7%), and likewise had nearly identical secondary structure between the Type 1 and Type 2 operons. Furthermore, the 16S-23S ITS showed marked differences consistent between operons among numerous strains. Both operons have promoter sequences that satisfy consensus requirements for functional prokaryotic transcription initiation. Horizontal gene transfer from another unknown heterocytous cyanobacterium is considered the most likely explanation for the origin of this molecule, but does not explain the ultimate origin of this sequence, which is very divergent from all 16S rRNA sequences found thus far in cyanobacteria. The divergent sequence is highly conserved among numerous strains of S. hyalinum, suggesting adaptive advantage and selective constraint of the divergent sequence.

Removal of Exogenous Materials from the Outer Portion of Frozen Cores to Investigate the Ancient Biological Communities Harbored Inside

The cryosphere offers access to preserved organisms that persisted under past environmental conditions. In fact, these frozen materials could reflect conditions over vast time periods and investigation of biological materials harbored inside could provide insight of ancient environments. To appropriately analyze these ecosystems and extract meaningful biological information from frozen soils and ice, proper collection and processing of the frozen samples is necessary. This is especially critical for microbial and DNA analyses since the communities present may be so uniquely different from modern ones. Here, a protocol is presented to successfully collect and decontaminate frozen cores. Both the absence of the colonies used to dope the outer surface and exogenous DNA suggest that we successfully decontaminated the frozen cores and that the microorganisms detected were from the material, rather than contamination from drilling or processing the cores.

Microbiome of HIV-infected people

Consistent interactions between the gut microbiome and adaptive immunity recently led several research groups to evaluate modifications of human gut microbiota composition during HIV infection. Herein we propose to review the shifts reported in infected individuals, as their correlation to disease progression. Though the gut microbiota is consistently altered in HIV individuals, the literature reveals several discrepancies, such as changes in microbial diversity associated with HIV status, taxa modified in infected subjects or influence of ART on gut flora restoration. Similarly, mechanisms involved in interactions between gut bacteria and immunity are to date poorly elucidated, emphasizing the importance of understanding how microbes can promote HIV replication. Further research is needed to propose adjuvant therapeutics dedicated to controlling disease progression through gut microbiome restoration.

The Road to Metagenomics: From Microbiology to DNA Sequencing Technologies and Bioinformatics

The study of microorganisms that pervade each and every part of this planet has encountered many challenges through time such as the discovery of unknown organisms and the understanding of how they interact with their environment. The aim of this review is to take the reader along the timeline and major milestones that led us to modern metagenomics. This new and thriving area is likely to be an important contributor to solve different problems. The transition from classical microbiology to modern metagenomics studies has required the development of new branches of knowledge and specialization. Here, we will review how the availability of high-throughput sequencing technologies has transformed microbiology and bioinformatics and how to tackle the inherent computational challenges that arise from the DNA sequencing revolution. New computational methods are constantly developed to collect, process, and extract useful biological information from a variety of samples and complex datasets, but metagenomics needs the integration of several of these computational methods. Despite the level of specialization needed in bioinformatics, it is important that life-scientists have a good understanding of it for a correct experimental design, which allows them to reveal the information in a metagenome.

Isolation of a diphenylamine-degrading bacterium and characterization of its metabolic capacities, bioremediation and bioaugmentation potential

The antioxidant diphenylamine (DPA) is used in fruit-packaging plants for the control of the physiological disorder apple scald. Its use results in the production of DPA-contaminated wastewater which should be treated before finally discharged. Biological treatment systems using tailored-made microbial inocula with specific catabolic activities comprise an appealing and sustainable solution. This study aimed to isolate DPA-degrading bacteria, identify the metabolic pathway of DPA and evaluate their potential for future implementation in bioremediation and biodepuration applications. A Pseudomonas putida strain named DPA1 able to rapidly degrade and utilize DPA as the sole C and N source was enriched from a DPA-contaminated soil. The isolated strain degraded spillage-level concentrations of DPA in liquid culture (2000 mg L(-1)) and in contaminated soil (1000 mg kg(-1)) and metabolized DPA via the transient formation of aniline and catechol. Further evidence for the bioremediation and biodepuration potential of the P. putida strain DPA1 was provided by its capacity to degrade the post-harvest fungicide ortho-phenylphenol (OPP), concurrently used by the fruit-packaging plants, although at slower rates and DPA in a wide range of pH (4.5-9) and temperatures (15-37 °C). These findings revealed the high potential of the P. putida strain DPA1 for use in future soil bioremediation strategies and/or as start-up inocula in wastewater biodepuration systems.

Genomic insights into the taxonomic status of the Bacillus cereus group

The identification and phylogenetic relationships of bacteria within the Bacillus cereus group are controversial. This study aimed at determining the taxonomic affiliations of these strains using the whole-genome sequence-based Genome BLAST Distance Phylogeny (GBDP) approach. The GBDP analysis clearly separated 224 strains into 30 clusters, representing eleven known, partially merged species and accordingly 19–20 putative novel species. Additionally, 16S rRNA gene analysis, a novel variant of multi-locus sequence analysis (nMLSA) and screening of virulence genes were performed. The 16S rRNA gene sequence was not sufficient to differentiate the bacteria within this group due to its high conservation. The nMLSA results were consistent with GBDP. Moreover, a fast typing method was proposed using the pycA gene, and where necessary, the ccpA gene. The pXO plasmids and cry genes were widely distributed, suggesting little correlation with the phylogenetic positions of the host bacteria. This might explain why classifications based on virulence characteristics proved unsatisfactory in the past. In summary, this is the first large-scale and systematic study of the taxonomic status of the bacteria within the B. cereus group using whole-genome sequences, and is likely to contribute to further insights into their pathogenicity, phylogeny and adaptation to diverse environments.

Fecal microbiota analysis: an overview of sample collection methods and sequencing strategies

Despite huge interest, there are still no universally accepted standards to conduct clinical studies in the field of gut microbiota analysis. Stool material is frequently used as a proxy of gut microbiota, but many different protocols can be used for collection and DNA extraction. Whereas 16S rRNA encoding gene amplification and sequencing has been widely used to study the composition of bacterial populations, it is now being challenged by the random, shotgun approach that brings far more information, although at a higher cost. In this review we give an overview of existing methods and important points to consider when conducting gut microbiota studies, with the objective to provide recommendations to those who would like to conduct such research.

Extending the conserved phylogenetic core of archaea disentangles the evolution of the third domain of life

Initial studies of the archaeal phylogeny relied mainly on the analysis of the RNA component of the small subunit of the ribosome (SSU rRNA). The resulting phylogenies have provided interesting but partial information on the evolutionary history of the third domain of life because SSU rRNA sequences do not contain enough phylogenetic signal to resolve all nodes of the archaeal tree. Thus, many relationships, and especially the most ancient ones, remained elusive. Moreover, SSU rRNA phylogenies can be heavily biased by tree reconstruction artifacts. The sequencing of complete genomes allows using a variety of protein markers as an alternative to SSU rRNA. Taking advantage of the recent burst of archaeal complete genome sequences, we have carried out an in-depth phylogenomic analysis of this domain. We have identified 200 new protein families that, in addition to the ribosomal proteins and the subunits of the RNA polymerase, form a conserved phylogenetic core of archaeal genes. The accurate analysis of these markers combined with desaturation approaches shed new light on the evolutionary history of Archaea and reveals that several relationships recovered in recent analyses are likely the consequence of tree reconstruction artifacts. Among others, we resolve a number of important relationships, such as those among methanogens Class I, and we propose the definition of two new superclasses within the Euryarchaeota: Methanomada and Diaforarchaea.

Flow cytometry and real-time quantitative PCR as tools for assessing plasmid persistence

The maintenance of a plasmid in the absence of selection for plasmid-borne genes is not guaranteed. However, plasmid persistence can evolve under selective conditions. Studying the molecular mechanisms behind the evolution of plasmid persistence is key to understanding how plasmids are maintained under nonselective conditions. Given the current crisis of rapid antibiotic resistance spread by multidrug resistance plasmids, this insight is of high medical relevance. The conventional method for monitoring plasmid persistence (i.e., the fraction of plasmid-containing cells in a population over time) is based on cultivation and involves differentiating colonies of plasmid-containing and plasmid-free cells on agar plates. However, this technique is time-consuming and does not easily lend itself to high-throughput applications. Here, we present flow cytometry (FCM) and real-time quantitative PCR (qPCR) as alternative tools for monitoring plasmid persistence. For this, we measured the persistence of a model plasmid, pB10::gfp, in three Pseudomonas hosts and in known mixtures of plasmid-containing and -free cells. We also compared three performance criteria: dynamic range, resolution, and variance. Although not without exceptions, both techniques generated estimates of overall plasmid loss rates that were rather similar to those generated by the conventional plate count (PC) method. They also were able to resolve differences in loss rates between artificial plasmid persistence assays. Finally, we briefly discuss the advantages and disadvantages for each technique and conclude that, overall, both FCM and real-time qPCR are suitable alternatives to cultivation-based methods for routine measurement of plasmid persistence, thereby opening avenues for high-throughput analyses.

Core functional traits of bacterial communities in the Upper Mississippi River show limited variation in response to land cover

Taxonomic characterization of environmental microbial communities via high-throughput DNA sequencing has revealed that patterns in microbial biogeography affect community structure. However, shifts in functional diversity related to variation in taxonomic composition are poorly understood. To overcome limitations due to the prohibitive cost of high-depth metagenomic sequencing, tools to infer functional diversity based on phylogenetic distributions of functional traits have been developed. In this study we characterized functional microbial diversity at 11 sites along the Mississippi River in Minnesota using both metagenomic sequencing and functional-inference-based (PICRUSt) approaches. This allowed us to determine how distance and variation in land cover throughout the river influenced the distribution of functional traits, as well as to validate PICRUSt inferences. The distribution and abundance of functional traits, by metagenomic analysis, were similar among sites, with a median standard deviation of 0.0002% among tier 3 functions in KEGG. Overall inferred functional variation was significantly different (P ≤ 0.035) between two water basins surrounded by agricultural vs. developed land cover, and abundances of bacterial orders that correlated with functional traits by metagenomic analysis were greater where abundances of the trait were inferred to be higher. PICRUSt inferences were significantly correlated (r = 0.147, P = 1.80 × 10(-30)) with metagenomic annotations. Discrepancies between metagenomic and PICRUSt taxonomic-functional relationships, however, suggested potential functional redundancy among abundant and rare taxa that impeded the ability to accurately assess unique functional traits among rare taxa at this sequencing depth. Results of this study suggest that a suite of "core functional traits" is conserved throughout the river and distributions of functional traits, rather than specific taxa, may shift in response to environmental heterogeneity.

Biodegradation of chlorobenzene, 1,2-dichlorobenzene, and 1,4-dichlorobenzene in the vadose zone

Much of the microbial activity in nature takes place at interfaces, which are often associated with redox discontinuities. One example is the oxic/anoxic interface where polluted groundwater interacts with the overlying vadose zone. We tested whether microbes in the vadose zone can use synthetic chemicals as electron donors and thus protect the overlying air and buildings from groundwater pollutants. Samples from the vadose zone of a site contaminated with chlorobenzene (CB), 1,2-dichlorobenzene (12DCB), and 1,4-dichlorobenzene (14DCB) were packed in a multiport column to simulate the interface of the vadose zone with an underlying groundwater plume. A mixture of CB, 12DCB, and 14DCB in anoxic water was pumped continuously through the bottom of column to an outlet below the first sampling port to create an oxic/anoxic interface and a capillary fringe. Removal to below the detection limits by rapid biodegradation with rates of 21 ± 1 mg of CB • m(-2) • d(-1), 3.7 ± 0.5 mg of 12DCB • m(-2) • d(-1), and 7.4 ± 0.7 mg of 1.4 DCB • m(-2) • d(-1) indicated that natural attenuation in the capillary fringe can prevent the migration of CB, 12DCB, and 14DCB vapors. Enumeration of bacteria capable of degrading chlorobenzenes suggested that most of the biodegradation takes place within the first 10 cm above the saturated zone. Biodegradation also increased the upward flux of contaminants and thus enhanced their elimination from the underlying water. The results revealed a substantial biodegradation capacity for chlorinated aromatic compounds at the oxic/anoxic interface and illustrate the role of microbes in creating steep redox gradients.

Intra- and intergenomic variation of ribosomal RNA operons in concurrent Alteromonas macleodii strains

Biodiversity estimates based on ribosomal operon sequence diversity rely on the premise that a sequence is characteristic of a single specific taxon or operational taxonomic unit (OTU). Here, we have studied the sequence diversity of 14 ribosomal RNA operons (rrn) contained in the genomes of two isolates (five operons in each genome) and four metagenomic fosmids, all from the same seawater sample. Complete sequencing of the isolate genomes and the fosmids establish that they represent strains of the same species, Alteromonas macleodii, with average nucleotide identity (ANI) values >97 %. Nonetheless, we observed high levels of intragenomic heterogeneity (i.e., variability between operons of a single genome) affecting multiple regions of the 16S and 23S rRNA genes as well as the internally transcribed spacer 1 (ITS-1) region. Furthermore, the ribosomal operons exhibited intergenomic heterogeneity (i.e., variability between operons located in separate genomes) in each of these regions, compounding the variability. Our data reveal the extensive heterogeneity observed in natural populations of A. macleodii at a single point in time and support the idea that distinct lineages of A. macleodii exist in the deep Mediterranean. These findings highlight the potential of rRNA fingerprinting methods to misrepresent species diversity while simultaneously failing to recognize the ecological significance of individual strains.

Legionellaand non-Legionellabacteria in a biological treatment plant

Legionella pneumophila were previously identified in the aeration ponds (up to 1010CFU/L) of a biological wastewater treatment plant at Borregaard Ind. Ltd., Sarpsborg, Norway, and in air samples (up to 3300 CFU/m3) collected above the aeration ponds. After 3 outbreaks of Legionnaires’ disease reported in this area in 2005 and 2008, the aeration ponds of the plant were shut down by the Norwegian authorities in September 2008. The aim of the present work was to analyze the Legionella and non-Legionella bacterial communities in the aeration ponds before and during the shutdown process and to identify potential human pathogens. The non-Legionella bacterial community was investigated in selected samples during the shutdown process by 16S rDNA sequencing of clone libraries (400 clones) and growth analysis. The concentration of L. pneumophila and Pseudomonas spp. DNA were monitored by quantitative PCR. Results showed a decrease in the concentration of L. pneumophila and Pseudomonas spp. during the shutdown. This was accompanied by a significant change in the composition of the bacterial community in the aeration ponds. This study demonstrated that several advanced analytical methods are necessary to characterize the bacterial population in complex environments, such as the industrial aeration ponds.