Assessment of Biological Contribution to Natural Recovery of Anthropized Freshwater Sediments From Argentina: Autochthonous Microbiome Structure and Functional Prediction

Monitored natural recovery (MNR) is an in situ technique of conventional remediation for the treatment of contaminated sediments that relies on natural processes to reduce the bioavailability or toxicity of contaminants. Metabarcoding and bioinformatics approaches to infer functional prediction were applied in bottom sediments of a tributary drainage channel of Río de La Plata estuary, in order to assess the biological contribution to MNR. Hydrocarbon concentration in water samples and surface sediments was below the detection limit. Surface sediments were represented with high available phosphorous, alkaline pH, and the bacterial classes Anaerolineae, Planctomycetia, and Deltaproteobacteria. The functional prediction in surface sediments showed an increase of metabolic activity, carbon fixation, methanogenesis, and synergistic relationships between Archaeas, Syntrophobacterales, and Desulfobacterales. The prediction in non-surface sediments suggested the capacity to respond to different kinds of environmental stresses (oxidative, osmotic, heat, acid pH, and heavy metals), predicted mostly in Lactobacillales order, and the capacity of Alphaproteobacteria, Betaproteobacteria, Gammaproteobacteria, and Actinomyces classes to degrade xenobiotic compounds. Canonical correspondence analysis (CCA) suggests that depth, phosphate content, redox potential, and pH were the variables that structured the bacterial community and not the hydrocarbons. The characterization of sediments by metabarcoding and functional prediction approaches, allowed to assess how the microbial activity would contribute to the recovery of the site.


INTRODUCTION
Petroleum exploitation in Argentina began in 1907 with the discovery of crude oil in the region of Patagonia. Throughout the 20th century, it registered a sustained growth that was accompanied by the installation of refineries to supply the domestic market and even export oil by products. The crude oil produced in Patagonia is transported to the refineries by trunk lines and sea. Basins have a local pipeline network that carries the extracted crude to the loading terminals of ships in order to transport it mainly to the operating refineries of Buenos Aires region (Riccardi, 2015). It is not uncommon that the exploitation and transportation of the crude oil generate and discharge into aquatic environments. When hydrocarbons are released into water courses, they are scavenged from the water column to the bottom through a number of processes, such as flocculation, sedimentation, and/or coagulation. This leads to concentrations in the sediment many orders of magnitude higher than in the water column, being the sediments, therefore, a potential sink for petroleum hydrocarbons (Adeniji et al., 2017). Once on the sedimentary column, hydrocarbons suffer from a wide variety of processes including: dilution, volatilization, chemical transformation, and biodegradation (Rothermich et al., 2002;Magar and Wenning, 2006;Morasch et al., 2007), as well as solid-association with sediments (Bao et al., 2013). Besides, another possible fate is the migration through the dynamic environment to other compartments, ultimately being propagated through the food chain (Xia et al., 2012).
One of the conventional remediation options for the treatment of contaminated sediments is the monitored natural recovery. MNR is an in situ approach that relies on biological, chemical, and physical natural processes that contribute to reducing the concentration, mobility, and/or toxicity of the contaminant (Foote et al., 2014). Biological and chemical processes can transform and remove contaminants. In addition, physical processes, like the deposition of fresh and clean sediment can bury and confine contaminated sediments through the formation of a natural cap. Physical processes typically occur faster than chemical and biological transformations (Lin et al., 2017), but altogether the time required for risk reduction can be relatively long. On the other hand, one of the main implications of the use of natural processes to reduce risks is the absence of building costs, often associated with other technologies.
Before selecting MNR as the remediation strategy, several lines of evidence are required in order to identify and evaluate the processes that could effectively contribute to pollutant risk reduction. For this purpose, is essential a detailed characterization involving biological, chemical, and physical analyses that support the contaminant source control, the chemical degradation, and the occurrence of physical processes. Likewise, biological trials that provide the proof of biodegradation are also essential. It results from the above that the application of MNR is the result of a thoughtful decision-making process following careful site assessment and characterization (National Research Council, 2001;US Environmental Protection Agency, 2005;Magar and Wenning, 2006). Sediments of lakes, rivers, streams, and ponds are environments with a high microbial diversity. Bacteria, archaea, and eukaryotic microorganisms play important roles in natural biogeochemical processes, like organic matter decomposition, nutrient cycling, and food webs (Gadd and Raven, 2010;Pernthaler, 2013). Particularly, bacteria influence several processes, such as the consumption of oxygen and the turnover of nutrients (Tšertova et al., 2013). Microbes have the adaptive advantage of using hydrocarbons as an optional energy source in hydrocarbon contaminated ecosystems, regulating the nutrient flux with their surrounding environment. However, microbial dynamics in association with hydrocarbon contaminated environments are still poorly studied (Reid et al., 2018).
Recently, the development of environmental DNA metabarcoding techniques have revolutionized modern biodiversity surveys increasing the confidence and the potential of biomonitoring programs in the risk assessment of pollutants (Goodwin et al., 2017;Wang et al., 2020). Although information derived from DNA sequencing is yet to be formally included in assessment programs (Bourlat et al., 2013;Aylagas et al., 2020), nowadays this methodology is receiving major attention, and several projects start to involve evaluations based in metabarcoding techniques and bioinformatics approaches (Goodwin et al., 2017;Mukherjee et al., 2017;Aylagas et al., 2020). Sequencing studies generating data on autochthonous microbial communities are vital to develop ecosystem management programs and to design efficient clean-up strategies to cope with pollutants.
Before selecting MNR as an option for site management is necessary to provide evidence that natural remediation processes are occurring. As a result, the present work was conducted with the aim to establish the microbiological processes that could be involved in an ecosystem from a tributary drainage channel of Río de La Plata estuary. The central hypothesis of this work is that the study of the diversity and functional potential of sediment bacterial community using metabarcoding and bioinformatics tools as well as its interpretation considering a detailed physicochemical and geochemical characterization of the site, will allow to identify the hydrocarbon bioremediation potential and support the selection of MNR recovery strategy.

Sampling Sites
The selected site was a tributary drainage channel of Río de La Plata estuary. Several industries, mainly petrochemicals and oil refineries, are located at both sides of this channel, close to its mouth and in a relatively small area. While discharges occurred in the past, for the last 20 years, no new wastes or discharges were created. Sediment cores and water samples were collected in February 2017 at four sites along this area near of Río de la Plata estuary. During sampling, two active urban discharges were detected on the left margin of the stream. Sampling and description of sediment cores was carried out by the Laboratory of Continental and Marine Cores (SACMa, Instituto de Estudios Andinos Don Pablo Groeber, UBA CONICET, Argentina). Undisturbed cores were retrieved from the selected sampling sites from a platform with a manually operated hammer corer with transparent PVC liners devoid of core catcher. The single core tubes measured up to 2,000 mm in length with an internal diameter of 68 mm. The coring procedure was the same for all sampling points: the corer was lowered into the sediment where it penetrated initially by its own weight. Thereafter, a driving weight (the "hammer") was lifted from the top of the corer and then released to drive the core tube into the sediment. Visual inspection after sampling allowed confirming that the water-sediment interface was kept at an undisturbed state in the liner, even after the corer was lifted to the surface. Because of the highly fluid near-surface sediments, the core was transported from the field in the core tube and split in the laboratory by means of a longitudinal section.
Core surfaces were then cleaned, properly oriented, and measured. Lithology was visually described considering a number of parameters: the nature of the sediment, grain size, wet and dry coloring according to Munsell Color (Org;2000), sedimentary structure, water content, and the presence of biogenic material. Any unusual features were noted. Major horizons or lithological units were defined for each core. Sediment subsamples were taken out from each horizon immediately after core opening. Samples were named according to the core number (1-4) and the horizon (H), being H0 the uppermost unit defined, i.e., the top core, and H6 the deepest (Supplementary Table S1). Samples were collected from the middle part of each horizon avoiding the contact between adjacent units but also the tube wall. This was particularly important to avoid the potential dragged-in contamination due to the plastic deformation of the sediment column. Samples were collected in sterile glass bottles or polypropylene tubes. According to lithological units, a total of 22 sediment subsamples were taken for subsequent analysis (Supplementary Table S1). Samples were then stored at −20°C in order to perform DNA extraction.

DNA Extraction, 16S rRNA Gene Amplification, and Illumina MiSeq Sequencing
The DNA of sediment samples was extracted using E.Z.N.A Soil DNA kit (OMEGA, United States) according to the manufacturer protocol and stored at −20°C until further analysis. DNA extracts were checked for quantity and purity with NanoDrop™ 2000/2000c Spectrophotometers (Thermo Fisher Scientific) and subsequently submitted to the sequencing platform Mr. DNA (Molecular Research LP, Shallowater, TX). The 16S ribosomal RNA (rRNA) gene V4 variable region PCR primers 515/806 for Eubacteria and Archaea with barcode on the forward primer were used in a 30 cycle PCR using the HotStarTaq Plus Master Mix Kit (Qiagen, United States) under the following conditions: 94°C for 3 min, followed by 28 cycles of 94°C for 30 s, 53°C for 40 s, and 72°C for 1 min, after which a final elongation step at 72°C for 5 min was performed. After amplification, PCR products were checked in 2% agarose gel to determine the success of amplification and the relative intensity of bands. PCR products were used to prepare a DNA library following Illumina TruSeq DNA library preparation protocol. Sequencing was performed on Illumina MiSeq Next Generation Sequencer in accordance with the manufacturer guidelines.

16S rRNA Amplicon Processing and Bioinformatic Analysis
The sequences generated in the present study were deposited in the National Center of biotechnology information 1 as the Bioproject "Natural recovery of anthropized fresh-water sediments", and are available under the accesion number PRJNA592525.

Statistics
Differences in diversity indexes and bacteria relative abundances among sediment samples were compared via statistical ANOVA and post hoc Tukey-Kramer tests, using online calculator SciStatCalc version 1.4 2 and the Statistical Analysis of Metagenomic Profiles 3 (STAMP v2.1.3) software (Parks et al., 2014), respectively. Additionally, multiple test correction Benjamini-Hochberg FDR with a q value filter of 0.05 was applied in bacteria relative abundances using STAMP software. The differences in bacterial community structure among samples based on weighted UniFrac distance metric in PCoA were assessed through the PERMANOVA test with 999 permutations in QIIME 1.

Canonical Correspondence Analysis
The data of physicochemical analysis and taxonomic assignment were used for multivariate analysis. Canonical correspondence analysis (CCA) between microbial diversity (Species data) and physicochemical parameters (Environment data) was performed using the software Canoco 4.5 (Microcomputer Power, Ithaca, NY, United States) in association with CanoDraw ® . Relative abundances (x) at class level (Species data) were transformed by arcsin (x/100) 1/2 , while the environmental variables were transformed using log 10 (x + 1). Chloroplast, unassigned sequences, and classes below 5% of relative abundance were removed for the species data matrix for further analysis. To obtain tests of independent effects of physicochemical variables, co-variables were eliminated in CCA. The effect of each environmental variable and their interaction on the sum of all canonical eigenvalues was tested by Monte Carlo permutation tests available in Canoco using 999 replicate runs to assess significance on the trace values.

PICRUSt Functional Prediction
PICRUSt (v1.1.0) was used to infer the functional potential of prokaryotic communities in terms of Kyoto Encyclopedia of Genes and Genomes 4 (KEGG) orthology and pathways, and to associate taxonomic changes with functional differences (Langille et al., 2013 into KEGG Pathways defined at levels 1-3 (Ogata et al., 1999;Kanehisa and Goto, 2000), and (v) connecting the OTUs that are contributing to each K0 of interest. STAMP (v2.1.3) was used for data visualization and statistical analyses (Parks et al., 2014). The Nearest Sequenced Taxon Index (NSTI) score, which is the sum of phylogenetic distances for each OTU between its nearest relative with a sequenced reference genome, measured in terms of substitutions per site in the 16S rRNA gene and weighted according to the frequency of that OTU, was used as an indicator for the accuracy of PICRUSt.

Core Description
Sediment cores showed negligible shortening and no evidence of vertical mixing. Water-sediment interface and top layers of sediment were conveniently preserved. Core lengths varied between 74.5 cm and 194 cm, and the number of lithological units/horizons described varied between 3 and 7. Notably, the number of horizons recognized does not necessarily correlate with core length (Figure 1). Cores 1 and 2 showed similar stratigraphic sequences. According to visual inspection, the topmost part of the sediment column was composed of massive rich organic muds. The color of these cores was black or dark greenish gray when the sample was moist and turned to olive/greenishgray after drying. Toward the base, the organic matter rich layers rich organic mud layers were interbedded with centimeter-scale layers of hydrocarbon-rich muds. Dark gray, black, and bluish-black, faint-laminated hydrocarbon-rich muds underlie these levels. The base of the cores was composed of laminated light gray silty clays with fragments of molluscs valves.
Cores 3 and 4 were stratigraphically more complex. They were composed of several horizons of hydrocarbon-rich mud layers evident because of differences in color, consistency, and proportions of silt and clay. These horizons usually presented abundant, well-preserved particulate organic matter composed mainly by plant fragments, rubble, and debris derived from domestic activities. In core 4, a carbon-rich, semi-consolidated horizon, probably derived from oil refining, was recognized. The base of the sequence of core 4 was characterized by finegrained sediments, mainly dark gray clayed silt with faint lamination and some discrete levels of sandy silt. The core length and the number of horizons of each core are summarized in Supplementary Table S1.

Physicochemical Characteristics of the Sediments
The physical and chemical characteristics [i.e., PAH, BTEX, aliphatic and aromatic hydrocarbons, ratios nC17/pristane (nC17/Pr), nC18/phytane (nC17/Ph), total nitrogen, soluble sulfate, total sulfide, available phosphorus, pH, redox potential, organic matter, cation exchange capacity, humidity, and heavy metals] that were studied in sediment samples are Frontiers in Microbiology | www.frontiersin.org 5 March 2021 | Volume 12 | Article 601705 summarized in Table 1. From these data, it can be seen that, independently of the core, the lowest values of hydrocarbons (i.e., PAH, BTEX, aliphatic hydrocarbons, aromatic hydrocarbons, total hydrocarbons) were found in samples from surface horizons (1H0, 2H0, 3H0, and 4H0). Also, most of sediment samples presented values of the nC17/Pr and nC18/Ph ratios lower than 1 (0.01-0.47) and (0.02-0.79), respectively. The ratios of surface samples 1H0, 2H0, and 4H0 were not possible to calculate due to nC17 and nC18 were below to detection limit (<50 mg/kg). pH was within the range 6.7-8.3, with the exception of samples from horizons H1, H2, H3, and H4 from core 4 (in the point of the urban discharge), in which pH values were downright acid (3.0-4.1). Redox potential was in average 200 mV for all the sediments, being indicative of a rather oxidizing environment. The content of organic matter was higher in middle sediments compared to surface and base sediments. Total nitrogen values were in average 1,500 mg/ kg for sediments samples, with the highest value in 4H2 (2,610 mg/kg). Both soluble sulfate and total sulfide tended to increase with core depth. However, total sulfide was below the detection limit (50 mg/kg) in surface sediments of cores 1 and 2 and in almost sediments samples from core 4. Available phosphorus showed the highest values in H0 and H1 sediments of core 3 (52.3 and 58.2 mg/kg), while in the other samples it was found an average of 8 mg/kg. Arsenic, cadmium, chromium, and lead, were below the permitted sediment values according to the Interim Sediment Quality Guideline [ISQG from Canadian Environmental Quality Guideline (CEQG); 5.9, 0.6, 37.3, and 35 mg/kg, respectively; data not shown].

Diversity Indexes
Illumina sequencing results showed a total of 1,322,571 filtered sequences, ranging from 20,259 to 110,528 sequences per sample. The minor number of reads (20,259) was used to normalize the OTU table to the same sample depth. A total of 1,756 observed OTUs were generated after clustering at 97% similarity level. Setting this sample depth was enough to allow the coverage of all the samples, as observed by the plateau in the rarefaction curves (Supplementary Figure S1) and the calculated values of Good's coverage index close to 1 ( Table 2). Diversity indexes (richness Shannon, Simpson, and Chao) are presented in Table 2.

Beta Diversity
Principal coordinate analysis based on weighted UniFrac distance metric was performed in order to study the differences in bacterial community structure among samples. Figure 2 shows the dissimilarities of bacterial community structures across the samples. Variances between bacterial community structures were explained by the three axes: PCo1 (35.68%), PCo2 (15.25%), and PCo3 (10.57%). From the chart, it can be seen that the samples from the surface sediments (H0) were grouped together (marked in a black circle), evidencing the similarity in their bacterial community structures (PERMANOVA p = 0,001) and the differences in comparison with the others samples (PERMANOVA p = 0,001). The samples 1H1, 2H1, 2H2, 2H3, 2H4, 2H5, 4H1, 4H2, 4H3 4H4, 4H5, and 4H6 formed a second cluster in the Euclidean space, suggesting the presence of a similar community among them. The most divergent samples were 1H2, 1H3, 3H4, and 4H4.

Canonical Correspondence Analysis of the Differences in the Taxonomy Among Sediment Samples
The relationship between the structure of bacterial communities and physicochemical variables in sediment samples were studied with a CCA, which is presented in Figure 4.  together in an area of the euclidean space represented by high available phosphorus, pH, and cation exchange capacity (group 1). The classes Deltaproteobacteria, Planctomycetia, Anaerolineae Spirochaetes, and Bacteroidia represent this group of samples. ANOVA and post hoc Tukey-Kramer statistical tests between the three groups showed that only Deltaproteobacteria, Planctomycetia, and Anaerolineae were significantly enhanced in group 1 (p < 0.05; Table 3B). Table 3C showed that Syntrophobacterales and Desulfobacterales, belonging to Deltaproteobacteria class, Anaerolineales and Methanomicrobiales, belonging to Anaerolinea and Methanomicrobia classes, respectively, were significantly enhanced in group 1. Samples 1H2, 1H3, 3H2, 3H3, and 3H4 were clustered in group 2, which correlated with high redox potential and THC. Dehalococcoidetes, Actinobacteria, Caldisericia, OP8, OPB14, and Methanomicrobia, were the most representative classes of this group. However, only the classes OPB14, with unclasiffied order_OPB14, and Actinobacteria, including members of Actynomycetales order, were found in significant (p < 0.05) and major relative abundances in this group of samples (Table 3C).
Group 3 gathered together, the deep samples of core 2 (2H2, 2H3, 2H4, and 2H5) and all hydrocarbon-contaminated samples of core 4 (4H1, 4H2, 4H3, 4H4, 4H5, and 4H6), characterized by an acid pH and by the fact of being taken from a point near to an urban discharge. Despite the fact that Alpha, Beta, Gammaproteobacteria, Bacilli, Fusobacteria, and Clostridia were the representative classes of this group, only Bacilli presented significantly higher relative abundance (p < 0.05), with almost exclusively members of Lactobacillales order. Among the orders affiliated to Gammaproteobacteria class, Pseudomonadales order was significantly more prevalent in groups 2 and 3 than in group 1.

Predicted Functional Profiles and Metagenome Contribution of the Bacterial Communities
In order to determine whether sediment samples were tractable for PICRUSt prediction, NSTI values were calculated ( Table 2). As described by Langille et al. (2013), the accuracy of PICRUSt prediction in general decreases with increasing NSTI values. For most samples, NSTI values were below 0.18, which agreed with the mid-range values calculated in soils by Langille et al. (2013). Conversely, NSTI values were greater (0.18-0.25) for 1H0, 1H2, 2H0, 3H0, 3H1, 3H2, and 4H0 samples. It is expected that these samples contain much unexplored diversity, and therefore PICRUSt results should be interpreted with caution.

Functional Prediction
Significant results of PICRUSt functional prediction are summarized in Figure 5. As indicated by the comparison between groups, sediment samples clustered in group 1 showed a significant enrichment in the number of sequences included in the category of cellular processes (Figure 5A), and in the two subcategories of cell division, cell motility, and secretion ( Figure 5B), protein folding and associated processing, translation, folding sorting and degradation (Figure 5C), and energy metabolism (Figure 5D), in relation with group 2 and 3. In this group also prevailed the pathways involved in: carbon fixation in prokaryotes and in photosynthetic organisms (Figure 5F), methane metabolism (Figure 5E), and citrate cycle (Figure 5G). Glycolysis and gluconeogenesis were enriched in sediments of group 2 ( Figure 5G) in comparison with samples from groups 1 and 3. Sediments of group 3 were significantly enriched in lipid metabolism and xenobiotics biodegradation and metabolism subcategories (Figure 5D), being the degradation pathways of benzoate, bisphenol, chlorocyclohexane and chlorobenzene, chloroalkane and chloroalkene, dioxin, ethylbenzene, fluorobenzoate, limonene and pinene, naphthalene, polycyclic aromatic hydrocarbons, and the metabolism of xenobiotics by cytocrome P450 enhanced in group 3, in comparison with group 1 and 2 ( Figure 5F). Table 4 provides an overview of several KEGG orthologs (K0s), which were involved and significantly increased in the most relevant metabolic pathways of every group. In addition, Table 4 also shows the number of sequences predicted in most relevant classes for each K0.

DISCUSSION
Monitored natural recovery is a sustainable remedy approach highly site-specific that, and, in many cases, may require years or decades to establish the amount of risk reduction desired by stakeholders (Magar and Wenning, 2006). For this reason, MNR combined with a thin-layer cap of clean sand or clean sediment over the contaminated sediment is generally selected. The capping to exceed the depth of bioactive zone, typically 5-15 cm, so that immediate risk reduction can be achieved (Fetters et al., 2020). The chemical analyses in the present work showed that the hydrocarbon concentration was below the detection limit (< 50 mg/kg) both in water and surface sediments samples. These results are likely to be related to a capping process, in which natural sedimentation permitted the physical isolation of the hydrocarbons. Interestingly this allowed the constitution of an effective hydrocarbon primary source control barrier that led to pollutant-risk reduction. In addition, the ratios nC17/pristane and nC18/phytane were used as indicators of biodegradation (Christensen and Larsen, 1993;Companioni-Damas et al., 2011). Isoprenoids (i.e., pristane, phytane) present a greater resistance to biodegradation than the n-alkanes. The fact that such ratios were generally less than one was indicative of the occurrence of advanced FIGURE 4 | Canonical correspondence analysis of microbial community composition and environmental variables. Content of available phosphorus (P), redox potential (RP), pH, polycyclic aromatic hydrocarbons (PAH), total hydrocarbons (HCT), cation exchange capacity (IIC), total sulfide (TS), water content (H), and exchangeable ammonium (NH4). CCA1: first axis, CCA2: second axis. Samples belonging to group 1, group 2, and group 3 are gathered together in green, orange, and blue squares, respectively. Samples and classes are represented with black and blue dots, respectively.
Frontiers in Microbiology | www.frontiersin.org 11 March 2021 | Volume 12 | Article 601705 biodegradation process in non-surface sediments, suggesting that natural attenuation could have happened after industrial source control (20 years ago).
The bacterial diversity results showed very similar community structures among surface sediments, with a greater and distinct diversity than that observed for the non-surface and contaminated  (Figure 3). These results are in agreement with a previous study indicating that due to intensive resuspension events, the bottom bathymetry and the patchiness of submerged vegetation, the bacterial community composition in the upper 2-5 cm tends to remain homogenous (Tšertova et al., 2013). Furthermore, different reports of biomass and enzymatic activity measurements showed that bacterial abundance and total activity are generally higher near the sediment surface and tend to decrease with depth (Haglund et al., 2003;Wurzbacher et al., 2017). These observations are probably related to the process of natural sediment formation, in which a proportion of newly settled organic matter is rapidly recycled in the surface by biological activity and subsequently transformed into secondary compounds. Consequently, vertical gradients of respiratory electron acceptors are generated (Froelich et al., 1979), being organic matter and electron acceptors the main processes that structure the microbial community (Wurzbacher et al., 2017). Multivariable analysis was performed in order to put in evidence which is the main variables that structure the community composition in the present sediment samples (Figure 4). The analysis showed yet again similarities among surface samples (Group 1), whereas the bacterial community of the non-surface sediments samples showed to be more heterogeneous, forming two separate clusters (Groups 2 and 3). Depth, available phosphorous, redox potential, and pH, but not the hydrocarbon concentration, were the environmental variables that differentiated the community composition. These results are in agreement with previous observations, which showed that the selection pressure exerted by natural filter in sedimentary coastal environments is greater than that exerted by hydrocarbons (Jeanbille et al., 2016).
In silico metagenomes obtained in the present work further support the idea of a superficial enhanced metabolic activity, as demonstrated by the enrichment of cellular processes, like cell division, cell motility, and genetic information processing in surface sediments (Figure 5). The major contribution to the diversity and extensive metabolic capacity of the surface sediments would be delivered by the presence of the orders: Syntrophobacterales, Desulfobacterales, Anaerolineales, and Methanomicrobiales (affiliated to Deltaproteobacteria, Planctomycetia, Anaerolineae, and Methanomicrobia classes, respectively; Figures 3, 4; Table 3). The enhanced carbon fixation pathways in surface samples suggest the prevalence of autotrophic processes, which could be achieved through three different pathways of carbon assimilation by microbial communities. First, Deltaproteobacteria and Anaerolineae probably perform the carbon fixation through the reductive citric acid cycle (Arnon-Buchanan cycle), since it has been predicted that these classes possess genes that encode the enzyme pyruvate-ferredoxin/flavodoxin oxidoreductase (K03737 ; Table 4). Second, Syntrophobacterales order has been demonstrated to possess the ability to fix carbon through the reductive acetyl-CoA pathway (Wood-Ljungdahl pathway; related to K00198: carbon monoxide dehydrogenase/ acetyl-CoA synthase). Finally, the classes Gammaproteobacteria, Alphaproteobacteria, and Anaerolineae could contribute to carbon fixation by the Calvin cycle as it has been predicted that these classes own genes that encode the large chain of the enzyme ribulose-bisphosphate carboxylase (EC:4.1.1.39, K01601; Table 4).
Consistent with the literature (Grabowski et al., 2005), heterotrophic and fermentative members of Deltaproteobacteria co-occurred with hydrogenotrophic methanogens (H 2 and CO 2 consumers) in surface sediments ( Table 3C). The presence of Syntrophobacterales and Desulfobacterales, together with Methanomicrobiales and the consequent enrichment in methane metabolism pathways (Figure 5E), suggests that methanogenesis could be performed by these taxa. Probably, this pathway occurs in surface sediments reducing CO 2 by formylmethanofuran dehydrogenase (EC:1.2.7.12, subunit A, K00200; Table 4). Particularly, members of the orders Desulfobacterales and Syntrophobacterales are sulfate-reducing bacteria able to grow as syntrophs in cooperation with methanogens in anoxic sediments (Morris et al., 2013). The decline of the relative abundance of Desulfobacterales with sample depth could be related to an important deterioration of "ecosystem health" (Jeanbille et al., 2016). The presence of autotrophic and syntrophic microorganisms in surface samples, involved in natural processes and commonly found in anoxic sediments suggest the establishment of a microbial community with the ability to restructure the environmental functions of the site.
Some authors have related diversity reduction to the persistence and dominance of certain taxa whose relative abundance increases with sediment depth. These persistent taxa have shown a number of advantageous features like metabolic versatility, mechanisms of resistance to environmental factors, and/or adaptations to energy limitations (e.g., substrate fermentation, carbon storage capacity, DNA repair mechanisms, sporulation ability, high ATP synthesis efficiency, and reduced genome size; Petro et al., 2017;Starnawski et al., 2017;Rissanen et al., 2019). In the present work, according to the metagenomic prediction and its abundance in acid-pH sediment samples, Bacilli could be a likely candidate for a persistent taxon (Tables 3, 4; Figure 5). Members of this class showed the genetic capacity to respond to different kinds of environmental stresses (e.g., oxidative, energy, and acid-pH stress), inducing the expression of an alternative sigma factor (σ B ) GST13 (K07550), which activates a large set of general stress genes, thereby conferring broad resistance in natural environments (Petersohn et al., 2001;Price et al., 2001). The heat shock response in Gram-positive bacteria is regulated by ctsR gene (K03708) and under osmotic stress, the presence of ABC transporter Opu A, B, C (K05845/6/7) could confer to Bacilli the capacity to resist the acid pH and hostile environment of deep sediments (Table 4). Besides, the presence of number of K0s that respond to metal ion stress were also found in Bacilli and Proteobacteria (ars genes K11811, K00537, and K06201). In agreement with this, genes implicated in the cross capacity to resist both acid pH and lead have also been reported for Lactobacillales (Giri et al., 2019). Therefore, all these genes may allow Bacilli to withstand heavy metal, acid-pH, and osmotic stress, characteristics that are actually present in deeper sediments from cores 2 and 4 (Table 4; Figure 5). Notably, in deep samples from cores 1 and 3, which are under the influence of a positive redox potential and a high hydrocarbon concentration, the presence of OPB14 class and sequences of acetyl-CoA synthase (K00198) involved in Wood-Ljungdahl pathway were enriched (Table 4). Although, due to the lack of subsurface cultured representatives, it is difficult to link the phylogenetic patterns in deep sediments, it is known that there are many groups that are capable of growing fermentatively or by acetogenesis in deep layers (Petro et al., 2017). In contrast to surface samples, the taxonomic contribution of acetyl-CoA synthase suggests that the class OPB14 could be responsible for carbon fixation in the non-surface samples, clustered in Group 2 ( Table 4).
Metabolic pathways for hydrocarbon degradation under anoxic conditions are fairly well understood (Sierra-Garcia and Oliveira, 2013). Phylogeny and metagenome prediction suggest that    Table 4]. In freshwater sediments, the oxygen is often limited to the surface and decreases within the first few millimeter or centimeter down the sediment core (Barker Jørgensen and Peter Revsbech, 1985). Below that depth, other electron acceptors, such as nitrate, iron, or sulfate become energetically favorable for the microorganisms (Rivett et al., 2008). The degradation of alkanes and aromatic hydrocarbons under nitrate-reducing conditions has been detected in in vitro assays as phenol-degrading enrichment cultures (Baek et al., 2003), in microcosms made with contaminated subsurface aquifer sediments (Wang et al., 2012a) and in the environment (Martirani-Von Abercron et al., 2016). In accordance with the present results, the bacterial diversity in the above-mentioned works was limited to subclasses of Proteobacteria, and some Gram positive bacteria (Baek et al., 2003;Wang et al., 2012b;Martirani-Von Abercron et al., 2016). In particular, members of Pseudomonadales have been demonstrated to degrade hydrocarbons using nitrate as electron donor (Mcnally et al., 1999), suggesting that the increment of this family in deeper sediments could be associated with this metabolic capacity (Table 3). Recently, some authors have described that during nitrate respiration some microorganisms can produce intracellular oxygen through a putative nitric oxide dismutase advantageous for aerobic pathways in anoxic environments (Zedelius et al., 2011;Ettwig et al., 2012). Likewise, other authors could quantify a high proportion of the corresponding nod gene in the environment (Zhu et al., 2017), suggesting that the occurrence of aerobic degradation pathways under anoxic, denitrifying conditions should not be excluded.
The characterization of sediments from a historically anthropized tributary drainage channel of Río de La Plata by environmental measures, metabarcoding, and functional prediction approaches, allowed to assess how the microbial activity would contribute to the recovery of the site. The evidence from this study suggests the formation of a natural and uncontaminated sediment cap, which allowed the exclusion of hydrocarbons from the water, thus possibly preventing pollutants entering the food chain. In addition, the findings of this research indicate the presence of microorganisms related with "ecosystem health" in natural and uncontaminated cap (H0); whereas, in deeper and contaminated sediments, hydrocarbon degradation pathways were predicted.
Overall, this study strengthens the idea that the use of metabarcoding and bioinformatics tools to assess the structure, diversity, and function of sediment bacterial communities in association with a deep analyses of site physicochemical and geochemical characteristics, proved to be useful to elucidate the possible contribution of biological processes to hydrocarbon degradation and site recovery. As MNR relies on natural biodegradation processes, this study could support the selection of MNR as a recovery strategy of the present anthropized freshwater course. This is the first time that such an approach is undertaken in the region.
Finally, this study intends to encourage the incorporation of metabarcoding techniques and bioinformatics approaches to provide new insights for MNR assessment.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/ Supplementary Material.

AUTHOR CONTRIBUTIONS
LM conceived and designed the experiments, analyzed the data and performed the computational analysis, contributed reagents, and materials, prepared figures and tables, and wrote and reviewed drafts and the final version of the paper. VS performed the laboratory analysis and sample extraction, contributed with reagents, materials, and analysis tools, prepared figures and tables, and reviewed drafts and the final version of the paper. AA analyzed DNA sequencing data, performed the computational analysis (i.e., 16S Bacteria and Archaea Standard Operating Procedure, PICRUSt pipeline), and reviewed (e.g., English language editing) drafts and the final version of the paper. BC analyzed the data and reviewed drafts of the paper. CL conducted environmental sampling, analyzed the data, and reviewed drafts of the paper. NV conceived the experiments and reviewed drafts of the paper. PM conceived the experiments and reviewed drafts of the paper. MO conceived and designed the experiments. MP conceived and designed the experiments, analyzed the data, and wrote and reviewed drafts of the paper. IM conceived and designed the experiments, analyzed the data, and wrote and reviewed drafts and the final version of the paper. All authors contributed to the article and approved the submitted version.

FUNDING
This work was financially supported by the Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET) and YPF Tecnología (Y-TEC) under the Agreement: COL 01662016-2016. The funder, Y-TEC, had the following involvement with the study: study design and sample collection.