Long-read sequencing revealed cooccurrence, host range, and potential mobility of antibiotic resistome in cow feces.

While it is well recognized that the environmental resistome is global, diverse, and augmented by human activities, it has been difficult to assess risk because of the inability to culture many environmental organisms, and it is difficult to evaluate risk from current sequence-based environmental methods. The four most important criteria to determine risk are whether the antibiotic-resistance genes (ARGs) are a complete, potentially functional complement; if they are linked with other resistances; whether they are mobile; and the identity of their host. Long-read sequencing fills this important gap between culture and short sequence-based methods. To address these criteria, we collected feces from a ceftiofur-treated cow, enriched the samples in the presence of antibiotics to favor ARG functionality, and sequenced long reads using Nanopore and PacBio technologies. Multidrug-resistance genes comprised 58% of resistome abundance, but only 0.8% of them were plasmid associated; fluroquinolone-, aminoglycoside-, macrolide-lincosamide-streptogramin (MLS)-, and ß-lactam-resistance genes accounted for 2.7 to 12.3% of resistome abundance but with 19 to 78% located on plasmids. A variety of plasmid types were assembled, some of which share low similarity to plasmids in current databases. Enterobacteriaceae were dominant hosts of antibiotic-resistant plasmids; physical linkage of extended-spectrum ß-lactamase genes (CTX-M, TEM, CMY, and CARB) was largely found with aminoglycoside-, MLS-, tetracycline-, trimethoprim-, phenicol-, sulfonamide-, and mercury-resistance genes. A draft circular chromosome of Vagococcus lutrae was assembled; it carries MLS-, tetracycline- (including tetM and tetL on an integrative conjugative element), and trimethoprim-resistance genes flanked by many transposase genes and insertion sequences, implying that they remain transferrable.

Permafrost thaw with warming reduces microbial metabolic capacities in subsurface soils.

Microorganisms are major constituents of the total biomass in permafrost regions, whose underlain soils are frozen for at least two consecutive years. To understand potential microbial responses to climate change, here we examined microbial community compositions and functional capacities across four soil depths in an Alaska tundra site. We showed that a 5-year warming treatment increased soil thaw depth by 25.7% (p = .011) within the deep organic layer (15-25 cm). Concurrently, warming reduced 37% of bacterial abundance and 64% of fungal abundances in the deep organic layer, while it did not affect microbial abundance in other soil layers (i.e., 0-5, 5-15, and 45-55 cm). Warming treatment altered fungal community composition and microbial functional structure (p < .050), but not bacterial community composition. Using a functional gene array, we found that the relative abundances of a variety of carbon (C)-decomposing, iron-reducing, and sulphate-reducing genes in the deep organic layer were decreased, which was not observed by the shotgun sequencing-based metagenomics analysis of those samples. To explain the reduced metabolic capacities, we found that warming treatment elicited higher deterministic environmental filtering, which could be linked to water-saturated time, soil moisture, and soil thaw duration. In contrast, plant factors showed little influence on microbial communities in subsurface soils below 15 cm, despite a 25.2% higher (p < .05) aboveground plant biomass by warming treatment. Collectively, we demonstrate that microbial metabolic capacities in subsurface soils are reduced, probably arising from enhanced thaw by warming.

Winter warming rapidly increases carbon degradation capacities of fungal communities in tundra soil: Potential consequences on carbon stability.

High-latitude tundra ecosystems are increasingly affected by climate warming. As an important fraction of soil microorganisms, fungi play essential roles in carbon degradation, especially the old, chemically recalcitrant carbon. However, it remains obscure how fungi respond to climate warming and whether fungi, in turn, affect carbon stability of tundra. In a 2-year winter soil warming experiment of 2°C by snow fences, we investigated responses of fungal communities to warming in the active layer of an Alaskan tundra. Although fungal community composition, revealed by the 28S rRNA gene amplicon sequencing, remained unchanged (p > .05), fungal functional gene composition, revealed by a microarray named GeoChip, was altered (p < .05). Changes in functional gene composition were linked to winter soil temperature, thaw depth, soil moisture, and gross primary productivity (canonical correlation analysis, p < .05). Specifically, relative abundances of fungal genes encoding invertase, xylose reductase and vanillin dehydrogenase significantly increased (p < .05), indicating higher carbon degradation capacities of fungal communities under warming. Accordingly, we detected changes in fungal gene networks under warming, including higher average path distance, lower average clustering coefficient and lower percentage of negative links, indicating that warming potentially changed fungal interactions. Together, our study reveals higher carbon degradation capacities of fungal communities under short-term warming and highlights the potential impacts of fungal communities on tundra ecosystem respiration, and consequently future carbon stability of high-latitude tundra.

The Microbial Genomes Atlas (MiGA) webserver: taxonomic and gene diversity analysis of Archaea and Bacteria at the whole genome level.

The small subunit ribosomal RNA gene (16S rRNA) has been successfully used to catalogue and study the diversity of prokaryotic species and communities but it offers limited resolution at the species and finer levels, and cannot represent the whole-genome diversity and fluidity. To overcome these limitations, we introduced the Microbial Genomes Atlas (MiGA), a webserver that allows the classification of an unknown query genomic sequence, complete or partial, against all taxonomically classified taxa with available genome sequences, as well as comparisons to other related genomes including uncultivated ones, based on the genome-aggregate Average Nucleotide and Amino Acid Identity (ANI/AAI) concepts. MiGA integrates best practices in sequence quality trimming and assembly and allows input to be raw reads or assemblies from isolate genomes, single-cell sequences, and metagenome-assembled genomes (MAGs). Further, MiGA can take as input hundreds of closely related genomes of the same or closely related species (a so-called 'Clade Project') to assess their gene content diversity and evolutionary relationships, and calculate important clade properties such as the pangenome and core gene sets. Therefore, MiGA is expected to facilitate a range of genome-based taxonomic and diversity studies, and quality assessment across environmental and clinical settings. MiGA is available at http://microbial-genomes.org/.

ARGs-OAP v2.0 with an expanded SARG database and Hidden Markov Models for enhancement characterization and quantification of antibiotic resistance genes in environmental metagenomes.

Motivation: Much global attention has been paid to antibiotic resistance in monitoring its emergence, accumulation and dissemination. For rapid characterization and quantification of antibiotic resistance genes (ARGs) in metagenomic datasets, an online analysis pipeline, ARGs-OAP has been developed consisting of a database termed Structured Antibiotic Resistance Genes (the SARG) with a hierarchical structure (ARGs type-subtype-reference sequence). Results: The new release of the database, termed SARG version 2.0, contains sequences not only from CARD and ARDB databases, but also carefully selected and curated sequences from the latest protein collection of the NCBI-NR database, to keep up to date with the increasing number of ARG deposited sequences. SARG v2.0 has tripled the sequences of the first version and demonstrated improved coverage of ARGs detection in metagenomes from various environmental samples. In addition to annotation of high-throughput raw reads using a similarity search strategy, ARGs-OAP v2.0 now provides model-based identification of assembled sequences using SARGfam, a high-quality profile Hidden Markov Model (HMM), containing profiles of ARG subtypes. Additionally, ARGs-OAP v2.0 improves cell number quantification by using the average coverage of essential single copy marker genes, as an option in addition to the previous method based on the 16S rRNA gene. Availability and implementation: ARGs-OAP can be accessed through http://smile.hku.hk/SARGs. The database could be downloaded from the same site. Source codes for this study can be downloaded from https://github.com/xiaole99/ARGs-OAP-v2.0. Supplementary information: Supplementary data are available at Bioinformatics online.

RNAcentral: a comprehensive database of non-coding RNA sequences.

RNAcentral is a database of non-coding RNA (ncRNA) sequences that aggregates data from specialised ncRNA resources and provides a single entry point for accessing ncRNA sequences of all ncRNA types from all organisms. Since its launch in 2014, RNAcentral has integrated twelve new resources, taking the total number of collaborating database to 22, and began importing new types of data, such as modified nucleotides from MODOMICS and PDB. We created new species-specific identifiers that refer to unique RNA sequences within a context of single species. The website has been subject to continuous improvements focusing on text and sequence similarity searches as well as genome browsing functionality. All RNAcentral data is provided for free and is available for browsing, bulk downloads, and programmatic access at http://rnacentral.org/.

How Much Do rRNA Gene Surveys Underestimate Extant Bacterial Diversity?

The most common practice in studying and cataloguing prokaryotic diversity involves the grouping of sequences into operational taxonomic units (OTUs) at the 97% 16S rRNA gene sequence identity level, often using partial gene sequences, such as PCR-generated amplicons. Due to the high sequence conservation of rRNA genes, organisms belonging to closely related yet distinct species may be grouped under the same OTU. However, it remains unclear how much diversity has been underestimated by this practice. To address this question, we compared the OTUs of genomes defined at the 97% or 98.5% 16S rRNA gene identity level against OTUs of the same genomes defined at the 95% whole-genome average nucleotide identity (ANI), which is a much more accurate proxy for species. Our results show that OTUs resulting from a 98.5% 16S rRNA gene identity cutoff are more accurate than 97% compared to 95% ANI (90.5% versus 89.9% accuracy) but indistinguishable from any other threshold in the 98.29 to 98.78% range. Even with the more stringent thresholds, however, the 16S rRNA gene-based approach commonly underestimates the number of OTUs by ∼12%, on average, compared to the ANI-based approach (∼14% underestimation when using the 97% identity threshold). More importantly, the degree of underestimation can become 50% or more for certain taxa, such as the genera Pseudomonas, Burkholderia, Escherichia, Campylobacter, and Citrobacter These results provide a quantitative view of the degree of underestimation of extant prokaryotic diversity by 16S rRNA gene-defined OTUs and suggest that genomic resolution is often necessary.IMPORTANCE Species diversity is one of the most fundamental pieces of information for community ecology and conservational biology. Therefore, employing accurate proxies for what a species or the unit of diversity is are cornerstones for a large set of microbial ecology and diversity studies. The most common proxies currently used rely on the clustering of 16S rRNA gene sequences at some threshold of nucleotide identity, typically 97% or 98.5%. Here, we explore how well this strategy reflects the more accurate whole-genome-based proxies and determine the frequency with which the high conservation of 16S rRNA sequences masks substantial species-level diversity.

Microbial functional diversity covaries with permafrost thaw-induced environmental heterogeneity in tundra soil.

Permafrost soil in high latitude tundra is one of the largest terrestrial carbon (C) stocks and is highly sensitive to climate warming. Understanding microbial responses to warming-induced environmental changes is critical to evaluating their influences on soil biogeochemical cycles. In this study, a functional gene array (i.e., geochip 4.2) was used to analyze the functional capacities of soil microbial communities collected from a naturally degrading permafrost region in Central Alaska. Varied thaw history was reported to be the main driver of soil and plant differences across a gradient of minimally, moderately, and extensively thawed sites. Compared with the minimally thawed site, the number of detected functional gene probes across the 15-65 cm depth profile at the moderately and extensively thawed sites decreased by 25% and 5%, while the community functional gene ß-diversity increased by 34% and 45%, respectively, revealing decreased functional gene richness but increased community heterogeneity along the thaw progression. Particularly, the moderately thawed site contained microbial communities with the highest abundances of many genes involved in prokaryotic C degradation, ammonification, and nitrification processes, but lower abundances of fungal C decomposition and anaerobic-related genes. Significant correlations were observed between functional gene abundance and vascular plant primary productivity, suggesting that plant growth and species composition could be co-evolving traits together with microbial community composition. Altogether, this study reveals the complex responses of microbial functional potentials to thaw-related soil and plant changes and provides information on potential microbially mediated biogeochemical cycles in tundra ecosystems.

Biotic responses buffer warming-induced soil organic carbon loss in Arctic tundra.

Climate warming can result in both abiotic (e.g., permafrost thaw) and biotic (e.g., microbial functional genes) changes in Arctic tundra. Recent research has incorporated dynamic permafrost thaw in Earth system models (ESMs) and indicates that Arctic tundra could be a significant future carbon (C) source due to the enhanced decomposition of thawed deep soil C. However, warming-induced biotic changes may influence biologically related parameters and the consequent projections in ESMs. How model parameters associated with biotic responses will change under warming and to what extent these changes affect projected C budgets have not been carefully examined. In this study, we synthesized six data sets over 5 years from a soil warming experiment at the Eight Mile Lake, Alaska, into the Terrestrial ECOsystem (TECO) model with a probabilistic inversion approach. The TECO model used multiple soil layers to track dynamics of thawed soil under different treatments. Our results show that warming increased light use efficiency of vegetation photosynthesis but decreased baseline (i.e., environment-corrected) turnover rates of SOC in both the fast and slow pools in comparison with those under control. Moreover, the parameter changes generally amplified over time, suggesting processes of gradual physiological acclimation and functional gene shifts of both plants and microbes. The TECO model predicted that field warming from 2009 to 2013 resulted in cumulative C losses of 224 or 87 g/m2 , respectively, without or with changes in those parameters. Thus, warming-induced parameter changes reduced predicted soil C loss by 61%. Our study suggests that it is critical to incorporate biotic changes in ESMs to improve the model performance in predicting C dynamics in permafrost regions.

ARGs-OAP: online analysis pipeline for antibiotic resistance genes detection from metagenomic data using an integrated structured ARG-database.

MOTIVATION: Environmental dissemination of antibiotic resistance genes (ARGs) has become an increasing concern for public health. Metagenomics approaches can effectively detect broad profiles of ARGs in environmental samples; however, the detection and subsequent classification of ARG-like sequences are time consuming and have been severe obstacles in employing metagenomic methods. We sought to accelerate quantification of ARGs in metagenomic data from environmental samples. RESULTS: A Structured ARG reference database (SARG) was constructed by integrating ARDB and CARD, the two most commonly used databases. SARG was curated to remove redundant sequences and optimized to facilitate query sequence identification by similarity. A database with a hierarchical structure (type-subtype-reference sequence) was then constructed to facilitate classification (assigning ARG-like sequence to type, subtype and reference sequence) of sequences identified through similarity search. Utilizing SARG and a previously proposed hybrid functional gene annotation pipeline, we developed an online pipeline called ARGs-OAP for fast annotation and classification of ARG-like sequences from metagenomic data. We also evaluated and proposed a set of criteria important for efficiently conducting metagenomic analysis of ARGs using ARGs-OAP. AVAILABILITY AND IMPLEMENTATION: Perl script for ARGs-OAP can be downloaded from https://github.com/biofuture/Ublastx_stageone ARGs-OAP can be accessed through http://smile.hku.hk/SARGs CONTACT: zhangt@hku.hk or tiedjej@msu.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

RNAcentral: an international database of ncRNA sequences.

The field of non-coding RNA biology has been hampered by the lack of availability of a comprehensive, up-to-date collection of accessioned RNA sequences. Here we present the first release of RNAcentral, a database that collates and integrates information from an international consortium of established RNA sequence databases. The initial release contains over 8.1 million sequences, including representatives of all major functional classes. A web portal (http://rnacentral.org) provides free access to data, search functionality, cross-references, source code and an integrated genome browser for selected species.

Isothermal assay targeting class 1 integrase gene for environmental surveillance of antibiotic resistance markers.

Antimicrobial resistance genes (ARGs) present in the environment pose a risk to human health due to potential for transfer to human pathogens. Surveillance is an integral part of mitigating environmental dissemination. Quantification of the mobile genetic element class 1 integron-integrase gene (intI1) has been proposed as a surrogate to measuring multiple ARGs. Measurement of such indicator genes can be further simplified by adopting emerging nucleic acids methods such as loop mediated isothermal amplification (LAMP). In this study, LAMP assays were designed and tested for estimating relative abundance of the intI1 gene, which included design of a universal bacteria 16S rRNA gene assay. Following validation of sensitivity and specificity with known bacterial strains, the assays were tested using DNA extracted from river and lake samples. Results showed a significant Pearson correlation (R2 = 0.8) between the intI1 gene LAMP assay and ARG relative abundance (measured via qPCR). To demonstrate the ruggedness of the LAMP assays, experiments were also run in the hands of relatively "untrained" personnel by volunteer undergraduate students at a local community college using a hand-held real-time DNA analysis device - Gene-Z. Overall, results support use of the intI1 gene as an indicator of ARGs and the LAMP assays exhibit the opportunity for volunteers to monitor environmental samples for anthropogenic pollution outside of a specialized laboratory.

Microbial Community Analysis with Ribosomal Gene Fragments from Shotgun Metagenomes.

Shotgun metagenomic sequencing does not depend on gene-targeted primers or PCR amplification; thus, it is not affected by primer bias or chimeras. However, searching rRNA genes from large shotgun Illumina data sets is computationally expensive, and no approach exists for unsupervised community analysis of small-subunit (SSU) rRNA gene fragments retrieved from shotgun data. We present a pipeline, SSUsearch, to achieve the faster identification of short-subunit rRNA gene fragments and enabled unsupervised community analysis with shotgun data. It also includes classification and copy number correction, and the output can be used by traditional amplicon analysis platforms. Shotgun metagenome data using this pipeline yielded higher diversity estimates than amplicon data but retained the grouping of samples in ordination analyses. We applied this pipeline to soil samples with paired shotgun and amplicon data and confirmed bias against Verrucomicrobia in a commonly used V6-V8 primer set, as well as discovering likely bias against Actinobacteria and for Verrucomicrobia in a commonly used V4 primer set. This pipeline can utilize all variable regions in SSU rRNA and also can be applied to large-subunit (LSU) rRNA genes for confirmation of community structure. The pipeline can scale to handle large amounts of soil metagenomic data (5 Gb memory and 5 central processing unit hours to process 38 Gb [1 lane] of trimmed Illumina HiSeq2500 data) and is freely available at https://github.com/dib-lab/SSUsearch under a BSD license.

Ribosomal Database Project: data and tools for high throughput rRNA analysis.

Ribosomal Database Project (RDP; http://rdp.cme.msu.edu/) provides the research community with aligned and annotated rRNA gene sequence data, along with tools to allow researchers to analyze their own rRNA gene sequences in the RDP framework. RDP data and tools are utilized in fields as diverse as human health, microbial ecology, environmental microbiology, nucleic acid chemistry, taxonomy and phylogenetics. In addition to aligned and annotated collections of bacterial and archaeal small subunit rRNA genes, RDP now includes a collection of fungal large subunit rRNA genes. RDP tools, including Classifier and Aligner, have been updated to work with this new fungal collection. The use of high-throughput sequencing to characterize environmental microbial populations has exploded in the past several years, and as sequence technologies have improved, the sizes of environmental datasets have increased. With release 11, RDP is providing an expanded set of tools to facilitate analysis of high-throughput data, including both single-stranded and paired-end reads. In addition, most tools are now available as open source packages for download and local use by researchers with high-volume needs or who would like to develop custom analysis pipelines.

Fungal identification using a Bayesian classifier and the Warcup training set of internal transcribed spacer sequences.

Fungi are key organisms in many ecological processes and communities. Rapid and low cost surveys of the fungal members of a community can be undertaken by isolating and sequencing a taxonomically informative genomic region, such as the ITS (internal transcribed spacer), from DNA extracted from a metagenomic sample, and then classifying these sequences to determine which organisms are present. This paper announces the availability of the Warcup ITS training set and shows how it can be used with the Ribosomal Database Project (RDP) Bayesian Classifier to rapidly and accurately identify fungi using ITS sequences. The classifications can be down to species level and use conventional literature-based mycological nomenclature and taxonomic assignments.

A scalable and accurate targeted gene assembly tool (SAT-Assembler) for next-generation sequencing data.

Gene assembly, which recovers gene segments from short reads, is an important step in functional analysis of next-generation sequencing data. Lacking quality reference genomes, de novo assembly is commonly used for RNA-Seq data of non-model organisms and metagenomic data. However, heterogeneous sequence coverage caused by heterogeneous expression or species abundance, similarity between isoforms or homologous genes, and large data size all pose challenges to de novo assembly. As a result, existing assembly tools tend to output fragmented contigs or chimeric contigs, or have high memory footprint. In this work, we introduce a targeted gene assembly program SAT-Assembler, which aims to recover gene families of particular interest to biologists. It addresses the above challenges by conducting family-specific homology search, homology-guided overlap graph construction, and careful graph traversal. It can be applied to both RNA-Seq and metagenomic data. Our experimental results on an Arabidopsis RNA-Seq data set and two metagenomic data sets show that SAT-Assembler has smaller memory usage, comparable or better gene coverage, and lower chimera rate for assembling a set of genes from one or multiple pathways compared with other assembly tools. Moreover, the family-specific design and rapid homology search allow SAT-Assembler to be naturally compatible with parallel computing platforms. The source code of SAT-Assembler is available at https://sourceforge.net/projects/sat-assembler/. The data sets and experimental settings can be found in supplementary material.

In-feed antibiotic effects on the swine intestinal microbiome.

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.

A Sensitive and Accurate protein domain cLassification Tool (SALT) for short reads.

MOTIVATION: Protein domain classification is an important step in functional annotation for next-generation sequencing data. For RNA-Seq data of non-model organisms that lack quality or complete reference genomes, existing protein domain analysis pipelines are applied to short reads directly or to contigs that are generated using de novo sequence assembly tools. However, these strategies do not provide satisfactory performance in classifying short reads into their native domain families. RESULTS: We introduce SALT, a protein domain classification tool based on profile hidden Markov models and graph algorithms. SALT carefully incorporates the characteristics of reads that are sequenced from the domain regions and assembles them into contigs based on a supervised graph construction algorithm. We applied SALT to two RNA-Seq datasets of different read lengths and quantified its performance using the available protein domain annotations and the reference genomes. Compared with existing strategies, SALT showed better sensitivity and accuracy. In the third experiment, we applied SALT to a non-model organism. The experimental results demonstrated that it identified more transcribed protein domain families than other tested classifiers. AVAILABILITY: The source code and supplementary data are available at https://sourceforge.net/projects/salt1/

The Genomic Standards Consortium.

A vast and rich body of information has grown up as a result of the world's enthusiasm for 'omics technologies. Finding ways to describe and make available this information that maximise its usefulness has become a major effort across the 'omics world. At the heart of this effort is the Genomic Standards Consortium (GSC), an open-membership organization that drives community-based standardization activities, Here we provide a short history of the GSC, provide an overview of its range of current activities, and make a call for the scientific community to join forces to improve the quality and quantity of contextual information about our public collections of genomes, metagenomes, and marker gene sequences.

History and impact of RDP: a legacy from Carl Woese to microbiology.

The Ribosomal Database Project (RDP) grew out of Carl Woese's vision of how rRNA comparative methods could transform biology. First at the University of Illinois Urbana-Champaign, and later at Michigan State University's Center for Microbial Ecology, the project has grown from a few hundred to several million rRNA gene sequences. In the years since Woese started the RDP, publications describing the database and related tools have been cited over 11,000 times in journals spanning a wide range of disciplines, while the RDP website is accessed by 10,000 researchers in over 20,000 analysis sessions each month. This article describes the history of RDP's development over the last two decades.