Characterization of a biogas-producing microbial community by short-read next generation DNA sequencing

ABSTRACT: BACKGROUND: Renewable energy production is currently a major issue worldwide. Biogas is a promising renewable energy carrier as the technology of its production combines the elimination of organic waste with the formation of a versatile energy carrier, methane. In consequence of the complexity of the microbial communities and metabolic pathways involved the biotechnology of the microbiological process leading to biogas production is poorly understood. Metagenomic approaches are suitable means of addressing related questions. In the present work a novel high-throughput technique was tested for its benefits in resolving the functional and taxonomical complexity of such microbial consortia. RESULTS: It was demonstrated that the extremely parallel SOLiDTM short-read DNA sequencing platform is capable of providing sufficient useful information to decipher the systematic and functional contexts within a biogas-producing community. Although this technology has not been employed to address such problems previously, the data obtained compare well with those from similar high-throughput approaches such as 454-pyrosequencing GS FLX or Titanium. The predominant microbes contributing to the decomposition of organic matter include members of the Eubacteria, class Clostridia, order Clostridiales, family Clostridiaceae. Bacteria belonging in other systematic groups contribute to the diversity of the microbial consortium. Archaea comprise a remarkably small minority in this community, given their crucial role in biogas production. Among the Archaea, the predominant order is the Methanomicrobiales and the most abundant species is Methanoculleus marisnigri. The Methanomicrobiales are hydrogenotrophic methanogens. Besides corroborating earlier findings on the significance of the contribution of the Clostridia to organic substrate decomposition, the results demonstrate the importance of the metabolism of hydrogen within the biogas producing microbial community. CONCLUSIONS: Both microbiological diversity and the regulatory role of the hydrogen metabolism appear to be the driving forces optimizing biogas-producing microbial communities. The findings may allow a rational design of these communities to promote greater efficacy in large-scale practical systems. The composition of an optimal biogas-producing consortium can be determined through the use of this approach, and this systematic methodology allows the design of the optimal microbial community structure for any biogas plant. In this way, metagenomic studies can contribute to significant progress in the efficacy and economic improvement of biogas production.     

Complete genome sequence of Methanoculleus marisnigri Romesser et al. 1981 type strain JR1

Methanoculleus marisnigri Romesser et al. 1981 is a methanogen belonging to the order Methanomicrobiales within the archaeal phylum Euryarchaeota. The type strain, JR1, was isolated from anoxic sediments of the Black Sea. M. marisnigri is of phylogenetic interest because at the time the sequencing project began only one genome had previously been sequenced from the order Methanomicrobiales. We report here the complete genome sequence of M. marisnigri type strain JR1 and its annotation. This is part of a Joint Genome Institute 2006 Community Sequencing Program to sequence genomes of diverse Archaea.     

Complete Genome Sequence of the Hydrogenotrophic, Methanogenic Archaeon Methanoculleus bourgensis Strain MS2T, Isolated from a Sewage Sludge Digester

Methanoculleus bourgensis, of the order Methanomicrobiales, is a dominant methanogenic archaeon in many biogas-producing reactor systems fed with renewable primary products. It is capable of synthesizing methane via the hydrogenotrophic pathway utilizing hydrogen and carbon dioxide or formate as the substrates. Here we report the complete and finished genome sequence of M. bourgensis strain MS2(T), isolated from a sewage sludge digester.     

Characterization of microbial biofilms in a thermophilic biogas system by high-throughput metagenome sequencing

DNAs of two biofilms of a thermophilic two-phase leach-bed biogas reactor fed with rye silage and winter barley straw were sequenced by 454-pyrosequencing technology to assess the biofilm-based microbial community and their genetic potential for anaerobic digestion. The studied biofilms matured on the surface of the substrates in the hydrolysis reactor (HR) and on the packing in the anaerobic filter reactor (AF). The classification of metagenome reads showed Clostridium as most prevalent bacteria in the HR, indicating a predominant role for plant material digestion. Notably, insights into the genetic potential of plant-degrading bacteria were determined as well as further bacterial groups, which may assist Clostridium in carbohydrate degradation. Methanosarcina and Methanothermobacter were determined as most prevalent methanogenic archaea. In consequence, the biofilm-based methanogenesis in this system might be driven by the hydrogenotrophic pathway but also by the aceticlastic methanogenesis depending on metabolite concentrations such as the acetic acid concentration. Moreover, bacteria, which are capable of acetate oxidation in syntrophic interaction with methanogens, were also predicted. Finally, the metagenome analysis unveiled a large number of reads with unidentified microbial origin, indicating that the anaerobic degradation process may also be conducted by up to now unknown species.     

Taxonomic composition and gene content of a methane-producing microbial community isolated from a biogas reactor

A total community DNA sample from an agricultural biogas reactor continuously fed with maize silage, green rye, and small proportions of chicken manure has recently been sequenced using massively parallel pyrosequencing. In this study, the sample was computationally characterized without a prior assembly step, providing quantitative insights into the taxonomic composition and gene content of the underlying microbial community. Clostridiales from the phylum Firmicutes is the most prevalent phylogenetic order, Methanomicrobiales are dominant among methanogenic archaea. An analysis of Operational Taxonomic Units (OTUs) revealed that the entire microbial community is only partially covered by the sequenced sample, despite that estimates suggest only a moderate overall diversity of the community. Furthermore, the results strongly indicate that archaea related to the genus Methanoculleus, using CO(2) as electron acceptor and H(2) as electron donor, are the main producers of methane in the analyzed biogas reactor sample. A phylogenetic analysis of glycosyl hydrolase protein families suggests that Clostridia play an important role in the digestion of polysaccharides and oligosaccharides. Finally, the results unveiled that most of the organisms constituting the sample are still unexplored.     

Genome‐wide characterization and expression analysis of genetic variants in sweet orange

Heterozyosity is an important feature of many plant genomes, and is related to heterosis. Sweet orange, a highly heterozygous species, is thought to have originated from an inter‐species hybrid between pummelo and mandarin. To investigate the heterozygosity of the sweet orange genome and examine how this heterozygosity affects gene expression, we characterized the genome of Valencia orange for single nucleotide variations (SNVs), small insertions and deletions (InDels) and structural variations (SVs), and determined their functional effects on protein‐coding genes and non‐coding sequences. Almost half of the genes containing large‐effect SNVs and InDels were expressed in a tissue‐specific manner. We identified 3542 large SVs (>50 bp), including deletions, insertions and inversions. Most of the 296 genes located in large‐deletion regions showed low expression levels. RNA‐Seq reads and DNA sequencing reads revealed that the alleles of 1062 genes were differentially expressed. In addition, we detected approximately 42 Mb of contigs that were not found in the reference genome of a haploid sweet orange by de novo assembly of unmapped reads, and annotated 134 protein‐coding genes within these contigs. We discuss how this heterozygosity affects the quality of genome assembly. This study advances our understanding of the genome architecture of sweet orange, and provides a global view of gene expression at heterozygous loci.     

Local de novo assembly of RAD paired-end contigs using short sequencing reads

Despite the power of massively parallel sequencing platforms, a drawback is the short length of the sequence reads produced. We demonstrate that short reads can be locally assembled into longer contigs using paired-end sequencing of restriction-site associated DNA (RAD-PE) fragments. We use this RAD-PE contig approach to identify single nucleotide polymorphisms (SNPs) and determine haplotype structure in threespine stickleback and to sequence E. coli and stickleback genomic DNA with overlapping contigs of several hundred nucleotides. We also demonstrate that adding a circularization step allows the local assembly of contigs up to 5 kilobases (kb) in length. The ease of assembly and accuracy of the individual contigs produced from each RAD site sequence suggests RAD-PE sequencing is a useful way to convert genome-wide short reads into individually-assembled sequences hundreds or thousands of nucleotides long.     

Isolation of a methanogen from deep marine sediments that contain methane hydrates,and description of Methanoculleus submarinus sp. nov

We isolated a methanogen from deep in the sediments of the Nankai Trough off the eastern coast of Japan. At the sampling site, the water was 950 m deep and the sediment core was collected at 247 m below the sediment surface. The isolated methanogen was named Nankai-1. Cells of Nankai-1 were nonmotile and highly irregular coccoids (average diameter, 0.8 to 2 micro m) and grew with hydrogen or formate as a catabolic substrate. Cells required acetate as a carbon source. Yeast extract and peptones were not required but increased the growth rate. The cells were mesophilic, growing most rapidly at 45 degrees C (no growth at /=55 degrees C). Cells grew with a maximum specific growth rate of 2.43 day(-1) at 45 degrees C. Cells grew at pH values between 5.0 and 8.7 but did not grow at pH 4.7 or 9.0. Strain Nankai-1 grew in a wide range of salinities, from 0.1 to 1.5 M Na(+). The described phenotypic characteristics of this novel isolate were consistent with the in situ environment of the Nankai Trough. This is the first report of a methanogenic isolate from methane hydrate-bearing sediments. Phylogenetic analysis of its 16S rRNA gene sequence indicated that it is most closely related to Methanoculleus marisnigri (99.1% sequence similarity), but DNA hybridization experiments indicated a DNA sequence similarity of only 49%. Strain Nankai-1 was also found to be phenotypically similar to M. marisnigri, but two major phenotypic differences were found: strain Nankai-1 does not require peptones, and it grows fastest at a much higher temperature. We propose a new species, Methanoculleus submarinus, with strain Nankai-1 as the type strain.     

High-throughput pyrosequencing used for the discovery of a novel cellulase from a thermophilic cellulose-degrading microbial consortium

Objectives

To analyze the microbial diversity and gene content of a thermophilic cellulose-degrading consortium from hot springs in Xiamen, China using 454 pyrosequencing for discovering cellulolytic enzyme resources.

Results

A thermophilic cellulose-degrading consortium, XM70 that was isolated from a hot spring, used sugarcane bagasse as sole carbon and energy source. DNA sequencing of the XM70 sample resulted in 349,978 reads with an average read length of 380 bases, accounting for 133,896,867 bases of sequence information. The characterization of sequencing reads and assembled contigs revealed that most microbes were derived from four phyla: Geobacillus (Firmicutes), Thermus, Bacillus, and Anoxybacillus. Twenty-eight homologous genes belonging to 15 glycoside hydrolase families were detected, including several cellulase genes. A novel hot spring metagenome-derived thermophilic cellulase was expressed and characterized.

Conclusions

The application value of thermostable sugarcane bagasse-degrading enzymes is shown for production of cellulosic biofuel. The practical power of using a short-read-based metagenomic approach for harvesting novel microbial genes is also demonstrated.

     

Taxonomic Profiling and Metagenome Analysis of a Microbial Community from a Habitat Contaminated with Industrial Discharges

Industrial units, manufacturing dyes, chemicals, solvents, and xenobiotic compounds, produce liquid and solid wastes, which upon conventional treatment are released in the nearby environment and thus are the major cause of pollution. Soil collected from contaminated Kharicut Canal bank (N 22°57.878′; E 072°38.478′), Ahmedabad, Gujarat, India was used for metagenomic DNA preparation to study the capabilities of intrinsic microbial community in dealing with xenobiotics. Sequencing of metagenomic DNA on the Genome Sequencer FLX System using titanium chemistry resulted in 409,782 reads accounting for 133,529,997 bases of sequence information. Taxonomic analyses and gene annotations were carried out using the bioinformatics platform Sequence Analysis and Management System for Metagenomic Datasets. Taxonomic profiling was carried out by three different complementary approaches: (a) 16S rDNA, (b) environmental gene tags, and (c) lowest common ancestor. The most abundant phylum and genus were found to be “Proteobacteria” and “Pseudomonas,” respectively. Metagenome reads were mapped on sequenced microbial genomes and the highest numbers of reads were allocated to Pseudomonas stutzeri A1501. Assignment of obtained metagenome reads to Gene Ontology terms, Clusters of Orthologous Groups of protein categories, protein family numbers, and Kyoto Encyclopedia of Genes and Genomes hits revealed genomic potential of indigenous microbial community. In total, 157,024 reads corresponded to 37,028 different KEGG hits, and amongst them, 11,574 reads corresponded to 131 different enzymes potentially involved in xenobiotic biodegradation. These enzymes were mapped on biodegradation pathways of xenobiotics to elucidate their roles in possible catalytic reactions. Consequently, information obtained from the present study will act as a baseline which, subsequently along with other “-omic” studies, will help in designing future bioremediation strategies in effluent treatment plants and environmental clean-up projects.     

Quantitative metagenomic analyses based on average genome size normalization

Over the past quarter-century, microbiologists have used DNA sequence information to aid in the characterization of microbial communities. During the last decade, this has expanded from single genes to microbial community genomics, or metagenomics, in which the gene content of an environment can provide not just a census of the community members but direct information on metabolic capabilities and potential interactions among community members. Here we introduce a method for the quantitative characterization and comparison of microbial communities based on the normalization of metagenomic data by estimating average genome sizes. This normalization can relieve comparative biases introduced by differences in community structure, number of sequencing reads, and sequencing read lengths between different metagenomes. We demonstrate the utility of this approach by comparing metagenomes from two different marine sources using both conventional small-subunit (SSU) rRNA gene analyses and our quantitative method to calculate the proportion of genomes in each sample that are capable of a particular metabolic trait. With both environments, to determine what proportion of each community they make up and how differences in environment affect their abundances, we characterize three different types of autotrophic organisms: aerobic, photosynthetic carbon fixers (the Cyanobacteria); anaerobic, photosynthetic carbon fixers (the Chlorobi); and anaerobic, nonphotosynthetic carbon fixers (the Desulfobacteraceae). These analyses demonstrate how genome proportionality compares to SSU rRNA gene relative abundance and how factors such as average genome size and SSU rRNA gene copy number affect sampling probability and therefore both types of community analysis.     

Direct Chloroplast Sequencing: Comparison of Sequencing Platforms and Analysis Tools for Whole Chloroplast Barcoding

Direct sequencing of total plant DNA using next generation sequencing technologies generates a whole chloroplast genome sequence that has the potential to provide a barcode for use in plant and food identification. Advances in DNA sequencing platforms may make this an attractive approach for routine plant identification. The HiSeq (Illumina) and Ion Torrent (Life Technology) sequencing platforms were used to sequence total DNA from rice to identify polymorphisms in the whole chloroplast genome sequence of a wild rice plant relative to cultivated rice (cv. Nipponbare). Consensus chloroplast sequences were produced by mapping sequence reads to the reference rice chloroplast genome or by de novo assembly and mapping of the resulting contigs to the reference sequence. A total of 122 polymorphisms (SNPs and indels) between the wild and cultivated rice chloroplasts were predicted by these different sequencing and analysis methods. Of these, a total of 102 polymorphisms including 90 SNPs were predicted by both platforms. Indels were more variable with different sequencing methods, with almost all discrepancies found in homopolymers. The Ion Torrent platform gave no apparent false SNP but was less reliable for indels. The methods should be suitable for routine barcoding using appropriate combinations of sequencing platform and data analysis.     

Insight into Dominant Cellulolytic Bacteria from Two Biogas Digesters and Their Glycoside Hydrolase Genes

Diverse cellulolytic bacteria are essential for maintaining high lignocellulose degradation ability in biogas digesters. However, little was known about functional genes and gene clusters of dominant cellulolytic bacteria in biogas digesters. This is the foundation to understand lignocellulose degradation mechanisms of biogas digesters and apply these gene resource for optimizing biofuel production. A combination of metagenomic and 16S rRNA gene clone library methods was used to investigate the dominant cellulolytic bacteria and their glycoside hydrolase (GH) genes in two biogas digesters. The 16S rRNA gene analysis revealed that the dominant cellulolytic bacteria were strains closely related to Clostridium straminisolvens and an uncultured cellulolytic bacterium designated BG-1. To recover GH genes from cellulolytic bacteria in general, and BG-1 in particular, a refined assembly approach developed in this study was used to assemble GH genes from metagenomic reads; 163 GH-containing contigs ≥ 1 kb in length were obtained. Six recovered GH5 genes that were expressed in E. coli demonstrated multiple lignocellulase activities and one had high mannanase activity (1255 U/mg). Eleven fosmid clones harboring the recovered GH-containing contigs were sequenced and assembled into 10 fosmid contigs. The composition of GH genes in the 163 assembled metagenomic contigs and 10 fosmid contigs indicated that diverse GHs and lignocellulose degradation mechanisms were present in the biogas digesters. In particular, a small portion of BG-1 genome information was recovered by PhyloPythiaS analysis. The lignocellulase gene clusters in BG-1 suggested that it might use a possible novel lignocellulose degradation mechanism to efficiently degrade lignocellulose. Dominant cellulolytic bacteria of biogas digester possess diverse GH genes, not only in sequences but also in their functions, which may be applied for production of biofuel in the future.     

Methanoculleus spp. as a biomarker of methanogenic activity in swine manure storage tanks

Greenhouse gas emissions represent a major problem associated with manure management in the livestock industry. A prerequisite to mitigate methane emissions occurring during manure storage is a clearer understanding of how the microbial consortia involved in methanogenesis function. Here, we have examined manure stored in outdoor tanks from two different farms, at different locations and depths. Physico-chemical and microbiological characterization of these samples indicated differences between each tank, as well as differences within each tank dependent on the depth of sampling. The dynamics of both the bacterial and archaeal communities within these samples were monitored over a 150-day period of anaerobic incubation to identify and track emerging microorganisms, which may be temporally important in the methanogenesis process. Analyses based on DNA fingerprinting of microbial communities identified trends common among all samples as well as trends specific to certain samples. All archaeal communities became enriched with Methanoculleus spp. over time, indicating that the hydrogenotrophic pathway of methanogenesis predominated. Although the emerging species differed in samples obtained from shallow depths compared to deep samples, the temporal enrichment of Methanoculleus suggests that this genus may represent a relevant indicator of methanogenic activity in swine manure storage tanks.     

Rapid Genome Mapping in Nanochannel Arrays for Highly Complete and Accurate De Novo Sequence Assembly of the Complex Aegilops tauschii Genome

Next-generation sequencing (NGS) technologies have enabled high-throughput and low-cost generation of sequence data; however, de novo genome assembly remains a great challenge, particularly for large genomes. NGS short reads are often insufficient to create large contigs that span repeat sequences and to facilitate unambiguous assembly. Plant genomes are notorious for containing high quantities of repetitive elements, which combined with huge genome sizes, makes accurate assembly of these large and complex genomes intractable thus far. Using two-color genome mapping of tiling bacterial artificial chromosomes (BAC) clones on nanochannel arrays, we completed high-confidence assembly of a 2.1-Mb, highly repetitive region in the large and complex genome of Aegilops tauschii, the D-genome donor of hexaploid wheat (Triticum aestivum). Genome mapping is based on direct visualization of sequence motifs on single DNA molecules hundreds of kilobases in length. With the genome map as a scaffold, we anchored unplaced sequence contigs, validated the initial draft assembly, and resolved instances of misassembly, some involving contigs <2 kb long, to dramatically improve the assembly from 75% to 95% complete.