www.nature.com/scientificreports OPEN Enhanced carbon-sulfur cycling in the sediments of Arabian Sea oxygen minimum zone center Received: 10 January 2018 1 1 2 1 Svetlana Fernandes , Aninda Mazumdar , Sabyasachi Bhattacharya , Aditya Peketi , 3,4 2 1 2 5 Accepted: 21 May 2018 Tarunendu Mapder , Rimi Roy , Mary Ann Carvalho , Chayan Roy , P. Mahalakshmi , 1 6 3 2 Published: xx xx xxxx Rheane Da Silva , P. L. Srinivasa Rao , Suman Kumar Banik & Wriddhiman Ghosh Biogeochemistry of oxygen minimum zone (OMZ) sediments, which are characterized by high input of labile organic matter, have crucial bearings on the benthic biota, gas and metal fluxes across the sediment-water interface, and carbon-sulfur cycling. Here we couple pore-fluid chemistry and comprehensive microbial diversity data to reveal the sedimentary carbon-sulfur cycle across a water- depth transect covering the entire thickness of eastern Arabian Sea OMZ, off the west coast of India. Geochemical data show remarkable increase in average total organic carbon content and aerial sulfate 2− reduction rate (J ) in the sediments of the OMZ center coupled with shallowing of sulfate methane SO4 transition zone and hydrogen sulfide and ammonium build–up. Total bacterial diversity, including those of complex organic matter degraders, fermentative and exoelectrogenic bacteria, and sulfate-reducers (that utilize only simple carbon compounds) were also found to be highest in the same region. The above findings indicate that higher organic carbon sequestration from the water-columns (apparently due to lower benthic consumption, biodegradation and biotransformation) and greater bioavailability of simple organic carbon compounds (apparently produced by fermetative microflora of the sediments) are instrumental in intensifying the carbon-sulfur cycle in the sediments of the OMZ center. Oxygen minimum zones (OMZs) of the global ocean encompass perennially oxygen-depleted (<20 µ M dissolved O ) subsurface water masses spanning approximately 100 and 1000 mbsl , and have critical implications for the 2–4 biodiversity/productivity of the ocean biome . The Eastern tropical North Pacific, Eastern tropical South Pacific and Arabian Sea harbor some of the most prominent OMZs defined by secondary nitrite maxima (> 0.5 µ M − 5,6 7 NO ) , arising out of active heterotrophic denitrification in the water-column. Among the OMZs, the central 8 9 Arabian Sea contains the thickest (~1.2 km vertical expanse) and most intensely oxygen-depleted water mass , 6 2 4 with the total OMZ area covering approximately 3.2 × 10 km . The Arabian Sea OMZ originates from high productivity in the euphotic zone juxtaposed with high respiratory oxygen demand for organic matter degra- dation and poor intermediate water ventilation caused by the landlocked geography of the northern portion of the sea . High productivity, in turn, is driven by upwelling of nitrate-rich water mass due to monsoonal cur- rent and winter convective mixing . Even though the most intense upwelling occurs along the western bound- ary of Arabian Sea , the zone of denitrification, remarkably, occurs in the central and eastern regions of the 5,8,13,14 sea . This is potentially attributable to the eastward advection of organic matter from the western side. As 6 2 15 the oxygen-deficient water mass impinges upon a large area (1.15 × 10 km ) of the upper continental slope of the Arabian Sea, high primary productivity, lack of oxic respiration and zooplankton grazing result in organic matter enrichment in the underlying sediments, which in turn become a hotspot of carbon-sulfur sequestra- 17–22 16 tion and cycling . Enhanced preservation of labile organic matter in anoxic sediments fuels an intricate network of microbe-mediated anaerobic processes such as organic compounds fermentation, methanogene- 23–26 sis/methanotrophy, sulfate reduction and sulfide oxidation . However, compared to the other major OMZs 27–32 of the world relatively less information is available on the sediment carbon-sulfur biogeochemistry of the 1 2 CSIR-National Institute of Oceanography, Dona Paula, Goa, 403004, India. Department of Microbiology, Bose Institute, Kolkata, 700054, West Bengal, India. Department of Chemistry, Bose Institute, 93/1 APC Road, Kolkata, 700009, India. Indian Institute of Engineering Science and Technology, Shibpur, Howrah, 711103, West Bengal, 5 6 India. Prima Dona Homes, Dona Paula, Goa, 403004, India. Gujarat Energy Research and Management Institute, Gujarat, 382421, India. Correspondence and requests for materials should be addressed to A.M. (email: maninda@ nio.org) or W.G. (email: email@example.com) Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 1 www.nature.com/scientificreports/ Figure 1. Bathymetry map showing the coring sites (red dots) of this study. The sampling locations are magnified in the inset. 24,33–36 Arabian Sea (moreover, most of the latter studies are centered off the Oman and Pakistan margins) . So far as the overall geomicrobiology of Arabian Sea OMZ sediments are concerned, apart from few recent studies on reflection of climatic oscillation within the microbiome and microbial characterization coupled with enzymatic diversity , such studies are rare that integrate solid-phase/pore-fluid chemistry with comprehensive microbial diversity data obtained using high-throughput metagenomics. Furthermore, reports of this kind of holistic study of in situ carbon-sulfur cycling, involving many discretely located sediment cores from across a transect cov- ering the entire thickness of the OMZ, are practically nil. In this scenario, the present cross-disciplinary study characterizes sedimentary carbon-sulfur chemistry and total bacterial diversity along eight discrete gravity cores (∼3 m long) collected from water-depths spanning 225 and 1275 mbsl, covering the entire thickness of the Eastern Arabian Sea OMZ. These cores collected on-board RV Sindhu Sankalp (SSK42), and designated as SSK42/1 through 8, were located across a ∼30 km transect spread across the upper continental slope, off the West coast of India (see Fig. 1 and Supplementary Table 1). Results of this geomicrobiological investigation revealed marked water-depth-dependent spatial variations in (i) depth-integrated sulfate reduction rate; (ii) total bac- terial diversity, as well as the diversities of complex organic matter hydrolyzing bacteria, fermentative bacteria and sulfate-reducing bacteria; (iii) organic carbon content; and (iv) shallowing of the depth of sulfate-methane transition zone; and (v) pore-water sulfide buildup. In view of the potential expansion of global OMZs in the recent times, and in the light of the ecologically damaging H S expulsion events reported from sediments off the 39–41 Peru-Chile and Namibian shelves , our findings assume critical oceanographic significance. Results and Discussion High organic carbon sequestration boosts microbial diversity in sediments underlying the OMZ center. Consistency of the (TOC/TN) ratios (13.6 ± 3.4) and δ C (−20.6 ± 1‰) values measured molar TOC across the water-depth transect (Supplementary Fig. 1 and Supplementary Table 2) suggests that marine produc- tivity is the dominant source of organic matter in the OMZ sediments. Total organic carbon (TOC) content in the investigated cores vary from 0.32 to 5.0 (wt %) (Supplementary Fig. 2b and Supplementary Table 2). Average TOC contents of the cores exhibited marked water-depth dependent variation: cores located close to the upper (SSK42/8; 225 mbsl) and lower (SSK42/1; 1275 mbsl) boundaries of the OMZ have lower average TOC content than those within the OMZ (Fig. 2a). This is attributable to higher consumption and remineralization of organic 43,44 45,46 matter due to relatively higher availability of dissolved oxygen and abundance of benthic fauna near the edges of the OMZ. In contrast, higher TOC content in the sediments underlying the OMZ center (SSK42/5 and 6, located at 580 and 530 mbsl, respectively) could be resultant of greater organic matter flux from the perennially oxygen-depleted water column and diminutive activity of the benthic biota , which not only leaves organic mat- 22,47–49 ter largely unutilized but also contributes to the preservation of its biochemical composition . Biodiversity rich, chemoorganoheterotroph dominated microbial communities were found to occur till the bottoms of all SSK42 sediment cores studied for their microbiology (Fig. 3). A maximum of 3374, and a minimum of 332, operational taxonomic units (OTUs, or microbial species-level entities) were identified by PCR-amplification and sequencing (up to a data-throughput level that ensured plateau in the rarefaction curve) Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 2 www.nature.com/scientificreports/ Figure 2. Trends of water-depth dependent variation of key geochemical and microbiological parameters of the Arabian Sea OMZ sediments. (a) Variations in the average TOC content of the studied sediment cores (light blue curve), and the average OTU-count of all Bacteria (purple curve) in the sediment cores. The error bars in the average TOC plot indicate standard deviation. The range of OTU-count of all Bacteria in the relevant sediment cores can be seen from Supplementary Table 9. (b) Variations in the bottom-water oxygen concentration (red triangle), sulfate flux J (black circles), and total SRB-OTU-count (blue dots) in the studied SO4 sediment cores. The range of SRB-OTU-count in the relevant sediment cores can be seen from Supplementary Table 9. The light blue shaded region represents the thickness of oxygen minimum zone. of the V3 regions of all bacterial 16S rRNA genes present in the individual metagenomes isolated from the total 87 sediment datapoints studied across SSK42/1, 3, 5, 6, 7 and 8 (Supplementary Tables 3–8). Whereas, Bacteria-specific V3 primers generated PCR products of desired size (∼200 bp) from all the metagenomes inves- tigated, Archaea–specific primers failed to generate any amplicon from 45 of them, plausibly due to the presence of very low archaeal cell number in the corresponding sediment samples. Notably, deep shotgun sequencing and analysis of whole metagenomes isolated from the sediment datapoints of SSK42/1, 5 and 6 showed that in all the metagenomes for which Archaea–specific primers failed to amplify V3 amplicons, archaeal protein-coding and 16S rRNA-encoding reads were ≤0.5% and ≤0.0001% of all annotatable reads, respectively. In contrast, in the same metagenomes, bacterial protein-coding and 16S rRNA-encoding reads were ≥99.5% and 0.1–0.5% of all annotatable reads, respectively. In this way, more than half of the sedimentary communities under investigation had no matching information on their archaeal diversity; the present study therefore focused only on the bacterial component of the communities. Taxonomic composition of communities were found to vary considerably along the sediment-depths of indi- vidual cores, while a number of core-specific variations (for instance, high diversity of Chloroflexi in SSK42/5, Firmicutes in SSK42/5 and 135-270 cmbsf of SSK42/8) were also apparent across the water-depth transect. However, amidst the variabilities, certain features of global consistency unified microbial community architec- tures across the sediment cores; these included the presence of diverse Gammaproteobacteria, Alphaproteobacteria and Bacteroidetes, many members of which are chemoorganoheterotrophic and capable of transforming wide 50–52 range of low to high molecular weight carbon compounds . Although the metagenome isolation technique used in the present study did not discriminate between cellular DNA from living and dead microbial cells, or for that matter between cellular DNA and extracellular DNA adhering to the sediment particles, the absence of down-depth decrease in the OTU-count for most of the phyla/classes present suggests that major fractions of the bacterial diversities revealed across the SSK42 cores are viable in situ. In the same way as that observed for the average TOC contents of the cores, average OTU count of the sediment cores also exhibited remarkable variation with bottom water oxygen concentration: across the water-depth transect, average OTU count is highest in the Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 3 www.nature.com/scientificreports/ Figure 3. Phylum-/class-level assortment of the taxonomically-classifiable bacterial OTUs identified at the different sediment-datapoints of SSK42/1, 3, 5, 6, 7 and 8. (A) SSK42/1; (B) SSK42/3; (C) SSK42/5; (D) SSK42/6; (E) SSK42/7; (F) SSK42/8. Color-code for the identification of taxa is common for A through F, and is given at the bottom of the figure. sediments of the OMZ center (i.e. for SSK42/5) and relatively lower in the two extremities (Fig. 2a). These findings indicate that dissolved oxygen in the bottom water, and consequentially the TOC content of the sediments, are the key determinants of microbial species richness in the OMZ-underlain sediments. Diversity of sulfate-reducing bacteria (SRB) is higher in the sediment horizon of the OMZ center. Almost every sediment sample explored for microbiology across SSK42/1, 3, 5, 6, 7 and 8 were found to encom- pass at least one OTU ascribable to known sulfate-reducing bacterial taxa (Supplementary Table 9). In this context it is noteworthy that the actual SRB diversity may well be higher for some of the sample sites because the present analyses were based exclusively on bacterial diversity data (see the Methods section “Species-richness estimation for sulfate-reducing bacteria” below) and the sulfate-reducing genera such as Desulfurococcus, Desulfurolobus (phylum Crenarchaeota), and Archaeoglobus (phylum Euryarchaeota) of the domain Archaea may also be pres- ent. In all the sediment cores, maximum SRB-OTUs were identified within 30 cmbsf; subsequently, however, multiple phases of rise and fall in SRB diversity were detected down the sediment-depths (Fig. 4; Table 1). Trends of variation in SRB-OTU distribution along the sediment-depth of each core were analyzed by determining the probability density function for every individual zone of fluctuation. In doing so, a common pattern of SRB diversity distribution along the surface-sediment to core bottom trajectory was revealed for the entire region. This involved (i) an upper zone of exponential decay, (ii) a middle zone encompassing a Gaussian distribution, and (iii) a lower zone of exponential rise, in SRB-OTU count (Fig. 4), defined by the probability density functions given in equations (1), (2) and (3), respectively. In these equations, the constants A , t , w, and x were computed 1 1 c simultaneously from the iterative fitting of the functions with the OTU-count data (to optimize the probability density functions they were sampled for up to 4000 iterations till the χ values reached their minima). P(X) = Ae 1 (1) Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 4 www.nature.com/scientificreports/ Figure 4. Zones of SRB-OTU-count fluctuation along the sediment cores of SSK42/1 (row A ), SSK42/3 (row B), SSK42/5 (row C), SSK42/6 (row D), SSK42/7 (row E) and SSK42/8 (row F), as defined by probability density functions. The histograms shown in the six panels of column I were derived based on the normalized values of the SRB-OTU-counts obtained at the individual sediment-depths of the six coring-stations. The differently- colored lines represent the probability density functions fitting to the different zones of SRB-OTU-count uc fl tuation. Expanse of the different zones along the sediment-depths of the six coring-station have been shown separately in columns II through IV, where all panels under Column II depict Zone 1 (red lines, which follow exponentially decaying distribution); all panels under column II depict Zone 2 (blue lines, which follow Gaussian distribution; and all panels under column IV depict Zone 3 (green lines, which follow exponential growth). Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 5 www.nature.com/scientificreports/ Zone 2 Zone 1 Zone 3 () xx − Name of the c x x −2 − A t t (yA =+ y) e 1 (ye =+ y) w (yA =+ y) e 1 coring-station 1 1 0 0 0 w π /2 D = 0–84 cm D = 84–125 cm SSK42/1 — 2 2 χ = 185.031 χ = 0.134 D = 8–65 cm D = 65–180 cm D = 180–285 cm SSK42/3 2 2 2 χ = 1190.057 χ = 62.253 χ = 6.376 D = 0–45 cm D = 15–160 cm D = 140–260 cm SSK42/5 2 2 2 χ = 0.624 χ = 373.333 χ = 151.524 D = 2–135 cm D = 60–220 cm D = 220–275 cm SSK42/6 2 2 2 χ = 69.652 χ = 73.105 χ = 77.225 D = 30–165 cm D = 45–165 cm D = 150–240 cm SSK42/7 2 2 2 χ = 11.919 χ = 6.355 χ = 0.447 D = 2–105 cm D = 15–180 cm D = 105–210 cm SSK42/8 2 2 2 χ = 16.333 χ = 4.704 χ = 6.001 a b c Table 1. Segmentation and characterization of the three zones of SRB-OTU distribution along the sediment- depths of the six SSK42 cores. Expanse of a given zone over a sediment core (i.e. the range of sediment-depth b 2 covered by the zone) is denoted by D. e minim Th al χ value that optimized the distribution function of a particular zone in a particular core segment is given next to information on D. e O Th TU-distribution function that characterizes a particular zone is shown in parenthesis below the zone name. () xx − A −2 P(X) = e w w π/2 (2) P(X) = Ae 1 (3) Notwithstanding the down-depth fluctuations revealed for the total OTU-count of SRB, genus-level clas- sification of the SRB OTUs showed that members of the deltaproteobacterial orders Desulfobacterales and Desulfovibrionales are present in all the sediment communities explored for their microbiology. Furthermore, overall SRB diversity for individual sediment cores, across the sampling transect, was found to increase from SSK42/1 through SSK42/5, and then decline SSK42/6 onward (Fig. 2b). The mean, median, maximum, as well as minimum, SRB-OTU-count for a sediment core (as determined from the data presented in Supplementary Table 9) were all found to be highest for SSK42/5. Microbial sulfate reduction intensifies in the sediments underlying the OMZ center. Consumption of pore-water sulfate via organoclastic sulfate reduction (OSR) and anaerobic oxidation of meth- ane (AOM), together with fermentation and methanogenesis, are the major drivers of carbon-sulfur cycle in 53,54 anoxic sediments . In all the SSK42 cores, pore-water sulfate concentration was found to decrease linearly 2− with sediment-depth (Fig. 5; Supplementary Table 2) with the SO concentration gradients ranging between −1 0.015 (for SSK42/8) and 0.11 (for SSK42/7) mM cm . In this context it is noteworthy that in the gravity cores of the present study, potential loss of core-top sediment layers (that have highest sulfate reduction rates) is unlikely to have any major influence on the sulfate concentration gradients because over the top few centim- 2− eters of any sediment package there is virtually no decrease in SO concentration due to the short diffusion 31,54 2− distance involved . Sediment-depth-integrated sulfate reduction rates (J ) calculated for the individual SO4 cores from their sulfate concentration profiles were found to range between 0.0008 (for SSK42/8) and 0.0113 (for −2 −1 SSK42/7) mmol cm yr (Fig. 2b). These values represent net sulfate reduction rates that are generally lower 35 2− − 27 than the radiotracer based ( S ) gross sulfate reduction rates , which, unlike the net rates, incorporate SO4 /HS the high sulfate reduction rates of the core-tops. Potential underestimations relative to the expected gross rates 2− notwithstanding, the present J values are well within the range of net sulfate reduction rates reported for SO4 sediments underlying the eastern Arabian Sea OMZ, and are also comparable with those obtained for other 20 2− prominent OMZ sediments . Notably, the J values calculated along the present water-depth transect reach SO4 their maxima within the center of the OMZ (Fig. 2b) and decrease towards the boundaries of its vertical expanse due to the relatively higher oxygen enrichment at the sediment water interface, a phenomenon that has parallels 55,56 reported from in vitro studies on SRB . On the other hand, the enhanced sulfate reduction rate observed in the sediments of the OMZ center could be attributable to potentially greater availability of labile organic matter sup- porting OSR and methanogenesis . While the abundance of labile organic matter (simple short chain fatty acids, such as acetate, lactate, succinate and formate) and H in sediments are known to determine the extent of OSR, 58,59 methanogenesis and consequently the depth of sulfate-methane transition zones (SMTZs), it is noteworthy 60 61 that the frequencies of OTUs affiliated to fermentative and exoelectrogenic bacteria, which are potent sources of in situ lactate, acetate, CO and/or H , are also highest along the sediment-depth of SSK42/6, as compared to 2 2 the other cores (Supplementary Tables 10–15). Additional factors such as sedimentation rate may also be respon- sible for the spatial variation in sulfate reduction rate; for instance, it is known that slower sedimentation rate may 62,63 enhance the oxygen exposure time of organic matter at the sediment water interface and in doing to preserve the lability of organic matter crucial for sulfate reduction . Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 6 www.nature.com/scientificreports/ 2− − Figure 5. Concentrations of pore-water SO (hollow circle), ΣHS (hollow square) and CH (filled triangle) 4 4 in sediment sub-samples in cores SSK42/1-8. e s Th ulfate-methane transition zone (SMTZ) in SSK42-6 is within 250 to 275 cmbsf. The maximum hydrogen sulfide concentration is observed within the SMTZ. 2− Coupled influence of organic matter remineralization and OSR on the SO concentration gradients of the 2− + SSK42 cores is evident from the down-depth decrease in TOC and SO , and increase in dissolved NH and 4 4 − 65 ΣHS concentrations (Fig. 5, Supplementary Fig. 2a, Supplementary Table 2). Across the eight cores studied, + − dissolved NH ranged between 59 and 2214 µ M, whereas ΣHS was detectable only along SSK42/5, 6, 7 and 8. Across these sediment cores, dissolved ΣHS concentrations were found to vary between non-detectable and 2− − 2010 µ M. Down-depth variations in the pore-water SO and ΣHS concentrations are also accompanied by 2− − 34 corresponding shifts in sulfur isotope ratios of SO and ΣHS . In SSK42/1 through 8, the core-top δ S 4 SO4 varies between 22.1 and 28.7‰ VCDT while the down-depth profiles of δ S show enrichment trends for all SO4 the eight cores (Fig. 6A). e Th δ S maxima for the individual sediment cores were found to vary between 30.2 SO4 34 − and 51‰ VCDT. δ S values recorded along SSK42/6, 7 and 8 vary between −22.4 and +32.5‰ VCDT and HS 34 2− 34 − exhibit down-depth S enrichment, as observed for the source SO (Fig. 6B). Whereas in SSK42/8, δ S 4 HS ranges from −29.4 to −21.4‰ VCDT with no definite enrichment trend. e Th down-depth enrichment trend of dissolved pore-water sulfate is attributed to the OSR by obligate anaerobes which preferentially utilize the lighter 32 34 66 sulfur isotope S, leaving pore fluids progressively enriched in S . In sediment cores studied for sulfur isotopic 2− − fractionation (i.e. SSK42/6, 7 and 8), Δ values range between 22 and 59.8 (Avg: 46.7 ± 10.4‰ VCDT) SO4 −HS and are well within the range of S-fractionations reported for OSR, with or without disproportionation of inter- 2− 67–70 mediate sulfur species (S°, S O ) . Notably, genus-level classification of the OTUs identified in the SSK42 2 3 cores showed the SRB diversities across the sample-sites were dominated by Desulfovibrio species that are known to be capable of sulfur isotopic fractionation up to 46‰ . Other sulfate-reducers prevalent in the SSK42 cores include Desulfococcus, Desulfobacteria and Desulfotalea spp., which can render S-isotope fractionation between 16.1 and 32.7‰ VCDT . In this context it is noteworthy that OTUs belonging to the anaerobic/facultatively anaerobic genera Arcobacter, Paracoccus, Sulfurimonas, i Th ohalomonas and Thiohalophilus that are capable of chemolithotrophically-oxidizing reduced sulfur species including sulfide , were detected intermittently along the sediment-depth of every coring station (Supplementary Tables 10–15). However, more information is required to ascertain whether they actually render oxidation of H S (or other intermediate sulfur species) in situ, and thereby influence net S-isotope fractionation in the pore-waters. 54,58 AOM is known to play a major role in shaping the sulfate concentration profile of sediments . AOM, 74,75 mediated by a syntrophic consortium of anaerobic methanotrophic archaea (ANME) and SRBs , takes place − − typically in the SMTZs of sediments, and results in the in situ enrichment of HCO and ΣHS . The carbon isotope ratio of dissolved inorganic carbon (DIC) in organic-matter–rich, anoxic marine sediments depends on 53,76 13 the mixing of DIC produced via OSR and AOM pathways . δ C of DIC produced via OSR approaches the 13 13 δ C of the in situ TOC, whereas δ C of AOM derived DIC is generally much lower than that of the in situ TOC. In the SSK42 cores, DIC concentrations increase with depth (DIC varies between 2.4 and 14.9 mM, across the Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 7 www.nature.com/scientificreports/ 34 2− − Figure 6. δ S profiles of dissolved sulfate (SO ) and hydrogen sulfide (ΣHS ) in SSK42 pore water 34 34 samples. δ S and δ S profiles show enrichment trends with depth. SO4 HS eight cores studied) (Fig. 7). It is noteworthy that below 60–90 cmbsf of SSK42/5, 6 and 7 (the cores underlying 13 13 the OMZ center), δ C gets significantly depleted (− 22 to −36.4‰ VPDB) in comparison to δ C (Fig. 7), DIC TOC 12 77,78 thereby indicating strong upward diffusion of C-enriched DIC produced via AOM at the SMTZ . On the 13 13 other hand, δ C remains much higher than δ C in the other sediment cores investigated, thereby suggest- DIC TOC ing the deeper occurrence of SMTZs in those areas. Shallow depths (250–275 cmbsf ) of SMTZs within the OMZ center (SK42/6 and 7) corroborate the intensification of AOM within this territory. Carbon isotopic composition of the methane detected in SSK42/6 and 7 (δ C values between −87.8 and −62.2‰ VCDT; Supplementary Table 2) is within the range of δ C values recorded previously for microbial methane . Furthermore, deep shotgun sequencing and analysis of whole metagenomes isolated from the sed- iment datapoints of SSK42/6 showed that the predominant methanogen in these sediments is Methanosarcina 13 80,81 (S. Bhattacharya et al., Unpublished), which produces methane having −108 to −76‰ δ C . Bioavailable organic carbon is the potential driver of intense biogeochemical activity in the sediments underlying the OMZ center. Characteristic features such as (i) enhanced TOC content (ii) 2− − tandem increase in J and SRB-diversity, (iii) shallowing of SMTZ, and (iv) buildup of high ΣHS , are all SO4 reflective of heightened carbon-sulfur cycling in the sediments underlying the center of the eastern Arabian Sea OMZ. In addition, increase in total bacterial diversity within this region implicates high organic matter input as the potential trigger for the intensification of the entire biogeochemical cascade. Collectively, these features sug- gest that high as well as low molecular weight carbon compounds, utilizable by complex organic matter degraders and sulfate-reducers respectively, are abundant under the OMZ center. Greater preservation of labile organic matter within the OMZ center has also been insinuated in previous studies that measured chlorophyll degrada- tion products in OMZ sediments across the eastern Arabian Sea . In this context it is further noteworthy that although TOC and SRB diversity reach their maxima within the zone of high J , water-depth dependent varia- SO4 tion profiles of the first two parameters do not exactly overlap with that of the third (Fig. 2b). We therefore suggest that the potential nature of the bioavailable organic-compounds is a key determinant of in situ sulfate reduc- 2− tion rate, over and above TOC content and SRB diversity. Corroborating this conclusion, J drops sharply SO4 beyond the OMZ center, concomitant to the rise in bottom-water oxygen, thereby pointing towards the progres- sive reduction in the burial of bioavailable organic matter (from the water columns) owing to oxygen influenced increase in benthic consumption, biodegradation and biotransformation of organic matter . In a nutshell, our study identifies the sediments impinged by the hypoxic waters of the OMZ center as major sinks of marine organic matter sequestration, and hotspots of carbon-sulfur cycling. The present findings hold implications for the predictive modeling of potential loci of hydrogen sulfide and methane expulsions, which can damage the pelagic and benthic biota severely. Solid phase analyses of these sediments might oer m ff ore insights 84,85 into the past variations in AOM intensity, SMTZ depth, and overall expanse of the OMZ . Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 8 www.nature.com/scientificreports/ Figure 7. Concentrations of pore-water DIC (filled black circle), and stable carbon isotopic ratios (δ C) of DIC (hollow black cicle) and TOC (filled blue triangle), down the sediment-depths of SSK42/1-8. Methods Coring and CTD profiling. Eight gravity cores were collected onboard RV Sindhu Sankalp (SSK-42) in December 2012 across the upper continental slope along the western continental margin of India (Supplementary Table 1 and Supplementary Fig. 1). Barring core SSK42/1 (water-depth: 1275 mbsl), all other cores are located within the perennial oxygen minimum zone. Sea water-depth, temperature, salinity and oxygen concentrations were measured using a SBE 911plus CTD profiler (Sea Bird Electronics, Bellevue, Washington, USA) in the vicin - ity of the coring locations. (Supplementary Table 1 and Supplementary Fig. 1). Sampling for geochemistry, and porosity measurement. Cores were sub-sampled on-board for hydrocarbon gas analysis and pore-water extraction. For hydrocarbon gas analyses, sediment was extracted by inserting 10 ml cut syringe deep into the core at an interval of 15 cm. Sediment was transferred to 28 ml glass vials containing 10 ml of sodium azide to arrest microbial activity. The vials were u fl shed with helium, crimped and stored at 4 °C until further analysis. For extraction of pore-water, 50 ml cut syringes were inserted deep into the core, at an interval of 15 cm. Sediment was transferred from the syringe into 50 ml centrifuge tubes (Tarsons Products Private Limited, Kolkata, India) under a stream of nitrogen to avoid oxidation of dissolved sulfide. The tubes were centrifuged at 4700 g and 4 °C for 15 minutes in a C30 cryocentrifuge (REMI Laboratory Instruments, Mumbai, India). The supernatant was filtered through 0.22 µ m Whatman syringe filter and the aliquots were transferred into glass vials (Supelco Inc, Bellefonte, PA, USA) for measurement of concentrations and isotope ratios of various constit- uents. Dissolved ΣHS was precipitated as CdS by the addition of CdNO . e Th vials were crimped immediately following nitrogen flushing and stored at 4 °C until further analysis. Sulfate concentration was measured in the ΣHS free supernatant pore-water. Sub-sampling of cores for solid phase analysis was carried out as 1 cm slabs. The sub samples were stored at 4 °C until further analysis. For porosity measurements, measured volume of sediment was dried at 105 °C. Moisture content in the sediment was calculated from the difference in wet and dry weight of sediments. Porosity was calculated as: (volume of sediment pore-water/wet sediment volume weight) × 100. Measurement of concentrations and isotope ratios. Pore-water sulfate concentrations were measured by a Basic IC plus 883 ion chromatograph (Metrohm AG, Herisau, Switzerland) equipped with a suppressed con- ductivity detector (Metrohm, IC detector 1.850.9010) and a MetrosepASupp 5 (150/4.0) anion exchange column (Metrohm AG, Herisau, Switzerland). A mixed solution of 1 mM NaHCO and 3.2 mM Na CO was used as the 3 2 3 eluent and 0.2 N H SO as the regeneration fluid. Pore-water samples were diluted 1000-fold with de-ionized 2 4 Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 9 www.nature.com/scientificreports/ water (Siemens, <0.06 µ S) prior to analysis. A 10 ppm Fluka mixed anion standard was used to prepare the calibration curve for quantification. The sample reproducibility was ±0.2 ppm. Dissolved sulfide concentrations were measured from the CdS precipitate following methylene blue complexation method. The absorbance of methylene blue complex was measured at 670 nm using a spectrophotometer (Spectrascan UV 2700, Chemito Instruments, Mumbai, India). Carbon coulometer (CM5130) was used to determine the concentration of dis- solved inorganic carbon. Ammonia was determined using the indophenol blue method. Absorbances were meas- ured at 640 nm using the spectrophotometer mentioned above. Concentration and carbon isotopic composition of hydrocarbon gases present in the head space samples was determined using a continuous flow isotope ratio mass spectrometer (Delta V Plus, Thermo Fisher Scientific, Germany) coupled with a gas chromatograph (Thermo Fisher Scientific, Germany). Peak area calibration curve was prepared using standard hydrocarbon gas mixtures for concentration measurement. Carbon isotope ratios are reported in standard format relative to the Vienna Peedee Belemnite (VPDB) standard. The external pre- 13 13 cision calculated for δ C is typically 0.07–0.09‰ (VPDB). δ C of DIC was analyzed on a isotope-ratio mass Plus spectrometer (Thermo Finnigan Delta XP) coupled to a Thermo Finnigan Gasbench II at the Stable Isotope Geochemistry Laboratory, University of Kentucky, USA. Pore-water sulfate was precipitated as BaSO by the addition of 1 ml of 1 M BaCl solution following acidification. Dissolved sulfide precipitated as CdS was filtered through a 0.22 µ m nitrocellulose filter paper. The precipitate was converted to silver sulfide (Ag S) by the addition of silver nitrate (AgNO ) solution. The dried and homogenized BaSO and Ag S precipitates were mixed with 3 4 2 V O and flash combusted at 1150 °C in an elemental analyzer (EA1112, Thermo Fisher Scientific, Germany) and 2 5 subsequently analysed for stable sulfur isotope ratios using a continuous flow isotope ratio mass spectrometer (Thermo Delta V plus, Thermo Fisher Scientific, Germany). All results are reported in standard delta notation (δ S) as per mil deviations from the VCDT (Vienna Canyon Diablo Troilite). IAEA standards S-1, S-2, S0-5 and S0-6 were used for instrument calibration. Sample reproducibility of ±0.3‰ is reported here. TIC content of standard reference material (Ultrapure CaCO from Sigma-Aldrich) was 12.0 ± 0.25%. Total carbon (TC) and nitrogen (TN) content of freeze dried and desalinated sediment samples were measured by EA1112 elemental analyzer (Thermo Fisher Scientific, Germany). Total organic Carbon (TOC) content was cal- culated by subtracting TIC from TC. 2,4-DNP was used as a calibration standard for TC. Reproducibility for TC in NIST-SRM1944 sediment standard was found to be 4.4 ± 0.2% and TN in MAG-1 sediment standard was found to be 0.24 ± 0.03%. δ C was measured on decarbonated sediments using a Delta-V-plus isotope ratio TOC mass spectrometer coupled with an EA1112 elemental analyzer (Thermo Fisher Scientific, Germany). The exter - nal standard reproducibility calculated for δ C using cellulose is typically −24.7 ± 0.1‰ (VPDB). TOC Depth-integrated Sulfate ux fl calculations. Diffusive sulfate flux (J ) was calculated from concentra- SO4 tion profiles using Fick’s first law assuming steady state conditions. JD =ϕ s(dC/dX) (4) SO4 −2 −1 Where J, C and ϕ represent the depth-integrated sulfate reduction rate (mmol cm yr ), sulfate concentration (mM) in the pore-water and average sediment porosity respectively. dC/dX is the sulfate concentration gradient 2 −1 and Ds (cm s ) is the molecular diffusivity corrected for tortuosity. Ds is calculated by the formula: Ds =+ Do/[1n(1 −ϕ)] (5) where n = 3 for clays and silt . Do is the sulfate diffusivity in the absence of particles. Since Do varies with tem- perature (water-depth), it was calculated for the varying bottom water temperatures using the empirical relation based on previous data . Sampling for studies of in situ microbial ecology, and isolation of metagenomic DNA from the sediment samples. For microbiological investigations, sediment samples were collected onboard, asepti- cally from various sediment-depths of SSK42/1, 3, 5, 6, 7 and 8, immediately after cutting and removal of lon- gitudinal halves of the PVC core-liners. Identities of the 87 sediment-datapoints investigated for microbiology across the six coring stations are given in Supplementary Tables 3–8. As mentioned above, the top one-cm of the exposed surface was first removed so as to avoid atmospheric and sea water contamination; sampling was carried out strictly under nitrogen flow. In order to sample a given sediment-depth from a given core (i.e. a particular sediment-datapoint), an approximately 5-mm-thick sediment-slice (spanning equally on either side of the core-height marking) was scooped out with sterile scalpels and put into a sterile polypropylene bottle. The head-space of every sample-containing bottle was flushed with nitrogen gas, following which it was sealed with Parafilm (Neenah, WI, USA) and immediately placed under refrigeration: samples meant for metagen- ome analysis and culture-based study (the latter is not a part of this paper) were stored at −20 °C and 4 °C, respectively. From the laboratory onboard SSK-42, en route to the laboratory at Bose Institute (Kolkata) these sample-preservation temperatures were maintained all-through. Upon reaching the laboratory, extraction of total community DNA (metagenome) from every sediment sample was completed within a week. Metagenomes were isolated using PowerSoil DNA Isolation Kit (MoBio, Carlsbad, CA, USA), as per manufacturer’s protocol. Taxonomic diversity estimation in the sedimentary microbial communities. Taxonomic com- position of the microbial communities occurring at various sediment-depths of SSK42/1, 3, 5, 6, 7 and 8 was determined by high-throughput sequencing and analyses of the V3 regions of all 16S rRNA genes present in the 88,89 respective metagenomes. Following methods described earlier , V3 regions of all Bacteria- or Archaea-specific 16S rRNA genes that were present in the metagenome of a given community were PCR-amplified separately using domain-specific oligonucleotide primer-pairs and sequenced by the Ion Torrent Personal Genome Machine (Ion Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 10 www.nature.com/scientificreports/ PGM; Thermo Fisher Scientific, Waltham, Massachusetts, USA) up to such a data-throughput level that ensured plateau in the rarefaction curve. Subsequently, all the V3 reads present in a given sequence dataset were clustered (at 97% sequence similarity level) into operational taxonomic units (OTUs) or potential microbial species-level entities (see below); taxonomic classification of these OTUs revealed the total diversity of the community. Amplification of 16S rRNA gene fragments and sequencing by Ion PGM. Amplification of 16S rRNA gene fragments and sequencing by Ion PGM was carried out using the Fusion Primer protocol, as described 88,89 previously . V3 regions of all potential bacterial or archaeal 16S rRNA genes present in the studied sample were amplified by polymerase chain reaction (PCR) using domain-specific oligonucleotide primers. In order to enable tandem sequencing of multiple samples on a single chip each sample DNA was subjected to PCR using a 16S forward primer prefixed with an Ion Torrent adapter and a unique sample-specific barcode or multiplex identi- fier in the following order in the 5′ to 3′ direction: (i) a 26-mer A-linker followed by a 4-mer A-linker key (bases represented in bold fonts), common for all sample primers, (ii) a 10-mer barcode unique to each sample primer followed by a common 3-mer barcode adaptor (all marked as stars), and finally (iii) the relevant domain-specific universal forward primer in its 5′ to 3′ direction (underlined bases). Reverse primers, in their turn, had (i) a com- mon trP1 adapter (bases represented in italics), followed by (ii) the relevant domain-specific universal reverse primer in its 5′ to 3′ direction (underlined bases). V3 regions of all bacterial 16S rRNA genes present in a metagenome were amplified using the forward primer 5′-CCA TCT CAT CCC TGC GTG TCT CCG ACT CAG *** *** *** *** *CC TAC GGG AGG CAG CAG-3′ (where the underlined portion represents the universal primer 341 f ) and the reverse primer 5′-CCT CTC TAT GGG CAG TCG GTG ATA TTA CCG CGG CTG CTG G-3′ (where the underlined portion represents the uni- versal primer 515r). V3 regions of all archaeal 16S rRNA genes present in a metagenome were amplified using the forward primer 5′-CCA TCT CAT CCC TGC GTG TCT CCG ACT CAG *** *** *** *AA TTG GAN TCA ACG CCG G-3′ (where the underlined portion represents the universal primer 344 f ) and the reverse primer 5′-CCT CTC TAT GGG CAG TCG GTG ATC GRC GGC CAT GCA CCW C-3′ (where the underlined portion represents the universal primer 522r). Each 50 µ l PCR reaction contained 10 µ l template (corresponding to ∼100 ng metagenomic DNA), 5 µ l 10 × KOD DNA polymerase buffer, 5 µl dNTP (0.25 mM each), 2 µl MgCl (25 mM), 1.5 µl (3%) DMSO, 3 µl each of the forward and reverse primer (0.3 µM each), 19.5 µl dH O and 1 µl KOD Hot Start DNA Polymerase enzyme (Novagen, Merck Biosciences, Darmstadt, Germany). PCR products were amplified for 30 cycles as follows: 95 °C for 20S, 65 °C for 30 s and 70 °C for 10 s. After amplification, all PCR products were electrophoresed on 2.5% −1 w/v agarose gel, purified by size selection, and adjusted to final concentrations of 10 ng μl using molecular grade water. PCR products from all the samples were pooled up at equal concentrations for subsequent Ion PGM sequencing. Before Ion PGM sequencing, size distribution and DNA concentration in the pooled-up amplicon mixture was examined using a Bioanalyzer 2100 (Agilent Technologies, Santa Clara, California, USA). The mixed sample was adjusted to a final concentration of 26 pM and attached to the surface of Ion Sphere Particles (ISPs) using an Ion Onetouch 200 Template kit (Thermo Fisher Scientific, Waltham, Massachusetts, USA) according to the manufacturer’s instructions. Manual enrichment of the resulting ISPs resulted in >95% templated-ISPs, which were then sequenced on Ion 316 Chips using the Ion PGM (Ion Express Template 200 chemistry) for 500 flows that gives an expected average read length of >220 bp. Post sequencing, and before retrieval from the sequencing machine, all reads were filtered by the inbuilt PGM software to remove low-quality- and polyclonal-sequences; sequences matching the PGM 3′ adaptor were also automatically trimmed. Every V3 sequence dataset ae ft r this quality-filtration was exported as a fastq file for downstream analyses; at the same time all were deposited to the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI) with individual run accession numbers (given in Supplementary Tables 3–8) under the BioProject accession number PRJNA309469. Clustering of V3 reads into OTUs. Every V3 sequence dataset retrieved from the Ion PGM was further filtered for quality value 20 and read-length threshold of 100 nucleotides. Selected reads were converted to fasta using Fastx_toolkit 0.0.13.2 (http://hannonlab.cshl.edu/fastx_toolkit/download.html). Subsequently preproc- essing was done using the “preproc” module of ESPRIT ; this involved the screening of only such reads that had intact 16S rDNA universal primer sequence plus the barcode adaptor and the barcode. Adaptor sequences were clipped off every individual read, subsequent to which all such reads were removed from the dataset that contained unidentified bases or were identical to another read, or were sheer subsets of a longer read. After preprocessing, denoising of reads was carried out as follows: (i) pairwise-alignment of reads was carried out following Needleman algorithm using the module “pairwise.seqs” available in MOTHUR ; (ii) the alignment 92 93 matrix created in this way was utilized for Single Linkage Preclustering using a script provided by the VAMPS (Visualization and Analysis of Microbial Population Structure project) integrated collection of tools http://vamps. mbl.edu/overview.php. Pairwise distance criterion of 0.02 (98% identity) was used to select reads in precluster- ing, which were then saved in fasta format. The purpose of this preclustering was to preempt the possibility of diversity-overestimation emanating from the potential 2% sequencing error of the Ion PGM platform. Hierarchical clustering of reads was carried out using various modules available in ESPRIT . During this pro- cess, reads were again aligned pairwise following Needleman algorithm, clusters or operational taxonomic units (OTUs) created with a pairwise distance criterion of 0.03 (97% identity), and finally the total number of OTUs or species-level entities determined. OTUs formed in this way were further used to draw statistical inferences. OTUs were filtered using a Perl script to remove all the singletons. Singletons were removed from* .Cluster as well as*. Cluster_List to make new_Cluster & new_Cluster_List files. The new Cluster files were run in the statistical mod- ule of ESPRIT to get ACE and Rarefaction analyses. The rarefaction information was further used in R package to plot for the number of reads taking part in OTU-formation against the number of OTUs formed. The clusters Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 11 www.nature.com/scientificreports/ (minus singletons) were used together with the fasta and frequency files created during the clustering operations to create consensus sequences, which in turn were taxonomically classified with the help of the “RDP Classifier” located at http://rdp.cme.msu.edu/classifier/classifier.jsp . Species-richness estimation for sulfate-reducing bacteria. To estimate the species-richness (or total OTU-count) of sulfate-reducing bacteria within a given sediment-community, the total number of OTUs ali ffi ated to such taxa were enumerated (from the RDP-classified OTU-set of that community), every member of which are known to reduce sulfate dissimilatorily. es Th e included the genera Desulfacinum, Desulfobacca, Desulfobaculum, Desulfocurvus, Desulfoglaeba, Desulfomonas, Desulfomonile, Desulforhabdus, Desulfosoma, Desulfovibrio, Desulfovirga, Desulfurella, Desulfuromonas, Desulfuromusa and Thermodesulforhabdus ; plus the entire fam- ilies Desulfarculaceae, Desulfobacteraceae, Desulfobulbaceae, Desulfohalobiaceae, Desulfomicrobiaceae and Desulfonatronaceae of Deltaproteobacteria . Notably, the OTU-set of none of the explored sediment-communities encompassed any consensus sequence ascribable to the SRB genera of the other bacterial phyla Chrysiogenetes, Firmicutes, Nitrospiraea, Synergistetes or e Th rmodesulfobacteria . Mathematical analysis of SRB diversity distribution. Probability density functions of the trends of variation observed in the distribution of SRB-OTUs along each sediment core were determined via the following steps. (I) First, the SRB-OTU-count profile of the core was segmented into apparent zones of u fl ctuation along the sediment-surface to core-bottom trajectory. This was done by referring to the relevant histogram within column I of Fig. 3 (each of these histograms was derived based on the normalized values of the SRB-OTU-counts for the various sediment-datapoints of the core; normalization of the SRB-OTU-count for a given datapoint was carried out by taking its ratio with the highest SRB-OTU-count obtained within that core). (II) In the second step, inde- pendent prediction of an approximate probability density function was done for every fluctuation-zone identified. (III) Finally, optimization of every probability density function was done with reference to the histogram through 2 2 96,97 minimization of the χ value. Levenberg-Marquardt algorithm was used for χ minimization and the work was done using the software OriginPro 9 in a computer having an Intel(R) Core(TM) i5-2450 64-bit CPU and 8 GB RAM. In order to ensure that the theoretical density functions were best fitted, up to 4000 independent iterations were performed on every dataset. For each calculation, tolerance level of the reduced χ value was −9 considered at 10 . Data availability. All data supporting the findings of this paper are available within the Article and Supplementary Information files. DNA sequence data pertaining to the microbial diversity study are available from the Sequence Read Archive (SRA) of the National Center for Biotechnology Information (NCBI), a part of the US National Library of Medicine (NLM), a branch of the National Institutes of Health, USA: relevant acces- sion numbers are all listed in Supplementary Tables 3–8 within the Supplementary Information files. References 1. Lam, P. & Kuypers, M. M. M. Microbial nitrogen cycling processes in oxygen minimum zones. Ann. Rev. Mar. Sci. 3, 317–345 (2011). 2. Stramma, L., Johnson, G. C., Firing, E. & Schmidtko, S. Eastern Pacific oxygen minimum zones: Supply paths and multidecadal changes. J. Geophys. Res. 115, 1–12 (2010). 3. Wright, J. J., Konwar, K. M. & Hallam, S. J. Microbial ecology of expanding oxygen minimum zones. Nat. Rev. Microbiol. 10, 381–394 (2012). 4. Mot ffi t, S. E. et al . Response of seafloor ecosystems to abrupt global climate change. Proc. Natl. Acad. Sci. USA 112, 4684–4689 (2015). 5. Naqvi, W. A. Geographical extent of denitrification in the Arabian Sea in relation to some physical processes. Oceanol. Acta 14, 281–290 (1991). 6. Ulloa, O., Canfield, D. E., De Long, E. F., Letelier, R. M. & Stewart, F. J. Microbial oceanography of anoxic oxygen minimum zones. Proc. Natl. Acad. Sci. USA 109, 15996–16003 (2012). 7. Ward, B. B. et al. Denitrification as the dominant nitrogen loss process in the Arabian Sea. Nature 461, 78–81 (2009). 8. Acharya, S. S. & Panigrahi, M. K. Eastward shift and maintenance of Arabian Sea oxygen minimum zone: understanding the paradox. Deep Sea Res. Part I Oceanogr. Res. Pap. 115, 240–252 (2016). 9. Bange, H. W., Naqvi, S. W. A. & Codispoti, L. A. The nitrogen cycle in the Arabian Sea. Prog. Oceanogr. 65, 145–158 (2005). 10. Naqvi, S. W. A. et al. Marine hypoxia/anoxia as a source of CH and N O. Biogeosciences 7, 2159–2190 (2010). 4 2 11. McCreary, J. P. et al. Dynamics of the Indian–Ocean oxygen minimum zones. Prog. Oceanogr. 112–113, 15–37 (2013). 12. Madhupratap, M. et al. Mechanism of the biological response to winter cooling in the northeastern Arabian Sea. Nature 384, 549–552 (1996). 13. Naqvi, S. W. A. et al. Coastal versus open–ocean denitrification in the Arabian Sea. Biogeosciences 3, 621–633 (2006). 14. Banse, K., Naqvi, S. W. A., Narvekar, P. V., Postel, J. R. & Jayakumar, D. A. Oxygen minimum zone of the open Arabian Sea: variability of oxygen and nitrite from daily to decadal timescales. Biogeosciences 11, 2237–2261 (2014). 15. Helly, J. J. & Levin, L. A. Global distribution of naturally occurring marine hypoxia on continental margins. Deep Sea Res. Part I Oceanogr. Res. Pap. 51, 1159–1168 (2004). 16. Cavan, E. L., Trimmer, M., Shelley, F. & Sanders, R. Remineralization of particulate organic carbon in an ocean oxygen minimum zone. Nat. Commun. 8, 14847 (2017). 17. Van der Weijden, C. H., Reichart, G. J. & Visser, H. J. Enhanced preservation of organic matter in sediments deposited within the oxygen minimum zone in the northeastern Arabian Sea. Deep Sea Res. Part I Oceanogr. Res. Pap. 46, 807–830 (1999). 18. Schulte, S., Mangelsdorf, K. & Rullkötter, J. Organic matter preservation on the Pakistan continental margin as revealed by biomarker geochemistry. Org. Geochem. 31, 1005–1022 (2000). 19. Seiter, K., Hensen, C. & Zabel, M. Benthic carbon mineralization on a global scale. Global Biogeochem. Cycles 19, 1–26 (2005). 20. Bowles, M. W., Mogollón, J. M., Kasten, S., Zabel, M. & Hinrichs, K.-U. Global rates of marine sulfate reduction and implications for sub–sea–floor metabolic activities. Science 344, 889–891 (2014). 21. Cowie, G. et al. Comparative organic geochemistry of Indian margin (Arabian Sea) sediments: estuary to continental slope. Biogeosciences 11, 6683–6696 (2014). 22. Jessen, G. L. et al. Hypoxia causes preservation of labile organic matter and changes seafloor microbial community composition (Black Sea). Sci. Adv. 3, e1601897 (2017). Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 12 www.nature.com/scientificreports/ 23. Luckge, A., Ercegovac, M., Strauss, H. & Littke, R. Early diagenetic alteration of organic matter by sulfate reduction in Quaternary sediments from the northeastern Arabian Sea. Mar. Geol. 158, 1–13 (1999). 24. Schmaljohann, R. et al. Oxygen–minimum zone sediments in the northeastern Arabian Sea o P ff akistan: a habitat for the bacterium Thioploca. Mar. Ecol. Prog. Ser. 211, 27–42 (2001). 25. Jørgensen, B. B. & Nelson, D. C. Sulfide oxidation in marine sediments: geochemistry meets microbiology in Sulfur Biogeochemistry – Past and Present (ed. Amend, J. P., Edwards, K. J. & Lyons, T. W.) 379, 63–81 (The Geological Society of America Inc., 2004). 26. Chronopoulou, P. M., Shelley, F., Pritchard, W. J., Maanoja, S. T. & Trimmer, M. Origin and fate of methane in the Eastern Tropical North Pacific oxygen minimum zone. ISME J. 11(6), 1386–1399 (2017). 27. Fossing, H., Ferdelman, T. G. & Berg, P. Sulfate reduction and methane oxidation in continental margin sediments influenced by irrigation (South–East Atlantic off Namibia). Geochim. Cosmochim. Acta 64, 897–910 (2000). 28. Kuwabara, J. S., van Geen, A., McCorkle, D. C. & Bernhard, J. M. Dissolved sulfide distributions in the water column and sediment pore-waters of the Santa Barbara basin. Geochim. Cosmochim. Acta 63, 2199–2209 (1999). 29. Schubert, C. J., Ferdelman, T. G. & Strotmann, B. Organic matter composition and sulfate reduction rates in sediments off Chile. Org. Geochem. 31, 351–361 (2000). 30. Schulz, H. N. & Schulz, H. D. Large sulfur bacteria and the formation of phosphorite. Science 307, 416–418 (2005). 31. Treude, T. et al. Anaerobic oxidation of methane and sulfate reduction along the Chilean continental margin. Geochim. Cosmochim. Acta 69, 2767–2779 (2005). 32. Dale, A. W., Graco, M. & Wallmann, K. Strong and dynamic benthic-pelagic coupling and feedbacks in a coastal upwelling system (Peruvian shelf ). Front. Mar. Sci. 4, 29 (2017). 33. Passier, H. F., Luther, G. W. & de Lange, G. J. Early diagenesis and sulphur speciation in sediments of the Oman Margin, northwestern Arabian Sea. Deep Sea Res. Part II Top. Stud. Oceanogr. 44(6), 1361–1380 (1997). 34. Schenau, S. J., Passier, H. F., Reichart, G. J. & De Lange, G. J. Sedimentary pyrite formation in the Arabian Sea. Mar. Geol. 185(3), 393–402 (2002). 35. Lückge, A., Horsfield, B., Littke, R. & Scheeder, G. Organic matter preservation and sulfur uptake in sediments from the continental margin off Pakistan. Org. Geochem. 33(4), 477–488 (2002). 36. Yoshinaga, M. Y. et al. Possible roles of uncultured archaea in carbon cycling in methane–seep sediments. Geochim. Cosmochim. Acta 164, 35–52 (2015). 37. Orsi, W. D. et al. Climate oscillations reflected within the microbiome of Arabian Sea sediments. Sci. Rep. 7(1), 6040 (2017). 38. Divya, B., Soumya, K. V. & Nair, S. 16SrRNA and enzymatic diversity of culturable bacteria from the sediments of oxygen minimum zone in the Arabian Sea. Antonie Leeuwenhoek 98, 9–18 (2010). 39. Weeks, S. J., Currie, B., Bakun, A. & Peard, K. R. Hydrogen sulphide eruptions in the Atlantic Ocean off southernAfrica: implications of a new view based on SeaWiFS satellite imagery. Deep Sea Res. Part I Oceanogr. Res. Pap. 51(2), 153–172 (2004). 40. Lavik, G. et al. Detoxification of sulphidic African shelf waters by blooming chemolithotrophs. Nature 457, 581–584 (2009). 41. Schunck, H. et al. Giant hydrogen sulfide plume in the oxygen minimum zone off Peru supports chemolithoautotrophy. PLoS One 8(8), e68661 (2013). 42. Tyson, R. V. Sedimentary Organic Matter (Chapman and Hall, New York, 615 (1995) 43. Vandewiele, S., Cowie, G., Soetaert, K. & Middelburg, J. J. Amino acid biogeochemistry and organic matter degradation state across the Pakistan margin oxygen minimum zone. Deep Sea Res. Part II Top. Stud. Oceanogr. 56, 376–392 (2009). 44. Resplandy, L. et al. Controlling factors of the oxygen balance in the Arabian Sea’s OMZ. Biogeosciences 9, 5095–5109 (2012). 45. Levin, L. A. Oxygen minimum zone benthos: adaptation and community response to hypoxia. Oceanogr. Mar. Biol. 41, 1–45 (2003). 46. Ingole, B. S., Sautya, S., Sivadas, S., Singh, R. & Nanajkar, M. Macrofaunal community structure in the western Indian continental margin including the oxygen minimum zone. Mar. Ecol. 31, 148–166 (2010). 47. Hedges, J. I. et al. Evidence for non–selective preservation of organic matter in sinking marine particles. Nature 409, 801–804 (2001). 48. Keil, R. G., Neibauer, J. A., Biladeau, C., van der Elst, K. & Devol, A. H. A multiproxy approach to understanding the “enhanced” flux of organic matter through the oxygen–deficient waters of the Arabian Sea. Biogeosciences 13, 2077–2092 (2016). 49. Koho, K. A. et al. Microbial bioavailability regulates organic matter preservation in marine sediments. Biogeosciences 10, 1131–1141 (2013). 50. Cottrell, M. T. & Kirchman, D. L. Natural assemblages of marine proteobacteria and members of the Cytophaga-Flavobacter cluster consuming low–and high–molecular–weight dissolved organic matter. Appl. Environ. Microbiol. 66, 1692–1697 (2000). 51. Sezenna M. L. Proteobacteria: Phylogeny, Metabolic Diversity and Ecological Ee ff cts. Nova Science Publishers Inc.: Hauppauge NY, USA. (2011). 52. Liu, W. et al. Microbial network for waste activated sludge cascade utilization in an integrated system of microbial electrolysis and anaerobic fermentation. Biotechnol. biofuels 9, 83 (2016). 53. Komada, T. et al. Organic matter cycling across the sulfate-methane transition zone of the Santa Barbara Basin, California Borderland. Geochim. Cosmochim. Acta 176, 259–278 (2016). 54. Jørgensen, B. B., Weber, A. & Zop, J fi . Sulfate reduction and anaerobic methane oxidation in Black Sea sediments. Deep Sea Res. Part I Oceanogr. Res. Pap. 48, 2097–2120 (2001). 55. Sigalevich, P., Meshorer, E., Helman, Y. & Cohen, Y. Transition from anaerobic to aerobic growth conditions for the sulfate-reducing bacterium Desulfovibrio oxyclinae results in flocculation. Appl. Environ. Microbiol. 66, 5005–5012 (2000a). 56. Sigalevich, P., Baev, M. V., Teske, A. & Cohen, Y. Sulfate reduction and possible aerobic metabolism of the sulfate-reducing bacterium Desulfovibrio oxyclinae in a chemostat coculture with Marinobacter sp. strain MB under exposure to increasing oxygen concentrations. Appl. Environ. Microbiol. 66, 5013–5018 (2000b). 57. Maltby, J., Sommer, S., Dale, A. W. & Treude, T. Microbial methanogenesis in the sulfate–reducing zone of surface sediments traversing the Peruvian margin. Biogeosciences 13, 283–299 (2016). 58. Jørgensen, B. B. & Kasten. S. Sulfur cycling and methane oxidation in Marine geochemistry 2nd ed. (ed. Schulz, H. D. & Zabel, M.) 271–309 (Springer–Verlag Berlin Heidelberg, 2006). 59. Drake H. L., Küsel K. & Matthies C. Acetogenic Prokaryotes in e Th Prokaryotes Vol. 2. (ed. Dworkin, M., Falkow, S., Rosenberg, E., Schleifer, K.-H. & Stackebrandt, E.). 354–420 (Springer, New York, 2006). 60. Tamang, J. P., Watanabe, K. & Holzapfel, W. H. Review: Diversity of microorganisms in global fermented foods and beverages. Front. Microbiol. 7, 377 (2016). 61. Logan, B. E. Exoelectrogenic bacteria that power microbial fuel cells. Nat. Rev. Microbiol. 7(5), 375–381 (2009). 62. Hedges, J. & Keil, R. Sedimentary organic matter preservation: an assessment and speculative synthesis. Mar. Chem. 49, 81–115 (1995). 63. Hartnett, H. E., Keil, R. G., Hedges, J. I. & Devol, A. H. Influence of oxygen exposure time on organic carbon preservation in continental margin sediments. Nature 391, 572–575 (1998). 64. Canfield, D. E. Sulfate reduction and oxic respiration in marine sediments: implications from organic carbon preservation in euxinic environments. Deep Sea Res. 36(1), 121–138 (1989). 65. Froehlich, P. N. et al. Early oxidation of organic matter in pelagic sediments of the eastern equatorial Atlantic: suboxic diagenesis. Geochim. Cosmochim. Acta 43, 1075–1090 (1979). 66. Mazumdar, A., Peketi, A., Joao, H., Dewangan, P., Borole, D. V. & Kocherla, M. Sulfidization in a shallow coastal depositional setting: diagenetic and paleoclimatic implications. Chem. Geol. 322-323, 68–78 (2012). Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 13 www.nature.com/scientificreports/ 67. Habicht, K. S. & Canfield, D. E. Sulfur isotope fractionation during bacterial sulfate reduction in organic-rich sediments. Geochim. Cosmochim. Acta 61, 5351–5361 (1997). 68. Canfield, D. E. Biogeochemistry of sulfur isotopes. Rev. Mineral. Geochem. 43, 607–636 (2001). 69. Habicht, K. S. & Canfield, D. E. Isotope fractionation by sulfate–reducing natural populations and the isotopic composition of sulfide in marine sediments. Geology 29, 55–558 (2001). 70. Sim, M. S., Bosak, T. & Ono, S. Large sulfur isotope fractionation does not require disproportionation. Science 333, 74–77 (2011). 71. Kaplan, I. R. & Rittenberg, S. C. Microbiological fractionation of sulphur isotopes. Microbiology 34(2), 195–212 (1964). 72. Detmers, J., Bruchert, V., Habicht, K. S. & Kuever, J. Diversity of sulfur isotope fractionations by sulfate–reducing prokaryotes. Appl. Env. Microbiol. 67(2), 888–894 (2001). 73. Ghosh, W. & Dam, B. Biochemistry and molecular biology of lithotrophic sulfur-oxidation by taxonomically and ecologically diverse bacteria and archaea. FEMS Microbiol. Rev. 33, 999–1043 (2009). 74. Orphan, V. J., House, C. H., Hinrichs, K.-U., McKeegan, K. D. & DeLong, E. F. Methane– consuming archaea revealed by directly coupled isotopic and phylogenetic analysis. Science 293, 484–487 (2001). 75. Knittel, K. & Boetius, A. Anaerobic oxidation of methane: progress with an unknown process. Ann. Rev. Microbol. 63, 311–334 (2009). 76. Holler, T. et al. Carbon and sulfur back flux during anaerobic microbial oxidation of methane and coupled sulfate reduction. Proc. Natl. Acad. Sci. USA 108, 1484–1490 (2011). 77. Yoshinaga, M. Y. et al. Carbon isotope equilibration during sulphate-limited anaerobic oxidation of methane. Nat. Geosci. 7, 190–194 (2014). 78. Pohlman, J. W. et al. Anaerobic methane oxidation in low–organic content methane seep sediments. Geochim. Cosmochim. Acta 108, 184–201 (2013). 79. Whiticar, M. J. Carbon and hydrogen isotope systematics of bacterial formation and oxidation of methane. Chem. Geol. 161, 291–314 (1999). 80. Londry, K. L. & des Marais, D. J. Stable carbon isotope fractionation by sulfate–reducing bacteria. Appl. Environ. Microbiol. 69, 2942–2949 (2003). st 81. Canfield, D. E., Kristensen, E & Thamdrup, B. Aquatic Geomicrobiology, 1 Ed., Adv. in Mar. Biol. 48, 640 (Elsevier Academic Press, 2005). 82. Rasiq, K. T., Kurian, S., Karapurkar, S. G. & Naqvi, S. W. A. Sedimentary pigments and nature of organic matter within the oxygen minimum zone (OMZ) of the Eastern Arabian Sea (Indian margin). Estuar. Coast. Shelf Sci. 176, 91–101 (2016). 83. Levin, L. A., Whitcraft, C. R., Mendoza, G. F., Gonzalez, J. P. & Cowie, G. Oxygen and organic matter thresholds for benthic faunal activity on the Pakistan margin oxygen minimum zone (700–1100m). Deep Sea Res. Part II Top. Stud. Oceanogr. 56, 449–471 (2009). 84. Peketi, A. et al. Coupled C–S–Fe geochemistry in a rapidly accumulating marine sedimentary system: diagenetic and depositional implications. Geochem. Geophys. Geosyst. 16(9), 2865–2883 (2015). 85. Lin, Q., Wang, J., Algeo, T. J., Sun, F. & Lin, R. Enhanced framboidal pyrite formation related to anaerobic oxidation of methane in the sulfate–methane transition zone of the northern South China Sea. Mar. Geol. 379, 100–108 (2016). 86. Iversen, N. & Jørgensen, B. B. Diffusion coefficients of sulfate and methane in marine sediments: influence of porosity. Geochim. Cosmochim. Acta 57, 571–578 (1993). 87. Li, Y.-H. & Gregory, S. Diffusion of ions in sea water and in deep sea sediments. Geochim. Cosmochim. Acta 38, 703–714 (1974). 88. Ghosh, W. et al. Resilience and receptivity worked in tandem to sustain a geothermal mat community amidst erratic environmental conditions. Sci. Rep. 5, 12179 (2015). 89. Roy, C. et al. Global association between thermophilicity and vancomycin susceptibility in bacteria. Front. Microbiol. 7, 412 (2016). 90. Sun, Y. et al. ESPRIT: estimating species richness using large collections of 16S rRNA pyrosequences. Nucleic Acids Res. 37, e76 (2009). 91. Schloss, P. D. et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–41 (2009). 92. Junemann, S. et al. Bacterial community shift in treated periodontitis patients revealed by ion torrent 16S rRNA gene amplicon sequencing. PLoS One 7, e41606 (2012). 93. Huse, S. M., Welch, D. M., Morrison, H. G. & Sogin, M. L. Ironing out the wrinkles in the rare biosphere through improved OTU clustering. Environ. Microbiol. 12, 1889–1898 (2010). 94. Dixon, P. VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930 (2003). 95. Rabus, R., Hansen, T. A. & Widdel, F. Dissimilatory sulfate–and sulfur–reducing prokaryotes in The Prokaryotes Vol. 2 (eds Dworkin, M. et al.) 659–768 (Springer, 2006). 96. Marquardt, D. W. An Algorithm for Least–Squares Estimation of Nonlinear Parameters. SIAM J. Appl. Math. 11, 431–441 (1963). 97. Moré, J. J. The Levenberg-Marquardt algorithm: implementation and theory in Numerical Analysis (ed. Watson, G. A.) 105-116 (Springer, 1978). Acknowledgements We thank the Director CSIR-National Institute of Oceanography for supporting this study. Funding from Ministry of Earth Sciences, Government of India (GoI) (GAP2303) for the acquisition of sediment cores and data is sincerely acknowledged. We acknowledge the financial support (MoES/36/00IS/Extra/19/2013) given by the Earth System Science Organization, Ministry of Earth Sciences, GoI, for conducting the microbiological studies. Sincere thanks are due to ship cell and crew members for on board support. S.F. and R.R. were endowed with fellowships by the Council of Scientific and Industrial Research, GoI. S.B. and C.R. received fellowship from the Bose Institute and University Grants Comision, GoI, respectively. Sincere thanks to Dr. P. Babu for providing coulometry facility. We are also extremely thankful to our friends Mr. Pravin Nilawe, Dr. Neeraj Kumar Chauhan, Dr. Atima Agarwal and Mr. Ashish Kumar George of Invitrogen BioServices India Pvt Ltd, Thermo Fisher Scientific, for constantly tutoring us on current metagenomic methodologies. Comments from Dr. D. Shenoy (CSIR-NIO) and Dr. Siby Kurian (CSIR-NIO) have benefited the manuscript. Author Contributions A.M. and W.G. designed the study. A.M., W.G. and S.F. compiled and interpreted all data, and wrote the paper. S.B. and A.P. provided inputs in analyses and interpretations. S.F., A.P., M.A.C., P.M., R.D.S. and P.L.S.R. performed the geochemical experiments. S.B. R.R. and C.R. performed the microbiological experiments. T.M. and S.B. , , carried out statistical analysis of microbiological data. All authors reviewed the manuscript. Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 14 www.nature.com/scientificreports/ Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-27002-2. Competing Interests: The authors declare no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2018 Scientific REPO R TS | (2018) 8:8665 | DOI:10.1038/s41598-018-27002-2 15
Scientific Reports – Springer Journals
Published: Jun 6, 2018
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.
Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.
All the latest content is available, no embargo periods.
“Hi guys, I cannot tell you how much I love this resource. Incredible. I really believe you've hit the nail on the head with this site in regards to solving the research-purchase issue.”Daniel C.
“Whoa! It’s like Spotify but for academic articles.”@Phil_Robichaud
“I must say, @deepdyve is a fabulous solution to the independent researcher's problem of #access to #information.”@deepthiw
“My last article couldn't be possible without the platform @deepdyve that makes journal papers cheaper.”@JoseServera