Biogeochemical Responses and Seasonal Dynamics of the Benthic Boundary Layer Microbial Communities during the El Niño 2015 in an Eastern Boundary Upwelling System

The Eastern South Pacific coastal zone is characterized by seasonal and interannual variability, driven by upwelling and El Niño Southern Oscillation (ENSO), respectively. These oceanographical conditions influence microbial communities and their contribution to nutrient and greenhouse gases recycling, especially in bottom waters due to oxygenation. This article addresses the seasonal hydrographic and biogeochemical conditions in the water and sediments during El Niño 2015. Bottom water active microbial communities, including nitrifiers, were studied using amplicon sequencing of 16S rRNA (cDNA) and RT-qPCR, respectively. The results of the hydrographic analysis showed changes in the water column associated with the predominance of sub-Antarctic Waters characterized by warmed and low nutrients in the surface and more oxygenated conditions at the bottom in comparison with El Niño 2014. The organic matter quantity and quality decreased during fall and winter. The bottom water active microbial assemblages were dominated by archaea (Ca. Poseidoniales) and putative ammonia oxidizing archaea. Active bacteria affiliated to SAR11, Marinimicrobia and Nitrospina, and oxygen deficient realms (Desulfobacterales, SUP05 clade and anammox) suffered variations, possibly associated with oxygen and redox conditions in the benthic boundary layer. Nitrifying functional groups contributed significantly more during late fall and winter which was consistent with higher bottom water oxygenation. Relationships between apparent oxygen utilization nitrate and nitrous oxide in the water support the contribution of nitrification to this greenhouse gas distribution in the water. In general, our study suggests that seasonal oceanographic variability during an El Niño year influences the microbial community and thus remineralization potential, which supports the need to carry out longer time series to identify the relevance of seasonality under ENSO in Eastern Boundary Upwelling Systems (EBUS) areas.


Introduction
The coastal zone of central Chile is marked by high productivity, associated with upwelling events driven by the anticyclonic gyre seasonality, characteristic of Eastern Boundary Upwelling Systems (EBUS). In this area, upwelling occurs seasonally fertilizing the surface layers with high nutrient and oxygen deficient subsurface waters due to the presence of the oxygen minimum zone (OMZ) on the continental shelf [1]. Additionally, the upwelled water is characterized by the accumulation of greenhouse gases, such as CO 2 and N 2 O, usually associated with pelagic and benthic anaerobic microbial activities [2][3][4]. However, the EBUS coastal areas in the eastern South Pacific (ESP) are also subject to great spatiotemporal dynamics, associated, for example, with the variability of upwelling-favorable winds, causing meso-and sub-mesoscale structures (eddies, filament meanders) [5]. Among other sources of oceanographic variability, coastal-trapped waves could be a relevant driver of interannual variability, influencing the oceanographic dynamics in this coastal area, for example, during "El Niño" when a deepening of the thermocline, and warmer and more oxygenated waters are characteristics at the surface layer [6,7]. During 2015-2016, strong temperature anomalies were observed in the equatorial region during the "Godzilla El Niño" event, in comparison with previous events since 1950 [8], impacting the biogeochemical conditions in the water, from moderate to strong, between early 2015 and mid-2015 at Callao coasts off Perú [9]. These authors reported higher temperatures, oxygen, pH and nitrate but lower silicate and phosphate concentrations at the surface layer during "El Niño". The shifts in the nutrient reservoir in the surface water and oxygenation influences the phytoplankton community structure and thus primary productivity at different upwelling areas along the ESP coast [10][11][12]. Changes in the OMZ and nutrient availability are key factors in shaping the microbial communities' structure and activities [2], including key microbial communities and their contribution to organic matter respiration, producing CO 2 and N 2 O greenhouse gases. Aerobic versus anaerobic processes influence the N 2 O accumulation in the water column, which is aerobically produced through nitrification as a by-product of the ammonia oxidation to nitrite, carried out by bacteria and archaea in marine ecosystems [13,14]. N 2 O could also be anaerobically produced or consumed by denitrification, since it is an intermediate product of the dissimilative reduction of nitrate, generating molecular nitrogen, a significant process producing greenhouse gases and their potential exchange in the atmosphere in OMZs [15]. In addition, changes in oxygenation and in substrate availability generate profound microbial composition effects as determined in coastal areas of the ESP which are also subject to seasonal and interannual variability [16][17][18]. El Niño was found to influence benthic marine communities like macro and meiofauna [19,20] and sulfur oxidizing large filamentous bacteria Ca. Marithioploca. [21], due to changes in bottom water oxygenation, shifting sediment biogeochemistry [22,23]. However, little is known about the impact of seasonality and other sources of variability present during El Niño on the active microbial community in the central area of the ESP and its potential connection with water column biogeochemistry. Herein, we hypothesize that water column biogeochemical conditions will change during selected representative months per season during El Niño 2015, influencing benthic boundary layer microbial community composition and active aerobic ammonia oxidizing microbes in the continental shelf off Valparaíso, central Chile.

Materials and Methods
Valparaíso Bay in central Chile (32 • 9 S, 71 • 6 W) is situated in an Eastern Boundary Upwelling System influenced mainly by the seasonal upwelling of Eastern Subsurface Waters (ESSW) characterized by cold, nutrient rich and low oxygen waters, associated with the southern portion of the Oxygen Minimum Zone (OMZ) [24]. The alongshore southerly winds field versus the northerly winds field drives the upwelling dynamics and water column biogeochemical conditions, including the OMZ spatiotemporal variability, which is also affected by interannual conditions, such as El Niño, that deepens the thermocline [6]. In fact, ENSO cycles affect coastal primary productivity in the Valparaíso Bay [11]. Moreover, the Valparaíso Bay area is characterized by the influence of the Aconcagua River and anthropogenic effects associated with urban sewage outfalls and harbor activities, contributing to organic matter enrichment and benthic habitat heterogeneity in the coastal zone [25,26].
The oceanographic station is situated in the middle of the continental shelf area ( Figure S1). This Station corresponds to the Oceanographic Time Series off Concón Water & Sediment (SOCAS, in Spanish) visited seasonally between 2012 and 2015 [25,26]. In this study, the SOCAS oceanographic station time series was visited during selected months between 2014 and 2015. The first stage of the sampling consisted in identifying the location and depth using a GPS Echo-sounder Garmin ® . Temperature, salinity, and oxygen profiles were obtained by a CTDO (Seabird-19 plus). Discrete water sampling was carried out using Niskin bottles at four-six depth levels (0, 25, 50, 75, 100 and 125 m). Water was collected for oxygen (Winkler bottles) and ammonium (40 mL Schott Glass Bottles, DURAN) determinations in triplicate, immediately incubated with respective reagents following standard procedures and recommendations [27,28]. Seawater for nitrate, nitrite and phosphate analysis was filtered through a GF/F filter (Whatman) using syringes and filter holders of 25 mm diameter, into acid-clean Nalgene bottles (125 mL). These bottles were kept on board the ship in the dark and with ice packs until arrival at the lab (<5 h), and then they were frozen at −20 • C in the freezer until they were analyzed. Seawater for Chlorophyll-a (1 L) analyses was collected and kept in the dark until filtration in the laboratory. Seawater for greenhouse gases (GHG = CO 2 , CH 4 and N 2 O) determination was collected in triplicates in 20 mL vials. The GHG samples were poisoned to avoid further biological activity, using 50 µL of saturated mercuric chloride. These samples were stored in the dark at room temperature until they were analyzed. Seawater (60 mL) for microbial cells was concentrated by filtration onto 0.22 µm hydrophilic PVDF filters, (GVWP02500, Millipore, Burlington, MA, USA) using sterile syringes and filter holders of 25 mm diameter during January, April, August, and December 2015. The filters were preserved with RNAlater ® solution (Ambion, Austin, TX, USA), stored in ice packs until arrival at the lab and transferred to a −80 • C freezer until RNA extraction.
Non-disturbed sediment samples for biochemical analysis were collected from independent replicates using a gravity corer (internal diameter of 50 mm and 1000 mm length). Redox potential (Ehnhe) and pH measurements at surface sediment layer (0-2 cm) were determined on board using a platinum standard combination electrode with a calomel internal reference (SGTM, Mettler Toledo). Additionally, surface sediments for total organic matter, Chlorophyll-a, phaeopigments, carbon and nitrogen stable isotope ratios were collected and stored in Ziploc bags, kept with ice packs until arrival at the lab and finally they transferred to a −20 • C freezer until sediment chemical analyses. Sediments for DNA analysis were collected for benthic microbial community characterization. Approximately 0.5 g of surface sediment was subsampled from the gravity cores using plastic vials (duplicates) and preserved with RNAlater ® solution (Ambion, Austin, TX, USA), kept with ice packs until arrival at the lab, and transferred to a −80 • C freezer until RNA extraction.

Water Column Chemical Analyses
Oxygen concentrations were determined by a modified Winkler method and using a DOSIMAT for microtitration [29]. This data was also used to calibrate oxygen data registered using a CTDO (Conductivity-Temperature-Depth-Optical). Ammonium was determined fluorometrically [28] using a Turner design ® fluorometer (Sunnyvale, CA, USA). Nutrients were analyzed by spectrophotometry based on standard methods with a nutrient automatic analyzer [30]. Chlorophyll-a was determined in discrete seawater samples according to the method described by Holm-Hansen et al. [31]. GHG were determined based on the headspace method [32] and analyzed by a gas chromatograph (GC-2014 Greenhouse, Shimadzu, Kyoto, Japan) equipped with an electron capture detector (ECD) and a flame ionization detector (FID). Three-point calibration curves using helium, air and a standard of 600, 5 and 1 for CO 2 , CH 4 and N 2 O mixture were run for absolute quantification (Scotty gas mixture; Air Liquid Co., Paris, France).

Sediment Chemical Analyses
The total organic matter (TOM) was determined using weight differences before and after ignition at 475-500 • C (four hours) [33]. TN and carbon (δ13C) and nitrogen (δ15N) isotopes in surface sediments were analyzed by mass spectrometry (Thermo Scientific Delta V Advantage IRMS with EA-2000 Flash Elemental analyzer). Carbon and nitrogen isotope values were reported according to Pee Dee Belemnite (VPDB) and air standard and analyzed at the Laboratory of Biogeochemistry and Applied Stable Isotopes (LABASI, PUC), Chile. Sediment-bound chlorophyll-a and phaeopigments were extracted from 1 g duplicate subsamples of wet sediment using 90% acetone following a procedure described by Gutiérrez et al. [34]. The total sulfides content was determined using colorimetric methods following Cline [35].

Microbial Communities' Analyses
Sediment DNA extraction was conducted with the PowerSoil DNA Isolation Kit (Mo-Bio Laboratories, Carlsbad, CA, USA) in accordance with the manufacturer's specifications. Water RNA was extracted using Ambion ® RNA extraction kit (AM1560) according to manufacturer's protocol with the addition of a mechanical disruption step with 200-µm-diameter zirconium beads (Low Binding Zirconium Beads, OPS Diagnostics, Lebanon, NJ, USA) for two rounds of 30 s (~3000 rpm) using a Mini-Beadbeater-8™ (Biospec Products, Bartlesville, OK, USA). DNA and RNA quantification were determined by using fluorescence (Qbit 2.0). For the analyses of RNA, the DNA was removed using the TURBO DNA-free™ kit (Applied Biosystems, Foster City, CA, USA) and RNA was tested for residual DNA by bacterial 16S rRNA gene PCR amplification following standard procedures [16]. RNA was reversetranscribed with random primers provided by the ImProm-II™ Reverse Transcription System (Promega Corp., Madison, WI, USA).
Bacterial and archaeal 16S rRNA gene were amplified using primers 28F (5 -GAGTT TGATCNTGGCTCAG) and 519R (5 -GTNTTACNGCGGCKGCTG) and Arch349F (5 -CCCTAYGGGGYGCASCAG) and Arch958R (5 -YCCGGCGTTGAMTCCAATT), respectively, and sequenced using high-throughput Illumina Miseq by a commercial laboratory (MrDNA, Shallowater, TX, USA). The 16S rRNA and rDNA gene sequences were checked, and quality standards were revised using the Mothur software v1.35.1 [36]. Sequences with short reads <200 bp, with more than one undetermined nucleotide (N) and with a maximum homopolymer length of eight nucleotides were deleted. Chimeric sequences were also identified using the Chimera UCHIME algorithm [37] and removed. The 16S rRNA and rDNA gene sequences retrieved were taxonomically classified using the automatic software pipeline SILVAngs available from https://www.arb-silva.de/ [38]. The sequences were deposited in the European Nucleotide Archive (ENA) under the project accession PRJEB31703 (experiments accession ID ERX3252743-ERX3252756).
Alpha diversity indexes were analyzed based on Operational Taxonomic Units (OTUs) using PRIMER (7.0.11) with the PERMANOVA add on PRIMER [39]. Rarefaction curves, Heatmap, and Venn Diagrams were analyzed with R using the ampvis2 package [40]. The comparison of the active bottom water versus surface sediment microbial community (16S rRNA versus 16S rDNA) was performed to identify a core active and present microbiome from the benthic boundary layer using dissimilarity analyses between categorical factors (sediment/water) considering the Similarity Percentages-species contributions (SIMPER). In addition, the comparison between active bottom water microbial communities from Upwelling (January and December) versus non-Upwelling seasons (April and August) was conducted to identify dissimilarities between seasons using SIMPER. Principal Coordinated Analyses (PCoA) were generated based on the Bray Curtis similarity matrix using square root transformed reads of the active bottom water microbial communities using OTUs. The relationships between water and sediment environmental conditions with bottom water active microbial community structure were visualized using correlation analyses (Pearson) overlaid on PCoA. Distance based linear models (DistLM) were estimated using a similarity matrix to evaluate the relative contribution of the bottom water conditions in the microbial community structure variability.
Aerobic ammonia oxidizing archaea Nitrosopumilus-like and ammonia oxidizing bacteria from Betaproteobacteria (βAOB) were quantified using RT-qPCR based on the ammonia monooxygenase subunit A (amoA) gene and the nitrite oxidizer, Nitrospina-like bacteria using 16S rRNA gene. First, cDNA was generated using the reversed primer Arch-amoAR [41], amoA-2R [42] and NitSSU_286R [43] with the ImProm-II Reverse Transcription System (Promega, Madison, WI, USA). Quantification was performed following previous studies [44], using the Brilliant II Sybr Green qPCR master mix (Stratagene) kit and PCR reactions consisting of a 20-µL, with approximately 5 ng of cDNA (2 µL) as a template and 0.4 µM of primer (423F-Arch-amoAR, amoA 1F-2R and NitSSU_130F-NitSSU_286R primers), in triplicate for samples, standards and negative controls. The PCR amplification protocol consisted in an initial activation and denaturation for 10 min at 95 • C, then 40 cycles of denaturation (30 s at 95 • C), annealing (30 s at 56 • C for amoA and 57 • C for 16S rRNA gene), and a final extension (40 s at 72 • C). A melting curve was also included. The standard curve was carried out using a PCR product (4-4 × 10 7 copies/µL) with the following efficiency (109-114%) and r 2 (>0.997). The qPCR was run on a MX3000P Stratagene Real-Time PCR System.

Oceanographic Features Observed
Hydrological conditions determined during time series are shown in Figure 1. The water column off Valparaíso was characterized by the presence of two water masses, the sub-Antarctic Water (SAAW) and the Equatorial Subsurface Water (ESSW) according to Silva et al. [24], the first usually located at the surface layer (<25 m depth) and the ESSW at deep and bottom waters (>50 m depth). During our study, waters with high salinity (34.6 PSU) and density (>26.2 Sigma-t), but with low temperature (<12 • C) and oxygen conditions (<2 mL/L) were associated with ESSW water mass mainly during 2014 and January 2015 ( Figure 1, Table S1). However, during April, August, and December 2015 an increment of SAAW contribution in comparison with ESSW was observed in the water column using temperature versus salinity (T/S) diagram analyses ( Figure 2).   Figure S2), one of the strongest El Niño events registered in the ESP [8]. Moreover, besides the El Niño event of 2015, an extremely strong storm and meteo-tsunami event [45] were reported twelve days before our sampling collection in August 2015.
Phytoplanktonic biomass analyses based on discrete chlorophyll-a surface water determinations, (Table S1 Supplementary Material), and satellite image analyses are depicted in Figure S3. Chlorophyll-a maxima were observed in the spring months of (September and November) 2014 and between spring and summer (October and December) 2015. Higher nutrient availability for phytoplankton in the study area was reported to be correlated with upwelling and other allochthonous sources, for example from the Aconcagua river [11]. In fact, during our sampling, the Aconcagua River showed a discharge increment in late 2015, the highest of the three-year analysis ( Figure S4, Table S1 Supplementary Material).
However, only slight salinity decreases were observed at the surface waters ( Figure 1). Our results are consistent with the primary productivity observed during the last decade at a shallower coastal area off Valparaíso associated with the interannual upwelling dynamics, including the influence of ENSO warm and cold phases in longer time series during 1986 and 1996 [11].

Biogeochemical Conditions in the Water Column
The water column biogeochemical conditions are shown in Figure 3. High nitrate (>20 µM) and phosphate (>2.5 µM) were predominant mostly during 2014 in subsurface waters (>15 m), whereas 2015 sampling periods presented high values restricted to the bottom waters during December. Lower nutrient concentrations were observed between April and August 2015. These latter months were also characterized by higher nitrite and ammonium concentration in subsurface and bottom waters, supporting a potential enhancement of remineralization in the water column. Greenhouse distribution in the water column is present in Figure S5. Nitrous oxide presented a higher concentration at the oxyclines and deep subsurface waters during August and December 2015. Carbon dioxide presented higher concentration in bottom waters, whereas low methane concentrations were detected throughout the water column. The apparent oxygen utilization (AOU) was estimated using the difference between the expected oxygen saturation in the water, based on salinity and temperature, and the observed oxygen concentration ( Figure 3). The AOU determination indicates a higher organic matter respiration potential during 2014 in comparison with the 2015 sampling months, in particular during late fall (April). In addition, nitrate and phosphate relationships were estimated to evaluate potential shifts from Redfield ratio~16:1 and potential nitrogen loss processes based on N* Nitrogen deficit (Figures 3 and 4). During August 2015, low N/P and N* were determined in subsurface waters. Negative N* values are associated with anaerobic processes such as denitrification and anammox [46,47], and both nitrogen loss processes are known to be active in the eastern South Pacific (ESP) as well as in the central-south coastal upwelling areas [17,48,49] and other oxygen-deficient regions [15,50]. Moreover, the nitrous oxide emissions contribute to nitrogen loss from the ocean, associated with intermediate products of denitrification and of aerobic ammonia oxidation [51]. The contribution of aerobic ammonia oxidation, the first step of nitrification, can be inferred from the relationship between nitrous oxide and AOU and nitrate [52]. In our study, nitrous oxide presented a positive relationship with nitrate and AOU, and nitrate was also positively related to AOU (Figure 4, Table S2).  (Table S2) indicate the predominance of nitrification as a source of nitrous oxide associated with more oxygenated waters. However, partial denitrification could also be active in high nitrate and low oxygen concentration ESSW (Figure 4). Our results were supported by nitrous oxide spatiotemporal dynamics in the upwelling area time series off Concepción (Centro de Investigación Oceanográfica en el Pacífico Sur-Oriental, COPAS time series) situated at central-southern Chile [3,4].

Sediment Conditions during the Selected Months
The sediments of the study area were mainly muddy (60%-67%) and characterized by high organic matter content (TOM 6%-8%) during January and April; meanwhile, low values (TOM < 4.9%) during the following months were registered ( Figure 5, Table S3 Supplementary Material). In addition, a decrease in the %TOM associated with shifts in the organic matter composition was observed characterized by low C:N values in summer (<8.8, January and December) in comparison with the other months ( Figure 5A). Total nitrogen varied between 23.74-33.89, showing higher values during April (Table S3). A decrease in the chlorophyll-a content was observed from January towards December ( Figure 5A). In contrast, phaeopigments were variable and increased towards the last months sampled, having a maximum during winter (August). Higher redox and pH values, and lower total sulfides were observed during the winter months, and parameters were associated with organic matter composition ( Figure 5B). The carbon and nitrogen stable isotopes composition of the sediment ranged between −24.45-−23.85 for δ13C and between 6.59-6.88 for δ15N (Table S3). These isotopic signals were consistent with those values reported in other sites in central Chile [53], suggesting an important contribution of marine organic matter derived from primary production. Indeed, the higher δ15N values in August and lower values in January (Table S3) could be associated with sediment transport and the resuspension process of organic matter in coastal areas [54]. It is important to mention that the sedimentary conditions were highly variable during the study period considering oceanographic condition changes. Previous reports in the study area indicate a strong relationship between the food supply from primary production and the physical properties of the seafloor structuring the macrobenthos biodiversity and distribution [26]. Chlorophyll-a, used as a proxy for fresh organic matter content versus phaeopigments as an indicator of degradation products [55], indicate that January and April presented higher, and fresher, organic matter content, in comparison with the other months. The C:N values determined during our study were lower (<9.5) than previously reported for the same study area during 2013 (>10.9, [25]), suggesting a slightly higher lability during 2015. Total sulfide indicates that the potential of sulfate reduction, the predominant microbial remineralization process in marine sediments [56], was active during our sampling time, but its potential accumulation decreased during August. This latter month presented a higher oxygenation and greater redox potential suggesting an increase in aerobic remineralization processes and a potential bioturbation [34]. On the other hand, the δ13C and δ15N signatures indicate that during August a higher contribution of allochthonous material (terrigenous), more typical of shore stations closer to Aconcagua river discharge areas, are expected according with previous studies [25,53]. During August, the distinct isotopic signal in the sediments could be accounted for by modified oceanographic conditions on the water column and bottom due to the occurrence of meteo-tsunamis recorded 0n the coast of central Chile [45].

Microbial Community Dynamics in the Different Biomes
Amplicon 16S rRNA sequencing results indicate that sediments were characterized by a rich and diverse microbial community in comparison with bottom water (Table S4), especially during January and December (upwelling period). Rarefaction curves indicate that a higher richness was expected in the sediments, since a plateau was not reached in our sequence results ( Figure S6A). In general, most of the microbial communities present in the surface sediments were active in the bottom water since >90.7% OTUs were shared in both biomes ( Figure S6B), indicating a significant bottom boundary layer exchange. A total of 64 phyla were detected, but only 52 were distributed in more than three months and were listed according to SILVA 132 classification (Table S4). A heatmap analyses helped us to visualize the top abundant phyla that were variable in the different biomes and sampling months ( Figure 6). Moreover, the specific high frequent OTUs are plotted in Figure 7.  The Marine Group II, currently affiliated within the order Ca. Poseidoniales, accounted for 25% of the total reads and was predominant during January and December (>42.4%, Figure 6). Putative ammonia oxidizing archaea related to Ca. Nitrosopumilus and Ca. Nitrosopelagicus from the order Nitrosophaerales contributed 22% of the total reads and were enriched in the water mainly during August, representative of winter ( Figure 6). These predominant archaeal taxa were comparable with other marine areas globally, characterized by the abundant Ca. Poseidoniales known as a key heterotrophic order that peaks following phytoplankton blooms, and by Nitrosophaerales as putative relevant ammonia oxidizers in the ocean [14]. For example, our results agreed with Tara ocean metagenomic expedition reports showing the ubiquity of these group and the high frequency of Ca. Poseidoniales supported by its genetic versatility, including its ability to harvest light [57]. SIMPER analyses indicate that many archaeal OTUs, specially associated with putative ammonia oxidizers, were favored during April and August linked to non-upwelling period (Table S5).
On the other hand, Proteobacteria was the predominant phylum of the bacteria domain (Table S4), and was characterized by 17 classes with Gammaproteobacteria, Deltaproteobacteria and Alphaproteobacteria dominant in terms of sequences ( Figure S7A). Uncultured Proteobacteria associated with single OTUs presented a significant contribution in the different biomes ( Figure S7B). Many of these taxa, corresponded to cosmopolitan marine pelagic microbial groups, like SAR11 (Ca. Pelagibacter), ZD0405 and SAR86 (Oceanospirillales), and Marine Group B (SAR324). The SIMPER analyses indicate that these groups contributed greatly. mainly during the non-upwelling period (Table S5), supporting previous reports of longer historical time series studies such as at the North Sea, during winter [58], or oligotrophic conditions in different marine areas, such as The Bermuda Atlantic Time-series Study (BATS) and The San Pedro Ocean Time-series (SPOT) (reviewed in [59]).
In general, our predominant bottom water microbial groups were comparable with assemblages associated with SPOT deeper waters (>150 m depth), e.g., Thaumarchaeota, Marine Group A, SAR86 (Oceanospirillales), Nitrospina spp [59,60]. On the other hand, our results support other relevant functional microbial groups typically found in the OMZ of the Eastern South Pacific (ESP) [61,62], such as sulfur oxidizing bacteria from the SUP05 cluster ( Figure S7B) and the anaerobic ammonia oxidizers (anammox bacteria) Ca. Scalindua, ( Figure S8). In addition, Deltaproteobacteria OTU sva0081 (Desulfobacterales) presented a high number of reads in the sediment ( Figure S7B). These marine bacteria are a H 2 oxidizer key functional group of anaerobic organic matter degradation in marine sediment [63], known for rapid growth under micro-aerophilic conditions [64]. In our study, Desulfobacterales were characterized by various OTUs (22) and is an order known to be enriched during low-oxygen conditions further south in the ESP area off Concepción [16]. Besides hydrogenotrophic, other bacteria of this order are sulfate reducers [65], a significant remineralization process in marine sediments [56,61].

Nitrifying Microbial Groups and Their Predominance in the Bottom Water during Non-Upwelling
Keystone microorganisms associated with nitrogen cycling, such as putative archaea ammonia oxidizers (Ca. Nitrosopumilus, Ca. Nitrosopelagicus, and others) and Nitrite oxidizers affiliated with Nitrospina spp. presented a higher contribution in the water biome, and during non-upwelling conditions based in SIMPER analyses (Table S5, Figure 7 and Figure S8). The specific quantification of active putative ammonia oxidizers associated with Nitrosupumilus-like and βAOB, and the nitrite oxidizers Nitrospina-like (RT-qPCR amoA and 16S rRNA genes, respectively) supports the results based on amplicon 16S rRNA sequencing. Higher RT-qPCR counts indicate that the nitrifying functional groups presented a higher activity during the non-upwelling period, i.e., April and August, in comparison with the upwelling period January and December ( Figure 8).
In addition, the greater activities determined during April are consistent with the increment of oxygen ( Figure 5) and ammonium concentration in bottom waters (Figure 8). Nitrospina-like was the most active nitrifying microbial group followed by Nitrosupumiluslike and then βAOB, as observed in the time series upwelling area off Concepción in the ESP [18,44]. Nitrospina spp. was considered as a key marine microbial group usually associated with putative ammonia oxidizing archaea in network analyses in the coastal ocean water column, such as Monterey Bay [66,67]. Moreover, both functional groups are known to be active at extremely low oxygen concentration in deep waters of the ESP upwelling area off Concepción [68]. Besides the biogeochemical contribution of specific microbial taxa, the benthic boundary layer microbial communities have a relevant ecological interaction with benthic fauna, such as meiofauna (>18 µm and <1 mm), which are widespread in marine sediments around the world [69]. Meiofauna are dominated by free-living nematodes usually grazing on microorganisms and reported to influence organic matter remineralization processes and benthic bacteria composition [70]. These communities have been reported in high abundances in the study area [25], showing a variability controlled by organic matter content and oxygen availability [19]. Bioturbation by meiofauna and macrofauna could drive benthic microbial processes changes associated with nitrification, denitrification, and the contribution of greenhouse gas production [71]. Oxygenation could have profound changes in the benthic boundary layer conditions and benthic fauna responses.
The influence of environmental conditions on the active bottom water microbial community structure were explored using Principal Coordinated Analyses (PCoA, Figure 9). The PCoA analyses indicate that the bottom water active microbial community variability at OTU level was differentiated during our sampling months. This is evident in the PCO1 axes, accounting for 57.6% of the variability. Bottom water temperature and δ15N signature in the sediment were the variables having a greater correlation with the microbial community variability represented by the PCO1 and PCO2 axes, respectively (Table S6). Moreover, temperature was a significant (p = 0.04) explanatory variable of the variability based on DistLM (Table S7). Other variables showing correlation with changes in the microbial community structure, were organic matter composition and remineralization indicator parameters, such as organic matter % and chlorophyll in the sediments, particulate organic carbon in the water, and remineralization conditions such as redox (Table S5). Our results are consistent with previously published results indicating that temperature was a predominant explanatory variable influencing microbial community structure [60].

Conclusions
The coastal area of the eastern South Pacific presented highly dynamic oceanographic conditions that influenced the water column and sediment biogeochemistry during El Niño 2015. Benthic boundary layer microbial communities and their contribution to organic matter remineralization were seasonally variable, favoring aerobic groups such as nitrifying assemblages during fall and winter, related to positive redox conditions. The stimulation of nitrification could also be evidenced in the water column resulting in the accumulation of nitrous oxide, a potent greenhouse gas. Further studies are needed in coastal areas to determine regulatory factors of marine microbiome biodiversity and their resulting contribution to biogeochemical cycles and benthic ecology during warm and cold phases of El Niño.
Supplementary Materials: The following are available online at https://www.mdpi.com/2073-4 441/13/2/180/s1, Figure S1. Map of the study area., Figure S2: Sea surface temperature (SST) plotted together with the Oceanic Niño Index (ONI), Figure S3: Chlorophyll-a analyses based on satellite images and detected biomass in-situ during our time series, Figure S4: Aconcagua river discharge monthly average, Figure S5: Greenhouse gases distribution in the water column, Figure S6 (Table S1: General conditions in the water column and biogeochemical conditions at the surface and bottom waters, Table S2: Spearman rank correlation analyses of the environmental variables. Significant correlations (p < 0.05) were depicted (pink positive, green negative), Table S3: Specific analyses associated with sediment conditions studied during 2015., Table S4: Table S4. Alpha diversity analyses associated with the 16S rRNA sequencing reads, ENA experiment ID and Archaea and Bacteria phyla (% associated with each domain)., Table S5: Similarity Percentages-species contributions (SIMPER) analyses associated with the 16S rRNA gene iTag sequencing reads using different categorical factors, Table S6: Variables correlations (Pearson) associated with the first axes of the Principal Coordinate Analyses derived from OTUs Bray Curtis similarity, Table S7: Distance based linear models (DistLM) results based on pelagic microbial community OTUs).