Plasmids do not consistently stabilize cooperation across bacteria but may promote broad pathogen host-range

Bacterial lifestyle shapes pangenomes

Pangenomes vary across bacteria. Some species have fluid pangenomes, with a high proportion of genes varying between individual genomes. Other species have less fluid pangenomes, with different genomes tending to contain the same genes. Two main hypotheses have been suggested to explain this variation: differences in species' bacterial lifestyle and effective population size. However, previous studies have not been able to test between these hypotheses because the different features of lifestyle and effective population size are highly correlated with each other, and phylogenetically conserved, making it hard to disentangle their relative importance. We used phylogeny-based analyses, across 126 bacterial species, to tease apart the causal role of different factors. We found that pangenome fluidity was lower in i) host-associated compared with free-living species and ii) host-associated species that are obligately dependent on a host, live inside cells, and are more pathogenic and less motile. In contrast, we found no support for the competing hypothesis that larger effective population sizes lead to more fluid pangenomes. Effective population size appears to correlate with pangenome variation because it is also driven by bacterial lifestyle, rather than because of a causal relationship.

Co-localization of antibiotic resistance genes is widespread in the infant gut microbiome and associates with an immature gut microbial composition

BACKGROUND

In environmental bacteria, the selective advantage of antibiotic resistance genes (ARGs) can be increased through co-localization with genes such as other ARGs, biocide resistance genes, metal resistance genes, and virulence genes (VGs). The gut microbiome of infants has been shown to contain numerous ARGs, however, co-localization related to ARGs is unknown during early life despite frequent exposures to biocides and metals from an early age.

RESULTS

We conducted a comprehensive analysis of genetic co-localization of resistance genes in a cohort of 662 Danish children and examined the association between such co-localization and environmental factors as well as gut microbial maturation. Our study showed that co-localization of ARGs with other resistance and virulence genes is common in the early gut microbiome and is associated with gut bacteria that are indicative of low maturity. Statistical models showed that co-localization occurred mainly in the phylum Proteobacteria independent of high ARG content and contig length. We evaluated the stochasticity of co-localization occurrence using enrichment scores. The most common forms of co-localization involved tetracycline and fluoroquinolone resistance genes, and, on plasmids, co-localization predominantly occurred in the form of class 1 integrons. Antibiotic use caused a short-term increase in mobile ARGs, while non-mobile ARGs showed no significant change. Finally, we found that a high abundance of VGs was associated with low gut microbial maturity and that VGs showed even higher potential for mobility than ARGs.

CONCLUSIONS

We found that the phenomenon of co-localization between ARGs and other resistance and VGs was prevalent in the gut at the beginning of life. It reveals the diversity that sustains antibiotic resistance and therefore indirectly emphasizes the need to apply caution in the use of antimicrobial agents in clinical practice, animal husbandry, and daily life to mitigate the escalation of resistance. Video Abstract.

A pervasive large conjugative plasmid mediates multispecies biofilm formation in the intestinal microbiota increasing resilience to perturbations

Although horizontal gene transfer is pervasive in the intestinal microbiota, we understand only superficially the roles of most exchanged genes and how the mobile repertoire affects community dynamics. Similarly, little is known about the mechanisms underlying the ability of a community to recover after a perturbation. Here, we identified and functionally characterized a large conjugative plasmid that is one of the most frequently transferred elements among Bacteroidales species and is ubiquitous in diverse human populations. This plasmid encodes both an extracellular polysaccharide and fimbriae, which promote the formation of multispecies biofilms in the mammalian gut. We use a hybridization-based approach to visualize biofilms in clarified whole colon tissue with unprecedented 3D spatial resolution. These biofilms increase bacterial survival to common stressors encountered in the gut, increasing strain resiliency, and providing a rationale for the plasmid's recent spread and high worldwide prevalence.

Microbial Primer: Cooperation in bacteria

The growth and success of many bacteria appear to rely on a stunning range of cooperative behaviours. But what is cooperation and how is it studied?

Genes for cooperation are not more likely to be carried by plasmids

Cooperation is prevalent across bacteria, but risks being exploited by non-cooperative cheats. Horizontal gene transfer, particularly via plasmids, has been suggested as a mechanism to stabilize cooperation. A key prediction of this hypothesis is that genes which are more likely to be transferred, such as those on plasmids, should be more likely to code for cooperative traits. Testing this prediction requires identifying all genes for cooperation in bacterial genomes. However, previous studies used a method which likely misses some of these genes for cooperation. To solve this, we used a new genomics tool, SOCfinder, which uses three distinct modules to identify all kinds of genes for cooperation. We compared where these genes were located across 4648 genomes from 146 bacterial species. In contrast to the prediction of the hypothesis, we found no evidence that plasmid genes are more likely to code for cooperative traits. Instead, we found the opposite-that genes for cooperation were more likely to be carried on chromosomes. Overall, the vast majority of genes for cooperation are not located on plasmids, suggesting that the more general mechanism of kin selection is sufficient to explain the prevalence of cooperation across bacteria.

SOCfinder: a genomic tool for identifying social genes in bacteria

Bacteria cooperate by working collaboratively to defend their colonies, share nutrients, and resist antibiotics. Nevertheless, our understanding of these remarkable behaviours primarily comes from studying a few well-characterized species. Consequently, there is a significant gap in our understanding of microbial social traits, particularly in natural environments. To address this gap, we can use bioinformatic tools to identify genes that control cooperative or otherwise social traits. Existing tools address this challenge through two approaches. One approach is to identify genes that encode extracellular proteins, which can provide benefits to neighbouring cells. An alternative approach is to predict gene function using annotation tools. However, these tools have several limitations. Not all extracellular proteins are cooperative, and not all cooperative behaviours are controlled by extracellular proteins. Furthermore, existing functional annotation methods frequently miss known cooperative genes. We introduce SOCfinder as a new tool to find bacterial genes that control cooperative or otherwise social traits. SOCfinder combines information from several methods, considering if a gene is likely to [1] code for an extracellular protein [2], have a cooperative functional annotation, or [3] be part of the biosynthesis of a cooperative secondary metabolite. We use data on two extensively-studied species (P. aeruginosa and B. subtilis) to show that SOCfinder is better at finding known cooperative genes than existing tools. We also use theory from population genetics to identify a signature of kin selection in SOCfinder cooperative genes, which is lacking in genes identified by existing tools. SOCfinder opens up a number of exciting directions for future research, and is available to download from https://github.com/lauriebelch/SOCfinder.

Signatures of kin selection in a natural population of the bacteria Bacillus subtilis

Laboratory experiments have suggested that bacteria perform a range of cooperative behaviors, which are favored because they are directed toward relatives (kin selection). However, there is a lack of evidence for cooperation and kin selection in natural bacterial populations. Molecular population genetics offers a promising method to study natural populations because the theory predicts that kin selection will lead to relaxed selection, which will result in increased polymorphism and divergence at cooperative genes. Examining a natural population of Bacillus subtilis, we found consistent evidence that putatively cooperative traits have higher polymorphism and greater divergence than putatively private traits expressed at the same rate. In addition, we were able to eliminate alternative explanations for these patterns and found more deleterious mutations in genes controlling putatively cooperative traits. Overall, our results suggest that cooperation is favored by kin selection, with an average relatedness of r = .79 between interacting individuals.

Interspecific competition can drive plasmid loss from a focal species in a microbial community

Plasmids are key disseminators of antimicrobial resistance genes and virulence factors, and it is therefore critical to predict and reduce plasmid spread within microbial communities. The cost of plasmid carriage is a key metric that can be used to predict plasmids' ecological fate, and it is unclear whether plasmid costs are affected by growth partners in a microbial community. We carried out competition experiments and tracked plasmid maintenance using a model system consisting of a synthetic and stable five-species community and a broad host-range plasmid, engineered to carry different payloads. We report that both the cost of plasmid carriage and its long-term maintenance in a focal strain depended on the presence of competitors, and that these interactions were species specific. Addition of growth partners increased the cost of a high-payload plasmid to a focal strain, and accordingly, plasmid loss from the focal species occurred over a shorter time frame. We propose that the destabilising effect of interspecific competition on plasmid maintenance may be leveraged in clinical and natural environments to cure plasmids from focal strains.

PlasBin-flow: a flow-based MILP algorithm for plasmid contigs binning

MOTIVATION

The analysis of bacterial isolates to detect plasmids is important due to their role in the propagation of antimicrobial resistance. In short-read sequence assemblies, both plasmids and bacterial chromosomes are typically split into several contigs of various lengths, making identification of plasmids a challenging problem. In plasmid contig binning, the goal is to distinguish short-read assembly contigs based on their origin into plasmid and chromosomal contigs and subsequently sort plasmid contigs into bins, each bin corresponding to a single plasmid. Previous works on this problem consist of de novo approaches and reference-based approaches. De novo methods rely on contig features such as length, circularity, read coverage, or GC content. Reference-based approaches compare contigs to databases of known plasmids or plasmid markers from finished bacterial genomes.

RESULTS

Recent developments suggest that leveraging information contained in the assembly graph improves the accuracy of plasmid binning. We present PlasBin-flow, a hybrid method that defines contig bins as subgraphs of the assembly graph. PlasBin-flow identifies such plasmid subgraphs through a mixed integer linear programming model that relies on the concept of network flow to account for sequencing coverage, while also accounting for the presence of plasmid genes and the GC content that often distinguishes plasmids from chromosomes. We demonstrate the performance of PlasBin-flow on a real dataset of bacterial samples.

AVAILABILITY AND IMPLEMENTATION

https://github.com/cchauve/PlasBin-flow.

Cooperative antibiotic resistance facilitates horizontal gene transfer

The rise of β-lactam resistance among pathogenic bacteria, due to the horizontal transfer of plasmid-encoded β-lactamases, is a current global health crisis. Importantly, β-lactam hydrolyzation by β-lactamases, not only protects the producing cells but also sensitive neighboring cells cooperatively. Yet, how such cooperative traits affect plasmid transmission and maintenance is currently poorly understood. Here we experimentally show that KPC-2 β-lactamase expression and extracellular activity were higher when encoded on plasmids compared with the chromosome, resulting in the elevated rescue of sensitive non-producers. This facilitated efficient plasmid transfer to the rescued non-producers and expanded the potential plasmid recipient pool and the probability of plasmid transfer to new genotypes. Social conversion of non-producers by conjugation was efficient yet not absolute. Non-cooperative plasmids, not encoding KPC-2, were moderately more competitive than cooperative plasmids when β-lactam antibiotics were absent. However, in the presence of a β-lactam antibiotic, strains with non-cooperative plasmids were efficiently outcompeted. Moreover, plasmid-free non-producers were more competitive than non-producers imposed with the metabolic burden of a plasmid. Our results suggest that cooperative antibiotic resistance especially promotes the fitness of replicons that transfer horizontally such as conjugative plasmids.

Is cooperation favored by horizontal gene transfer?

It has been hypothesized that horizontal gene transfer on plasmids can facilitate the evolution of cooperation, by allowing genes to jump between bacteria, and hence increase genetic relatedness at the cooperative loci. However, we show theoretically that horizontal gene transfer only appreciably increases relatedness when plasmids are rare, where there are many plasmid-free cells available to infect (many opportunities for horizontal gene transfer). In contrast, when plasmids are common, there are few opportunities for horizontal gene transfer, meaning relatedness is not appreciably increased, and so cooperation is not favored. Plasmids, therefore, evolve to be rare and cooperative, or common and noncooperative, meaning plasmid frequency and cooperativeness are never simultaneously high. The overall level of plasmid-mediated cooperation, given by the product of plasmid frequency and cooperativeness, is therefore consistently negligible or low.

Diffusion and enrichment of high-risk antibiotic resistance genes (ARGs) via the transmission chain (mulberry leave, guts and feces of silkworm, and soil) in an ecological restoration area of manganese mining, China: Role of heavy metals

This study investigated the diffusion and enrichment of antibiotic resistance genes (ARGs) and pathogens via the transmission chain (mulberry leaves - silkworm guts - silkworm feces - soil) near a manganese mine restoration area (RA) and control area (CA, away from RA). Horizontal gene transfer (HGT) of ARGs was testified by an IncP a-type broad host range plasmid RP4 harboring ARGs (tetA) and conjugative genes (e.g., korB, trbA, and trbB) as an indicator. Compared to leaves, the abundances of ARGs and pathogens in feces after silkworms ingested leaves from RA increased by 10.8% and 52.3%, respectively, whereas their abundance in feces from CA dropped by 17.1% and 97.7%, respectively. The predominant ARG types in feces involved the resistances to β-lactam, quinolone, multidrug, peptide, and rifamycin. Therein, several high-risk ARGs (e.g., qnrB, oqxA, and rpoB) carried by pathogens were more enriched in feces. However, HGT mediated by plasmid RP4 in this transmission chain was not a main factor to promote the enrichment of ARGs due to the harsh survival environment of silkworm guts for the plasmid RP4 host E. coli. Notably, Zn, Mn, and As in feces and guts promoted the enrichment of qnrB and oqxA. Worriedly, the abundance of qnrB and oqxA in soil increased by over 4-fold after feces from RA were added into soil for 30 days regardless of feces with or without E. coli RP4. Overall, ARGs and pathogens could diffuse and enrich in environment via the sericulture transmission chain developed at RA, especially some high-risk ARGs carried by pathogens. Thus, greater attentions should be paid to dispel such high-risk ARGs to support benign development of sericulture industry in the safe utilization of some RAs.

Social Diversification Driven by Mobile Genetic Elements

Social diversification in microbes is an evolutionary process where lineages bifurcate into distinct populations that cooperate with themselves but not with other groups. In bacteria, this is frequently driven by horizontal transfer of mobile genetic elements (MGEs). Here, the resulting acquisition of new genes changes the recipient's social traits and consequently how they interact with kin. These changes include discriminating behaviors mediated by newly acquired effectors. Since the producing cell is protected by cognate immunity factors, these selfish elements benefit from selective discrimination against recent ancestors, thus facilitating their proliferation and benefiting the host. Whether social diversification benefits the population at large is less obvious. The widespread use of next-generation sequencing has recently provided new insights into population dynamics in natural habitats and the roles MGEs play. MGEs belong to accessory genomes, which often constitute the majority of the pangenome of a taxon, and contain most of the kin-discriminating loci that fuel rapid social diversification. We further discuss mechanisms of diversification and its consequences to populations and conclude with a case study involving myxobacteria.

Plasmid-Encoded Traits Vary across Environments

Plasmids are key mobile genetic elements in bacterial evolution and ecology as they allow the rapid adaptation of bacteria under selective environmental changes. However, the genetic information associated with plasmids is usually considered separately from information about their environmental origin. To broadly understand what kinds of traits may become mobilized by plasmids in different environments, we analyzed the properties and accessory traits of 9,725 unique plasmid sequences from a publicly available database with known bacterial hosts and isolation sources. Although most plasmid research focuses on resistance traits, such genes made up <1% of the total genetic information carried by plasmids. Similar to traits encoded on the bacterial chromosome, plasmid accessory trait compositions (including general Clusters of Orthologous Genes [COG] functions, resistance genes, and carbon and nitrogen genes) varied across seven broadly defined environment types (human, animal, wastewater, plant, soil, marine, and freshwater). Despite their potential for horizontal gene transfer, plasmid traits strongly varied with their host's taxonomic assignment. However, the trait differences across environments of broad COG categories could not be entirely explained by plasmid host taxonomy, suggesting that environmental selection acts on the plasmid traits themselves. Finally, some plasmid traits and environments (e.g., resistance genes in human-related environments) were more often associated with mobilizable plasmids (those having at least one detected relaxase) than others. Overall, these findings underscore the high level of diversity of traits encoded by plasmids and provide a baseline to investigate the potential of plasmids to serve as reservoirs of adaptive traits for microbial communities. IMPORTANCE Plasmids are well known for their role in the transmission of antibiotic resistance-conferring genes. Beyond human and clinical settings, however, they disseminate many other types of genes, including those that contribute to microbially driven ecosystem processes. In this study, we identified the distribution of traits genetically encoded by plasmids isolated from seven broadly categorized environments. We find that plasmid trait content varied with both bacterial host taxonomy and environment and that, on average, half of the plasmids were potentially mobilizable. As anthropogenic activities impact ecosystems and the climate, investigating and identifying the mechanisms of how microbial communities can adapt will be imperative for predicting the impacts on ecosystem functioning.

Comparative Genomics of the Genus Pseudomonas Reveals Host- and Environment-Specific Evolution

Each Earth ecosystem has unique microbial communities. Pseudomonas bacteria have evolved to occupy a plethora of different ecological niches, including living hosts, such as animals and plants. Many genes necessary for the Pseudomonas-niche interaction and their encoded functions remain unknown. Here, we describe a comparative genomic study of 3,274 genomes with 19,056,667 protein-coding sequences from Pseudomonas strains isolated from diverse environments. We detected functional divergence of Pseudomonas that depends on the niche. Each group of strains from a certain environment harbored a distinctive set of metabolic pathways or functions. The horizontal transfer of genes, which mainly proceeded between closely related taxa, was dependent on the isolation source. Finally, we detected thousands of undescribed proteins and functions associated with each Pseudomonas lifestyle. This research represents an effort to reveal the mechanisms underlying the ecology, pathogenicity, and evolution of Pseudomonas, and it will enable clinical, ecological, and biotechnological advances. IMPORTANCE Microbes play important roles in the health of living beings and in the environment. The knowledge of these functions may be useful for the development of new clinical and biotechnological applications and the restoration and preservation of natural ecosystems. However, most mechanisms implicated in the interaction of microbes with the environment remain poorly understood; thus, this field of research is very important. Here, we try to understand the mechanisms that facilitate the differential adaptation of Pseudomonas-a large and ubiquitous bacterial genus-to the environment. We analyzed more than 3,000 Pseudomonas genomes and searched for genetic patterns that can be related with their coevolution with different hosts (animals, plants, or fungi) and environments. Our results revealed that thousands of genes and genetic features are associated with each niche. Our data may be useful to develop new technical and theoretical advances in the fields of ecology, health, and industry.

Gene transferability and sociality do not correlate with gene connectivity

The connectivity of a gene, defined as the number of interactions a gene's product has with other genes' products, is a key characteristic of a gene. In prokaryotes, the complexity hypothesis predicts that genes which undergo more frequent horizontal transfer will be less connected than genes which are only very rarely transferred. We tested the role of horizontal gene transfer, and other potentially important factors, by examining the connectivity of chromosomal and plasmid genes, across 134 diverse prokaryotic species. We found that (i) genes on plasmids were less connected than genes on chromosomes; (ii) connectivity of plasmid genes was not correlated with plasmid mobility; and (iii) the sociality of genes (cooperative or private) was not correlated with gene connectivity.

Bacterial plasmid-associated and chromosomal proteins have fundamentally different properties in protein interaction networks

Plasmids facilitate horizontal gene transfer, which enables the diversification of pathogens into new anatomical and environmental niches, implying that plasmid-encoded genes can cooperate well with chromosomal genes. We hypothesise that such mobile genes are functionally different to chromosomal ones due to this ability to encode proteins performing non-essential functions like antimicrobial resistance and traverse distinct host cells. The effect of plasmid-driven gene gain on protein-protein interaction network topology is an important question in this area. Moreover, the extent to which these chromosomally- and plasmid-encoded proteins interact with proteins from their own groups compared to the levels with the other group remains unclear. Here, we examined the incidence and protein-protein interactions of all known plasmid-encoded proteins across representative specimens from most bacteria using all available plasmids. We found that plasmid-encoded genes constitute ~ 0.65% of the total number of genes per bacterial sample, and that plasmid genes are preferentially associated with different species but had limited taxonomical power beyond this. Surprisingly, plasmid-encoded proteins had both more protein-protein interactions compared to chromosomal proteins, countering the hypothesis that genes with higher mobility rates should have fewer protein-level interactions. Nonetheless, topological analysis and investigation of the protein-protein interaction networks' connectivity and change in the number of independent components demonstrated that the plasmid-encoded proteins had limited overall impact in > 96% of samples. This paper assembled extensive data on plasmid-encoded proteins, their interactions and associations with diverse bacterial specimens that is available for the community to investigate in more detail.

Estimating the transfer rates of bacterial plasmids with an adapted Luria–Delbrück fluctuation analysis

To increase our basic understanding of the ecology and evolution of conjugative plasmids, we need reliable estimates of their rate of transfer between bacterial cells. Current assays to measure transfer rate are based on deterministic modeling frameworks. However, some cell numbers in these assays can be very small, making estimates that rely on these numbers prone to noise. Here, we take a different approach to estimate plasmid transfer rate, which explicitly embraces this noise. Inspired by the classic fluctuation analysis of Luria and Delbrück, our method is grounded in a stochastic modeling framework. In addition to capturing the random nature of plasmid conjugation, our new methodology, the Luria–Delbrück method (“LDM”), can be used on a diverse set of bacterial systems, including cases for which current approaches are inaccurate. A notable example involves plasmid transfer between different strains or species where the rate that one type of cell donates the plasmid is not equal to the rate at which the other cell type donates. Asymmetry in these rates has the potential to bias or constrain current transfer estimates, thereby limiting our capabilities for estimating transfer in microbial communities. In contrast, the LDM overcomes obstacles of traditional methods by avoiding restrictive assumptions about growth and transfer rates for each population within the assay. Using stochastic simulations and experiments, we show that the LDM has high accuracy and precision for estimation of transfer rates compared to the most widely used methods, which can produce estimates that differ from the LDM estimate by orders of magnitude.

Plasmid transfer can often spread resistance between important clinical pathogens. This study shows that widely used methods can lead to biased estimates of plasmid transfer rate by several orders of magnitude, and presents a new approach, inspired by the classic Luria-Delbrück approach, for accurately assessing this fundamental rate parameter

Impact of horizontal gene transfer on emergence and stability of cooperative virulence in Salmonella Typhimurium

Intestinal inflammation fuels the transmission of Salmonella Typhimurium (S.Tm). However, a substantial fitness cost is associated with virulence expression. Mutations inactivating transcriptional virulence regulators generate attenuated variants profiting from inflammation without enduring virulence cost. Such variants interfere with the transmission of fully virulent clones. Horizontal transfer of functional regulatory genes (HGT) into attenuated variants could nevertheless favor virulence evolution. To address this hypothesis, we cloned hilD, coding for the master regulator of virulence, into a conjugative plasmid that is highly transferrable during intestinal colonization. The resulting mobile hilD allele allows virulence to emerge from avirulent populations, and to be restored in attenuated mutants competing against virulent clones within-host. However, mutations inactivating the mobile hilD allele quickly arise. The stability of virulence mediated by HGT is strongly limited by its cost, which depends on the hilD expression level, and by the timing of transmission. We conclude that robust evolution of costly virulence expression requires additional selective forces such as narrow population bottlenecks during transmission.

Salmonella Typhimurium virulence is costly and can be lost by mutation during infection. Bakkeren et al. show that virulence restoration via horizontal gene transfer is only transient while transmission bottlenecks promote long-term virulence stability.

Kin selection for cooperation in natural bacterial populations

Bacteria secrete many molecules outside the cell, where they provide benefits to other cells. One potential reason for producing these “public goods” is that they benefit closely related cells that share the gene for cooperation (kin selection). While many laboratory studies have supported this hypothesis, there is a lack of evidence that kin selection favors cooperation in natural populations. We examined bacterial genomes from the environment and used population genetics theory to analyze the DNA sequences. Our analyses suggest that public goods cooperation has indeed been favored by kin selection in natural populations.

Bacteria produce a range of molecules that are secreted from the cell and can provide a benefit to the local population of cells. Laboratory experiments have suggested that these “public goods” molecules represent a form of cooperation, favored because they benefit closely related cells (kin selection). However, there is a relative lack of data demonstrating kin selection for cooperation in natural populations of bacteria. We used molecular population genetics to test for signatures of kin selection at the genomic level in natural populations of the opportunistic pathogen Pseudomonas aeruginosa. We found consistent evidence from multiple traits that genes controlling putatively cooperative traits have higher polymorphism and greater divergence and are more likely to harbor deleterious mutations relative to genes controlling putatively private traits, which are expressed at similar rates. These patterns suggest that cooperative traits are controlled by kin selection, and we estimate that the relatedness for social interactions in P. aeruginosa is r = 0.84. More generally, our results demonstrate how molecular population genetics can be used to study the evolution of cooperation in natural populations.

Legacy and Emerging Pollutants in an Urban River Stretch and Effects on the Bacterioplankton Community

River contamination is due to a chemical mixture of point and diffuse pollution, which can compromise water quality. Polycyclic Aromatic Hydrocarbons (PAHs) and emerging compounds such as pharmaceuticals and antibiotics are frequently found in rivers flowing through big cities. This work evaluated the presence of fifteen priority PAHs, eight pharmaceuticals including the antibiotics ciprofloxacin (CIP) and sulfamethoxazole (SMX), together with their main antibiotic resistant genes (ARGs) and the structure of the natural bacterioplankton community, in an urbanized stretch of the river Danube. SMX and diclofenac were the most abundant chemicals found (up to 20 ng/L). ARGs were also found to be detected as ubiquitous contaminants. A principal component analysis of the overall microbiological and chemical data revealed which contaminants were correlated with the presence of certain bacterial groups. The highest concentrations of naphthalene were associated with Deltaproteobacteria and intI1 gene. Overall, the most contaminated site was inside the city and located immediately downstream of a wastewater treatment plant. However, both the sampling points before the river reached the city and in its southern suburban area were still affected by emerging and legacy contamination. The diffuse presence of antibiotics and ARGs causes particular concern because the river water is used for drinking purposes.