The effect of anthropogenic arsenic contamination on the earthworm microbiome

Earthworms are globally distributed and perform essential roles for soil health and microbial structure. We have investigated the effect of an anthropogenic contamination gradient on the bacterial community of the keystone ecological species Lumbricus rubellus through utilizing 16S rRNA pyrosequencing for the first time to establish the microbiome of the host and surrounding soil. The earthworm‐associated microbiome differs from the surrounding environment which appears to be a result of both filtering and stimulation likely linked to the altered environment associated with the gut micro‐habitat (neutral pH, anoxia and increased carbon substrates). We identified a core earthworm community comprising Proteobacteria (～50%) and Actinobacteria (～30%), with lower abundances of Bacteroidetes (～6%) and Acidobacteria (～3%). In addition to the known earthworm symbiont (Verminephrobacter sp.), we identified a potential host‐associated Gammaproteobacteria species (Serratia sp.) that was absent from soil yet observed in most earthworms. Although a distinct bacterial community defines these earthworms, clear family‐ and species‐level modification were observed along an arsenic and iron contamination gradient. Several taxa observed in uncontaminated control microbiomes are suppressed by metal/metalloid field exposure, including eradication of the hereto ubiquitously associated Verminephrobacter symbiont, which raises implications to its functional role in the earthworm microbiome.     

The core microbiome bonds the Alpine bog vegetation to a transkingdom metacommunity

Bog ecosystems fulfil important functions in Earth's carbon and water turnover. While plant communities and their keystone species Sphagnum have been well studied, less is known about the microbial communities associated with them. To study our hypothesis that bog plants share an essential core of their microbiome despite their different phylogenetic origins, we analysed four plant community plots with 24 bryophytes, vascular plants and lichen species in two Alpine bogs in Austria by 16S rDNA amplicon sequencing followed by bioinformatic analyses. The overall bog microbiome was classified into 32 microbial phyla, while Proteobacteria (30.8%), Verrucomicrobia (20.3%) and Planctomycetes (15.1%) belonged to the most abundant groups. Interestingly, the archaeal phylum Euryarcheota represented 7.2% of total microbial abundance. However, a high portion of micro‐organisms remained unassigned at phylum and class level, respectively. The core microbiome of the bog vegetation contained 177 operational taxonomic units (OTUs) (150 526 seq.) and contributed to 49.5% of the total microbial abundance. Only a minor portion of associated core micro‐organisms was host specific for examined plant groups (5.9–11.6%). Using our new approach to analyse plant–microbial communities in an integral framework of ecosystem, vegetation and microbiome, we demonstrated that bog vegetation harboured a core microbiome that is shared between plants and lichens over the whole ecosystem and formed a transkingdom metacommunity. All micro‐ and macro‐organisms are connected to keystone Sphagnum mosses via set of microbial species, for example Burkholderia bryophila which was found associated with a wide spectrum of host plants and is known for a beneficial plant–microbe interaction.     

Environmental switching during biofilm development in a cold seep system and functional determinants of species sorting

The functional basis for species sorting theory remains elusive, especially for microbial community assembly in deep‐sea environments. Using artificial surface‐based biofilm models, our recent work revealed taxonomic succession during biofilm development in a newly defined cold seep system, the Thuwal cold seeps II, which comprises a brine pool and the adjacent normal bottom water (NBW) to form a metacommunity via the potential immigration of organisms from one patch to another. Here, we designed an experiment to investigate the effects of environmental switching between the brine pool and the NBW on biofilm assembly, which could reflect environmental filtering effects during bacterial immigration to new environments. Analyses of 16S rRNA genes of 71 biofilm samples suggested that the microbial composition of biofilms established in new environments was determined by both the source community and the incubation conditions. Moreover, a comparison of 18 metagenomes provided evidence for biofilm community assembly that was based primarily on functional features rather than taxonomic identities; metal ion resistance and amino acid metabolism were the major species sorting determinants for the succession of biofilm communities. Genome binning and pathway reconstruction of two bacterial species (Marinobacter sp. and Oleispira sp.) further demonstrated metal ion resistance and amino acid metabolism as functional traits conferring the survival of habitat generalists in both the brine pool and NBW. The results of this study shed new light on microbial community assembly in special habitats and bridge a gap in species sorting theory.     

Networking in the Plant Microbiome

Almost all higher organisms, including plants, insects, and mammals, are colonized by complex microbial communities and harbor a microbiome. Emerging studies with plants reveal that these microbiomes are structured and form complex, interconnected microbial networks. Within these networks, different taxa have different roles, and keystone species have been identified that could be crucial for plant health and ecosystem functioning. A new paper in this issue of PLOS Biology by Agler et al. highlights the presence of microbial hubs in these networks that may act as mediators between the plant and its microbiome. A next major frontier is now to link microbiome composition to function. In order to do this, we present a number of hypothetical examples of how microbiome diversity and function potentially influence host performance.     

Assembly of the amphibian microbiome is influenced by the effects of land-use change on environmental reservoirs

A growing focus in microbial ecology is understanding how beneficial microbiome function is created and maintained through various assembly mechanisms. This study explores the role of both the environment and disease in regulating the composition of microbial species in the soil and on amphibian hosts. We compared the microbial communities of Plethodon cinereus salamanders along a land-use gradient in the New York metropolitan area and paired these with associated soil cores. Additionally, we characterized the diversity of bacterial and fungal symbionts that putatively inhibit the pathogenic fungus Batrachochytrium dendrobatidis. We predicted that variation in skin microbial community composition would correlate with changes seen in the soil which functions as the regional species pool. We found that salamanders and soil share many microbial taxa but that these two communities exhibit differences in the relative abundances of the bacterial phyla Acidobacteria, Actinobacteria, and Proteobacteria and the fungal phyla Ascomycota and genus Basidiobolus. Microbial community composition varies with changes in land-use associated factors creating site-specific compositions. By employing a quantitative, null-based assembly model, we identified that dispersal limitation, variable selection, and drift guide assembly of microbes onto their skin, creating high dissimilarity between individuals with likely consequences in disease preventative function.     

Identifying Keystone Species in the Human Gut Microbiome from Metagenomic Timeseries Using Sparse Linear Regression

Human associated microbial communities exert tremendous influence over human health and disease. With modern metagenomic sequencing methods it is now possible to follow the relative abundance of microbes in a community over time. These microbial communities exhibit rich ecological dynamics and an important goal of microbial ecology is to infer the ecological interactions between species directly from sequence data. Any algorithm for inferring ecological interactions must overcome three major obstacles: 1) a correlation between the abundances of two species does not imply that those species are interacting, 2) the sum constraint on the relative abundances obtained from metagenomic studies makes it difficult to infer the parameters in timeseries models, and 3) errors due to experimental uncertainty, or mis-assignment of sequencing reads into operational taxonomic units, bias inferences of species interactions due to a statistical problem called “errors-in-variables”. Here we introduce an approach, Learning Interactions from MIcrobial Time Series (LIMITS), that overcomes these obstacles. LIMITS uses sparse linear regression with boostrap aggregation to infer a discrete-time Lotka-Volterra model for microbial dynamics. We tested LIMITS on synthetic data and showed that it could reliably infer the topology of the inter-species ecological interactions. We then used LIMITS to characterize the species interactions in the gut microbiomes of two individuals and found that the interaction networks varied significantly between individuals. Furthermore, we found that the interaction networks of the two individuals are dominated by distinct “keystone species”, Bacteroides fragilis and Bacteroided stercosis, that have a disproportionate influence on the structure of the gut microbiome even though they are only found in moderate abundance. Based on our results, we hypothesize that the abundances of certain keystone species may be responsible for individuality in the human gut microbiome.     

Microbiome Profiles in Periodontitis in Relation to Host and Disease Characteristics

Periodontitis is an inflammatory condition that affects the supporting tissues surrounding teeth. The occurrence of periodontitis is associated with shifts in the structure of the communities that inhabit the gingival sulcus. Although great inter-subject variability in the subgingival microbiome has been observed in subjects with periodontitis, it is unclear whether distinct community types exist and if differences in microbial signatures correlate with host characteristics or with the variable clinical presentations of periodontitis. Therefore, in this study we explored the existence of different community types in periodontitis and their relationship with host demographic, medical and disease-related clinical characteristics. Clustering analyses of microbial abundance profiles suggested two types of communities (A and B) existed in the 34 subjects with periodontitis evaluated. Type B communities harbored greater proportions of certain periodontitis-associated taxa, including species historically associated with the disease, such as Porphyromonas gingivalis, Tannerella forsythia and Treponema denticola, and taxa recently linked to periodontitis. In contrast, subjects with type A communities had increased proportions of different periodontitis-associated species, and were also enriched for health-associated species and core taxa (those equally prevalent in health and periodontitis). Periodontitis subgingival clusters were not associated with demographic, medical or disease-specific clinical parameters other than periodontitis extent (proportion of sites affected), which positively correlated with the total proportion of cluster B signature taxa. In conclusion, two types of microbial communities were detected in subjects with periodontitis. Host demographics and underlying medical conditions did not correlate with these profiles, which instead appeared to be related to periodontitis extent, with type B communities present in more widespread disease cases. The two identified periodontitis profiles may represent distinct dysbiotic processes potentially requiring community-tailored therapeutic interventions.     

Elevated nutrients change bacterial community composition and connectivity: high throughput sequencing of young marine biofilms

Biofilms are integral to many marine processes but their formation and function may be affected by anthropogenic inputs that alter environmental conditions, including fertilisers that increase nutrients. Density composition and connectivity of biofilms developed in situ (under ambient and elevated nutrients) were compared using 454-pyrosequencing of the 16S gene. Elevated nutrients shifted community composition from bacteria involved in higher processes (eg Pseudoalteromonas spp. invertebrate recruitment) towards more nutrient-tolerant bacterial species (eg Terendinibacter sp.). This may enable the persistence of biofilm communities by increasing resistance to nutrient inputs. A core biofilm microbiome was identified (predominantly Alteromonadales and Oceanospirillales) and revealed shifts in abundances of core microbes that could indicate enrichment by fertilisers. Fertiliser decreased density and connectivity within biofilms indicating that associations were disrupted perhaps via changes to energetic allocations within the core microbiome. Density composition and connectivity changes suggest nutrients can affect the stability and function of these important marine communities.     

A metabarcoding analysis of the wrackbed microbiome indicates a phylogeographic break along the North Sea–Baltic Sea transition zone

Sandy beaches are biogeochemical hotspots that bridge marine and terrestrial ecosystems via the transfer of organic matter, such as seaweed (termed wrack). A keystone of this unique ecosystem is the microbial community, which helps to degrade wrack and re-mineralize nutrients. However, little is known about this community. Here, we characterize the wrackbed microbiome as well as the microbiome of a primary consumer, the seaweed fly Coelopa frigida, and examine how they change along one of the most studied ecological gradients in the world, the transition from the marine North Sea to the brackish Baltic Sea. We found that polysaccharide degraders dominated both microbiomes, but there were still consistent differences between wrackbed and fly samples. Furthermore, we observed a shift in both microbial communities and functionality between the North and Baltic Sea driven by changes in the frequency of different groups of known polysaccharide degraders. We hypothesize that microbes were selected for their abilities to degrade different polysaccharides corresponding to a shift in polysaccharide content in the different seaweed communities. Our results reveal the complexities of both the wrackbed microbial community, with different groups specialized to different roles, and the cascading trophic consequences of shifts in the near shore algal community.     

Phylogenetic Profiles of In-House Microflora in Drains at a Food Production Facility: Comparison and Biocontrol Implications of Listeria-Positive and -Negative Bacterial Populations

Listeria species experience complex interactions with other microorganisms, which may promote growth and colonization of the organism in local environments or negatively affect them. This study investigated the microbial community at a food production facility, examining interactions between Listeria and the associated microbiome. Listeria species can be transferred between zones in the production environment by individuals or equipment, and drains may act as a reservoir for the organism, reflecting the microbial flora potentially in the production environment. Drains that were colonized by Listeria species and those determined to be free of Listeria were examined. In each case, 16S rRNA gene analysis was performed using the PhyloChip platform. Some general similarities in bacterial population structure were observed when Listeria-negative and -positive drain communities were compared, with some distinct differences also noted. These included increased populations of the genera Prevotella and Janthinobacterium associated with the absence of Listeria species, whereas Enterococcus and Rhodococcus were in higher abundance in drains colonized by Listeria species. Based on these results, a selection of bacterial species were grown in coculture biofilm with a Listeria monocytogenes strain identified as having colonized a drain at the facility. Mixed-species biofilm experiments showed that Janthinobacterium inhibited attachment and subsequent biofilm formation of L. monocytogenes; however, Enterococcus gallinarum significantly increased it. The results of this study suggest the microbial community in food processing facilities can impact the colonization of Listeria species and that influencing the microbiome in favor of antilisterial species may reduce the colonization of Listeria species and limit the likelihood of product/process contamination.     

Fungal and Bacterial Communities in Indoor Dust Follow Different Environmental Determinants

People spend most of their time inside buildings and the indoor microbiome is a major part of our everyday environment. It affects humans’ wellbeing and therefore its composition is important for use in inferring human health impacts. It is still not well understood how environmental conditions affect indoor microbial communities. Existing studies have mostly focussed on the local (e.g., building units) or continental scale and rarely on the regional scale, e.g. a specific metropolitan area. Therefore, we wanted to identify key environmental determinants for the house dust microbiome from an existing collection of spatially (area of Munich, Germany) and temporally (301 days) distributed samples and to determine changes in the community as a function of time. To that end, dust samples that had been collected once from the living room floors of 286 individual households, were profiled for fungal and bacterial community variation and diversity using microbial fingerprinting techniques. The profiles were tested for their association with occupant behaviour, building characteristics, outdoor pollution, vegetation, and urbanization. Our results showed that more environmental and particularly outdoor factors (vegetation, urbanization, airborne particulate matter) affected the community composition of indoor fungi than of bacteria. The passage of time affected fungi and, surprisingly, also strongly affected bacteria. We inferred that fungal communities in indoor dust changed semi-annually, whereas bacterial communities paralleled outdoor plant phenological periods. These differences in temporal dynamics cannot be fully explained and should be further investigated in future studies on indoor microbiomes.     

Examining the fish microbiome: vertebrate-derived bacteria as an environmental niche for the discovery of unique marine natural products

Historically, marine invertebrates have been a prolific source of unique natural products, with a diverse array of biological activities. Recent studies of invertebrate-associated microbial communities are revealing microorganisms as the true producers of many of these compounds. Inspired by the human microbiome project, which has highlighted the human intestine as a unique microenvironment in terms of microbial diversity, we elected to examine the bacterial communities of fish intestines (which we have termed the fish microbiome) as a new source of microbial and biosynthetic diversity for natural products discovery. To test the hypothesis that the fish microbiome contains microorganisms with unique capacity for biosynthesizing natural products, we examined six species of fish through a combination of dissection and culture-dependent evaluation of intestinal microbial communities. Using isolation media designed to enrich for marine Actinobacteria, we have found three main clades that show taxonomic divergence from known strains, several of which are previously uncultured. Extracts from these strains exhibit a wide range of activities against both gram-positive and gram-negative human pathogens, as well as several fish pathogens. Exploration of one of these extracts has identified the novel bioactive lipid sebastenoic acid as an anti-microbial agent, with activity against Staphylococcus aureus, Bacillus subtilis, Enterococcus faecium, and Vibrio mimicus.     

Rhizosphere microbiome functional diversity and pathogen invasion resistance build up during plant development

The rhizosphere microbiome is essential for plant growth and health, and numerous studies have attempted to link microbiome functionality to species and trait composition. However, to date little is known about the actual ecological processes shaping community composition, complicating attempts to steer microbiome functionality. Here, we assess the development of microbial life history and community-level species interaction patterns that emerge during plant development. We use microbial phenotyping to experimentally test the development of niche complementarity and life history traits linked to microbiome performance. We show that the rhizosphere microbiome assembles from pioneer assemblages of species with random resource overlap into high-density, functionally complementary climax communities at later stages. During plant growth, fast-growing species were further replaced by antagonistic and stress-tolerant ones. Using synthetic consortia isolated from different plant growth stages, we demonstrate that the high functional diversity of ‘climax’ microbiomes leads to a better resistance to bacterial pathogen invasion. By demonstrating that different life-history strategies prevail at different plant growth stages and that community-level processes may supersede the importance of single species, we provide a new toolbox to understand microbiome assembly and steer its functionality at a community level.     

The rice foot rot pathogen Dickeya zeae alters the in-field plant microbiome

Studies on bacterial plant diseases have thus far been focused on the single bacterial species causing the disease, with very little attention given to the many other microorganisms present in the microbiome. This study intends to use pathobiome analysis of the rice foot rot disease, caused by Dickeya zeae, as a case study to investigate the effects of this bacterial pathogen to the total resident microbiome and to highlight possible interactions between the pathogen and the members of the community involved in the disease process. The microbiome of asymptomatic and the pathobiome of foot-rot symptomatic field-grown rice plants over two growing periods and belonging to two rice cultivars were determined via 16S rRNA gene amplicon sequencing. Results showed that the presence of D. zeae is associated with an alteration of the resident bacterial community in terms of species composition, abundance and richness, leading to the formation of microbial consortia linked to the disease state. Several bacterial species were significantly co-presented with the pathogen in the two growing periods suggesting that they could be involved in the disease process. Besides, culture-dependent isolation and in planta inoculation studies of a bacterial member of the pathobiome, identified as positive correlated with the pathogen in our in silico analysis, indicated that it benefits from the presence of D. zeae. A similar microbiome/pathobiome experiment was also performed in a symptomatically different rice disease evidencing that not all plant diseases have the same consequence/relationship with the plant microbiome. This study moves away from a pathogen-focused stance and goes towards a more ecological perception considering the effect of the entire microbial community which could be involved in the pathogenesis, persistence, transmission and evolution of plant pathogens.     

Enrichment versus biofilm culture: a functional and phylogenetic comparison of polycyclic aromatic hydrocarbon-degrading microbial communities

The effect that culture methods have on the diversity of degradative microbial communities is not well understood. We compared conventional batch enrichment with a biofilm culture method for the isolation of polycyclic aromatic hydrocarbon (PAH)-degrading microbial communities from a PAH-contaminated soil. The two methods were assessed by comparing: (i) the diversity of culturable bacteria; (ii) the diversity of PAH-catabolic genes in isolated bacteria; (iii) the inter- and intraspecific diversity of active PAH-catabolic gene classes; (iv) the diversity of bacteria present in 16S rRNA gene libraries generated from RNA extracted from the two communities and soil; and (v) the estimated diversity of active bacteria in the soil and culture systems. Single-strand conformation polymorphism analysis showed that the biofilm culture yielded 36 bacterial and two fungal species compared with 12 bacterial species from the enrichment culture. Application of accumulation and non-parametric estimators to clone libraries generated from 16S rRNA confirmed that the biofilm community contained greater diversity. Sequencing of clones showed that only species from the Proteobacteria were active in the enrichment culture, and that these species were expressing an identical nahAc-like naphthalene dioxygenase. 16S rRNA clones generated from the biofilm community indicated that species from the Cytophaga/Flavobacterium, high G+C bacteria and Proteobacteria were active at the time of sampling, expressing cndA-, nahAc- and phnAc-like naphthalene dioxygenases. The diversity of active species in the biofilm culture system closely matched that in the PAH-contaminated source soil. The results of this study showed that biofilm culture methods are more appropriate for the study of community-level interactions in PAH-degrading microbial communities. The study also indicated that cultivation of microbial communities on solid media might be the primary source of bias in the recovery of diverse species.     

Metabolic reconstruction for metagenomic data and its application to the human microbiome

Microbial communities carry out the majority of the biochemical activity on the planet, and they play integral roles in processes including metabolism and immune homeostasis in the human microbiome. Shotgun sequencing of such communities' metagenomes provides information complementary to organismal abundances from taxonomic markers, but the resulting data typically comprise short reads from hundreds of different organisms and are at best challenging to assemble comparably to single-organism genomes. Here, we describe an alternative approach to infer the functional and metabolic potential of a microbial community metagenome. We determined the gene families and pathways present or absent within a community, as well as their relative abundances, directly from short sequence reads. We validated this methodology using a collection of synthetic metagenomes, recovering the presence and abundance both of large pathways and of small functional modules with high accuracy. We subsequently applied this method, HUMAnN, to the microbial communities of 649 metagenomes drawn from seven primary body sites on 102 individuals as part of the Human Microbiome Project (HMP). This provided a means to compare functional diversity and organismal ecology in the human microbiome, and we determined a core of 24 ubiquitously present modules. Core pathways were often implemented by different enzyme families within different body sites, and 168 functional modules and 196 metabolic pathways varied in metagenomic abundance specifically to one or more niches within the microbiome. These included glycosaminoglycan degradation in the gut, as well as phosphate and amino acid transport linked to host phenotype (vaginal pH) in the posterior fornix. An implementation of our methodology is available at http://huttenhower.sph.harvard.edu/humann. This provides a means to accurately and efficiently characterize microbial metabolic pathways and functional modules directly from high-throughput sequencing reads, enabling the determination of community roles in the HMP cohort and in future metagenomic studies.     

Spatiotemporal development of the bacterial community in a tubular longitudinal microbial fuel cell

The spatiotemporal development of a bacterial community in an exoelectrogenic biofilm was investigated in sucrose-fed longitudinal tubular microbial fuel cell reactors, consisting of two serially connected modules. The proportional changes in the microbial community composition were assessed by polymerase chain reaction–denaturing gradient gel electrophoresis (DGGE) and DNA sequencing in order to relate them to the performance and stability of the bioelectrochemical system. The reproducibility of duplicated reactors, evaluated by cluster analysis and Jaccard’s coefficient, shows 80–90% similarity in species composition. Biofilm development through fed-batch start-up and subsequent stable continuous operation results in a population shift from γ-Proteobacteria- and Bacteroidetes- to Firmicutes-dominated communities, with other diverse species present at much lower relative proportions. DGGE patterns were analysed by range-weighted richness (Rr) and Pareto–Lorenz evenness distribution curves to investigate the evolution of the bacterial community. The first modules shifted from dominance by species closely related to Bacteroides graminisolvens, Raoultella ornithinolytica and Klebsiella sp. BM21 at the start of continuous-mode operation to a community dominated by Paludibacter propionicigenes-, Lactococcus sp.-, Pantoea agglomerans- and Klebsiella oxytoca-related species with stable power generation (6.0 W/m³) at day 97. Operational strategies that consider the dynamics of the population will provide useful parameters for evaluating system performance in the practical application of microbial fuel cells.     

Impact of Substratum Surface on Microbial Community Structure and Treatment Performance in Biological Aerated Filters

The impact of substratum surface property change on biofilm community structure was investigated using laboratory biological aerated filter (BAF) reactors and molecular microbial community analysis. Two substratum surfaces that differed in surface properties were created via surface coating and used to develop biofilms in test (modified surface) and control (original surface) BAF reactors. Microbial community analysis by 16S rRNA gene-based PCR-denaturing gradient gel electrophoresis (DGGE) showed that the surface property change consistently resulted in distinct profiles of microbial populations during replicate reactor start-ups. Pyrosequencing of the bar-coded 16S rRNA gene amplicons surveyed more than 90% of the microbial diversity in the microbial communities and identified 72 unique bacterial species within 19 bacterial orders. Among the 19 orders of bacteria detected, Burkholderiales and Rhodocyclales of the Betaproteobacteria class were numerically dominant and accounted for 90.5 to 97.4% of the sequence reads, and their relative abundances in the test and control BAF reactors were different in consistent patterns during the two reactor start-ups. Three of the five dominant bacterial species also showed consistent relative abundance changes between the test and control BAF reactors. The different biofilm microbial communities led to different treatment efficiencies, with consistently higher total organic carbon (TOC) removal in the test reactor than in the control reactor. Further understanding of how surface properties affect biofilm microbial communities and functional performance would enable the rational design of new generations of substrata for the improvement of biofilm-based biological treatment processes.