The metagenomic analysis show that b-lactamases were present in each metagenome analyzed, ranged from 0.0003% (percent of reads assigned to b-lactamases from the total number of analyzed reads) in a ocean metagenome (accession number SRS582462) to 0.0335% in a human gut metagenome (accession number SRS056259). The abundance of b-lactamases per environment is in average a 0,01% in soils (n=80 metagenomes, including agricultural and non-agricultural soils), 0.0068% in glacier (n=11), 0.0046% in fresh water (n=11), 0.002% in ocean metagenomes (n=22), 0.0121 in human gut (n=63), 0.0049% in cow gut (n=27 metagenomes, including feces and rumen samples) and 0.0092% in wastewater treatment plant environment (n=19). Soils and glaciers metagenomes show the higher diversity of b-lactamases genes (number of different b-lactamase genes), followed by wastewater metagenomes. Is important indicate that some b-lactamase genes can present different variants, for instance, blaOXA gen present 257 different variants in the EX-B database, but all the variants hits obtained in a given sample/environment are assigned to blaOXA gene, which explains the differences between the number of b-lactamases in the EX-B database and the number of b-lactamase genes in the richness graph. The total b-lactamase relative abundance show that soils, human gut and wastewater metagenomes are the most b-lactamase enriched environments.
The presence of b-lactamases, grouped according molecular classes, show differences according environment (Fig. 2); thus, class A b-lactamases are dominant in non-agricultural soils, cow (including feces and rumen metagenomes) and human gut environments (more than 70% in each case), class B b-lactamases are dominant in agricultural soils, fresh water, oceans and wastewater treatment plant environments (between 40 to 60% in each environment), class C b-lactamases are more represented in agricultural soils and fresh water environments (20% approx. in each environment) and class D b-lactamases being more abundant in fresh water, wastewater treatment plant and glacier environments (from 15% to 30% approx.). A more detailed picture of b-lactamases in the environment is presented in table 1; where b-lactamase genes were analyzed according their occurrence in a particular environment throughout an indicator species analysis (Dufrene 1997
). According these results, 16 b-lactamase genes show a high faithfulness of occurrence in non-agricultural soils, other 16 b-lactamase genes show a high faithfulness of occurrence in agricultural soils and 14 b-laxctamase genes show a higher probability to be present in the wastewater treatment plant environment. Only four b-lactamases (blaEBR, CfxA, HGI and mecA) were found to be highly present in the human gut environment (p<0.005) instead of other analyzed environments. Interestingly, when the environments are grouped according level of anthropogenic impact (anthropogenic impact level 1= wastewater treatment plant; level 2= human gut; level 3= agricultural soils, fresh water, oceans and cow gut; level 4= non-agricultural soils and glaciers), the less impacted environments show a high level of faithfulness of occurrence of b-lactamase genes (50 genes for anthropogenic impact level 3 and 4) than the observed in the more anthropogenic impacted environments, where only 13 b-lactamase genes show faithfulness of occurrence.
b-lactamase hits obtained by BLAST were used to construct metagenome distance matrices and b-lactamase gene networks. When all the b-lactamase hits were used to construct the gene network (Fig. S1), different cluster were clearly observable; the same trend was obtained when genes poorly described at gene level (i.e. identified only such b-latamase, ESBL gene or class A/B/C/D) were removed from the analysis (Figure 3). The network analysis indicate that each clusters harbor metagenomes almost exclusively related to a given environment; thus five cluster representing soils (including a differentiation between agricultural and non-agricultural soils), human gut, animal gut and wastewater were identified. Cluster analysis performed on b-lactamase gene networks, indicated that samples from the same environment are more tightly connected than samples from different environments (Table S2). This is clearly visible in Fig. 3, where nodes from a given environment (soils, human gut, cow gut and wastewater treatment plant environment) show more connections between them than the number of connections with nodes related to other environments.
Explicar asortividad, clustering coeficient, y algun otro..............
In order to test if sample geography influences on b-lactamase gene content, nodes present in Fig. 3 were presented according their geographic origin (Fig. 4). Clearly, the results obtained here indicated that geography is not correlated with both b-lactamase gene content and diversity.
The BLAST analysis performed on metagenomic reads indicate the presence of b-lactamase genes in different environments, but also can indicates how close are the metagenomic reads related to the reference database. To elucidate it, we choose four clinically important b-lactamases such as blaCTX-M, blaGES, blaOXA and blaTEM and the hits obtained by BLAST were grouped in five groups according the percent of identity to the same genes in X-B database (identities from 50 to 59, 60 to 69, 70 to 79, 80 to 89 and 90 to 100 percent) (Fig. 5). The results indicated that in each case, most of the hits show a low percent of identity; thus, in average, between 43.7 (in human gut metagenomes) to 69.9% (bovine feces) of blaOXA hits show a percent of identity between 50 to 59 percent. Hits with higher percent of identity to our database were obtained from human gut, wastewater treatment plant environment, fresh water and glacier metagenomes. In the case of blaGES, 44.2% (non-agricultural soils) of hits until 81.2 % (glaciers) of hits show a percent of identity between 50 to 59 percent. In addition all the obtained hits in human gut show the same range of identity, but it is based only in six hits. In the case of blaCTX-M and blaTEM, close to 40% of the hits show a percent of identity of 50 to 59 percent. In some cases, a higher percent of identity can be observed for most of the hits in a given environment (blaCTX-M gene in human and bovine environment; blaTEM in ocean, human gut or bovine rumen metagenomes) but in those cases, the total number of hits obtained was very low to be taken into account. Real higher percent of identity (between 90 to 100%) were only detected in some specific cases that include a high number of hits; thus, blaCTX-M gene shows hits with the high percent of identity in samples obtained from human gut (n=14) and hits with high percent of identity to blaGES gene were found in wastewater treatment plant environment (n=105). Hits related to blaTEM show a percent of identity to our EX-B database, higher than 79%, in human gut, wastewater treatment plants, oceans, fresh water, glacier and agricultural soil metagenomes.
The target of antibiotic resistance genes in natural environments should be considered, specially if we consider the evidence that these genes have a long evolutionary history in natural environments (Aminov 2009
, Hall 2004
, Garau 2005
, Baltz 2008). Thus, natural environments appear such as a reservoir of potential ARGs to pathogenic bacteria (Berglund 2015
, Versluis 2015
) that cannot be ignored.
The results obtained in this study show the high content of b-lactamase genes in each analyzed environments, but a deep analysis indicate important differences in the diversity of those genes. Thus, soil and glacier metagenomes exhibit the higher diversity of b-lactamase genes, followed by wastewater metagenomes. The high diversity of b-lactamases in soils is in accordance with a previous metagenomic-based study of antibiotic resistance genes in the environment (Nesme 2014
), in which soil metagenomes showed the most diverse pool of ARGs than human feces, ocean, cow and chicken gut and artic snow metagenomes. In relation to glaciar metagenomes, a recent study showed a big diversity of ARGs in samples obtained from glacier from different geographic locations, including Antarctica, Alaska, Himalaya, Greenland and others (