Clades of huge phages from across Earth's ecosystems.

Bacteriophages typically have small genomes1 and depend on their bacterial hosts for replication2. Here we sequenced DNA from diverse ecosystems and found hundreds of phage genomes with lengths of more than 200 kilobases (kb), including a genome of 735 kb, which is-to our knowledge-the largest phage genome to be described to date. Thirty-five genomes were manually curated to completion (circular and no gaps). Expanded genetic repertoires include diverse and previously undescribed CRISPR-Cas systems, transfer RNAs (tRNAs), tRNA synthetases, tRNA-modification enzymes, translation-initiation and elongation factors, and ribosomal proteins. The CRISPR-Cas systems of phages have the capacity to silence host transcription factors and translational genes, potentially as part of a larger interaction network that intercepts translation to redirect biosynthesis to phage-encoded functions. In addition, some phages may repurpose bacterial CRISPR-Cas systems to eliminate competing phages. We phylogenetically define the major clades of huge phages from human and other animal microbiomes, as well as from oceans, lakes, sediments, soils and the built environment. We conclude that the large gene inventories of huge phages reflect a conserved biological strategy, and that the phages are distributed across a broad bacterial host range and across Earth's ecosystems.

Poisson hurdle model-based method for clustering microbiome features.

MOTIVATION: High-throughput sequencing technologies have greatly facilitated microbiome research and have generated a large volume of microbiome data with the potential to answer key questions regarding microbiome assembly, structure and function. Cluster analysis aims to group features that behave similarly across treatments, and such grouping helps to highlight the functional relationships among features and may provide biological insights into microbiome networks. However, clustering microbiome data are challenging due to the sparsity and high dimensionality. RESULTS: We propose a model-based clustering method based on Poisson hurdle models for sparse microbiome count data. We describe an expectation-maximization algorithm and a modified version using simulated annealing to conduct the cluster analysis. Moreover, we provide algorithms for initialization and choosing the number of clusters. Simulation results demonstrate that our proposed methods provide better clustering results than alternative methods under a variety of settings. We also apply the proposed method to a sorghum rhizosphere microbiome dataset that results in interesting biological findings. AVAILABILITY AND IMPLEMENTATION: R package is freely available for download at https://cran.r-project.org/package=PHclust. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Conservation of beneficial microbes between the rhizosphere and the cyanosphere.

Biocrusts are phototroph-driven communities inhabiting arid soil surfaces. Like plants, most photoautotrophs (largely cyanobacteria) in biocrusts are thought to exchange fixed carbon for essential nutrients like nitrogen with cyanosphere bacteria. Here, we aim to compare beneficial interactions in rhizosphere and cyanosphere environments, including finding growth-promoting strains for hosts from both environments. To examine this, we performed a retrospective analysis of 16S rRNA gene sequencing datasets, host-microbe co-culture experiments between biocrust communities/biocrust isolates and a model grass (Brachypodium distachyon) or a dominant biocrust cyanobacterium (Microcoleus vaginatus), and metabolomic analysis. All 18 microbial phyla in the cyanosphere were also present in the rhizosphere, with additional 17 phyla uniquely found in the rhizosphere. The biocrust microbes promoted the growth of the model grass, and three biocrust isolates (Bosea sp._L1B56, Pseudarthrobacter sp._L1D14 and Pseudarthrobacter picheli_L1D33) significantly promoted the growth of both hosts. Moreover, pantothenic acid was produced by Pseudarthrobacter sp._L1D14 when grown on B. distachyon exudates, and supplementation of plant growth medium with this metabolite increased B. distachyon biomass by over 60%. These findings suggest that cyanobacteria and other diverse photoautotrophic hosts can be a source for new plant growth-promoting microbes and metabolites.

Exploring the roles of microbes in facilitating plant adaptation to climate change.

Plants benefit from their close association with soil microbes which assist in their response to abiotic and biotic stressors. Yet much of what we know about plant stress responses is based on studies where the microbial partners were uncontrolled and unknown. Under climate change, the soil microbial community will also be sensitive to and respond to abiotic and biotic stressors. Thus, facilitating plant adaptation to climate change will require a systems-based approach that accounts for the multi-dimensional nature of plant-microbe-environment interactions. In this perspective, we highlight some of the key factors influencing plant-microbe interactions under stress as well as new tools to facilitate the controlled study of their molecular complexity, such as fabricated ecosystems and synthetic communities. When paired with genomic and biochemical methods, these tools provide researchers with more precision, reproducibility, and manipulability for exploring plant-microbe-environment interactions under a changing climate.

Division of labor in honey bee gut microbiota for plant polysaccharide digestion.

Bees acquire carbohydrates from nectar and lipids; and amino acids from pollen, which also contains polysaccharides including cellulose, hemicellulose, and pectin. These potential energy sources could be degraded and fermented through microbial enzymatic activity, resulting in short chain fatty acids available to hosts. However, the contributions of individual microbiota members to polysaccharide digestion have remained unclear. Through analysis of bacterial isolate genomes and a metagenome of the honey bee gut microbiota, we identify that Bifidobacterium and Gilliamella are the principal degraders of hemicellulose and pectin. Both Bifidobacterium and Gilliamella show extensive strain-level diversity in gene repertoires linked to polysaccharide digestion. Strains from honey bees possess more such genes than strains from bumble bees. In Bifidobacterium, genes encoding carbohydrate-active enzymes are colocated within loci devoted to polysaccharide utilization, as in Bacteroides from the human gut. Carbohydrate-active enzyme-encoding gene expressions are up-regulated in response to particular hemicelluloses both in vitro and in vivo. Metabolomic analyses document that bees experimentally colonized by different strains generate distinctive gut metabolomic profiles, with enrichment for specific monosaccharides, corresponding to predictions from genomic data. The other 3 core gut species clusters (Snodgrassella and 2 Lactobacillus clusters) possess few or no genes for polysaccharide digestion. Together, these findings indicate that strain composition within individual hosts determines the metabolic capabilities and potentially affects host nutrition. Furthermore, the niche specialization revealed by our study may promote overall community stability in the gut microbiomes of bees.

Feature selection and causal analysis for microbiome studies in the presence of confounding using standardization.

BACKGROUND: Microbiome studies have uncovered associations between microbes and human, animal, and plant health outcomes. This has led to an interest in developing microbial interventions for treatment of disease and optimization of crop yields which requires identification of microbiome features that impact the outcome in the population of interest. That task is challenging because of the high dimensionality of microbiome data and the confounding that results from the complex and dynamic interactions among host, environment, and microbiome. In the presence of such confounding, variable selection and estimation procedures may have unsatisfactory performance in identifying microbial features with an effect on the outcome. RESULTS: In this manuscript, we aim to estimate population-level effects of individual microbiome features while controlling for confounding by a categorical variable. Due to the high dimensionality and confounding-induced correlation between features, we propose feature screening, selection, and estimation conditional on each stratum of the confounder followed by a standardization approach to estimation of population-level effects of individual features. Comprehensive simulation studies demonstrate the advantages of our approach in recovering relevant features. Utilizing a potential-outcomes framework, we outline assumptions required to ascribe causal, rather than associational, interpretations to the identified microbiome effects. We conducted an agricultural study of the rhizosphere microbiome of sorghum in which nitrogen fertilizer application is a confounding variable. In this study, the proposed approach identified microbial taxa that are consistent with biological understanding of potential plant-microbe interactions. CONCLUSIONS: Standardization enables more accurate identification of individual microbiome features with an effect on the outcome of interest compared to other variable selection and estimation procedures when there is confounding by a categorical variable.

Peatland microbial community responses to plant functional group and drought are depth-dependent.

Peatlands store one-third of Earth's soil carbon, the stability of which is uncertain due to climate change-driven shifts in hydrology and vegetation, and consequent impacts on microbial communities that mediate decomposition. Peatland carbon cycling varies over steep physicochemical gradients characterizing vertical peat profiles. However, it is unclear how drought-mediated changes in plant functional groups (PFGs) and water table (WT) levels affect microbial communities at different depths. We combined a multiyear mesocosm experiment with community sequencing across a 70-cm depth gradient, to test the hypotheses that vascular PFGs (Ericaceae vs. sedges) and WT (high vs. low) structure peatland microbial communities in depth-dependent ways. Several key results emerged. (i) Both fungal and prokaryote (bacteria and archaea) community structure shifted with WT and PFG manipulation, but fungi were much more sensitive to PFG whereas prokaryotes were much more sensitive to WT. (ii) PFG effects were largely driven by Ericaceae, although sedge effects were evident in specific cases (e.g., methanotrophs). (iii) Treatment effects varied with depth: the influence of PFG was strongest in shallow peat (0-10, 10-20 cm), whereas WT effects were strongest at the surface and middle depths (0-10, 30-40 cm), and all treatment effects waned in the deepest peat (60-70 cm). Our results underline the depth-dependent and taxon-specific ways that plant communities and hydrologic variability shape peatland microbial communities, pointing to the importance of understanding how these factors integrate across soil profiles when examining peatland responses to climate change.

Large-scale replicated field study of maize rhizosphere identifies heritable microbes.

Soil microbes that colonize plant roots and are responsive to differences in plant genotype remain to be ascertained for agronomically important crops. From a very large-scale longitudinal field study of 27 maize inbred lines planted in three fields, with partial replication 5 y later, we identify root-associated microbiota exhibiting reproducible associations with plant genotype. Analysis of 4,866 samples identified 143 operational taxonomic units (OTUs) whose variation in relative abundances across the samples was significantly regulated by plant genotype, and included five of seven core OTUs present in all samples. Plant genetic effects were significant amid the large effects of plant age on the rhizosphere microbiome, regardless of the specific community of each field, and despite microbiome responses to climate events. Seasonal patterns showed that the plant root microbiome is locally seeded, changes with plant growth, and responds to weather events. However, against this background of variation, specific taxa responded to differences in host genotype. If shown to have beneficial functions, microbes may be considered candidate traits for selective breeding.

Restoring wetlands on intensive agricultural lands modifies nitrogen cycling microbial communities and reduces N2O production potential.

The concentration of nitrous oxide (N2O), an ozone-depleting greenhouse gas, is rapidly increasing in the atmosphere. Most atmospheric N2O originates in terrestrial ecosystems, of which the majority can be attributed to microbial cycling of nitrogen in agricultural soils. Here, we demonstrate how the abundance of nitrogen cycling genes vary across intensively managed agricultural fields and adjacent restored wetlands in the Sacramento-San Joaquin Delta in California, USA. We found that the abundances of nirS and nirK genes were highest at the intensively managed organic-rich cornfield and significantly outnumber any other gene abundances, suggesting very high N2O production potential. The quantity of nitrogen transforming genes, particularly those responsible for denitrification, nitrification and DNRA, were highest in the agricultural sites, whereas nitrogen fixation and ANAMMOX was strongly associated with the wetland sites. Although the abundance of nosZ genes was also high at the agricultural sites, the ratio of nosZ genes to nir genes was significantly higher in wetland sites indicating that these sites could act as a sink of N2O. These findings suggest that wetland restoration could be a promising natural climate solution not only for carbon sequestration but also for reduced N2O emissions.

Discovery of enzymes for toluene synthesis from anoxic microbial communities.

Microbial toluene biosynthesis was reported in anoxic lake sediments more than three decades ago, but the enzyme catalyzing this biochemically challenging reaction has never been identified. Here we report the toluene-producing enzyme PhdB, a glycyl radical enzyme of bacterial origin that catalyzes phenylacetate decarboxylation, and its cognate activating enzyme PhdA, a radical S-adenosylmethionine enzyme, discovered in two distinct anoxic microbial communities that produce toluene. The unconventional process of enzyme discovery from a complex microbial community (>300,000 genes), rather than from a microbial isolate, involved metagenomics- and metaproteomics-enabled biochemistry, as well as in vitro confirmation of activity with recombinant enzymes. This work expands the known catalytic range of glycyl radical enzymes (only seven reaction types had been characterized previously) and aromatic-hydrocarbon-producing enzymes, and will enable first-time biochemical synthesis of an aromatic fuel hydrocarbon from renewable resources, such as lignocellulosic biomass, rather than from petroleum.

Interactions between plants and soil shaping the root microbiome under abiotic stress.

Plants growing in soil develop close associations with soil microorganisms, which inhabit the areas around, on, and inside their roots. These microbial communities and their associated genes - collectively termed the root microbiome - are diverse and have been shown to play an important role in conferring abiotic stress tolerance to their plant hosts. In light of growing concerns over the threat of water and nutrient stress facing terrestrial ecosystems, especially those used for agricultural production, increased emphasis has been placed on understanding how abiotic stress conditions influence the composition and functioning of the root microbiome and the ultimate consequences for plant health. However, the composition of the root microbiome under abiotic stress conditions will not only reflect shifts in the greater bulk soil microbial community from which plants recruit their root microbiome but also plant responses to abiotic stress, which include changes in root exudate profiles and morphology. Exploring the relative contributions of these direct and plant-mediated effects on the root microbiome has been the focus of many studies in recent years. Here, we review the impacts of abiotic stress affecting terrestrial ecosystems, specifically flooding, drought, and changes in nitrogen and phosphorus availability, on bulk soil microbial communities and plants that interact to ultimately shape the root microbiome. We conclude with a perspective outlining possible directions for future research needed to advance our understanding of the complex molecular and biochemical interactions between soil, plants, and microbes that ultimately determine the composition of the root microbiome under abiotic stress.

Insight into the Bacterial Endophytic Communities of Peach Cultivars Related to Crown Gall Disease Resistance.

Crown gall disease caused by Agrobacterium tumefaciens severely impacts the production of peach and other fruit trees. Several peach cultivars are partially resistant to A. tumefaciens, but little is known about the roles of endophytic microbiota in disease resistance. In the present study, the endophytic bacterial communities of resistant and susceptible peach cultivars "Honggengansutao" and "Okinawa" were analyzed using universal 16S rRNA gene amplicon sequencing in parallel with the cultivation and characterization of bacterial isolates. A total of 1,357,088 high-quality sequences representing 3,160 distinct operational taxonomic units (OTUs; Proteobacteria, Actinobacteria, Bacteroidetes, and Firmicutes) and 1,200 isolates of 20 genera and 305 distinct ribotypes were collected from peach roots and twigs. It was found that factors including plant developmental stage, cultivar, and A. tumefaciens invasion strongly influenced the peach endophytic communities. The community diversity of endophytic bacteria and the abundance of culturable bacteria were both higher in the roots of the resistant cultivar, particularly after inoculation. Strikingly, the pathogen antagonists Streptomyces and Pseudomonas in roots and Rhizobium in twigs were most frequently detected in resistant plants. Our results suggest that the higher abundance and diversity of endophytic bacteria and increased proportions of antagonistic bacteria might contribute to the natural defense of the resistant cultivar against A. tumefaciens This work reveals the relationships between endophytic bacteria and disease resistance in peach plants and provides important information for microbiome-based biocontrol of crown gall disease in fruit trees.IMPORTANCEAgrobacterium tumefaciens as the causal agent of peach crown gall disease can be controlled by planting resistant cultivars. This study profiles the endophytic bacteria in susceptible and resistant peach cultivars, advancing our understanding of the relationships between endophytic bacterial communities and peach crown gall disease, with potential implications for other complex microbiome-plant-pathogen interactions. The resistant cultivar may defend itself by increasing the diversity and abundance of beneficial endophytic bacteria. The antagonists identified among the genera Streptomyces, Pseudomonas, and Rhizobium may have application potential for biocontrol of crown gall disease in fruit trees.

Accurate, multi-kb reads resolve complex populations and detect rare microorganisms.

Accurate evaluation of microbial communities is essential for understanding global biogeochemical processes and can guide bioremediation and medical treatments. Metagenomics is most commonly used to analyze microbial diversity and metabolic potential, but assemblies of the short reads generated by current sequencing platforms may fail to recover heterogeneous strain populations and rare organisms. Here we used short (150-bp) and long (multi-kb) synthetic reads to evaluate strain heterogeneity and study microorganisms at low abundance in complex microbial communities from terrestrial sediments. The long-read data revealed multiple (probably dozens of) closely related species and strains from previously undescribed Deltaproteobacteria and Aminicenantes (candidate phylum OP8). Notably, these are the most abundant organisms in the communities, yet short-read assemblies achieved only partial genome coverage, mostly in the form of short scaffolds (N50 = ∼ 2200 bp). Genome architecture and metabolic potential for these lineages were reconstructed using a new synteny-based method. Analysis of long-read data also revealed thousands of species whose abundances were <0.1% in all samples. Most of the organisms in this "long tail" of rare organisms belong to phyla that are also represented by abundant organisms. Genes encoding glycosyl hydrolases are significantly more abundant than expected in rare genomes, suggesting that rare species may augment the capability for carbon turnover and confer resilience to changing environmental conditions. Overall, the study showed that a diversity of closely related strains and rare organisms account for a major portion of the communities. These are probably common features of many microbial communities and can be effectively studied using a combination of long and short reads.

Defining the core Arabidopsis thaliana root microbiome.

Land plants associate with a root microbiota distinct from the complex microbial community present in surrounding soil. The microbiota colonizing the rhizosphere (immediately surrounding the root) and the endophytic compartment (within the root) contribute to plant growth, productivity, carbon sequestration and phytoremediation. Colonization of the root occurs despite a sophisticated plant immune system, suggesting finely tuned discrimination of mutualists and commensals from pathogens. Genetic principles governing the derivation of host-specific endophyte communities from soil communities are poorly understood. Here we report the pyrosequencing of the bacterial 16S ribosomal RNA gene of more than 600 Arabidopsis thaliana plants to test the hypotheses that the root rhizosphere and endophytic compartment microbiota of plants grown under controlled conditions in natural soils are sufficiently dependent on the host to remain consistent across different soil types and developmental stages, and sufficiently dependent on host genotype to vary between inbred Arabidopsis accessions. We describe different bacterial communities in two geochemically distinct bulk soils and in rhizosphere and endophytic compartments prepared from roots grown in these soils. The communities in each compartment are strongly influenced by soil type. Endophytic compartments from both soils feature overlapping, low-complexity communities that are markedly enriched in Actinobacteria and specific families from other phyla, notably Proteobacteria. Some bacteria vary quantitatively between plants of different developmental stage and genotype. Our rigorous definition of an endophytic compartment microbiome should facilitate controlled dissection of plant-microbe interactions derived from complex soil communities.

The Tale of a Neglected Energy Source: Elevated Hydrogen Exposure Affects both Microbial Diversity and Function in Soil.

The enrichment of H2-oxidizing bacteria (HOB) by H2 generated by nitrogen-fixing nodules has been shown to have a fertilization effect on several different crops. The benefit of HOB is attributed to their production of plant growth-promoting factors, yet their interactions with other members of soil microbial communities have received little attention. Here we report that the energy potential of H2, when supplied to soil, alters ecological niche partitioning of bacteria and fungi, with multifaceted consequences for both generalist and specialist microbial functions. We used dynamic microcosms to expose soil to the typical atmospheric H2 mixing ratio (0.5 ppmv) permeating soils, as well as mixing ratios comparable to those found at the soil-nodule interface (10,000 ppmv). Elevated H2 exposure exerted direct effects on two HOB subpopulations distinguished by their affinity for H2 while enhancing community level carbon substrate utilization potential and lowering CH4 uptake activity in soil. We found that H2 triggered changes in the abundance of microorganisms that were reproducible yet inconsistent across soils at the taxonomic level and even among HOB. Overall, H2 exposure altered microbial process rates at an intensity that depends upon soil abiotic and biotic features. We argue that further examination of direct and indirect effects of H2 on soil microbial communities will lead to a better understanding of the H2 fertilization effect and soil biogeochemical processes.IMPORTANCE An innovative dynamic microcosm chamber system was used to demonstrate that H2 diffusing in soil triggers changes in the distribution of HOB and non-HOB. Although the response was uneven at the taxonomic level, an unexpected coordinated response of microbial functions was observed, including abatement of CH4 oxidation activity and stimulation of carbon turnover. Our work suggests that elevated H2 rewires soil biogeochemical structure through a combination of direct effects on the growth and persistence of HOB and indirect effects on a variety of microbial processes involving HOB and non-HOB.

Metagenomic and Metatranscriptomic Analyses Reveal the Structure and Dynamics of a Dechlorinating Community Containing Dehalococcoides mccartyi and Corrinoid-Providing Microorganisms under Cobalamin-Limited Conditions.

The aim of this study is to obtain a systems-level understanding of the interactions between Dehalococcoides and corrinoid-supplying microorganisms by analyzing community structures and functional compositions, activities, and dynamics in trichloroethene (TCE)-dechlorinating enrichments. Metagenomes and metatranscriptomes of the dechlorinating enrichments with and without exogenous cobalamin were compared. Seven putative draft genomes were binned from the metagenomes. At an early stage (2 days), more transcripts of genes in the Veillonellaceae bin-genome were detected in the metatranscriptome of the enrichment without exogenous cobalamin than in the one with the addition of cobalamin. Among these genes, sporulation-related genes exhibited the highest differential expression when cobalamin was not added, suggesting a possible release route of corrinoids from corrinoid producers. Other differentially expressed genes include those involved in energy conservation and nutrient transport (including cobalt transport). The most highly expressed corrinoid de novo biosynthesis pathway was also assigned to the Veillonellaceae bin-genome. Targeted quantitative PCR (qPCR) analyses confirmed higher transcript abundances of those corrinoid biosynthesis genes in the enrichment without exogenous cobalamin than in the enrichment with cobalamin. Furthermore, the corrinoid salvaging and modification pathway of Dehalococcoides was upregulated in response to the cobalamin stress. This study provides important insights into the microbial interactions and roles played by members of dechlorinating communities under cobalamin-limited conditions.IMPORTANCE The key chloroethene-dechlorinating bacterium Dehalococcoides mccartyi is a cobalamin auxotroph, thus acquiring corrinoids from other community members. Therefore, it is important to investigate the microbe-microbe interactions between Dehalococcoides and the corrinoid-providing microorganisms in a community. This study provides systems-level information, i.e., taxonomic and functional compositions and dynamics of the supportive microorganisms in dechlorinating communities under different cobalamin conditions. The findings shed light on the important roles of Veillonellaceae species in the communities compared to other coexisting community members in producing and providing corrinoids for Dehalococcoides species under cobalamin-limited conditions.

Conversion of Amazon rainforest to agriculture alters community traits of methane-cycling organisms.

Land use change is one of the greatest environmental impacts worldwide, especially to tropical forests. The Amazon rainforest has been subject to particularly high rates of land use change, primarily to cattle pasture. A commonly observed response to cattle pasture establishment in the Amazon is the conversion of soil from a methane sink in rainforest, to a methane source in pasture. However, it is not known how the microorganisms that mediate methane flux are altered by land use change. Here, we use the deepest metagenomic sequencing of Amazonian soil to date to investigate differences in methane-cycling microorganisms and their traits across rainforest and cattle pasture soils. We found that methane-cycling microorganisms responded to land use change, with the strongest responses exhibited by methane-consuming, rather than methane-producing, microorganisms. These responses included a reduction in the relative abundance of methanotrophs and a significant decrease in the abundance of genes encoding particulate methane monooxygenase. We also observed compositional changes to methanotroph and methanogen communities as well as changes to methanotroph life history strategies. Our observations suggest that methane-cycling microorganisms are vulnerable to land use change, and this vulnerability may underlie the response of methane flux to land use change in Amazon soils.

Genome-Resolved Meta-Omics Ties Microbial Dynamics to Process Performance in Biotechnology for Thiocyanate Degradation.

Remediation of industrial wastewater is important for preventing environmental contamination and enabling water reuse. Biological treatment for one industrial contaminant, thiocyanate (SCN-), relies upon microbial hydrolysis, but this process is sensitive to high loadings. To examine the activity and stability of a microbial community over increasing SCN- loadings, we established and operated a continuous-flow bioreactor fed increasing loadings of SCN-. A second reactor was fed ammonium sulfate to mimic breakdown products of SCN-. Biomass was sampled from both reactors for metagenomics and metaproteomics, yielding a set of genomes for 144 bacteria and one rotifer that constituted the abundant community in both reactors. We analyzed the metabolic potential and temporal dynamics of these organisms across the increasing loadings. In the SCN- reactor, Thiobacillus strains capable of SCN- degradation were highly abundant, whereas the ammonium sulfate reactor contained nitrifiers and heterotrophs capable of nitrate reduction. Key organisms in the SCN- reactor expressed proteins involved in SCN- degradation, sulfur oxidation, carbon fixation, and nitrogen removal. Lower performance at higher loadings was linked to changes in microbial community composition. This work provides an example of how meta-omics can increase our understanding of industrial wastewater treatment and inform iterative process design and development.

Tackling soil diversity with the assembly of large, complex metagenomes.

The large volumes of sequencing data required to sample deeply the microbial communities of complex environments pose new challenges to sequence analysis. De novo metagenomic assembly effectively reduces the total amount of data to be analyzed but requires substantial computational resources. We combine two preassembly filtering approaches--digital normalization and partitioning--to generate previously intractable large metagenome assemblies. Using a human-gut mock community dataset, we demonstrate that these methods result in assemblies nearly identical to assemblies from unprocessed data. We then assemble two large soil metagenomes totaling 398 billion bp (equivalent to 88,000 Escherichia coli genomes) from matched Iowa corn and native prairie soils. The resulting assembled contigs could be used to identify molecular interactions and reaction networks of known metabolic pathways using the Kyoto Encyclopedia of Genes and Genomes Orthology database. Nonetheless, more than 60% of predicted proteins in assemblies could not be annotated against known databases. Many of these unknown proteins were abundant in both corn and prairie soils, highlighting the benefits of assembly for the discovery and characterization of novelty in soil biodiversity. Moreover, 80% of the sequencing data could not be assembled because of low coverage, suggesting that considerably more sequencing data are needed to characterize the functional content of soil.

Gill bacteria enable a novel digestive strategy in a wood-feeding mollusk.

Bacteria play many important roles in animal digestive systems, including the provision of enzymes critical to digestion. Typically, complex communities of bacteria reside in the gut lumen in direct contact with the ingested materials they help to digest. Here, we demonstrate a previously undescribed digestive strategy in the wood-eating marine bivalve Bankia setacea, wherein digestive bacteria are housed in a location remote from the gut. These bivalves, commonly known as shipworms, lack a resident microbiota in the gut compartment where wood is digested but harbor endosymbiotic bacteria within specialized cells in their gills. We show that this comparatively simple bacterial community produces wood-degrading enzymes that are selectively translocated from gill to gut. These enzymes, which include just a small subset of the predicted wood-degrading enzymes encoded in the endosymbiont genomes, accumulate in the gut to the near exclusion of other endosymbiont-made proteins. This strategy of remote enzyme production provides the shipworm with a mechanism to capture liberated sugars from wood without competition from an endogenous gut microbiota. Because only those proteins required for wood digestion are translocated to the gut, this newly described system reveals which of many possible enzymes and enzyme combinations are minimally required for wood degradation. Thus, although it has historically had negative impacts on human welfare, the shipworm digestive process now has the potential to have a positive impact on industries that convert wood and other plant biomass to renewable fuels, fine chemicals, food, feeds, textiles, and paper products.