Taxator-tk: precise taxonomic assignment of metagenomes by fast approximation of evolutionary neighborhoods

Spontaneous fermentation mitigates the frequency of genes encoding antimicrobial resistance spreading from the phyllosphere reservoir to the diet

The phyllosphere microbiome of vegetable products constitutes an important reservoir for multidrug resistant bacteria and Antibiotic Resistance Genes (ARG). Vegetable products including fermented products such as Paocai therefore may serve as a shuttle for extrinsic microorganisms with ARGs into the gut of consumers. Here we study the effect of fermentation on Paocai ARG dissemination by metagenomic analysis. Microbial abundance and diversity of the Paocai microbiome were diminished during fermentation, which correlated with the reduction of abundance in ARGs. Specifically, as fermentation progressed, Enterobacterales overtook Pseudomonadales as the predominant ARG carriers, and Lactobacillales and Enterobacteriales became the determinants of Paocai resistome variation. Moreover, the dual effect of microbes and metal resistance genes (MRGs) was the major contributor driving Paocai resistome dynamics. We recovered several metagenome-assembled genomes (MAGs) carrying acquired ARGs in the phyllosphere microbiome. ARGs of potential clinical and epidemiological relevance such as tet M and emrB-qacA, were mainly hosted by non-dominant bacterial genera. Overall, our study provides evidence that changes in microbial community composition by fermentation aid in constraining ARG dispersal from raw ingredients to the human microbiome but does not eliminate them.

Metagenomes from microbial populations beneath a chromium waste tip give insight into the mechanism of Cr (VI) reduction

Dumped Chromium Ore Processing Residue (COPR) at legacy sites poses a threat to health through leaching of toxic Cr(VI) into groundwater. Previous work implicates microbial activity in reducing Cr(VI) to less mobile and toxic Cr(III), but the mechanism has not been explored. To address this question a combined metagenomic and geochemical study was undertaken. Soil samples from below the COPR waste were used to establish anaerobic microcosms which were challenged with Cr(VI), with or without acetate as an electron donor, and incubated for 70 days. Cr was rapidly reduced in both systems, which also reduced nitrate, nitrite then sulfate, but this sequence was accelerated in the acetate amended microcosms. 16S rRNA gene sequencing revealed that the original soil sample was diverse but both microcosm systems became less diverse by the end of the experiment. A high proportion of 16S rRNA gene reads and metagenome-assembled genomes (MAGs) with high completeness could not be taxonomically classified, highlighting the distinctiveness of these alkaline Cr impacted systems. Examination of the coding capacity revealed widespread capability for metal tolerance and Fe uptake and storage, and both populations possessed metabolic capability to degrade a wide range of organic molecules. The relative abundance of genes for fatty acid degradation was 4× higher in the unamended compared to the acetate amended system, whereas the capacity for dissimilatory sulfate metabolism was 3× higher in the acetate amended system. We demonstrate that naturally occurring in situ bacterial populations have the metabolic capability to couple acetate oxidation to sequential reduction of electron acceptors which can reduce Cr(VI) to less mobile and toxic Cr(III), and that microbially produced sulfide may be important in reductive precipitation of chromate. This capability could be harnessed to create a Cr(VI) trap-zone beneath COPR tips without the need to disturb the waste.

Rich microbial and depolymerising diversity in Antarctic krill gut

With almost a quadrillion individuals, the Antarctic krill processes five million tons of organic carbon every day during austral summer. This high carbon flux requires a broad range of hydrolytic enzymes to decompose the diverse food-derived biopolymers. While krill itself possesses numerous such enzymes, it is unclear, to what extent the endogenous microbiota contribute to the hydrolytic potential of the gut environment. Here we applied amplicon sequencing, shotgun metagenomics, cultivation, and physiological assays to characterize the krill gut microbiota. The broad bacterial diversity (273 families, 919 genera, and 2,309 species) also included a complex potentially anaerobic sub-community. Plate-based assays with 198 isolated pure cultures revealed widespread capacities to utilize lipids (e.g., tributyrin), followed by proteins (casein) and to a lesser extent by polysaccharides (e.g., alginate and chitin). While most isolates affiliated with the genera Pseudoalteromonas and Psychrobacter, also Rubritalea spp. (Verrucomicrobia) were observed. The krill gut microbiota growing on marine broth agar plates possess 13,012 predicted hydrolyses; 15-fold more than previously predicted from a transcriptome-proteome compendium of krill. Cultivation-independent and -dependent approaches indicated members of the families Flavobacteriaceae and Pseudoalteromonadaceae to dominate the capacities for lipid/protein hydrolysis and to provide a plethora of carbohydrate-active enzymes, sulfatases, and laminarin- or porphyrin-depolymerizing hydrolases. Notably, also the potential to hydrolyze plastics such as polyethylene terephthalate and polylactatide was observed, affiliating mostly with Moraxellaceae. Overall, this study shows extensive microbial diversity in the krill gut, and suggests that the microbiota likely play a significant role in the nutrient acquisition of the krill by enriching its hydrolytic enzyme repertoire.IMPORTANCEThe Antarctic krill (Euphausia superba) is a keystone species of the Antarctic marine food web, connecting the productivity of phyto- and zooplankton with the nutrition of the higher trophic levels. Accordingly, krill significantly contributes to biomass turnover, requiring the decomposition of seasonally varying plankton-derived biopolymers. This study highlights the likely role of the krill gut microbiota in this ecosystem function by revealing the great number of diverse hydrolases that microbes contribute to the krill gut environment. The here resolved repertoire of hydrolytic enzymes could contribute to the overall nutritional resilience of krill and to the general organic matter cycling under changing environmental conditions in the Antarctic sea water. Furthermore, the krill gut microbiome could serve as a valuable resource of cold-adapted hydrolytic enzymes for diverse biotechnological applications.

Selective pressure of arsenic and antimony co-contamination on microbial community in alkaline sediments

Although response of microbial community to arsenic (As) and antimony (Sb) co-contamination has been investigated in neutral and acidic environments, little is known in alkaline environment. Herein, the microbial response and survival strategies under the stress of As and Sb co-contamination were determined in the alkaline sediments. Elevated concentrations of As (13700 ± 5012 mg/kg) and Sb (10222 ± 1619 mg/kg) were introduced into the alkaline sediments by the mine drainage, which was partially adopted in the aquatic environment and resulted in a relatively lower contamination (As, 6633 ± 1707 mg/kg; Sb, 6108 ± 1095 mg/kg) in the downstream sediments. The microbial richness was significantly damaged and the microbial compositions were dramatically shifted by the As and Sb co-contamination. Metagenomic analysis shed light on the survival strategies of the microbes under the pressure of As and Sb co-contamination including metal oxidation coupled with denitrification, metal reduction, and metal resistance. The representative microbes were revealed in the sediments with higher (Halomonas) and lower (Thiobacillus, Hydrogenophaga and Flavihumibacter) As and Sb concentration, respectively. In addition, antibiotic resistance genes were found to co-occur with metal resistance genes in the assembled bins. These findings might provide theoretical guidance for bioremediation of As and Sb co-contamination in alkaline environment.

Viral Communities Contribute More to the Lysis of Antibiotic-Resistant Bacteria than the Transduction of Antibiotic Resistance Genes in Anaerobic Digestion Revealed by Metagenomics

Ecological role of the viral community on the fate of antibiotic resistance genes (ARGs) (reduction vs proliferation) remains unclear in anaerobic digestion (AD). Metagenomics revealed a dominance of Siphoviridae and Podoviridae among 13,895 identified viral operational taxonomic units (vOTUs) within AD, and only 21 of the vOTUs carried ARGs, which only accounted for 0.57 ± 0.43% of AD antibiotic resistome. Conversely, ARGs locating on plasmids and integrative and conjugative elements accounted for above 61.0%, indicating a substantial potential for conjugation in driving horizontal gene transfer of ARGs within AD. Virus-host prediction based on CRISPR spacer, tRNA, and homology matches indicated that most viruses (80.2%) could not infect across genera. Among 480 high-quality metagenome assembly genomes, 95 carried ARGs and were considered as putative antibiotic-resistant bacteria (pARB). Furthermore, lytic phages of 66 pARBs were identified and devoid of ARGs, and virus/host abundance ratios with an average value of 71.7 indicated extensive viral activity and lysis. The infectivity of lytic phage was also elucidated through laboratory experiments concerning changes of the phage-to-host ratio, pH, and temperature. Although metagenomic evidence for dissemination of ARGs by phage transduction was found, the higher proportion of lytic phages infecting pARBs suggested that the viral community played a greater role in reducing ARB numbers than spreading ARGs in AD.

Soil Humic Acid Stimulates Potentially Active Dissimilatory Arsenate-Reducing Bacteria in Flooded Paddy Soil as Revealed by Metagenomic Stable Isotope Probing

Dissimilatory arsenate reduction contributes a large proportion of arsenic flux from flooded paddy soil, which is closely linked to soil organic carbon input and efflux. Humic acid (HA) represents a natural ingredient in soil and is shown to enhance microbial arsenate respiration to promote arsenic mobility. However, the community and function profiles of metabolically active arsenate-respiring bacteria and their interactions with HA in paddy soil remain unclear. To probe this linkage, we performed a genome-centric comparison of potentially active arsenate-respiring bacteria in anaerobic microcosms amended with 13C-lactate and HA by combining stable-isotope probing with genome-resolved metagenomics. Indeed, HA greatly accelerated the microbial reduction of arsenate to arsenite. Enrichment of bacteria that harbor arsenate-respiring reductase genes (arrA) in HA-enriched 13C-DNA was confirmed by metagenomic binning, which are affiliated with Firmicutes (mainly Desulfitobacterium, Bacillus, Brevibacillus, and Clostridia) and Acidobacteria. Characterization of reference extracellular electron transfer (EET)-related genes in these arrA-harboring bacteria supports the presence of EET-like genes, with partial electron-transport chain genes identified. This suggests that Gram-positive Firmicutes- and Acidobacteria-related members may harbor unspecified EET-associated genes involved in metal reduction. Our findings highlight the link between soil HA and potentially active arsenate-respiring bacteria, which can be considered when using HA for arsenic removal.

Fecal microbial gene transfer contributes to the high-grain diet-induced augmentation of aminoglycoside resistance in dairy cattle

A high-grain (HG) diet can rapidly lower the rumen pH and thus modify the gastrointestinal microbiome in dairy cattle. Although the prevalence of antibiotic resistance is strongly linked with the gut microbiome, the influences of HG diet on animals' gut resistome remain largely unexplored. Here, we examined the impact and mechanism of an HG diet on the fecal resistome in dairy cattle by metagenomically characterizing the gut microbiome. Eight lactating Holstein cattle were randomly allocated into two groups and fed either a conventional (CON) or HG diet for 3 weeks. The fecal microbiome and resistome were significantly altered in dairy cattle from HG, demonstrating an adaptive response that peaks at day 14 after the dietary transition. Importantly, we determined that feeding an HG diet specifically elevated the prevalence of resistance to aminoglycosides (0.11 vs 0.24 RPKG, P < 0.05). This diet-induced resistance increase is interrelated with the disproportional propagation of microbes in Lachnospiraceae, indicating a potential reservoir of aminoglycosides resistance. We further showed that the prevalence of acquired resistance genes was also modified by introducing a different diet, likely due to the augmented frequency of lateral gene transfer (LGT) in microbes (CON vs HG: 254 vs 287 taxa) such as Lachnospiraceae. Consequently, we present that diet transition is associated with fecal resistome modification in dairy cattle and an HG diet specifically enriched aminoglycosides resistance that is likely by stimulating microbial LGT.IMPORTANCEThe increasing prevalence of antimicrobial resistance is one of the most severe threats to public health, and developing novel mitigation strategies deserves our top priority. High-grain (HG) diet is commonly applied in dairy cattle to enhance animals' performance to produce more high-quality milk. We present that despite such benefits, the application of an HG diet is correlated with an elevated prevalence of resistance to aminoglycosides, and this is a combined effect of the expansion of antibiotic-resistant bacteria and increased frequency of lateral gene transfer in the fecal microbiome of dairy cattle. Our results provided new knowledge in a typically ignored area by showing an unexpected enrichment of antibiotic resistance under an HG diet. Importantly, our findings laid the foundation for designing potential dietary intervention strategies to lower the prevalence of antibiotic resistance in dairy production.

Potential of Microbial Communities to Perform Dehalogenation Processes in Natural and Anthropogenically Modified Environments-A Metagenomic Study

Halogenated organic compounds (HOCs) pose a serious problem for the environment. Many are highly toxic and accumulate both in soil and in organisms. Their biological transformation takes place by dehalogenation, in which the halogen substituents are detached from the carbon in the organic compound by enzymes produced by microorganisms. This increases the compounds' water solubility and bioavailability, reduces toxicity, and allows the resulting compound to become more susceptible to biodegradation. The microbial halogen cycle in soil is an important part of global dehalogenation processes. The aim of the study was to examine the potential of microbial communities inhabiting natural and anthropogenically modified environments to carry out the dehalogenation process. The potential of microorganisms was assessed by analyzing the metagenomes from a natural environment (forest soils) and from environments subjected to anthropopression (agricultural soil and sludge from wastewater treatment plants). Thirteen genes encoding enzymes with dehalogenase activity were identified in the metagenomes of both environments, among which, 2-haloacid dehalogenase and catechol 2,3-dioxygenase were the most abundant genes. Comparative analysis, based on comparing taxonomy, identified genes, total halogens content and content of DDT derivatives, demonstrated the ability of microorganisms to transform HOCs in both environments, indicating the presence of these compounds in the environment for a long period of time and the adaptive need to develop mechanisms for their detoxification. Metagenome analyses and comparative analyses indicate the genetic potential of microorganisms of both environments to carry out dehalogenation processes, including dehalogenation of anthropogenic HOCs.

Virus impacted community adaptation in oligotrophic groundwater environment revealed by Hi-C coupled metagenomic and viromic study

Viruses play a crucial role in microbial mortality, diversity and biogeochemical cycles. Groundwater is the largest global freshwater and one of the most oligotrophic aquatic systems on Earth, but how microbial and viral communities are shaped in this special habitat is largely unexplored. In this study, we collected groundwater samples from 23 to 60 m aquifers at Yinchuan Plain, China. In total, 1920 non-reductant viral contigs were retrieved from metagenomes and viromes constructed by Illumina and Nanopore hybrid sequencing. Only 3% of them could be clustered with known viruses, most of which were Caudoviricetes. Coupling 1.2 Tb Hi-C sequencing with CRISPR matching and homology search, we connected 469 viruses with their hosts while some viral clusters presented a broad-host-range trait. Meanwhile, a large proportion of biosynthesis related auxiliary metabolism genes were identified. Those characteristics might benefit viruses for a better survival in this special oligotrophic environment. Additionally, the groundwater virome showed genomic features distinct from those of the open ocean and wastewater treatment facilities in GC distribution and unannotated gene compositions. This paper expands the current knowledge of the global viromic records and serves as a foundation for a more thorough understanding of viruses in groundwater.

Gut microbiota enhance energy accumulation of black-necked crane to cope with impending migration

Less is known about the role of gut microbiota in overwintering environmental adaptation in migratory birds. Here, we performed metagenomic sequencing on fresh fecal samples (n = 24) collected during 4 periods of overwintering (Dec: early; Jan: middle I; Feb: middle II; Mar: late) to characterize gut microbial taxonomic and functional characteristics of black-necked crane (Grus nigricollis). The results demonstrated no significant change in microbial diversity among overwintering periods. Analysis of compositions of microbiomes with bias correction (ANCOM-BC) determined 15 Proteobacteria species enriched in late overwintering period. Based on previous reports, these species are associated with degradation of chitin, cellulose, and lipids. Meanwhile, fatty acid degradation and betalain biosynthesis pathways are enriched in late overwintering period. Furthermore, metagenomic binning obtained 91 high-quality bins (completeness >70% and contamination <10%), 5 of which enriched in late overwintering period. Carnobacterium maltaromaticum, unknown Enterobacteriaceae, and Yersinia frederiksenii have genes for chitin and cellulose degradation, acetate, and glutamate production. Unknown Enterobacteriaceae and Y. frederiksenii hold genes for synthesis of 10 essential amino acids required by birds, and the latter has genes for γ-aminobutyrate production. C. maltaromaticum has genes for pyridoxal synthesis. These results implied the gut microbiota is adapted to the host diet and may help black-necked cranes in pre-migratory energy accumulation by degrading the complex polysaccharide in their diet, supplying essential amino acids and vitamin pyridoxal, and producing acetate, glutamate, and γ-aminobutyrate that could stimulate host feeding. Additionally, enriched Proteobacteria also encoded more carbohydrate-active enzymes (CAZymes) and antibiotic resistance genes (ARGs) in late overwintering period. KEY POINTS: • Differences in gut microbiota function during overwintering period of black-necked cranes depend mainly on changes in core microbiota abundance • Gut microbiota of black-necked crane adapted to the diet during overwintering period • Gut microbiota could help black-necked cranes to accumulate more energy in the late overwintering period.

Microbial biogeochemical cycling reveals the sustainability of the rice-crayfish co-culture model

Aquaculture has great potential in nourishing the global growing population, while such staggering yields are coupled with environmental pollution. Rice-crayfish co-culture models (RCFP) have been widely adopted in China due to their eco-friendliness. However, little is known about RCFP's microbiome pattern, which hinders our understanding of its sustainability. This study has conducted metagenomic analysis across aquaculture models and habitats, which revealed aquaculture model-specific biogeochemical cycling pattern (e.g., nitrogen (N), sulfur (S), and carbon (C)): RCFP is advantageous in N-assimilation, N-contamination, and S-pollutants removal, while non-RCFP features N denitrification process and higher S metabolism ability, producing several hazardous pollutants in non-RCFP (e.g., nitric oxide, nitrogen monoxide, and sulfide). Moreover, RCFP has greater capacity for carbohydrate enzyme metabolism compared with non-RCFP in environmental habitats, but not in crayfish gut. Collectively, RCFP plays an indispensable role in balancing aquaculture productivity and environmental protection, which might be applied to the blue transformation of aquaculture.

DNA virome composition of two sympatric wild felids, bobcat (Lynx rufus) and puma (Puma concolor) in Sonora, Mexico

With viruses often having devastating effects on wildlife population fitness and wild mammals serving as pathogen reservoirs for potentially zoonotic diseases, determining the viral diversity present in wild mammals is both a conservation and One Health priority. Additionally, transmission from more abundant hosts could increase the extinction risk of threatened sympatric species. We leveraged an existing circular DNA enriched metagenomic dataset generated from bobcat (Lynx rufus, n = 9) and puma (Puma concolor, n = 13) scat samples non-invasively collected from Sonora, Mexico, to characterize fecal DNA viromes of each species and determine the extent that viruses are shared between them. Using the metaWRAP pipeline to co-assemble viral genomes for comparative metagenomic analysis, we observed diverse circular DNA viruses in both species, including circoviruses, genomoviruses, and anelloviruses. We found that differences in DNA virome composition were partly attributed to host species, although there was overlap between viruses in bobcats and pumas. Pumas exhibited greater levels of alpha diversity, possibly due to bioaccumulation of pathogens in apex predators. Shared viral taxa may reflect dietary overlap, shared environmental resources, or transmission through host interactions, although we cannot rule out species-specific host-virus coevolution for the taxa detected through co-assembly. However, our detection of integrated feline foamy virus (FFV) suggests Sonoran pumas may interact with domestic cats. Our results contribute to the growing baseline knowledge of wild felid viral diversity. Future research including samples from additional sources (e.g., prey items, tissues) may help to clarify host associations and determine the pathogenicity of detected viruses.

Recovery of 197 eukaryotic bins reveals major challenges for eukaryote genome reconstruction from terrestrial metagenomes

As most eukaryotic genomes are yet to be sequenced, the mechanisms underlying their contribution to different ecosystem processes remain untapped. Although approaches to recovering Prokaryotic genomes have become common in genome biology, few studies have tackled the recovery of eukaryotic genomes from metagenomes. This study assessed the reconstruction of microbial eukaryotic genomes using 6000 metagenomes from terrestrial and some transition environments using the EukRep pipeline. Only 215 metagenomic libraries yielded eukaryotic bins. From a total of 447 eukaryotic bins recovered 197 were classified at the phylum level. Streptophytes and fungi were the most represented clades with 83 and 73 bins, respectively. More than 78% of the obtained eukaryotic bins were recovered from samples whose biomes were classified as host-associated, aquatic, and anthropogenic terrestrial. However, only 93 bins were taxonomically assigned at the genus level and 17 bins at the species level. Completeness and contamination estimates were obtained for a total of 193 bins and consisted of 44.64% (σ = 27.41%) and 3.97% (σ = 6.53%), respectively. Micromonas commoda was the most frequent taxon found while Saccharomyces cerevisiae presented the highest completeness, probably because more reference genomes are available. Current measures of completeness are based on the presence of single-copy genes. However, mapping of the contigs from the recovered eukaryotic bins to the chromosomes of the reference genomes showed many gaps, suggesting that completeness measures should also include chromosome coverage. Recovering eukaryotic genomes will benefit significantly from long-read sequencing, development of tools for dealing with repeat-rich genomes, and improved reference genomes databases.

Monthly dynamics of microbial communities and variation of nitrogen-cycling genes in an industrial-scale expanded granular sludge bed reactor

Introduction

The expanded granular sludge bed (EGSB) is a major form of anaerobic digestion system during wastewater treatment. Yet, the dynamics of microbial and viral communities and members functioning in nitrogen cycling along with monthly changing physicochemical properties have not been well elucidated.

Methods

Here, by collecting the anaerobic activated sludge samples from a continuously operating industrial-scale EGSB reactor, we conducted 16S rRNA gene amplicon sequencing and metagenome sequencing to reveal the microbial community structure and variation with the ever-changing physicochemical properties along within a year.

Results

We observed a clear monthly variation of microbial community structures, while COD, the ratio of volatile suspended solids (VSS) to total suspended solids (TSS) (VSS/TSS ratio), and temperature were predominant factors in shaping community dissimilarities examined by generalized boosted regression modeling (GBM) analysis. Meanwhile, a significant correlation was found between the changing physicochemical properties and microbial communities (p <0.05). The alpha diversity (Chao1 and Shannon) was significantly higher (p <0.05) in both winter (December, January, and February) and autumn (September, October, and November) with higher organic loading rate (OLR), higher VSS/TSS ratio, and lower temperature, resulting higher biogas production and nutrition removal efficiency. Further, 18 key genes covering nitrate reduction, denitrification, nitrification, and nitrogen fixation pathways were discovered, the total abundance of which was significantly associated with the changing environmental factors (p <0.05). Among these pathways, the dissimilatory nitrate reduction to ammonia (DNRA) and denitrification had the higher abundance contributed by the top highly abundant genes narGH, nrfABCDH, and hcp. The COD, OLR, and temperature were primary factors in affecting DNRA and denitrification by GBM evaluation. Moreover, by metagenome binning, we found the DNRA populations mainly belonged to Proteobacteria, Planctomycetota, and Nitrospirae, while the denitrifying bacteria with complete denitrification performance were all Proteobacteria. Besides, we detected 3,360 non-redundant viral sequences with great novelty, in which Siphoviridae, Podoviridae, and Myoviridae were dominant viral families. Interestingly, viral communities likewise depicted clear monthly variation and had significant associations with the recovered populations (p <0.05).

Discussion

Our work highlights the monthly variation of microbial and viral communities during the continuous operation of EGSB affected by the predominant changing COD, OLR, and temperature, while DNRA and denitrification pathways dominated in this anaerobic system. The results also provide a theoretical basis for the optimization of the engineered system.

Gut microbiota contributes to lignocellulose deconstruction and nitrogen fixation of the larva of Apriona swainsoni

Apriona swainsoni is a vital forest pest prevalent in China. The larvae of A. swainsoni live solely in the branches of trees and rely entirely on the xylem for nutrition. However, there is still a lack of in-depth research on the gut microbiota's use of almost nitrogen-free wood components to provide bio-organic macromolecular components needed for their growth. Thus, in this study, the metagenome, metaproteome, and metabolome of the A. swainsoni larvae in four gut segments (foregut; midgut; anterior hindgut; posterior hindgut) were analyzed by the multi-omics combined technology, to explore the metabolic utilization mechanism of the corresponding gut microbiota of A. swainsoni. Firstly, we found that the metagenome of different gut segments was not significantly different in general, but there were different combinations of dominant bacteria and genes in different gut segments, and the metaproteome and metabolome of four gut segments were significantly different in general. Secondly, the multi-omics results showed that there were significant gradient differences in the contents of cellulose and hemicellulose in different segments of A. swainsoni, and the expression of corresponding metabolic proteins was the highest in the midgut, suggesting the metabolic characteristics of these lignocellulose components in A. swainsoni gut segments. Finally, we found that the C/N ratio of woody food was significantly lower than that of frass, and metagenomic results showed that nitrogen fixation genes mainly existed in the foregut and two hindgut segments. The expression of the key nitrogen fixing gene nifH occurred in two hindgut parts, indicating the feature of nitrogen fixation of A. swainsoni. In conclusion, our results provide direct evidence that the larvae of A. swainsoni can adapt to the relatively harsh niche conditions through the highly organized gut microbiome in four gut segments, and may play a major role in their growth.

Binning Metagenomic Contigs Using Unsupervised Clustering and Reference Databases

Metagenomics can directly extract the genetic material of all microorganisms from the environment, and obtain metagenomic samples with a large number of unknown DNA sequences. Binning of metagenomic contigs is a hot topic in metagenomics research. There are two key challenges for the current unsupervised metagenomic clustering algorithms. First, unsupervised metagenomic clustering methods rarely use reference databases, causing a certain waste of resources. Second, unsupervised metagenomic clustering methods are restricted by the characteristics of the sequences and the clustering algorithms, and the binning effect is limited. Therefore, a new binning method for metagenomic contigs using unsupervised clustering methods and reference databases is proposed to address these challenges, to make full use of the advantages of unsupervised clustering methods and reference databases constructed by scientists to improve the overall binning effect. This method uses the integrated SVM classification model to further bin the unsupervised clustering parts that do not perform well. Our proposed method was tested on simulated datasets and a real dataset and compared with other state-of-the-art metagenomic clustering methods including CONCOCT, Metabin2.0, Autometa, and MetaBAT. The results show that our method can achieve higher precision rate and improve the binning effect.

A preliminary study on the feasibility of industrialization for n-caproic acid recovery from food wastewater: From lab to pilot

Food wastewater is associated with greenhouse gas emission and has a significant water footprint. Here, the platform chemical n-caproate was recovered from liquor brewing wastewater at maximum and mean concentrations of 26.4 g/L and 17.0 ± 4.3 g/L, respectively, after 377 d operation. Laboratory-scale lactate-driven chain elongation (CE) process was implemented first. Taxonomic composition and metagenomic data analyses revealed that Caproiciproducens (e.g., Ruminococcaceae bacterium CPB6) and bacteria affiliated with Lachnospiraceae transformed lactate to n-caproate by reverse β-oxidation and/or fatty acid biosynthesis. The lactate-driven CE process was then scaled up from 2.5 L to 500 L and achieved a n-caproate production of 14.5 ± 0.6 g/L within 96 h. n-Caproic acid was extracted at a concentration and purity of 815.9 ± 8.3 g/L and 88.6 ± 8.9 %, respectively. The present study demonstrated a commercially viable strategy for resource recovery and carbon fixation from food waste streams.

Characterization of arsenic-metabolizing bacteria in an alkaline soil

Arsenite (As(III)) is more toxic, mobilizable and bioavailable than arsenate (As(V)). Hence, the transformations between As(III) and As(V) are crucial for the toxicity and mobility of arsenic (As). However, As transformation and microbial communities involved in alkaline soils are largely unknown. Here we investigate two major pathways of As transformation, i.e., As(III) oxidation and As(V) reduction, and identify the bacteria involved in the alkaline soil by combining stable isotope probing with shotgun metagenomic sequencing. As(III) oxidation and significant increase of the aioA genes copies were observed in the treatments amended with As(III) and NO₃^-, suggesting that As(III) oxidation can couple with nitrate reduction and was mainly catalyzed by the microorganisms containing aioA genes. As(V) reduction was detected in the treatments amended with As(V) and acetate where the abundance of arrA gene significantly increased, indicating that microorganisms with arrA genes were the key As(V) reducers. Acidovorax, Hydrogenophaga, and Ramlibacter were the putative nitrate-dependent As(III) oxidizers, and Deinococcus and Serratia were the putative respiratory As(V) reducers. These findings will improve our understanding of As metabolism and are meaningful for mapping out bioremediation strategies of As contamination in alkaline environment.

Understanding of mercury and methylmercury transformation in sludge composting by metagenomic analysis

Municipal sewage especially the produced sewage sludge is a significant source releasing mercury (Hg) to the environment. However, the Hg speciation especially methylmercury (MeHg) transformation in sewage sludge treatment process remains poorly understood. This study investigated the transformation of Hg speciation especially MeHg in sludge composting. The distribution of Hg transformation related gene pairs hgcAB and merAB, and their putative microbial hosts were comprehensively analyzed. Both Hg (from 3.16±0.22 mg/kg to 3.20±0.19 mg/kg) and MeHg content (from 4.77±0.64 ng/g to 4.36±0.37 ng/g) were not obviously changed before and after composting, but about 19.69% of Hg and 27.36% of MeHg were lost according to mass balance calculation. The metagenomic analysis further revealed that anaerobes (Desulfobacterota and Euryarchaeota) were the mainly putative Hg methylators especially carrying high abundance of hgcA gene in the initial periods of composting. Among the 151 reconstructed metagenome-assembled genomes (MAGs), only 4 hgcA gene carriers (Myxococcota, Firmicutes, Cyclobacteriaceae, and Methanothermobacter) and 16 merB gene carriers were identified. But almost all of the MAGs carried hgcB gene and merA gene. The merA gene was widely distributed in genomes, which indicated the widespread functionality of microbes for reducing Hg(II) to Hg(0). The hgcA carrying microbes tends to present the similar metabolic pathways including methanogenesis and sulfur metabolism. Besides, both the irregular distribution of hgcA in various species (including Actinobacteria, Archaea, Bacteroidetes, Desulfobacterota, Euryarchaeota, and Nitrospirae, etc.) and opposite evolution trends between hgcA gene abundance and its host genome abundance can be an indication of horizontal gene transfer or gene deletions of hgcA during composting. Our findings thus revealed that sludge composting is not only a hotspot for Hg speciation transformation, but also a potential hotspot for MeHg transformation.

Shotgun metagenomics of fecal samples from children in Peru reveals frequent complex co-infections with multiple Campylobacter species

Campylobacter spp. are a major cause of bacterial diarrhea worldwide and are associated with high rates of mortality and linear growth faltering in children living in low- to middle-income countries (LMICs). Campylobacter jejuni and Campylobacter coli are most often the causative agents of enteric disease among children in LMICs. However, previous work on a collection of stool samples from children under 2 years of age, living in a low resource community in Peru with either acute diarrheal disease or asymptomatic, were found to be qPCR positive for Campylobacter species but qPCR negative for C. jejuni and C. coli. The goal of this study was to determine if whole-genome shotgun metagenomic sequencing (WSMS) could identify the Campylobacter species within these samples. The Campylobacter species identified in these stool samples included C. jejuni, C. coli, C. upsaliensis, C. concisus, and the potential new species of Campylobacter, "Candidatus Campylobacter infans". Moreover, WSMS results demonstrate that over 65% of the samples represented co-infections with multiple Campylobacter species present in a single stool sample, a novel finding in human populations.

Analysis of shotgun metagenomic data obtained from fecal samples of children living in a low resource tropical community of Peru revealed multiple Campylobacter species. Co-infections with more than one Campylobacter species within the same sample was a common finding. A potential new species of Campylobacter was also detected within these samples.

Acetate reprograms gut microbiota during alcohol consumption

Liver damage due to chronic alcohol use is among the most prevalent liver diseases. Alcohol consumption frequency is a strong factor of microbiota variance. Here we use isotope labeled [1-13C] ethanol, metagenomics, and metatranscriptomics in ethanol-feeding and intragastric mouse models to investigate the metabolic impacts of alcohol consumption on the gut microbiota. First, we show that although stable isotope labeled [1-13C] ethanol contributes to fatty acid pools in the liver, plasma, and cecum contents of mice, there is no evidence of ethanol metabolism by gut microbiota ex vivo under anaerobic conditions. Next, we observe through metatranscriptomics that the gut microbiota responds to ethanol-feeding by activating acetate dissimilation, not by metabolizing ethanol directly. We demonstrate that blood acetate concentrations are elevated during ethanol consumption. Finally, by increasing systemic acetate levels with glyceryl triacetate supplementation, we do not observe any impact on liver disease, but do induce similar gut microbiota alterations as chronic ethanol-feeding in mice. Our results show that ethanol is not directly metabolized by the gut microbiota, and changes in the gut microbiota linked to ethanol are a side effect of elevated acetate levels. De-trending for these acetate effects may be critical for understanding gut microbiota changes that cause alcohol-related liver disease.

Metagenomics: A Tool for Exploring Key Microbiome With the Potentials for Improving Sustainable Agriculture

Microorganisms are immense in nature and exist in every imaginable ecological niche, performing a wide range of metabolic processes. Unfortunately, using traditional microbiological methods, most microorganisms remain unculturable. The emergence of metagenomics has resolved the challenge of capturing the entire microbial community in an environmental sample by enabling the analysis of whole genomes without requiring culturing. Metagenomics as a non-culture approach encompasses a greater amount of genetic information than traditional approaches. The plant root-associated microbial community is essential for plant growth and development, hence the interactions between microorganisms, soil, and plants is essential to understand and improve crop yields in rural and urban agriculture. Although some of these microorganisms are currently unculturable in the laboratory, metagenomic techniques may nevertheless be used to identify the microorganisms and their functional traits. A detailed understanding of these organisms and their interactions should facilitate an improvement of plant growth and sustainable crop production in soil and soilless agriculture. Therefore, the objective of this review is to provide insights into metagenomic techniques to study plant root-associated microbiota and microbial ecology. In addition, the different DNA-based techniques and their role in elaborating plant microbiomes are discussed. As an understanding of these microorganisms and their biotechnological potentials are unlocked through metagenomics, they can be used to develop new, useful and unique bio-fertilizers and bio-pesticides that are not harmful to the environment.

Analysis of antibiotic resistance genes reveals their important roles in influencing the community structure of ocean microbiome

Antibiotic resistance gene (ARG) content is a well-established driver of microbial abundance and diversity in an environment. By reanalyzing 132 metagenomic datasets from the Tara Oceans project, we aim to unveil the associations between environmental factors, the ocean microbial community structure and ARG contents. We first investigated the structural patterns of microbial communities including both prokaryotes such as bacteria and eukaryotes such as protists. Additionally, several ARG-dominant horizontal gene transfer events between Protist and Prokaryote have been identified, indicating the potential roles of ARG in shaping the ocean microbial communities. For a deeper insight into the role of ARGs in ocean microbial communities on a global scale, we identified 1926 unique types of ARGs and discovered that the ARGs are more abundant and diverse in the mesopelagic zone than other water layers, potentially caused by limited resources. Finally, we found that ARG-enriched genera were often more abundant compared to their ARG-less neighbors in the same environment (e.g. coastal oceans). A deeper understanding of the ARG-microbiome relationships could help in the conservation of the oceanic ecosystem.

Potential Use of Microbial Community Genomes in Various Dimensions of Agriculture Productivity and Its Management: A Review

Agricultural productivity is highly influenced by its associated microbial community. With advancements in omics technology, metagenomics is known to play a vital role in microbial world studies by unlocking the uncultured microbial populations present in the environment. Metagenomics is a diagnostic tool to target unique signature loci of plant and animal pathogens as well as beneficial microorganisms from samples. Here, we reviewed various aspects of metagenomics from experimental methods to techniques used for sequencing, as well as diversified computational resources, including databases and software tools. Exhaustive focus and study are conducted on the application of metagenomics in agriculture, deciphering various areas, including pathogen and plant disease identification, disease resistance breeding, plant pest control, weed management, abiotic stress management, post-harvest management, discoveries in agriculture, source of novel molecules/compounds, biosurfactants and natural product, identification of biosynthetic molecules, use in genetically modified crops, and antibiotic-resistant genes. Metagenomics-wide association studies study in agriculture on crop productivity rates, intercropping analysis, and agronomic field is analyzed. This article is the first of its comprehensive study and prospects from an agriculture perspective, focusing on a wider range of applications of metagenomics and its association studies.

Metagenomics insights into the profiles of antibiotic resistome in combined sewage overflows from reads to metagenome assembly genomes

Combined sewage overflows (CSOs) have become an important source of antibiotic resistance genes (ARGs) in the environment, while the distribution and dynamics of antibiotic resistome in the CSOs events have not been well understood. This study deciphered the profiles of antibiotic resitome in the CSOs based on metagenomics analysis from reads to metagenome assembly genomes (MAGs), and the dynamical changes of ARGs were clarified through continuous monitoring of the CSO event. Results showed that antibiotic inactivation was the dominant resistance mechanism, and sulfonamide, aminoglycoside along with multidrug resistance were the dominant antibiotic resistance types. It was speculated that the antibiotic resistome were generally determined by sewer sediment flushed out along with the CSOs not domestic sewage in the pipes. The host range and mobility of the antibiotic resistome were determined at contigs level, and the hosts mainly belonged to the Proteobacteria with the genus of Pseudomonas, Escherichia, Enterobacter and Aeromonas being dominant. The transposase (tnpA), IS91 and integrons were mobile genetic elements (MGEs) located together with ARGs, and a MAG carrying 32 ARGs and 140 VFGs was assembled. Although microbial community contributed most to the changes of antibiotic resistome in the CSOs directly, the risks caused by the MGEs should be paid more attention.

Metagenomics revealed the mobility and hosts of antibiotic resistance genes in typical pesticide wastewater treatment plants

Pesticide showed a crucial selective pressure of antibiotic resistance genes (ARGs) in the environmental dimension, especially in the pesticide wastewater treatment process, where the information on the mobility and hosts of ARGs was very important but limited. This study tried to clarify the mobile antibiotic resistome and ARG hosts in three typical pesticide wastewater treatment plants (PWWTPs) through metagenomics. Results showed that ARGs associated with antibiotic efflux and multi-drug resistance generally dominated in the PWWTPs, and the relative abundance of ARGs was generally higher in the water phase than that in sludge phase. The mobile antibiotic resistome accounted for 43.6% ± 16.2% and 44.8% ± 18.0% of the total relative abundance of ARGs in the water phase and sludge phase, respectively. The tnpA, IS91 and intI1 were the dominant mobile genetic elements (MGEs) closely associated with ARGs. MCR-5 and MCR-9 were first identified in the PWWTPs and located together with the tnpA, tnpA2 and int2. The potential human pathogens belonging to Citrobacter, Pseudomonas, Enterobacter, Acinetobacter, and Kluyvern were the major ARG hosts in the PWWTPs. Statistical analysis indicated that microbial community contributed the most to the occurrence of antibiotic resistome, and the reduction of the major ARG hosts was crucial from the perspective of ARGs control.

The phyllosphere microbiome shifts toward combating melanose pathogen

BACKGROUND

Plants can recruit beneficial microbes to enhance their ability to defend against pathogens. However, in contrast to the intensively studied roles of the rhizosphere microbiome in suppressing plant pathogens, the collective community-level change and effect of the phyllosphere microbiome in response to pathogen invasion remains largely elusive.

RESULTS

Here, we integrated 16S metabarcoding, shotgun metagenomics and culture-dependent methods to systematically investigate the changes in phyllosphere microbiome between infected and uninfected citrus leaves by Diaporthe citri, a fungal pathogen causing melanose disease worldwide. Multiple microbiome features suggested a shift in phyllosphere microbiome upon D. citri infection, highlighted by the marked reduction of community evenness, the emergence of large numbers of new microbes, and the intense microbial network. We also identified the microbiome features from functional perspectives in infected leaves, such as enriched microbial functions for iron competition and potential antifungal traits, and enriched microbes with beneficial genomic characteristics. Glasshouse experiments demonstrated that several bacteria associated with the microbiome shift could positively affect plant performance under D. citri challenge, with reductions in disease index ranging from 65.7 to 88.4%. Among them, Pantoea asv90 and Methylobacterium asv41 identified as "recruited new microbes" in the infected leaves, exhibited antagonistic activities to D. citri both in vitro and in vivo, including inhibition of spore germination and/or mycelium growth. Sphingomonas spp. presented beneficial genomic characteristics and were found to be the main contributor for the functional enrichment of iron complex outer membrane receptor protein in the infected leaves. Moreover, Sphingomonas asv20 showed a stronger suppression ability against D. citri in iron-deficient conditions than iron-sufficient conditions, suggesting a role of iron competition during their antagonistic action.

CONCLUSIONS

Overall, our study revealed how phyllosphere microbiomes differed between infected and uninfected citrus leaves by melanose pathogen, and identified potential mechanisms for how the observed microbiome shift might have helped plants cope with pathogen pressure. Our findings provide novel insights into understanding the roles of phyllosphere microbiome responses during pathogen challenge. Video abstract.

Complete Genome Sequence of Aggregatibacter actinomycetemcomitans Strain CU1000N

Here, we report the complete genome sequence of Aggregatibacter actinomycetemcomitans strain CU1000N. This rough strain is used extensively as a model organism to characterize localized aggressive periodontitis pathogenesis, the basic biology and oral cavity colonization of A. actinomycetemcomitans, and its interactions with other members of the oral microbiome.

Evidence for Long-Term Anthropogenic Pollution: The Hadal Trench as a Depository and Indicator for Dissemination of Antibiotic Resistance Genes

Knowledge of the distribution and dissemination of antibiotic resistance genes (ARGs) is essential for understanding anthropogenic impacts on natural ecosystems. The transportation of ARGs via aquatic environments is significant and has received great attention, but whether there has been anthropogenic ARG pollution to the hadal ocean ecosystem has not been well explored. For investigating ecological health concerns, we profiled the ARG occurrence in sediments of the Mariana Trench (MT) (10 890 m), the deepest region of the ocean. Metagenomic-based ARG profiles showed a sudden increase of abundance and diversity in the surface layer of MT sediments reaching 2.73 × 10^-2 copy/cell and 81 subtypes, and a high percentage of ∼63.6% anthropogenic pollution sources was predicted by the Bayesian-modeling classification method. These together suggested that ARG accumulation and anthropogenic impacts have already permeated into the bottom of the deepest corner on the earth. Moreover, six ARG-carrying draft genomes were retrieved using a metagenomic binning strategy, one of which assigned as Streptococcus was identified as a potential bacterial host to contribute to the ARG accumulation in MT, carrying ermF, tetM, tetQ, cfxA2, PBP-2X, and PBP-1A. We propose that the MT ecosystem needs further long-term monitoring for the assessment of human impacts, and our identified three biomarkers (cfxA2, ermF, and mefA) could be used for the rapid monitoring of anthropogenic pollution. Together our findings imply that anthropogenic pollution has penetrated into the deepest region of the ocean and urge for better pollution control to reduce the risk of ARG dissemination to prevent the consistent accumulation and potential threat to the natural environment.

Transition of Serotype 35B Pneumococci From Commensal to Prevalent Virulent Strain in Children

In our community-based prospective cohort study in young children, we observed a significant increase in pneumococcal serotype 35B nasopharyngeal (NP) commensal colonization during the 2011–2014 timeframe, but these strains were not associated with disease. Beginning in 2015 and continuing through to the present, the serotype 35B virulence changed, and it became the dominant bacteria isolated and associated with pneumococcal acute otitis-media (AOM) in our cohort. We performed comparative analyses of 250 35B isolates obtained from 140 children collected between 2006 and 2019. Changes in prevalence, clonal-complex composition, and antibiotic resistance were analyzed. Seventy-two (29%) of 35B isolates underwent whole-genome sequencing to investigate genomic changes associated with the shift in virulence that resulted in increased rates of 35B-associated AOM disease. 35B strains that were commensals and AOM disease-causing were mainly associated with sequence type (ST) 558. Antibiotic concentrations of β-lactams and ofloxacin necessary to inhibit growth of 35B strains rose significantly (2006–2019) (p<0.005). However, only isolates from the 35B/ST558 showed significant increases in MIC₅₀ of penicillin and ofloxacin between the years 2006–2014 and 2015–2019 (p=0.007 and p<0.0001). One hundred thirty-eight SNPs located in 34 different genes were significantly associated with post-2015 strains. SNPs were found in nrdG (metal binding, 10%); metP and metN (ABC transporter, 9%); corA (Mg²⁺ transporter, 6%); priA (DNA replication, 5%); and on the enzymic gene ldcB (LD-carboxypeptidase, 3%). Pneumococcal serotype 35B strains was a common NP commensal during 2010–2014. In 2015, a shift in increasing number of AOM cases occurred in young children caused by 35B, that was associated with changes in genetic composition and antibiotic susceptibility.

Tiara: deep learning-based classification system for eukaryotic sequences

MOTIVATION

With a large number of metagenomic datasets becoming available, eukaryotic metagenomics emerged as a new challenge. The proper classification of eukaryotic nuclear and organellar genomes is an essential step toward a better understanding of eukaryotic diversity.

RESULTS

We developed Tiara, a deep-learning-based approach for the identification of eukaryotic sequences in the metagenomic datasets. Its two-step classification process enables the classification of nuclear and organellar eukaryotic fractions and subsequently divides organellar sequences into plastidial and mitochondrial. Using the test dataset, we have shown that Tiara performed similarly to EukRep for prokaryotes classification and outperformed it for eukaryotes classification with lower calculation time. In the tests on the real data, Tiara performed better than EukRep in analyzing the small dataset representing eukaryotic cell microbiome and large dataset from the pelagic zone of oceans. Tiara is also the only available tool correctly classifying organellar sequences, which was confirmed by the recovery of nearly complete plastid and mitochondrial genomes from the test data and real metagenomic data.

AVAILABILITY AND IMPLEMENTATION

Tiara is implemented in python 3.8, available at https://github.com/ibe-uw/tiara and tested on Unix-based systems. It is released under an open-source MIT license and documentation is available at https://ibe-uw.github.io/tiara. Version 1.0.1 of Tiara has been used for all benchmarks.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-Resolved Metagenomic Analyses Reveal the Presence of a Putative Bacterial Endosymbiont in an Avian Nasal Mite (Rhinonyssidae; Mesostigmata)

Rhinonyssidae (Mesostigmata) is a family of nasal mites only found in birds. All species are hematophagous endoparasites, which may damage the nasal cavities of birds, and also could be potential reservoirs or vectors of other infections. However, the role of members of Rhinonyssidae as disease vectors in wild bird populations remains uninvestigated, with studies of the microbiomes of Rhinonyssidae being almost non-existent. In the nasal mite (Tinaminyssus melloi) from rock doves (Columba livia), a previous study found evidence of a highly abundant putatively endosymbiotic bacteria from Class Alphaproteobacteria. Here, we expanded the sample size of this species (two different hosts- ten nasal mites from two independent samples per host), incorporated contamination controls, and increased sequencing depth in shotgun sequencing and genome-resolved metagenomic analyses. Our goal was to increase the information regarding this mite species and its putative endosymbiont. We obtained a metagenome assembled genome (MAG) that was estimated to be 98.1% complete and containing only 0.9% possible contamination. Moreover, the MAG has characteristics typical of endosymbionts (namely, small genome size an AT bias). Overall, our results support the presence of a potential endosymbiont, which is the first described for avian nasal mites to date, and improve the overall understanding of the microbiota inhabiting these mites.

Patterns of Microbiome Variation Among Infrapopulations of Permanent Bloodsucking Parasites

While interspecific variation in microbiome composition can often be readily explained by factors such as host species identity, there is still limited knowledge of how microbiomes vary at scales lower than the species level (e.g., between individuals or populations). Here, we evaluated variation in microbiome composition of individual parasites among infrapopulations (i.e., populations of parasites of the same species living on a single host individual). To address this question, we used genome-resolved and shotgun metagenomic data of 17 infrapopulations (balanced design) of the permanent, bloodsucking seal louse Echinophthirius horridus sampled from individual Saimaa ringed seals Pusa hispida saimensis. Both genome-resolved and read-based metagenomic classification approaches consistently show that parasite infrapopulation identity is a significant factor that explains both qualitative and quantitative patterns of microbiome variation at the intraspecific level. This study contributes to the general understanding of the factors driving patterns of intraspecific variation in microbiome composition, especially of bloodsucking parasites, and has implications for understanding how well-known processes occurring at higher taxonomic levels, such as phylosymbiosis, might arise in these systems.

Bacteria responsible for antimonite oxidation in antimony-contaminated soil revealed by DNA-SIP coupled to metagenomics

Antimony (Sb), the analog of arsenic (As), is a toxic metalloid that poses risks to the environment and human health. Antimonite (Sb(III)) oxidation can decrease Sb toxicity, which contributes to the bioremediation of Sb contamination. Bacteria can oxidize Sb(III), but the current knowledge regarding Sb(III)-oxidizing bacteria (SbOB) is limited to pure culture studies, thus underestimating the diversity of SbOB. In this study, Sb(III)-oxidizing microcosms were set up using Sb-contaminated rice paddies as inocula. Sb(III) oxidation driven by microorganisms was observed in the microcosms. The increasing copies and transcription of the arsenate-oxidizing gene, aioA, in the microcosms during biotic Sb(III) oxidation indicated that microorganisms mediated Sb(III) oxidation via the aioA genes. Furthermore, a novel combination of DNA-SIP and shotgun metagenomic was applied to identify the SbOB and predict their metabolic potential. Several putative SbOB were identified, including Paracoccus, Rhizobium, Achromobacter and Hydrogenophaga. Furthermore, the metagenomic analysis indicated that all of these putative SbOB contained aioA genes, confirming their roles in Sb(III) oxidation. These results suggested the concept of proof of combining DNA-SIP and shotgun metagenomics directly. In addition, the identification of the novel putative SbOB expands the current knowledge regarding the diversity of SbOB.

Prevalence of antibiotic resistance genes and bacterial pathogens along the soil-mangrove root continuum

Plants roots are colonised by soil bacteria that are known to be the reservoir of antibiotic resistance genes (ARGs). ARGs can transfer between these microorganisms and pathogens, but to what extent these ARGs and pathogens disseminate from soil into plant is poorly understood. Here, we examined a high-resolution resistome profile along the soil-root continuum of mangrove saplings using amplicon and metagenomic sequencing. Data revealed that 91.4% of total ARGs were shared across four root-associated compartments (endosphere, episphere, rhizosphere and unplanted soil). Rather than compartment-selective dynamics of microbiota, the resistome was disseminated in a continuous fashion along the soil-root continuum. Such dissemination was independent of underlying root-associated bacterial and fungal microbiota, but might be facilitated by a multiplicity of mobile genetic elements. As the multiple-drug resistant pathogens, Vibrio vulnificus, pathogenic Escherichia coli and Klebsiella pneumoniae consistently predominated across four compartments, indicating the potential dissemination of antibiotic pathogens along the soil-root continuum. Through deciphering the profile and dynamics of the root-associated resistome and pathogens, our study identified the soil-root continuum as an interconnected sink through which certain ARGs and pathogens can flow from soil into the plant.

Impact of oxytetracycline on anaerobic wastewater treatment and mitigation using enhanced hydrolysis pretreatment

In this study, two parallel-operated up-flow anaerobic sludge bed reactors, one used to treat synthetic wastewater spiked with oxytetracycline and the other used to treat the same wastewater after enhanced hydrolysis, were used to evaluate the impact of oxytetracycline on anaerobic digestion and resistance development and the efficacy of enhanced hydrolysis pretreatment on the elimination of adverse effects. The reactors were operated under a constant organic-loading rate (10 g/L/d) with increasing oxytetracycline doses (0 mg/L to 200 mg/L) over a period of 15 months. For the reactor without pretreatment, the chemical oxygen demand removal reached up to 89.5%%at oxytetracycline doses ranging from 0 mg/L to 100 mg/L, which collapsed at higher oxytetracycline doses. Miseq sequencing showed that a diverse hydrolysis/fermentation/acetogenesis bacterial community was maintained as the oxytetracycline dose was increased from 0 mg/L to 100 mg/L, while extreme dominance of Macellibacteroides (65.70%%- 71.56%) was found to occur at higher oxytetracycline doses. The total abundance of antibiotic resistance genes increased from 1.3 × 10^-1 copies per cell to 2.6 × 10^-1 copies per cell with increasing oxytetracycline dose from 0 mg/L to 5 mg/L, remained unchanged at oxytetracycline doses ranging from 25 mg/L to 100 mg/L, and then increased to 4.8 × 10^-1 copies per cell and 1.3 copies per cell at oxytetracycline doses of 150 mg/L and 200 mg/L, respectively. Multidrug resistance developed in response to oxytetracycline treatment at 200 mg/L. Poor chemical oxygen demand removal and a marked enrichment in antibiotic resistance genes was validated using a full-scale up-flow anaerobic sludge bed system fed with an influent oxytetracycline concentration of approximately 200 mg/L. For the reactor treating wastewater pretreated with enhanced hydrolysis (85 °C for 6 h), the chemical oxygen demand removal rate and antibiotic resistance genes level over the whole oxytetracycline dose range were found to be similar to those achieved with zero oxytetracycline treatment. These results demonstrated that the control of conventional pollutants and ARGs could be achieved simultaneously in the UASB reactor by employing enhanced hydrolysis pretreatment.

Bacteria are important dimethylsulfoniopropionate producers in marine aphotic and high-pressure environments

Dimethylsulfoniopropionate (DMSP) is an important marine osmolyte. Aphotic environments are only recently being considered as potential contributors to global DMSP production. Here, our Mariana Trench study reveals a typical seawater DMSP/dimethylsulfide (DMS) profile, with highest concentrations in the euphotic zone and decreased but consistent levels below. The genetic potential for bacterial DMSP synthesis via the dsyB gene and its transcription is greater in the deep ocean, and is highest in the sediment.s DMSP catabolic potential is present throughout the trench waters, but is less prominent below 8000 m, perhaps indicating a preference to store DMSP in the deep for stress protection. Deep ocean bacterial isolates show enhanced DMSP production under increased hydrostatic pressure. Furthermore, bacterial dsyB mutants are less tolerant of deep ocean pressures than wild-type strains. Thus, we propose a physiological function for DMSP in hydrostatic pressure protection, and that bacteria are key DMSP producers in deep seawater and sediment.

Dimethylsulfoniopropionate (DMSP) is an osmolyte produced by marine microbes that plays an important role in nutrient cycling and atmospheric chemistry. Here the authors go to the Mariana Trench—the deepest point in the ocean—and find bacteria are key DMSP producers, and that DMSP has a role in protection against high pressure.

Evaluation of computational methods for human microbiome analysis using simulated data

BACKGROUND

Our understanding of the composition, function, and health implications of human microbiota has been advanced by high-throughput sequencing and the development of new genomic analyses. However, trade-offs among alternative strategies for the acquisition and analysis of sequence data remain understudied.

METHODS

We assessed eight popular taxonomic profiling pipelines; MetaPhlAn2, metaMix, PathoScope 2.0, Sigma, Kraken, ConStrains, Centrifuge and Taxator-tk, against a battery of metagenomic datasets simulated from real data. The metagenomic datasets were modeled on 426 complete or permanent draft genomes stored in the Human Oral Microbiome Database and were designed to simulate various experimental conditions, both in the design of a putative experiment; read length (75-1,000 bp reads), sequence depth (100K-10M), and in metagenomic composition; number of species present (10, 100, 426), species distribution. The sensitivity and specificity of each of the pipelines under various scenarios were measured. We also estimated the relative root mean square error and average relative error to assess the abundance estimates produced by different methods. Additional datasets were generated for five of the pipelines to simulate the presence within a metagenome of an unreferenced species, closely related to other referenced species. Additional datasets were also generated in order to measure computational time on datasets of ever-increasing sequencing depth (up to 6 × 107).

RESULTS

Testing of eight pipelines against 144 simulated metagenomic datasets initially produced 1,104 discrete results. Pipelines using a marker gene strategy; MetaPhlAn2 and ConStrains, were overall less sensitive, than other pipelines; with the notable exception of Taxator-tk. This difference in sensitivity was largely made up in terms of runtime, significantly lower than more sensitive pipelines that rely on whole-genome alignments such as PathoScope2.0. However, pipelines that used strategies to speed-up alignment between genomic references and metagenomic reads, such as kmerization, were able to combine both high sensitivity and low run time, as is the case with Kraken and Centrifuge. Absent species genomes in the database mostly led to assignment of reads to the most closely related species available in all pipelines. Our results therefore suggest that taxonomic profilers that use kmerization have largely superseded those that use gene markers, coupling low run times with high sensitivity and specificity. Taxonomic profilers using more time-consuming read reassignment, such as PathoScope 2.0, provided the most sensitive profiles under common metagenomic sequencing scenarios. All the results described and discussed in this paper can be visualized using the dedicated R Shiny application (https://github.com/microgenomics/HumanMicrobiomeAnalysis). All of our datasets, pipelines and results are made available through the GitHub repository for future benchmarking.

Metagenomic approaches in microbial ecology: an update on whole-genome and marker gene sequencing analyses

Metagenomics and marker gene approaches, coupled with high-throughput sequencing technologies, have revolutionized the field of microbial ecology. Metagenomics is a culture-independent method that allows the identification and characterization of organisms from all kinds of samples. Whole-genome shotgun sequencing analyses the total DNA of a chosen sample to determine the presence of micro-organisms from all domains of life and their genomic content. Importantly, the whole-genome shotgun sequencing approach reveals the genomic diversity present, but can also give insights into the functional potential of the micro-organisms identified. The marker gene approach is based on the sequencing of a specific gene region. It allows one to describe the microbial composition based on the taxonomic groups present in the sample. It is frequently used to analyse the biodiversity of microbial ecosystems. Despite its importance, the analysis of metagenomic sequencing and marker gene data is quite a challenge. Here we review the primary workflows and software used for both approaches and discuss the current challenges in the field.

Reservoirs of antimicrobial resistance genes in retail raw milk

BACKGROUND

It has been estimated that at least 3% of the USA population consumes unpasteurized (raw) milk from animal sources, and the demand to legalize raw milk sales continues to increase. However, consumption of raw milk can cause foodborne illness and be a source of bacteria containing transferrable antimicrobial resistance genes (ARGs). To obtain a comprehensive understanding of the microbiome and antibiotic resistome in both raw and processed milk, we systematically analyzed 2034 retail milk samples including unpasteurized milk and pasteurized milk via vat pasteurization, high-temperature-short-time pasteurization, and ultra-pasteurization from the United States using complementary culture-based, 16S rRNA gene, and metagenomic sequencing techniques.

RESULTS

Raw milk samples had the highest prevalence of viable bacteria which were measured as all aerobic bacteria, coliform, and Escherichia coli counts, and their microbiota was distinct from other types of milk. 16S rRNA gene sequencing revealed that Pseudomonadaceae dominated raw milk with limited levels of lactic acid bacteria. Among all milk samples, the microbiota remained stable with constant bacterial populations when stored at 4 °C. In contrast, storage at room temperature dramatically enriched the bacterial populations present in raw milk samples and, in parallel, significantly increased the richness and abundance of ARGs. Metagenomic sequencing indicated raw milk possessed dramatically more ARGs than pasteurized milk, and a conjugation assay documented the active transfer of blaCMY-2, one ceftazidime resistance gene present in raw milk-borne E. coli, across bacterial species. The room temperature-enriched resistome differed in raw milk from distinct geographic locations, a difference likely associated with regionally distinct milk microbiota.

CONCLUSION

Despite advertised "probiotic" effects, our results indicate that raw milk microbiota has minimal lactic acid bacteria. In addition, retail raw milk serves as a reservoir of ARGs, populations of which are readily amplified by spontaneous fermentation. There is an increased need to understand potential food safety risks from improper transportation and storage of raw milk with regard to ARGs. Video Abstract.

Characterization of Nitrate-Dependent As(III)-Oxidizing Communities in Arsenic-Contaminated Soil and Investigation of Their Metabolic Potentials by the Combination of DNA-Stable Isotope Probing and Metagenomics

Arsenite (As(III)) oxidation has important environmental implications by decreasing both the mobility and toxicity of As in the environment. Microbe-mediated nitrate-dependent As(III) oxidation (NDAO) may be an important process for As(III) oxidation in anoxic environments. Our current knowledge of nitrate-dependent As(III)-oxidizing bacteria (NDAB), however, is largely based on isolates, and thus, the diversity of NDAB may be underestimated. In this study, DNA-stable isotope probing (SIP) with ¹³C-labeled NaHCO₃ as the sole carbon source, amplicon sequencing, and shotgun metagenomics were combined to identify NDAB and investigate their NDAO metabolism. As(III) oxidation was observed in the treatment amended with nitrate, while no obvious As(III) oxidation was observed without nitrate addition. The increase in the gene copies of aioA in the nitrate-amended treatment suggested that As(III) oxidation was mediated by microorganisms containing the aioA genes. Furthermore, diverse putative NDAB were identified in the As-contaminated soil cultures, such as Azoarcus, Rhodanobacter, Pseudomonas, and Burkholderiales-related bacteria. Metagenomic analysis further indicated that most of these putative NDAB contained genes for As(III) oxidation and nitrate reduction, confirming their roles in NDAO. The identification of novel putative NDAB expands current knowledge regarding the diversity of NDAB. The current study also suggests the proof of concept of using DNA-SIP to identify the slow-growing NDAB.

A review of methods and databases for metagenomic classification and assembly

Microbiome research has grown rapidly over the past decade, with a proliferation of new methods that seek to make sense of large, complex data sets. Here, we survey two of the primary types of methods for analyzing microbiome data: read classification and metagenomic assembly, and we review some of the challenges facing these methods. All of the methods rely on public genome databases, and we also discuss the content of these databases and how their quality has a direct impact on our ability to interpret a microbiome sample.

Longitudinal survey of microbiome associated with particulate matter in a megacity

BACKGROUND

While the physical and chemical properties of airborne particulate matter (PM) have been extensively studied, their associated microbiome remains largely unexplored. Here, we performed a longitudinal metagenomic survey of 106 samples of airborne PM_2.5 and PM₁₀ in Beijing over a period of 6 months in 2012 and 2013, including those from several historically severe smog events.

RESULTS

We observed that the microbiome composition and functional potential were conserved between PM_2.5 and PM₁₀, although considerable temporal variations existed. Among the airborne microorganisms, Propionibacterium acnes, Escherichia coli, Acinetobacter lwoffii, Lactobacillus amylovorus, and Lactobacillus reuteri dominated, along with several viral species. We further identified an extensive repertoire of genes involved in antibiotic resistance and detoxification, including transporters, transpeptidases, and thioredoxins. Sample stratification based on Air Quality Index (AQI) demonstrated that many microbial species, including those associated with human, dog, and mouse feces, exhibit AQI-dependent incidence dynamics. The phylogenetic and functional diversity of air microbiome is comparable to those of soil and water environments, as its composition likely derives from a wide variety of sources.

CONCLUSIONS

Airborne particulate matter accommodates rich and dynamic microbial communities, including a range of microbial elements that are associated with potential health consequences.

Metagenomics methods for the study of plant-associated microbial communities: A review

Plant microbiota have different effects on the plant which can be beneficial or pathogenic. In this study, we concentrated on beneficial microbes associated with plants using endophytic microbes as a case study. Detailed knowledge of the microbial diversity, abundance, composition, functional genes patterns, and metabolic pathways at genome level could assist in understanding the contributions of microbial community towards plant growth and health. Recently, the study of microbial community has improved greatly with the discovery of next-generation sequencing and bioinformatics technologies. Analysis of next generation sequencing data and a proper computational method plays a key role in examining microbial metagenome. This review presents the general metagenomics and computational methods used in processing plant associated metagenomes with concentration on endophytes. This includes 1) introduction of plant-associated microbiota and the factors driving their diversity. 2) plant metagenome focusing on DNA extraction, verification and quality control. 3) metagenomics methods used in community analysis of endophytes focusing on maize plant and, 4) computational methods used in the study of endophytic microbiomes. Limitations and future prospects of metagenomics and computational methods for the analysis of plant-associated metagenome (endophytic metagenome) were also discussed with the aim of fostering its development. We conclude that there is need to adopt advanced genomic features such as k-mers of random size, which do not depend on annotation and can represent other sequence alternatives.

Current challenges and best-practice protocols for microbiome analysis

Analyzing the microbiome of diverse species and environments using next-generation sequencing techniques has significantly enhanced our understanding on metabolic, physiological and ecological roles of environmental microorganisms. However, the analysis of the microbiome is affected by experimental conditions (e.g. sequencing errors and genomic repeats) and computationally intensive and cumbersome downstream analysis (e.g. quality control, assembly, binning and statistical analyses). Moreover, the introduction of new sequencing technologies and protocols led to a flood of new methodologies, which also have an immediate effect on the results of the analyses. The aim of this work is to review the most important workflows for 16S rRNA sequencing and shotgun and long-read metagenomics, as well as to provide best-practice protocols on experimental design, sample processing, sequencing, assembly, binning, annotation and visualization. To simplify and standardize the computational analysis, we provide a set of best-practice workflows for 16S rRNA and metagenomic sequencing data (available at https://github.com/grimmlab/MicrobiomeBestPracticeReview).

Optimizing sequencing protocols for leaderboard metagenomics by combining long and short reads

As metagenomic studies move to increasing numbers of samples, communities like the human gut may benefit more from the assembly of abundant microbes in many samples, rather than the exhaustive assembly of fewer samples. We term this approach leaderboard metagenome sequencing. To explore protocol optimization for leaderboard metagenomics in real samples, we introduce a benchmark of library prep and sequencing using internal references generated by synthetic long-read technology, allowing us to evaluate high-throughput library preparation methods against gold-standard reference genomes derived from the samples themselves. We introduce a low-cost protocol for high-throughput library preparation and sequencing.

Robust taxonomic classification of uncharted microbial sequences and bins with CAT and BAT

Current-day metagenomics analyses increasingly involve de novo taxonomic classification of long DNA sequences and metagenome-assembled genomes. Here, we show that the conventional best-hit approach often leads to classifications that are too specific, especially when the sequences represent novel deep lineages. We present a classification method that integrates multiple signals to classify sequences (Contig Annotation Tool, CAT) and metagenome-assembled genomes (Bin Annotation Tool, BAT). Classifications are automatically made at low taxonomic ranks if closely related organisms are present in the reference database and at higher ranks otherwise. The result is a high classification precision even for sequences from considerably unknown organisms.

Robust taxonomic classification of uncharted microbial sequences and bins with CAT and BAT

Current-day metagenomics analyses increasingly involve de novo taxonomic classification of long DNA sequences and metagenome-assembled genomes. Here, we show that the conventional best-hit approach often leads to classifications that are too specific, especially when the sequences represent novel deep lineages. We present a classification method that integrates multiple signals to classify sequences (Contig Annotation Tool, CAT) and metagenome-assembled genomes (Bin Annotation Tool, BAT). Classifications are automatically made at low taxonomic ranks if closely related organisms are present in the reference database and at higher ranks otherwise. The result is a high classification precision even for sequences from considerably unknown organisms.

The fecal resistome of dairy cattle is associated with diet during nursing

Antimicrobial resistance is a global public health concern, and livestock play a significant role in selecting for resistance and maintaining such reservoirs. Here we study the succession of dairy cattle resistome during early life using metagenomic sequencing, as well as the relationship between resistome, gut microbiota, and diet. In our dataset, the gut of dairy calves serves as a reservoir of 329 antimicrobial resistance genes (ARGs) presumably conferring resistance to 17 classes of antibiotics, and the abundance of ARGs declines gradually during nursing. ARGs appear to co-occur with antibacterial biocide or metal resistance genes. Colostrum is a potential source of ARGs observed in calves at day 2. The dynamic changes in the resistome are likely a result of gut microbiota assembly, which is closely associated with diet transition in dairy calves. Modifications in the resistome may be possible via early-life dietary interventions to reduce overall antimicrobial resistance.