Konstantinos T. Konstantinidis

The bacterial species definition in the genomic era

The bacterial species definition, despite its eminent practical significance for identification, diagnosis, quarantine and diversity surveys, remains a very difficult issue to advance. Genomics now offers novel insights into intra-species diversity and the potential for emergence of a more soundly based system. Although we share the excitement, we argue that it is premature for a universal change to the definition because current knowledge is based on too few phylogenetic groups and too few samples of natural populations. Our analysis of five important bacterial groups suggests, however, that more stringent standards for species may be justifiable when a solid understanding of gene content and ecological distinctiveness becomes available. Our analysis also reveals what is actually encompassed in a species according to the current standards, in terms of whole-genome sequence and gene-content diversity, and shows that this does not correspond to coherent clusters for the environmental Burkholderia and Shewanella genera examined. In contrast, the obligatory pathogens, which have a very restricted ecological niche, do exhibit clusters. Therefore, the idea of biologically meaningful clusters of diversity that applies to most eukaryotes may not be universally applicable in the microbial world, or if such clusters exist, they may be found at different levels of distinction.

Genomic analysis of the uncultivated marine crenarchaeote Cenarchaeum symbiosum

Crenarchaeota are ubiquitous and abundant microbial constituents of soils, sediments, lakes, and ocean waters. To further describe the cosmopolitan nonthermophilic Crenarchaeota, we analyzed the genome sequence of one representative, the uncultivated sponge symbiont Cenarchaeum symbiosum. C. symbiosum genotypes coinhabiting the same host partitioned into two dominant populations, corresponding to previously described a- and b-type ribosomal RNA variants. Although they were syntenic, overlapping a- and b-type ribotype genomes harbored significant variability. A single tiling path comprising the dominant a-type genotype was assembled and used to explore the genomic properties of C. symbiosum and its planktonic relatives. Of 2,066 ORFs, 55.6% matched genes with predicted function from previously sequenced genomes. The remaining genes partitioned between functional RNAs (2.4%) and hypotheticals (42%) with limited homology to known functional genes. The latter category included some genes likely involved in the archaeal–sponge symbiotic association. Conversely, 525 C. symbiosum ORFs were most highly similar to sequences from marine environmental genomic surveys, and they apparently represent orthologous genes from free-living planktonic Crenarchaeota. In total, the C. symbiosum genome was remarkably distinct from those of other known Archaea and shared many core metabolic features in common with its free-living planktonic relatives.

Towards a Genome-Based Taxonomy for Prokaryotes

The ranks higher than the species in the prokaryotic taxonomy are primarily designated based on phylogenetic analysis of the 16S rRNA gene sequences, but no definite standards exist for the absolute relatedness (measured by 16S rRNA or other means) between the ranks. Accordingly, it remains unknown how comparable the ranks are between different organisms. To gain insights into this question, we studied the relationship between shared gene content and genetic relatedness for 175 fully sequenced strains, using as a robust measure of relatedness the average amino acid identity (AAI) of the shared genes. Our results reveal that adjacent ranks (e.g., phylum versus class) frequently show extensive overlap in terms of genetic and gene content relatedness of the grouped organisms, and hence, the current system is of limited predictive power in this respect. The overlap between nonadjacent ranks (e.g., phylum versus family) is generally limited and attributable to clear inconsistencies of the taxonomy. In addition to providing means for standardizing taxonomy, our AAI-based approach provides a means to evaluate the robustness of alternative genetic markers for phylogenetic purposes. For instance, the 23S rRNA gene was found to be as good a marker as the 16S rRNA gene, while several of the widely distributed protein-coding genes, such as the RNA polymerase and gyrase subunits, show a strong phylogenetic signal, albeit less strong than the rRNA genes (0.78 > R2 > 0.69 for the protein-coding genes versus R2 = 0.84 for the rRNA genes). The AAI approach outlined here could contribute significantly to a genome-based taxonomy for all microbial organisms.

Direct comparisons of Illumina vs. Roche 454 sequencing technologies on the same microbial community DNA sample.

Next-generation sequencing (NGS) is commonly used in metagenomic studies of complex microbial communities but whether or not different NGS platforms recover the same diversity from a sample and their assembled sequences are of comparable quality remain unclear. We compared the two most frequently used platforms, the Roche 454 FLX Titanium and the Illumina Genome Analyzer (GA) II, on the same DNA sample obtained from a complex freshwater planktonic community. Despite the substantial differences in read length and sequencing protocols, the platforms provided a comparable view of the community sampled. For instance, derived assemblies overlapped in ∼90% of their total sequences and in situ abundances of genes and genotypes (estimated based on sequence coverage) correlated highly between the two platforms (R2>0.9). Evaluation of base-call error, frameshift frequency, and contig length suggested that Illumina offered equivalent, if not better, assemblies than Roche 454. The results from metagenomic samples were further validated against DNA samples of eighteen isolate genomes, which showed a range of genome sizes and G+C% content. We also provide quantitative estimates of the errors in gene and contig sequences assembled from datasets characterized by different levels of complexity and G+C% content. For instance, we noted that homopolymer-associated, single-base errors affected ∼1% of the protein sequences recovered in Illumina contigs of 10× coverage and 50% G+C; this frequency increased to ∼3% when non-homopolymer errors were also considered. Collectively, our results should serve as a useful practical guide for choosing proper sampling strategies and data possessing protocols for future metagenomic studies.

Dehalococcoides mccartyi gen. nov., sp. nov., obligately organohalide-respiring anaerobic bacteria relevant to halogen cycling and bioremediation, belong to a novel bacterial class, Dehalococcoidia classis nov., order Dehalococcoidales ord. nov. and family Dehalococcoidaceae fam. nov., within the phylum Chloroflexi.

Six obligately anaerobic bacterial isolates (195(T), CBDB1, BAV1, VS, FL2 and GT) with strictly organohalide-respiring metabolisms were obtained from chlorinated solvent-contaminated aquifers, contaminated and uncontaminated river sediments or anoxic digester sludge. Cells were non-motile with a disc-shaped morphology, 0.3-1 µm in diameter and 0.1-0.2 µm thick, and characteristic indentations on opposite flat sides of the cell. Growth occurred in completely synthetic, reduced medium amended with a haloorganic electron acceptor (mostly chlorinated but also some brominated compounds), hydrogen as electron donor, acetate as carbon source, and vitamins. No other growth-supporting redox couples were identified. Aqueous hydrogen consumption threshold concentrations were <1 nM. Growth ceased when vitamin B(12) was omitted from the medium. Addition of sterile cell-free supernatant of Dehalococcoides-containing enrichment cultures enhanced dechlorination and growth of strains 195 and FL2, suggesting the existence of so-far unidentified stimulants. Dechlorination occurred between pH 6.5 and 8.0 and over a temperature range of 15-35 °C, with an optimum growth temperature between 25 and 30 °C. The major phospholipid fatty acids were 14 : 0 (15.7 mol%), br15 : 0 (6.2 mol%), 16 : 0 (22.7 mol%), 10-methyl 16 : 0 (25.8 mol%) and 18 : 0 (16.6 mol%). Unusual furan fatty acids including 9-(5-pentyl-2-furyl)-nonanoate and 8-(5-hexyl-2-furyl)-octanoate were detected in strains FL2, BAV1 and GT, but not in strains 195(T) and CBDB1. The 16S rRNA gene sequences of the six isolates shared more than 98 % identity, and phylogenetic analysis revealed an affiliation with the phylum Chloroflexi and more than 10 % sequence divergence from other described isolates. The genome sizes and G+C contents ranged from 1.34 to 1.47 Mbp and 47 to 48.9 mol% G+C, respectively. Based on 16S rRNA gene sequence comparisons, genome-wide average nucleotide identity and phenotypic characteristics, the organohalide-respiring isolates represent a new genus and species, for which the name Dehalococcoides mccartyi gen. nov., sp. nov. is proposed. Isolates BAV1 ( = ATCC BAA-2100  = JCM 16839  = KCTC 5957), FL2 ( = ATCC BAA-2098  = DSM 23585  = JCM 16840  = KCTC 5959), GT ( = ATCC BAA-2099  = JCM 16841  = KCTC 5958), CBDB1, 195(T) ( = ATCC BAA-2266(T)  = KCTC 15142(T)) and VS are considered strains of Dehalococcoides mccartyi, with strain 195(T) as the type strain. The new class Dehalococcoidia classis nov., order Dehalococcoidales ord. nov. and family Dehalococcoidaceae fam. nov. are described to accommodate the new taxon.

Burkholderia Xenovorans LB400 Harbors a Multi-Replicon, 9.73-Mbp Genome Shaped for Versatility

Burkholderia xenovorans LB400 (LB400), a well studied, effective polychlorinated biphenyl-degrader, has one of the two largest known bacterial genomes and is the first nonpathogenic Burkholderia isolate sequenced. From an evolutionary perspective, we find significant differences in functional specialization between the three replicons of LB400, as well as a more relaxed selective pressure for genes located on the two smaller vs. the largest replicon. High genomic plasticity, diversity, and specialization within the Burkholderia genus are exemplified by the conservation of only 44% of the genes between LB400 and Burkholderia cepacia complex strain 383. Even among four B. xenovorans strains, genome size varies from 7.4 to 9.73 Mbp. The latter is largely explained by our findings that >20% of the LB400 sequence was recently acquired by means of lateral gene transfer. Although a range of genetic factors associated with in vivo survival and intercellular interactions are present, these genetic factors are likely related to niche breadth rather than determinants of pathogenicity. The presence of at least eleven “central aromatic” and twenty “peripheral aromatic” pathways in LB400, among the highest in any sequenced bacterial genome, supports this hypothesis. Finally, in addition to the experimentally observed redundancy in benzoate degradation and formaldehyde oxidation pathways, the fact that 17.6% of proteins have a better LB400 paralog than an ortholog in a different genome highlights the importance of gene duplication and repeated acquirement, which, coupled with their divergence, raises questions regarding the role of paralogs and potential functional redundancies in large-genome microbes.

Unexpected nondenitrifier nitrous oxide reductase gene diversity and abundance in soils

Agricultural and industrial practices more than doubled the intrinsic rate of terrestrial N fixation over the past century with drastic consequences, including increased atmospheric nitrous oxide (N2O) concentrations. N2O is a potent greenhouse gas and contributor to ozone layer destruction, and its release from fixed N is almost entirely controlled by microbial activities. Mitigation of N2O emissions to the atmosphere has been attributed exclusively to denitrifiers possessing NosZ, the enzyme system catalyzing N2O to N2 reduction. We demonstrate that diverse microbial taxa possess divergent nos clusters with genes that are related yet evolutionarily distinct from the typical nos genes of denitirifers. nos clusters with atypical nosZ occur in Bacteria and Archaea that denitrify (44% of genomes), do not possess other denitrification genes (56%), or perform dissimilatory nitrate reduction to ammonium (DNRA; (31%). Experiments with the DNRA soil bacterium Anaeromyxobacter dehalogenans demonstrated that the atypical NosZ is an effective N2O reductase, and PCR-based surveys suggested that atypical nosZ are abundant in terrestrial environments. Bioinformatic analyses revealed that atypical nos clusters possess distinctive regulatory and functional components (e.g., Sec vs. Tat secretion pathway in typical nos), and that previous nosZ-targeted PCR primers do not capture the atypical nosZ diversity. Collectively, our results suggest that nondenitrifying populations with a broad range of metabolisms and habitats are potentially significant contributors to N2O consumption. Apparently, a large, previously unrecognized group of environmental nosZ has not been accounted for, and characterizing their contributions to N2O consumption will advance understanding of the ecological controls on N2O emissions and lead to refined greenhouse gas flux models.

Genome sequencing of environmental Escherichia coli expands understanding of the ecology and speciation of the model bacterial species

Defining bacterial species remains a challenging problem even for the model bacterium Escherichia coli and has major practical consequences for reliable diagnosis of infectious disease agents and regulations for transport and possession of organisms of economic importance. E. coli traditionally is thought to live within the gastrointestinal tract of humans and other warm-blooded animals and not to survive for extended periods outside its host; this understanding is the basis for its widespread use as a fecal contamination indicator. Here, we report the genome sequences of nine environmentally adapted strains that are phenotypically and taxonomically indistinguishable from typical E. coli (commensal or pathogenic). We find, however, that the commensal genomes encode for more functions that are important for fitness in the human gut, do not exchange genetic material with their environmental counterparts, and hence do not evolve according to the recently proposed fragmented speciation model. These findings are consistent with a more stringent and ecologic definition for bacterial species than the current definition and provide means to start replacing traditional approaches of defining distinctive phenotypes for new species with omics-based procedures. They also have important implications for reliable diagnosis and regulation of pathogenic E. coli and for the coliform cell-counting test.

Microbiome of the upper troposphere: Species composition and prevalence, effects of tropical storms, and atmospheric implications

The composition and prevalence of microorganisms in the middle-to-upper troposphere (8–15 km altitude) and their role in aerosol-cloud-precipitation interactions represent important, unresolved questions for biological and atmospheric science. In particular, airborne microorganisms above the oceans remain essentially uncharacterized, as most work to date is restricted to samples taken near the Earth’s surface. Here we report on the microbiome of low- and high-altitude air masses sampled onboard the National Aeronautics and Space Administration DC-8 platform during the 2010 Genesis and Rapid Intensification Processes campaign in the Caribbean Sea. The samples were collected in cloudy and cloud-free air masses before, during, and after two major tropical hurricanes, Earl and Karl. Quantitative PCR and microscopy revealed that viable bacterial cells represented on average around 20% of the total particles in the 0.25- to 1-μm diameter range and were at least an order of magnitude more abundant than fungal cells, suggesting that bacteria represent an important and underestimated fraction of micrometer-sized atmospheric aerosols. The samples from the two hurricanes were characterized by significantly different bacterial communities, revealing that hurricanes aerosolize a large amount of new cells. Nonetheless, 17 bacterial taxa, including taxa that are known to use C1–C4 carbon compounds present in the atmosphere, were found in all samples, indicating that these organisms possess traits that allow survival in the troposphere. The findings presented here suggest that the microbiome is a dynamic and underappreciated aspect of the upper troposphere with potentially important impacts on the hydrological cycle, clouds, and climate.

Microbial species delineation using whole genome sequences

Increased sequencing of microbial genomes has revealed that prevailing prokaryotic species assignments can be inconsistent with whole genome information for a significant number of species. The long-standing need for a systematic and scalable species assignment technique can be met by the genome-wide Average Nucleotide Identity (gANI) metric, which is widely acknowledged as a robust measure of genomic relatedness. In this work, we demonstrate that the combination of gANI and the alignment fraction (AF) between two genomes accurately reflects their genomic relatedness. We introduce an efficient implementation of AF,gANI and discuss its successful application to 86.5M genome pairs between 13,151 prokaryotic genomes assigned to 3032 species. Subsequently, by comparing the genome clusters obtained from complete linkage clustering of these pairs to existing taxonomy, we observed that nearly 18% of all prokaryotic species suffer from anomalies in species definition. Our results can be used to explore central questions such as whether microorganisms form a continuum of genetic diversity or distinct species represented by distinct genetic signatures. We propose that this precise and objective AF,gANI-based species definition: the MiSI (Microbial Species Identifier) method, be used to address previous inconsistencies in species classification and as the primary guide for new taxonomic species assignment, supplemented by the traditional polyphasic approach, as required.