Molecular Nature of Marine Particulate Organic Iron-Carrying Moieties Revealed by Electrospray Ionization Fourier-Transform Ion Cyclotron Resonance Mass Spectrometry (ESI-FTICRMS)

Marine sinking particulate organic matter (POM), acting as a link between surface primary production and burial of organic matter in marine sediments, undergoes a variety of physical and biochemical alterations on its way to the deep ocean, resulting in an increase in its un-characterizable proportion with diagenesis. Further, the binding ligands in POM for iron, an essential nutrient to marine life and tightly coupled with organic matter, has rarely been studied. In the current study, we employed an approach combining sequential extraction with ultrahigh resolution mass spectrometry (ESI-FTICRMS), in order to explore and unravel the chemical characteristics of organic matter compounds relevant to marine particle flux within the mesopelagic and deep ocean, with a focus on the potential iron-carrying molecules. With increasing depth, POM increases in aliphaticity, and decreases in intensity-normalized O/C ratios, aromatics, and carboxylic-rich alicyclic molecules (CRAM)-like compounds. The potential iron-carrying molecules account for ∼14% of total identified molecules, and appear to have been incorporated into the marine particles via ion complexation, hydrophobic interaction, and/or interlayered “occlusion.” The relative abundance of iron-binding organic molecules in these three operationally-defined categories changes with depth: “surficially-complexed” fraction decreases with depth, the “interlayered-occluded” fraction increases to a comparable extent and “hydrophobic interaction” fraction occurs at all depths. Collectively, the potential iron-carrying organic molecules exhibit a set of unique molecular characteristics: a relatively lower average H/C ratio and a higher O/C ratio compared to bulk POM, a dominance of aromatics, black carbon-like compounds and CRAM-like compounds, and minor amounts of aliphatics. These molecules exhibit partial similar molecular features as precursors formed from photochemical reactions in the surface ocean, but they have been greatly modified by flux processes. Noticeably, a minor fraction of these iron-carrying molecules (<1%) was identified to contain hydroxamate-like moieties [N(OH)-COOH], the key functionality of one of the strongest iron-binding ligands in the dissolved phase. This agrees with improved spectrophotometric results and corroborates their presence in the POM. These hydroxamate-like moieties play an important role in controlling the distributions and fluxes of Fe and particle-reactive radionuclides with similar chemical complexing properties as Fe (e.g., thorium) in the ocean.


INTRODUCTION
Sinking particles, composed of organic (sinking particulate organic matter, POM) and inorganic material, are the critical conduit for delivering photosynthetically-fixed organic carbon, food and energy from the surface ocean to the deep sea (Verdugo et al., 2004;Boeuf et al., 2019). They vary in size and shape, consisting of a suit of dead phytoplankton, zooplankton carcasses, larvacean houses, fecal pellet, etc. (Omand et al., 2020). Typically accounting for <20% of the surface primary production, sinking POM forms through a variety of biological, physical and chemical processes that alter its composition and attenuate the flux on its way to the sea floor (Volkman and Tanoue, 2002). Despite its paramount role in global carbon cycling, sinking POM is still a poorly understood organic matter pool in terms of the physical and chemical processes it undergoes while sinking (Lee, 2002).
Iron, as a trace element in the dissolved phase (dFe: 0.02-1.0 nmol L −1 ; Bruland et al., 1991) at the surface ocean, is tightly coupled with organic carbon in both of their biogeochemical cycles and also regulates the export of POM flux from the euphotic zone to the deep ocean. For example, in high-nutrient, low-chlorophyll (HNLC) regions, the low surface ocean dFe concentration can limit primary productivity by as much as 30-40% (Marchetti and Maldonado, 2016). In response to the scarcity of Fe, microorganisms produce a variety of lowmolecular-weight (LMW, <1 kDa) organic compounds, i.e., siderophores, with strong binding affinity to metals, especially Fe(III) (Reid et al., 1993;Neilands, 1995;Butler, 1998Butler, , 2005Winkelmann, 2002;Vraspir and Butler, 2009;Butler and Theisen, 2010). Siderophores are typically multidentate, oxygen-donor ligands that include hydroxamate, catecholate and α-hydroxycarboxylate types. Out of these three types of siderophores, hydroxamate has been the focus of many studies on marine dFe binding ligands (Macrellis et al., 2001;Mawji et al., 2008;Velasquez et al., 2011;Boiteau et al., 2016). First, dissolved hydroxamate siderophores (300-1000 Da) are ubiquitous in oceanic systems. For example, in surface waters and water below the photic zone of the Central California upwelling, more than half of the Fe-binding compounds have been characterized as either hydroxamate type or catecholate type siderophores (Macrellis et al., 2001). In contrast, in the non-Fe-limiting regions (e.g., the North Atlantic Ocean), the characterizable siderophores exclusively consist of hydroxamate siderophores, which are produced by heterotrophic bacteria and fungi (Raines et al., 2015) in the concentration range of 3-20 pM (Mawji et al., 2008). In the eastern Pacific Ocean, two types of hydroxamate siderophores (hydrophilic siderophore Ferrioxamine B and amphiphilic amphibactins) dominate the siderophore pool and alternate in their abundance across a transect from coastal waters to the HNLC region to the oligotrophic ocean (Boiteau et al., 2016). Second, dissolve hydroxamate siderophores have the highest stability constants (as logK) that range from 16 to 62 at 25 • C for complexation of Fe(H 2 O) 6 3+ with the fully deprotonated ligand of hydroxamate functionality (N(OH)-COOH) (Boukhalfa and Crumbliss, 2002). The tri-hydroxamate siderophores have an average conditional stability constant for Fe of 30.6 (Kraemer et al., 2005;Ahmed and Holmstrom, 2014). Third, hydroxamate siderophores are hydrophilic or amphiphilic, and they have relatively higher resistance towards photolysis compared to other types of Fe-binding ligands (Barbeau et al., 2001(Barbeau et al., , 2003. These properties likely enhance their transport and stability in the upper ocean.
Most current studies have focused on the iron-carrying ligands in the dissolved phase (Macrellis et al., 2001;Hunter and Boyd, 2007;Vraspir and Butler, 2009;Bundy et al., 2016), and only recently it was suggested that hydroxamate-like moieties also existed in the particles/colloids (Chuang et al., 2013(Chuang et al., , 2015a). Yet the mechanisms by which these hydroxamatelike moieties are formed/incorporated in the particulate or colloidal phases, their molecular characteristics, their micromolecular environment, and whether these potential Fe-binding sites play a significant role in Fe biogeochemical cycling remain largely unknown. The reason lies in the analytical difficulty of qualitatively and quantitatively assessing these iron-carrying ligands in the particulate phase, especially those hydroxamatelike moieties, the latter of which are a minor component of the particulate organic carbon (POC) pool [determined as <1.5% of POC in Chuang et al. (2013Chuang et al. ( , 2015a], and thus "masked" when applying techniques like NMR, FTIR or C K-edge X-ray absorption spectroscopy, and/or gas and liquid chromatography. In the present study, electrospray ionization (ESI) Fouriertransform ion cyclotron resonance mass spectrometry (FTICRMS) was applied to elucidate the composition of organic matter from marine settling particles, using a sequential extraction method. Our objectives were to: (1) examine the overall molecular information of the settling particulate organic matter (POM) with depth, i.e., what the micromolecular environment is like for the iron-complexing moieties; (2) characterize the molecular features of iron-complexing sites in fractions of the marine particles that were probed by different extraction reagents, thus representing possible mechanisms for their incorporation into the settling particles; and (3) investigate their distributions at three depths (500, 1500, and 3200 m) in the northwestern Sargasso Sea, where there have been several studies in the particle and elemental fluxes with the Bermuda Time-Series Site (Conte et al., 2019).

MATERIALS AND METHODS
Here we apply an optimized sequential extraction method (Chuang et al., 2015b) followed by ESI-FTICRMS, to characterize the molecular features of the Fe-complexing sites in sinking particles (500, 1500, and 3200 m) in the northwestern Sargasso Sea near Bermuda. Using Fe(III) as the tracer, we identified potentially iron-complexing ligands in different operationallydefined fractions: an EDTA (Ethylenediaminetetraacetic acid disodium salt dihydrate) fraction of "surficial metal-complexed" POM (EDTA is used to chelate the surficial cations, dissociate the POM-metal complex, and thus the OM will be released into the solution); a SDS (sodium dodecyl sulfate) fraction of molecules that are bound to particles via hydrophobic interactions, i.e., "hydrophobically-sorbed POM;" and a HF (hydrofluoric acid) fraction of residual "interlayered and interstitial" mineral complexed/occluded POM.

Sample Collection
Particle flux material was collected by the Bermuda Oceanic Flux Program (OFP) time-series, located at 31 • 50 N, 64 • 10 W. The OFP mooring utilizes McLane Parflux sediment traps (0.5 m 2 surface area; McLane Labs, Falmouth, MA, United States) at 500, 1500, and 3200 m depths that are programmed at a nominal 2-week sampling interval. Sample cups are filled with high purity Sargasso Sea deep water (∼3000 m) brine (40 ppt) and amended with ultra-high purity HgCl 2 (200 mg/L) to minimize microbial degradation of the POM during sample collection. Details of sample collection and analytical methods are described elsewhere (Conte et al., 2019). Samples analyzed herein were collected in 2005 and 2006, and stored refrigerated. Additional sample information is provided in Supplementary Table S1.

Sequential Extraction of Marine Particles
The bulk hydroxamate concentration was determined using samples collected at different time points (n = 3 for each depth, see below section "Hydroxamate-Like Functionalities Determination in Marine Particles" and Supplementary  Table S1). Due to the limited sample material, the sequential extraction required pooling of three discrete samples collected at different time points of the same depth (Supplementary Table S1). The sequential extraction procedure is shown in Figure 1. Briefly, ∼10 mg of particles from each depth were first incubated with 1 mL 100 mM EDTA solution (pH 8.0) at 4 • C overnight on an end-to-end orbital shaker. The extractant soluble phase was separated from the particles via centrifugation and filtration with a syringe filter (PTFE membrane, diameter 25 mm, pore size 0.45 µm). The filtrate was extensively diafiltered (10-12 times) against ultrapure water with an Amicon Ultra-4 centrifugal filter unit with a membrane cutoff of 3 kDa. This retentate (>3 kDa and <0.45 µm) is defined as "EDTA extract." The residual particles were then re-suspended with 1 mL 1% SDS in 0.01 M Tris buffer (pH 6.8) at 95 • C for 1 h. This extractant was processed the same way as described for "EDTA extract" and defined as the "SDS extract." The residual pellet from SDS extraction was then washed three times with ultrapure water and subsequently treated with 1 mL 52% HF on ice for 1 h. The slurry was evaporated to dryness under a N 2 stream. After this, the residual particle material was re-suspended in 1 mL Tris-HCl (250 mM, pH 8.0) and centrifuged at 5000 × g for 15 min to separate the particles from the solution, the latter of which was processed as described previously by ultrafiltration to isolate the "HF extract" (Chuang et al., 2015b).
In order to better understand the stepwise extraction efficiency of POM from the marine particles, particles collected at different time points from the three depths (Supplementary Table S1) were well mixed (∼5 mg each) with a spatula and four replicate particles were extracted in parallel to obtain the particles from each sequential step for elemental analysis, using a Perkin-Elmer CHNS/O 2400 analyzer. The mass loss and organic carbon contents of the residual particles after each extraction step were determined and the analytical error of these four replicates were reported. For elemental analysis, the pre-weighed sample was directly acidified in a silver capsule three times and was heated at 60 • C overnight to remove any residual acid, prior to the analysis. Acetanilide (71.09% C, 6.71% H, 10.36% N) was used as the calibration standard.

Incubation of Fe (III) and Sequentially-Extracted POM Fractions
Sequential extractants, in the matrix of ultrapure water, from marine particles (section "Sequential Extraction of Marine Particles") were split into two aliquots, one amended with FeCl 3 (10 µM) and one without (control). The pH difference between these two treatments [i.e., Fe (+) and Fe (-)] were within 0.5 and both were around neutral. These were incubated at 4 • C on an end-to-end orbital shaker in the dark for 3 days. After that, excess FeCl 3 was removed by extensively diafiltering the sample against ultrapure water with an Amicon Ultra-4 centrifugal filter unit with a membrane cutoff of 3 kDa (Figure 1), and the retentate (>3 kDa) was kept for further ESI-FTICRMS analysis. The sample of the Fe (-) treatment was processed the same manner (with diafiltration) for comparison purpose.

ESI-FTICRMS Analysis
Marine particle extracts were diluted with methanol (1:1, v/v) immediately before injection into a Bruker Daltonics 12 Tesla Apex Qe ESI-FTICRMS (College of Sciences Major Instrumentation Cluster, COSMIC, Old Dominion University, Virginia). Positive ion mode was used, as the molecules of interest are iron-carrying molecules, which have better ionization efficiency in this mode (Waska et al., 2016). Samples were continuously injected by a syringe pump (infusion rate of 120 µL/h), and electrospray voltages were optimized for each sample. Ions were accumulated in a hexapole for 2.0 s before being transferred to the ICR cell and 300 transients were co-added. Fourier transformation and magnitude calculation of the free induction decay signal (FID) was accomplished using the Bruker Daltonics Data Analysis software (Version 4.0 SP 2). A total of 18 samples (three depths × three extractions × two treatments), as well as six procedure blanks (three extractions × two treatments) following the same steps as the samples (one for Fe addition and the other for no Fe addition), were included for FTICRMS analysis. Selected samples were analyzed multiple times to ensure the analytical reproducibility and all samples were analyzed within the same week to ensure stable instrumentation parameters and conditions for all samples (Sleighter et al., 2012).
Mass spectra were externally calibrated with a polyethylene glycol standard and then internally calibrated using a set of m/z values within the samples. Several criteria were met for picking out the calibration peak list: (1) peaks have a signal-to-noise (S/N) ratio ≥5, (2) peaks that have intermediate intensities to minimize instrumental artifacts that yield peaks with too high or too low intensities, (3) peaks have an error (i.e., the difference between the exact mass of calculated formula and the measured m/z) of <0.5 ppm, and (4) peaks can generate only one possible formula containing only C, H, O and/or Na elements, and all formulas have their corresponding structures found in Pubchem database. There was at least one calibration point every 15-20 m/z units (i.e., every other CH 2 ), across the entire m/z region (200-1400 m/z). Linear mode calibration was performed on each dataset.
All m/z lists, with an S/N ≥ 4, were first searched for 13 C peaks and confirmed they were mostly singly-charged. Those 13 C peaks were not included in the formula calculations. Moreover, the m/z values present in the procedure blank spectra were excluded from the sample mass list. A molecular formula calculator (Molecular Formula Calc Version 1.2.3 NHMFL, 2016) generated empirical formula matches within 1.0 ppm using elemental ranges of . All formulas were screened to eliminate those that are unlikely to occur in OM, according to a list of selection criteria (Neubert et al., 2002;Kind and Fiehn, 2007;Stubbins et al., 2010;Walker et al., 2016). Kendrick mass defect (KMD) analyses were then performed to look for homologous series of CH 2 , CH 2 O, COO, O, H 2 O, H 2 , etc., respectively, formed by at least two unique formulas (i.e., only one possible formula corresponds to a single m/z value after the formula rule's screening). If there are multiple possible formulas for a single m/z value (usually those with m/z values >400 Da), the one that fell within a homologous series was kept as the more likely formula (Sleighter and Hatcher, 2007). Subsequently, multiple formulas of one m/z value were screened by searching chemical structures in the Pubchem database. Formulas with no chemical structures identified were excluded from the final formula assignment. Fe-carrying molecules were also examined if their 54 Fe isotopic peak could be found (i.e., the singly-charged 56 Fe and 54 Fe-carrying peak m/z difference is 1.99533 ± 0.00070 and the isotopic ratio of 54 Fe/ 56 Fe as 0.063 ± 0.030) (Boiteau and Repeta, 2015;Walker et al., 2016). Adduct of iron hydroxide with methanol is likely irrelevant, since the mass list of the procedure blank (with the addition of FeCl 3 yet without the addition of particle extracts) had been subtracted Frontiers in Earth Science | www.frontiersin.org from the sample mass lists. Besides, this issue should occur much less likely at µM FeCl 3 concentrations (Waska et al., 2016).
To unambiguously identify a specific functional group, a tandem mass spectrometry analysis (MS/MS) would be required Hatcher, 2007, 2008). However, MS/MS analysis would have great difficulty to identify a trace component, i.e., a hydroxamate functional group (<1% of POC, see section "Bulk Hydroxamate Concentrations in Marine Particles") using the bulk OM without any extensive pre-separation and concentration. Instead, a KMD (Kendrick Mass Defect) analysis was conducted to search the whole mass list for molecules that fall in some homologous series of the key functionality of hydroxamate, N(OH)-COOH (with the formula total H number being corrected for Na adducts or Fe complexation). A homologous series of N(OH)-COOH was formed, if at least two formulas differ only by [N(OH)-COOH] n .

Hydroxamate-Like Functionalities Determination in Marine Particles
Bulk hydroxamate concentration of the particles was quantified using the classical "Csaky" spectrophotometric method (Gillam et al., 1981), which required a hydrolysis step to first release the bound hydroxamic acid. The un-hydrolyzed particle yielded no response to the reactions, suggesting it is hydroxamic acid rather than hydroxylamine that is being detected. We noticed a rapid "browning" effect of this method when analyzing particle samples, likely caused by sulfuric acid. This effect was corrected by subtracting the absorbance of a control sample without the addition of the final coloring reagent from that of a sample aliquot with the addition of the final coloring reagent. The net absorbance was used to calculate hydroxamate concentration from the calibration curve using acetohydroxamic acid (AHA) as the standard. This modified method was validated with siderophore standards added to marine particles or soil particles (Xu et al., 2015). In addition, we also critically scrutinized if there could be any other interferences by the following biomolecule classes: (1) neutral polysaccharides (e.g., dextran); (2) acidic polysaccharides and monosaccharides (e.g., alginic acid and glucuronic acid); (3) proteins/amino acid (e.g., cytochrome c and phenylalanine); (4) fatty acids (e.g., palmitic acid); (5) humic acid (e.g., 1S104H from IHSS); (6) fulvic acid (2S101F from IHSS); (7) deoxyribonuclease. We verified that none of these compounds caused any significant interferences to this method (i.e., their absorbance is the same as in ultrapure water) (Xu et al., data not shown).
The presence of N(OH)-COOH in the marine particles, either detected as "positive" with the "Csaky" method, or calculated by KMD analysis to fall in the homologous N(OH)-COOH series, are defined as "hydroxamate-like" molecules below.

Relative Distribution of POM in Sequential Extracts With Depth
The mass and organic carbon changes after each sequential extraction step are shown in Supplementary Table S2. The organic carbon (%, residual mass) generally increased with each extraction step, from the original value of (5.55% ± 0.05) to (45.48% ± 1.90) after HF extraction. EDTA extraction released (67% ± 2) of the particle mass and (39.74% ± 7.52) of the total original organic carbon into the solution. SDS and HF extraction released similar amounts of the original mass (13% ± 2) and (17% ± 3), yet quite different amounts of original organic carbon: (30.17% ± 2.35) and (7.46% ± 0.95), respectively. Their respective extraction efficiencies, which were calculated as the loss of particulate mass or organic carbon (the part being extracted into the solution) to those in the previous particles/residues being extracted, change in the following order: HF (86% ± 3) > EDTA (67% ± 2) > SDS (41% ± 7) by mass, and SDS: (49.92% ± 2.34) > EDTA (39.74% ± 7.52) > HF (24.68% ± 2.74) based on organic carbon. These rankings are likely to be different if one changes the sequential extraction order, which is not the focus of the current study. These three consecutive steps removed substantial amounts of original mass (98%) and of original organic carbon (77%) from the marine sinking particles, supporting the idea that ESI-FTICRMS analysis of these fractions represents the overall composition of OM in the marine particles.
The ESI-FTICRMS intensity percentages of the three extracts with depth is given in Figure 3. For each depth, the mass lists of the three extract (EDTA, SDS and HF) of the non-FeCl 3amended samples were combined. If there are common m/z values shared by different samples within the same depth, the average intensity was adopted. At 500 m, the POM consists mostly of SDS-extractable material (87%), followed by EDTAextractable compounds (12%), and the HF-extractable material  is only a very minor fraction (1%). At 1500 m depth, the SDS fraction decreases to 51% of the total extractable POM, while the EDTA and HF fractions increase (22 and 27%, respectively). At 3200 m, POM compounds are more evenly distributed and comprise 34, 46, and 20% for the EDTA, SDS and HF fractions, respectively.

Molecular Characteristics of Sequential Extracts of Marine Particles
Each sequential extractant (the non-FeCl 3 -amended) recovered a distinct set of compounds (in terms of elemental groupings and compound classifications) from the particles. The formula elemental groupings are listed in Table 1. CHO, CHON, and CHONS formulas are three major groups, all together accounting for 96, 91, and 85% of the total formulas for 500, 1500, and 3200 m, respectively. The number of peaks shared among POM at the different depths is shown in Supplementary Table S3 (773 formulas between 500 and 1500 m; 790 formulas between 1500 and 3200 m; 877 formulas between 500 and 3200 m). CHO formulas increased with depth, from 27% at 500 m to 35% at 1500 m, and to 41% at 3200 m. In contrast, the CHON formulas slightly decreased with depth, from 38% at 500 m to 39% at 1500 m, and to 34% at 3200 m. CHONS formulas decreased more rapidly with depth, from 31% at 500 m to 17% at 1500 m, and to 9% at 3200 m. CHOS, a minor fraction of the total identified formulas (within 5%) increased by twofold from 1.5% at 500 m to 3.7% at 1500 m, and stayed steady as 3.4% at 3200 m. Notably, deoxycompounds increased from 2% at 500 m to 5% at 1500 m, and to 12% at 3200 m.  The compound classifications are displayed in van Krevelen diagrams (Figure 4) and also listed in Table 2. From 500 to 3200 m, aliphatics steadily increased from 41.8 to 67.9%, while the N-containing aliphatics, which mostly consist of protein, peptides, and amino sugars, first increased from 15.4% at 500 m to 33.1% at 1500 m, and then decreased to 17.6% at 3200 m. Aromatics, condensed aromatics and black-carbon-like compounds decreased from 27.2 to 8.6%, 7.5 to 3.7%, and 3.1 to 0.7%, respectively. Similarly, with increasing depth the intensitynormalized H/C values (Sleighter and Hatcher, 2008) (calculated by summing up the product of individual intensity and the corresponding H/C ratio of each formula, and then being divided by the total intensity) increased from 1.31 to 1.62 whereas O/C values decreased from 0.20 to 0.16. CRAM-like compounds were only a minor fraction (5.0% at 500 m to 4.4% at 3200 m).

Iron-Carrying Molecules
In total, 745 iron-carrying molecules were identifiable, accounting for ∼14% of total assignable molecules including all depths, FeCl 3 -amended and non-FeCl 3 -amended treatments ( Figure 5A). The majority of the Fe-carrying molecules were found in the FeCl 3 -amended samples, and only 33 formulas were found in the non-FeCl 3 -amended samples. These ironcarrying molecules could have been corroborated by adding another trivalent metal (e.g., Ga) at the same concentration (10 µM) and identifying if the same apo molecule forms the similar complex with the other metal, which can be applied in future study.
The absolute number of iron-carrying molecules in the "EDTA extract" increased from 103 at 500 m to 587 at 1500 m and then dropped to 101 at 3200 m; those in the "SDS extract" increased from 31 at 500 m to 106 at 1500 m and then dropped to 37 at 3200 m; those in the "HF extract" increased from 35 at 500 m to 74 at 1500 m and to 99 at 3200 m ( Figure 5A). The EDTA fraction ranked highest in the number of iron-carrying molecules at both the 500 m and 1500 m depths, while at 3200 m depth the EDTA and HF fractions have a similar number of Fe-carrying molecules. CHON and CHONS were the dominant elemental groups and accounted for 86.3% of total assignable Fecarrying formulas (Figure 5B). CHO and CHOS were both minor elemental groups, and accounted for 7.5 and 4.4% of total assignable Fe-carrying formulas, respectively. The compound classification of the Fe-carrying formulas is distinct from that of the bulk POM in that these Fecarrying molecules consist mainly of aromatics (57.6% if using AI, and 61.7% if using modified AI), black carbon (14.8%), CRAM-like compounds (11.5%), and very minor amount of aliphatics (1.7%) (Table 3 and Figure 6). The intensity-normalized H/C and O/C ratios are 0.75 and 0.44 ( Figure 6), respectively, indicating an unsaturated/aromatic and highly-oxygenated composition.

Direct Evidence for the Presence of Hydroxamate-Like Molecules
Eleven formulas were found that belong to some N(OH)-COOH homologous series (Supplementary Table S4). Their m/z ranged from 431.98069 to 684.99753, mostly falling into the regions of aromatics and condensed aromatics (Figure 6). All these hydroxamate-like molecules accounted for 0.20 and 0.41% of all the assignable formulas in number and intensity, respectively, and 1.48 and 0.95% of the total Fe-carrying formulas in number and intensity, respectively ( Table 3). These values provide a lower limit of hydroxamate-like moieties because not all these molecules may necessarily fall into any N(OH)-COOH homologous series. In the CH 2 -Kendrick-Mass-Defect (KMD)-Kendrick-Mass (KM) plot, which categorizes the molecules according to their mass defect (or the elemental composition), Kendrick mass defect of the Fe-containing molecules shifted below that of their corresponding apo molecules, the latter of which were found in the Fe (-) treatment (Supplementary Figure S1).

Characteristics of the Molecular Micro-Environment of the Iron-Carrying Molecules in the Settling Particles
The changing molecular characteristics of the sinking particles with depth can be contrasted with trends observed for dissolved organic matter (DOM) (Chen et al., 2014b;Timko et al., 2015). First, higher average H/C values of DOM have been reported for the surface water of both North Pacific and Atlantic oceans, compared to those values of their respective deep water (Chen et al., 2014b;Timko et al., 2015), whereas the sinking particles exhibit an opposite trend with depth in our study (i.e., increasing H/C values with depth, Table 2). Second, the contents of aliphatic and carbohydrate-like compounds in DOM decrease with depth, likely due to the selective biodegradation of these molecules with increasing depth (Chen et al., 2014b). In samples collected at the same Sargasso Sea location, Timko et al. (2015) found higher concentrations of aliphatic compounds in the surface DOM pool than deep waters, while conversely those of aromatic/unsaturated and polyphenolic-like compounds were higher in deep waters. In contrast, we observe that particulate aliphatics (Figure 4 and Table 2) contents increased with depth. In addition, though carbohydrates are not well ionized in ESI (Shen and Perreault, 1998;Chen et al., 2011), spectrophotometric measurements by Chuang et al. (2013) provide evidence for the carbohydrate preservation in the sinking particles: 30% of OC at 500 m, 33% at 1500 m, and 26% at 3200 m, compared to the percentage of POC in particle mass as 11% at 500 m, 7% at 1500 m, and 5% at 3200 m. Third, higher percentages of CRAM-like formulas were observed in the intermediate and deep waters compared to surface water DOM in both North Pacific and Atlantic oceans, generally accounting for a substantial portion of the total mass spectrum intensity (40-65%) (Chen et al., 2014b). However, the particulate CRAM-like compounds are not a significant component of the POM (<6%) (Figure 4 and Table 2).
The results show that aliphatic compounds, and carbohydratelike molecules become more predominant in sinking particles with increasing depth while aromatics, condensed aromatics, black carbon-like molecules, and CRAM-like compounds are preferentially "removed." One hypothetical explanation for this pattern is the microbial community attached on particles actively modifies POM and thus produces labile material (Boeuf et al., 2019;Omand et al., 2020). The other alternative explanation is these compounds might be condensation products resulting from a Maillard-type reaction (Brandes et al., 2004) and thus become less biodegradable (Hedges et al., 2001). The differences in the molecular data for the three sequential extraction fractions provide some hints for understanding the molecular interactions within the organic phase of the sinking particles. Hydrophobic interactions (as probed by SDS extraction) appear to be the dominant mechanism for surface-active organic materials to associate with particles at mesopelagic depths (500 m) (Figure 3). With increasing depth there is, however, a greater importance of metal complexation and cation bridging. It could be due to the fact that more intact (dead or senescent) cells are present in the sinking particles at the 500 m than the other two deeper depths (1500 and 3200 m), and thus more cellular and membrane material became available for SDS extraction. However, at the other two depths, the cells could have mostly broken down, and the resulting OM either has been degraded or resulted in the metal complexed fraction. Lastly, more organic materials appear to partition into particle interstices (as probed by HF extraction) with depth (1500 and 3200 m vs. 500 m, Figure 3), which could provide a physical protective mechanism for organic molecules to become less accessible to microbial activity.

Refractory Nature of Particulate Hydroxamate-Like Moieties Revealed by Their Depth Distribution
In general, the absolute hydroxamate concentration in the settling particles decreased with depth in the Sargasso Sea while the POC-normalized concentration did not. This could indicate that hydroxamate-like moieties in the settling particles might be relatively less labile, compared to the bulk POC (Figures 2A,B). These hydroxamate concentrations are 10 to 20-fold lower than what was reported by Chuang et al. (2013). The source of this discrepancy is likely because Chuang et al. (2013) did not correct for the "browning effect" caused by the sulfuric acid step (see above).
The particulate hydroxamate concentrations determined at different seasons in the upper mesopelagic depth (500 m) has a larger variation compared to those in the intermediate (1500 m) and bathypelagic (3200 m) depths (Figure 2 and Supplementary  Table S1). This suggests that abundances of iron and microbiallyderived iron-binding ligands are more seasonally-driven in upper mesopelagic (e.g., seasonality of atmospheric dust input to surface waters) than at the intermediate and deep depths (Sedwick et al., 2005;Bundy et al., 2016). Once incorporated into the settling FIGURE 6 | Compound classification of Fe-carrying compounds from three extractants (blue circle, Fe-containing formulas; red square, Fe, hydroxamate-like containing formulas).
particles and transported to depth, hydroxamate-like molecules could have been stabilized via metal complexation (as probed by EDTA) (Iskrenova-Tchoukova et al., 2010), or hydrophobic sorption (as probed by SDS) or mineral interlayered occlusion (as probed by HF), and thus they may have been protected against enzymatic reaction (see below discussion).
While one needs to be cautious about applying the ESI-FTICRMS data for any quantitative purpose given the ionization efficiencies of ESI are variable for different types of compound groups, our results are in good agreement with the spectrophotometric results (a range of 0.2-0.6% of POC contributed by the carbon of N(OH)-COOH, in a comparison to 0.41% of total peak intensity contributed by hydroxamatelike molecules, Table 3 and Supplementary Table S1). Both results show that the N(OH)-COOH functionality is only a trace component compared to the bulk POM. These two independent lines of evidence collectively support the presence of N(OH)-COOH functional group in the sinking particles of the open ocean, at a level < 1% of the POC.
Using the lower limit of particle mass concentration and flux [10 µg/L and 10 mg/m 2 /d reported in Chuang et al. (2013)] and a median value of 10 µM-AHA eq./g-particle for particulate hydroxamate concentration, the particulate hydroxamate concentration (in overall volume) and flux at the studied site can be estimated to be 100 pM and 0.1 µmol-AHA eq./m 2 /d, respectively. This is on the same order of magnitude as the dissolved equivalent hydroxamate flux in the same area which is estimated to be 0.2 µmol-AHA eq. /m 2 /d, assuming a dissolved organic carbon flux of 15 mg/m 2 /d (Giering et al., 2014) and a hydroxamate-C to dissolved organic carbon ratio of 10 −6 . Meanwhile, Mawji et al. (2008) reported the dissolved hydroxamate siderophore in the surface water of the Atlantic Ocean was 3-20 pM. Thus, our results indicate that while hydroxamate-like moieties appear to be a trace component of POM, they could be comparable to their dissolved counterparts in the concentrations and fluxes. However, one has to be cautious that a direct comparison between the hydroxamate-like moieties in the particulate phase and the dissolved hydroxamate siderophores deserves more scrutiny. The reason is that there are still quantitative uncertainties based on the assumptions of relative concentrations of hydroxamate moieties and organic carbon in particles and water, and lack of evidence to connect dissolved hydroxamate siderophore as the source of the hydroxamate-like moieties in the particulate phase.

Possible Mechanisms for Incorporation and Preservation of Iron-Carrying Molecules and the Subset Hydroxamate-Like Molecules
We observed distinct molecular differences between Fecarrying molecules with likely Fe-chelating properties and the bulk POM. This could reflect differences in biogeochemical sources and respective biosynthetic pathways. For example, Chen et al. (2014a) showed that POM produced from photoflocculated DOM was enriched with both aliphatic and aromatic components, yet no direct evidence supported whether the co-precipitated Fe was complexed with both types of compounds or preferentially complexed with one or the other. This Fe-rich photo-produced black carbon pool (∼10% of the POC) was estimated to be nearly equivalent to the estimated global flux of dissolved black carbon from land to the ocean (Chen et al., 2014a). In a study that was not conducted in the oceanic context, plutonium (an A-type or hard metal (i.e., ion with low polarizability) that was often referred to having similar chemical complexing properties as Fe due to their similar ionic potential Z/r) (Chuang et al., 2013) was found to be predominantly associated with nitrogen-enriched, low average H/C (∼0.64), high DBE (∼16), aromatic and condensed aromatic compounds (DiDonato et al., 2017). In the present study, we did observe the Fe-carrying molecules are enriched in aromatic and condensed aromatics, and also contain black carbon with a content similar to what was reported in Chen et al. (2014a), however, these particulate Fe-carrying molecules were deplete in aliphatic (only 1.7%, Table 3), and also distinctly different from the bulk POM (41.8-67.9%, Table 2). A hypothetical explanation is as follows: Below 500 m where photochemical reactions are minor, the aliphatic moieties of these molecules are preferentially removed via microbial activity. Meanwhile, both the black carbon-like and CRAM-like compounds are preserved, compared to the bulk of the POM, and the latter exhibit a totally opposite pattern (i.e., aliphatic-enriched, aromatic and CRAM-depleted and deoxygenation) with increasing depth (Table 2 vs. Table 3).
While the hydroxamate-like moieties are only a minor fraction of the potential Fe-binding ligands of the POM, they may play an important role in both carbon and Fe cycling. As shown here, a substantial portion of the POM (∼77% of POC) can be removed by our sequential extraction method and saturated with dissolved Fe 3+ , which would otherwise not likely be exposed and accessible to Fe 3+ under natural conditions. Furthermore, the flux of the hydroxamate-like moieties in the sinking particles is the same order of magnitude as their dissolved counterparts. We adopt here the hydroxamate flux of 0.1 µmol-AHA eq./m 2 /d as we calculated in section "Refractory Nature of Particulate Hydroxamate-Like Moieties Revealed by Their Depth Distribution." Considering hydroxamate and Fe usually forms hexadentate or octahedral complexes via three bidentate ligands, this particulate hydroxamate flux has a theoretical maximum capacity of carrying an Fe flux of 0.03 µmol-Fe/m 2 /d. This is to be compared to an average particulate Fe flux of 4 µmol-Fe/m 2 /d (Chuang et al., 2013). Uncertainties in this assessment are possibly large, yet these rough calculations suggest that particulate hydroxamate could be one of the strongest organic Fe-binding ligands in natural sinking particles, considering their high binding affinity to Fe(III) over other known ligands (Reid et al., 1993;Neilands, 1995;Winkelmann, 2002;Butler, 2005;Butler and Theisen, 2010).
In addition to the above mentioned general molecular features of the iron-carrying molecules (Figure 6 and Table 3), the prevalence of sulfur with these hydroxamatelike functionalities (Supplementary Table S4) is noticeable, as dissolved S-containing hydroxamates were previously reported FIGURE 7 | Conceptual scheme showing the possible mechanisms for hydroxamate-like moieties incorporation into marine settling particle.
Frontiers in Earth Science | www.frontiersin.org yet in a low proportion to the total dissolved hydroxamate siderophores in natural environments (Gledhill et al., 2019). How sulfur was incorporated into these hydroxamate-like moieties in the marine particles is equivocal, with one option being via abiotic pathways (Sleighter and Hatcher, 2008;Raven et al., 2016;Pohlabeln et al., 2017), or biological processes (Boeuf et al., 2019). Another molecular feature is that these hydroxamate-like molecules have relatively low H/C ratio (average H/C as 0.83, Supplementary Table S4) compared to those of the dissolved hydroxamate siderophores which usually exceed 1.0 (e.g., average H/C in the dissolved hydroxamate siderophore is 1.49; Gledhill et al., 2019). It is possible that dissolved hydroxamate siderophores produced by particle-attached bacteria or fungi are a direct source to be incorporated into the particles. However, these molecular characteristics suggest that these hydroxamate-like moieties are by no means similar in chemical composition or structure to their dissolved counterparts, thus they likely have experienced complex biogeochemical processes yet to be defined (Figure 7).
First, hydroxamate-like moieties might be incorporated into the particles via sunlight-induced reactive oxygen species (ROS) mediated protein aggregation (Sun et al., 2017(Sun et al., , 2018(Sun et al., , 2019 in surface waters. In support of this mechanism, many previous studies have shown, although indirectly, that more nitrogen than carbon could accumulate into particles during irradiation (Helms et al., 2013;Chen et al., 2014a;Sun et al., 2019). We have observed N-containing aliphatics, consisting mostly of proteins, peptides, and amino sugars (Chen et al., 2018), accounting for 15.4-33.1% of the POM (in intensity %, Table 2). In comparison Chuang et al. (2013) found that about 45% of POC in particles consisted of protein (Chuang et al., 2013). Thus, it is compelling to believe that proteinaceous components are a significant component of the settling particles. The photo-reaction-induced and ROSmediated crosslinking process in surface waters would generate precursors including aliphatic and condensed aromatic-like molecules (Chen et al., 2014a) that would be further modified as particles are processed in the water column (see above in this section).
A protein encapsulating mechanism could be an alternative pathway for hydroxamate-like moieties incorporation into marine particles. For example, ferritin, a large protein primarily storing iron mineral internally (Frankel, 1991) by marine microbes (Marchetti et al., 2009;Shire and Kustka, 2015), has been demonstrated to also encapsulate siderophores into cavitylike structures within the apoferritin (Dominguez-Vera, 2004). This ferritin-encapsulating mechanism has been exploited in health sciences to design slow-release drugs (Dominguez-Vera, 2004;Zhang et al., 2006;He and Marles-Wright, 2015).
Finally, some labile organic molecules, e.g., amino acids or simple sugars, have been shown to be preserved in marine sediments via sorption mechanisms (Keil et al., 1994;Barber et al., 2017). This likely involves relatively "weak" forces such as hydrophobic, ionic interactions or hydrogen bonding that are associated with mineral phases. Some bacteria (e.g., Pseudomonas aeruginosa) have been shown to utilize pyoverdine (a siderophore that contains both catecholateand hydroxamate-like ligands) to complex the metal oxide surface for attachment (Upritchard et al., 2007). Thus, it is possible that some hydroxamate-like moieties in particles might be located with bacteria. In addition, the apparent LMW hydroxamate-like moieties observed in ESI-FTICRMS (Supplementary Table S4) could also be ionically bound as they were reversibly desorbed by the EDTA reagent. LMW molecules that are held together via hydrophobic interactions, hydrogen bonding and ionic interactions (complexation) can form supramolecular ("physical cross-linking") associations (Sutton and Sposito, 2005) with an apparent molecular weight over 3 kDa (either present in the marine particles captured by sediment trap or a 0.45 µm membrane, or an ultrafiltration membrane). These supramolecular associations can dissociate within the electrospray interface, releasing LMW molecules.
Regardless of the exact mechanism of incorporation, it is clearly evident from our work that hydroxamate-like moieties exist in the particulate phase, potentially equivalent in concentrations and fluxes to their counterparts in the dissolved phase. The Fe-carrying compounds which represent all the potential organic Fe-binding ligands of the settling particles, have molecular features in a contrast to the bulk POM, the latter of which is then different from the vast DOM pool. The next step would be to identify the exact mechanism of hydroxamate-like molecule formation/incorporation, as this would be of great value for better understanding of the cycling of Fe, organic matter and organic nitrogen in the ocean. While our study only analyzed particulate Fe-carrying molecules in the oligotrophic North Atlantic gyre, our results strongly suggest the relative importance of hydroxamate-like molecules, which possibly have been incorporated via various physical and chemical mechanisms. These, along with the major fractions of Fe-binding compounds identified here, likely play an important role in controlling Fe and other A-type particle-reactive radionuclides (e.g., thorium) distributions and flux in the ocean.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/Supplementary Material.

AUTHOR CONTRIBUTIONS
CX and PS conceived and designed the experiment. CX and WX performed the particle extraction experiment. CX performed the iron amendment incubation experiment. PL, CX, and HC performed the ESI-FTICRMS experiment. CX, PL, LS, HC, WX, and MK performed the ESI-FTICRMS data analysis. CX, PL, HC, WX, MK, PH, MC, AQ, and PS wrote the manuscript. All authors contributed to the article and approved the submitted version.