Microdiversity of extracellular enzyme genes among sequenced prokaryotic genomes

Coupled high-throughput functional screening and next generation sequencing for identification of plant polymer decomposing enzymes in metagenomic libraries

Recent advances in sequencing technologies generate new predictions and hypotheses about the functional roles of environmental microorganisms. Yet, until we can test these predictions at a scale that matches our ability to generate them, most of them will remain as hypotheses. Function-based mining of metagenomic libraries can provide direct linkages between genes, metabolic traits and microbial taxa and thus bridge this gap between sequence data generation and functional predictions. Here we developed high-throughput screening assays for function-based characterization of activities involved in plant polymer decomposition from environmental metagenomic libraries. The multiplexed assays use fluorogenic and chromogenic substrates, combine automated liquid handling and use a genetically modified expression host to enable simultaneous screening of 12,160 clones for 14 activities in a total of 170,240 reactions. Using this platform we identified 374 (0.26%) cellulose, hemicellulose, chitin, starch, phosphate and protein hydrolyzing clones from fosmid libraries prepared from decomposing leaf litter. Sequencing on the Illumina MiSeq platform, followed by assembly and gene prediction of a subset of 95 fosmid clones, identified a broad range of bacterial phyla, including Actinobacteria, Bacteroidetes, multiple Proteobacteria sub-phyla in addition to some Fungi. Carbohydrate-active enzyme genes from 20 different glycoside hydrolase (GH) families were detected. Using tetranucleotide frequency (TNF) binning of fosmid sequences, multiple enzyme activities from distinct fosmids were linked, demonstrating how biochemically-confirmed functional traits in environmental metagenomes may be attributed to groups of specific organisms. Overall, our results demonstrate how functional screening of metagenomic libraries can be used to connect microbial functionality to community composition and, as a result, complement large-scale metagenomic sequencing efforts.

Microbial community structure of a pilot-scale thermophilic anaerobic digester treating poultry litter

The microbial community structure of a stable pilot-scale thermophilic continuous stirred tank reactor digester stabilized on poultry litter was investigated. This 40-m(3) digester produced biogas with 57% methane, and chemical oxygen demand removal of 54%. Bacterial and archaeal diversity were examined using both cloning and pyrosequencing that targeted 16S rRNA genes. The bacterial community was dominated by phylum Firmicutes, constituting 93% of the clones and 76% of the pyrotags. Of the Firmicutes, class Clostridia (52% pyrotags) was most abundant followed by class Bacilli (13% pyrotags). The bacterial libraries identified 94 operational taxonomic units (OTUs) and pyrosequencing identified 577 OTUs at the 97% minimum similarity level. Fifteen OTUs were dominant (≥2% abundance), and nine of these were novel unclassified Firmicutes. Several of the dominant OTUs could not be classified more specifically than Clostridiales, but were most similar to plant biomass degraders, including Clostridium thermocellum. Of the rare pyrotag OTUs (<0.5% abundance), 75% were Firmicutes. The dominant methanogen was Methanothermobacter which has hydrogenotrophic metabolism, and accounted for >99% of the archaeal clones. Based on the primary methanogen, as well as digester chemistry (high VA and ammonia levels), we propose that bacterial acetate oxidation is the primary pathway in this digester for the control of acetate levels.

Phylogenetic distribution of potential cellulases in bacteria

Many microorganisms contain cellulases that are important for plant cell wall degradation and overall soil ecosystem functioning. At present, we have extensive biochemical knowledge of cellulases but little is known about the phylogenetic distribution of these enzymes. To address this, we analyzed the distribution of 21,985 genes encoding proteins related to cellulose utilization in 5,123 sequenced bacterial genomes. First, we identified the distribution of glycoside hydrolases involved in cellulose utilization and synthesis at different taxonomic levels, from the phylum to the strain. Cellulose degradation/utilization capabilities appeared in nearly all major groups and resulted in strains displaying various enzyme gene combinations. Potential cellulose degraders, having both cellulases and β-glucosidases, constituted 24% of all genomes whereas potential opportunistic strains, having β-glucosidases only, accounted for 56%. Finally, 20% of the bacteria have no relevant enzymes and do not rely on cellulose utilization. The latter group was primarily connected to specific bacterial lifestyles like autotrophy and parasitism. Cellulose degraders, as well as opportunists, have multiple enzymes with similar functions. However, the potential degraders systematically harbor about twice more β-glucosidases than their potential opportunistic relatives. Although scattered, the distribution of functional types, in bacterial lineages, is not random but mostly follows a Brownian motion evolution model. Degraders form clusters of relatives at the species level, whereas opportunists are clustered at the genus level. This information can form a mechanistic basis for the linking of changes in microbial community composition to soil ecosystem processes.