Nonpareil 3: Fast Estimation of Metagenomic Coverage and Sequence Diversity

Insights into Ammonia Adaptation and Methanogenic Precursor Oxidation by Genome-Centric Analysis

Ammonia released from the degradation of protein and/or urea usually leads to suboptimal anaerobic digestion (AD) when N-rich organic waste is used. However, the insights behind the differential ammonia tolerance of anaerobic microbiomes remain an enigma. In this study, the cultivation in synthetic medium with different carbon sources (acetate, methanol, formate, and H2/CO2) shaped a common initial inoculum into four unique ammonia-tolerant syntrophic populations. Specifically, various levels of ammonia tolerance were observed: consortia fed with methanol and H2/CO2 could grow at ammonia levels up to 7.25 g NH+-N/L, whereas the other two groups (formate and acetate) only thrived at 5.25 and 4.25 g NH+-N/L, respectively. Metabolic reconstruction highlighted that this divergent microbiome might be achieved by complementary metabolisms to maximize biomethane recovery from carbon sources, thus indicating the importance of the syntrophic community in the AD of N-rich substrates. Besides, sodium/proton antiporter operon, osmoprotectant/K+ regulator, and osmoprotectant synthesis operon may function as the main drivers of adaptation to the ammonia stress. Moreover, energy from the substrate-level phosphorylation and multiple energy-converting hydrogenases (e.g., Ech and Eha) could aid methanogens to balance the energy request for anabolic activities and contribute to thriving when exposed to high ammonia levels.

Microbiome-Informed Food Safety and Quality: Longitudinal Consistency and Cross-Sectional Distinctiveness of Retail Chicken Breast Microbiomes

Chicken has recently overtaken beef as the most-consumed meat in the United States. The growing popularity of chicken is accompanied by frequent occurrences of foodborne pathogens and increasing concerns over antibiotic usage. Our study represents a proof-of-concept investigation into the possibility and practicality of leveraging microbiome-informed food safety and quality. Through a longitudinal and cross-sectional survey, we established the chicken microbiome as a robust and multifaceted food microbiology attribute that could provide a variety of safety and quality information and retain systematic signals characteristic of overall processing environments.

Microorganisms and their communities on foods are important determinants and indicators of food safety and quality. Despite growing interests in studying food and food-related microbiomes, how effective and practical it is to glean various food safety and quality information from food commodity microbiomes remains underinvestigated. Microbiomes of retail chicken breast from 4 processing establishments in 3 major U.S. broiler production states displayed longitudinal consistency over 7 months and cross-sectional distinctiveness associated with individual processing environments. Packaging type and processing environment but not antibiotic usage and seasonality affected composition and diversity of the microbiomes. Low abundances of antimicrobial resistance genes were found on chicken breasts, and no significant resistome difference was observed between antibiotic-free and conventional products. Benchmarked by culture enrichment, shotgun metagenomics sequencing delivered sensitive and specific detection of Salmonella enterica from chicken breasts.

IMPORTANCE Chicken has recently overtaken beef as the most-consumed meat in the United States. The growing popularity of chicken is accompanied by frequent occurrences of foodborne pathogens and increasing concerns over antibiotic usage. Our study represents a proof-of-concept investigation into the possibility and practicality of leveraging microbiome-informed food safety and quality. Through a longitudinal and cross-sectional survey, we established the chicken microbiome as a robust and multifaceted food microbiology attribute that could provide a variety of safety and quality information and retain systematic signals characteristic of overall processing environments.

Nitrite Oxidizer Activity and Community Are More Responsive Than Their Abundance to Ammonium-Based Fertilizer in an Agricultural Soil

Autotrophic nitrification is mediated by ammonia oxidizing bacteria (AOB) or ammonia oxidizing archaea (AOA) and nitrite oxidizing bacteria (NOB). Mounting studies have examined the impact of nitrogen (N) fertilization on the dynamic and diversity of AOA and AOB, while we have limited information on the response of the activity, abundance, and diversity of NOB to N fertilization. We investigated the influence of organic and inorganic N fertilizers on soil NOB in silage corn field plots that received contrasting nitrogen (N) treatments: control (no additional N), ammonium sulfate (AS 100 and 200 kg N ha-1), and compost (200 kg N ha-1). Nitrifying community was examined using a universal marker (16S rRNA gene), functional gene markers (AOB amoA and Nitrospira nxrB), and metagenomics. The overall nitrifying community was not altered after the first fertilization but was significantly shifted by 4-year repeated application of ammonium fertilizers. Nitrospira were the dominant NOB (>99.7%) in our agricultural soil. Both community compositions of AOB and Nitrospira were significantly changed by ammonium fertilizers but not by compost after 4 years of repeated applications. All nitrifiers, including comammox, were recovered in soil metagenomes based on a gene-targeted assembly, but their sequence counts were very low. Although N treatment did not affect the abundance of Nitrospira nxrB determined by real-time quantitative PCR, ammonium fertilizers significantly promoted rates of potential nitrite oxidation determined at 0.15 mM nitrite in soil slurries. Understanding the response of both ammonia oxidizers and nitrite oxidizers to N fertilization may initiate or improve strategies for mitigating potential environmental impacts of nitrate production in agricultural ecosystems.

Iterative subtractive binning of freshwater chronoseries metagenomes identifies over 400 novel species and their ecologic preferences

Recent advances in sequencing technology and bioinformatic pipelines have allowed unprecedented access to the genomes of yet-uncultivated microorganisms from diverse environments. However, the catalogue of freshwater genomes remains limited, and most genome recovery attempts in freshwater ecosystems have only targeted specific taxa. Here, we present a genome recovery pipeline incorporating iterative subtractive binning, and apply it to a time series of 100 metagenomic datasets from seven connected lakes and estuaries along the Chattahoochee River (Southeastern USA). Our set of metagenome-assembled genomes (MAGs) represents >400 yet-unnamed genomospecies, substantially increasing the number of high-quality MAGs from freshwater lakes. We propose names for two novel species: 'Candidatus Elulimicrobium humile' ('Ca. Elulimicrobiota', 'Patescibacteria') and 'Candidatus Aquidulcis frankliniae' ('Chloroflexi'). Collectively, our MAGs represented about half of the total microbial community at any sampling point. To evaluate the prevalence of these genomospecies in the chronoseries, we introduce methodologies to estimate relative abundance and habitat preference that control for uneven genome quality and sample representation. We demonstrate high degrees of habitat-specialization and endemicity for most genomospecies in the Chattahoochee lakes. Wider ecological ranges characterized smaller genomes with higher coding densities, indicating an overall advantage of smaller, more compact genomes for cosmopolitan distributions.

Terabase-scale metagenome coassembly with MetaHipMer

Metagenome sequence datasets can contain terabytes of reads, too many to be coassembled together on a single shared-memory computer; consequently, they have only been assembled sample by sample (multiassembly) and combining the results is challenging. We can now perform coassembly of the largest datasets using MetaHipMer, a metagenome assembler designed to run on supercomputers and large clusters of compute nodes. We have reported on the implementation of MetaHipMer previously; in this paper we focus on analyzing the impact of very large coassembly. In particular, we show that coassembly recovers a larger genome fraction than multiassembly and enables the discovery of more complete genomes, with lower error rates, whereas multiassembly recovers more dominant strain variation. Being able to coassemble a large dataset does not preclude one from multiassembly; rather, having a fast, scalable metagenome assembler enables a user to more easily perform coassembly and multiassembly, and assemble both abundant, high strain variation genomes, and low-abundance, rare genomes. We present several assemblies of terabyte datasets that could never be coassembled before, demonstrating MetaHipMer’s scaling power. MetaHipMer is available for public use under an open source license and all datasets used in the paper are available for public download.

Multimodal Approach to Assessment of Fecal Microbiota Donors based on Three Complementary Methods

Methods of stool assessment are mostly focused on next-generation sequencing (NGS) or classical culturing, but only rarely both. We conducted a series of experiments using a multi-method approach to trace the stability of gut microbiota in various donors over time, to find the best method for the proper selection of fecal donors and to find "super-donor" indicators. Ten consecutive stools donated by each of three donors were used for the experiments (30 stools in total). The experiments assessed bacterial viability measured by flow cytometry, stool culturing on different media and in various conditions, and NGS (90 samples in total). There were no statistically significant differences between live and dead cell numbers; however, we found a group of cells classified as not-dead-not-alive, which may be possibly important in selection of "good" donors. Donor C, being a regular stool donor, was characterized by the largest number of cultivable species (64). Cultivable core microbiota (shared by all donors) was composed of only 16 species. ANCOM analysis of NGS data highlighted particular genera to be more abundant in one donor vs. the others. There was a correlation between the not-dead-not-alive group found in flow cytometry and Anaeroplasma found by NGS, and we could distinguish a regular stool donor from the others. In this work, we showed that combining various methods of microbiota assessment gives more information than each method separately.

Phages Actively Challenge Niche Communities in Antarctic Soils

By modulating the structure, diversity, and trophic outputs of microbial communities, phages play crucial roles in many biomes. In oligotrophic polar deserts, the effects of katabatic winds, constrained nutrients, and low water availability are known to limit microbial activity. Although phages may substantially govern trophic interactions in cold deserts, relatively little is known regarding the precise ecological mechanisms. Here, we provide the first evidence of widespread antiphage innate immunity in Antarctic environments using metagenomic sequence data from hypolith communities as model systems. In particular, immunity systems such as DISARM and BREX are shown to be dominant systems in these communities. Additionally, we show a direct correlation between the CRISPR-Cas adaptive immunity and the metavirome of hypolith communities, suggesting the existence of dynamic host-phage interactions. In addition to providing the first exploration of immune systems in cold deserts, our results suggest that phages actively challenge niche communities in Antarctic polar deserts. We provide evidence suggesting that the regulatory role played by phages in this system is an important determinant of bacterial host interactions in this environment.IMPORTANCE In Antarctic environments, the combination of both abiotic and biotic stressors results in simple trophic levels dominated by microbiomes. Although the past two decades have revealed substantial insights regarding the diversity and structure of microbiomes, we lack mechanistic insights regarding community interactions and how phages may affect these. By providing the first evidence of widespread antiphage innate immunity, we shed light on phage-host dynamics in Antarctic niche communities. Our analyses reveal several antiphage defense systems, including DISARM and BREX, which appear to dominate in cold desert niche communities. In contrast, our analyses revealed that genes which encode antiphage adaptive immunity were underrepresented in these communities, suggesting lower infection frequencies in cold edaphic environments. We propose that by actively challenging niche communities, phages play crucial roles in the diversification of Antarctic communities.

Metagenomic- and Cultivation-Based Exploration of Anaerobic Chloroform Biotransformation in Hypersaline Sediments as Natural Source of Chloromethanes

Chloroform (CF) is an environmental contaminant that can be naturally formed in various environments ranging from forest soils to salt lakes. Here we investigated CF removal potential in sediments obtained from hypersaline lakes in Western Australia. Reductive dechlorination of CF to dichloromethane (DCM) was observed in enrichment cultures derived from sediments of Lake Strawbridge, which has been reported as a natural source of CF. No CF removal was observed in abiotic control cultures without artificial electron donors, indicating biotic CF dechlorination in the enrichment cultures. Increasing vitamin B₁₂ concentration from 0.04 to 4 µM in enrichment cultures enhanced CF removal and reduced DCM formation. In cultures amended with 4 µM vitamin B₁₂ and ¹³C labelled CF, formation of ¹³CO₂ was detected. Known organohalide-respiring bacteria and reductive dehalogenase genes were neither detected using quantitative PCR nor metagenomic analysis of the enrichment cultures. Rather, members of the order Clostridiales, known to co-metabolically transform CF to DCM and CO₂, were detected. Accordingly, metagenome-assembled genomes of Clostridiales encoded enzymatic repertoires for the Wood-Ljungdahl pathway and cobalamin biosynthesis, which are known to be involved in fortuitous and nonspecific CF transformation. This study indicates that hypersaline lake microbiomes may act as a filter to reduce CF emission to the atmosphere.

The preceding root system drives the composition and function of the rhizosphere microbiome

BACKGROUND

The soil environment is responsible for sustaining most terrestrial plant life, yet we know surprisingly little about the important functions carried out by diverse microbial communities in soil. Soil microbes that inhabit the channels of decaying root systems, the detritusphere, are likely to be essential for plant growth and health, as these channels are the preferred locations of new root growth. Understanding the microbial metagenome of the detritusphere, and how it responds to agricultural management such as crop rotations and soil tillage, is vital for improving global food production.

RESULTS

This study establishes an in-depth soil microbial gene catalogue based on the living-decaying rhizosphere niches in a cropping soil. The detritusphere microbiome regulates the composition and function of the rhizosphere microbiome to a greater extent than plant type: rhizosphere microbiomes of wheat and chickpea were homogenous (65-87% similarity) in the presence of decaying root (DR) systems but were heterogeneous (3-24% similarity) where DR was disrupted by tillage. When the microbiomes of the rhizosphere and the detritusphere interact in the presence of DR, there is significant degradation of plant root exudates by the rhizosphere microbiome, and genes associated with membrane transporters, carbohydrate and amino acid metabolism are enriched.

CONCLUSIONS

The study describes the diversity and functional capacity of a high-quality soil microbial metagenome. The results demonstrate the contribution of the detritusphere microbiome in determining the metagenome of developing root systems. Modifications in root microbial function through soil management can ultimately govern plant health, productivity and food security.

Genome-resolved metagenomics to study co-occurrence patterns and intraspecific heterogeneity among plant pathogen metapopulations

Assessment of pathogen diversity in agricultural fields is essential for informing management decisions and the development of resistant plant varieties. However, many population genomic studies have relied on culture-based approaches that do not provide quantitative assessment of pathogen populations at the field-level or the associated host microbiome. Here, we applied whole-genome shotgun sequencing of microbial DNA extracted directly from the washings of pooled leaf samples, collected from individual tomato and pepper fields in Alabama that displayed the classical symptoms of bacterial spot disease caused by Xanthomonas spp. Our results revealed that while the occurrence of both X. perforans and X. euvesicatoria within fields was limited, evidence of co-occurrence of up to three distinct X. perforans genotypes was obtained in 7 of 10 tomato fields sampled. These population dynamics were accompanied by the corresponding type 3 secreted effector repertoires associated with the co-occurring X. perforans genotypes, indicating that metapopulation structure within fields should be considered when assessing the adaptive potential of X. perforans. Finally, analysis of microbial community composition revealed that co-occurrence of the bacterial spot pathogens Pseudomonas cichorii and Xanthomonas spp. is common in Alabama fields and provided evidence for the non-random association of several other human and plant opportunists.

Disinfection exhibits systematic impacts on the drinking water microbiome

Limiting microbial growth during drinking water distribution is achieved either by maintaining a disinfectant residual or through nutrient limitation without using a disinfectant. The impact of these contrasting approaches on the drinking water microbiome is not systematically understood. We use genome-resolved metagenomics to compare the structure, metabolic traits, and population genomes of drinking water microbiome samples from bulk drinking water across multiple full-scale disinfected and non-disinfected drinking water systems. Microbial communities cluster at the structural- and functional potential-level based on the presence/absence of a disinfectant residual. Disinfectant residual alone explained 17 and 6.5% of the variance in structure and functional potential of the drinking water microbiome, respectively, despite including multiple drinking water systems with variable source waters and source water communities and treatment strategies. The drinking water microbiome is structurally and functionally less diverse and variable across disinfected compared to non-disinfected systems. While bacteria were the most abundant domain, archaea and eukaryota were more abundant in non-disinfected and disinfected systems, respectively. Community-level differences in functional potential were driven by enrichment of genes associated with carbon and nitrogen fixation in non-disinfected systems and γ-aminobutyrate metabolism in disinfected systems likely associated with the recycling of amino acids. Genome-level analyses for a subset of phylogenetically-related microorganisms suggests that disinfection selects for microorganisms capable of using fatty acids, presumably from microbial decay products, via the glyoxylate cycle. Overall, we find that disinfection exhibits systematic selective pressures on the drinking water microbiome and may select for microorganisms able to utilize microbial decay products originating from disinfection-inactivated microorganisms. Video abstract.

Characteristics of Wetting-Induced Bacteriophage Blooms in Biological Soil Crust

Biological soil crusts (biocrusts) are photosynthetic "hot spots" in deserts and cover ∼12% of the Earth's terrestrial surface, and yet they face an uncertain future given expected shifts in rainfall events. Laboratory wetting of biocrust communities is known to cause a bloom of Firmicutes which rapidly become dominant community members within 2 days after emerging from a sporulated state. We hypothesized that their bacteriophages (phages) would respond to such a dramatic increase in their host's abundance. In our experiment, wetting caused Firmicutes to bloom and triggered a significant depletion of cyanobacterial diversity. We used genome-resolved metagenomics to link phage to their hosts and found that the bloom of the genus Bacillus correlated with a dramatic increase in the number of Caudovirales phages targeting these diverse spore-formers (r = 0.762). After 2 days, we observed dramatic reductions in the relative abundances of Bacillus, while the number of Bacillus phages continued to increase, suggestive of a predator-prey relationship. We found predicted auxiliary metabolic genes (AMGs) associated with sporulation in several Caudovirales genomes, suggesting that phages may influence and even benefit from sporulation dynamics in biocrusts. Prophage elements and CRISPR-Cas repeats in Firmicutes metagenome-assembled genomes (MAGs) provide evidence of recent infection events by phages, which were corroborated by mapping viral contigs to their host MAGs. Combined, these findings suggest that the blooming Firmicutes become primary targets for biocrust Caudovirales phages, consistent with the classical "kill-the-winner" hypothesis.IMPORTANCE This work forms part of an overarching research theme studying the effects of a changing climate on biological soil crust (biocrust) in the Southwestern United States. To our knowledge, this study was the first to characterize bacteriophages in biocrust and offers a view into the ecology of phages in response to a laboratory wetting experiment. The phages identified here represent lineages of Caudovirales, and we found that the dynamics of their interactions with their Firmicutes hosts explain the collapse of a bacterial bloom that was induced by wetting. Moreover, we show that phages carried host-altering metabolic genes and found evidence of proviral infection and CRISPR-Cas repeats within host genomes. Our results suggest that phages exert controls on population density by lysing dominant bacterial hosts and that they further impact biocrust by acquiring host genes for sporulation. Future research should explore how dominant these phages are in other biocrust communities and quantify how much the control and lysis of blooming populations contributes to nutrient cycling in biocrusts.

Metagenomic Signatures of Gut Infections Caused by Different Escherichia coli Pathotypes

Escherichia coli is a leading contributor to infectious diarrhea and child mortality worldwide, but it remains unknown how alterations in the gut microbiome vary for distinct E. coli pathotype infections and whether these signatures can be used for diagnostic purposes. Further, the majority of enteric diarrheal infections are not diagnosed with respect to their etiological agent(s) due to technical challenges. To address these issues, we devised a novel approach that combined traditional, isolate-based and molecular-biology techniques with metagenomics analysis of stool samples and epidemiological data. Application of this pipeline to children enrolled in a case-control study of diarrhea in Ecuador showed that, in about half of the cases where an E. coli pathotype was detected by culture and PCR, E. coli was likely not the causative agent based on the metagenome-derived low relative abundance, the level of clonality, and/or the virulence gene content. Our results also showed that diffuse adherent E. coli (DAEC), a pathotype that is generally underrepresented in previous studies of diarrhea and thus, thought not to be highly virulent, caused several small-scale diarrheal outbreaks across a rural to urban gradient in Ecuador. DAEC infections were uniquely accompanied by coelution of large amounts of human DNA and conferred significant shifts in the gut microbiome composition relative to controls or infections caused by other E. coli pathotypes. Our study shows that diarrheal infections can be efficiently diagnosed for their etiological agent and categorized based on their effects on the gut microbiome using metagenomic tools, which opens new possibilities for diagnostics and treatment.IMPORTANCEE. coli infectious diarrhea is an important contributor to child mortality worldwide. However, diagnosing and thus treating E. coli infections remain challenging due to technical and other reasons associated with the limitations of the traditional culture-based techniques and the requirement to apply Koch's postulates. In this study, we integrated traditional microbiology techniques with metagenomics and epidemiological data in order to identify cases of diarrhea where E. coli was most likely the causative disease agent and evaluate specific signatures in the disease-state gut microbiome that distinguish between diffuse adherent, enterotoxigenic, and enteropathogenic E. coli pathotypes. Therefore, our methodology and results should be highly relevant for diagnosing and treating diarrheal infections and have important applications in public health.

A Bioinformatics Guide to Plant Microbiome Analysis

Recent evidence for intimate relationship of plants with their microbiota shows that plants host individual and diverse microbial communities that are essential for their survival. Understanding their relatedness using genome-based and high-throughput techniques remains a hot topic in microbiome research. Molecular analysis of the plant holobiont necessitates the application of specific sampling and preparatory steps that also consider sources of unwanted information, such as soil, co-amplified plant organelles, human DNA, and other contaminations. Here, we review state-of-the-art and present practical guidelines regarding experimental and computational aspects to be considered in molecular plant-microbiome studies. We discuss sequencing and "omics" techniques with a focus on the requirements needed to adapt these methods to individual research approaches. The choice of primers and sequence databases is of utmost importance for amplicon sequencing, while the assembly and binning of shotgun metagenomic sequences is crucial to obtain quality data. We discuss specific bioinformatic workflows to overcome the limitation of genome database resources and for covering large eukaryotic genomes such as fungi. In transcriptomics, it is necessary to account for the separation of host mRNA or dual-RNAseq data. Metaproteomics approaches provide a snapshot of the protein abundances within a plant tissue which requires the knowledge of complete and well-annotated plant genomes, as well as microbial genomes. Metabolomics offers a powerful tool to detect and quantify small molecules and molecular changes at the plant-bacteria interface if the necessary requirements with regard to (secondary) metabolite databases are considered. We highlight data integration and complementarity which should help to widen our understanding of the interactions among individual players of the plant holobiont in the future.

Discovery of hyperstable carbohydrate-active enzymes through metagenomics of extreme environments

The enzymes from hyperthermophilic microorganisms populating volcanic sites represent interesting cases of protein adaptation and biotransformations under conditions where conventional enzymes quickly denature. The difficulties in cultivating extremophiles severely limit access to this class of biocatalysts. To circumvent this problem, we embarked on the exploration of the biodiversity of the solfatara Pisciarelli, Agnano (Naples, Italy), to discover hyperthermophilic carbohydrate-active enzymes (CAZymes) and to characterize the entire set of such enzymes in this environment (CAZome). Here, we report the results of the metagenomic analysis of two mud/water pools that greatly differ in both temperature and pH (T = 85 °C and pH 5.5; T = 92 °C and pH 1.5, for Pool1 and Pool2, respectively). DNA deep sequencing and following in silico analysis led to 14 934 and 17 652 complete ORFs in Pool1 and Pool2, respectively. They exclusively belonged to archaeal cells and viruses with great genera variance within the phylum Crenarchaeota, which reflected the difference in temperature and pH of the two Pools. Surprisingly, 30% and 62% of all of the reads obtained from Pool1 and 2, respectively, had no match in nucleotide databanks. Genes associated with carbohydrate metabolism were 15% and 16% of the total in the two Pools, with 278 and 308 putative CAZymes in Pool1 and 2, corresponding to ~ 2.0% of all ORFs. Biochemical characterization of two CAZymes of a previously unknown archaeon revealed a novel subfamily GH5_19 β-mannanase/β-1,3-glucanase whose hemicellulose specificity correlates with the vegetation surrounding the sampling site, and a novel NAD⁺ -dependent GH109 with a previously unreported β-N-acetylglucosaminide/β-glucoside specificity. DATABASES: The sequencing reads are available in the NCBI Sequence Read Archive (SRA) database under the accession numbers SRR7545549 (Pool1) and SRR7545550 (Pool2). The sequences of GH5_Pool2 and GH109_Pool2 are available in GenBank database under the accession numbers MK869723 and MK86972, respectively. The environmental data relative to Pool1 and Pool2 (NCBI BioProject PRJNA481947) are available in the Biosamples database under the accession numbers SAMN09692669 (Pool1) and SAMN09692670 (Pool2).

Review, Evaluation, and Directions for Gene-Targeted Assembly for Ecological Analyses of Metagenomes

Shotgun metagenomics has greatly advanced our understanding of microbial communities over the last decade. Metagenomic analyses often include assembly and genome binning, computationally daunting tasks especially for big data from complex environments such as soil and sediments. In many studies, however, only a subset of genes and pathways involved in specific functions are of interest; thus, it is not necessary to attempt global assembly. In addition, methods that target genes can be computationally more efficient and produce more accurate assembly by leveraging rich databases, especially for those genes that are of broad interest such as those involved in biogeochemical cycles, biodegradation, and antibiotic resistance or used as phylogenetic markers. Here, we review six gene-targeted assemblers with unique algorithms for extracting and/or assembling targeted genes: Xander, MegaGTA, SAT-Assembler, HMM-GRASPx, GenSeed-HMM, and MEGAN. We tested these tools using two datasets with known genomes, a synthetic community of artificial reads derived from the genomes of 17 bacteria, shotgun sequence data from a mock community with 48 bacteria and 16 archaea genomes, and a large soil shotgun metagenomic dataset. We compared assemblies of a universal single copy gene (rplB) and two N cycle genes (nifH and nirK). We measured their computational efficiency, sensitivity, specificity, and chimera rate and found Xander and MegaGTA, which both use a probabilistic graph structure to model the genes, have the best overall performance with all three datasets, although MEGAN, a reference matching assembler, had better sensitivity with synthetic and mock community members chosen from its reference collection. Also, Xander and MegaGTA are the only tools that include post-assembly scripts tuned for common molecular ecology and diversity analyses. Additionally, we provide a mathematical model for estimating the probability of assembling targeted genes in a metagenome for estimating required sequencing depth.

Prevalence of viral photosynthesis genes along a freshwater to saltwater transect in Southeast USA

Bacteriophages encode host-acquired functional genes known as auxiliary metabolic genes (AMGs). Photosynthesis AMGs are commonly found in marine cyanobacteria-infecting Myoviridae and Podoviridae cyanophages, but their ecology remains understudied in freshwater environments. To advance knowledge of this issue, we analysed viral metagenomes collected in the summertime for four years from five lakes and two estuarine locations interconnected by the Chattahoochee River, Southeast USA. Sequences representing ten different AMGs were recovered and found to be prevalent in all sites. Most freshwater AMGs were 10-fold less abundant than estuarine and marine AMGs and were encoded by novel Myoviridae and Podoviridae cyanophage genera. Notably, several of the corresponding viral genomes showed endemism to a specific province along the river. This translated into psbA gene phylogenetic clustering patterns that matched a marine vs. freshwater origin indicating that psbA may serve as a robust classification and source-tracking biomarker. Genomes classified in a novel viral lineage represented by isolate S-EIVl contained psbA, which is unprecedented for this lineage. Collectively, our findings indicated that the acquisition of photosynthesis AMGs is a widespread strategy used by cyanophages in aquatic ecosystems, and further indicated the existence of viral provinces in which certain viral species and/or genotypes are locally abundant.

Literature-based, manually-curated database of PCR primers for the detection of antibiotic resistance genes in various environments

A literature-based, manually-curated database of PCR primers for the detection of antibiotic resistance genes in various environments was constructed (LCPDb-ARG; lcpdb.ddg.biol.uw.edu.pl and lcpdb.ddlemb.com). Currently, this database is comprised of 607 PCR primer pairs designed for the amplification of various genes conferring resistance to antibiotics representing 10 classes of antimicrobial agents. Three parameters were assigned for each primer pair: specificity, efficacy and taxonomic efficacy. These parameters were evaluated using a novel bioinformatic tool, UniPriVal, developed for this study. UniPriVal was used to validate each primer pair against various databases, including the Bacterial Antimicrobial Resistance Reference Gene Database (BARRGDB) and those of the National Center of Biotechnology Information (NCBI). Primer pairs specific for each antibiotic resistance gene were ranked based on their model success metric value. To validate the utility and correctness of the information collected in the LCPDb-ARG, selected primer pairs were tested in bioinformatic and experimental PCR surveys. To our knowledge, this is the first database designed to facilitate PCR monitoring of the occurrence and diversity of antibiotic resistance genes in environmental and clinical samples. The internal validation system of this user-friendly application enables the quantified ranking of PCR primer pairs, which assists selection of the best primers for each application.

Responses of tundra soil microbial communities to half a decade of experimental warming at two critical depths

Northern-latitude tundra soils harbor substantial carbon (C) stocks that are highly susceptible to microbial degradation with rising global temperatures. Understanding the magnitude and direction (e.g., C release or sequestration) of the microbial responses to warming is necessary to accurately model climate change. In this study, Alaskan tundra soils were subjected to experimental in situ warming by ∼1.1 °C above ambient temperature, and the microbial communities were evaluated using metagenomics after 4.5 years, at 2 depths: 15 to 25 cm (active layer at outset of the experiment) and 45 to 55 cm (transition zone at the permafrost/active layer boundary at the outset of the experiment). In contrast to small or insignificant shifts after 1.5 years of warming, 4.5 years of warming resulted in significant changes to the abundances of functional traits and the corresponding taxa relative to control plots (no warming), and microbial shifts differed qualitatively between the two soil depths. At 15 to 25 cm, increased abundances of carbohydrate utilization genes were observed that correlated with (increased) measured ecosystem carbon respiration. At the 45- to 55-cm layer, increased methanogenesis potential was observed, which corresponded with a 3-fold increase in abundance of a single archaeal clade of the Methanosarcinales order, increased annual thaw duration (45.3 vs. 79.3 days), and increased CH₄ emissions. Collectively, these data demonstrate that the microbial responses to warming in tundra soil are rapid and markedly different between the 2 critical soil layers evaluated, and identify potential biomarkers for the corresponding microbial processes that could be important in modeling.

Evaluation of Sequencing Library Preparation Protocols for Viral Metagenomic Analysis from Pristine Aquifer Groundwaters

Viral ecology of terrestrial habitats is yet-to be extensively explored, in particular the terrestrial subsurface. One problem in obtaining viral sequences from groundwater aquifer samples is the relatively low amount of virus particles. As a result, the amount of extracted DNA may not be sufficient for direct sequencing of such samples. Here we compared three DNA amplification methods to enrich viral DNA from three pristine limestone aquifer assemblages of the Hainich Critical Zone Exploratory to evaluate potential bias created by the different amplification methods as determined by viral metagenomics. Linker amplification shotgun libraries resulted in lowest redundancy among the sequencing reads and showed the highest diversity, while multiple displacement amplification produced the highest number of contigs with the longest average contig size, suggesting a combination of these two methods is suitable for the successful enrichment of viral DNA from pristine groundwater samples. In total, we identified 27,173, 5,886 and 32,613 viral contigs from the three samples from which 11.92 to 18.65% could be assigned to taxonomy using blast. Among these, members of the Caudovirales order were the most abundant group (52.20 to 69.12%) dominated by Myoviridae and Siphoviridae. Those, and the high number of unknown viral sequences, substantially expand the known virosphere.

Comparative Metagenomics Reveals Enhanced Nutrient Cycling Potential after 2 Years of Biochar Amendment in a Tropical Oxisol

The complex structural and functional responses of agricultural soil microbial communities to the addition of carbonaceous compounds such as biochar remain poorly understood. This severely limits the predictive ability for both the potential enhancement of soil fertility and greenhouse gas mitigation. In this study, we utilized shotgun metagenomics in order to decipher changes in the microbial community in soil microcosms after 14 days of incubation at 23°C, which contained soils from biochar-amended and control plots cultivated with Napier grass. Our analyses revealed that biochar-amended soil microbiomes exhibited significant shifts in both community composition and predicted metabolism. Key metabolic pathways related to carbon turnover, such as the utilization of plant-derived carbohydrates as well as denitrification, were enriched under biochar amendment. These community shifts were in part associated with increased soil carbon, such as labile and aromatic carbon compounds, which was likely stimulated by the increased available nutrients associated with biochar amendment. These findings indicate that the soil microbiome response to the combination of biochar addition and to incubation conditions confers enhanced nutrient cycling and a small decrease in CO₂ emissions and potentially mitigates nitrous oxide emissions.IMPORTANCE The incorporation of biochar into soil is a promising management strategy for sustainable agriculture owing to its potential to sequester carbon and improve soil fertility. Expanding the addition of biochar to large-scale agriculture hinges on its lasting beneficial effects on the microbial community. However, there exists a significant knowledge gap regarding the specific role that biochar plays in altering the key biological soil processes that influence plant growth and carbon storage in soil. Previous studies that examined the soil microbiome under biochar amendment principally characterized only how the composition alters in response to biochar amendment. In the present study, we shed light on the functional alterations of the microbial community response 2 years after biochar amendment. Our results show that biochar increased the abundance of genes involved in denitrification and carbon turnover and that biochar-amended soil microcosms had a reduction in cumulative CO₂ production.

Sources and selection of snow-specific microbial communities in a Greenlandic sea ice snow cover

Sea ice and its snow cover are critical for global processes including climate regulation and biogeochemical cycles. Despite an increase in studies focused on snow microorganisms, the ecology of snow inhabitants remains unclear. In this study, we investigated sources and selection of a snowpack-specific microbial community by comparing metagenomes from samples collected in a Greenlandic fjord within a vertical profile including atmosphere, snowpack with four distinct layers of snow, sea ice brine and seawater. Microbial communities in all snow layers derived from mixed sources, both marine and terrestrial, and were more similar to atmospheric communities than to sea ice or seawater communities. The surface snow metagenomes were characterized by the occurrence of genes involved in photochemical stress resistance, primary production and metabolism of diverse carbon sources. The basal saline snow layer that was in direct contact with the sea ice surface harbored a higher abundance of cells than the overlying snow layers, with a predominance of Alteromonadales and a higher relative abundance of marine representatives. However, the overall taxonomic structure of the saline layer was more similar to that of other snow layers and the atmosphere than to underlying sea ice and seawater. The expulsion of relatively nutrient-rich sea ice brine into basal snow might have stimulated the growth of copiotrophic psychro- and halotolerant snow members. Our study indicates that the size, composition and function of snowpack microbial communities over sea ice were influenced primarily by atmospheric deposition and inflow of sea ice brine and that they form a snow-specific assemblage reflecting the particular environmental conditions of the snowpack habitat.

A metagenomic survey of soil microbial communities along a rehabilitation chronosequence after iron ore mining

Microorganisms are useful environmental indicators, able to deliver essential insights to processes regarding mine land rehabilitation. To compare microbial communities from a chronosequence of mine land rehabilitation to pre-disturbance levels from references sites covered by native vegetation, we sampled non-rehabilitated, rehabilitating and reference study sites from the Urucum Massif, Southwestern Brazil. From each study site, three composed soil samples were collected for chemical, physical, and metagenomics analysis. We used a paired-end library sequencing technology (NextSeq 500 Illumina); the reads were assembled using MEGAHIT. Coding DNA sequences (CDS) were identified using Kaiju in combination with non-redundant NCBI BLAST reference sequences containing archaea, bacteria, and viruses. Additionally, a functional classification was performed by EMG v2.3.2. Here, we provide the raw data and assembly (reads and contigs), followed by initial functional and taxonomic analysis, as a base-line for further studies of this kind. Further investigation is needed to fully understand the mechanisms of environmental rehabilitation in tropical regions, inspiring further researchers to explore this collection for hypothesis testing.

Anaerobic degradation of hexadecane and phenanthrene coupled to sulfate reduction by enriched consortia from northern Gulf of Mexico seafloor sediment

To advance understanding of the fate of hydrocarbons released from the Deepwater Horizon oil spill and deposited in marine sediments, this study characterized the microbial populations capable of anaerobic hydrocarbon degradation coupled with sulfate reduction in non-seep sediments of the northern Gulf of Mexico. Anaerobic, sediment-free enrichment cultures were obtained with either hexadecane or phenanthrene as sole carbon source and sulfate as a terminal electron acceptor. Phylogenetic analysis revealed that enriched microbial populations differed by hydrocarbon substrate, with abundant SSU rRNA gene amplicon sequences from hexadecane cultures showing high sequence identity (up to 98%) to Desulfatibacillum alkenivorans (family Desulfobacteraceae), while phenanthrene-enriched populations were most closely related to Desulfatiglans spp. (up to 95% sequence identity; family Desulfarculaceae). Assuming complete oxidation to CO2, observed stoichiometric ratios closely resembled the theoretical ratios of 12.25:1 for hexadecane and 8.25:1 for phenanthrene degradation coupled to sulfate reduction. Phenanthrene carboxylic acid was detected in the phenanthrene-degrading enrichment cultures, providing evidence to indicate carboxylation as an activation mechanism for phenanthrene degradation. Metagenome-assembled genomes (MAGs) revealed that phenanthrene degradation is likely mediated by novel genera or families of sulfate-reducing bacteria along with their fermentative syntrophic partners, and candidate genes linked to the degradation of aromatic hydrocarbons were detected for future study.

Microbial Community Structure and Functional Potential in Cultivated and Native Tallgrass Prairie Soils of the Midwestern United States

The North American prairie covered about 3.6 million-km2 of the continent prior to European contact. Only 1-2% of the original prairie remains, but the soils that developed under these prairies are some of the most productive and fertile in the world, containing over 35% of the soil carbon in the continental United States. Cultivation may alter microbial diversity and composition, influencing the metabolism of carbon, nitrogen, and other elements. Here, we explored the structure and functional potential of the soil microbiome in paired cultivated-corn (at the time of sampling) and never-cultivated native prairie soils across a three-states transect (Wisconsin, Iowa, and Kansas) using metagenomic and 16S rRNA gene sequencing and lipid analysis. At the Wisconsin site, we also sampled adjacent restored prairie and switchgrass plots. We found that agricultural practices drove differences in community composition and diversity across the transect. Microbial biomass in prairie samples was twice that of cultivated soils, but alpha diversity was higher with cultivation. Metagenome analyses revealed denitrification and starch degradation genes were abundant across all soils, as were core genes involved in response to osmotic stress, resource transport, and environmental sensing. Together, these data indicate that cultivation shifted the microbiome in consistent ways across different regions of the prairie, but also suggest that many functions are resilient to changes caused by land management practices - perhaps reflecting adaptations to conditions common to tallgrass prairie soils in the region (e.g., soil type, parent material, development under grasses, temperature and rainfall patterns, and annual freeze-thaw cycles). These findings are important for understanding the long-term consequences of land management practices to prairie soil microbial communities and their genetic potential to carry out key functions.