Explore the vast range of titles available.

Release of carbon previously locked in permafrost is a potentially important positive climate feedback. Now metagenomics reveal the vulnerability of active-layer soil carbon to warming-induced microbial decomposition in Alaskan tundra. Microbial decomposition of soil carbon in high-latitude tundra underlain with permafrost is one of the most important, but poorly understood, potential positive feedbacks of greenhouse gas emissions from terrestrial ecosystems into the atmosphere in a warmer world 1 , 2 , 3 , 4 . Using integrated metagenomic technologies, we showed that the microbial functional community structure in the active layer of tundra soil was significantly altered after only 1.5 years of warming, a rapid response demonstrating the high sensitivity of this ecosystem to warming. The abundances of microbial functional genes involved in both aerobic and anaerobic carbon decomposition were also markedly increased by this short-term warming. Consistent with this, ecosystem respiration ( R eco ) increased up to 38%. In addition, warming enhanced genes involved in nutrient cycling, which very likely contributed to an observed increase (30%) in gross primary productivity (GPP). However, the GPP increase did not offset the extra R eco , resulting in significantly more net carbon loss in warmed plots compared with control plots. Altogether, our results demonstrate the vulnerability of active-layer soil carbon in this permafrost-based tundra ecosystem to climate warming and the importance of microbial communities in mediating such vulnerability.

Background Soil is an important reservoir of antibiotic resistance genes (ARGs), but their potential risk in different ecosystems as well as response to anthropogenic land use change is unknown. We used a metagenomic approach and datasets with well-characterized metadata to investigate ARG types and amounts in soil DNA of three native ecosystems: Alaskan tundra, US Midwestern prairie, and Amazon rainforest, as well as the effect of conversion of the latter two to agriculture and pasture, respectively. Results High diversity (242 ARG subtypes) and abundance (0.184–0.242 ARG copies per 16S rRNA gene copy) were observed irrespective of ecosystem, with multidrug resistance and efflux pump the dominant class and mechanism. Ten regulatory genes were identified and they accounted for 13–35% of resistome abundances in soils, among them arlR , cpxR , ompR , vanR , and vanS were dominant and observed in all studied soils. We identified 55 non-regulatory ARGs shared by all 26 soil metagenomes of the three ecosystems, which accounted for more than 81% of non-regulatory resistome abundance. Proteobacteria, Firmicutes, and Actinobacteria were primary ARG hosts, 7 of 10 most abundant ARGs were found in all of them. No significant differences in both ARG diversity and abundance were observed between native prairie soil and adjacent long-term cultivated agriculture soil. We chose 12 clinically important ARGs to evaluate at the sequence level and found them to be distinct from those in human pathogens, and when assembled they were even more dissimilar. Significant correlation was found between bacterial community structure and resistome profile, suggesting that variance in resistome profile was mainly driven by the bacterial community composition. Conclusions Our results identify candidate background ARGs (shared in all 26 soils), classify ARG hosts, quantify resistance classes, and provide quantitative and sequence information suggestive of very low risk but also revealing resistance gene variants that might emerge in the future. Aify4ac4FW_3-HKjdJXvGW Video abstract

Antibiotics have been administered to agricultural animals for disease treatment disease prevention, and growth promotion for over 50 y. The impact of such antibiotic use on the treatment of human diseases is hotly debated. We raised pigs in a highly controlled environment with one portion of the littermates receiving a diet containing performance-enhancing antibiotics [chlortetracycline, sulfamethazine, and penicillin (known as ASP250)] and the other portion receiving the same diet but without the antibiotics. We used phylogenetic, metagenomic, and quantitative PCR-based approaches to address the impact of antibiotics on the swine gut microbiota. Bacterial phylotypes shifted after 14 d of antibiotic treatment with the medicated pigs showing an increase in Proteobacteria (1-11%) compared with nonmedicated pigs at the same time point. This shift was driven by an increase in Escherichia coli populations. Analysis of the metagenomes showed that microbial functional genes relating to energy production and conversion were increased in the antibiotic-fed pigs. The results also indicate that antibiotic resistance genes increased in abundance and diversity in the medicated swine microbiome despite a high background of resistance genes in nonmedicated swine. Some enriched genes, such as aminoglycoside O-phosphotransferases, confer resistance to antibiotics that were not administered in this study, demonstrating the potential for indirect selection of resistance to classes of antibiotics not fed. The collateral effects of feeding subtherapeutic doses of antibiotics to agricultural animals are apparent and must be considered in cost-benefit analyses.

The high-throughput antibiotic resistance gene (ARG) qPCR array, initially published in 2012, is increasingly used to quantify resistance and mobile determinants in environmental matrices. Continued utility of the array; however, necessitates improvements such as removing or redesigning questionable primer sets, updating targeted genes and coverage of available sequences. Towards this goal, a new primer design tool (EcoFunPrimer) was used to aid in identification of conserved regions of diverse genes. The total number of assays used for diverse genes was reduced from 91 old primer sets to 52 new primer sets, with only a 10% loss in sequence coverage. While the old and new array both contain 384 primer sets, a reduction in old primer sets permitted 147 additional ARGs and mobile genetic elements to be targeted. Results of validating the updated array with a mock community of strains resulted in over 98% of tested instances incurring true positive/negative calls. Common queries related to sensitivity, quantification and conventional data analysis (e.g. Ct cutoff value, and estimated genomic copies without standard curves) were also explored. A combined list of new and previously used primer sets is provided with a recommended set based on redesign of primer sets and results of validation.

Soil microbial respiration is an important source of uncertainty in projecting future climate and carbon (C) cycle feedbacks. However, its feedbacks to climate warming and underlying microbial mechanisms are still poorly understood. Here we show that the temperature sensitivity of soil microbial respiration ( Q 10 ) in a temperate grassland ecosystem persistently decreases by 12.0 ± 3.7% across 7 years of warming. Also, the shifts of microbial communities play critical roles in regulating thermal adaptation of soil respiration. Incorporating microbial functional gene abundance data into a microbially-enabled ecosystem model significantly improves the modeling performance of soil microbial respiration by 5–19%, and reduces model parametric uncertainty by 55–71%. In addition, modeling analyses show that the microbial thermal adaptation can lead to considerably less heterotrophic respiration (11.6 ± 7.5%), and hence less soil C loss. If such microbially mediated dampening effects occur generally across different spatial and temporal scales, the potential positive feedback of soil microbial respiration in response to climate warming may be less than previously predicted. Mechanisms and consequences of the acclimation of soil respiration to warming are unclear. Here, the authors combine soil respiration, metagenomics, and functional gene results from a 7-year grassland warming experiment to a microbial-enzyme decomposition model, showing functional gene information to lower uncertainty and improve fit.

Increases in carbon (C) inputs to soil can replenish soil organic C (SOC) through various mechanisms. However, recent studies have suggested that the increased C input can also stimulate the decomposition of old SOC via priming. Whether the loss of old SOC by priming can override C replenishment has not been rigorously examined. Here we show, through data–model synthesis, that the magnitude of replenishment is greater than that of priming, resulting in a net increase in SOC by a mean of 32% of the added new C. The magnitude of the net increase in SOC is positively correlated with the nitrogen-to-C ratio of the added substrates. Additionally, model evaluation indicates that a two-pool interactive model is a parsimonious model to represent the SOC decomposition with priming and replenishment. Our findings suggest that increasing C input to soils likely promote SOC accumulation despite the enhanced decomposition of old C via priming. The magnitudes of replenishment and priming, two important but opposing fluxes in soil organic carbon (SOC) dynamics, have not been compared. Here the authors show that the magnitude of replenishment is greater than that of priming, resulting in a net SOC accumulation after additional carbon input to soils.

Antibiotic resistance is a worldwide health risk, but the influence of animal agriculture on the genetic context and enrichment of individual antibiotic resistance alleles remains unclear. Using quantitative PCR followed by amplicon sequencing, we quantified and sequenced 44 genes related to antibiotic resistance, mobile genetic elements, and bacterial phylogeny in microbiomes from U.S. laboratory swine and from swine farms from three Chinese regions. We identified highly abundant resistance clusters: groups of resistance and mobile genetic element alleles that cooccur. For example, the abundance of genes conferring resistance to six classes of antibiotics together with class 1 integrase and the abundance of IS 6100 -type transposons in three Chinese regions are directly correlated. These resistance cluster genes likely colocalize in microbial genomes in the farms. Resistance cluster alleles were dramatically enriched (up to 1 to 10% as abundant as 16S rRNA) and indicate that multidrug-resistant bacteria are likely the norm rather than an exception in these communities. This enrichment largely occurred independently of phylogenetic composition; thus, resistance clusters are likely present in many bacterial taxa. Furthermore, resistance clusters contain resistance genes that confer resistance to antibiotics independently of their particular use on the farms. Selection for these clusters is likely due to the use of only a subset of the broad range of chemicals to which the clusters confer resistance. The scale of animal agriculture and its wastes, the enrichment and horizontal gene transfer potential of the clusters, and the vicinity of large human populations suggest that managing this resistance reservoir is important for minimizing human risk. IMPORTANCE Agricultural antibiotic use results in clusters of cooccurring resistance genes that together confer resistance to multiple antibiotics. The use of a single antibiotic could select for an entire suite of resistance genes if they are genetically linked. No links to bacterial membership were observed for these clusters of resistance genes. These findings urge deeper understanding of colocalization of resistance genes and mobile genetic elements in resistance islands and their distribution throughout antibiotic-exposed microbiomes. As governments seek to combat the rise in antibiotic resistance, a balance is sought between ensuring proper animal health and welfare and preserving medically important antibiotics for therapeutic use. Metagenomic and genomic monitoring will be critical to determine if resistance genes can be reduced in animal microbiomes, or if these gene clusters will continue to be coselected by antibiotics not deemed medically important for human health but used for growth promotion or by medically important antibiotics used therapeutically. Agricultural antibiotic use results in clusters of cooccurring resistance genes that together confer resistance to multiple antibiotics. The use of a single antibiotic could select for an entire suite of resistance genes if they are genetically linked. No links to bacterial membership were observed for these clusters of resistance genes. These findings urge deeper understanding of colocalization of resistance genes and mobile genetic elements in resistance islands and their distribution throughout antibiotic-exposed microbiomes. As governments seek to combat the rise in antibiotic resistance, a balance is sought between ensuring proper animal health and welfare and preserving medically important antibiotics for therapeutic use. Metagenomic and genomic monitoring will be critical to determine if resistance genes can be reduced in animal microbiomes, or if these gene clusters will continue to be coselected by antibiotics not deemed medically important for human health but used for growth promotion or by medically important antibiotics used therapeutically.

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

Estimation of the coverage provided by a metagenomic data set, i.e., what fraction of the microbial community was sampled by DNA sequencing, represents an essential first step of every culture-independent genomic study that aims to robustly assess the sequence diversity present in a sample. However, estimation of coverage remains elusive because of several technical limitations associated with high computational requirements and limiting statistical approaches to quantify diversity. Here we described Nonpareil 3, a new bioinformatics algorithm that circumvents several of these limitations and thus can facilitate culture-independent studies in clinical or environmental settings, independent of the sequencing platform employed. In addition, we present a new metric of sequence diversity based on rarefied coverage and demonstrate its use in communities from diverse ecosystems. Estimations of microbial community diversity based on metagenomic data sets are affected, often to an unknown degree, by biases derived from insufficient coverage and reference database-dependent estimations of diversity. For instance, the completeness of reference databases cannot be generally estimated since it depends on the extant diversity sampled to date, which, with the exception of a few habitats such as the human gut, remains severely undersampled. Further, estimation of the degree of coverage of a microbial community by a metagenomic data set is prohibitively time-consuming for large data sets, and coverage values may not be directly comparable between data sets obtained with different sequencing technologies. Here, we extend Nonpareil, a database-independent tool for the estimation of coverage in metagenomic data sets, to a high-performance computing implementation that scales up to hundreds of cores and includes, in addition, a k -mer-based estimation as sensitive as the original alignment-based version but about three hundred times as fast. Further, we propose a metric of sequence diversity ( N d ) derived directly from Nonpareil curves that correlates well with alpha diversity assessed by traditional metrics. We use this metric in different experiments demonstrating the correlation with the Shannon index estimated on 16S rRNA gene profiles and show that N d additionally reveals seasonal patterns in marine samples that are not captured by the Shannon index and more precise rankings of the magnitude of diversity of microbial communities in different habitats. Therefore, the new version of Nonpareil, called Nonpareil 3, advances the toolbox for metagenomic analyses of microbiomes. IMPORTANCE Estimation of the coverage provided by a metagenomic data set, i.e., what fraction of the microbial community was sampled by DNA sequencing, represents an essential first step of every culture-independent genomic study that aims to robustly assess the sequence diversity present in a sample. However, estimation of coverage remains elusive because of several technical limitations associated with high computational requirements and limiting statistical approaches to quantify diversity. Here we described Nonpareil 3, a new bioinformatics algorithm that circumvents several of these limitations and thus can facilitate culture-independent studies in clinical or environmental settings, independent of the sequencing platform employed. In addition, we present a new metric of sequence diversity based on rarefied coverage and demonstrate its use in communities from diverse ecosystems.

Currently, the naïve Bayesian classifier provided by the Ribosomal Database Project (RDP) is one of the most widely used tools to classify 16S rRNA sequences, mainly collected from environmental samples. We show that RDP has 97+% assignment accuracy and is fast for 250 bp and longer reads when the read originates from a taxon known to the database. Because most environmental samples will contain organisms from taxa whose 16S rRNA genes have not been previously sequenced, we aim to benchmark how well the RDP classifier and other competing methods can discriminate these novel taxa from known taxa. Because each fragment is assigned a score (containing likelihood or confidence information such as the boostrap score in the RDP classifier), we \"train\" a threshold to discriminate between novel and known organisms and observe its performance on a test set. The threshold that we determine tends to be conservative (low sensitivity but high specificity) for naïve Bayesian methods. Nonetheless, our method performs better with the RDP classifier than the other methods tested, measured by the f-measure and the area-under-the-curve on the receiver operating characteristic of the test set. By constraining the database to well-represented genera, sensitivity improves 3-15%. Finally, we show that the detector is a good predictor to determine novel abundant taxa (especially for finer levels of taxonomy where novelty is more likely to be present). We conclude that selecting a read-length appropriate RDP bootstrap score can significantly reduce the search space for identifying novel genera and higher levels in taxonomy. In addition, having a well-represented database significantly improves performance while having genera that are \"highly\" similar does not make a significant improvement. On a real dataset from an Amazon Terra Preta soil sample, we show that the detector can predict (or correlates to) whether novel sequences will be assigned to new taxa when the RDP database \"doubles\" in the future.