Time-series sewage metagenomics distinguishes seasonal, human-derived and environmental microbial communities potentially allowing source-attributed surveillance.

Sewage metagenomics has risen to prominence in urban population surveillance of pathogens and antimicrobial resistance (AMR). Unknown species with similarity to known genomes cause database bias in reference-based metagenomics. To improve surveillance, we seek to recover sewage genomes and develop a quantification and correlation workflow for these genomes and AMR over time. We use longitudinal sewage sampling in seven treatment plants from five major European cities to explore the utility of catch-all sequencing of these population-level samples. Using metagenomic assembly methods, we recover 2332 metagenome-assembled genomes (MAGs) from prokaryotic species, 1334 of which were previously undescribed. These genomes account for ~69% of sequenced DNA and provide insight into sewage microbial dynamics. Rotterdam (Netherlands) and Copenhagen (Denmark) show strong seasonal microbial community shifts, while Bologna, Rome, (Italy) and Budapest (Hungary) have occasional blooms of Pseudomonas-dominated communities, accounting for up to ~95% of sample DNA. Seasonal shifts and blooms present challenges for effective sewage surveillance. We find that bacteria of known shared origin, like human gut microbiota, form communities, suggesting the potential for source-attributing novel species and their ARGs through network community analysis. This could significantly improve AMR tracking in urban environments.

Discovery of Vibrio cholerae in Urban Sewage in Copenhagen, Denmark.

We report the discovery of a persistent presence of Vibrio cholerae at very low abundance in the inlet of a single wastewater treatment plant in Copenhagen, Denmark at least since 2015. Remarkably, no environmental or locally transmitted clinical case of V. cholerae has been reported in Denmark for more than 100 years. We, however, have recovered a near-complete genome out of 115 metagenomic sewage samples taken over the past 8 years, despite the extremely low relative abundance of one V. cholerae read out of 500,000 sequenced reads. Due to the very low relative abundance, routine screening of the individual samples did not reveal V. cholerae. The recovered genome lacks the gene responsible for cholerae toxin production, but although this strain may not pose an immediate public health risk, our finding illustrates the importance, challenges, and effectiveness of wastewater-based pathogen surveillance.

Utilizing co-abundances of antimicrobial resistance genes to identify potential co-selection in the resistome.

The rapid spread of antimicrobial resistance (AMR) is a threat to global health, and the nature of co-occurring antimicrobial resistance genes (ARGs) may cause collateral AMR effects once antimicrobial agents are used. Therefore, it is essential to identify which pairs of ARGs co-occur. Given the wealth of next-generation sequencing data available in public repositories, we have investigated the correlation between ARG abundances in a collection of 214,095 metagenomic data sets. Using more than 6.76∙108 read fragments aligned to acquired ARGs to infer pairwise correlation coefficients, we found that more ARGs correlated with each other in human and animal sampling origins than in soil and water environments. Furthermore, we argued that the correlations could serve as risk profiles of resistance co-occurring to critically important antimicrobials (CIAs). Using these profiles, we found evidence of several ARGs conferring resistance for CIAs being co-abundant, such as tetracycline ARGs correlating with most other forms of resistance. In conclusion, this study highlights the important ARG players indirectly involved in shaping the resistomes of various environments that can serve as monitoring targets in AMR surveillance programs. IMPORTANCE: Understanding the collateral effects happening in a resistome can reveal previously unknown links between antimicrobial resistance genes (ARGs). Through the analysis of pairwise ARG abundances in 214K metagenomic samples, we observed that the co-abundance is highly dependent on the environmental context and argue that these correlations can be used to show the risk of co-selection occurring in different settings.

The European livestock resistome.

Metagenomic sequencing has proven to be a powerful tool in the monitoring of antimicrobial resistance (AMR). Here, we provide a comparative analysis of the resistome from pigs, poultry, veal calves, turkey, and rainbow trout, for a total of 538 herds across nine European countries. We calculated the effects of per-farm management practices and antimicrobial usage (AMU) on the resistome in pigs, broilers, and veal calves. We also provide an in-depth study of the associations between bacterial diversity, resistome diversity, and AMR abundances as well as co-occurrence analysis of bacterial taxa and antimicrobial resistance genes (ARGs) and the universality of the latter. The resistomes of veal calves and pigs clustered together, as did those of avian origin, while the rainbow trout resistome was different. Moreover, we identified clear core resistomes for each specific food-producing animal species. We identified positive associations between bacterial alpha diversity and both resistome alpha diversity and abundance. Network analyses revealed very few taxa-ARG associations in pigs but a large number for the avian species. Using updated reference databases and optimized bioinformatics, previously reported significant associations between AMU, biosecurity, and AMR in pig and poultry farms were validated. AMU is an important driver for AMR; however, our integrated analyses suggest that factors contributing to increased bacterial diversity might also be associated with higher AMR load. We also found that dispersal limitations of ARGs are shaping livestock resistomes, and future efforts to fight AMR should continue to emphasize biosecurity measures.IMPORTANCEUnderstanding the occurrence, diversity, and drivers for antimicrobial resistance (AMR) is important to focus future control efforts. So far, almost all attempts to limit AMR in livestock have addressed antimicrobial consumption. We here performed an integrated analysis of the resistomes of five important farmed animal populations across Europe finding that the resistome and AMR levels are also shaped by factors related to bacterial diversity, as well as dispersal limitations. Thus, future studies and interventions aimed at reducing AMR should not only address antimicrobial usage but also consider other epidemiological and ecological factors.

Transmission of antimicrobial resistance in the gut microbiome of gregarious cockroaches: the importance of interaction between antibiotic exposed and non-exposed populations.

Antimicrobial resistance (AMR) is a major global health concern, further complicated by its spread via the microbiome bacterial members. While mathematical models discuss AMR transmission through the symbiotic microbiome, experimental studies are scarce. Herein, we used a gregarious cockroach, Pycnoscelus surinamensis, as an in vivo animal model for AMR transmission investigations. We explored whether the effect of antimicrobial treatment is detectable with metagenomic sequencing, and whether AMR genes can be spread and established in unchallenged (not treated with antibiotics) individuals following contact with treated donors, and under various frequencies of interaction. Gut and soil substrate microbiomes were investigated by metagenomic sequencing for bacterial community composition and resistome profiling. We found that tetracycline treatment altered the treated gut microbiome by decreasing bacterial diversity and increasing the abundance of tetracycline resistance genes. Untreated cockroaches that interacted with treated donors also had elevated tetracycline resistance. The levels of resistance differed depending on the magnitude and frequency of donor transfer. Additionally, treated donors showed signs of microbiome recovery due to their interaction with the untreated ones. Similar patterns were also recorded in the soil substrate microbiomes. Our results shed light on how interacting microbiomes facilitate AMR gene transmission to previously unchallenged hosts, a dynamic influenced by the interaction frequencies, using an in vivo model to validate theoretical AMR transmission models.IMPORTANCEAntimicrobial resistance is a rising threat to human and animal health. The spread of resistance through the transmission of the symbiotic gut microbiome is of concern and has been explored in theoretical modeling studies. In this study, we employ gregarious insect populations to examine the emergence and transmission of antimicrobial resistance in vivo and validate modeling hypotheses. We find that antimicrobial treatment increases the levels of resistance in treated populations. Most importantly, we show that resistance increased in untreated populations after interacting with the treated ones. The level of resistance transmission was affected by the magnitude and frequency of population mixing. Our results highlight the importance of microbial transmission in the spread of antimicrobial resistance.

Genomic analysis of sewage from 101 countries reveals global landscape of antimicrobial resistance.

Antimicrobial resistance (AMR) is a major threat to global health. Understanding the emergence, evolution, and transmission of individual antibiotic resistance genes (ARGs) is essential to develop sustainable strategies combatting this threat. Here, we use metagenomic sequencing to analyse ARGs in 757 sewage samples from 243 cities in 101 countries, collected from 2016 to 2019. We find regional patterns in resistomes, and these differ between subsets corresponding to drug classes and are partly driven by taxonomic variation. The genetic environments of 49 common ARGs are highly diverse, with most common ARGs carried by multiple distinct genomic contexts globally and sometimes on plasmids. Analysis of flanking sequence revealed ARG-specific patterns of dispersal limitation and global transmission. Our data furthermore suggest certain geographies are more prone to transmission events and should receive additional attention.

A curated data resource of 214K metagenomes for characterization of the global antimicrobial resistome.

The growing threat of antimicrobial resistance (AMR) calls for new epidemiological surveillance methods, as well as a deeper understanding of how antimicrobial resistance genes (ARGs) have been transmitted around the world. The large pool of sequencing data available in public repositories provides an excellent resource for monitoring the temporal and spatial dissemination of AMR in different ecological settings. However, only a limited number of research groups globally have the computational resources to analyze such data. We retrieved 442 Tbp of sequencing reads from 214,095 metagenomic samples from the European Nucleotide Archive (ENA) and aligned them using a uniform approach against ARGs and 16S/18S rRNA genes. Here, we present the results of this extensive computational analysis and share the counts of reads aligned. Over 6.76∙108 read fragments were assigned to ARGs and 3.21∙109 to rRNA genes, where we observed distinct differences in both the abundance of ARGs and the link between microbiome and resistome compositions across various sampling types. This collection is another step towards establishing global surveillance of AMR and can serve as a resource for further research into the environmental spread and dynamic changes of ARGs.

Global Distribution of mcr Gene Variants in 214K Metagenomic Samples.

Since the initial discovery of a mobilized colistin resistance gene (mcr-1), several other variants have been reported, some of which might have circulated a while beforehand. Publicly available metagenomic data provide an opportunity to reanalyze samples to understand the evolutionary history of recently discovered antimicrobial resistance genes (ARGs). Here, we present a large-scale metagenomic study of 442 Tbp of sequencing reads from 214,095 samples to describe the dissemination and emergence of nine mcr gene variants (mcr-1 to mcr-9). Our results show that the dissemination of each variant is not uniform. Instead, the source and location play a role in the spread. However, the genomic context and the genes themselves remain primarily unchanged. We report evidence of new subvariants occurring in specific environments, such as a highly prevalent and new variant of mcr-9. This work emphasizes the importance of sharing genomic data for the surveillance of ARGs in our understanding of antimicrobial resistance. IMPORTANCE The ever-growing collection of metagenomic samples available in public data repositories has the potential to reveal new details on the emergence and dissemination of mobilized colistin resistance genes. Our analysis of metagenomes deposited online in the last 10 years shows that the environmental distribution of mcr gene variants depends on sampling source and location, possibly leading to the emergence of new variants, although the contig on which the mcr genes were found remained consistent.

Long-Term Temporal Stability of the Resistome in Sewage from Copenhagen.

Antimicrobial resistance (AMR) is a major threat to global health, and it is crucial to understand the epidemiological aspects in order to predict the emergence and propagation of AMR genes. The aim of this study was to assess the variability and medium-term AMR trends within the mostly healthy human population of a single city. We monitored over 36 months (November 2015 to November 2018) the AMR level in the city of Copenhagen, Denmark, by taking bi-weekly sewage samples from the inlets of the three main water treatment plants, extracting the DNA, performing metagenomic sequencing, and read-mapping against a database of known AMR genes. We found that the AMR level was surprisingly stable with no periodic variability and no signs of drift over the measured period. We found, however, that the seemingly random variations at each site correlate in time with each other, suggesting that the variations we see are due to real environmental changes in the occurrence of AMR.IMPORTANCE The Copenhagen sewage resistome is surprisingly stable in time. The implication is that, at least for cities that are comparable to Copenhagen in terms of sewer infrastructure, few or even single samples provide a robust picture of the resistome within a city.

The gut microbiome but not the resistome is associated with urogenital schistosomiasis in preschool-aged children.

Helminth parasites have been shown to have systemic effects in the host. Using shotgun metagenomic sequencing, we characterise the gut microbiome and resistome of 113 Zimbabwean preschool-aged children (1-5 years). We test the hypothesis that infection with the human helminth parasite, Schistosoma haematobium, is associated with changes in gut microbial and antimicrobial resistance gene abundance/diversity. Here, we show that bacteria phyla Bacteroidetes, Firmicutes, Proteobacteria, and fungi phyla Ascomycota, Microsporidia, Zoopagomycota dominate the microbiome. The abundance of Proteobacteria, Ascomycota, and Basidiomycota differ between schistosome-infected versus uninfected children. Specifically, infection is associated with increases in Pseudomonas, Stenotrophomonas, Derxia, Thalassospira, Aspergillus, Tricholoma, and Periglandula, with a decrease in Azospirillum. We find 262 AMR genes, from 12 functional drug classes, but no association with individual-specific data. To our knowledge, we describe a novel metagenomic dataset of Zimbabwean preschool-aged children, indicating an association between urogenital schistosome infection and changes in the gut microbiome.

Schooling in adolescence raises IQ scores.

Although some scholars maintain that education has little effect on intelligence quotient (IQ) scores, others claim that IQ scores are indeed malleable, primarily through intervention in early childhood. The causal effect of education on IQ at later ages is often difficult to uncover because analyses based on observational data are plagued by problems of reverse causation and self-selection into further education. We exploit a reform that increased compulsory schooling from 7 to 9 y in Norway in the 1960s to estimate the effect of education on IQ. We find that this schooling reform, which primarily affected education in the middle teenage years, had a substantial effect on IQ scores measured at the age of 19 y.

Detection of the water reservoir in a forming planetary system.

Icy bodies may have delivered the oceans to the early Earth, yet little is known about water in the ice-dominated regions of extrasolar planet-forming disks. The Heterodyne Instrument for the Far-Infrared on board the Herschel Space Observatory has detected emission lines from both spin isomers of cold water vapor from the disk around the young star TW Hydrae. This water vapor likely originates from ice-coated solids near the disk surface, hinting at a water ice reservoir equivalent to several thousand Earth oceans in mass. The water's ortho-to-para ratio falls well below that of solar system comets, suggesting that comets contain heterogeneous ice mixtures collected across the entire solar nebula during the early stages of planetary birth.

Rapid advance of spring arrival dates in long-distance migratory birds.

Several bird species have advanced the timing of their spring migration in response to recent climate change. European short-distance migrants, wintering in temperate areas, have been assumed to be more affected by change in the European climate than long-distance migrants wintering in the tropics. However, we show that long-distance migrants have advanced their spring arrival in Scandinavia more than short-distance migrants. By analyzing a long-term data set from southern Italy, we show that long-distance migrants also pass through the Mediterranean region earlier. We argue that this may reflect a climate-driven evolutionary change in the timing of spring migration.