Review, Evaluation, and Directions for Gene-Targeted Assembly for Ecological Analyses of Metagenomes.

Shotgun metagenomics has greatly advanced our understanding of microbial communities over the last decade. Metagenomic analyses often include assembly and genome binning, computationally daunting tasks especially for big data from complex environments such as soil and sediments. In many studies, however, only a subset of genes and pathways involved in specific functions are of interest; thus, it is not necessary to attempt global assembly. In addition, methods that target genes can be computationally more efficient and produce more accurate assembly by leveraging rich databases, especially for those genes that are of broad interest such as those involved in biogeochemical cycles, biodegradation, and antibiotic resistance or used as phylogenetic markers. Here, we review six gene-targeted assemblers with unique algorithms for extracting and/or assembling targeted genes: Xander, MegaGTA, SAT-Assembler, HMM-GRASPx, GenSeed-HMM, and MEGAN. We tested these tools using two datasets with known genomes, a synthetic community of artificial reads derived from the genomes of 17 bacteria, shotgun sequence data from a mock community with 48 bacteria and 16 archaea genomes, and a large soil shotgun metagenomic dataset. We compared assemblies of a universal single copy gene (rplB) and two N cycle genes (nifH and nirK). We measured their computational efficiency, sensitivity, specificity, and chimera rate and found Xander and MegaGTA, which both use a probabilistic graph structure to model the genes, have the best overall performance with all three datasets, although MEGAN, a reference matching assembler, had better sensitivity with synthetic and mock community members chosen from its reference collection. Also, Xander and MegaGTA are the only tools that include post-assembly scripts tuned for common molecular ecology and diversity analyses. Additionally, we provide a mathematical model for estimating the probability of assembling targeted genes in a metagenome for estimating required sequencing depth.

Ecological patterns of nifH genes in four terrestrial climatic zones explored with targeted metagenomics using FrameBot, a new informatics tool.

UNLABELLED: Biological nitrogen fixation is an important component of sustainable soil fertility and a key component of the nitrogen cycle. We used targeted metagenomics to study the nitrogen fixation-capable terrestrial bacterial community by targeting the gene for nitrogenase reductase (nifH). We obtained 1.1 million nifH 454 amplicon sequences from 222 soil samples collected from 4 National Ecological Observatory Network (NEON) sites in Alaska, Hawaii, Utah, and Florida. To accurately detect and correct frameshifts caused by indel sequencing errors, we developed FrameBot, a tool for frameshift correction and nearest-neighbor classification, and compared its accuracy to that of two other rapid frameshift correction tools. We found FrameBot was, in general, more accurate as long as a reference protein sequence with 80% or greater identity to a query was available, as was the case for virtually all nifH reads for the 4 NEON sites. Frameshifts were present in 12.7% of the reads. Those nifH sequences related to the Proteobacteria phylum were most abundant, followed by those for Cyanobacteria in the Alaska and Utah sites. Predominant genera with nifH sequences similar to reads included Azospirillum, Bradyrhizobium, and Rhizobium, the latter two without obvious plant hosts at the sites. Surprisingly, 80% of the sequences had greater than 95% amino acid identity to known nifH gene sequences. These samples were grouped by site and correlated with soil environmental factors, especially drainage, light intensity, mean annual temperature, and mean annual precipitation. FrameBot was tested successfully on three ecofunctional genes but should be applicable to any. IMPORTANCE: High-throughput phylogenetic analysis of microbial communities using rRNA-targeted sequencing is now commonplace; however, such data often allow little inference with respect to either the presence or the diversity of genes involved in most important ecological processes. To study the gene pool for these processes, it is more straightforward to assess the genes directly responsible for the ecological function (ecofunctional genes). However, analyzing these genes involves technical challenges beyond those seen for rRNA. In particular, frameshift errors cause garbled downstream protein translations. Our FrameBot tool described here both corrects frameshift errors in query reads and determines their closest matching protein sequences in a set of reference sequences. We validated this new tool with sequences from defined communities and demonstrated the tool's utility on nifH gene fragments sequenced from soils in well-characterized and major terrestrial ecosystem types.

Effect of storage conditions on the stability and fermentability of enzymatic lignocellulosic hydrolysate.

To minimize the change of lignocellulosic hydrolysate composition during storage, the effects of storage conditions (temperature, pH and time) on the composition and fermentability of hydrolysate prepared from AFEX™ (Ammonia Fiber Expansion - a trademark of MBI, Lansing, MI) pretreated corn stover were investigated. Precipitates formed during hydrolysate storage increased with increasing storage pH and time. The precipitate amount was the least when hydrolysate was stored at 4 °C and pH 4.8, accounting for only 0.02% of the total hydrolysate weight after 3-month storage. No significant changes of NMR (Nuclear Magnetic Resonance) spectra and concentrations of sugars, minerals and heavy metals were observed after storage under this condition. When pH was adjusted higher before fermentation, precipitates also formed, consisting of mostly struvite (MgNH4PO4·6H2O) and brushite (CaHPO4·2H2O). Escherichia coli and Saccharomyces cerevisiae fermentation studies and yeast cell growth assays showed no significant difference in fermentability between fresh hydrolysate and stored hydrolysate.

Functional genes to assess nitrogen cycling and aromatic hydrocarbon degradation: primers and processing matter.

Targeting sequencing to genes involved in key environmental processes, i.e., ecofunctional genes, provides an opportunity to sample nature's gene guilds to greater depth and help link community structure to process-level outcomes. Vastly different approaches have been implemented for sequence processing and, ultimately, for taxonomic placement of these gene reads. The overall quality of next generation sequence analysis of functional genes is dependent on multiple steps and assumptions of unknown diversity. To illustrate current issues surrounding amplicon read processing we provide examples for three ecofunctional gene groups. A combination of in silico, environmental and cultured strain sequences was used to test new primers targeting the dioxin and dibenzofuran degrading genes dxnA1, dbfA1, and carAa. The majority of obtained environmental sequences were classified into novel sequence clusters, illustrating the discovery value of the approach. For the nitrite reductase step in denitrification, the well-known nirK primers exhibited deficiencies in reference database coverage, illustrating the need to refine primer-binding sites and/or to design multiple primers, while nirS primers exhibited bias against five phyla. Amino acid-based OTU clustering of these two N-cycle genes from soil samples yielded only 114 unique nirK and 45 unique nirS genus-level groupings, likely a reflection of constricted primer coverage. Finally, supervised and non-supervised OTU analysis methods were compared using the nifH gene of nitrogen fixation, with generally similar outcomes, but the clustering (non-supervised) method yielded higher diversity estimates and stronger site-based differences. High throughput amplicon sequencing can provide inexpensive and rapid access to nature's related sequences by circumventing the culturing barrier, but each unique gene requires individual considerations in terms of primer design and sequence processing and classification.

Microbial Communities Associated with Potato Common Scab-Suppressive Soil Determined by Pyrosequencing Analyses.

Potato common scab, caused by Streptomyces spp., is an annual production problem for potato growers, and not effectively controlled by current methods. A field with naturally occurring common scab suppression has been identified in Michigan, and confirmed to have a biological basis for this disease suppression. This field and an adjacent scab nursery conducive to disease were studied using pyrosequencing to compare the two microbial communities. Total DNA was extracted from both the disease-conducive and -suppressive soils. A phylogenetically taxon-informative region of the 16S rRNA gene was used to establish operational taxonomic units (OTUs) to characterize bacterial community richness and diversity. In total, 1,124 OTUs were detected and 565 OTUs (10% dissimilarity) were identified in disease-conducive soil and 859 in disease-suppressive soil, including 300 shared both between sites. Common phyla based on relative sequence abundance were Acidobacteria, Proteobacteria, and Firmicutes. Sequences of Lysobacter were found in significantly higher numbers in the disease-suppressive soil, as were sequences of group 4 and group 6 Acidobacteria. The relative abundance of sequences identified as the genus Bacillus was significantly higher by an order of magnitude in the disease-conducive soil.

DNA-stable isotope probing integrated with metagenomics for retrieval of biphenyl dioxygenase genes from polychlorinated biphenyl-contaminated river sediment.

Stable isotope probing with [(13)C]biphenyl was used to explore the genetic properties of indigenous bacteria able to grow on biphenyl in PCB-contaminated River Raisin sediment. A bacterial 16S rRNA gene clone library generated from [(13)C]DNA after a 14-day incubation with [(13)C]biphenyl revealed the dominant organisms to be members of the genera Achromobacter and Pseudomonas. A library built from PCR amplification of genes for aromatic-ring-hydroxylating dioxygenases from the [(13)C]DNA fraction revealed two sequence groups similar to bphA (encoding biphenyl dioxygenase) of Comamonas testosteroni strain B-356 and of Rhodococcus sp. RHA1. A library of 1,568 cosmid clones was produced from the [(13)C]DNA fraction. A 31.8-kb cosmid clone, detected by aromatic dioxygenase primers, contained genes of biphenyl dioxygenase subunits bphAE, while the rest of the clone's sequence was similar to that of an unknown member of the Gammaproteobacteria. A discrepancy in G+C content near the bphAE genes implies their recent acquisition, possibly by horizontal transfer. The biphenyl dioxygenase from the cosmid clone oxidized biphenyl and unsubstituted and para-only-substituted rings of polychlorinated biphenyl (PCB) congeners. A DNA-stable isotope probing-based cosmid library enabled the retrieval of functional genes from an uncultivated organism capable of PCB metabolism and suggest dispersed dioxygenase gene organization in nature.

Stimulatory effects of a microbially dechlorinated polychlorinated biphenyl (PCB) mixture on rat uterine contraction in vitro.

Microbially mediated reductive dechlorination has been advocated as the first part of a two-stage (anaerobic/aerobic) biotreatment process for polychlorinated biphenyls (PCBs) in sediments, and is generally viewed as a detoxication process. However, previous studies suggest that microbial dechlorination increases the ability to stimulate uterine contractions compared with the original PCB mixtures. Here, we investigate the composition and uterotonic activity of the commercial PCB mixture Aroclor 1260 before and after incubation with microorganisms eluted from PCB-contaminated sediment of the Hudson River. Incubation with microorganisms resulted in a partially dechlorinated mixture (HR1260) dominated by ortho-substituted PCBs with four or fewer chlorines per biphenyl. Aroclor 1260 that had not been incubated with microorganisms had no significant effect on contraction frequency of rat uterine strips (gestation day 10) in vitro, whereas HR1260 dramatically increased contraction frequency to 718+/-134% of the basal rate at a total PCB concentration of 70 microM (p<0.05). The microbial dechlorination increased 2,2',4,4'-tetrachlorobiphenyl and one or more of four congeners that co-eluted during chromatography (2,3,3',5-tetrachlorobiphenyl, 2,3',4,5-tetrachlorobiphenyl, 2,2',4,4',6-pentachlorobiphenyl and 2,2',4,5',6-pentachlorobiphenyl) to 24 and 8 mol%, respectively. However, the uterotonic activities of the latter congeners were modest when evaluated either solely or in a reconstituted mixture and could not fully account for the uterotonic activity of HR1260. Nonetheless, the relative abundance of congeners with three or fewer chlorines increased to 14 mol% as a group in HR1260, suggesting that these congeners collectively contribute to the uterotonic activity even though the abundance of any one congener in this group was less than 5 mol%.

Degradation of aroclor 1242 dechlorination products in sediments by Burkholderia xenovorans LB400(ohb) and Rhodococcus sp. strain RHA1(fcb).

Burkholderia xenovorans strain LB400, which possesses the biphenyl pathway, was engineered to contain the oxygenolytic ortho dehalogenation (ohb) operon, allowing it to grow on 2-chlorobenzoate and to completely mineralize 2-chlorobiphenyl. A two-stage anaerobic/aerobic biotreatment process for Aroclor 1242-contaminated sediment was simulated, and the degradation activities and genetic stabilities of LB400(ohb) and the previously constructed strain RHA1(fcb), capable of growth on 4-chlorobenzoate, were monitored during the aerobic phase. The population dynamics of both strains were also followed by selective plating and real-time PCR, with comparable results; populations of both recombinants increased in the contaminated sediment. Inoculation at different cell densities (10(4) or 10(6) cells g(-1) sediment) did not affect the extent of polychlorinated biphenyl (PCB) biodegradation. After 30 days, PCB removal rates for high and low inoculation densities were 57% and 54%, respectively, during the aerobic phase.

Long-term recovery of PCB-contaminated sediments at the Lake Hartwell superfund site: PCB dechlorination. 1. End-member characterization.

Under anaerobic conditions, such as those typically found in buried sediments, the primary metabolic pathway for polychlorinated biphenyls (PCBs) is reductive dechlorination in which chlorine removal and substitution with hydrogen by bacteria result in a reduced organic compound with fewer chlorines. Vertical sediment cores were collected from Lake Hartwell (Pickens County, SC) and analyzed in 5-cm intervals for 107 PCB congeners in a total of more than 280 samples from 18 sediment cores and surface samples. This paper reports on extensive PCB dechlorination measured in Lake Hartwell sediments and the characterization of dechlorination end-member (EM) patterns using chemical forensic methods. PCB congener fingerprinting and a multivariate receptor modeling method, polytopic vector analysis (PVA), were used for identification and characterization of weathered and dechlorinated PCB congener patterns. Dechlorination resulted in a substantial shift in buried sediments from tetra- through decachlorobiphenyl congeners to mono- through trichlorobiphenyl congeners. Mono- through trichlorobiphenyls comprised approximately 80% of the PCBs in buried sediments that underwent maximum dechlorination as compared to approximately 20% in surface sediments. The major concentration decreases were seen in the tetra- through hexachlorobiphenyl homologues, which accounted for over 90% of the dechlorination. Octa- through decachlorobiphenyl congeners also were dechlorinated, but their overall contribution to dechlorination was relatively small due to their low initial concentrations (< 5%). The net accumulation of 2-CB, 2,2'/2,6-DCBs, 2,4'-DCB, 2,2',4-TCB, and 2,2',6-TCB at Lake Hartwell matched characteristic PCB dechlorination products reported in the literature, such as those for Processes M, Q, and C; and the persistence of tetrachlorobiphenyls (TeCBs) that contained 24- and 25-congener groups resembled dechlorination Processes H or H'. Although dechlorination tended to be very extensive in most of the cores, it was not always consistent from core to core or at various depth intervals within a single core. The reason for this variability in dechlorination extent could not be determined from the existing data and did not appear to correlate with such factors as PCB concentration, total organic carbon, or age. The authors used fingerprinting analysis and a PVA multivariate receptor model as exploratory data analysis tools to characterize PCB sources and their alteration patterns. Dominant sources and alteration patterns were determined in this large data set by comparing PVA EM patterns with known source patterns (i.e., Aroclors or Aroclor mixtures) and literature-reported alteration patterns. PVA also afforded an opportunity to characterize the vertical and lateral distributions of the weathered and unweathered PCB source patterns and dechlorination patterns, a task that would have been much more difficult to accomplish through comparison of chromatograms alone.

Long-term recovery of PCB-contaminated sediments at the Lake Hartwell superfund site: PCB dechlorination. 2. Rates and extent.

This paper reports on extensive polychlorinated biphenyl (PCB) dechlorination measured in Lake Hartwell (Pickens County, SC) sediments. Vertical sediment cores were collected from 18 locations in Lake Hartwell (Pickens County, SC) and analyzed in 5-cm increments for PCB congeners. The preferential loss of meta and para chlorines with sediment depth demonstrated that PCBs in the sediments underwent reductive dechlorination after burial. Notably, ortho chlorines were highly conserved for more than 5 decades; since the first appearance of PCBs, ca. 1950-1955. These dechlorination characteristics resulted in the accumulation of lower chlorinated congeners dominated by ortho chlorine substituents. Dechlorination rates were determined by plotting the numbers of meta plus para chlorines per biphenyl molecule (mol of chlorine/mol of PCB) with sediment age. Regression analyses showed linear correlations between meta plus para chlorine concentrations with time. The average dechlorination rate was 0.094 +/- 0.063 mol of Cl/mol of PCB/yr. The rates measured using the 2001 cores were approximately twice those measured using the 2000 cores, most likely because the 2001 cores were collected only at transects O, L, and I, which had the highest rates measured in 2000. An inverse of the dechlorination rates indicated that 16.4 +/- 11.6 yr was required per meta plus para chlorine removal (ranging from 4.3 to 43.5 yr per chlorine removal). The rates determined from this study were 1-2 orders of magnitude lower than rates reported from laboratory microcosm studies using Hudson River and St. Lawrence River sediments, suggesting that dechlorination rates reported for laboratory experiments are much higher than those occurring in situ.

The inhibition of LPS-induced splenocyte proliferation by ortho-substituted and microbially dechlorinated polychlorinated biphenyls is associated with a decreased expression of cyclin D2.

Immunological effects of polychlorinated biphenyls (PCBs) have been demonstrated in our laboratories with the peferential inhibition of lipopolysaccharide (LPS)-induced splenocyte proliferation by ortho-substituted PCB congeners. An investigation of the mechanism behind this immunotoxicity revealed an interruption in the progression of murine lymphocytes from G0/G1 into S phase by Aroclor 1242 and the di-ortho-substituted congener, 2,2'-chlorobiphenyl (CB), whereas, a non-ortho-substituted congener, 4,4'-CB, did not affect cell cycle progression. This interruption of cell cycle progression by 2,2'-CB and Aroclor 1242 was associated with a decreased expression of the cell cycle regulatory protein, cyclin D2, while expression was not affected by exposure to the non-ortho-substituted 4,4'-CB. These results suggest the preferential inhibition of LPS-induced splenocyte proliferation by ortho-substituted congeners is a result of a decreased expression of cyclin D2, which leads to an interruption in cell cycle progression. In addition, PCB mixtures with an increased percentage of chlorines in the ortho position following an environmentally occurring degradation process inhibited LPS-induced proliferation, interrupted cell cycle progression, and decreased cyclin D2 expression. This study provides evidence for a mechanism of action of the immunological effects of ortho-substituted individual congeners as well as environmentally relevant mixtures enriched in congeners with this substitution pattern.

Classification of the biphenyl- and polychlorinated biphenyl-degrading strain LB400T and relatives as Burkholderia xenovorans sp. nov.

Strain LB400T is the best-studied polychlorinated biphenyl (PCB) degrader. This organism has previously been allocated in the genus Burkholderia, since its 16S rRNA gene sequence shows 98.6 % sequence similarity to the type strains of Burkholderia graminis and Burkholderia terricola. A polyphasic study was undertaken to clarify the actual taxonomic position of this biotechnologically important organism and of two strains, one recovered from a blood culture vial and one from a coffee plant rhizosphere, both of which resembled strain LB400T in their whole-cell protein patterns. DNA-DNA hybridization experiments revealed that the three strains represented a single novel species, for which the name Burkholderia xenovorans sp. nov. is proposed. Strains of this novel species can be differentiated phenotypically from nearly all other Burkholderia species by their inability to assimilate L-arabinose. The whole-cell fatty acid profile of B. xenovorans strains is consistent with their classification in the genus Burkholderia, with 18 : 1omega7c, 16 : 1omega7c, 16 : 0, 14 : 0 3OH, 16 : 0 3OH, 17 : 0 cyclo and 14 : 0 being the most abundant fatty acids. The G + C content of the species varies between 62.4 and 62.9 mol%. The type strain of B. xenovorans is LB400T (= LMG 21463T = CCUG 46959T = NRRL B-18064T).

Inhibition of LPS-induced splenocyte proliferation by ortho-substituted polychlorinated biphenyl congeners.

Polychlorinated biphenyls (PCBs) are persistent environmental contaminants, and their ubiquitous nature has prompted studies of their potential health hazards. As a result of their lipophilic nature, PCBs accumulate in breast milk and subsequently affect the health of offspring of exposed individuals. Biological effects of PCBs in animals have mostly been attributed to coplanar congeners, although effects of ortho congeners also have been demonstrated. To investigate the relationship of immunotoxicity and chlorine substitution pattern, the effects of PCB congeners and mixtures of ortho and non-ortho-substituted constituents of Aroclor 1242 on splenocytes from C57B1/6 mice were examined. The immunotoxic endpoints investigated included splenocyte viability, lipopolysaccharide (LPS)-induced splenocyte proliferation, and LPS-induced antibody secretion. Congeners with multiple ortho chlorines preferentially inhibited splenocyte proliferation as compared with non- or mono-ortho-substituted congeners. However, mixtures of non- and mono-ortho-substituted congeners and multi-ortho-substituted congeners inhibited LPS-induced splenocyte proliferation and antibody secretion at similar concentrations. Exposure of splenocytes to these mixtures did not activate the aryl hydrocarbon receptor (AhR) signal transduction pathway. These results suggest individual multi-ortho-substituted congeners preferentially inhibit LPS-induced splenocyte proliferation, while congeners not exhibiting an effect individually may have additive effects in a mixture to produce an immunotoxic response through an AhR-independent pathway.

Dechlorination of Four Commercial Polychlorinated Biphenyl Mixtures (Aroclors) by Anaerobic Microorganisms from Sediments.

The rate, extent, and pattern of dechlorination of four Aroclors by inocula prepared from two polychlorinated biphenyl (PCB)-contaminated sediments were compared. The four mixtures used, Aroclors 1242, 1248, 1254, and 1260, average approximately three, four, five, and six chlorines, respectively, per biphenyl molecule. All four Aroclors were dechlorinated with the loss of meta plus para chlorines ranging from 15 to 85%. Microorganisms from an Aroclor 1242-contaminated site in the upper Hudson River dechlorinated Aroclor 1242 to a greater extent than did microorganisms from Aroclor 1260-contaminated sediments from Silver Lake, Mass. The Silver Lake inoculum dechlorinated Aroclor 1260 more rapidly than the Hudson River inoculum did and showed a preferential removal of meta chlorines. For each inoculum the rate and extent of dechlorination tended to decrease as the degree of chlorination of the Aroclor increased, especially for Aroclor 1260. The maximal observed dechlorination rates were 0.3, 0.3, and 0.2 mug-atoms of Cl removed per g of sediment per week for Aroclors 1242, 1248, and 1254, respectively. The maximal observed dechlorination rates for Hudson River and Silver Lake organisms for Aroclor 1260 were 0.04 and 0.21 mug-atoms of Cl removed per g of sediment per week, respectively. The dechlorination patterns obtained suggested that the Hudson River microorganisms were more capable than the Silver Lake organisms of removing the last para chlorine. These results suggest that there are different PCB-dechlorinating microorganisms at different sites, with characteristic specificities for PCB dechlorination.