Assessment of genome annotation using gene function similarity within the gene neighborhood.

BACKGROUND: Functional annotation of bacterial genomes is an obligatory and crucially important step of information processing from the genome sequences into cellular mechanisms. However, there is a lack of computational methods to evaluate the quality of functional assignments. RESULTS: We developed a genome-scale model that assigns Bayesian probability to each gene utilizing a known property of functional similarity between neighboring genes in bacteria. CONCLUSIONS: Our model clearly distinguished true annotation from random annotation with Bayesian annotation probability >0.95. Our model will provide a useful guide to quantitatively evaluate functional annotation methods and to detect gene sets with reliable annotations.

Diversity of Pseudomonas Genomes, Including Populus-Associated Isolates, as Revealed by Comparative Genome Analysis.

The Pseudomonas genus contains a metabolically versatile group of organisms that are known to occupy numerous ecological niches, including the rhizosphere and endosphere of many plants. Their diversity influences the phylogenetic diversity and heterogeneity of these communities. On the basis of average amino acid identity, comparative genome analysis of >1,000 Pseudomonas genomes, including 21 Pseudomonas strains isolated from the roots of native Populus deltoides (eastern cottonwood) trees resulted in consistent and robust genomic clusters with phylogenetic homogeneity. All Pseudomonas aeruginosa genomes clustered together, and these were clearly distinct from other Pseudomonas species groups on the basis of pangenome and core genome analyses. In contrast, the genomes of Pseudomonas fluorescens were organized into 20 distinct genomic clusters, representing enormous diversity and heterogeneity. Most of our 21 Populus-associated isolates formed three distinct subgroups within the major P. fluorescens group, supported by pathway profile analysis, while two isolates were more closely related to Pseudomonas chlororaphis and Pseudomonas putida. Genes specific to Populus-associated subgroups were identified. Genes specific to subgroup 1 include several sensory systems that act in two-component signal transduction, a TonB-dependent receptor, and a phosphorelay sensor. Genes specific to subgroup 2 contain hypothetical genes, and genes specific to subgroup 3 were annotated with hydrolase activity. This study justifies the need to sequence multiple isolates, especially from P. fluorescens, which displays the most genetic variation, in order to study functional capabilities from a pangenomic perspective. This information will prove useful when choosing Pseudomonas strains for use to promote growth and increase disease resistance in plants.

Insights from 20 years of bacterial genome sequencing.

Since the first two complete bacterial genome sequences were published in 1995, the science of bacteria has dramatically changed. Using third-generation DNA sequencing, it is possible to completely sequence a bacterial genome in a few hours and identify some types of methylation sites along the genome as well. Sequencing of bacterial genome sequences is now a standard procedure, and the information from tens of thousands of bacterial genomes has had a major impact on our views of the bacterial world. In this review, we explore a series of questions to highlight some insights that comparative genomics has produced. To date, there are genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. However, the distribution is quite skewed towards a few phyla that contain model organisms. But the breadth is continuing to improve, with projects dedicated to filling in less characterized taxonomic groups. The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system provides bacteria with immunity against viruses, which outnumber bacteria by tenfold. How fast can we go? Second-generation sequencing has produced a large number of draft genomes (close to 90 % of bacterial genomes in GenBank are currently not complete); third-generation sequencing can potentially produce a finished genome in a few hours, and at the same time provide methlylation sites along the entire chromosome. The diversity of bacterial communities is extensive as is evident from the genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. Genome sequencing can help in classifying an organism, and in the case where multiple genomes of the same species are available, it is possible to calculate the pan- and core genomes; comparison of more than 2000 Escherichia coli genomes finds an E. coli core genome of about 3100 gene families and a total of about 89,000 different gene families. Why do we care about bacterial genome sequencing? There are many practical applications, such as genome-scale metabolic modeling, biosurveillance, bioforensics, and infectious disease epidemiology. In the near future, high-throughput sequencing of patient metagenomic samples could revolutionize medicine in terms of speed and accuracy of finding pathogens and knowing how to treat them.

The evolution of host specialization in the vertebrate gut symbiont Lactobacillus reuteri.

Recent research has provided mechanistic insight into the important contributions of the gut microbiota to vertebrate biology, but questions remain about the evolutionary processes that have shaped this symbiosis. In the present study, we showed in experiments with gnotobiotic mice that the evolution of Lactobacillus reuteri with rodents resulted in the emergence of host specialization. To identify genomic events marking adaptations to the murine host, we compared the genome of the rodent isolate L. reuteri 100-23 with that of the human isolate L. reuteri F275, and we identified hundreds of genes that were specific to each strain. In order to differentiate true host-specific genome content from strain-level differences, comparative genome hybridizations were performed to query 57 L. reuteri strains originating from six different vertebrate hosts in combination with genome sequence comparisons of nine strains encompassing five phylogenetic lineages of the species. This approach revealed that rodent strains, although showing a high degree of genomic plasticity, possessed a specific genome inventory that was rare or absent in strains from other vertebrate hosts. The distinct genome content of L. reuteri lineages reflected the niche characteristics in the gastrointestinal tracts of their respective hosts, and inactivation of seven out of eight representative rodent-specific genes in L. reuteri 100-23 resulted in impaired ecological performance in the gut of mice. The comparative genomic analyses suggested fundamentally different trends of genome evolution in rodent and human L. reuteri populations, with the former possessing a large and adaptable pan-genome while the latter being subjected to a process of reductive evolution. In conclusion, this study provided experimental evidence and a molecular basis for the evolution of host specificity in a vertebrate gut symbiont, and it identified genomic events that have shaped this process.

Gene and translation initiation site prediction in metagenomic sequences.

MOTIVATION: Gene prediction in metagenomic sequences remains a difficult problem. Current sequencing technologies do not achieve sufficient coverage to assemble the individual genomes in a typical sample; consequently, sequencing runs produce a large number of short sequences whose exact origin is unknown. Since these sequences are usually smaller than the average length of a gene, algorithms must make predictions based on very little data. RESULTS: We present MetaProdigal, a metagenomic version of the gene prediction program Prodigal, that can identify genes in short, anonymous coding sequences with a high degree of accuracy. The novel value of the method consists of enhanced translation initiation site identification, ability to identify sequences that use alternate genetic codes and confidence values for each gene call. We compare the results of MetaProdigal with other methods and conclude with a discussion of future improvements. AVAILABILITY: The Prodigal software is freely available under the General Public License from http://code.google.com/p/prodigal/.

Enigmatic, ultrasmall, uncultivated Archaea.

Metagenomics has provided access to genomes of as yet uncultivated microorganisms in natural environments, yet there are gaps in our knowledge-particularly for Archaea-that occur at relatively low abundance and in extreme environments. Ultrasmall cells (<500 nm in diameter) from lineages without cultivated representatives that branch near the crenarchaeal/euryarchaeal divide have been detected in a variety of acidic ecosystems. We reconstructed composite, near-complete approximately 1-Mb genomes for three lineages, referred to as ARMAN (archaeal Richmond Mine acidophilic nanoorganisms), from environmental samples and a biofilm filtrate. Genes of two lineages are among the smallest yet described, enabling a 10% higher coding density than found genomes of the same size, and there are noncontiguous genes. No biological function could be inferred for up to 45% of genes and no more than 63% of the predicted proteins could be assigned to a revised set of archaeal clusters of orthologous groups. Some core metabolic genes are more common in Crenarchaeota than Euryarchaeota, up to 21% of genes have the highest sequence identity to bacterial genes, and 12 belong to clusters of orthologous groups that were previously exclusive to bacteria. A small subset of 3D cryo-electron tomographic reconstructions clearly show penetration of the ARMAN cell wall and cytoplasmic membranes by protuberances extended from cells of the archaeal order Thermoplasmatales. Interspecies interactions, the presence of a unique internal tubular organelle [Comolli, et al. (2009) ISME J 3:159-167], and many genes previously only affiliated with Crenarchaea or Bacteria indicate extensive unique physiology in organisms that branched close to the time that Cren- and Euryarchaeotal lineages diverged.

Complete genome sequence of Desulfurococcus fermentans, a hyperthermophilic cellulolytic crenarchaeon isolated from a freshwater hot spring in Kamchatka, Russia.

Desulfurococcus fermentans is the first known cellulolytic archaeon. This hyperthermophilic and strictly anaerobic crenarchaeon produces hydrogen from fermentation of various carbohydrates and peptides without inhibition by accumulating hydrogen. The complete genome sequence reported here suggested that D. fermentans employs membrane-bound hydrogenases and novel glycohydrolases for hydrogen production from cellulose.

Caldicellulosiruptor core and pangenomes reveal determinants for noncellulosomal thermophilic deconstruction of plant biomass.

Extremely thermophilic bacteria of the genus Caldicellulosiruptor utilize carbohydrate components of plant cell walls, including cellulose and hemicellulose, facilitated by a diverse set of glycoside hydrolases (GHs). From a biofuel perspective, this capability is crucial for deconstruction of plant biomass into fermentable sugars. While all species from the genus grow on xylan and acid-pretreated switchgrass, growth on crystalline cellulose is variable. The basis for this variability was examined using microbiological, genomic, and proteomic analyses of eight globally diverse Caldicellulosiruptor species. The open Caldicellulosiruptor pangenome (4,009 open reading frames [ORFs]) encodes 106 GHs, representing 43 GH families, but only 26 GHs from 17 families are included in the core (noncellulosic) genome (1,543 ORFs). Differentiating the strongly cellulolytic Caldicellulosiruptor species from the others is a specific genomic locus that encodes multidomain cellulases from GH families 9 and 48, which are associated with cellulose-binding modules. This locus also encodes a novel adhesin associated with type IV pili, which was identified in the exoproteome bound to crystalline cellulose. Taking into account the core genomes, pangenomes, and individual genomes, the ancestral Caldicellulosiruptor was likely cellulolytic and evolved, in some cases, into species that lost the ability to degrade crystalline cellulose while maintaining the capacity to hydrolyze amorphous cellulose and hemicellulose.

Complete genome sequence of the thermophilic, piezophilic, heterotrophic bacterium Marinitoga piezophila KA3.

Marinitoga piezophila KA3 is a thermophilic, anaerobic, chemoorganotrophic, sulfur-reducing bacterium isolated from the Grandbonum deep-sea hydrothermal vent site at the East Pacific Rise (13°N, 2,630-m depth). The genome of M. piezophila KA3 comprises a 2,231,407-bp circular chromosome and a 13,386-bp circular plasmid. This genome was sequenced within Department of Energy Joint Genome Institute CSP 2010.

Complete genome sequences of Desulfosporosinus orientis DSM765T, Desulfosporosinus youngiae DSM17734T, Desulfosporosinus meridiei DSM13257T, and Desulfosporosinus acidiphilus DSM22704T.

Desulfosporosinus species are sulfate-reducing bacteria belonging to the Firmicutes. Their genomes will give insights into the genetic repertoire and evolution of sulfate reducers typically thriving in terrestrial environments and able to degrade toluene (Desulfosporosinus youngiae), to reduce Fe(III) (Desulfosporosinus meridiei, Desulfosporosinus orientis), and to grow under acidic conditions (Desulfosporosinus acidiphilus).

Deep Learning Encoding for Rapid Sequence Identification on Microbiome Data.

We present a novel approach for rapidly identifying sequences that leverages the representational power of Deep Learning techniques and is applied to the analysis of microbiome data. The method involves the creation of a latent sequence space, training a convolutional neural network to rapidly identify sequences by mapping them into that space, and we leverage the novel encoded latent space for denoising to correct sequencing errors. Using mock bacterial communities of known composition, we show that this approach achieves single nucleotide resolution, generating results for sequence identification and abundance estimation that match the best available microbiome algorithms in terms of accuracy while vastly increasing the speed of accurate processing. We further show the ability of this approach to support phenotypic prediction at the sample level on an experimental data set for which the ground truth for sequence identities and abundances is unknown, but the expected phenotypes of the samples are definitive. Moreover, this approach offers a potential solution for the analysis of data from other types of experiments that currently rely on computationally intensive sequence identification.

Complete genome sequence of the anaerobic, halophilic alkalithermophile Natranaerobius thermophilus JW/NM-WN-LF.

The genome of the anaerobic halophilic alkalithermophile Natranaerobius thermophilus consists of one 3,165,557-bp chromosome and two plasmids (17,207 bp and 8,689 bp). The present study is the first to report the completely sequenced genome of an anaerobic polyextremophile and genes associated with roles in regulation of intracellular osmotic pressure, pH homeostasis, and growth at elevated temperatures.

Complete genome sequence of the haloalkaliphilic, hydrogen-producing bacterium Halanaerobium hydrogeniformans.

Halanaerobium hydrogenoformans is an alkaliphilic bacterium capable of biohydrogen production at pH 11 and 7% (wt/vol) salt. We present the 2.6-Mb genome sequence to provide insights into its physiology and potential for bioenergy applications.

Complete genome sequence of the cellulolytic thermophile Clostridium thermocellum DSM1313.

Clostridium thermocellum DSM1313 is a thermophilic, anaerobic bacterium with some of the highest rates of cellulose hydrolysis reported. The complete genome sequence reveals a suite of carbohydrate-active enzymes and demonstrates a level of diversity at the species level distinguishing it from the type strain ATCC 27405.

Genome sequence of Chthoniobacter flavus Ellin428, an aerobic heterotrophic soil bacterium.

Chthoniobacter flavus Ellin428 is the first isolate from the class Spartobacteria of the bacterial phylum Verrucomicrobia. C. flavus Ellin428 can metabolize many of the saccharide components of plant biomass but is incapable of growth on amino acids or organic acids other than pyruvate.

Genome sequences of Alicycliphilus denitrificans strains BC and K601T.

Alicycliphilus denitrificans strain BC and A. denitrificans strain K601(T) degrade cyclic hydrocarbons. These strains have been isolated from a mixture of wastewater treatment plant material and benzene-polluted soil and from a wastewater treatment plant, respectively, suggesting their role in bioremediation of soil and water. Although the strains are phylogenetically closely related, there are some clear physiological differences. The hydrocarbon cyclohexanol, for example, can be degraded by strain K601(T) but not by strain BC. Furthermore, both strains can use nitrate and oxygen as an electron acceptor, but only strain BC can use chlorate as electron acceptor. To better understand the nitrate and chlorate reduction mechanisms coupled to the oxidation of cyclic compounds, the genomes of A. denitrificans strains BC and K601(T) were sequenced. Here, we report the complete genome sequences of A. denitrificans strains BC and K601(T).

Genome sequence of the Arctic methanotroph Methylobacter tundripaludum SV96.

Methylobacter tundripaludum SV96(T) (ATCC BAA-1195) is a psychrotolerant aerobic methane-oxidizing gammaproteobacterium (Methylococcales, Methylococcaceae) living in High Arctic wetland soil. The strain was isolated from soil harvested in July 1996 close to the settlement Ny-Ålesund, Svalbard, Norway (78°56'N, 11°53'E), and described as a novel species in 2006. The genome includes pmo and pxm operons encoding copper membrane monooxygenases (Cu-MMOs), genes required for nitrogen fixation, and the nirS gene implicated in dissimilatory nitrite reduction to NO but no identifiable inventory for further processing of nitrogen oxides. These genome data provide the basis to investigate M. tundripaludum SV96, identified as a major player in the biogeochemistry of Arctic environments.

Complete genome sequences of Krokinobacter sp. strain 4H-3-7-5 and Lacinutrix sp. strain 5H-3-7-4, polysaccharide-degrading members of the family Flavobacteriaceae.

Two members of the family Flavobacteriaceae were isolated from subseafloor sediments using artificial seawater with cellulose, xylan, and chitin as the sole carbon and energy sources. Here, we present the complete genome sequences of Krokinobacter sp. strain 4H-3-7-5 and Lacinutrix sp. strain 5H-3-7-4, which both encode putatively novel enzymes involved in cellulose, hemicellulose, and chitin metabolism.

Genome sequence of the verrucomicrobium Opitutus terrae PB90-1, an abundant inhabitant of rice paddy soil ecosystems.

Bacteria of the deeply branching phylum Verrucomicrobia are rarely cultured yet commonly detected in metagenomic libraries from aquatic, terrestrial, and intestinal environments. We have sequenced the genome of Opitutus terrae PB90-1, a fermentative anaerobe within this phylum, isolated from rice paddy soil and capable of propionate production from plant-derived polysaccharides.

Complete genome sequence of the cellulose-degrading bacterium Cellulosilyticum lentocellum.

Cellulosilyticum lentocellum DSM 5427 is an anaerobic, endospore-forming member of the Firmicutes. We describe the complete genome sequence of this cellulose-degrading bacterium, which was originally isolated from estuarine sediment of a river that received both domestic and paper mill waste. Comparative genomics of cellulolytic clostridia will provide insight into factors that influence degradation rates.