Geology and land use shape nitrogen and sulfur cycling groundwater microbial communities in Pacific Island aquifers

ISME Communications volume 3, Article number: 58 (2023) Cite this article

Metrics details

Resource-constrained island populations have thrived in Hawai’i for over a millennium, but now face aggressive new challenges to fundamental resources, including the security and sustainability of water resources. Characterizing the microbial community in groundwater ecosystems is a powerful approach to infer changes from human impacts due to land management in hydrogeological complex aquifers. In this study, we investigate how geology and land management influence geochemistry, microbial diversity and metabolic functions. We sampled a total of 19 wells over 2-years across the Hualālai watershed of Kona, Hawai’i analyzing geochemistry, and microbial communities by 16S rRNA amplicon sequencing. Geochemical analysis revealed significantly higher sulfate along the northwest volcanic rift zone, and high nitrogen (N) correlated with high on-site sewage disposal systems (OSDS) density. A total of 12,973 Amplicon Sequence Variants (ASV) were identified in 220 samples, including 865 ASVs classified as putative N and sulfur (S) cyclers. The N and S cyclers were dominated by a putative S-oxidizer coupled to complete denitrification (Acinetobacter), significantly enriched up to 4-times comparatively amongst samples grouped by geochemistry. The significant presence of Acinetobacter infers the bioremediation potential of volcanic groundwater for microbial-driven coupled S-oxidation and denitrification providing an ecosystem service for island populations dependent upon groundwater aquifers.

Microbes play a regulatory role in the convoluted interactions between hydrology, geology, and land use [1, 2] resulting in highly variable groundwater biogeochemistry. Microbes in subsurface groundwater aquifers represent roughly 40% of Earth's microbial life [3, 4], and are crucial for ecosystem services including providing clean drinking water. Subsurface microbial communities improve water quality by breaking down organic material and cycling nitrogen (N), sulfur (S), and iron (Fe) [5]. Unlike the redox structuring of sedimentary systems, opposing redox reactions co-occur in unique microenvironments (e.g., biofilms, particles) resulting in a diverse groundwater microbiome in large subsurface aquatic habitats [6, 7]. As a result, the use of inorganic compounds as energy sources (chemolithoautotrophy) are very common metabolic pathways in groundwater ecosystems, often linked to carbon (C), hydrogen (H), and sulfur (S) cycles with oxygen(O) and nitrate (NO3−) as the most widely used electron acceptors [1]. In one of the few microbiome studies conducted in deep Hawaiian groundwater aquifers, analysis revealed high functional diversity characterized primarily by chemolithotrophic metabolisms [8].

Characterizing the microbial community of groundwater ecosystems is a powerful approach to infer geochemical changes resulting from human impacts and land use in hydrogeologically complex aquifers [9, 10]. Changes in subsurface microbial community structure and diversity can alert users to decreases in water quality, because microbes in the subsurface react rapidly to alterations in their environment due to land use changes and/or contaminants, resulting in decreased diversity and shortened biogeochemical pathways [11, 12]. This is particularly important in Hawai’i as it depends almost entirely on groundwater for all water resources. Approximately 89% of Hawai’i's potable water supply is located in subsurface water resources, with minimal contributions from surface water (i.e., reservoirs, rain catchment) [13], and small-scale desalination activities. Hawai’i has the greatest number of cesspools per capita in the United States [14]. Non-point source pollution from personal on-site sewage disposal systems (OSDS, cesspools) and agriculture are the two major sources of human-derived nutrients that impact coastal ecosystems and have been recognized as major environmental problems in Hawai’i [15,16,17,18].

Nitrogen is one of the most common nonpoint source contaminants in groundwater [12, 19, 20]. Inputs of N are important determinants of major ion groundwater geochemistry [21], and subsurface microbial communities are responsible for regulating the effects of anthropogenic N and its transition into surface water or coastal ecosystems [22,23,24]. Autotrophic denitrification is the dominant N-removal process in oligotrophic groundwater where inorganic compounds (Fe, S, manganese, Mn; hydrogen gas, H2) are oxidized by reducing nitrite (NO2−) or nitrate (NO3−) to gaseous nitrogen (N2) [25]. [23, 26]. In addition, microbial-driven denitrification can occur both anaerobically, which was thought to be the only process until the 1980's, and aerobically, first described in a sulfide-oxidizing wastewater treatment plant [27, 28]. Human inputs of N into groundwater and its detrimental effects on coastal ecosystems has been heavily studied in Hawai’i [29,30,31], however the effect of excess N in volcanic, deep, freshwater aquifers and their subsurface microbial communities has yet to be characterized.

In addition to N, the presence of sulfate drives community structure in groundwater systems, and often points to a link between substrates and availability of electron donors [5, 32]. The presence of S species such as sulfate (SO42−), sulfides (S2−), polysulfides (SxO62−), thiosulfate (S2O32−), and sulfite (SO32−) are indicative of inputs from geothermal activity or seawater in groundwater systems [33,34,35]. The presence of SO42− in groundwater is often attributed to products of abiotic reactions, however, both anaerobic and aerobic biological sulfide oxidation is thermodynamically favorable in most environments [36]. Further, the relative availability of reduced S, C and Fe are also key determinants of nitrate removal pathways [26]. When S is present in carbon-limited systems, denitrification is powered by microbial mediated oxidation of reduced sulfur to sulfate [23], and occurs much faster than abiotic production of sulfate.

The Hualālai watershed (Kona, Hawai’i) has no surface runoff or drainage [37], is a semi-arid climate and is acutely vulnerable to stress from urban development due to the complete dependence on groundwater for municipal and agricultural water use [38]. Further, economic demands to support development in the Keauhou aquifer seek to exploit these aquifers and put undue strain on the quality of groundwater supplies. The watershed has the potential for large inputs of N due to personal OSDS, as well as S inputs from salt-water intrusion, aerosol deposition, and geothermal activity. In addition, fluctuating rainfall patterns reveal an overall drying trend, and rising temperatures is predicted to result in decreasing recharge and groundwater storage threatening the future of secure and safe freshwater supplies [39, 40].

In this study, we investigated subsurface microbial community structure and function with regard to N and S cycling in deep, volcanic Hawaiian island aquifers. We explored 1) how geology and land management practices in the Hualālai watershed influenced aquifer geochemistry, 2) the factors driving variability in microbial communities including geothermal inputs (S) along the rift zone and OSDS (N) inputs in the southern Keauhou aquifer, and 3) the microbial functional capacity for N and S cycling and how it varies with geochemical spatial pattern in these aquifers. The large and complex groundwater volcanic aquifers in Hawaii provide an opportunity to characterize the N and S cycling functions of microbial communities across different geologies and land management practices. This research has implications for land and water managers that must address decreased water quality in these drinking water resources affecting public health and future development in this watershed.

Hawai’i Island is the largest (10,464 km2) and youngest island in the Hawaiian archipelago [41]. The island was formed by five shield volcanoes with overlapping spatial and temporal lava flows, and geothermally active rift zones, which created complex hydrogeological aquifer connections [42, 43]. Hualālai, with a peak elevation of 2500 m above mean sea level, is the third most active volcano and consists of the Kīholo and Keauhou aquifer systems [44, 45]. Basalt flows and cinder cones build Hualālai and lie primarily along three geothermally-active rift zones trending northwest, northeast and southeast from the volcano summit (Fig. 1) [46]. The Hualālai watershed is located on the leeward (dry) side of Hawai’i Island with annual rainfall along the Kona coast ranging from 204 to 750 mm, however, the Kona rain belt upslope receives higher annual rainfall between 750 and 1350 mm [47].

Sample sites are colored by groups (N1, S, N2, and N3), and the Hualālai rift zones are named by direction (North, Northwest, and Southeast). Agriculture lands are shaded light blue, and OSDS sites are indicated by black circles.

Groundwater samples were collected from 19 private and publicly owned vertical or inclined shaft production wells in quarterly intervals from August 2017 to March 2019 (Fig. 1; Supplemental Information (SI)). Raw groundwater samples were collected prior to chlorination after sample rinsing 10-L cubitainers 3 times at well pump stations following a 15–20 min well purge [48] and stored at 4 °C until arrival in the field laboratory. Groundwater was subsampled into acid-washed 250-mL HDPE (Nalgene) bottles and stored (15–30 days) at 4 °C until geochemical analysis. Triplicate 2-L samples were filtered first through a 47-mm diameter 0.8 µm (GH Polypro, Pall Gelman Inc., MI) and then through a 0.2 µm pore-sized hydrophilic polypropylene membrane (Pall Gelman Inc., MI) filter and both were frozen at −20 °C prior to DNA extraction.

Groundwater was separated into acid-washed 250 mL HDPE (Nalgene) bottles and stored (15–30 days) at 4 °C until nutrient analysis. Standard methods for nutrient analysis were followed according to manufacture protocols for a four channel autoanalyzer (Astoria-Pacific International, Astoria Pacific, Clackamas, OR). Total N and P was analyzed as PO43− and NOx− following alkaline persulfate digestion [49]. Nutrient concentrations were determined for NO3− + NO2− = NOx−; orthophosphate, (PO43−); silicate (Si); ammonium (NH4+); total nitrogen, (TN); and total phosphorus, (TP). A YSI Pro Plus multiparameter water quality meter was used to measure groundwater temperature (°C), pH, dissolved oxygen (DO, mgL−1), and specific conductance (SPC, µS cm−1) after pump flushing and prior to the collection of water for DNA analysis (YSI, Inc., Yellow Springs, OH). Samples for major ions (SO42−; magnesium, Mg2+; chloride, Cl−) and select trace metals (chromium, Cr; iron, Fe; and manganese, Mn) were filtered through a 47-mm, 0.2 µm pore-sized hydrophilic polypropylene membrane (Pall Gelman Inc., Ann Arbor, MI) in a reusable filter holder and receiver (Thermo Fisher Scientific, Waltham, MA) and separated into 60-mL HDPE acid-washed bottles. Further information regarding geochemical analysis is provided in Supplementary Information.

Geographic spatial land use data was analyzed using the sf package [50] for R software to perform spatial analysis. Sample sites were converted to vector data as a 1 km circle around each site and land use data was aggregated and compared by site. Land use types include probability of high temperature (inferred caldera boundary/geothermal activity), agriculture (cropland and pasture, orchards, vineyards, ornamental horticulture), land use planning and allocation (urban high human density, rural human density, conserved lands), acres as golf courses, and resorts, as well as quantity of OSDS and total effluent flux in millions gallons per day [51].

DNA was extracted from all filters according to the PowerWater DNA Extraction protocol using the DNeasy PowerWater Kit (QIAGEN, Germantown, MD). Library preparation for 16S rRNA gene sequencing was performed in 4 separate libraries due to the length of the study and quarterly sampling between 2018 and 2019 using dual-indexed primers for the MiSeq Illumina platform by protocols described in Supplementary Information [52, 53]. Sequencing was performed on an Illumina MiSeq (300 cycle, V3 chemistry kit) at the UCI Genomics High-Throughput Facility, at U of California, Irvine. Methods for Quantitative PCR (qPCR) analysis of 16S rRNA, dsrA and nirS genes are provided in Supplementary Information.

Sequence data was processed through MetaFlow|mics pipeline [54, 55]. The dada2 R package [56] was used to filter, denoise and merge the raw reads into Amplicon Sequence Variants (ASVs). Reads were truncated at positions 250F and 170R and discarded if they contained one or more bases with quality scores <2 or >3. Sequencing error probabilities were modeled with dada2's iterative learn errors function and denoised with an iterative partitioning algorithm with default parameters. Reads were merged and any pairs with an overlap of fewer than 20 bases, or with more than one mismatch, were discarded. Both the mothur v1.44.1 [57] and the Silva database (version 138) [58] were used to align and annotate sequences. Sequences with a start or stop position outside the 5th-95th percentile range (over all sequences) were discarded. Potential chimeras were removed with VSEARCH [59] as implemented in Mothur and assigned taxonomy using the RDP classifier [60]. ASVs with no taxonomic information at phylum level, or matching mitochondria or chloroplasts were discarded. Sequences can be referenced at BioProject ID: PRJNA819449. Statistical analysis was performed using the Phyloseq [61], Vegan [62], DESeq2 [63] factoextra [64], sf [50], and base statistical packages for R [65]. Nutrient and major ion concentrations were z-score transformed (mg L−1) prior to Principal Component Analysis (PCA). ASVs were subsampled to 2000 for Beta diversity (CCA, NMDS, PERMANOVA) analysis, and PERMANOVA tests were performed with 99999 permutations. An initial beta-diversity analysis was performed comparing the 0.8 and 0.2 µm 16S community and no significant difference (Permanova p > 0.5) was observed and therefore samples were combined for further analysis.

Putative functional annotations based on 16S rRNA genes were made utilizing version 1.2.4 of FAPROTAX (Functional Annotation of Prokaryotic Taxa; [66]. Due to recent taxonomic updates in the classification of certain clades of microorganisms [67] that are reflected in the Silva database (Silva138) but not in the FAPROTAX database, we manually re-annotated the FAPROTAX taxa of putative N and S cyclers to match Silva138 using the Genome Taxonomy Database (GTDB). A table of these taxa reclassifications along with their associated metabolisms are provided in Supplementary Information.

Differential abundance analysis (DESeq2) was designed for gene-level expression analysis of RNA-seq data, but is adapted and recommended for detection of differentially abundant species in high-throughput data as utilized in our analysis [68]. Differential abundance testing was performed on genera with N and S cycling capabilities with ASV prevalence greater than 3. A pseudo count of 1 was used to prevent null geometric means in the normalized process. P-values were computed using a Wald test and multiple hypothesis testing was adjusted with the Benjamini–Hochberg method [69]. Differentially abundant ASVs were extracted from all pairwise group comparisons, and selected for display in the heat map when the adjusted p-value was below 0.1 and the log2 fold-change above 1. The abundance values displayed were summed for each site, then converted to relative abundance. To improve heat map visualization, relative abundances values were converted to z-scores per ASVs. A list of displayed ASVs from DESeq2 analysis is provided in the SI.

Based on principal component analysis (PCA), groundwater geochemistry across sampling sites clustered into two distinct groupings (PCA, Fig. 2A). Samples were categorized into an S group including five wells along the Hualālai rift zone (Fig. 1), and an N group consisting of 11 wells. Magnesium was the dominant driver (PC1 - 38.2%) of the differences between groups, followed by ammonium (NH4+; PC2 - 16.0%), which together explained 54.2% of the variation observed in geochemistry. The N groupings were further differentiated based on surrounding land use classifications. A PCA based on land use revealed four distinct groups explaining 81.2% of the observed variation (Fig. 2B). Group S formed a unique cluster based on an area with high probability of geothermal influence, while samples from group N were further clustered into N1 (rural land use), N2 (urban and OSDS), and N3 groups (agriculture; Fig. 2B). PCA revealed OSDS quantities were the dominant driver (PC1 - 48.8%) of variability in samples followed by conserved lands (PC2 - 32.4%).

A principal component analysis (PCA) demonstrating associations between groundwater and A geochemistry variables, and B land-use, and geology. Summary variables are colored by sample groups (N1-square, S-circle, N2-triangle, and N3-diamond) with the larger shape representing the mean of principle variables, and black vectors representing strength of association inferred from vector length. Geochemistry includes pH; Cr Chromium, DO dissolved oxygen, NOx nitrate+nitrite, TN total nitrogen, NH4 ammonium, Mn Manganese, Fe Iron, Temp temperature, Cl Chlorine, SPC specific conductance, PO4 orthophosphate, TP total orthophosphate, SO4 sulfate, Mg Magnesium. Land use includes conservation lands (LU_Consv), agriculture (Ag_Acres), high dentistry urban (LU_HDUrban), rural (LU_Rural), golf courses (Glf_Acres), resorts (LUResorts), OSDS quantities (OSDS_Qty), and OSDS N flux (OSDS_NFlux), and 35% probability of geothermal activity (Prob_HT35).

Geochemistry of the S group was significantly different from N2 and N3 (Tukey, p < 0.001) for DO, NOx−, SO4−2, PO4−3, Si, SPC, salinity, and pH (Fig. 3). Group S was characterized by a significantly higher mean, and variable sulfate concentrations (179.02 mg L−1; Fig. 3). Several samples from one well in group S had SO42− concentrations greater than the EPA secondary maximum contaminant level threshold based on taste considerations (250 mg L−1) [70]. The S group had lower mean concentrations of DO (5.03 mg L−1), NOx− (0.88 mg L−1), pH (6.93), and high mean SPC (1.57 mS cm−1). In addition, mean averages of PO43− (0.29 mg L−1) and Si (33.91 mg L−1) were significantly higher for group S. In comparison, NOx− was significantly higher in N2 and N3 (Tukey, p < 0.001) than group S. Samples from N3 had the highest mean concentrations of NOx− (1.21 mg L−1) but are not significantly different from N2 (1.06 mg L−1). Two outlier samples of NOx− occurred from N2 (5.73 mg L−1) and N3 (8.28 mg L−1), but were still below the EPA maximum contaminant level (mcl) for nitrates in groundwater (10 mg L−1) [21]. Groups in N also have significantly lower SPC (Tukey, p < 0.001), and higher pH (Tukey, p < 0.001; Fig. 3). Significant differences in geochemistry were not observed for sample times (months/year; Supplementary Fig. 1). Water quality parameters are published and available at the Hydroshare database [71].

The Mean (horizontal line), and individual sample values (dark circles) are indicated for each variable by group. Vertical y-axis for DO, NH4, NOx, and SO4−2, PO43−, and Si are measured in mg L−1, SPC (µS cm−1), and Temp (°C). Post-hoc analysis (Tukey-post hoc) was performed to determine significant differences between groups and are labeled by letters to indicate significant differences.

A mantel test revealed significant spatial relationships between geochemistry and agriculture, urban density, golf course, and resort land use types (Table 1). Geochemical variables including DO (R2 = 0.66, p = 0.001), pH (R2 = 0.41, p < 0.05), SPC (R2 = 0.31, p = 0.01), and SO4−2 (R2 = 0.90, p = 0.001) were all strongly correlated with golf course and resort land use. The S group has the largest average land use as golf course and resort lands (1679 km2) and areas of high geothermal probability (Fig. 1). In contrast, NOx− concentrations were strongly correlated with agricultural lands (R2 = 0.45, p < 0.01), and urban density (R2 = 0.32, p < 0.01; Table 1). The greatest average acreage (5160 km2) classified as urban occurred in group N3, while agricultural land use (pasture, conserved land) was greatest in group S and N1 with an average of 1.1 × 106 km2. Mean quantities of OSDS were greatest in N2 at 5146 units per parcel, compared with 4009 in N3, and the smallest quantities occurring in S at 143.7 units per parcel. Rural land use was greatest by far in group S at 7080 km2 by comparison to all N groups with an average size of 1597 km2.

Microbial density, as estimated from qPCR 16S rRNA gene copies mL−1, was low and significantly different across geochemistry groups (ANOVA, p < 0.001). The largest mean density was measured in N3 (9.7 × 104 copies mL−1; Fig. 4) and the smallest abundances were measured in S (2.5 × 103 copies mL−1). Group N3 had a significantly greater density than S (Tukey, p < 0.01), and N2 (Tukey, p < 0.001), but was not different from N1.

Post-hoc analysis (Tukey-post hoc) was performed between groups and are labeled by letters to indicate significant difference.

A total of 12,973 ASVs were identified in 220 samples collected between May 2017 and March 2019. Mean sample size was 12,966 reads per sample with a range of 2079–42,804. All samples were dominated by Proteobacteria (Supplementary Fig. 2) with mean relative abundances of 43.6% for Gammaproteobacteria (Ɣ) and 12.9% for Alphaproteobacteria (α). The top 10 most abundant orders (Supplementary Fig. 2) made up >50% of overall relative abundances of groundwater samples. Beta-diversity of the microbial community was significantly different between PCA (S, N1, N2, N3) groups (PERMANOVA, p < 0.01; Supplementary Fig. 3). However, alpha-diversity (Pielou and Shannon indices) was not significantly different between PCA groups (Supplementary Fig. 4).

Canonical correspondence analysis (CCA) identified significant associations between groundwater microbial community structure, geochemistry, and land-use (ANOVA of CCA model, p < 0.001; Fig. 5). Agricultural land use (CCA1), pH, OSDS quantity, and NOx− were significantly associated with microbial communities from groups N2 and N3, compared to resort land use and SO4−2 concentrations (CCA2) associations in groups N1 and S. Furthermore, mantel tests (Table 1) reveal that microbial communities were correlated with OSDS density (R2 = 0.43, p < 0.05), and rural densities (R2 = 0.32, p < 0.05). The mantel tests did not detect any significant correlation between community structure and agriculture, urban densities, rift zone, golf courses or resort lands, conservation lands, or any geochemical variable (Table 1).

Circles represent samples for each well colored by group, and black vectors represent association between communities, geochemistry and land use.

A total of 865 ASVs across groups (6.7% of total relative abundances) were classified as putative N and/or S cyclers between genus and order level based on Functional Annotation of Prokaryotic Taxa (FAPROTAX) [66]. Relative sequence abundances were calculated for each PCA group to account for differences in 16S rRNA gene abundances across groups.

Relative abundances of putative S-cyclers (353 ASVs) were significantly different across groups (p < 0.001; Fig. 6A). Group S had the highest mean relative abundances of putative S-cycling organisms (26.5%) and were significantly different from N1 (11.8%; Tukey, p < 0.01). Group N2 (25.2%) had significantly higher mean relative abundances of S-cyclers than N3 (18.8%; Tukey, p < 0.05) and N1 (Tukey, p = 0.0012). Putative S-oxidizers (159 ASVs) were dominated by the genus Acinetobacter with the highest mean relative abundances in group S (20.5%), compared with N2 (18.7%), N3 (15.5%), and lowest in N1 (7.9%; Fig. 6D). Sulfur reducers were dominated by an uncultured genus classified to order Desulfobulbales with low mean relative abundances (1.0%) across groups, and with the largest relative abundances occurring in group S (1.5%; Fig. 6E). Abundance of the sulfate reducing gene (dsrA) was also very low, but significantly different across all PCA groups (ANOVA, p < 0.001; Fig. 5). Group N1 had significantly higher mean abundances of dsrA (2.2 × 103 copies mL−1) than N2 (p < 0.01). Group S (1.89 × 102 copies mL−1) was also significantly different (p < 0.01) than N1 (2.15 × 103 copies mL−1). If we assume an average of one dsrA gene copy per genome [72, 73], 7.7% of the microbial community of group S has the potential for S reduction, compared to 5.3% at N1, 1.1% at N2, and 0.94% at N3.

Relative abundance (%) of ASVs are displayed by functional PCA groups including A N and S Metabolisms, B Denitrification included nitrate, nitrite, and nitrous oxide denitrification, C N Reduction includes N or S reduction, and nitrate reduction, D S oxidation includes photoautotrophic S oxidation, and E S Reduction.

Putative N-cycling microbes (512 ASVs) also had significantly different relative abundances across groups (ANOVA, p < 0.001; Fig. 6A). Group S had the highest mean relative abundance of N-cyclers (28.9%), and was significantly different from N3 (20.91%; Tukey, p < 0.05) and N1 (15.27%; Tukey p < 0.01). Groups S and N1 had similar relative abundances of nitrifying organisms compared to lower abundances in groups N2 and N3. Most nitrifying microbes belonged to genus Candidatus Nitrosotenuis (Archaea) in N groups, compared to genus Candidatus Nitrotoga (Bacteria) in group S. Microorganisms capable of dissimilatory N reduction (via any part of the reduction cascade; 396 ASVs) had high mean relative abundances (18.5%) across all groups and were dominated by five ASVs of Ɣ-proteobacteria in genus Acinetobacter, with the largest relative abundance occurring in group S (11.1%; Fig. 6C). The next most abundant N reducer was classified to genus Aquabacterium with the highest relative abundance in group N3 (1.1%; Fig. 6C).

Denitrifying (capable of reducing nitrate, nitrite, or nitrous oxide to N2) microbes (212 ASVs) had low mean abundances (3.14%) and were not significantly different by groups. The largest denitrifying taxa contributions were from genus Denitratasoma with the largest relative abundance in group N2 (0.45%) and S (0.42%) and genus Flavobacterium in group S (0.38%; Fig. 6B). A more accurate measurement of denitrification potential measured by qPCR of the denitrifying gene (nirS) demonstrated high and significant differences between all groups (ANOVA, p < 0.001; Fig. 5). Group N3 had the highest mean nirS abundances (4.5 × 104 copies mL−1), followed by S (3.88 × 103 copies mL−1), and N1 (3.76 × 103 copies mL−1) and N2 (1.0 × 103 copies mL−1). The N3 group samples were significantly different from N2 (p < 0.001) and S (p < 0.01) for nirS gene abundances. Estimating an average one nirS gene copy per genome [74, 75] the majority of group S microbial community (>100%) has the potential for denitrification, compared to 9.3% at N1, 30.1% at N2 and 46.3% at N3.

Differential abundance analysis (DESeq2) revealed unique clustering of significantly differentiated ASVs (adjusted p value < 0.1) across groups with similar functional potential for N and S metabolisms (Fig. 7). Group N1 and S had the greatest number of differentially enriched ASVs (68 and 64 respectively), compared with fewer in groups N2, and N3 (35 and 38 respectively). Putative N functional metabolisms (e.g., aerobic nitrification, denitrification, and dissimilatory N-reduction) are present across all groups, but differ by taxa. Putative S functional metabolisms (e.g., S-oxidation, and S-reduction) are highly enriched in groups S and N1, dominated by a putative S-oxidizing, N-reducer (Acinetobacter) in group S. A summary of the significantly different ASVs with the highest enrichment (mean fold change) organized by N and S metabolism is discussed in Supplementary Information and Supplementary Table 1. Results presented here are focused on the dominant putative S-oxidizing, N-reducing taxa differentiated across PCA groups.

Abundance values for each ASV selected were summed by site and converted to relative abundance and transformed to z-scores for each ASV. ASV numbers identified as Acinetobacter are highlighted.

Putative S-oxidation was highly enriched in group S (27 ASVs) and N1 (20 ASVs), dominated by taxa capable of both S-oxidation, and N-reducing taxa (Fig. 7; Supplementary Table 1). The highest enriched ASV in group S belonged to a putative S-oxidizing and N-reducer, genus Acinetobacter (12 ASVs). However, another putative S-oxidizing, N-reducer, Thiobacillus (3 ASVs) had the highest enrichment across all taxa in group N1. Group N2 taxa was highly enriched in Acinetobacter, but also Pseudomonas, while group N3 had unique contributions from Rhodobacter, and Dechloromonas. In addition, taxa capable of putative S-oxidation not coupled with N reduction was diverse across groups, with unique contributions in groups N1, S and N2. Group S had contributions from Sulfuricurvum (2 ASVs) and Methylobacterium-Methylorubrum (3 ASVs), compared to putative S-oxidizers in group N1 including Sulfurifustis (2 ASVs), Magnetovibrio (1 ASV), and Meiothermus (1 ASV) and contributions from Chromatiaceae (1 ASV), and Chlorobium (1 ASV) in group N2.

Groundwater samples from group S, which are located near the northwest rift zone (Fig. 1), were biogeochemically unique and characterized by high SO42−. Samples from group S have a distinct volcanic CO2 signal suggesting geothermal outgassing [18], which can include sulfur gasses (e.g., H2S, SO22−) [76], as well as elevated concentrations of PO43− and Si as in our samples, indicating increased weathering of rock minerals. Groundwater age dating from the rift zone (C14) indicates much older water (>5000 years) due to either excess geothermal-derived inorganic C, or isolated flow paths resulting in decreased connectivity with the rest of the aquifer [77, 78]. Decreased recharge, and older age indicates a diminished likelihood of abiotic production of oxidized S species such as SO42−, and S2O32− from volcanic outgassing. However, biotic H2S oxidation rates have been shown to far exceed abiotic oxidation in both aerobic and anaerobic conditions [36], further suggesting that SO42− production in these samples may be produced by microbial activity.

Multivariate analysis demonstrated that SO42− concentrations are significant drivers of microbial community structure in both groups S and N1 (Fig. 4 CCA), consistent with other groundwater microbiome studies [5]. The S group was dominated by a high diversity (45 of the 159 ASVs) of putative S-oxidizers belonging to the genus Acinetobacter. The genus has been characterized in oligotrophic groundwater with a versatile genome and a great potential for water treatment [79, 80]. FAPROTAX classified Acinetobacter as a N-reducer [66], however more recent work has shown members have capabilities for nitrification (oxidation of ammonia to nitrate) using organic substrates (heterotrophic) and aerobic denitrification [80]. Several species of Acinetobacter have been described in several industrial systems including wastewater treatment of sulfur, Mn [81], and H2 [79]. Acinetobacter strains were first identified in the process of denitrifying sulfide removal using microbial communities capable of both heterotrophy and autotrophy to gain energy from S-oxidation with a complete set of genes to reduce NO32− to N2 [81].

Taxa classified to Acinetobacter were 4-times more enriched in group S than any other N groups. The presence of Acinetobacter in our samples may indicate a direct linking of S and N functional metabolisms [82, 83], which is also supported by high nirS gene abundances throughout our samples. Acinetobacter contains one copy of the nirS gene [80] and high nirS relative abundances suggests that there is great potential for microbially-driven S-oxidation and N-removal in group S. This study identifies the functional capacity in Hawaiian volcanic aquifers for an Acinetobacter-type groundwater microbial community that has the potential to utilize energy from S to drive complete N removal. Acinetobacter has demonstrated high N removal efficiency under aerobic conditions (>2 mg L−1 DO) in low carbon, low temperature systems where it can remove up to 40.2% of N as gas at a rate of 0.203 mg L−1 h−1 [84]. The high sulfate concentrations coupled with significantly lower NOx− in group S compared to groups N2 and N3, further supports our hypothesis that the electrons required for N-reduction originates from microbial oxidation of sulfur producing SO42− by Acinetobacter (S group) and possibly Thiobacillus (N1 group).

Land management influences the N geochemistry in the Hualālai watershed primarily from impacts of urbanization (OSDS quantity) in groups N2 and N3 (Keauhou). Inputs of N have a major effect on microbial community structure in groundwater, and may differ depending upon the anthropogenic source of the N [6, 23, 85]. In this study, quantities of OSDS are the primary driver of differences (48%; Fig. 2B) and were significantly associated with microbial community structure (Fig. 4). Results demonstrating OSDS as a driver of groundwater geochemistry are consistent with a groundwater flow and nutrient transport model developed for the coastal region of the Keauhou aquifer (N3) [35]. The model shows that OSDS contributed the largest proportion of nutrients (54% of total N) to the aquifer and had the greatest effect on water quality relative to other nonpoint source contaminants [86]. In this study, areas of high urban densities with high OSDS quantities (N2, N3) had significantly higher NOx− concentrations than comparable areas dominated by passive agricultural land use (N1).

Human impacts from agriculture and wastewater are important determinants of geochemistry, water quality, and microbial community structure (especially N-cycling functional groups) in groundwater ecosystems [12, 21]. In this study, the highly enriched putative microbial N-removal function is complete denitrification (heterotrophic nitrification coupled with aerobic denitrification) via sulfur oxidation by Acinetobacter. We hypothesize a greater potential for microbial-mediated N-removal in group S compared to the N groups due to the presence of S compounds, elevated sulfate concentrations, and highly enriched Acinetobacter populations.

The availability of labile organic carbon is thought to limit the viability of heterotrophic N-removal in groundwater [24, 87]. Sewage effluent that is well oxidized also tends to contain smaller amounts of labile carbon [23] limiting microbial N-reduction processes in some environments. The potential for microbial mediated N-reduction is greatest in group N3 compared to group S based on taxa comparisons and abundances (nirS), however, group N3 may be limited by electron donor availability which is greater in group S due to the presence of S compounds. Further studies are necessary to delineate the microbial N-reducution capabilities in these groundwater aquifers and the differences that appear across our groups, especially in areas heavily influenced by OSDS as observed in groups N2 and N3.

Findings for the enriched potential of linked chemolithoautrophic sulfur-oxidizing and denitrification in the Hualālai watershed has implications for water managers. Although denitrification by sulfide oxidation leads to decreased N loading, the metabolism is potentially detrimental to well or pump operations. Sulfate concentrations, water hardness, and corrosion may increase, thereby causing ions to precipitate out of solution when oxygen is encountered as groundwater is pumped to the surface [6, 23]. Personal communication with private well owners along the rift zone confirmed that some wells do experience rancid, foul-smelling water and precipitation. However, the significant presence of Acinetobacter infers the bio-remediation potential of volcanic groundwater microbial communities to simultaneously remove N and S providing an ecosystem service for resource-constrained groundwater aquifers.

The datasets generated and analyzed during the current study are available at the Hydroshare repository, https://www.hydroshare.org/resource/d812bbb7c93348999371c9f1f517297f/ and the National Center for Biotechnology Information (NCBI) repository, ID PRJNA819449. All further data and analysis discussed in this study is provided in manuscript supplemental information.

Anantharaman K, Brown CT, Hug LA, Sharon I, Castelle CJ, Probst AJ, et al. Thousands of microbial genomes shed light on interconnected biogeochemical processes in an aquifer system. Nat Commun. 2016;7:13219.

Article CAS PubMed PubMed Central Google Scholar

Korbel KL, Hancock PJ, Serov P, Lim RP, Hose GC. Groundwater ecosystems vary with land use across a mixed agricultural landscape. J Environ Qual. 2013;42:380–90. https://doi.org/10.2134/jeq2012.0018.

Article CAS PubMed Google Scholar

Whitman WB, Coleman DC, Wiebe WJ. Prokaryotes: the unseen majority. Proc Natl Acad Sci USA. 1998;95:6578–83.

Article CAS PubMed PubMed Central Google Scholar

Griebler C, Lueders T. Microbial biodiversity in groundwater ecosystems. Freshwater Biol. 2009;54:649–77. https://doi.org/10.1111/j.1365-2427.2008.02013.x.

Article Google Scholar

Flynn TM, Sanford RA, Santo Domingo JW, Ashbolt NJ, Levine AD, Bethke CM. The active bacterial community in a pristine confined aquifer. Water Resour Res. 2012;48. https://onlinelibrary.wiley.com/doi/10.1029/2011WR011568.

Seitzinger S, Harrison JA, Böhlke JK, Bouwman AF, Lowrance R, Peterson B, et al. Denitrification across landscapes and waterscapes: a synthesis. Ecol Appl. 2006;16:2064–90.

Article CAS PubMed Google Scholar

Bethke CM, Sanford RA, Kirk MF, Jin Q, Flynn TM. The thermodynamic ladder in geomicrobiology. Am J Sci. 2011;311:183–210.

Article CAS Google Scholar

Kirs M, Kisand V, Nelson CE, Dudoit T, Moravcik PS. Distinct bacterial communities in tropical island aquifers. PLoS One. 2020;15:e0232265.

Article CAS PubMed PubMed Central Google Scholar

Sirisena KA, Daughney CJ, Moreau-Fournier M, Ryan KG, Chambers GK. National survey of molecular bacterial diversity of New Zealand groundwater: relationships between biodiversity, groundwater chemistry and aquifer characteristics. FEMS Microbiol Ecol. 2013;86:490–504.

Article CAS PubMed Google Scholar

Hubalek V, Wu X, Eiler A, Buck M, Heim C, Dopson M, et al. Connectivity to the surface determines diversity patterns in subsurface aquifers of the Fennoscandian shield. ISME J. 2016;10:2556.

Article PubMed PubMed Central Google Scholar

Hemme CL, Tu Q, Shi Z, Qin Y, Gao W, Deng Y, et al. Comparative metagenomics reveals impact of contaminants on groundwater microbiomes. Front Microbiol. 2015;6:1205.

Article PubMed PubMed Central Google Scholar

Korbel KL, Greenfield P, Hose GC. Agricultural practices linked to shifts in groundwater microbial structure and denitrifying bacteria. Sci Total Environ. 2022;807:150870.

Article CAS PubMed Google Scholar

Hawaiʻi Fresh Water Initiative. A blueprint for action: water security for an uncertain future. Honolulu, HI, USA: Hawaiʻi Community Foundation. 2016. https://www.hawaiicommunityfoundation.org/file/cat/Fresh_Water_Blueprint_FINAL_062215_small.pdf.

Whittier RB, El-Kadi A. Human health and environmental risk ranking of on-site sewage disposal systems for the Hawaiian Islands of Kauai, Molokai, Maui, and Hawaii. Honolulu, HI: Final report prepared for State of Hawai’i Department of Health, Safe Drinking Water Branch; 2014.

Amato DW, Bishop JM, Glenn CR, Dulai H, Smith CM. Impact of submarine groundwater discharge on marine water quality and reef Biota of Maui. PLoS One. 2016;11:e0165825.

Article PubMed PubMed Central Google Scholar

Delevaux JMS, Whittier R, Stamoulis KA, Bremer LL, Jupiter S, Friedlander AM, et al. A linked land-sea modeling framework to inform ridge-to-reef management in high oceanic islands. PLoS One. 2018;13:e0193230.

Article PubMed PubMed Central Google Scholar

McKenzie T, Habel SL, Dulai H. Increased coastal pollution expected under future sea level stands: chemical evidence for tidal groundwater inundation of coastal wastewater infrastructure. In: Geological Society of America Abstracts with Programs. 2019. https://doi.org/10.1130/abs/2019am-330663.

Okuhata BK, El-Kadi AI, Dulai H, Lee J, Wada CA, Bremer LL, et al. A density-dependent multi-species model to assess groundwater flow and nutrient transport in the coastal Keauhou aquifer, Hawai ‘i, USA. Hydrogeol J. 2022;30:231–50.

Article Google Scholar

Menció A, Mas-Pla J, Otero N, Regàs O, Boy-Roura M, Puig R, et al. Nitrate pollution of groundwater; all right, but nothing else? Sci Total Environ. 2016;539:241–51.

Article PubMed Google Scholar

Kim HR, Yu S, Oh J, Kim KH, Lee JH. Moniruzzaman, et al. Nitrate contamination and subsequent hydrogeochemical processes of shallow groundwater in agro-livestock farming districts in South Korea. Agric Ecosyst Environ. 2019;273:50–61.

Article CAS Google Scholar

Bouwman AF, Beusen AHW, Griffioen J, Van Groenigen JW, Hefting MM, Oenema O, et al. Global trends and uncertainties in terrestrial denitrification and N2O emissions. Philos Trans R Soc Lond B Biol Sci. 2013;368:20130112.

Article CAS PubMed PubMed Central Google Scholar

Slomp CP, Van, Cappellen P. Nutrient inputs to the coastal ocean through submarine groundwater discharge: controls and potential impact. J Hydrol. 2004;295:64–86.

Article CAS Google Scholar

Rivett MO, Buss SR, Morgan P, Smith JWN, Bemment CD. Nitrate attenuation in groundwater: a review of biogeochemical controlling processes. Water Res. 2008;42:4215–32.

Article CAS PubMed Google Scholar

Santoro AE, Boehm AB, Francis CA. Denitrifier community composition along a nitrate and salinity gradient in a coastal aquifer. Appl Environ Microbiol. 2006;72:2102–9.

Article CAS PubMed PubMed Central Google Scholar

Hashimoto S, Furukawa K, Shioyama M. Autotrophic denitrification using elemental sulfur. J Ferment Technol. 1987;65:683–92.

Article CAS Google Scholar

Burgin AJ, Hamilton SK. Have we overemphasized the role of denitrification in aquatic ecosystems? A review of nitrate removal pathways. Front Ecol Environ. 2007;5:89–96.

Article Google Scholar

Robertson LA, Kuenen JG. Thiosphaera pantotropha gen. nov. sp. nov., a Facultatively Anaerobic, Facultatively Autotrophic Sulphur Bacterium. Microbiology. 1983;129:2847–55.

Article CAS Google Scholar

Hiscock KM, Lloyd JW, Lerner DN. Review of natural and artificial denitrification of groundwater. Water Res. 1991;25:1099–111.

Article CAS Google Scholar

Howarth RW, Marino R. Nitrogen as the limiting nutrient for eutrophication in coastal marine ecosystems: Evolving views over three decades. Limnol Oceanogr. 2006;51:364–76.

Article CAS Google Scholar

Street JH, Knee KL, Grossman EE, Paytan A. Submarine groundwater discharge and nutrient addition to the coastal zone and coral reefs of leeward Hawai’i. Mar Chem. 2008;109:355–76.

Article CAS Google Scholar

Duarte TK, Pongkijvorasin S, Roumasset J, Amato D, Burnett K. Optimal management of a Hawaiian Coastal aquifer with nearshore marine ecological interactions. Water Resour Res. 2010;46. http://doi.wiley.com/10.1029/2010WR009094.

Di Capua F, Pirozzi F, Lens PNL, Esposito G. Electron donors for autotrophic denitrification. Chem Eng J. 2019;362:922–37.

Article Google Scholar

Xu Y, Schoonen MAA, Nordstrom DK, Cunningham KM, Ball JW. Sulfur geochemistry of hydrothermal waters in Yellowstone National Park: I. the origin of thiosulfate in hot spring waters. Geochim Cosmochim Acta. 1998;62:3729–43.

Article CAS Google Scholar

Stefánsson A, Arnórsson S, Gunnarsson I, Kaasalainen H, Gunnlaugsson E. The geochemistry and sequestration of H2S into the geothermal system at Hellisheidi, Iceland. J Volcanol Geotherm Res. 2011;202:179–88.

Article Google Scholar

Kaasalainen H, Stefánsson A. The chemistry of trace elements in surface geothermal waters and steam, Iceland. Chem Geol. 2012;330-331:60–85.

Article CAS Google Scholar

Luther GW, Findlay A, MacDonald D, Owings S, Hanson T, Beinart R, et al. Thermodynamics and kinetics of sulfide oxidation by oxygen: a look at inorganically controlled reactions and biologically mediated processes in the environment. Front Microbiol. 2011;2. https://www.frontiersin.org/article/10.3389/fmicb.2011.00062.

Oki DS. Geohydrology and Numerical Simulation of the Ground-water Flow System of Kona, Island of Hawaii. U.S. Department of the Interior, U.S. Geological Survey. Honolulu, Hawai'i: Pacific Island Water Science Center; 1999. 70 p.

Fukunaga and Associates, Inc. Hawai’i county water use and development plan update. Department of Water Supply, Hawaii County; 2017. https://www.hawaiidws.org/wp-content/uploads/2018/06/Combined-Ph-1-2-Keauhou-20170510_w-Appendix-final.pdf.

Marston L, Konar M, Cai X, Troy TJ. Virtual groundwater transfers from overexploited aquifers in the United States. Proc Natl Acad Sci USA. 2015;112:8561–6.

Article CAS PubMed PubMed Central Google Scholar

Owuor SO, Butterbach-Bahl K, Guzha AC, Rufino MC, Pelster DE, Díaz-Pinés E, et al. Groundwater recharge rates and surface runoff response to land use and land cover changes in semi-arid environments. Ecol Process. 2016;5. http://ecologicalprocesses.springeropen.com/articles/10.1186/s13717-016-0060-6.

Clague DA, Dalrymple GB. The Hawaiian-Emporer volcanic chain, part 1, geologic evolution. 1987. https://evols.library.manoa.hawaii.edu/handle/10524/33604.

Izuka, SK, Engott, JA, Rotzoll, K, Bassiouni, M, Johnson, AG, Miller, LD, et al. Volcanic aquifers of Hawai’i-Hydrogeology, water budgets, and conceptual models. U.S. Geological Survey; 2018. Report No.: 2015-5164. https://doi.org/10.3133/sir20155164.

Wolfe EW, Wise WS, Brent Dalrymple G. The geology and petrology of Mauna Kea Volcano, Hawaii; a study of postshield volcanism. Professional Paper; 1997. https://doi.org/10.3133/pp1557.

Moore JG, Clague D. Coastal lava flows from Mauna Loa and Hualalai volcanoes, Kona, Hawaii. Bull. Volcanol. 1987;49:752–64. https://doi.org/10.1007/bf01079826.

Article CAS Google Scholar

Moore JG, Clague DA, Holcomb RT, Lipman PW, Normark WR, Torresan ME. Prodigious submarine landslides on the Hawaiian Ridge. J Geophys Res. 1989;94:17465 https://doi.org/10.1029/jb094ib12p17465.

Article Google Scholar

Stearns HT, Macdonald GA. Geology and ground-water resources of the island of Hawaii. Honolulu Advertising; 1946. https://pubs.er.usgs.gov/publication/70160867.

Giambelluca TW, Chen Q, Frazier AG, Price JP, Chen YL, Chu PS, et al. Online rainfall Atlas of Hawai’i. Bull Am Meteorol Soc. 2013;94:313–6.

Article Google Scholar

Korbel K, Chariton A, Stephenson S, Greenfield P, Hose GC. Wells provide a distorted view of life in the aquifer: implications for sampling, monitoring and assessment of groundwater ecosystems. Sci Rep. 2017;7:40702.

Article CAS PubMed PubMed Central Google Scholar

U.S. Geological Survey. Methods of analysis by the U.S. Geological Survey National Water Quality Laboratory: evaluation of alkaline persulfate digestion as an alternative to Kjeldahl digestion for determination of total and dissolved nitrogen and phosphorus in water. Denver, CO: U.S. Geological Survey; 2003.

Pebesma E, Bivand R. sf: Simple Features for R. R package version 0.6-3. https://CRAN.R-project.org/package=sf.

Office of Planning and Sustainable Development. 2022. https://planning.hawaii.gov/gis/download-gis-data-expanded/.

Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20.

Article CAS PubMed PubMed Central Google Scholar

Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Huntley J, Fierer N, et al. Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms. ISME J. 2012;6:1621–4.

Article CAS PubMed PubMed Central Google Scholar

Arisdakessian C, Cleveland SB, Belcaid M. MetaFlow|mics: scalable and reproducible nextflow pipelines for the analysis of microbiome marker data. In: Practice and experience in advanced research computing. New York, NY, USA: ACM; 2020. https://dl.acm.org/doi/10.1145/3311790.3396664.

Cleveland S, Arisdakessian C, Nelson C, Belcaid M, Frank K, Jacobs G. The C-MĀIKI gateway: a modern science platform for analyzing microbiome data. In: Practice and experience in advanced research computing. New York, NY, USA: Association for Computing Machinery; 2022. p. 1–7. (PEARC ’22).

Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13:581–3.

Article CAS PubMed PubMed Central Google Scholar

Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.

Article CAS PubMed PubMed Central Google Scholar

Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.

Article CAS PubMed Google Scholar

Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584.

Article PubMed PubMed Central Google Scholar

Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.

Article CAS PubMed PubMed Central Google Scholar

McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217 https://doi.org/10.1371/journal.pone.0061217.

Article CAS PubMed PubMed Central Google Scholar

Dixon P. VEGAN, a package of R functions for community ecology. J Veg Sci. 2003;14:927–30. https://doi.org/10.1111/j.1654-1103.2003.tb02228.x.

Article Google Scholar

Anders S, Huber W. Differential expression analysis for sequence count data. Nat Prec. 2010. https://doi.org/10.1038/npre.2010.4282.2.

Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25:1–18.

Article Google Scholar

Team RC. R: 2019. A Language and Environment for Statistical Computing version. 2020;3.

Louca S, Parfrey LW, Doebeli M. Decoupling function and taxonomy in the global ocean microbiome. Science. 2016;353:1272–7.

Article CAS PubMed Google Scholar

Waite DW, Chuvochina M, Pelikan C, Parks DH, Yilmaz P, Wagner M, et al. Proposal to reclassify the proteobacterial classes Deltaproteobacteria and Oligoflexia, and the phylum Thermodesulfobacteria into four phyla reflecting major functional capabilities. Int J Syst Evol Microbiol. 2020;70:5972–6016.

Article CAS PubMed Google Scholar

Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550.

Article PubMed PubMed Central Google Scholar

McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10:e1003531.

Article PubMed PubMed Central Google Scholar

USEPA. Drinking water advisory: consumer acceptability advice and health effects analysis on sulfate. 2003. Report No.: EPA 822-R-03-007. https://www.epa.gov/sites/default/files/2014-09/documents/support_cc1_sulfate_healtheffects.pdf.

Tachera D. Groundwater chemistry: nutrient data. 2021. https://www.hydroshare.org/resource/d812bbb7c93348999371c9f1f517297f/.

Klein M, Friedrich M, Roger AJ, Hugenholtz P, Fishbain S, Abicht H, et al. Multiple lateral transfers of dissimilatory sulfite reductase genes between major lineages of sulfate-reducing prokaryotes. J Bacteriol. 2001;183:6028–35.

Article CAS PubMed PubMed Central Google Scholar

Kondo R, Nedwell DB, Purdy KJ, Silva SQ. Detection and enumeration of sulphate-reducing bacteria in estuarine sediments by competitive PCR. Geomicrobiol J. 2004;21:145–57.

Article CAS Google Scholar

Kandeler E, Deiglmayr K, Tscherko D, Bru D, Philippot L. Abundance of narG, nirS, nirK, and nosZ genes of denitrifying bacteria during primary successions of a glacier foreland. Appl Environ Microbiol. 2006;72:5957–62.

Article CAS PubMed PubMed Central Google Scholar

Chon K, Chang JS, Lee E, Lee J, Ryu J, Cho J. Abundance of denitrifying genes coding for nitrate (narG), nitrite (nirS), and nitrous oxide (nosZ) reductases in estuarine versus wastewater effluent-fed constructed wetlands. Ecol Eng. 2011;37:64–9.

Article Google Scholar

Kilbride BM, Edmonds M, Biggs J. Observing eruptions of gas-rich compressible magmas from space. Nat Commun. 2016;7:13744.

Article CAS PubMed PubMed Central Google Scholar

Attias E, Constable S, Sherman D, Ismail K, Shuler C, Dulai H. Marine electromagnetic imaging and volumetric estimation of freshwater plumes offshore Hawai’i. Geophys Res Lett. 2021;16:48. https://onlinelibrary.wiley.com/doi/10.1029/2020GL091249.

Okuhata BK, Thomas DM, Dulai H, Popp BN, Lee J, El-Kadi AI. Inference of young groundwater ages and modern groundwater proportions using chlorofluorocarbon and tritium/helium-3 tracers from West Hawai’i Island. J Hydrol. 2022;609:127755 https://doi.org/10.1016/j.jhydrol.2022.127755.

Article CAS Google Scholar

Jung J, Park W. Acinetobacter species as model microorganisms in environmental microbiology: current state and perspectives. Appl Microbiol Biotechnol. 2015;99:2533–48.

Article CAS PubMed Google Scholar

Su JF, Zheng SC, Huang TL, Ma F, Shao SC, Yang SF, et al. Characterization of the anaerobic denitrification bacterium Acinetobacter sp. SZ28 and its application for groundwater treatment. Bioresour Technol. 2015;192:654–9.

Article CAS PubMed Google Scholar

Lee DJ, Pan X, Wang A, Ho KL. Facultative autotrophic denitrifiers in denitrifying sulfide removal granules. Bioresour Technol. 2013;132:356–60.

Article CAS PubMed Google Scholar

Cardoso RB, Sierra-Alvarez R, Rowlette P, Flores ER, Gómez J, Field JA. Sulfide oxidation under chemolithoautotrophic denitrifying conditions. Biotechnol Bioeng. 2006;95:1148–57.

Article CAS PubMed Google Scholar

Chen C, Ren N, Wang A, Yu Z, Lee DJ. Simultaneous biological removal of sulfur, nitrogen and carbon using EGSB reactor. Appl Microbiol Biotechnol. 2008;78:1057–63. https://doi.org/10.1007/s00253-008-1396-3.

Article CAS PubMed Google Scholar

Wen G, Wang T, Li K, Wang H, Wang J, Huang T. Aerobic denitrification performance of strain Acinetobacter johnsonii WGX-9 using different natural organic matter as carbon source: effect of molecular weight. Water Res. 2019;164:114956 https://doi.org/10.1016/j.watres.2019.114956.

Article CAS PubMed Google Scholar

Ben Maamar S, Aquilina L, Quaiser A, Pauwels H, Michon-Coudouel S, Vergnaud-Ayraud V, et al. Groundwater isolation governs chemistry and microbial community structure along hydrologic flowpaths. Front Microbiol. 2015;6:1457.22.

PubMed PubMed Central Google Scholar

Article Google Scholar

Bellini MI, Kumaresan D, Tarlera S, Murrell JC, Fernández-Scavino A. Identification of active denitrifiers by DNA-stable isotope probing and amplicon sequencing reveals Betaproteobacteria as responsible for attenuation of nitrate. FEMS Microbiol Ecol. https://academic.oup.com/femsec/article-abstract/94/2/fix181/4757058.

Download references

We would like to thank the numerous Kona, Hawai’i community groups for generously giving their time, site access, and knowledge to help us conduct this work. Thanks go to the Hawai’i County Department of Water Supply, Hawai’i Water Service, Natural Energy Laboratory of Hawai’i Authority including Pam Madden, and Keith Olson, Commission on Water Resource Management, Kohanaiki Club, and Greg Chun for working with us to access wells and helping us to collect groundwater samples. Mahalo nui to Sean Cleveland, ITS and the Hawaii Data Science Institute for technical assistance and cyberinfrastructure support. This paper is SOEST Contribution No. (TBD after publication), University of Hawai’i at Mānoa.

This research was supported by the Hawai’i EPSCoR Program funded by the National Science Foundation Research Infrastructure Improvement Awards (RII) Track-1: ‘Ike Wai: Securing Hawaii's Water Future Award #OIA-1557349 and Advanced CyberInfrastructure - Tapis Framework #1931439 and #1931575. This work was made possible, in part, through access to the Genomics High Throughput Facility Shared Resource of the Cancer Center Support Grant (P30CA-062203) at the University of California, Irvine and NIH shared instrumentation grants 1S10RR025496-01, 1S10OD010794-01, and 1S10OD021718-01.

University of Hawai’i at Mānoa, Pacific Biosciences Research Center, Honolulu, HI, USA

Sheree J. Watson, Cédric Arisdakessian, Maria Petelo, Kekuʻiapōiula Keliipuleole & Kiana L. Frank

University of Hawai’i at Mānoa, Department of Information and Computer Sciences, Honolulu, HI, USA

Cédric Arisdakessian

University of Hawai’i at Mānoa, Marine Biology Graduate Program, Honolulu, HI, USA

Kekuʻiapōiula Keliipuleole

University of Hawai’i at Mānoa, Department of Earth Sciences, Honolulu, HI, USA

Diamond K. Tachera, Brytne K. Okuhata & Henrietta Dulai

You can also search for this author in PubMed Google Scholar

SW and KF conceived and designed the work that led to the submission of the manuscript including acquiring samples, data analysis and interpreting the results. CA was instrumental in data analysis and interpreting the results. KF, KK, BO, DT, and SW played an important role in collection of the data, MP, SW and KK were instrumental in lab preparation of the samples. The manuscript was drafted by SW, KF, and CA. Co-authors BO, and HD helped revise the manuscript and helped approve the final version. All co-authors agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Correspondence to Kiana L. Frank.

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 Creative 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 permitted 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/.

Reprints and Permissions

Watson, S.J., Arisdakessian, C., Petelo, M. et al. Geology and land use shape nitrogen and sulfur cycling groundwater microbial communities in Pacific Island aquifers. ISME COMMUN. 3, 58 (2023). https://doi.org/10.1038/s43705-023-00261-5

Download citation

Received: 09 September 2022

Revised: 18 April 2023

Accepted: 24 May 2023

Published: 07 June 2023

DOI: https://doi.org/10.1038/s43705-023-00261-5

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Blog

Geology and land use shape nitrogen and sulfur cycling groundwater microbial communities in Pacific Island aquifers