imGLAD: accurate detection and quantification of target organisms in metagenomes

Metagenomic insights into the impact of litter from poultry Concentrated Animal Feeding Operations (CAFOs) to adjacent soil and water microbial communities

In recent decades, human food consumption has led to an increased demand for animal-based foods, particularly chicken meat production. The state of Georgia, USA is one of the top broiler chicken producers in the United States, where animals are raised in Concentrated Animal Feeding Operations (CAFOs). Without proper management, CAFOs could negatively impact the environment and become a public health risk as a source of water and air pollution and/or by spreading antimicrobial resistance genes. In this study, we used metagenome sequencing to investigate the impact of the application of the CAFO's litter on adjacent soils and downstream creek waters in terms of microbial diversity and antimicrobial resistance profile changes. Our data indicate that while a few microbial groups increased in abundance within a short period of time after litter application, these populations subsequently decreased to levels similar to those found prior to the litter application or to below the detection limit of our metagenome sequencing effort. Microbial taxonomic composition analyses, relative abundance of Metagenome-Assembled Genomes (MAGs) and detection of Antimicrobial Resistance Genes (ARGs) allow us to conclude that this practice of litter application had a negligible effect on the microbiome or resistome profile of these soils and nearby waterways, likely due to its dilution in the field and/or outcompetition by indigenous microbes, revealing a minimal impact of these poultry facilities on the natural microbial communities.

A user's guide to the bioinformatic analysis of shotgun metagenomic sequence data for bacterial pathogen detection

Metagenomics, i.e., shotgun sequencing of the total microbial community DNA from a sample, has become a mature technique but its application to pathogen detection in clinical, environmental, and food samples is far from common or standardized. In this review, we summarize ongoing developments in metagenomic sequence analysis that facilitate its wider application to pathogen detection. We examine theoretical frameworks for estimating the limit of detection for a particular level of sequencing effort, current approaches for achieving species and strain analytical resolution, and discuss some relevant modern tools for these tasks. While these recent advances are significant and establish metagenomics as a powerful tool to provide insights not easily attained by culture-based approaches, metagenomics is unlikely to emerge as a widespread, routine monitoring tool in the near future due to its inherently high detection limits, cost, and inability to easily distinguish between viable and non-viable cells. Instead, metagenomics seems best poised for applications involving special circumstances otherwise challenging for culture-based and molecular (e.g., PCR-based) approaches such as the de novo detection of novel pathogens, cases of co-infection by more than one pathogen, and situations where it is important to assess the genomic composition of the pathogenic population(s) and/or its impact on the indigenous microbiome.

Fecal microbiota transplantation promotes reduction of antimicrobial resistance by strain replacement

Multidrug-resistant organism (MDRO) colonization is a fundamental challenge in antimicrobial resistance. Limited studies have shown that fecal microbiota transplantation (FMT) can reduce MDRO colonization, but its mechanisms are poorly understood. We conducted a randomized, controlled trial of FMT for MDRO decolonization in renal transplant recipients called PREMIX (NCT02922816). Eleven participants were enrolled and randomized 1:1 to FMT or an observation period followed by delayed FMT if stool cultures were MDRO positive at day 36. Participants who were MDRO positive after one FMT were treated with a second FMT. At last visit, eight of nine patients who completed all treatments were MDRO culture negative. FMT-treated participants had longer time to recurrent MDRO infection versus PREMIX-eligible controls who were not treated with FMT. Key taxa (Akkermansia muciniphila, Alistipes putredinis, Phocaeicola dorei, Phascolarctobacterium faecium, Alistipes species, Mesosutterella massiliensis, Barnesiella intestinihominis, and Faecalibacterium prausnitzii) from the single feces donor used in the study that engrafted in recipients and metabolites such as short-chain fatty acids and bile acids in FMT-responding participants uncovered leads for rational microbiome therapeutic and diagnostic development. Metagenomic analyses revealed a previously unobserved mechanism of MDRO eradication by conspecific strain competition in an FMT-treated subset. Susceptible Enterobacterales strains that replaced baseline extended-spectrum β-lactamase-producing strains were not detectable in donor microbiota manufactured as FMT doses but in one case were detectable in the recipient before FMT. These data suggest that FMT may provide a path to exploit strain competition to reduce MDRO colonization.

A genome catalogue of lake bacterial diversity and its drivers at continental scale

Lakes are heterogeneous ecosystems inhabited by a rich microbiome whose genomic diversity is poorly defined. We present a continental-scale study of metagenomes representing 6.5 million km2 of the most lake-rich landscape on Earth. Analysis of 308 Canadian lakes resulted in a metagenome-assembled genome (MAG) catalogue of 1,008 mostly novel bacterial genomospecies. Lake trophic state was a leading driver of taxonomic and functional diversity among MAG assemblages, reflecting the responses of communities profiled by 16S rRNA amplicons and gene-centric metagenomics. Coupling the MAG catalogue with watershed geomatics revealed terrestrial influences of soils and land use on assemblages. Agriculture and human population density were drivers of turnover, indicating detectable anthropogenic imprints on lake bacteria at the continental scale. The sensitivity of bacterial assemblages to human impact reinforces lakes as sentinels of environmental change. Overall, the LakePulse MAG catalogue greatly expands the freshwater genomic landscape, advancing an integrative view of diversity across Earth's microbiomes.

Pathogen Detection Using Metagenomic Next-Generation Sequencing of Plasma Samples from Patients with Sepsis in Uganda

Metagenomic sequencing is a promising new method for pathogen detection. We aimed to detect pathogens from archived plasma using metagenomic sequencing in a previously well-characterized cohort of 254 predominantly HIV-infected patients with sepsis in Uganda. We used Illumina sequencing and the Chan Zuckerberg ID metagenomics platform to sequence and identify pathogens. On average, each plasma sample yielded 3,404,737 ± 2,201,997 reads (mean ± standard deviation), of which 220,032 ± 416,691 (6.3% ± 8.6%) were identified as nonhuman reads. Using a background model filter, 414 genus-specific pathogen identifications were found in the 254 samples. Nineteen pathogens were previously detected positive by quantitative PCR (qPCR), compared to sequencing, which demonstrated 30.2% sensitivity and 99.5% specificity. Sensitivity was higher for viral pathogens than nonviral pathogens (37% versus 5%). For example, HIV viremia was detected in 69% of samples using qPCR, and sequencing revealed 70% sensitivity and 92% specificity. There were 75 genus-specific potential pathogens identified by sequencing in this cohort, including hepatitis B and Epstein-Barr virus (EBV), among several others. qPCR showed a prevalence of hepatitis B and EBV viremia of 17% and 45%, respectively. In-hospital mortality was associated with a lower qPCR threshold cycle value for EBV (adjusted odds ratio, 0.85; P < .001) but not for hepatitis B or HIV. In conclusion, a broad range of potential pathogens were identified by metagenomic sequencing in patients with sepsis in Uganda. Unexpectedly high rates of hepatitis B and EBV viremia were found. Whether these viral infections in HIV patients with sepsis are clinically important requires further study. IMPORTANCE The use of next-generation sequencing (NGS) in blood samples is an emerging technology for clinical microbiology labs. In this work, we performed NGS on plasma samples from a well-characterized cohort, where all samples had been previously tested by PCR for 43 pathogens. Therefore, we could compare sequencing performance against that of PCR and identify clinical correlates. A broad range of potential pathogens were identified by metagenomic sequencing in patients with sepsis in Uganda, particularly viruses, which we confirmed by PCR. In addition to HIV viremia, unexpectedly high rates of hepatitis B and EBV viremia were found, which may have important clinical implications.

Metagenomics indicate that public health risk may be higher from flooding following dry versus rainy periods

Urban floodwater could lead to significant risk for public and environmental health from mobilization of microbial pathogens and overflow of wastewater treatment systems. Here, we attempted to assess this risk by obtaining metagenomic profiles of antibiotic resistance genes (ARGs), virulence factors (VFs) and pathogens present in floodwater samples collected in urban Atlanta, GA that were categorized in two distinct groups: floods that occurred after periods of drought and those after regular (seasonal) rain events. Even though no major (known) pathogens were present at the limit of detection of our sequencing effort (~3 Gbp/sample), we observed that floodwaters after drought showed a 2.5-fold higher abundance of both ARGs and VFs compared to floodwater after rainy days. These differences were mainly derived by several novel species of the Pseudomonas genus, which were more dominant in the former versus the latter samples and carried several genes to cope with osmotic stress in addition to ARGs and VFs. These results revealed that there are previously undescribed species that become mobilized after flooding events in the Southeast US urban settings and could represent an increased public health risk, especially after periods of drought, which warrants further attention.

Toward shotgun metagenomic approaches for microbial source tracking sewage spills based on laboratory mesocosms

Little is known about the genomic diversity of the microbial communities associated with raw municipal wastewater (sewage), including whether microbial populations specific to sewage exist and how such populations could be used to improve source attribution and apportioning in contaminated waters. Herein, we used the influent of three wastewater treatment plants in Atlanta, Georgia (USA) to perturb laboratory freshwater mesocosms, simulating sewage contamination events, and followed these mesocosms with shotgun metagenomics over a 7-day observational period. We describe 15 abundant non-redundant bacterial metagenome-assembled genomes (MAGs) ubiquitous within all sewage inocula yet absent from the unperturbed freshwater control at our analytical limit of detection. Tracking the dynamics of the populations represented by these MAGs revealed varied decay kinetics, depending on (inferred) phenotypes, e.g., anaerobes decayed faster than aerobes under the well-aerated incubation conditions. Notably, a portion of these populations showed decay patterns similar to those of common markers, Enterococcus and HF183. Despite the apparent decay of these populations, the abundance of β-lactamase encoding genes remained high throughout incubation relative to the control. Lastly, we constructed genomic libraries representing several different fecal sources and outline a bioinformatic approach which leverages these libraries for identifying and apportioning contamination signal among multiple probable sources using shotgun metagenomic data.

The search for Schizosaccharomyces fission yeasts in environmental meta-transcriptomes

Fission yeast is an important model organism in evolutionary genetics and cell biology research. Nevertheless, most research is limited to a single laboratory strain and knowledge of its natural occurrence is limited, which reduces our understanding of its life history and hinders isolation of new strains from nature. Understanding the natural diversity of fission yeast can provide insight into its genetic and phenotypic diversity and the evolutionary processes that shaped these. Here we aimed to identify candidate natural habitats of fission yeasts by searching through a large collection of publicly available environmental metatranscriptomic datasets. Using a custom pipeline, we processed over 13,000 NCBI SRA accessions, from a wide range of 34 different environmental categories. Overall, we found a very low abundance of putative yeast transcripts, with most fission yeast signatures coming from the categories of 'food' and 'terrestrial arthropods'. Additionally, a signal could be found in a variety of marine and fresh aquatic habitats. Our results do not provide a conclusive answer on the natural habitat of fission yeasts, but our analysis further narrows the range of locations where fission yeasts naturally occur.

Longitudinal study on the effects of growth-promoting and therapeutic antibiotics on the dynamics of chicken cloacal and litter microbiomes and resistomes

BACKGROUND

Therapeutic and growth-promoting antibiotics are frequently used in broiler production. Indirect evidence indicates that these practices are linked to the proliferation of antimicrobial resistance (AMR), the spread of antibiotic-resistant bacteria from food animals to humans, and the environment, but there is a lack of comprehensive experimental data supporting this. We investigated the effects of growth promotor (bacitracin) and therapeutic (enrofloxacin) antibiotic administration on AMR in broilers for the duration of a production cycle, using a holistic approach that integrated both culture-dependent and culture-independent methods. We specifically focused on pathogen-harboring families (Enterobacteriaceae, Enterococcaceae, and Staphylococcaceae).

RESULTS

Antibiotic-resistant bacteria and antibiotic resistance genes were ubiquitous in chicken cloaca and litter regardless of antibiotic administration. Environment (cloaca vs. litter) and growth stage were the primary drivers of variation in the microbiomes and resistomes, with increased bacterial diversity and a general decrease in abundance of the pathogen-harboring families with age. Bacitracin-fed groups had higher levels of bacitracin resistance genes and of vancomycin-resistant Enterococcaceae (total Enterococcaceae counts were not higher). Although metagenomic analyses classified 28-76% of the Enterococcaceae as the commensal human pathogens E. faecalis and E. faecium, culture-based analysis suggested that approximately 98% of the vancomycin-resistant Enterococcaceae were avian and not human-associated, suggesting differences in the taxonomic profiles of the resistant and non-resistant strains. Enrofloxacin treatments had varying effects, but generally facilitated increased relative abundance of multidrug-resistant Enterobacteriaceae strains, which were primarily E. coli. Metagenomic approaches revealed a diverse array of Staphylococcus spp., but the opportunistic pathogen S. aureus and methicillin resistance genes were not detected in culture-based or metagenomic analyses. Camphylobacteriaceae were significantly more abundant in the cloacal samples, especially in enrofloxacin-treated chickens, where a metagenome-assembled C. jejuni genome harboring fluoroquinolone and β-lactam resistance genes was identified.

CONCLUSIONS

Within a "farm-to-fork, one health" perspective, considering the evidence that bacitracin and enrofloxacin used in poultry production can select for resistance, we recommend their use be regulated. Furthermore, we suggest routine surveillance of ESBL E. coli, vancomycin-resistant E. faecalis and E. faecium, and fluoroquinolone-resistant C. jejuni strains considering their pathogenic nature and capacity to disseminate AMR to the environment. Video Abstract.

Description of Candidatus Mesopelagibacter carboxydoxydans and Candidatus Anoxipelagibacter denitrificans: Nitrate-reducing SAR11 genera that dominate mesopelagic and anoxic marine zones

The diverse and ubiquitous members of the SAR11 lineage (Alphaproteobacteria) represent up to 30-40% of the surface and mesopelagic oceanic microbial communities. However, the molecular and ecological mechanisms that differentiate closely related, yet distinct, SAR11 members that often co-occur under similar environmental conditions remain speculative. Recently, two mesopelagic and oxygen minimum zone (OMZ)-associated subclades of SAR11 (Ic and IIa.A) were described using single-cell amplified genomes (SAGs) linked to nitrate reduction in OMZs. In this current study, the collection of genomes belonging to these two subclades was expanded with thirteen new metagenome-assembled genomes (MAGs), thus providing a more detailed phylogenetic and functional characterization of these subclades. Gene content-based predictions of metabolic functions revealed similarities in central carbon metabolism between subclades Ic and IIa.A and surface SAR11 clades, with small variations in central pathways. These variations included more versatile sulfur assimilation pathways, as well as a previously predicted capacity for nitrate reduction that conferred unique versatility on mesopelagic-adapted clades compared to their surface counterparts. Finally, consistent with previously reported abundances of carbon monoxide (CO) in surface and mesopelagic waters, subclades Ia (surface) and Ic (mesopelagic) have the genetic potential to oxidize carbon monoxide (CO), presumably taking advantage of this abundant compound as an electron donor. Based on genomic analyses, environmental distribution and metabolic reconstruction, we propose two new SAR11 genera, Ca. Mesopelagibacter carboxydoxydans (subclade Ic) and Ca. Anoxipelagibacter denitrificans (subclade IIa.A), which represent members of the mesopelagic and OMZ-adapted SAR11 clades.

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.

Genome repository of oil systems: An interactive and searchable database that expands the catalogued diversity of crude oil-associated microbes

Microbial communities ultimately control the fate of petroleum hydrocarbons (PHCs) that enter the natural environment, but the interactions of microbes with PHCs and the environment are highly complex and poorly understood. Genome-resolved metagenomics can help unravel these complex interactions. However, the lack of a comprehensive database that integrates existing genomic/metagenomic data from oil environments with physicochemical parameters known to regulate the fate of PHCs currently limits data analysis and interpretations. Here, we curated a comprehensive, searchable database that documents microbial populations in natural oil ecosystems and oil spills, along with available underlying physicochemical data, geocoded via geographic information system to reveal their geographic distribution patterns. Analysis of the ~2000 metagenome-assembled genomes (MAGs) available in the database revealed strong ecological niche specialization within habitats. Over 95% of the recovered MAGs represented novel taxa underscoring the limited representation of cultured organisms from oil-contaminated and oil reservoir ecosystems. The majority of MAGs linked to oil-contaminated ecosystems were detectable in non-oiled samples from the Gulf of Mexico but not in comparable samples from elsewhere, indicating that the Gulf is primed for oil biodegradation. The repository should facilitate future work toward a predictive understanding of the microbial taxa and their activities that control the fate of oil spills.

Metagenomics as a Public Health Risk Assessment Tool in a Study of Natural Creek Sediments Influenced by Agricultural and Livestock Runoff: Potential and Limitations

Little is known about the public health risks associated with natural creek sediments that are affected by runoff and fecal pollution from agricultural and livestock practices. For instance, the persistence of foodborne pathogens such as Shiga toxin-producing Escherichia coli (STEC) originating from these practices remains poorly quantified. Towards closing these knowledge gaps, the water-sediment interface of two creeks in the Salinas River Valley of California was sampled over a 9-month period using metagenomics and traditional culture-based tests for STEC. Our results revealed that these sediment communities are extremely diverse and have functional and taxonomic diversity comparable to that observed in soils. With our sequencing effort (∼4 Gbp per library), we were unable to detect any pathogenic E. coli in the metagenomes of 11 samples that had tested positive using culture-based methods, apparently due to relatively low abundance. Furthermore, there were no significant differences in the abundance of human- or cow-specific gut microbiome sequences in the downstream impacted sites compared to that in upstream more pristine (control) sites, indicating natural dilution of anthropogenic inputs. Notably, the high number of metagenomic reads carrying antibiotic resistance genes (ARGs) found in all samples was significantly higher than ARG reads in other available freshwater and soil metagenomes, suggesting that these communities may be natural reservoirs of ARGs. The work presented here should serve as a guide for sampling volumes, amount of sequencing to apply, and what bioinformatics analyses to perform when using metagenomics for public health risk studies of environmental samples such as sediments.IMPORTANCE Current agricultural and livestock practices contribute to fecal contamination in the environment and the spread of food- and waterborne disease and antibiotic resistance genes (ARGs). Traditionally, the level of pollution and risk to public health are assessed by culture-based tests for the intestinal bacterium Escherichia coli However, the accuracy of these traditional methods (e.g., low accuracy in quantification, and false-positive signal when PCR based) and their suitability for sediments remain unclear. We collected sediments for a time series metagenomics study from one of the most highly productive agricultural regions in the United States in order to assess how agricultural runoff affects the native microbial communities and if the presence of Shiga toxin-producing Escherichia coli (STEC) in sediment samples can be detected directly by sequencing. Our study provided important information on the potential for using metagenomics as a tool for assessment of public health risk in natural environments.

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.

Ecogenomic characterization of widespread, closely-related SAR11 clades of the freshwater genus "Candidatus Fonsibacter" and proposal of Ca. Fonsibacter lacus sp. nov

The ubiquitous alpha-proteobacteria of the order "Candidatus Pelagibacterales" (SAR11) are highly abundant in aquatic environments, and among them, members of the monophyletic lineage LD12 (also known as SAR11 clade IIIb) are specifically found in lacustrine ecosystems. Clade IIIb bacteria are some of the most prominent members of freshwater environments, but little is known about their biology due to the lack of genome representatives. Only recently, the first non-marine isolate was cultured and described as "Candidatus Fonsibacter ubiquis". Here, we expand the collection of freshwater IIIb representatives and describe a new IIIb species of the genus "Ca. Fonsibacter". Specifically, we assembled a collection of 67 freshwater metagenomic datasets from the interconnected lakes of the Chattahoochee River basin (GA, USA) and obtained nearly complete metagenome-assembled genomes (MAGs) representing 5 distinct IIIb subclades, roughly equivalent to species based on genomic standards, including the previously described "Ca. F. ubiquis". Genomic comparisons between members of the IIIb species revealed high similarity in gene content. However, when comparing their abundance profiles in the Chattahoochee basin and various aquatic environments, differences in temporal and spatial distributions among the distinct species were observed implying niche differentiation might be underlying the coexistence of the highly functionally similar representatives. The name Ca. Fonsibacter lacus sp. nov. is proposed for the most abundant and widespread species in the Chattahoochee River basin and various freshwater ecosystems.