Labour sharing promotes coexistence in atrazine degrading bacterial communities

Microbial communities are pivotal in the biodegradation of xenobiotics including pesticides. In the case of atrazine, multiple studies have shown that its degradation involved a consortia rather than a single species, but little is known about how interdependency between the species composing the consortium is set up. The Black Queen Hypothesis (BQH) formalized theoretically the conditions leading to the evolution of dependency between species: members of the community called ‘helpers’ provide publicly common goods obtained from the costly degradation of a compound, while others called ‘beneficiaries’ take advantage of the public goods, but lose access to the primary resource through adaptive degrading gene loss. Here, we test whether liquid media supplemented with the herbicide atrazine could support coexistence of bacterial species through BQH mechanisms. We observed the establishment of dependencies between species through atrazine degrading gene loss. Labour sharing between members of the consortium led to coexistence of multiple species on a single resource and improved atrazine degradation potential. Until now, pesticide degradation has not been approached from an evolutionary perspective under the BQH framework. We provide here an evolutionary explanation that might invite researchers to consider microbial consortia, rather than single isolated species, as an optimal strategy for isolation of xenobiotics degraders.

Microorganisms in nature usually co-exist as communities whose complexity is under the influence of local environmental conditions and interactions between individuals. Spatially structured environments, such as biofilms, are more prone to support coexistence of multiple species consuming the same resource because access to the resource and its metabolic by-products is conditioned by the structure of the environment 1,2 . Interactions of various natures between bacterial species that stabilize the community diversity can then arise and persist, but are most likely to disappear when the structure is disrupted 3,4 . In mixed liquid cultures, competition for a single resource is supposed to lead to competitive exclusion 5 , although very recent research has shown that cross-feeding may stabilize coexistence on a single carbon source 6 . If the breaking down of the resource includes both private and public goods, natural selection could favor the emergence of dependency between species, and therefore coexistence, through adaptive gene loss according to the Black Queen Hypothesis (BQH) 7,8 .
Morris et al. formulated three key characteristics for a function to follow the BQH. The function allowing consumption of the resource must be i) costly, ii) essential and iii) leaky. Also, Estrela et al. 9 and Oliveira et al. 10 recently specified the conditions that may favor the emergence of interdependency, which include low privatization and low cost of secretion. Atrazine is one of the most heavily applied herbicides worldwide which is relatively persistent and mobile in soil, and whose degradation products can be traced in soils decades after application 11,12 . Atrazine can be mineralized by microorganisms and is a source of nitrogen [13][14][15][16] . Its entire biodegradation pathway is known (Fig. 1) and can be split into two parts: the upper part consisting in the dechlorination of atrazine and the removal of the two lateral chains from the s-triazinic cycle, followed by the lower part consisting in the complete degradation of cyanuric acid. Genes (atz and trz families) involved in its mineralization are most often located on plasmids or on catabolic cassettes delimited by insertion sequences 17,18 . Because of the cost of maintaining plasmids, atrazine genes can be easily lost after a small number of generations in culture media without atrazine 19 . However, in nitrogen limited environments, atrazine degrading capacities must be conserved because mineralization of atrazine leads to the essential delivery of nitrogen 20 . Also, metabolites, such as desisopropylamine, aminoethanol and desethylamine, are produced and released in the environment during atrazine degradation 21 . Consequently, atrazine biodegradation is (i) conferred by genes costly to maintain, (ii) essential under N-limited conditions because it provides nitrogen supply, and (iii) benefits to other members of the community because metabolites leak in the environment. Atrazine biodegradation is therefore a good candidate function to test the BQH predictions, and in particular the emergence of dependency between species in a spatially unstructured environment.
Chelatobacter sp. SR38 22 , Pseudomonas sp. ADPe 19 , Arthrobacter sp. TES 23 and Variovorax sp. 38R 24 are all originally soil bacteria isolated from arable soils exposed to atrazine. They are all able to at least partly degrade either atrazine or its intermediary products (Fig. 1), and represent four different genetic repertoires for atrazine degradation. They have all been modified to resist to a different combination of two antibiotics in order to evaluate their frequency in mixed culture on double-antibiotics plates. Based on their genetic repertoire for atrazine degradation, it is quite straightforward to predict that in a liquid culture medium with atrazine as the unique nitrogen source, Arthrobacter sp. TES and Pseudomonas sp. ADPe must coexist because Arthrobacter sp. TES will degrade atrazine to supply its nitrogen needs and produce cyanuric acid that will be used by Pseudomonas sp. ADPe as its nitrogen reservoir. However, predicting what would happen in a more complex scenario where multiple competing strains evolve in environments containing various nitrogen sources is not trivial. Here, we experimentally questioned whether spatially unstructured atrazine-containing media could support coexistence of multiple atrazine degrading species through BQH mechanisms.
We propagated for ~100 generations four-species consortia, in triplicate, in seven minimal media supplemented with citrate in excess, as the carbon source, and either one among three nitrogen sources: atrazine, cyanuric acid: the metabolic product of the upper part of the atrazine degradation pathway, or ammonium sulfate; or 2-way and 3-way combinations of the above mentioned nitrogen sources, keeping the N molarity constant. We found that out of the seven media, all species coexisted in the four ones containing atrazine, whereas Arthrobacter sp. TES went extinct in the three atrazine-free media (Fig. 2a). The extinction of Arthrobacter sp. TES in cyanuric acid supplemented media was expected as it does not possess the metabolic pathway to degrade cyanuric acid, however its extinction also in ammonium sulfatesupplemented media, is more surprising since it can use ammonia as sole nitrogen source. Intriguingly, the equilibrium frequencies of the three other species, besides Arthrobacter sp. TES, are quite similar after ~100 generations in six out of the seven media, with Pseudomonas sp. ADPe dominant at ~10 7 CFU.mL −1 and the two remaining species at ~10 6 CFU.mL −1 , the exception being the ammonium sulfate -cyanuric acid supplemented medium where Variovorax sp. 38R appeared to be dominant at ~10 7 CFU.mL −1 .
We then assessed using multiplex PCRs the presence of the atrazine degradation genes in the evolved consortia across the seven different media. We found that in the four media containing atrazine, all initially present atz or trz genes were retained in the evolved consortia, preluding the conservation of the atrazine mineralisation potential in these environments (Fig. 2b). In the atrazine-free media, only atzA, responsible for the dechlorination of atrazine, a step that releases chlorine but does not provide any nitrogen input to cells, and atzD, responsible for the first step of the degradation of cyanuric acid, were conserved in the three environments while atzC and trzD were specifically kept in the ammonium sulfate -cyanuric acid supplemented medium. Since in that specific medium, atzC does not confer any advantage to its carrier, Chelatobacter sp. SR38, it has likely been carried along because of its colocalization with trzD on the same plasmid 25 , responsible for the degradation of cyanuric acid. Rapid evolution, probably through the loss of atz genes containing plasmids, therefore led to the loss of the atrazine degradation function in atrazine-free environments, confirming the genetic load of maintaining atrazine genes in atrazine uncontaminated environments.
We then evaluated the atrazine mineralization potential, using 14 C-atrazine either labelled on the ethylamino chain or on the s-triazinic cycle, of the consortia that have evolved in the four atrazine containing environments (Fig. 3) and compared them to the ancestral ones. Interestingly, we observed a significant increase of the mineralization potential of the ethylamino chain in evolved consortia compared to ancestors in all but the medium which contained all three nitrogen sources. The observed gain is however much stronger in the atrazine only supplemented medium (from 48%; CI 95% = [41.9-54.1] to 71.4%; CI 95% = [69. 4-73.4]). When evaluating the mineralization potential of the s-triazinic cycle, results are quite contrasted with a significant increase of the mineralization potential in the atrazine only supplemented medium (from 48.3%; CI 95% = [36.5-60] to 86.9%; CI 95% = [81.6-92.2]), and significant decreases or no change in media supplemented with combinations of two or all three nitrogen sources. Altogether, these results indicate that the upper part of the atrazine degradation pathway (atzA, trzN, atzB and atzC) has been improved in the evolved consortia, presumably leading to increased intermediary metabolite production, such as aminoethanol, ethylamine or hypoxanthine 21 . While the upper pathway is constitutively expressed, studies have shown that the lower pathway was repressed by the presence of ammonium in the environment 26 . It was therefore quite likely that the stronger accumulation of intermediary metabolites and their consumption by consortium members led to an increased concentration of ammonium in the environment explaining at least partly the decreased mineralization potential of the s-triazinic cycle in evolved consortia compared to ancestors.
The impressive boost of atrazine mineralization potential observed in the medium supplemented only with atrazine led us to focus on phenotypic and genotypic changes that occurred in this specific environment. We therefore isolated members of the evolved consortia, evaluated individually their genetic repertoire using multiplex PCRs, and assessed their corresponding growth characteristics using Bioscreen TM assays. We found that only Arthrobacter sp. TES was able to grow in isolated culture, while the three other members of the evolved consortia depend on Arthrobacter sp. TES to grow in this medium (Fig. 4a). This acquired dependency for Chelatobacter sp. SR38 is explained by loss of part of its atz genes repertoire, namely the upper pathway genes atzA and atzB (Fig. 4b). We then wanted to characterize the interactions between members of the evolved consortia. We therefore reconstructed all possible duos and trios, and compared growth performances of each member to their performance in the evolved four-species consortia and in monoculture over 3.5 days (Fig. 4c). A relatively simple and straightforward interaction scheme can be drawn from these results (Fig. 4d): Variovorax sp. 38R, Chelatobacter sp. SR38 and Pseudomonas sp. ADPe growths are clearly supported by Arthrobacter sp. TES. However, Arthrobacter sp. TES is penalized in duos or trios with Chelatobacter sp. SR38 and Pseudomonas sp. ADPe, in the absence of Variovorax sp. 38R. Interestingly, Variovorax sp. 38R by itself has no positive effect on Arthrobacter sp. TES, but it counteracts the negative impacts of the two other consortium members on Arthrobacter sp. TES. Therefore, coexistence of the four members might be favoured, even though two of the three direct beneficiaries exert antagonistic relationships against the public goods provider, because the third one acts as a buffering intermediary.
Here, we have witnessed the establishment of dependencies in a four species artificial consortium after evolution in an unstructured, atrazine only supplemented medium during ~100 generations. According to the BQH, dependency was set up via adaptive gene loss, and one specific member of the consortium, Arthrobacter sp. TES  www.nature.com/scientificreports www.nature.com/scientificreports/ ensuring the costly part of the function, became the provider of public goods to the three other beneficiaries, Variovorax sp. 38R, Chelatobacter sp. SR38 and Pseudomonas sp. ADPe. In many cases, pesticide degradation in nature is thought to occur through the involvement of microbial consortia, rather than being the corollary of a single species 14 . Our study provides strong experimental evidences, using the herbicide atrazine as a case study, that the division of labour, sensu Giri et al. 27 , between populations in microbial communities could also be considered for the biodegradation of xenobiotics. Also, we anticipate that future research should focus on the bioaugmentation with stabilized and tightly structured microbial degrading consortia as an effective solution for in situ bioremediation of sites polluted with recalcitrant compounds. ). An aliquot of each cycle culture was mixed with glycerol (30% final concentration) and preserved at −80 °C.

Species composition.
We were able to easily identify the four strains in the ancestor and evolved "community" lines because each of them was double-resistant to a different combination of two of the following antibiotics: rifampicin 100 mg/L, kanamycin 50 mg/L, spectinomycin 100 mg/l, streptomycin 100 mg/L. To do so, each bacterial strain was grown in MS medium containing either atrazine or cyanuric acid as the nitrogen source and then spread on solid rich medium (TY) containing the antibiotic for which it was to be resistant. Resistant growing colonies were collected and cultivated in their specific liquid MS medium to ensure that they have kept their metabolic ability towards atrazine. This procedure was repeated to confer a second resistance to each strain. Pseudomonas sp. ADPe is Strep + /Rif + , Chelatobacter sp. SR38 is Kan + /Strep + , Arthrobacter sp. TES is Spec + /Rif + and Variovorax sp. 38R is Spec + /Kan + .
Ten-fold successive dilutions of T1 and T23 "community" lines were done and inoculated onto LB agar plates supplemented with different antibiotic combinations. After 3-days of growth at 28 °C, CFU were counted across 3 dilution levels.
Atrazine-degrading genetic repertoire. Multiplex PCR targeting atzA/atzD/trzD and atzC/trzN/atzB genes were used to determine the atrazine degradation repertoire of each strain of the ancestral and evolved "community" lines. 3 or 4 isolated colonies from the ancestral and evolved "community" lines were resuspended in water and served as DNA template for PCR reactions. PCR reactions were conducted using the following primers Atrazine mineralisation activity. Ancestors and "community" lines that have evolved on the four media containing atrazine were inoculated at ~0.001 OD 600nm in 200 µL fresh medium supplemented with 4000 dpm final of 14 C atrazine labelled either on the ethylamino chain or on the s-triazine cycle in triplicate. Cultures were then placed at 28 °C. A whatman ® 3 mm Chr paper soaked with barite (saturated solution 0.37 M) was placed and sealed on the top of the 96-wells plates and replaced periodically during 3.5 days. Dried papers were fluorography printed on Storage PhosphorScreen (Molecular Dynamics ® ) for 2 days. Screens were scanned by a phosphorimager (Storm 860 Molecular Imager). Indirect reading of released radioactivity intensity was done by ImageQuant 5.2 software, and referred to a standard curve. Mineralisation percentages of lines were determinate with respect to initial radioactivity dosage realized on 200 µL of each culture medium with a Beckman ® liquid scintillation counter. Significant differences between ancetral and evolved mineralization percentages were assessed using Tukey's test (p < 0.05). population dynamics characteristics of the ancestral and evolved line. The population dynamics of ancestral and evolved "community" lines were measured using Bioscreen TM . All lines were inoculated at ~0.001 considered as significant (p < 0,05). (d) Interactions scheme drawn from mono-and co-cultures of the four community-evolved strains in the atrazine only supplemented medium. Interactions types were defined based on growth performances in co-cultures experiments. Positive (+) or negative (−) interactions are represented by arrows. Arrows point towards the impacted strain. The dotted arrow expresses the conditional positive effect of Variovorax sp. 38R on Arthrobacter sp. TES in the presence of the two other strains.