Search

Multiple fungal species, including potential opportunistic pathogens have been previously identified in water systems. Here, we investigated over 250 restroom sink fungal communities across a university campus and evaluated their diversity and core taxa present. Remarkable similarity in mycobial community composition was observed across buildings with Ascomycota consistently dominating. We found a core mycobiome independent of the building sampled, that included Exophiala species, potential opportunistic pathogenic black yeasts. Other prevalent and dominant taxa included Saccharomyces and Fusarium, common built environment fungi. The frequent presence of Malassezia, a common skin commensal, showed the external influence of human activities as a source of fungi to sinks. The study represents a novel exploration of sink P-traps mycobial communities from a public area and highlights their importance as reservoirs of possible pathogenic fungi, as well as emphasizing the relevance of further research in this understudied ecosystem within the built environment.

This work aimed to characterize the gut and faeces bacterial communities (BC) of Porcellionides pruinosus using high-throughput sequencing. Isopods were collected from the field and kept in laboratory conditions similar to those normally applied in ecotoxicology tests. Faeces and purged guts of isopods (n = 3 × 30) were analysed by pyrosequencing the V3-V4 region of 16 S rRNA encoding gene. Results showed that gut and faecal BCs were dominated by Proteobacteria, particularly by an OTU (Operational Taxonomic Unit) affiliated to genus Coxiella. Diversity and richness values were statistically higher for faecal BC, mainly due to the occurrence of several low-abundance phylotypes. These results may reflect faecal carriage of bacterial groups that cannot settle in the gut. BCs of P. pruinosus comprised: (1) common members of the soil microbiota, (2) bacterial symbionts, (3) bacteria related to host metabolic/ecological features, and (4) bacterial etiological agents. Comparison of BC of this isopod species with the BC from other invertebrates revealed common bacterial groups across taxa. The baseline information provided by this work will assist the design and data interpretation of future ecotoxicological or biomonitoring assays where the analysis of P. pruinosus BC should be included as an additional indicator. Terrestrial isopods bacterial communities might support ecotoxicological assays and biomonitoring processes as a valuable tool.

Soils are subjected to multiple anthropogenic modifications, but the synergistic impacts of simultaneous environmental stressors on below-ground communities are poorly understood. We used a large-scale (1152 plots), long-term (26 years), multi-factorial grassland experiment to assess the impact of five common agricultural practises (pesticides, herbicide, liming, fertilizers and grazing exclusion) and their interactive effects on the composition and activity of soil microbial communities. We confirmed that pH strongly impacts belowground communities, but further demonstrate that pH strongly mediates the impacts of other management factors. Notably, there was a significant interaction between liming and the effect of pesticide application, with only half of the taxa responding to pesticide being shared in both limed and unlimed treatments. Likewise, nutrient amendments significantly altered bacterial community structure in acidic soils. Not only do these results highlight an hierarchy of effect of commonly used agricultural practices but also the widespread interactions between treatments: many taxa were significantly affected by interactions between treatments, even in the absence of significant main effects. Furthermore, the results demonstrated that chemical amendments may not percolate deeply into physically unperturbed soils with effects concentrated between 0 and 30 cm, despite 20+ years of treatment. The research shows that future changes to agricultural practices will need to consider interactions among multiple factors.

Microplastics (MP) are emerging contaminants with the capacity to bind and transport hydrophobic organic compounds of environmental concern, such as polybrominated diphenyl ethers (PBDEs). The aim of this study was to investigate the ingestion of nylon (polyamide) MP alone and when associated with PBDEs and their effects on Chironomus sancticaroli larvae survival and microbiome structure. Survival, PBDE uptake and microbial community composition were measured in fourth instar larvae exposed for 96 h to BDEs- 47, 99, 100 and 153 in the presence and absence of 1% w/w MP in sediment. Microbiome community structures were determined through high throughput sequencing of 16S small subunit ribosomal RNA gene (16S rRNA). Initial experiments showed that larvae ingested MP faster at 0.5% w/w MP, while depuration was more efficient at 1% w/w MP, although retention of MP was seen even after 168 h depuration. No mortality was observed as a result of PBDEs and MP exposure. MP had a negative effect on PBDE concentration within larvae (η2 = 0.94) and a negative effect on sediment concentrations (η2 = 0.48). In all samples, microbial communities were dominated by Alphaproteobacteria, Betaproteobacteria, Actinobacteria and Gammaproteobacteria. Bacterial alpha diversity was not significantly affected by PBDEs or MP exposure. However, the abundance of discrete bacterial taxa was more sensitive to MP (X2 = 45.81, p = 0.02), and PBDE exposure. Our results highlight that C. sancticaroli showed no acute response to MPs and PBDEs, but that MPs influenced bacterial microbiome structure even after only short-term acute exposure.

Patterns of species diversity provide fundamental insights into the underlying mechanisms and processes that regulate biodiversity. The species–time relationship (STR) has the potential to be one such pattern; in a comparable manner to its more extensively studied spatial analogue, the species–area relationship (SAR), which has been pivotal in the development of ecological models and theories. We sought to determine the mechanisms and processes that underpin STR patterns of temporal turnover by sampling bacterial communities within ten water-filled tree-holes on the same European beech tree through the course of a year. We took this natural model system to represent an archipelago of islands of varying sizes and with shared common immigration sources. We observed an inverse relationship between STR-derived turnover rates and island size. Further, turnover was related to island size and not island isolation within the study system as indicated by a low frequency of dispersal limitation and high homogenizing dispersal. Compared to SARs, STRs are understudied, as such, the findings from the current study should provide a renewed interest in STR-based patterns and processes.

Summary The heterogeneous nature of lotic habitats plays an important role in the complex ecological and evolutionary processes that structure the microbial communities within them. Due to such complexity, our understanding of lotic microbial ecology still lacks conceptual frameworks for the ecological processes that shape these communities. We explored how bacterial community composition and underlying ecological assembly processes differ between lotic habitats by examining community composition and inferring community assembly processes across four major habitat types (free‐living, particle‐associated, biofilm on benthic stones and rocks, and sediment). This was conducted at 12 river sites from headwater streams to the main river in the River Thames, UK. Our results indicate that there are distinct differences in the bacterial communities between four major habitat types, with contrasting ecological processes shaping their community assembly processes. While the mobile free‐living and particle‐associated communities were consistently less diverse than the fixed sediment and biofilm communities, the latter two communities displayed higher homogeneity across the sampling sites. This indicates that the relative influence of deterministic environmental filtering is elevated in sediment and biofilm communities compared with free‐living and particle‐associated communities, where stochastic processes play a larger role.

International audience; Aims:Heterogeneous canopies in temperate montane forests affect microclimate and soil characteristics, with important effects on soil microbial communities and related processes. Here, we studied the interactions between plant root traits and soil bacterial and fungal communities in closed forest and open gaps in a mixed forest along an elevational gradient in the French Alps (1400, 1700 and 2000 m).Methods: Samples were separated into three fractions (plant root, rhizosphere and bulk soil), to further investigate the influence of plant zones on microbial communities. Bacterial (16S) and fungal (ITS) biodiversity was determined using high throughput sequencing, along with standard measures of soil, litter and root traits.Results: We found that (i) microbial community diversity was higher in gaps than in closed forest because of increased root trait diversity and density; (ii) open versus closed forest patches affected phylogenetic dispersion despite differences in elevations with phylogenetic clustering in closed forest; (iii) the interaction between root traits and microbial communities was stronger for rhizosphere and root fractions than for bulk soil and (iv) bacterial community composition was better explained by root traits than for fungi.Conclusions: Our findings highlight the importance of open gaps versus closed forest patches and associated root traits affecting microbial community structure, particularly for bacterial assemblages that exhibited a stronger interaction with root traits than for fungi

Background Antimicrobial resistance (AMR) in Enterobacterales is a global health threat. Capacity for individual-level surveillance remains limited in many countries, whilst population-level surveillance approaches could inform empiric antibiotic treatment guidelines. Methods In this exploratory study, a novel approach to population-level prediction of AMR in Enterobacterales clinical isolates using metagenomic (Illumina) profiling of pooled DNA extracts from human faecal samples was developed and tested. Taxonomic and AMR gene profiles were used to derive taxonomy-adjusted population-level AMR metrics. Bayesian modelling, and model comparison based on cross-validation, were used to evaluate the capacity of each metric to predict the number of resistant Enterobacterales invasive infections at a population-level, using available bloodstream/cerebrospinal fluid infection data. Findings Population metagenomes comprised samples from 177, 157, and 156 individuals in Kenya, the UK, and Cambodia, respectively, collected between September 2014 and April 2016. Clinical data from independent populations included 910, 3356 and 197 bacterial isolates from blood/cerebrospinal fluid infections in Kenya, the UK and Cambodia, respectively (samples collected between January 2010 and May 2017). Enterobacterales were common colonisers and pathogens, and faecal taxonomic/AMR gene distributions and proportions of antimicrobial-resistant Enterobacterales infections differed by setting. A model including terms reflecting the metagenomic abundance of the commonest clinical Enterobacterales species, and of AMR genes known to either increase the minimum inhibitory concentration (MIC) or confer clinically-relevant resistance, had a higher predictive performance in determining population-level resistance in clinical Enterobacterales isolates compared to models considering only AMR gene information, only taxonomic information, or an intercept-only baseline model (difference in expected log predictive density compared to best model, estimated using leave-one-out cross-validation: intercept-only model = -223 [95% credible interval (CI): -330,-116]; model considering only AMR gene information = -186 [95% CI: -281,-91]; model considering only taxonomic information = -151 [95% CI: -232,-69]). Interpretation Whilst our findings are exploratory and require validation, intermittent metagenomics of pooled samples could represent an effective approach for AMR surveillance and to predict population-level AMR in clinical isolates, complementary to ongoing development of laboratory infrastructures processing individual samples. Funding The study was funded by Bill & Melinda Gates Foundation (grant agreement OPP1160974) and was sponsored by University of Oxford. The study was also supported by the National Institute for Health Research (NIHR) Oxford Biomedical Research Centre, and NIHR Health Protection Research Unit in Healthcare-associated Infections and Antimicrobial Resistance (a partnership between the University of Oxford and Public Health England [PHE]). Kenyan samples were collected in a study funded by the MRC/DfID/Wellcome Joint Clinical Trials scheme: MR/M007367/1.

Risks posed to bees from neonicotinoid seed treatments (clothianidin, thiamethoxam, imidacloprid) led in 2013 to the European Union instigating a moratorium for their use on mass-flowering crops, including oilseed rape in the UK. This restriction did allow for the continued use of these seed treatments, in particular clothianidin, on non-flowering crops like winter wheat. To determine the impacts of the moratorium, we assessed neonicotinoid concentrations pre- (2014) and post- (2015−17) moratorium in 347 honey samples collected across Great Britain. While the probability of detecting clothianidin declined immediately following the moratorium, detection rates remained constant over the following three years (mean = 0.10 ppb, maximum = 2.8 ppb). In contrast, after three years thiamethoxam residues entirely disappeared while detection of imidacloprid was infrequent but persistent over the whole period. For those hives where neonicotinoids were detected, there was no evidence that the concentrations in the honey declined over the three years following the ban. Using metabarcoding approaches, we identified plants foraged upon by honeybees during the production of honey. After the moratorium came into effect, the highest neonicotinoid residues were associated with honey produced by foraging on both oilseed rape and several wild plants found in arable field margins. Concerns about soil persistence and uptake by non-target flowering plants ultimately led to a full European Union ban in 2018. Our results suggest that before this full ban came into effect, the use of clothianidin on non-flowering crops maintained a low-level probability of encountering this neonicotinoid within honey. However, these concentrations were low and would have been unlikely to pose significant risks to honeybees.

F-type plasmids are diverse and of great clinical significance, often carrying genes conferring antimicrobial resistance (AMR) such as extended-spectrum β-lactamases, particularly in Enterobacterales. Organising this plasmid diversity is challenging, and current knowledge is largely based on plasmids from clinical settings. Here, we present a network community analysis of a large survey of F-type plasmids from environmental (influent, effluent and upstream/downstream waterways surrounding wastewater treatment works) and livestock settings. We use a tractable and scalable methodology to examine the relationship between plasmid metadata and network communities. This reveals how niche (sampling compartment and host genera) partition and shape plasmid diversity. We also perform pangenome-style analyses on network communities. We show that such communities define unique combinations of core genes, with limited overlap. Building plasmid phylogenies based on alignments of these core genes, we demonstrate that plasmid accessory function is closely linked to core gene content. Taken together, our results suggest that stable F-type plasmid backbone structures can persist in environmental settings while allowing dramatic variation in accessory gene content that may be linked to niche adaptation. The association of F-type plasmids with AMR may reflect their suitability for rapid niche adaptation.

Worldwide honeybees (Apis mellifera L.) are one of the most widely kept domesticated animals, supporting domestic and commercial livelihoods through the production of honey and wax, as well as in the delivery of pollination services to crops. Quantifying which plant species are foraged upon by honeybees provides insights into their nutritional status as well as patterns of landscape scale habitat utilization. Here we outline a rapid and reproducible methodology for identifying environmental DNA (eDNA) originating principally from pollen grains suspended within honey. The process is based on a DNA extraction incorporating vacuum filtration prior to universal eukaryotic internal transcribed spacer 2 region (ITS2) amplicon generation, sequencing and identification. To provide a pre-cursor to sequence phylotyping, we outline systems for error-corrected processing amplicon sequence variant abundance tables that removes chimeras. This methodology underpins the new UK National Honey Monitoring Scheme.•We compare the efficacy and speed of centrifugation and filtration systems for removing pollen from honey samples as a precursor to plant DNA barcoding.•We introduce the ‘HONEYPI’ informatics pipeline, an open access resource implemented in python 2.7, to ensure long-term reproducibility during the process of amplicon sequence variant classification., Graphical abstract Image, graphical abstract

SUMMARYThe prevalence and proliferation of antimicrobial resistant bacteria is considered one of the critical issues of our time. Wastewater is a habitat for complex microbial communities where bacteria share genes of antimicrobial resistance through horizontal gene transfer. Hospital wastewater drainage systems are an ideal reservoir for environmental and pathogenic bacteria to interface and exchange genes of antimicrobial resistance1. Highly resistant Enterobacterales are of the greatest concern because they not only cause human infections but can also thrive in the environment2,3 and have shown to transmit from sink drains to patients4. Replacement of contaminated plumbing may be the most intuitive and widely deployed response once highly resistant potentially pathogenic bacteria are identified in a sink drain. Here we show using shotgun metagenomic sequencing of sink drain biofilms (in 6 intensive care patient rooms, 36 subsamples) paired with microbial culture analysis, an evident shift of biofilm community structure towards increased abundance of Enterobacteriaceae in the newly replaced sink plumbing. Significant increase in resistome load and abundance of clinically relevant resistance and typically mobile genes in the newly replaced plumbing was also observed. Findings from this study demonstrate that exchanging contaminated plumbing for new plumbing may actually have the unexpected consequence of increased abundance of Enterobacterales and genes of antimicrobial resistance in the sink drains. Disruption of preexisting complex environmental biofilms may result in an unintended increase in the amount of antimicrobial resistant Enterobacterales which are targeted for elimination. Temporal variability is an important attribute of ecosystems affecting colonization dynamics and biotic interactions (composition of introduced and resident communities) are key to functional outcomes. Absence of competition within an ecosystem favors invasion.

Escherichia coli and Klebsiella spp. are important human pathogens that cause a wide spectrum of clinical disease. In healthcare settings, sinks and other wastewater sites have been shown to be reservoirs of antimicrobial-resistant E. coli and Klebsiella spp., particularly in the context of outbreaks of resistant strains amongst patients. Without focusing exclusively on resistance markers or a clinical outbreak, we demonstrate that many hospital sink drains are abundantly and persistently colonized with diverse populations of E. coli , Klebsiella pneumoniae and Klebsiella oxytoca , including both antimicrobial-resistant and susceptible strains. Using whole-genome sequencing of 439 isolates, we show that environmental bacterial populations are largely structured by ward and sink, with only a handful of lineages, such as E. coli ST635, being widely distributed, suggesting different prevailing ecologies, which may vary as a result of different inputs and selection pressures. Whole-genome sequencing of 46 contemporaneous patient isolates identified one (2 %; 95 % CI 0.05–11 %) E. coli urine infection-associated isolate with high similarity to a prior sink isolate, suggesting that sinks may contribute to up to 10 % of infections caused by these organisms in patients on the ward over the same timeframe. Using metagenomics from 20 sink-timepoints, we show that sinks also harbour many clinically relevant antimicrobial resistance genes including bla CTX-M, bla SHV and mcr, and may act as niches for the exchange and amplification of these genes. Our study reinforces the potential role of sinks in contributing to Enterobacterales infection and antimicrobial resistance in hospital patients, something that could be amenable to intervention. This article contains data hosted by Microreact.

Antimicrobial resistance (AMR) is a global health threat, especially in low-/middle-income countries (LMICs), where there is limited surveillance to inform empiric antibiotic treatment guidelines. Enterobacterales are amongst the most important causes of drug-resistant bacterial infections. We developed a novel AMR surveillance approach for Enterobacterales by profiling pooled human faecal metagenomes from three sites (n=563 individuals; Cambodia, Kenya, UK) to derive a taxonomy-adjusted AMR metric (“resistance potential”) which could be used to predict the aggregate percentage of resistant invasive Enterobacterales infections within each setting. Samples were sequenced (Illumina); taxonomic and resistance gene profiling was performed using ResPipe. Data on organisms causing bacteraemia and meningitis and antibiotic susceptibility test results from 2010-2017 were collated for each site. Bayesian generalised linear models with a binomial likelihood were fitted to determine the capacity of the resistance potential to predict AMR in Enterobacterales infections in each setting. The most informative model accurately predicted the numbers of resistant infections in the target populations for 14/14 of antibiotics in the UK, 12/12 in Kenya, and 9/12 in Cambodia. Intermittent metagenomics of pooled human samples could represent a powerful pragmatic and economical approach for determining and monitoring AMR in clinical infections, especially in resource-limited settings.

Microplastics attract widespread attention, including for their potential to transport toxic chemicals in the form of plasticisers and associated hydrophobic organic chemicals, such as polybrominated diphenyl ethers (PBDEs). The aims of this study were to investigate how nylon (polyamide) microplastics may affect PBDE accumulation in snails, and the acute effects of nylon particles and PBDEs on survival, weight change and inherent microbiome diversity and community composition of the pond snail Lymnaea stagnalis. Snails were exposed for 96 h to BDEs-47, 99, 100 and 153 in the presence and absence of 1% w/w nylon microplastics in quartz sand sediment. No mortality was observed over the exposure period. Snails not exposed to microplastics lost significantly more weight compared to those exposed to microplastics. Increasing PBDE concentration in the sediment resulted in an increased PBDE body burden in the snails, however microplastics did not significantly influence total PBDE uptake. Based on individual congeners, uptake of BDE 47 by snails was significantly reduced in the presence of microplastics. The diversity and composition of the snail microbiome was not significantly altered by the presence of PBDEs nor by the microplastics, singly or combined. Significant effects on a few individual operational taxonomic units (OTUs) occurred when comparing the highest PBDE concentration with the control treatment, but in the absence of microplastics only. Overall within these acute experiments, only subtle effects on weight loss and slight microbiome alterations occurred. These results therefore highlight that L. stagnalis are resilient to acute exposures to microplastics and PBDEs, and that microplastics are unlikely to influence HOC accumulation or the microbiome of this species over short timescales.

High-throughput sequencing 16S rRNA gene surveys have enabled new insights into the diversity of soil bacteria, and furthered understanding of the ecological drivers of abundances across landscapes. However, current analytical approaches are of limited use in formalising syntheses of the ecological attributes of taxa discovered, because derived taxonomic units are typically unique to individual studies and sequence identification databases only characterise taxonomy. To address this, we used sequences obtained from a large nationwide soil survey (GB Countryside Survey, henceforth “CS”) to create a comprehensive soil specific 16S reference database, with coupled ecological information derived from the survey metadata. Specifically, we modelled taxon responses to soil pH at the OTU level using hierarchical logistic regression (HOF) models, to provide information on putative landscape scale pH-abundance responses. We identify that most of the soil OTUs examined exhibit predictable abundance responses across soil pH gradients, though with the exception of known acidophilic lineages, the pH optima of OTU relative abundance was variable and could not be generalised by broad taxonomy. This highlights the need for tools and databases to predict ecological traits at finer taxonomic resolution. We further demonstrate the utility of the database by testing against geographically dispersed query 16S datasets; evaluating efficacy by quantifying matches, and accuracy in predicting pH responses of query sequences from a separate large soil survey. We found that the CS database provided good coverage of dominant taxa; and that the taxa indicating soil pH in a query dataset corresponded with the pH classifications of top matches in the CS database. Furthermore we were able to predict query dataset community structure, using predicted abundances of dominant taxa based on query soil pH data and the HOF models of matched CS database taxa. The database with associated HOF model outputs is released as an online portal for querying single sequences of interest (https://shiny-apps.ceh.ac.uk/ID-TaxER), and as a DADA2 database for use in bioinformatics pipelines. The further development of advanced informatics infrastructures incorporating modelled ecological attributes along with new functional genomic information will likely facilitate large scale exploration and prediction of soil microbial functional biodiversity under current and future environmental change scenarios.

Antimicrobial resistance has been recognized as a threat to human health. The role of hospital sinks acting as a reservoir for some of the most concerning antibiotic resistant organisms, carbapenemase producing Enterobacterales (CPE) is evident but not well understood. Strategies to prevent establishment, interventions to eliminate these reservoirs and factors which drive persistence of CPE are not well established. We use a uniquely designed sink lab to transplant CPE colonized hospital sink plumbing with an aim to understand CPE dynamics in a controlled setting, notably exploiting both molecular and culture techniques. After ex situ installation the CPE population in the sink plumbing drop from previously detectable to undetectable levels. The addition of nutrients is followed by a quick rebound in CPE detection in the sinks after as many as 37 days. We did not however detect a significant shift in microbial community structure or the overall resistance gene carriage in longitudinal samples from a subset of these transplanted sinks using whole shotgun metagenomic sequencing. Comparing nutrient types in a benchtop culture study model, protein rich nutrients appear to be the most supportive for CPE growth and biofilm formation ability. The role of nutrients exposure is determining factor for maintaining a high bioburden of CPE in the sink drains and P-traps. Therefore, limiting nutrient disposal into sinks has reasonable potential with regard to decreasing the CPE wastewater burden, especially in hospitals seeking to control an environmental reservoir.

Background Shotgun metagenomics is increasingly used to characterise microbial communities, particularly for the investigation of antimicrobial resistance (AMR) in different animal and environmental contexts. There are many different approaches for inferring the taxonomic composition and AMR gene content of complex community samples from shotgun metagenomic data, but there has been little work establishing the optimum sequencing depth, data processing and analysis methods for these samples. In this study we used shotgun metagenomics and sequencing of cultured isolates from the same samples to address these issues. We sampled three potential environmental AMR gene reservoirs (pig caeca, river sediment, effluent) and sequenced samples with shotgun metagenomics at high depth (~ 200 million reads per sample). Alongside this, we cultured single-colony isolates of Enterobacteriaceae from the same samples and used hybrid sequencing (short- and long-reads) to create high-quality assemblies for comparison to the metagenomic data. To automate data processing, we developed an open-source software pipeline, ‘ResPipe’. Results Taxonomic profiling was much more stable to sequencing depth than AMR gene content. 1 million reads per sample was sufficient to achieve Thermus thermophilus DNA substantially changed the estimated gene abundance distributions. While the majority of genomic content from cultured isolates from effluent was recoverable using shotgun metagenomics, this was not the case for pig caeca or river sediment. Conclusions Sequencing depth and profiling method can critically affect the profiling of polymicrobial animal and environmental samples with shotgun metagenomics. Both sequencing of cultured isolates and shotgun metagenomics can recover substantial diversity that is not identified using the other methods. Particular consideration is required when inferring AMR gene content or presence by mapping metagenomic reads to a database. ResPipe, the open-source software pipeline we have developed, is freely available (https://gitlab.com/hsgweon/ResPipe).

Soil aggregates are critical to soil functionality, but there remain many uncertainties with respect to the role of biotic factors in forming aggregates. Understanding the interacting effects of soil, land use type, vegetation and microbial communities is a major challenge that needs assessment in both field and controlled laboratory conditions, as well as in bulk and rhizosphere soils. To address these effects and their feedbacks, we first examined the influence of soil, root and litter characteristics along a land use gradient (ancient woodland, secondary woodland, grassland, pasture and arable land) on microbial community structure (in both bulk and rhizosphere soil), as well as on aggregate stability. Then, we performed an inoculation experiment where we extracted soil columns from the arable and secondary woodland and used a third unstructured loamy soil as a control. We sterilized these three soils to remove microbial communities, and then inoculated the tops of sterilized soil columns with soil from the secondary woodland or the arable field sites. Control columns of all soil types were not inoculated. In a fully-crossed design, we planted two species possessing distinct root system morphologies: Brachypodium sylvaticum (fibrous system with many thin and fine roots) and Urtica dioica (taproot system with few fine roots). After four months, microbial communities (in bulk and rhizospheric soil) and aggregate stability were measured, along with root traits. In both the field and laboratory experiments, bacterial (16S) and fungal (ITS) diversity was determined using high throughput sequencing. In the field study we found that: i) there were strong relationships between aggregate stability and microbial community composition that were driven by land use, ii) the relationship between aggregate stability along the land use gradient and the trophic nature of bacterial communities was not significant, but certain soil, root and litter parameters shaped bacterial phyla, with oligotrophic bacteria conditioned by the rhizosphere niche, and copiotrophic phyla more dependent on bulk soil conditions, iii) land use gradient (from woodland to arable), reduced the relative abundance of saprotrophic and ectomycorrhizal fungi with an increase in the relative abundance of Ascomycota and a reduction in the relative abundance of Basidiomycota. In the laboratory experiment we found that: i) the inoculation of sterilized soils with soils from the field significantly increased aggregate stability in control soil that was initially poorly structured, ii) the effects of inoculation on aggregate stability were similar when either secondary woodland or arable soils were used as inoculums and iii) these effects were impacted significantly by root length density. Our results show that microbial communities influence soil structure and that bacterial communities are intimately associated to rhizospheric conditions and root traits (of which root length density was the most pertinent).

Soil microorganisms act as gatekeepers for soil–atmosphere carbon exchange by balancing the accumulation and release of soil organic matter. However, poor understanding of the mechanisms responsible hinders the development of effective land management strategies to enhance soil carbon storage. Here we empirically test the link between microbial ecophysiological traits and topsoil carbon content across geographically distributed soils and land use contrasts. We discovered distinct pH controls on microbial mechanisms of carbon accumulation. Land use intensification in low-pH soils that increased the pH above a threshold (~6.2) leads to carbon loss through increased decomposition, following alleviation of acid retardation of microbial growth. However, loss of carbon with intensification in near-neutral pH soils was linked to decreased microbial biomass and reduced growth efficiency that was, in turn, related to trade-offs with stress alleviation and resource acquisition. Thus, less-intensive management practices in near-neutral pH soils have more potential for carbon storage through increased microbial growth efficiency, whereas in acidic soils, microbial growth is a bigger constraint on decomposition rates., Land use intensification could modify microbial activity and thus ecosystem function. Here, Malik et al. sample microbes and carbon-related functions across a land use gradient, demonstrating that microbial biomass and carbon use efficiency are reduced in human-impacted near-neutral pH soils.

Explaining the large-scale diversity of soil organisms that drive biogeochemical processes—and their responses to environmental change—is critical. However, identifying consistent drivers of belowground diversity and abundance for some soil organisms at large spatial scales remains problematic. Here we investigate a major guild, the ectomycorrhizal fungi, across European forests at a spatial scale and resolution that is—to our knowledge—unprecedented, to explore key biotic and abiotic predictors of ectomycorrhizal diversity and to identify dominant responses and thresholds for change across complex environmental gradients. We show the effect of 38 host, environment, climate and geographical variables on ectomycorrhizal diversity, and define thresholds of community change for key variables. We quantify host specificity and reveal plasticity in functional traits involved in soil foraging across gradients. We conclude that environmental and host factors explain most of the variation in ectomycorrhizal diversity, that the environmental thresholds used as major ecosystem assessment tools need adjustment and that the importance of belowground specificity and plasticity has previously been underappreciated.

Bumblebees are ecologically and economically important as pollinators of crop and wild plants, especially in temperate systems. Species, such as the buff-tailed bumblebee (Bombus terrestris), are reared commercially to pollinate high value crops. Their highly specific gut microbiota, characterised by low diversity, may affect nutrition and immunity and are likely to be important for fitness and colony health. However, little is known about how environmental factors affect bacterial community structure. We analyzed the gut microbiota from three groups of worker bumblebees (B. terrestris) from distinct colonies that varied in rearing and foraging characteristics: commercially reared with restricted foraging (RR); commercially reared with outside foraging (RF); and wild-caught workers (W). Contrary to previous studies, which indicate that bacterial communities are highly conserved across workers, we found that RF individuals had an intermediate community structure compared to RR and W types. Further, this was shaped by differences in the abundances of common OTUs and the diversity of rare OTUs present which we propose results from an increase in the variety of carbohydrates obtained through foraging.

Zinc oxide nanoparticles (ZnO NPs) are used in an array of products and processes, ranging from personal care products to antifouling paints, textiles, food additives, antibacterial agents and environmental remediation processes. Soils are an environment likely to be exposed to manmade nanoparticles due to the practice of applying sewage sludge as a fertiliser or as an organic soil improver. However, understanding on the interactions between soil properties, nanoparticles and the organisms that live within soil is lacking, especially with regards to soil bacterial communities. We studied the effects of nanoparticulate, non-nanoparticulate and ionic zinc (in the form of zinc chloride) on the composition of bacterial communities in soil with a modified pH range (from pH 4.5 to pH 7.2). We observed strong pH-dependent effects on the interaction between bacterial communities and all forms of zinc, with the largest changes in bacterial community composition occurring in soils with low and medium pH levels (pH 4.8 and 5.9). The high pH soil (pH 7.2) was less susceptible to the effects of zinc exposure. At the highest doses of zinc (2500 mg/kg dw soil), both nano and non-nano particulate zinc applications elicited a similar response in the soil bacterial community, and this differed significantly to the ionic zinc salt treatment. The results highlight the importance of considering soil pH in nanotoxicology studies, although further work is needed to determine the exact mechanisms controlling the toxicity and fate and interactions of nanoparticles with soil microbial communities.

Characterization of microbial communities in stressful conditions at a field level is rather scarce, especially when considering fungal communities from aboveground habitats. We aimed at characterizing fungal communities from different poplar habitats at a Hg-contaminated phytomanagement site by using Illumina-based sequencing, network analysis approach, and direct isolation of Hg-resistant fungal strains. The highest diversity estimated by the Shannon index was found for soil communities, which was negatively affected by soil Hg concentration. Among the significant correlations between soil operational taxonomic units (OTUs) in the co-occurrence network, 80% were negatively correlated revealing dominance of a pattern of mutual exclusion. The fungal communities associated with Populus roots mostly consisted of OTUs from the symbiotic guild, such as members of the Thelephoraceae, thus explaining the lowest diversity found for root communities. Additionally, root communities showed the highest network connectivity index, while rarely detected OTUs from the Glomeromycetes may have a central role in the root network. Unexpectedly high richness and diversity were found for aboveground habitats, compared to the root habitat. The aboveground habitats were dominated by yeasts from the Lalaria, Davidiella, and Bensingtonia genera, not detected in belowground habitats. Leaf and stem habitats were characterized by few dominant OTUs such as those from the Dothideomycete class producing mutual exclusion with other OTUs. Aureobasidium pullulans, one of the dominating OTUs, was further isolated from the leaf habitat, in addition to Nakazawaea populi species, which were found to be Hg resistant. Altogether, these findings will provide an improved point of reference for microbial research on inoculation-based programs of tailings dumps.

Bacterial genome sizes have previously been shown to exhibit a bimodal distribution. This phenomenon has prompted discussion regarding the evolutionary forces driving genome size in bacteria and its ecological significance. We investigated the level of inherent redundancy in the public database and the effect it has on the shape of the apparent bimodal distribution. Our study reveals that there is a significant bias in the genome sequencing efforts towards a certain group of species, and that correcting the bias using species nomenclature and clustering of the 16S rRNA gene, results in a unimodal rather than the previously published bimodal distribution. The true genome size distribution and its wider ecological implications will soon emerge as we are currently witnessing rapid growth in the number of sequenced genomes from diverse environmental niches across a range of habitats at an unprecedented rate.

The complete mitochondrial genome of the European blackfly, Simulium variegatum Meigen, 1818 was sequenced using a combined Illumina and Sanger sequencing approach. Using the known sequence of Chironomus tepperi Skuse, 1889 (Chironomidae) homologous NGS reads were identified and assembled. The genome is 15,367 bp in length and includes 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes and a control region. Gene order resembles that of the ancestral dipteran gene arrangement. The base composition of the genome is A (37.6%), T (35.3%), C (15.8%) and G (11.3%). The control region between 12S rRNA and tRNAIle is composed of 362 bp with no obvious repetitive motifs.

International audience; Recent research has highlighted strong correlations between soil edaphic parameters and bacterial biodiversity. Here we seek to explore these relationships across the European Union member states with respect to mapping bacterial biodiversity at the continental scale. As part of the EU FP7 EcoFINDERs project, bacterial communities from 76 soil samples taken across Europe were assessed from eleven countries encompassing Arctic to Southern Mediterranean climes, representing a diverse range of soil types and land uses (grassland, forest and arable land). We found predictable relationships between community biodiversity (ordination site scores) and land use factors as well as soil properties such as pH. Based on the modelled relationship between soil pH and bacterial biodiversity found for the surveyed soils, we were able to predict biodiversity in 1000 soils for which soil pH data had been collected as part of national scale monitoring. We then performed interpolative mapping utilising existing EU wide soil pH data to present the first map of bacterial biodiversity across the EU member states. The predictive accuracy of the map was assessed again using the national scale data, but this time contrasting the EU wide spatial predictions with point data on bacterial communities. Generally the maps were useful at predicting broad extremes of biodiversity reflective of low or high pH soils, though predictive accuracy was limited for Britain particularly for organic/acidic soil communities. Spatial accuracy could however be increased by utilising published maps of soil pH calculated using geostatistical approaches at both global and national scales. These findings will contribute to wider efforts to predict and understand the spatial distribution of soil biodiversity at global scales. Further work should focus on enhancing the predictive power of such maps, by harmonising global datasets on soil conditioning parameters, soil properties and biodiversity; and the continued efforts to advance the geostatistical modelling of specific components of soil biodiversity at local to global scales.

Change history: In the HTML version of this Article, author 'Filipa Cox' had no affiliation in the author list, although she was correctly associated with affiliation 3 in the PDF. In addition, the blue circles for 'oak' were missing from Extended Data Fig. 1. These errors have been corrected online.

Summary Studying fungal biodiversity using data generated from Illumina MiSeq sequencing platforms poses a number of bioinformatic challenges with the analysis typically involving a large number of tools for each analytical step from quality filtering to generating identified operational taxonomic unit (OTU) abundance tables.Here, we introduce PIPITS, an open‐source stand‐alone suite of software for automated processing of Illumina MiSeq sequences for fungal community analysis. PIPITS exploits a number of state of the art applications to process paired‐end reads from quality filtering to producing OTU abundance tables.We provide detailed descriptions of the pipeline and show its utility in the analysis of 9 396 092 sequences generated on the MiSeq platform from Illumina MiSeq. PIPITS is the first automated bioinformatics pipeline dedicated for fungal ITS sequences which incorporates ITSx to extract subregions of ITS and exploits the latest RDP Classifier to classify sequences against the curated UNITE fungal data set.

Bumblebees are ecologically and economically important as pollinators of crop and wild plants, especially in temperate systems. Species, such as the buff-tailed bumblebee (Bombus terrestris), are reared commercially to pollinate high-value crops. Their highly specific gut microbiota, characterized by low diversity, may affect nutrition and immunity and are likely to be important for fitness and colony health. However, little is known about how environmental factors affect bacterial community structure. We analysed the gut microbiota from three groups of worker bumblebees (B. terrestris) from distinct colonies that varied in rearing and foraging characteristics: commercially reared with restricted foraging (RR); commercially reared with outside foraging (RF); and wild-caught workers (W). Contrary to previous studies, which indicate that bacterial communities are highly conserved across workers, we found that RF individuals had an intermediate community structure compared with RR and W types. Further, this was shaped by differences in the abundances of common operational taxonomic units (OTUs) and the diversity of rare OTUs present, which we propose results from an increase in the variety of carbohydrates obtained through foraging.

Lotic ecosystems such as rivers and streams are unique in that they represent a continuum of both space and time during the transition from headwaters to the river mouth. As microbes have very different controls over their ecology, distribution and dispersion compared with macrobiota, we wished to explore biogeographical patterns within a river catchment and uncover the major drivers structuring bacterioplankton communities. Water samples collected across the River Thames Basin, UK, covering the transition from headwater tributaries to the lower reaches of the main river channel were characterised using 16S rRNA gene pyrosequencing. This approach revealed an ecological succession in the bacterial community composition along the river continuum, moving from a community dominated by Bacteroidetes in the headwaters to Actinobacteria-dominated downstream. Location of the sampling point in the river network (measured as the cumulative water channel distance upstream) was found to be the most predictive spatial feature; inferring that ecological processes pertaining to temporal community succession are of prime importance in driving the assemblages of riverine bacterioplankton communities. A decrease in bacterial activity rates and an increase in the abundance of low nucleic acid bacteria relative to high nucleic acid bacteria were found to correspond with these downstream changes in community structure, suggesting corresponding functional changes. Our findings show that bacterial communities across the Thames basin exhibit an ecological succession along the river continuum, and that this is primarily driven by water residence time rather than the physico-chemical status of the river.

Because of the ever-increasing application of next-generation sequencing (NGS) in research, and the expectation of faster experiment turn-around, it is becoming unfeasible and unscalable for analysis to be done exclusively by existing trained bioinformaticians. Instead, researchers and bench biologists are performing at least parts of most analyses. In order for this to be realized, two conditions must be satisfied: (1) well designed and accessible tools need to be made available, and (2) researchers and biologists need to be trained to use such tools in order to confidently handle high volumes of NGS data. Bio-Linux is a fully featured, powerful, configurable and easy to maintain bioinformatics workstation and helps on both counts by offering well over one hundred bioinformatics tools packaged into a single distribution, easily accessible and readily usable. Bio- Linux is also accessible in the form of virtual images or on the cloud, thus providing researchers with immediate access to scalable compute infrastructure required to run the analysis. Furthermore this paper discusses how bioinformatics training on Bio-Linux is helping to bridge the data production and analysis gap.

Pantothenate (vitamin B5) is the precursor for the biosynthesis of the phosphopantetheine moiety of coenzyme A and acyl carrier protein, and is synthesised in Escherichia coli by four enzymic reactions. Ketopantoate hydroxymethyltransferase (KPHMT) and pantothenate synthetase (PtS) catalyse the first and last steps, respectively. Two genes encoding KPHMT and one for PtS were identified in the Arabidopsis thaliana genome, and cDNAs for all three genes were amplified by PCR. The cDNAs were able to complement their respective E. coli auxotrophs, demonstrating that they encoded functional enzymes. Subcellular localisation of the proteins was investigated using green fluorescent protein (GFP) fusions and confocal microscopy. The two KPHMT-GFP fusion proteins were targeted exclusively to mitochondria, whereas PtS-GFP was found in the cytosol. This implies that there must be transporters for pathway intermediates. KPHMT enzyme activity could be measured in purified mitochondria from both pea leaves and Arabidopsis suspension cultures. We investigated whether Arabidopsis encoded homologues of the remaining two pantothenate biosynthesis enzymes from E. coli, l-aspartate-alpha-decarboxylase (ADC) and ketopantoate reductase (KPR). No homologue of ADC could be identified using either conventional blast or searches with the program fugue in which the structure of the E. coli ADC was compared to all the annotated proteins in Arabidopsis. ADC also appears to be absent from the genome of the yeast, Saccharomyces cerevisiae, by the same criteria. In contrast, a putative Arabidopsis oxidoreductase with some similarity to KPR was identified with fugue.

Soil extracellular enzymes released by microorganisms break down organic matter and are crucial in regulating C, N and P cycling. Soil pH is known to influence enzyme activity, and is also a strong driver of microbial community composition; but little is known about how alterations in soil pH affect enzymatic activity and how this is mediated by microbial communities. To assess long term enzymatic adaptation to soil pH, we conducted enzyme assays at buffered pH levels (2.5 to 10, 0.5 interval) on two historically managed soils maintained at either pH 5 or 7 from the Rothamsted’s Park Grass Long-term experiment ). The pH optima for a range of enzymes was found to differ between the two soils, the direction of the shift being toward the source soil pH, indicating the production of pH adapted isoenzymes by the soil microbial community. Soil bacterial and fungal communities determined by amplicon sequencing were found to be clearly distinct between pH 5 and soil pH 7 soils, possibly explaining differences in enzymatic responses. Furthermore, β-glucosidase sequences extracted from metagenomes revealed an increased abundance of Acidobacteria in the pH 5 soils, and increased abundance of Actinobacteria in pH 7 soils; these taxonomic shifts were more pronounced for enzymatic sequences when compared with a number of housekeeping gene sequences. Particularly for the Acidobacteria, this indicates that broad taxonomic groups at phylum level may possess enzymatic adaptations which underpin competitiveness in different pH soils. More generally our findings have implications for modelling the efficiency of different microbial enzymatic processes under changing environmental conditions; and future work is required to identify trade-offs with pH adaptations, which could result in different activity responses to other environmental perturbations.