Josh D. Neufeld

Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads.

ABSTRACT Microbial communities host unparalleled taxonomic diversity. Adequate characterization of environmental and host-associated samples remains a challenge for microbiologists, despite the advent of 16S rRNA gene sequencing. In order to increase the depth of sampling for diverse bacterial communities, we developed a method for sequencing and assembling millions of paired-end reads from the 16S rRNA gene (spanning the V3 region; ∼200 nucleotides) by using an Illumina genome analyzer. To confirm reproducibility and to identify a suitable computational pipeline for data analysis, sequence libraries were prepared in duplicate for both a defined mixture of DNAs from known cultured bacterial isolates (>1 million postassembly sequences) and an Arctic tundra soil sample (>6 million postassembly sequences). The Illumina 16S rRNA gene libraries represent a substantial increase in number of sequences over all extant next-generation sequencing approaches (e.g., 454 pyrosequencing), while the assembly of paired-end 125-base reads offers a methodological advantage by incorporating an initial quality control step for each 16S rRNA gene sequence. This method incorporates indexed primers to enable the characterization of multiple microbial communities in a single flow cell lane, may be modified readily to target other variable regions or genes, and demonstrates unprecedented and economical access to DNAs from organisms that exist at low relative abundances.

Bacterial biogeography of the human digestive tract

We present bacterial biogeography as sampled from the human gastrointestinal tract of four healthy subjects. This study generated >32 million paired-end sequences of bacterial 16S rRNA genes (V3 region) representing >95,000 unique operational taxonomic units (OTUs; 97% similarity clusters), with >99% Goods coverage for all samples. The highest OTU richness and phylogenetic diversity was found in the mouth samples. The microbial communities of multiple biopsy sites within the colon were highly similar within individuals and largely distinct from those in stool. Within an individual, OTU overlap among broad site definitions (mouth, stomach/duodenum, colon and stool) ranged from 32–110 OTUs, 25 of which were common to all individuals and included OTUs affiliated with Faecalibacterium prasnitzii and the TM7 phylum. This first comprehensive characterization of the abundant and rare microflora found along the healthy human digestive tract represents essential groundwork to investigate further how the human microbiome relates to health and disease.

Who eats what, where and when? Isotope-labelling experiments are coming of age

Isotope-labelling experiments have changed the way microbial ecologists investigate the ecophysiology of microbial populations and cells in the environment. Insight into the ‘uncultivated majority’ accompanies methodology that involves the incorporation of stable isotopes or radioisotopes into sub-populations of environmental samples. Subsequent analysis of labelled biomarkers of sub-populations with stable-isotope probing (DNA-SIP, RNA-SIP, phospholipid-derived fatty acid-SIP) or individual cells with a combination of fluorescence in situ hybridization and microautoradiography reveals linked phylogenetic and functional information about the organisms that assimilated these compounds. Here, we review some of the most recent literature, with an emphasis on methodological improvements to the sensitivity and utility of these methods. We also highlight related isotope techniques that are in continued development and hold promise to transform the way we link phylogeny and function in complex microbial communities.

Methodological Considerations for the Use of Stable Isotope Probing in Microbial Ecology

Stable isotope probing (SIP) is a method used for labeling uncultivated microorganisms in environmental samples or directly in field studies using substrate enriched with stable isotope (e.g., 13C). After consumption of the substrate, the cells of microorganisms that consumed the substrate become enriched in the isotope. Labeled biomarkers, such as phospholipid-derived fatty acid (PLFA), ribosomal RNA, and DNA can be analyzed with a range of molecular and analytical techniques, and used to identify and characterize the organisms that incorporated the substrate. The advantages and disadvantages of PLFA-SIP, RNA-SIP, and DNA-SIP are presented. Using examples from our laboratory and from the literature, we discuss important methodological considerations for a successful SIP experiment.

Something from (almost) nothing: the impact of multiple displacement amplification on microbial ecology

Microbial ecology is a field that applies molecular techniques to analyze genes and communities associated with a plethora of unique environments on this planet. In the past, low biomass and the predominance of a few abundant community members have impeded the application of techniques such as PCR, microarray analysis and metagenomics to complex microbial populations. In the absence of suitable cultivation methods, it was not possible to obtain DNA samples from individual microorganisms. Recently, a method called multiple displacement amplification (MDA) has been used to circumvent these limitations by amplifying DNA from microbial communities in low-biomass environments, individual cells from uncultivated microbial species and active organisms obtained through stable isotope probing incubations. This review describes the development and applications of MDA, discusses its strengths and limitations and highlights the impact of MDA on the field of microbial ecology. Whole genome amplification via MDA has increased access to the genomic DNA of uncultivated microorganisms and low-biomass environments and represents a ‘power tool’ in the molecular toolbox of microbial ecologists.

Unexpectedly High Bacterial Diversity in Arctic Tundra Relative to Boreal Forest Soils, Revealed by Serial Analysis of Ribosomal Sequence Tags

ABSTRACT Arctic tundra and boreal forest soils have globally relevant functions that affect atmospheric chemistry and climate, yet the bacterial composition and diversity of these soils have received little study. Serial analysis of ribosomal sequence tags (SARST) and denaturing gradient gel electrophoresis (DGGE) were used to compare composite soil samples taken from boreal and arctic biomes. This study comprises an extensive comparison of geographically distant soil bacterial communities, involving the analysis of 12,850 ribosomal sequence tags from six composite soil samples. Bacterial diversity estimates were greater for undisturbed arctic tundra soil samples than for boreal forest soil samples, with the highest diversity associated with a sample from an extreme northern location (82oN). The lowest diversity estimate was obtained from an arctic soil sample that was disturbed by compaction and sampled from a greater depth. Since samples from the two biomes did not form distinct clusters on the basis of SARST data and DGGE fingerprints, factors other than latitude likely influenced the phylogenetic compositions of these communities. The high number of ribosomal sequences analyzed enabled the identification of possible cosmopolitan and endemic bacterial distributions in particular soils.

Stable-isotope probing implicates Methylophaga spp and novel Gammaproteobacteria in marine methanol and methylamine metabolism

The metabolism of one-carbon (C1) compounds in the marine environment affects global warming, seawater ecology and atmospheric chemistry. Despite their global significance, marine microorganisms that consume C1 compounds in situ remain poorly characterized. Stable-isotope probing (SIP) is an ideal tool for linking the function and phylogeny of methylotrophic organisms by the metabolism and incorporation of stable-isotope-labelled substrates into nucleic acids. By combining DNA-SIP and time-series sampling, we characterized the organisms involved in the assimilation of methanol and methylamine in coastal sea water (Plymouth, UK). Labelled nucleic acids were analysed by denaturing gradient gel electrophoresis (DGGE) and clone libraries of 16S rRNA genes. In addition, we characterized the functional gene complement of labelled nucleic acids with an improved primer set targeting methanol dehydrogenase (mxaF) and newly designed primers for methylamine dehydrogenase (mauA). Predominant DGGE phylotypes, 16S rRNA, methanol and methylamine dehydrogenase gene sequences, and cultured isolates all implicated Methylophaga spp, moderately halophilic marine methylotrophs, in the consumption of both methanol and methylamine. Additionally, an mxaF sequence obtained from DNA extracted from sea water clustered with those detected in 13C-DNA, suggesting a predominance of Methylophaga spp among marine methylotrophs. Unexpectedly, most predominant 16S rRNA and functional gene sequences from 13C-DNA were clustered in distinct substrate-specific clades, with 16S rRNA genes clustering with sequences from the Gammaproteobacteria. These clades have no cultured representatives and reveal an ecological adaptation of particular uncultured methylotrophs to specific C1 compounds in the coastal marine environment.

Marine methylotrophs revealed by stable-isotope probing, multiple displacement amplification and metagenomics

The concentrations of one-carbon substrates that fuel methylotrophic microbial communities in the ocean are limited and the specialized guilds of bacteria that use these molecules may exist at low relative abundance. As a result, these organisms are difficult to identify and are often missed with existing cultivation and gene retrieval methods. Here, we demonstrate a novel proof of concept: using environmentally-relevant substrate concentrations in stable-isotope probing (SIP) incubations to yield sufficient DNA for large-insert metagenomic analysis through multiple displacement amplification (MDA). A marine surface-water sample was labelled sufficiently by incubation with near in situ concentrations of methanol. Picogram quantities of labelled (13)C-DNA were purified from caesium chloride gradients, amplified with MDA to produce microgram amounts of high-molecular-weight DNA (<or= 40 kb) and cloned to produce a fosmid library of > 10 000 clones. Denaturing gradient gel electrophoresis (DGGE) demonstrated minimal bias associated with the MDA step and implicated Methylophaga-like phylotypes with the marine metabolism of methanol. Polymerase chain reaction screening of 1500 clones revealed a methanol dehydrogenase (MDH) containing insert and shotgun sequencing of this insert resulted in the assembly of a 9-kb fragment of DNA encoding a cluster of enzymes involved in MDH biosynthesis, regulation and assembly. This novel combination of methodology enables future structure-function studies of microbial communities to achieve the long-desired goal of identifying active microbial populations using in situ conditions and performing a directed metagenomic analysis for these ecologically relevant microorganisms.

Revealing the uncultivated majority: combining DNA stable-isotope probing, multiple displacement amplification and metagenomic analyses of uncultivated Methylocystis in acidic peatlands

Peatlands represent an enormous carbon reservoir and have a potential impact on the global climate because of the active methanogenesis and methanotrophy in these soils. Uncultivated methanotrophs from seven European peatlands were studied using a combination of molecular methods. Screening for methanotroph diversity using a particulate methane monooxygenase-based diagnostic gene array revealed that Methylocystis-related species were dominant in six of the seven peatlands studied. The abundance and methane oxidation activity of Methylocystis spp. were further confirmed by DNA stable-isotope probing analysis of a sample taken from the Moor House peatland (England). After ultracentrifugation, (13)C-labelled DNA, containing genomic DNA of these Methylocystis spp., was separated from (12)C DNA and subjected to multiple displacement amplification (MDA) to generate sufficient DNA for the preparation of a fosmid metagenomic library. Potential bias of MDA was detected by fingerprint analysis of 16S rRNA using denaturing gradient gel electrophoresis for low-template amplification (0.01 ng template). Sufficient template (1-5 ng) was used in MDA to circumvent this bias and chimeric artefacts were minimized by using an enzymatic treatment of MDA-generated DNA with S1 nuclease and DNA polymerase I. Screening of the metagenomic library revealed one fosmid containing methanol dehydrogenase and two fosmids containing 16S rRNA genes from these Methylocystis-related species as well as one fosmid containing a 16S rRNA gene related to that of Methylocella/Methylocapsa. Sequencing of the 14 kb methanol dehydrogenase-containing fosmid allowed the assembly of a gene cluster encoding polypeptides involved in bacterial methanol utilization (mxaFJGIRSAC). This combination of DNA stable-isotope probing, MDA and metagenomics provided access to genomic information of a relatively large DNA fragment of these thus far uncultivated, predominant and active methanotrophs in peatland soil.

Characterization of Humus Microbial Communities in Adjacent Forest Types That Differ in Nitrogen Availability

To address the link between soil microbial community composition and soil processes, we investigated the microbial communities in forest floors of two forest types that differ substantially in nitrogen availability. Cedar-hemlock (CH) and hemlock-amabilis fir (HA) forests are both common on northern Vancouver Island, B.C., occurring adjacently across the landscape. CH forest floors have low nitrogen availability and HA high nitrogen availability. Total microbial biomass was assessed using chloroform fumigation-extraction and community composition was assessed using several cultivation-independent approaches: denaturing gradient gel electrophoresis (DGGE) of the bacterial communities, ribosomal intergenic spacer analysis (RISA) of the bacterial and fungal communities, and phospholipid fatty acid (PLFA) profiles of the whole microbial community. We did not detect differences in the bacterial communities of each forest type using DGGE and RISA, but differences in the fungal communities were detected using RISA. PLFA analysis detected subtle differences in overall composition of the microbial community between the forest types, as well as in particular groups of organisms. Fungal PLFAs were more abundant in the nitrogen-poor CH forests. Bacteria were proportionally more abundant in HA forests than CH in the lower humus layer, and Gram-positive bacteria were proportionally more abundant in HA forests irrespective of layer. Bacterial and fungal communities were distinct in the F, upper humus, and lower humus layers of the forest floor and total biomass decreased in deeper layers. These results indicate that there are distinct patterns in forest floor microbial community composition at the landscape scale, which may be important for understanding nutrient availability to forest vegetation.