Search

Soil microbiomes are highly diverse, and to improve their representation in biogeochemical models, microbial genome data can be leveraged to infer key functional traits. By integrating genome-inferred traits into a theory-based hierarchical framework, emergent behaviour arising from interactions of individual traits can be predicted. Here we combine theory-driven predictions of substrate uptake kinetics with a genome-informed trait-based dynamic energy budget model to predict emergent life-history traits and trade-offs in soil bacteria. When applied to a plant microbiome system, the model accurately predicted distinct substrate-acquisition strategies that aligned with observations, uncovering resource-dependent trade-offs between microbial growth rate and efficiency. For instance, inherently slower-growing microorganisms, favoured by organic acid exudation at later plant growth stages, exhibited enhanced carbon use efficiency (yield) without sacrificing growth rate (power). This insight has implications for retaining plant root-derived carbon in soils and highlights the power of data-driven, trait-based approaches for improving microbial representation in biogeochemical models.

BackgroundViruses impact nearly all organisms on Earth, including microbial communities and their associated biogeochemical processes. In soils, highly diverse viral communities have been identified, with a global distribution seemingly driven by multiple biotic and abiotic factors, especially soil temperature and moisture. However, our current understanding of the stability of soil viral communities across time and their response to strong seasonal changes in environmental parameters remains limited. Here, we investigated the diversity and activity of environmental soil DNA and RNA viruses, focusing especially on bacteriophages, across dynamics' seasonal changes in a snow-dominated mountainous watershed by examining paired metagenomes and metatranscriptomes.ResultsWe identified a large number of DNA and RNA viruses taxonomically divergent from existing environmental viruses, including a significant proportion of fungal RNA viruses, and a large and unsuspected diversity of positive single-stranded RNA phages (Leviviricetes), highlighting the under-characterization of the global soil virosphere. Among these, we were able to distinguish subsets of active DNA and RNA phages that changed across seasons, consistent with a "seed-bank" viral community structure in which new phage activity, for example, replication and host lysis, is sequentially triggered by changes in environmental conditions. At the population level, we further identified virus-host dynamics matching two existing ecological models: "Kill-The-Winner" which proposes that lytic phages are actively infecting abundant bacteria, and "Piggyback-The-Persistent" which argues that when the host is growing slowly, it is more beneficial to remain in a dormant state. The former was associated with summer months of high and rapid microbial activity, and the latter with winter months of limited and slow host growth.ConclusionTaken together, these results suggest that the high diversity of viruses in soils is likely associated with a broad range of host interaction types each adapted to specific host ecological strategies and environmental conditions. As our understanding of how environmental and host factors drive viral activity in soil ecosystems progresses, integrating these viral impacts in complex natural microbiome models will be key to accurately predict ecosystem biogeochemistry. Video Abstract.

BackgroundAdvances in microbiome science are being driven in large part due to our ability to study and infer microbial ecology from genomes reconstructed from mixed microbial communities using metagenomics and single-cell genomics. Such omics-based techniques allow us to read genomic blueprints of microorganisms, decipher their functional capacities and activities, and reconstruct their roles in biogeochemical processes. Currently available tools for analyses of genomic data can annotate and depict metabolic functions to some extent; however, no standardized approaches are currently available for the comprehensive characterization of metabolic predictions, metabolite exchanges, microbial interactions, and microbial contributions to biogeochemical cycling.ResultsWe present METABOLIC (METabolic And BiogeOchemistry anaLyses In miCrobes), a scalable software to advance microbial ecology and biogeochemistry studies using genomes at the resolution of individual organisms and/or microbial communities. The genome-scale workflow includes annotation of microbial genomes, motif validation of biochemically validated conserved protein residues, metabolic pathway analyses, and calculation of contributions to individual biogeochemical transformations and cycles. The community-scale workflow supplements genome-scale analyses with determination of genome abundance in the microbiome, potential microbial metabolic handoffs and metabolite exchange, reconstruction of functional networks, and determination of microbial contributions to biogeochemical cycles. METABOLIC can take input genomes from isolates, metagenome-assembled genomes, or single-cell genomes. Results are presented in the form of tables for metabolism and a variety of visualizations including biogeochemical cycling potential, representation of sequential metabolic transformations, community-scale microbial functional networks using a newly defined metric "MW-score" (metabolic weight score), and metabolic Sankey diagrams. METABOLIC takes ~ 3 h with 40 CPU threads to process ~ 100 genomes and corresponding metagenomic reads within which the most compute-demanding part of hmmsearch takes ~ 45 min, while it takes ~ 5 h to complete hmmsearch for ~ 3600 genomes. Tests of accuracy, robustness, and consistency suggest METABOLIC provides better performance compared to other software and online servers. To highlight the utility and versatility of METABOLIC, we demonstrate its capabilities on diverse metagenomic datasets from the marine subsurface, terrestrial subsurface, meadow soil, deep sea, freshwater lakes, wastewater, and the human gut.ConclusionMETABOLIC enables the consistent and reproducible study of microbial community ecology and biogeochemistry using a foundation of genome-informed microbial metabolism, and will advance the integration of uncultivated organisms into metabolic and biogeochemical models. METABOLIC is written in Perl and R and is freely available under GPLv3 at https://github.com/AnantharamanLab/METABOLIC . Video abstract.

Remote sensing approaches have revolutionized the study of macroorganisms, allowing theories of population and community ecology to be tested across increasingly larger scales without much compromise in resolution of biological complexity. In microbial ecology, our remote window into the ecology of microorganisms is through the lens of genome sequencing. For microbial organisms, recent evidence from genomes recovered from metagenomic samples corroborate a highly complex view of their metabolic diversity and other associated traits which map into high physiological complexity. Regardless, during the first decades of this omics era, microbial ecological research has primarily focused on taxa and functional genes as ecological units, favoring breadth of coverage over resolution of biological complexity manifested as physiological diversity. Recently, the rate at which provisional draft genomes are generated has increased substantially, giving new insights into ecological processes and interactions. From a genotype perspective, the wide availability of genome-centric data requires new data synthesis approaches that place organismal genomes center stage in the study of environmental roles and functional performance. Extraction of ecologically relevant traits from microbial genomes will be essential to the future of microbial ecological research. Here, we present microTrait, a computational pipeline that infers and distills ecologically relevant traits from microbial genome sequences. microTrait maps a genome sequence into a trait space, including discrete and continuous traits, as well as simple and composite. Traits are inferred from genes and pathways representing energetic, resource acquisition, and stress tolerance mechanisms, while genome-wide signatures are used to infer composite, or life history, traits of microorganisms. This approach is extensible to any microbial habitat, although we provide initial examples of this approach with reference to soil microbiomes.

Carbohydrate rich substrates such as lignocellulosic hydrolysates remain one of the primary sources of potentially renewable fuel and bulk chemicals. The pentose sugar D-xylose is often present in significant amounts along with hexoses. Saccharomyces cerevisiae can acquire the ability to metabolize D-xylose through expression of heterologous D-xylose isomerase (XI). This enzyme is notoriously difficult to express in S. cerevisiae and only fourteen XIs have been reported to be active so far. We cloned a new D-xylose isomerase derived from microorganisms in the gut of the wood-feeding beetle Odontotaenius disjunctus. Although somewhat homologous to the XI from Piromyces sp. E2, the new gene was identified as bacterial in origin and the host as a Parabacteroides sp. Expression of the new XI in S. cerevisiae resulted in faster aerobic growth than the XI from Piromyces on D-xylose media. The D-xylose isomerization rate conferred by the new XI was also 72% higher, while absolute xylitol production was identical in both strains. Interestingly, increasing concentrations of xylitol (up to 8 g L-1) appeared not to inhibit D-xylose consumption. The newly described XI displayed 2.6 times higher specific activity, 37% lower KM for D-xylose, and exhibited higher activity over a broader temperature range, retaining 51% of maximal activity at 30 °C compared with only 29% activity for the Piromyces XI.

The interaction between gut bacterial and viral microbiota is thought to be important in human health. While fluctuations in female genital tract (FGT) bacterial microbiota similarly determine sexual health, little is known about the presence, persistence, and function of vaginal bacteriophages. We conducted shotgun metagenome sequencing of cervicovaginal samples from South African adolescents collected longitudinally, who received no antibiotics. We annotated viral reads and circular bacteriophages, identified CRISPR loci and putative prophages, and assessed their diversity, persistence, and associations with bacterial microbiota composition. Siphoviridae was the most prevalent bacteriophage family, followed by Myoviridae, Podoviridae, Herelleviridae, and Inoviridae. Full-length siphoviruses targeting bacterial vaginosis (BV)-associated bacteria were identified, suggesting their presence in vivo. CRISPR loci and prophage-like elements were common, and genomic analysis suggested higher diversity among Gardnerella than Lactobacillus prophages. We found that some prophages were highly persistent within participants, and identical prophages were present in cervicovaginal secretions of multiple participants, suggesting that prophages, and thus bacterial strains, are shared between adolescents. The number of CRISPR loci and prophages were associated with vaginal microbiota stability and absence of BV. Our analysis suggests that (pro)phages are common in the FGT and vaginal bacteria and (pro)phages may interact.

Microorganisms have evolved several mechanisms to mobilize and mineralize occluded and insoluble phosphorus (P), thereby promoting plant growth in terrestrial ecosystems. However, the linkages between microbial P-solubilization traits and the preponderance of insoluble P in natural ecosystems are not well known. We tested the P solubilization traits of hundreds of culturable bacteria representative of the rhizosphere from a natural gradient where P concentration and bioavailability decline as soil becomes progressively more weathered. Aluminum, iron phosphate and organic P (phytate) were expected to dominate in more weathered soils. A defined cultivation medium with these chemical forms of P was used for isolation. A combination of soil chemical, spectroscopic analyses and 16S rRNA gene sequencing were used to understand the in situ ability for solubilization of these predominant forms of P. Locations with more occluded and organic P harbored the greatest abundance of P-mobilizing microorganisms, especially Burkholderiaceae (Caballeronia and Paraburkholderia spp.). Nearly all bacteria utilized aluminum phosphate, however fewer could subsist on iron phosphate (FePO4) or phytate. Microorganisms isolated from phytic acid were also most effective at solubilizing FePO4, suggesting that phytate solubilization may be linked to the ability to solubilize Fe. Significantly, we observed Fe to be co-located with P in organic patches in soil. Siderophore addition in lab experiments reinstated phytase mediated P-solubilization from Fe-phytate complexes. Taken together, these results indicate that metal-organic-P complex formation may limit enzymatic P solubilization from phytate in soil. Additionally, the linked traits of phytase and siderophore production were mostly restricted to specific clades within the Burkholderiaceae. We propose that Fe complexation of organic P (e.g., phytate) represents a major constraint on P turnover and availability in acidic soils, as only a limited subset of bacteria appear to possess the traits required to access this persistent pool of soil P.

Young women in sub-Saharan Africa are disproportionally affected by HIV infection and unintended pregnancies. However, hormonal contraceptive (HC) use may influence HIV risk through changes in genital tract microbiota and inflammatory cytokines. To investigate this, 130 HIV negative adolescent females aged 15-19 years were enrolled into a substudy of UChoose, an open-label randomized crossover study (NCT02404038), comparing acceptability and contraceptive product preference as a proxy for HIV prevention delivery methods. Participants were randomized to injectable norethisterone enanthate (Net-En), combined oral contraceptives (COC) or etonorgesterol/ethinyl estradiol combined contraceptive vaginal ring (CCVR) for 16 weeks, then crossed over to another HC for 16 weeks. Cervicovaginal samples were collected at baseline, crossover and exit for characterization of the microbiota and measurement of cytokine levels; primary endpoints were cervical T cell activation, vaginal microbial diversity and cytokine concentrations. Adolescents randomized to COCs had lower vaginal microbial diversity and relative abundance of HIV risk-associated taxa compared to Net-En or CCVR. Cervicovaginal inflammatory cytokine concentrations were significantly higher in adolescents randomized to CCVR compared to COC and Net-En. This suggests that COC use may induce an optimal vaginal ecosystem by decreasing bacterial diversity and inflammatory taxa, while CCVR use is associated with genital inflammation.

The rhizosphere is a hotspot for microbial carbon transformations, and is the entry point for root polysaccharides and polymeric carbohydrates that are important precursors to soil organic matter (SOM). However, the ecological mechanisms that underpin rhizosphere carbohydrate depolymerization are poorly understood. Using Avena fatua, a common annual grass, we analyzed time-resolved metatranscriptomes to compare microbial functions in rhizosphere, detritusphere, and combined rhizosphere-detritusphere habitats. Transcripts were binned using a unique reference database generated from soil isolate genomes, single-cell amplified genomes, metagenomes, and stable isotope probing metagenomes. While soil habitat significantly affected both community composition and overall gene expression, the succession of microbial functions occurred at a faster time scale than compositional changes. Using hierarchical clustering of upregulated decomposition genes, we identified four distinct microbial guilds populated by taxa whose functional succession patterns suggest specialization for substrates provided by fresh growing roots, decaying root detritus, the combination of live and decaying root biomass, or aging root material. Carbohydrate depolymerization genes were consistently upregulated in the rhizosphere, and both taxonomic and functional diversity were highest in the combined rhizosphere-detritusphere, suggesting coexistence of rhizosphere guilds is facilitated by niche differentiation. Metatranscriptome-defined guilds provide a framework to model rhizosphere succession and its consequences for soil carbon cycling.

Soil microbial biomass can reach its annual maximum pool size beneath the winter snowpack and is known to decline abruptly following snowmelt in seasonally snow-covered ecosystems. Observed differences in winter versus summer microbial taxonomic composition also suggests that phylogenetically conserved traits may permit winter- versus summer-adapted microorganisms to occupy distinct niches. In this study, we sought to identify archaea, bacteria, and fungi that are associated with the soil microbial bloom overwinter and the subsequent biomass collapse following snowmelt at a high-altitude watershed in central Colorado, United States. Archaea, bacteria, and fungi were categorized into three life strategies (Winter-Adapted, Snowmelt-Specialist, Spring-Adapted) based upon changes in abundance during winter, the snowmelt period, and after snowmelt in spring. We calculated indices of phylogenetic relatedness (archaea and bacteria) or assigned functional attributes (fungi) to organisms within life strategies to infer whether phylogenetically conserved traits differentiate Winter-Adapted, Snowmelt-Specialist, and Spring-Adapted groups. We observed that the soil microbial bloom was correlated in time with a pulse of snowmelt infiltration, which commenced 65 days prior to soils becoming snow-free. A pulse of nitrogen (N, as nitrate) occurred after snowmelt, along with a collapse in the microbial biomass pool size, and an increased abundance of nitrifying archaea and bacteria (e.g., Thaumarchaeota, Nitrospirae). Winter- and Spring-Adapted archaea and bacteria were phylogenetically clustered, suggesting that phylogenetically conserved traits allow Winter- and Spring-Adapted archaea and bacteria to occupy distinct niches. In contrast, Snowmelt-Specialist archaea and bacteria were phylogenetically overdispersed, suggesting that the key mechanism(s) of the microbial biomass crash are likely to be density-dependent (e.g., trophic interactions, competitive exclusion) and affect organisms across a broad phylogenetic spectrum. Saprotrophic fungi were the dominant functional group across fungal life strategies, however, ectomycorrhizal fungi experienced a large increase in abundance in spring. If well-coupled plant-mycorrhizal phenology currently buffers ecosystem N losses in spring, then changes in snowmelt timing may alter ecosystem N retention potential. Overall, we observed that snowmelt separates three distinct soil niches that are occupied by ecologically distinct groups of microorganisms. This ecological differentiation is of biogeochemical importance, particularly with respect to the mobilization of nitrogen during winter, before and after snowmelt.

For free-living bacteria and archaea, the equivalent of the biological species concept does not exist, creating several obstacles to the study of the processes contributing to microbial diversification. These obstacles are particularly high in soil, where high bacterial diversity inhibits the study of closely related genotypes and therefore the factors structuring microbial populations. Here, we isolated strains within a single Curtobacterium ecotype from surface soil (leaf litter) across a regional climate gradient and investigated the phylogenetic structure, recombination, and flexible gene content of this genomic diversity to infer patterns of gene flow. Our results indicate that microbial populations are delineated by gene flow discontinuities, with distinct populations cooccurring at multiple sites. Bacterial population structure was further delineated by genomic features allowing for the identification of candidate genes possibly contributing to local adaptation. These results suggest that the genetic structure within this bacterium is maintained both by ecological specialization in localized microenvironments (isolation by environment) and by dispersal limitation between geographic locations (isolation by distance).IMPORTANCE Due to the promiscuous exchange of genetic material and asexual reproduction, delineating microbial species (and, by extension, populations) remains challenging. Because of this, the vast majority of microbial studies assessing population structure often compare divergent strains from disparate environments under varied selective pressures. Here, we investigated the population structure within a single bacterial ecotype, a unit equivalent to a eukaryotic species, defined as highly clustered genotypic and phenotypic strains with the same ecological niche. Using a combination of genomic and computational analyses, we assessed the phylogenetic structure, extent of recombination, and flexible gene content of this genomic diversity to infer patterns of gene flow. To our knowledge, this study is the first to do so for a dominant soil bacterium. Our results indicate that bacterial soil populations, similarly to those in other environments, are structured by gene flow discontinuities and exhibit distributional patterns consistent with both isolation by distance and isolation by environment. Thus, both dispersal limitation and local environments contribute to the divergence among closely related soil bacteria as observed in macroorganisms.

Determining the carbon sources for active microbial populations in the subsurface is a challenging but highly informative component of subsurface microbial ecology. This work developed a method to provide ecological insights into groundwater microbial communities by characterizing community RNA through its radiocarbon and ribosomal RNA (rRNA) signatures. RNA was chosen as the biomolecule of interest because rRNA constitutes the majority of RNA in prokaryotes, represents recently active organisms, and yields detailed taxonomic information. The method was applied to a groundwater filter collected from a shallow alluvial aquifer in Colorado. RNA was extracted, radiometrically dated, and the 16S rRNA was analyzed by RNA-Seq. The RNA had a radiocarbon signature (Δ14C) of -193.4 ± 5.6‰. Comparison of the RNA radiocarbon signature to those of potential carbon pools in the aquifer indicated that at least 51% of the RNA was derived from autotrophy, in close agreement with the RNA-Seq data, which documented the prevalence of autotrophic taxa, such as Thiobacillus and Gallionellaceae. Overall, this hybrid method for RNA analysis provided cultivation-independent information on the in-situ carbon sources of active subsurface microbes and reinforced the importance of autotrophy and the preferential utilization of dissolved over sedimentary organic matter in alluvial aquifers.

Beneficial microbial associations enhance the fitness of most living organisms, and wood-feeding insects offer some of the most striking examples of this. Odontotaenius disjunctus is a wood-feeding beetle that possesses a digestive tract with four main compartments, each of which contains well-differentiated microbial populations, suggesting that anatomical properties and separation of these compartments may enhance energy extraction from woody biomass. Here, using integrated chemical analyses, we demonstrate that lignocellulose deconstruction and fermentation occur sequentially across compartments, and that selection for microbial groups and their metabolic pathways is facilitated by gut anatomical features. Metaproteogenomics showed that higher oxygen concentration in the midgut drives lignocellulose depolymerization, while a thicker gut wall in the anterior hindgut reduces oxygen diffusion and favours hydrogen accumulation, facilitating fermentation, homoacetogenesis and nitrogen fixation. We demonstrate that depolymerization continues in the posterior hindgut, and that the beetle excretes an energy- and nutrient-rich product on which its offspring subsist and develop. Our results show that the establishment of beneficial microbial partners within a host requires both the acquisition of the microorganisms and the formation of specific habitats within the host to promote key microbial metabolic functions. Together, gut anatomical properties and microbial functional assembly enable lignocellulose deconstruction and colony subsistence on an extremely nutrient-poor diet.

OBJECTIVES:Thousands of cases of multidrug-resistant tuberculosis (TB) have been observed in the Philippines, but studies on the Mycobacterium tuberculosis (MTB) genotypes that underlie the observed drug resistance profiles are lacking. This study aimed to analyse the whole genomes of clinical MTB isolates representing various resistance profiles to identify single nucleotide polymorphisms (SNPs) in resistance-associated genes. METHODS:The genomes of ten MTB isolates cultured from banked sputum sources were sequenced. Bioinformatics analysis consisted of assembly, annotation and SNP identification in genes reported to be associated with resistance to isoniazid (INH), rifampicin (RIF), ethambutol (ETH), streptomycin, pyrazinamide (PZA) and fluoroquinolones (FQs). RESULTS:The draft assemblies covered an average of 97.08% of the expected genome size. Seven of the ten isolates belonged to the Indo-Oceanic lineage/EA12-Manila clade. Two isolates were classified into the Euro-American lineage, whilst the pre-XDR (pre-extensively drug-resistant) isolate was classified under the East Asian/Beijing clade. The SNPs katG Ser315Thr, rpoB Ser450Leu and embB Met306Val were found in INH- (4/7), RIF- (3/6) and ETH-resistant (2/6) isolates, respectively, but not in susceptible isolates. Mutations in the inhA promoter and in the pncA and gyrA genes known to be involved in resistance to INH, PZA and FQs, respectively, were also identified. CONCLUSIONS:This study represents the first effort to investigate the whole genomes of Philippine clinical strains of MTB exhibiting various multidrug resistance profiles. Whole-genome data can provide valuable insights to the mechanistic and epidemiological qualities of TB in a high-burden setting such as the Philippines.

BackgroundExclusive breastfeeding reduces the rate of postnatal human immunodeficiency virus (HIV) transmission compared to nonexclusive breastfeeding; however, the mechanisms of this protection are unknown. Our study aimed to interrogate the mechanisms underlying the protective effect of exclusive breastfeeding.MethodsWe performed a prospective, longitudinal study of infants from a high-HIV-prevalence, low-income setting in South Africa. We evaluated the role of any non-breast milk feeds, excluding prescribed medicines on stool microbial communities via 16S rRNA gene sequencing, peripheral T-cell activation via flow cytometry, and buccal mucosal gene expression via quantitative polymerase chain reaction assay.ResultsA total of 155 infants were recruited at birth with mean gestational age of 38.9 weeks and mean birth weight of 3.2 kg. All infants were exclusively breastfed (EBF) at birth, but only 43.5% and 20% remained EBF at 6 or 14 weeks of age, respectively. We observed lower stool microbial diversity and distinct microbial composition in exclusively breastfed infants. These microbial communities, and the relative abundance of key taxa, were correlated with peripheral CD4+ T-cell activation, which was lower in EBF infants. In the oral mucosa, gene expression of chemokine and chemokine receptors involved in recruitment of HIV target cells to tissues, as well as epithelial cytoskeletal proteins, was lower in EBF infants.ConclusionsThese data suggest that nonexclusive breastfeeding alters the gut microbiota, increasing T-cell activation and, potentially, mucosal recruitment of HIV target cells. Study findings highlight a biologically plausible mechanistic explanation for the reduced postnatal HIV transmission observed in EBF infants.

Biological soil crusts (biocrusts) account for a substantial portion of primary production in dryland ecosystems. They successionally mature to deliver a suite of ecosystem services, such as carbon sequestration, water retention and nutrient cycling, and climate regulation. Biocrust assemblages are extremely well adapted to survive desiccation and to rapidly take advantage of the periodic precipitation events typical of arid ecosystems. Here we focus on the wetting response of incipient cyanobacterial crusts as they mature from "light" to "dark." We sampled a cyanobacterial biocrust chronosequence before (dry) and temporally following a controlled wetting event and used high-throughput 16S rRNA and rRNA gene sequencing to monitor the dynamics of microbial response. Overall, shorter-term changes in phylogenetic beta diversity attributable to periodic wetting were as large as those attributable to biocrust successional stage. Notably, more mature crusts showed significantly higher resistance to precipitation disturbance. A large bloom of a few taxa within the Firmicutes, primarily in the order Bacillales, emerged 18 h after wetting, while filamentous crust-forming cyanobacteria showed variable responses to wet-up across the successional gradient, with populations collapsing in less-developed light crusts but increasing in later-successional-stage dark crusts. Overall, the consistent Bacillales bloom accompanied by the variable collapse of pioneer cyanobacteria of the Oscillatoriales order across the successional gradient suggests that the strong response of few organisms to a hydration pulse with the mortality of the autotroph might have important implications for carbon (C) balance in semiarid ecosystems.IMPORTANCE Desert biological soil crusts are terrestrial topsoil microbial communities common to arid regions that comprise 40% of Earth's terrestrial surface. They successionally develop over years to decades to deliver a suite of ecosystem services of local and global significance. Ecosystem succession toward maturity has been associated with both resistance and resilience to disturbance. Recent work has shown that the impacts of both climate change and physical disturbance on biocrusts increase the potential for successional resetting. A larger proportion of biocrusts are expected to be at an early developmental stage, hence increasing susceptibility to changes in precipitation frequencies. Therefore, it is essential to characterize how biocrusts respond to wetting across early developmental stages. In this study, we document the wetting response of microbial communities from a biocrust chronosequence. Overall, our results suggest that the cumulative effects of altered precipitation frequencies on the stability of biocrusts will depend on biocrust maturity.

Like all higher organisms, plants have evolved in the context of a microbial world, shaping both their evolution and their contemporary ecology. Interactions between plant roots and soil microorganisms are critical for plant fitness in natural environments. Given this co-evolution and the pivotal importance of plant-microbial interactions, it has been hypothesized, and a growing body of literature suggests, that plants may regulate the composition of their rhizosphere to promote the growth of microorganisms that improve plant fitness in a given ecosystem. Here, using a combination of comparative genomics and exometabolomics, we show that pre-programmed developmental processes in plants (Avena barbata) result in consistent patterns in the chemical composition of root exudates. This chemical succession in the rhizosphere interacts with microbial metabolite substrate preferences that are predictable from genome sequences. Specifically, we observed a preference by rhizosphere bacteria for consumption of aromatic organic acids exuded by plants (nicotinic, shikimic, salicylic, cinnamic and indole-3-acetic). The combination of these plant exudation traits and microbial substrate uptake traits interact to yield the patterns of microbial community assembly observed in the rhizosphere of an annual grass. This discovery provides a mechanistic underpinning for the process of rhizosphere microbial community assembly and provides an attractive direction for the manipulation of the rhizosphere microbiome for beneficial outcomes.

Metagenomic sequencing provides a window into microbial community structure and metabolic potential; however, linking these data to exogenous metabolites that microorganisms process and produce (the exometabolome) remains challenging. Previously, we observed strong exometabolite niche partitioning among bacterial isolates from biological soil crust (biocrust). Here we examine native biocrust to determine if these patterns are reproduced in the environment. Overall, most soil metabolites display the expected relationship (positive or negative correlation) with four dominant bacteria following a wetting event and across biocrust developmental stages. For metabolites that were previously found to be consumed by an isolate, 70% are negatively correlated with the abundance of the isolate's closest matching environmental relative in situ, whereas for released metabolites, 67% were positively correlated. Our results demonstrate that metabolite profiling, shotgun sequencing and exometabolomics may be successfully integrated to functionally link microbial community structure with environmental chemistry in biocrust.

Hypersaline photosynthetic microbial mats are stratified microbial communities known for their taxonomic and metabolic diversity and strong light-driven day-night environmental gradients. In this study of the upper photosynthetic zone of hypersaline microbial mats of Elkhorn Slough, California (USA), we show how metagenome sequencing can be used to meaningfully assess microbial ecology and genetic partitioning in these complex microbial systems. Mapping of metagenome reads to the dominant Cyanobacteria observed in the system, Coleofasciculus (Microcoleus) chthonoplastes, was used to examine strain variants within these metagenomes. Highly conserved gene subsystems indicated a core genome for the species, and a number of variant genes and subsystems suggested strain level differentiation, especially for nutrient utilization and stress response. Metagenome sequence coverage binning was used to assess ecosystem partitioning of remaining microbes to both reconstruct the model organisms in silico and identify their ecosystem functions as well as to identify novel clades and propose their role in the biogeochemical cycling of mats. Functional gene annotation of these bins (primarily of Proteobacteria, Bacteroidetes, and Cyanobacteria) recapitulated the known biogeochemical functions in microbial mats using a genetic basis, and revealed significant diversity in the Bacteroidetes, presumably in heterotrophic cycling. This analysis also revealed evidence of putative phototrophs within the Gemmatimonadetes and Gammaproteobacteria residing in microbial mats. This study shows that metagenomic analysis can produce insights into the systems biology of microbial ecosystems from a genetic perspective and to suggest further studies of novel microbes.

Biogenesis of hydrogen sulfide (H2S) (microbial souring) has detrimental impacts on oil production operations and can cause health and safety problems. Understanding the processes that control the rates and patterns of sulfate reduction is crucial in developing a predictive understanding of reservoir souring and associated mitigation processes. This work demonstrates an approach to utilize genomic information to constrain the biological parameters needed for modeling souring, providing a pathway for using microbial data derived from oil reservoir studies. Minimum generation times were calculated based on codon usage bias and optimal growth temperatures based on the frequency of amino acids. We show how these derived parameters can be used in a simplified multiphase reactive transport model by simulating the injection of cold (30 °C) seawater into a 70 °C reservoir, modeling the shift in sulfate reducing microorganisms (SRM) community composition, sulfate and sulfide concentrations through time and space. Finally, we explore the question of necessary model complexity by comparing results using different numbers of SRM. Simulations showed that the kinetics of a SRM community consisting of twenty-five SRM could be adequately represented by a reduced community consisting of nine SRM with parameter values derived from the mean and standard deviations of the original SRM.

Much genetic diversity within a bacterial community is likely obscured by microdiversity within operational taxonomic units (OTUs) defined by 16S rRNA gene sequences. However, it is unclear how variation within this microdiversity influences ecologically relevant traits. Here, we employ a multifaceted approach to investigate microdiversity within the dominant leaf litter bacterium, Curtobacterium, which comprises 7.8% of the bacterial community at a grassland site undergoing global change manipulations. We use cultured bacterial isolates to interpret metagenomic data, collected in situ over 2 years, together with lab-based physiological assays to determine the extent of trait variation within this abundant OTU. The response of Curtobacterium to seasonal variability and the global change manipulations, specifically an increase in relative abundance under decreased water availability, appeared to be conserved across six Curtobacterium lineages identified at this site. Genomic and physiological analyses in the lab revealed that degradation of abundant polymeric carbohydrates within leaf litter, cellulose and xylan, is nearly universal across the genus, which may contribute to its high abundance in grassland leaf litter. However, the degree of carbohydrate utilization and temperature preference for this degradation varied greatly among clades. Overall, we find that traits within Curtobacterium are conserved at different phylogenetic depths. We speculate that similar to bacteria in marine systems, diverse microbes within this taxon may be structured in distinct ecotypes that are key to understanding Curtobacterium abundance and distribution in the environment.IMPORTANCE Despite the plummeting costs of sequencing, characterizing the fine-scale genetic diversity of a microbial community-and interpreting its functional importance-remains a challenge. Indeed, most studies, particularly studies of soil, assess community composition at a broad genetic level by classifying diversity into taxa (OTUs) defined by 16S rRNA sequence similarity. However, these classifications potentially obscure variation in traits that result in fine-scale ecological differentiation among closely related strains. Here, we investigated "microdiversity" in a highly diverse and poorly characterized soil system (leaf litter in a southern Californian grassland). We focused on the most abundant bacterium, Curtobacterium, which by standard methods is grouped into only one OTU. We find that the degree of carbohydrate usage and temperature preference vary within the OTU, whereas its responses to changes in precipitation are relatively uniform. These results suggest that microdiversity may be key to understanding how soil bacterial diversity is linked to ecosystem functioning.

Organic matter deposits in alluvial aquifers have been shown to result in the formation of naturally reduced zones (NRZs), which can modulate aquifer redox status and influence the speciation and mobility of metals, affecting groundwater geochemistry. In this study, we sought to better understand how natural organic matter fuels microbial communities within anoxic biogeochemical hot spots (NRZs) in a shallow alluvial aquifer at the Rifle (CO) site. We conducted a 20-day microcosm experiment in which NRZ sediments, which were enriched in buried woody plant material, served as the sole source of electron donors and microorganisms. The microcosms were constructed and incubated under anaerobic conditions in serum bottles with an initial N2 headspace and were sampled every 5 days for metagenome and metatranscriptome profiles in combination with biogeochemical measurements. Biogeochemical data indicated that the decomposition of native organic matter occurred in different phases, beginning with mineralization of dissolved organic matter (DOM) to CO2 during the first week of incubation, followed by a pulse of acetogenesis that dominated carbon flux after 2 weeks. A pulse of methanogenesis co-occurred with acetogenesis, but only accounted for a small fraction of carbon flux. The depletion of DOM over time was strongly correlated with increases in expression of many genes associated with heterotrophy (e.g., amino acid, fatty acid, and carbohydrate metabolism) belonging to a Hydrogenophaga strain that accounted for a relatively large percentage (~8%) of the metatranscriptome. This Hydrogenophaga strain also expressed genes indicative of chemolithoautotrophy, including CO2 fixation, H2 oxidation, S-compound oxidation, and denitrification. The pulse of acetogenesis appears to have been collectively catalyzed by a number of different organisms and metabolisms, most prominently pyruvate:ferredoxin oxidoreductase. Unexpected genes were identified among the most highly expressed (>98th percentile) transcripts, including acetone carboxylase and cell-wall-associated hydrolases with unknown substrates (numerous lesser expressed cell-wall-associated hydrolases targeted peptidoglycan). Many of the most highly expressed hydrolases belonged to a Ca. Bathyarchaeota strain and may have been associated with recycling of bacterial biomass. Overall, these results highlight the complex nature of organic matter transformation in NRZs and the microbial metabolic pathways that interact to mediate redox status and elemental cycling.

Cyanobacterial mats are laminated microbial ecosystems which occur in highly diverse environments and which may provide a possible model for early life on Earth. Their ability to produce hydrogen also makes them of interest from a biotechnological and bioenergy perspective. Samples of an intertidal microbial mat from the Elkhorn Slough estuary in Monterey Bay, California, were transplanted to a greenhouse at NASA Ames Research Center to study a 24-h diel cycle, in the presence or absence of molybdate (which inhibits biohydrogen consumption by sulfate reducers). Here, we present metagenomic analyses of four samples that will be used as references for future metatranscriptomic analyses of this diel time series.

Cyanobacterial mats are laminated microbial ecosystems which occur in highly diverse environments and which may provide a possible model for early life on Earth. Their ability to produce hydrogen also makes them of interest from a biotechnological and bioenergy perspective. Samples of an intertidal microbial mat from the Elkhorn Slough estuary in Monterey Bay, California, were transplanted to a greenhouse at NASA Ames Research Center to study a 24-h diel cycle, in the presence or absence of molybdate (which inhibits biohydrogen consumption by sulfate reducers). Here, we present metagenomic analyses of four samples that will be used as references for future metatranscriptomic analyses of this diel time series.

The subterranean world hosts up to one-fifth of all biomass, including microbial communities that drive transformations central to Earth's biogeochemical cycles. However, little is known about how complex microbial communities in such environments are structured, and how inter-organism interactions shape ecosystem function. Here we apply terabase-scale cultivation-independent metagenomics to aquifer sediments and groundwater, and reconstruct 2,540 draft-quality, near-complete and complete strain-resolved genomes that represent the majority of known bacterial phyla as well as 47 newly discovered phylum-level lineages. Metabolic analyses spanning this vast phylogenetic diversity and representing up to 36% of organisms detected in the system are used to document the distribution of pathways in coexisting organisms. Consistent with prior findings indicating metabolic handoffs in simple consortia, we find that few organisms within the community can conduct multiple sequential redox transformations. As environmental conditions change, different assemblages of organisms are selected for, altering linkages among the major biogeochemical cycles.

Groundwater ecosystems are conventionally thought to be fueled by surface-derived allochthonous organic matter and dominated by heterotrophic microbes living under often-oligotrophic conditions. However, in a 2-month study of nitrate amendment to a perennially suboxic aquifer in Rifle (CO), strain-resolved metatranscriptomic analysis revealed pervasive and diverse chemolithoautotrophic bacterial activity relevant to C, S, N and Fe cycling. Before nitrate injection, anaerobic ammonia-oxidizing (anammox) bacteria accounted for 16% of overall microbial community gene expression, whereas during the nitrate injection, two other groups of chemolithoautotrophic bacteria collectively accounted for 80% of the metatranscriptome: (1) members of the Fe(II)-oxidizing Gallionellaceae family and (2) strains of the S-oxidizing species, Sulfurimonas denitrificans. Notably, the proportion of the metatranscriptome accounted for by these three groups was considerably greater than the proportion of the metagenome coverage that they represented. Transcriptional analysis revealed some unexpected metabolic couplings, in particular, putative nitrate-dependent Fe(II) and S oxidation among nominally microaerophilic Gallionellaceae strains, including expression of periplasmic (NapAB) and membrane-bound (NarGHI) nitrate reductases. The three most active groups of chemolithoautotrophic bacteria in this study had overlapping metabolisms that allowed them to occupy different yet related metabolic niches throughout the study. Overall, these results highlight the important role that chemolithoautotrophy can have in aquifer biogeochemical cycling, a finding that has broad implications for understanding terrestrial carbon cycling and is supported by recent studies of geochemically diverse aquifers.

Soil surface temperature, an important driver of terrestrial biogeochemical processes, depends strongly on soil albedo, which can be significantly modified by factors such as plant cover. In sparsely vegetated lands, the soil surface can be colonized by photosynthetic microbes that build biocrust communities. Here we use concurrent physical, biochemical and microbiological analyses to show that mature biocrusts can increase surface soil temperature by as much as 10 °C through the accumulation of large quantities of a secondary metabolite, the microbial sunscreen scytonemin, produced by a group of late-successional cyanobacteria. Scytonemin accumulation decreases soil albedo significantly. Such localized warming has apparent and immediate consequences for the soil microbiome, inducing the replacement of thermosensitive bacterial species with more thermotolerant forms. These results reveal that not only vegetation but also microorganisms are a factor in modifying terrestrial albedo, potentially impacting biosphere feedbacks on past and future climate, and call for a direct assessment of such effects at larger scales.

Soils are arguably the most microbially diverse ecosystems. Physicochemical properties have been associated with the maintenance of this diversity. Yet, the role of microbial substrate specialization is largely unexplored since substrate utilization studies have focused on simple substrates, not the complex mixtures representative of the soil environment. Here we examine the exometabolite composition of desert biological soil crusts (biocrusts) and the substrate preferences of seven biocrust isolates. The biocrust's main primary producer releases a diverse array of metabolites, and isolates of physically associated taxa use unique subsets of the complex metabolite pool. Individual isolates use only 13-26% of available metabolites, with only 2 out of 470 used by all and 40% not used by any. An extension of this approach to a mesophilic soil environment also reveals high levels of microbial substrate specialization. These results suggest that exometabolite niche partitioning may be an important factor in the maintenance of microbial diversity.

The coffee berry borer (Hypothenemus hampei) is the most devastating insect pest of coffee worldwide with its infestations decreasing crop yield by up to 80%. Caffeine is an alkaloid that can be toxic to insects and is hypothesized to act as a defence mechanism to inhibit herbivory. Here we show that caffeine is degraded in the gut of H. hampei, and that experimental inactivation of the gut microbiota eliminates this activity. We demonstrate that gut microbiota in H. hampei specimens from seven major coffee-producing countries and laboratory-reared colonies share a core of microorganisms. Globally ubiquitous members of the gut microbiota, including prominent Pseudomonas species, subsist on caffeine as a sole source of carbon and nitrogen. Pseudomonas caffeine demethylase genes are expressed in vivo in the gut of H. hampei, and re-inoculation of antibiotic-treated insects with an isolated Pseudomonas strain reinstates caffeine-degradation ability confirming their key role.

Advancements in molecular technologies have provided new platforms that are being increasingly adopted for use in the clinical microbiology laboratory. Among these, microarray methods are particularly well suited for diagnostics as they allow multiplexing, or the ability to test for multiple targets simultaneously from the same specimen. Microarray technologies commonly used for the detection and identification of microbial targets include solid-state microarrays, electronic microarrays and bead suspension microarrays. Microarray methods have been applied to microbial detection, genotyping and antimicrobial resistance gene detection. Microarrays can offer a panel approach to diagnose specific patient presentations, such as respiratory or gastrointestinal infections, and can discriminate isolates by genotype for tracking epidemiology and outbreak investigations. And, as more information has become available on specific genes and pathways involved in antimicrobial resistance, we are beginning to be able to predict susceptibility patterns based on sequence detection for particular organisms. With further advances in automated microarray processing methods and genotype-phenotype prediction algorithms, these tests will become even more useful as an adjunct or replacement for conventional antimicrobial susceptibility testing, allowing for more rapid selection of targeted therapy for infectious diseases.

Biological Soil Crusts (BSCs) are organosedimentary assemblages comprised of microbes and minerals in topsoil of terrestrial environments. BSCs strongly impact soil quality in dryland ecosystems (e.g., soil structure and nutrient yields) due to pioneer species such as Microcoleus vaginatus; phototrophs that produce filaments that bind the soil together, and support an array of heterotrophic microorganisms. These microorganisms in turn contribute to soil stability and biogeochemistry of BSCs. Non-cyanobacterial populations of BSCs are less well known than cyanobacterial populations. Therefore, we attempted to isolate a broad range of numerically significant and phylogenetically representative BSC aerobic heterotrophs. Combining simple pre-treatments (hydration of BSCs under dark and light) and isolation strategies (media with varying nutrient availability and protection from oxidative stress) we recovered 402 bacterial and one fungal isolate in axenic culture, which comprised 116 phylotypes (at 97% 16S rRNA gene sequence homology), 115 bacterial and one fungal. Each medium enriched a mostly distinct subset of phylotypes, and cultivated phylotypes varied due to the BSC pre-treatment. The fraction of the total phylotype diversity isolated, weighted by relative abundance in the community, was determined by the overlap between isolate sequences and OTUs reconstructed from metagenome or metatranscriptome reads. Together, more than 8% of relative abundance of OTUs in the metagenome was represented by our isolates, a cultivation efficiency much larger than typically expected from most soils. We conclude that simple cultivation procedures combined with specific pre-treatment of samples afford a significant reduction in the culturability gap, enabling physiological and metabolic assays that rely on ecologically relevant axenic cultures.

Biological Soil Crusts (BSCs) are organosedimentary assemblages comprised of microbes and minerals in topsoil of terrestrial environments. BSCs strongly impact soil quality in dryland ecosystems (e.g., soil structure and nutrient yields) due to pioneer species such as Microcoleus vaginatus; phototrophs that produce filaments that bind the soil together, and support an array of heterotrophic microorganisms. These microorganisms in turn contribute to soil stability and biogeochemistry of BSCs. Non-cyanobacterial populations of BSCs are less well known than cyanobacterial populations. Therefore, we attempted to isolate a broad range of numerically significant and phylogenetically representative BSC aerobic heterotrophs. Combining simple pre-treatments (hydration of BSCs under dark and light) and isolation strategies (media with varying nutrient availability and protection from oxidative stress) we recovered 402 bacterial and one fungal isolate in axenic culture, which comprised 116 phylotypes (at 97% 16S rRNA gene sequence homology), 115 bacterial and one fungal. Each medium enriched a mostly distinct subset of phylotypes, and cultivated phylotypes varied due to the BSC pre-treatment. The fraction of the total phylotype diversity isolated, weighted by relative abundance in the community, was determined by the overlap between isolate sequences and OTUs reconstructed from metagenome or metatranscriptome reads. Together, more than 8% of relative abundance of OTUs in the metagenome was represented by our isolates, a cultivation efficiency much larger than typically expected from most soils. We conclude that simple cultivation procedures combined with specific pre-treatment of samples afford a significant reduction in the culturability gap, enabling physiological and metabolic assays that rely on ecologically relevant axenic cultures.

UnlabelledHuman papillomavirus (HPV) causes a number of neoplastic diseases in humans. Here, we show a complex normal HPV community in a cohort of 103 healthy human subjects, by metagenomics analysis of the shotgun sequencing data generated from the NIH Human Microbiome Project. The overall HPV prevalence was 68.9% and was highest in the skin (61.3%), followed by the vagina (41.5%), mouth (30%), and gut (17.3%). Of the 109 HPV types as well as additional unclassified types detected, most were undetectable by the widely used commercial kits targeting the vaginal/cervical HPV types. These HPVs likely represent true HPV infections rather than transitory exposure because of strong organ tropism and persistence of the same HPV types in repeat samples. Coexistence of multiple HPV types was found in 48.1% of the HPV-positive samples. Networking between HPV types, cooccurrence or exclusion, was detected in vaginal and skin samples. Large contigs assembled from short HPV reads were obtained from several samples, confirming their genuine HPV origin. This first large-scale survey of HPV using a shotgun sequencing approach yielded a comprehensive map of HPV infections among different body sites of healthy human subjects.ImportanceThis nonbiased survey indicates that the HPV community in healthy humans is much more complex than previously defined by widely used kits that are target selective for only a few high- and low-risk HPV types for cervical cancer. The importance of nononcogenic viruses in a mixed HPV infection could be for stimulating or inhibiting a coexisting oncogenic virus via viral interference or immune cross-reaction. Knowledge gained from this study will be helpful to guide the designing of epidemiological and clinical studies in the future to determine the impact of nononcogenic HPV types on the outcome of HPV infections.

Coarse woody debris is an important biomass pool in forest ecosystems that numerous groups of insects have evolved to take advantage of. These insects are ecologically important and represent useful natural analogs for biomass to biofuel conversion. Using a range of molecular approaches combined with microelectrode measurements of oxygen, we have characterized the gut microbiome and physiology of Odontotaenius disjunctus, a wood-feeding beetle native to the eastern United States. We hypothesized that morphological and physiological differences among gut regions would correspond to distinct microbial populations and activities. In fact, significantly different communities were found in the foregut (FG), midgut (MG)/posterior hindgut (PHG) and anterior hindgut (AHG), with Actinobacteria and Rhizobiales being more abundant toward the FG and PHG. Conversely, fermentative bacteria such as Bacteroidetes and Clostridia were more abundant in the AHG, and also the sole region where methanogenic Archaea were detected. Although each gut region possessed an anaerobic core, micron-scale profiling identified radial gradients in oxygen concentration in all regions. Nitrogen fixation was confirmed by (15)N2 incorporation, and nitrogenase gene (nifH) expression was greatest in the AHG. Phylogenetic analysis of nifH identified the most abundant transcript as related to Ni-Fe nitrogenase of a Bacteroidetes species, Paludibacter propionicigenes. Overall, we demonstrate not only a compartmentalized microbiome in this beetle digestive tract but also sharp oxygen gradients that may permit aerobic and anaerobic metabolism to occur within the same regions in close proximity. We provide evidence for the microbial fixation of N2 that is important for this beetle to subsist on woody biomass.

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

Abstract Background Terabyte-scale collections of string-encoded data are expected from consortia efforts such as the Human Microbiome Project http://nihroadmap.nih.gov/hmp. Intra- and inter-project data similarity searches are enabled by rapid k-mer matching strategies. Software applications for sequence database partitioning, guide tree estimation, molecular classification and alignment acceleration have benefited from embedded k-mer searches as sub-routines. However, a rapid, general-purpose, open-source, flexible, stand-alone k-mer tool has not been available. Results Here we present a stand-alone utility, Simrank, which allows users to rapidly identify database strings the most similar to query strings. Performance testing of Simrank and related tools against DNA, RNA, protein and human-languages found Simrank 10X to 928X faster depending on the dataset. Conclusions Simrank provides molecular ecologists with a high-throughput, open source choice for comparing large sequence sets to find similarity.

Soils are immensely diverse microbial habitats with thousands of co-existing bacterial, archaeal, and fungal species. Across broad spatial scales, factors such as pH and soil moisture appear to determine the diversity and structure of soil bacterial communities. Within any one site however, bacterial taxon diversity is high and factors maintaining this diversity are poorly resolved. Candidate factors include organic substrate availability and chemical recalcitrance, and given that they appear to structure bacterial communities at the phylum level, we examine whether these factors might structure bacterial communities at finer levels of taxonomic resolution. Analyzing 16S rRNA gene composition of nucleotide analog-labeled DNA by PhyloChip microarrays, we compare relative growth rates on organic substrates of increasing chemical recalcitrance of >2,200 bacterial taxa across 43 divisions/phyla. Taxa that increase in relative abundance with labile organic substrates (i.e., glycine, sucrose) are numerous (>500), phylogenetically clustered, and occur predominantly in two phyla (Proteobacteria and Actinobacteria) including orders Actinomycetales, Enterobacteriales, Burkholderiales, Rhodocyclales, Alteromonadales, and Pseudomonadales. Taxa increasing in relative abundance with more chemically recalcitrant substrates (i.e., cellulose, lignin, or tannin-protein) are fewer (168) but more phylogenetically dispersed, occurring across eight phyla and including Clostridiales, Sphingomonadalaes, Desulfovibrionales. Just over 6% of detected taxa, including many Burkholderiales increase in relative abundance with both labile and chemically recalcitrant substrates. Estimates of median rRNA copy number per genome of responding taxa demonstrate that these patterns are broadly consistent with bacterial growth strategies. Taken together, these data suggest that changes in availability of intrinsically labile substrates may result in predictable shifts in soil bacterial composition.

Bacterial communities in the airways of cystic fibrosis (CF) patients are, as in other ecological niches, influenced by autogenic and allogenic factors. However, our understanding of microbial colonization in younger versus older CF airways and the association with pulmonary function is rudimentary at best. Using a phylogenetic microarray, we examine the airway microbiota in age stratified CF patients ranging from neonates (9 months) to adults (72 years). From a cohort of clinically stable patients, we demonstrate that older CF patients who exhibit poorer pulmonary function possess more uneven, phylogenetically-clustered airway communities, compared to younger patients. Using longitudinal samples collected form a subset of these patients a pattern of initial bacterial community diversification was observed in younger patients compared with a progressive loss of diversity over time in older patients. We describe in detail the distinct bacterial community profiles associated with young and old CF patients with a particular focus on the differences between respective "early" and "late" colonizing organisms. Finally we assess the influence of Cystic Fibrosis Transmembrane Regulator (CFTR) mutation on bacterial abundance and identify genotype-specific communities involving members of the Pseudomonadaceae, Xanthomonadaceae, Moraxellaceae and Enterobacteriaceae amongst others. Data presented here provides insights into the CF airway microbiota, including initial diversification events in younger patients and establishment of specialized communities of pathogens associated with poor pulmonary function in older patient populations.