Molecular and Microscopic Insights into the Formation of Soil Organic Matter in a Red Pine Rhizosphere

Microbially-derived carbon inputs to soils play an important role in forming soil organic matter (SOM), but detailed knowledge of basic mechanisms of carbon (C) cycling, such as stabilization of organic C compounds originating from rhizodeposition, is scarce. This study aimed to investigate the stability of rhizosphere-produced carbon components in a model laboratory mesocosm of Pinus resinosa grown in a designed mineral soil mix with limited nutrients. We utilized a suite of advanced imaging and molecular techniques to obtain a molecular-level identification of newly-formed SOM compounds, and considered implications regarding their degree of long-term persistence. The microbes in this controlled, nutrient-limited system, without pre-existing organic matter, produced extracellular polymeric substances that formed associations with nutrient-bearing minerals and contributed to the microbial mineral weathering process. Electron microscopy revealed unique ultrastructural residual signatures of biogenic C compounds, and the increased presence of an amorphous organic phase associated with the mineral phase was evidenced by X-ray diffraction. These findings provide insight into the formation of SOM products in ecosystems, and show that the plantand microbially-derived material associated with mineral matrices may be important components in current soil carbon models.


Introduction
Soils play a critical role in global carbon cycling, having a dynamic carbon reservoir in the form of soil organic matter (SOM) that is estimated to be four times larger than the amount stored in the atmosphere [1,2].In 2010, a new conceptual framework for a microbial carbon pump was proposed for the microbial production of recalcitrant non-living organic matter and the storage of fixed C in oceans [3].The related idea of a C pump that operates in terrestrial systems was subsequently proposed in 2011 [4,5].Detailed investigations of a microbial carbon pump in soils presents an area of several knowledge gaps, and a benchmark for upcoming research.Recently, new research on microbial carbon stabilization in soils has emerged, focusing on SOM molecular structure characterization, recalcitrance, and persistence [6][7][8], as well as studies showing the importance of microbially-derived C for stabilized SOM [9] and a conceptual framework for SOM stabilization [9,10].Understanding the dynamics of C pools in soil systems is critical for controlling atmospheric carbon dioxide levels and maintaining soil health and agricultural productivity.Microbial contributions to long-lived soil C pools has often have been regarded as low to negligible, partly because the carbon in living microbial biomass accounts for less than 4% of soil organic carbon [11].However, microorganisms add carbon to soil through rapid and persistent cell division, community growth, exudation of biogenic products, and cell death [12].More recently, there has been a greater recognition that microbial necromass, i.e., dead microbial biomass, may contribute to SOM pools that have longer turnover times [13][14][15].Microbially-derived C inputs to soils are now increasingly recognized as playing a far greater role in sequestration of SOM than previously believed [16].A new concept has been introduced, contrasting the traditional view of recalcitrant humic substances against the emergent paradigm with a "continuum of progressively decomposing organic compounds" [17].The challenges of characterizing the SOM pool need to be addressed through the development of high spatial, chemical, and molecular resolution techniques.A comprehensive review by Oburger and Schmidt [18] discusses approaches for research on C input in soils.Mikutta et al. [19] and Pett-Ridge and Firestone [20] present stable isotope-enabled capabilities for tracing dynamic processes of organic molecules from roots and microorganisms in SOM.Employing multimodal analytical approaches that combine physical, chemical, and biological techniques may bring us closer to addressing these research questions.
The rhizosphere, the zone immediately adjacent to plant roots, represents a hotspot of carbon exudation, turnover, respiration, transformation, and accumulation [21].Depending on the plant, 5 to 21% of atmospheric C that is photosynthetically fixed is transformed into root exudates [22,23].The importance of roots for the formation of SOM is discussed in Rasse et al. [24], with the root contribution to soil C pools as a key factor of SOM stabilization.These rhizospheric C compounds consist of easily degradable sugars and organic acids, as well as more recalcitrant materials including lignified cellulose and hemicellulose from sloughed-off root cells [25].These plant-specific growth substrates readily recruit rhizospheric soil microorganisms, which subsequently metabolize and decompose these and other rhizodeposits high in labile C. Additionally, microbially-secreted extracellular polymeric substances (EPS) play a dynamic role in numerous soil processes, including soil C stabilization on mineral surfaces and mineral weathering [26,27].These rhizospheric compounds released by plants and microbes contribute a substantial amount of organic C to soils.Depending on the susceptibility of these rhizospheric compounds to degradation and mineralization, they may become stabilized, and a part of the newly formed refractory SOM of plant and microbial origin (Figure 1).[28], Kuzyakov [21], and von Lützow [29].
The objective of this research was to gain insight into the role soil microorganisms play in C stabilization.In addition, the research aimed to develop a molecular understanding of newly formed SOM macromolecules produced by root exudates and microbes in the rhizosphere of red pine seedlings, in a controlled conditions experiment without any previous humic substances in the designed mineral matrix.The system was maintained under base cation-limited conditions in order to drive utilization of minerals in the substrate as cations source.Laboratory and field experiments have shown that plants and their microbial affiliates can extract and utilize base cations such as K and Ca from soil minerals when exposed to soil solutes lacking these cations resulting in increased mineral weathering [30][31][32][33].We anticipated that (a) microbes would contribute to the new SOM pool by stabilizing root exudates and producing complex refractory C compounds.With our recent development of a solvent-based method for molecular SOM characterization by high resolution mass spectrometry [34], we also expected that (b) microbial residues will leave a stable signature of persistent OM on mineral surfaces that remain attached following a series of extractions with different polarity solvents.Furthermore, we hypothesized that (c) SOM composition would differ between rhizosphere soil and bulk soil.To achieve these objectives, we implemented a collective approach of imaging and molecular methods (Figure 2) allowing us to characterize the nature and extent of refractory microbial material in the SOM.We categorized the organic C compounds based on elemental composition, biochemical classes, and lability and persistence.[28], Kuzyakov [21], and von Lützow [29].
The objective of this research was to gain insight into the role soil microorganisms play in C stabilization.In addition, the research aimed to develop a molecular understanding of newly formed SOM macromolecules produced by root exudates and microbes in the rhizosphere of red pine seedlings, in a controlled conditions experiment without any previous humic substances in the designed mineral matrix.The system was maintained under base cation-limited conditions in order to drive utilization of minerals in the substrate as cations source.Laboratory and field experiments have shown that plants and their microbial affiliates can extract and utilize base cations such as K and Ca from soil minerals when exposed to soil solutes lacking these cations resulting in increased mineral weathering [30][31][32][33].We anticipated that (a) microbes would contribute to the new SOM pool by stabilizing root exudates and producing complex refractory C compounds.With our recent development of a solvent-based method for molecular SOM characterization by high resolution mass spectrometry [34], we also expected that (b) microbial residues will leave a stable signature of persistent OM on mineral surfaces that remain attached following a series of extractions with different polarity solvents.Furthermore, we hypothesized that (c) SOM composition would differ between rhizosphere soil and bulk soil.To achieve these objectives, we implemented a collective approach of imaging and molecular methods (Figure 2) allowing us to characterize the nature and extent of refractory microbial material in the SOM.We categorized the organic C compounds based on elemental composition, biochemical classes, and lability and persistence.

Experimental Approach
Thirty six seedlings of red pine (Pinus resinosa), a pine species native to North America, were planted into 164 mL Ray Leach plastic cone-tainers TM (Steuwe and Sons, Tangent, OR, USA) that had been closed at the bottom with 60 μm nylon mesh [32] and filled with 170 g of sterilized quartz sand amended with 1.5% biotite and 1% anorthite by weight (Figure 2, 1. Tree sample in a column).These minerals were the seedlings' sole source of Ca and K. See Shi et al. [32] for more details, and Supplementary Information for the sterile medium mineral constituency, growth conditions, and microbial inoculum.To establish a rhizosphere microbial community comparable to those occurring under pine forest conditions, an inoculum consisting of a microbial suspension was prepared by washing roots of young white pine trees collected about 20 miles North-East of Potlatch, ID (46.9°N, 116.9°W).Time zero for the start of the column experiment was defined when 5 mL of inoculum was added to each mesocosm.The trees were provided with 20 mL of distilled water three times per week.After nine months, the columns were destructively assessed using the suite of analytical and imaging methods listed below (Figure 2, 2. Mechanical separation).Twelve biological replicates (trees) were used for the FTICR and ion microscopy, and 4 separate trees with a total of 30 technical replicates for microbial DNA extraction.The material for XRD was combined from the above samples after the respective analysis or treatment (extraction).Briefly, seedlings were harvested by gently rolling the cone-tainer to loosen the roots, and the seedlings and the soil core were placed on a clean surface.To isolate the rhizosphere from the bulk soil, the bulk substrate (BS) was first gently shaken off and collected, followed by more vigorous shaking to free the soil that remained loosely associated with

Experimental Approach
Thirty six seedlings of red pine (Pinus resinosa), a pine species native to North America, were planted into 164 mL Ray Leach plastic cone-tainers TM (Steuwe and Sons, Tangent, OR, USA) that had been closed at the bottom with 60 µm nylon mesh [32] and filled with 170 g of sterilized quartz sand amended with 1.5% biotite and 1% anorthite by weight (Figure 2, 1. Tree sample in a column).These minerals were the seedlings' sole source of Ca and K. See Shi et al. [32] for more details, and Supplementary Information for the sterile medium mineral constituency, growth conditions, and microbial inoculum.To establish a rhizosphere microbial community comparable to those occurring under pine forest conditions, an inoculum consisting of a microbial suspension was prepared by washing roots of young white pine trees collected about 20 miles North-East of Potlatch, ID (46.9 • N, 116.9 • W).Time zero for the start of the column experiment was defined when 5 mL of inoculum was added to each mesocosm.The trees were provided with 20 mL of distilled water three times per week.After nine months, the columns were destructively assessed using the suite of analytical and imaging methods listed below (Figure 2, 2. Mechanical separation).Twelve biological replicates (trees) were used for the FTICR and ion microscopy, and 4 separate trees with a total of 30 technical replicates for microbial DNA extraction.The material for XRD was combined from the above samples after the respective analysis or treatment (extraction).Briefly, seedlings were harvested by gently rolling the cone-tainer to loosen the roots, and the seedlings and the soil core were placed on a clean surface.To isolate the rhizosphere from the bulk soil, the bulk substrate (BS) was first gently shaken off and collected, followed by more vigorous shaking to free the soil that remained loosely associated with the root.The material that remained closely adhered to the root was defined as the rhizospheric material (RS).

Microbial DNA Extraction, Sequencing, and Data Analysis
Microbial diversity and community composition were determined by 16S rRNA amplicon sequencing and subsequent taxonomic analysis.Four biological replicates of red pine seedlings were sampled at 9 months after inoculation, each in two respective profiles: bulk substrate (BS) and rhizospheric material (RS) (Figure 2, 3C.DNA extraction).However, only data from the RS was used in this study due to a lack of amplifiable material extracted from the BS.Microbial DNA was extracted using a MoBio PowerSoil ® DNA Isolation Kit (MoBio Laboratories, Inc., Carlsbad, CA, USA) following the manufacturer's instructions.Eluted DNA was quantified with a QuBit ® 2.0 Fluorometer and QuBit ® dsDNA HS (high sensitivity) Assay Kit (Invitrogen, Carlsbad, CA, USA).The quality of the DNA product was checked by PCR amplification of the V4 region (515f-806r) of the 16S rRNA gene by using the protocol developed by the Earth Microbiome Project (http://press.igsb.anl.gov/earthmicrobiome/emp-standard-protocols/16s/), with a Taq PCR Master Mix Kit (Qiagen, Hilden, GE, USA) and a DNA Engine Tetrad 2 Peltier thermocycler (Bio Rad, Hercules, CA, USA).PCR conditions and details pertaining to the primer set are described in Caporaso et al [35].Eluted DNA was stored at −20 • C and subsamples (20 µL) from each biological replicate were sent to Argonne National Laboratory (ANL) for sequencing on the Illumina HiSeq2000.
Reads were demultiplexed, joined, and quality filtered using ea-utils [36], then analyzed with QIIME [37].VSEARCH [38], an optimal, open-source implementation of USEARCH, was used to dereplicate, sort by abundance, remove single reads, and cluster at 97% identity.The UCHIME algorithm as implemented in VSEARCH was used to check these clusters for chimeras and construct an abundance table by mapping labeled reads to chimera-checked clusters [39].Taxonomy was assigned to the centroid of each operational taxonomic unit (OTU) by finding the last common ancestor of the top three hits to the Greengenes database [40].The centroids were aligned to the database and a phylogenetic tree was constructed using FastTree [41].

Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FTICR-MS)
Biological replicates of red pine seedlings were sampled at 9 months after inoculation.After the mechanical separation of the material from the columns (Figure 2, 2. Mechanical separation), samples of bulk substrate (BS) and rhizospheric material (RS) were extracted in either water (Optima ® LC/MS grade, Fisher Scientific, Suwanee, GA, USA) or 100% methanol (Optima ® HPLC grade, Fisher Scientific, Suwanee, GA, USA) according to Tfaily et al. [34] as follows: 0.5 g of BS or RS were weighed and divided evenly for water and methanol extraction, then mixed and shaken for 5 minutes.For RS, root samples were dipped in the solvent, manually stirred and gently shaken for 5 minutes, to minimize the portions of cut roots to be extracted in water or methanol.(Figure 2, 3A.H 2 O, MeOH extraction).Water was the solvent of choice to extract labile compounds such as sugars and amino acids, in addition to polar lignin-like and aromatic compounds.Methanol was used as a complementary solvent to water and was used to extract compounds with low O:C elemental ratios in addition to polar lipid-like compounds.After the extraction, the material was spun down, and solvents with extracted material were collected for FTICR-MS analysis and TEM imaging.To obtain high-resolution FTICR mass spectra of the organic matter in the extracts, the analysis was performed using a 12 Tesla 12T-FTICR-MS (Bruker SolariX, Billerica, MA, USA) outfitted with standard electrospray ionization (ESI) interface.The spectra were collected in negative ion mode.Water-extracted samples were diluted 1:2 with methanol and methanol extracted samples were run directly and infused using a 250 µL Hamilton syringe at a flow rate of 3 µL/min.The coated glass capillary temperature was set to 180 • C and the electrospray voltage set to +4.4 kV.The ESI signal was allowed to stabilize for 5 min and data were then collected from 100 to 1300 m/z, resolution was 4.5 × 10 5 at 451 Da, ion accumulation time was 0.1 s for 96 scan averages co-added, time of flight was set to 0.65 ms, and Q1 was set to 100 m/z.The syringe was then flushed with 50/50 methanol/water (v/v) to prepare for the next sample.Spectra were calibrated by two internal series of DOM homologous series separated by 14 Da (-CH 2 groups) and the mass accuracy was calculated to be <1 ppm for singly charged ions ranging across the mass spectra distribution.Molecular formula assignments were made using a modified version of the Compound Identification Algorithm (CIA) described by Kujawinski et al. [42] and updated as described subsequently [43].Biochemical compound classes are reported as relative abundance values based on counts of C, H, and O for the following H:C and O:C ranges; lipids (0 and condensed aromatics (0.5 ≤ O:C ≤ 0.9, 0.2 ≤ H:C < 0.8) ( [34], modified from Kim et al. [44]).Spectra from water-extracted samples were selective for carbohydrates with high O:C ratios, whereas spectra from methanol-extracted samples were selective for compounds with low O:C rations (O:C < 0.5) [34].The van Krevelen diagram is constructed using the elemental oxygen-to-carbon ratio (O/C ratio) on the x-axis and hydrogen to carbon (H/C ratio) on the y-axis.Major biogeochemical classes of compounds (such as lignin-like compounds, lipid-like compounds etc.) have their own characteristic H/C or O/C ratios.As a result, each class of compounds plots in a specific location on the diagram [44].

Transmission Electron Microscopy (TEM)
A whole mount approach was applied for imaging the material extracted from rhizospheric soils using the same solvents that were utilized for the FTICR-MS analysis (Figure 2, TEM imaging).A 5-µL drop of water or methanol with extracted material was applied to a 100-mesh Cu grid covered with formvar support film sputtered with carbon (Electron Microscopy Sciences, Hatfield, PA, USA).The material was allowed to adhere to the grid for 1 min before the liquid was gently blotted with a filter paper, and the material was negatively stained with a 5-uL drop of Nano-W (Nanoprobes, Yaphank, NY, USA).After 30 s, the excess liquid was removed by wicking and the sample was allowed to air dry.Samples were examined with a Tecnai T-12 TEM (FEI) operating at 120 kV with a LaB6 filament.Images were collected digitally with a 2 × 2 K UltraScan CCD (Gatan, Pleasanton, CA, USA).Representative images of at least 50 areas per sample were collected at 10,000× magnification and analyzed to distinguish the representative material ultrastructure.

Scanning Ion Microscopy-Helium Ion Microscopy (HeIM)
Mineral material from the same samples as for the FTICR was air-dried, then mounted on carbon tape-covered aluminum SEM stubs (Ted Pella, Redding, CA, USA), and sputter-coated with carbon.The mineral particles were imaged with a high-resolution Orion Helium ion microscope (HeIM) (Zeiss, Peabody, MA, USA) at 30 kV (Figure 2, 3B.Dehydration; SEM imaging).

XRD, Micro-XRD
Solid material from the same samples as for the FTICR and HeIM, (bulk substrate and rhizosphere, 0.5g each), were homogenized into fine powder by dry grinding, and loaded into 0.5-mm diameter quartz capillaries (Charles Supper Co., Natick, MA, USA) for micro-XRD measurements (Figure 2, 3A.H 2 O, MeOH extraction, Solid phase, Micro-XRD; and 3B.Grinding, Micro-XRD).Data were collected using a Rigaku D/Max Rapid II instrument equipped with a 2D image plate detector (The Woodlands, TX, USA).X-rays were generated with a MicroMax 007HF generator fitted with a rotating Cr anode (Lambda = 2.2897 Å), and focused on the specimen through a 300 µm diameter collimator.2DP, a Rigaku 2D Data Processing Software (Ver.1.0, Rigaku, 2007) was used to integrate the diffraction rings captured by the detector.Regular powder diffraction data were collected on a Phillips X'Pert Multipurpose X-ray diffractometer equipped with a fixed Cu anode operating at 45 kV and 40 mA.Analysis of the diffraction data for phase matching was carried out using JADE 9.5 (Materials Data, Inc., Livermore, CA, USA) and the PDF4 + 2014 database from ICDD.For the determination of amorphous content, a known amount of high-crystallinity corundum internal standard (a-Al2O3, NIST Standard Reference Material 674a) was mixed into the samples, and the amorphous content of the samples was calculated using formula [45]: Real amorphous fraction (Wi) as implemented in the program TOPAS v4.2 (Bruker AXS GmbH, Karlsruhe, Germany), where Ws is the experimental weight fraction of internal standard (corundum), and Wsc is the calculated weight fraction of internal standard from the refinement.

Column Growth Experiment
In nutrient-limited conditions, and in the absence of pre-existing humic soil components in the system, the experimental mesocosm sustained a diverse microbial community and produced a variety of organic compounds.The mineral matrix was initially sterile, thus all detected SOM was created by tree root exudates and microbial metabolism in situ.The plants exhibited stunted growth after 9 months of dependence on nutrients (Ca and K) released from minerals by the root system and its microbial affiliates.There was no leaf litter influence because the young trees did not undergo senescence.

Microbial Community Composition
Microbial community operational taxonomic units (OTUs) were solely observed in the rhizospheric material; the bulk substrate lacked an extracted signal that could be amplified.In the rhizospheric soil, the most common families of microorganisms accounted for >95% of all observations (Figure 3A).These represented ten bacterial phyla, with Proteobacteria being the most abundant, comprising 73% of the community.The most common orders of Proteobacteria observed were Burkholderiales, Rhizobiales, Pseudomonacterdales, and Nitrosomonadales (Figure 3A).We also observed abundant Saprospirales, Sphingobacteriales (from phylum Bacteriodetes), Actinomycetales (from phylum Actinobacteria), and Solibacterales (from phylum Acidobacteria) in RS.These taxa were present in rhizospheric soil of all four sampled columns with no significant variations of microbial community composition across biological and technical replicates.Family relatedness is shown using an inferred phylogenetic tree (Figure 3B).This phylogenetic tree shows that in the rhizospheric soil, Actinobacteria were closely related to each other and differed from the majority of other Proteobacteria, which also tended to be grouped together.The Rhizobiales were also clustered together.Out of the OTUs identified as being in the order Rhizobiales, the most abundant families included Bradyrhizobiaceae (5.4%), Phyllobacteriaceae (4.2%), and Hyphomicrobiaceae (4.1%) (Figure 4).

FTICR-MS
We characterized the chemical complexity and composition of the newly formed carbon compounds.Figure 5 summarizes these data from one of the biological replicates in a van Krevelen diagram which is a plot of the elemental ratios of O:C on the x-axis and H:C on the y-axis.All column samples showed a substantial presence of biochemical compounds extracted by water and methanol (Figure 5), especially in the RS.Additionally, each solvent was selective for the type of compounds that it extracted.In the water-extracted samples, a decrease in the % abundance of lignin-like, proteinlike, and amino sugar-like compounds was observed in the BS (Figure 6) compared to RS. Lignin-like compounds are of plant origin compared to protein-like and amino sugar compounds that could be a mix of both plant and microbial origin (Figure 1).The bulk substrate sample was mainly enriched in aromatic compounds with little to no labile compounds.In the methanol-extracted samples, both

FTICR-MS
We characterized the chemical complexity and composition of the newly formed carbon compounds.Figure 5 summarizes these data from one of the biological replicates in a van Krevelen diagram which is a plot of the elemental ratios of O:C on the x-axis and H:C on the y-axis.All column samples showed a substantial presence of biochemical compounds extracted by water and methanol (Figure 5), especially in the RS.Additionally, each solvent was selective for the type of compounds that it extracted.In the water-extracted samples, a decrease in the % abundance of lignin-like, proteinlike, and amino sugar-like compounds was observed in the BS (Figure 6) compared to RS. Lignin-like compounds are of plant origin compared to protein-like and amino sugar compounds that could be a mix of both plant and microbial origin (Figure 1).The bulk substrate sample was mainly enriched in aromatic compounds with little to no labile compounds.In the methanol-extracted samples, both

FTICR-MS
We characterized the chemical complexity and composition of the newly formed carbon compounds.Figure 5 summarizes these data from one of the biological replicates in a van Krevelen diagram which is a plot of the elemental ratios of O:C on the x-axis and H:C on the y-axis.All column samples showed a substantial presence of biochemical compounds extracted by water and methanol (Figure 5), especially in the RS.Additionally, each solvent was selective for the type of compounds that it extracted.In the water-extracted samples, a decrease in the % abundance of lignin-like, protein-like, and amino sugar-like compounds was observed in the BS (Figure 6) compared to RS. Lignin-like compounds are of plant origin compared to protein-like and amino sugar compounds that could be a mix of both plant and microbial origin (Figure 1).The bulk substrate sample was mainly enriched in aromatic compounds with little to no labile compounds.In the methanol-extracted samples, both the BS and RS appeared to be rich in lipid-like compounds that could be of plant and microbial origin, even though a higher abundance of these compounds was observed in RS samples.the BS and RS appeared to be rich in lipid-like compounds that could be of plant and microbial origin, even though a higher abundance of these compounds was observed in RS samples.

TEM
The high resolution TEM imaging revealed a clear difference in appearance in the water-and methanol-extracted material (Figure 7), consistent with FTICR-MS data.The material dissolved in the BS and RS appeared to be rich in lipid-like compounds that could be of plant and microbial origin, even though a higher abundance of these compounds was observed in RS samples.

TEM
The high resolution TEM imaging revealed a clear difference in appearance in the water-and methanol-extracted material (Figure 7), consistent with FTICR-MS data.The material dissolved in

TEM
The high resolution TEM imaging revealed a clear difference in appearance in the water-and methanol-extracted material (Figure 7), consistent with FTICR-MS data.The material dissolved in water showed ultrastructure indicative of morphologically indistinct precipitation of components dissolved in water upon air dying.In contrast, the SOM material extracted into the methanol phase exhibited circular features of size ranging from 50 nm to 1 µm suggestive of an emulsion of lipid-like compounds in methanol, consistent with FTICR-MS data (Van Krevelen diagrams in Figure 5).water showed ultrastructure indicative of morphologically indistinct precipitation of components dissolved in water upon air dying.In contrast, the SOM material extracted into the methanol phase exhibited circular features of size ranging from 50 nm to 1 µ m suggestive of an emulsion of lipid-like compounds in methanol, consistent with FTICR-MS data (Van Krevelen diagrams in Figure 5).

HeIM
Minerals from the rhizosphere showed adhesion of microbial EPS onto surfaces, as well as intercalation and deposition of OM into the interlayer space of biotite (Figure 8).
After extractions with H2O and methanol, mineral surfaces showed a residual OM signature when compared with the initial clean minerals (Figure 9).In addition, we also observed preferential microbial coverage on biotite edges and in mineral cracks and crevices.

HeIM
Minerals from the rhizosphere showed adhesion of microbial EPS onto surfaces, as well as intercalation and deposition of OM into the interlayer space of biotite (Figure 8).
water showed ultrastructure indicative of morphologically indistinct precipitation of components dissolved in water upon air dying.In contrast, the SOM material extracted into the methanol phase exhibited circular features of size ranging from 50 nm to 1 µ m suggestive of an emulsion of lipid-like compounds in methanol, consistent with FTICR-MS data (Van Krevelen diagrams in Figure 5).

HeIM
Minerals from the rhizosphere showed adhesion of microbial EPS onto surfaces, as well as intercalation and deposition of OM into the interlayer space of biotite (Figure 8).
After extractions with H2O and methanol, mineral surfaces showed a residual OM signature when compared with the initial clean minerals (Figure 9).In addition, we also observed preferential microbial coverage on biotite edges and in mineral cracks and crevices.After extractions with H 2 O and methanol, mineral surfaces showed a residual OM signature when compared with the initial clean minerals (Figure 9).In addition, we also observed preferential microbial coverage on biotite edges and in mineral cracks and crevices.

XRD and Micro-XRD
After 9 months, minerals did not show any significant sign of crystalline change, nor any secondary mineral formation (Figures 10 and 11).However, minerals increased the relative amorphous fraction.This phenomenon was particularly evident in biotite, with pronounced intercalation of the phyllosilicate layers as also revealed by HeIM (Figure 9).The level of amorphous material detected by micro-XRD decreased in the material after the water and methanol extraction.

XRD and Micro-XRD
After 9 months, minerals did not show any significant sign of crystalline change, nor any secondary mineral formation (Figures 10 and 11).However, minerals increased the relative amorphous fraction.This phenomenon was particularly evident in biotite, with pronounced intercalation of the phyllosilicate layers as also revealed by HeIM (Figure 9).The level of amorphous material detected by micro-XRD decreased in the material after the water and methanol extraction.

XRD and Micro-XRD
After 9 months, minerals did not show any significant sign of crystalline change, nor any secondary mineral formation (Figures 10 and 11).However, minerals increased the relative amorphous fraction.This phenomenon was particularly evident in biotite, with pronounced intercalation of the phyllosilicate layers as also revealed by HeIM (Figure 9).The level of amorphous material detected by micro-XRD decreased in the material after the water and methanol extraction.

Discussion
This study evaluated how SOM is formed, stored, and protected in a laboratory experiment using red pine seedlings grown in a designed mineral matrix with forest soil microbial inoculum.Our main findings were: (a) soil microorganisms in the initially OM-free column experiment were able to transform and stabilize red pine root exudates into a variety of C compounds with varying projected residence times, (b) a significant amount of plant-and microbially-derived biomass was contributed to the mineral matrix, and (c) residual microbial biomass stayed in the system even after extraction of soil minerals with different solvents.All of the main findings support our hypothesis that plant and microbial OM significantly contributes to the pool of stabilized C in soils.

Rhizospheric Microbiome
Despite the growing interest in characterizing the microbiome from the soil rhizosphere, there is less known about microbial communities adapted to survive solely on root exudates [46].In this nutrient-limited environment, the microbial community was preferentially located within the rhizosphere.The core community members in this ecological niche had to obtain energy from root exudates as their primary C source, and additionally, to access and mobilize cations (K, Ca) by mineral weathering.Our system contained a high percentage of Burkholderiales, Gram-negative proteobacteria, that are typically present in intimate association with fresh organic matter [47].Conversely, Actinomycetes, a group of Gram positive bacteria also present in our system, are adapted to nutrient poor environments [48].The greater abundance of Pseudomonadales and Acidobacter is attributed to these families being ubiquitous in soils.Also, nitrogen-fixing, plant root symbiotic Rhizobiales were present in our samples.Many genera within the relatively abundant Bradyrhizobiaceae are recognized for their nitrogen cycling capabilities, such as biological nitrogen fixation in the genus Bradyrhizobium or nitrogen mineralization from species in the genus Nitrobacter [49].Bradyrhizobium is also known to form symbiotic relationships with plants, mainly legumes, but have not been reported to form symbiotic relationships with pines.The majority of OTUs associated with Bradyrhizobiaceae observed in our samples had no assigned genus, and of the ones that did, Bradyrhizobium was the most abundant.Common carbon sources for microorganisms from the Bradyrhizobium genus include carbohydrates and organic acids, with nitrates being also a common nitrogen source [49].This suggests that in our samples Bradyrhizobium and other bacteria from the Rhizobiales were most likely in a free-living state and survived on organic acids from root exudates because there were very little carbohydrate-like compounds identified via FTICR in the rhizospheric soil.
The second most common family-level OTUs observed from Rhizobiales was Hyphomicrobiaceae, mainly from the genera Devosia, Hyphomicrobium, and Rhodoplanes.This diverse family is known for

Discussion
This study evaluated how SOM is formed, stored, and protected in a laboratory experiment using red pine seedlings grown in a designed mineral matrix with forest soil microbial inoculum.Our main findings were: (a) soil microorganisms in the initially OM-free column experiment were able to transform and stabilize red pine root exudates into a variety of C compounds with varying projected residence times, (b) a significant amount of plant-and microbially-derived biomass was contributed to the mineral matrix, and (c) residual microbial biomass stayed in the system even after extraction of soil minerals with different solvents.All of the main findings support our hypothesis that plant and microbial OM significantly contributes to the pool of stabilized C in soils.

Rhizospheric Microbiome
Despite the growing interest in characterizing the microbiome from the soil rhizosphere, there is less known about microbial communities adapted to survive solely on root exudates [46].In this nutrient-limited environment, the microbial community was preferentially located within the rhizosphere.The core community members in this ecological niche had to obtain energy from root exudates as their primary C source, and additionally, to access and mobilize cations (K, Ca) by mineral weathering.Our system contained a high percentage of Burkholderiales, Gram-negative proteobacteria, that are typically present in intimate association with fresh organic matter [47].Conversely, Actinomycetes, a group of Gram positive bacteria also present in our system, are adapted to nutrient poor environments [48].The greater abundance of Pseudomonadales and Acidobacter is attributed to these families being ubiquitous in soils.Also, nitrogen-fixing, plant root symbiotic Rhizobiales were present in our samples.Many genera within the relatively abundant Bradyrhizobiaceae are recognized for their nitrogen cycling capabilities, such as biological nitrogen fixation in the genus Bradyrhizobium or nitrogen mineralization from species in the genus Nitrobacter [49].Bradyrhizobium is also known to form symbiotic relationships with plants, mainly legumes, but have not been reported to form symbiotic relationships with pines.The majority of OTUs associated with Bradyrhizobiaceae observed in our samples had no assigned genus, and of the ones that did, Bradyrhizobium was the most abundant.Common carbon sources for microorganisms from the Bradyrhizobium genus include carbohydrates and organic acids, with nitrates being also a common nitrogen source [49].This suggests that in our samples Bradyrhizobium and other bacteria from the Rhizobiales were most likely in a free-living state and survived on organic acids from root exudates because there were very little carbohydrate-like compounds identified via FTICR in the rhizospheric soil.
The second most common family-level OTUs observed from Rhizobiales was Hyphomicrobiaceae, mainly from the genera Devosia, Hyphomicrobium, and Rhodoplanes.This diverse family is known for having oligocarbophilic organisms, i.e., organisms adapted to low carbon environments [50,51], which is consistent with the soil environment provided in this study.While many of Hyphomicrobiaceae genera are aerobic chemoheterotrophs (e.g., Devosia spp.) or chemoorganotrophs (e.g., Hyphomicrobium spp.); there are some that are facultative anaerobes that respire nitrate (e.g., Rhodoplanes spp.) [51,52].Hyphomicrobium spp.are widely distributed in soils, produce hyphae, and are especially adapted to growing in low resource conditions [51].Hyphomicrobium spp.can even thrive on methanol [51].Overall, Hyphomicrobium spp.has been observed growing autotrophically using CO 2 , H 2 , and O 2 [53].As such, it is not surprising that we observed microorganisms that thrive on low C sources and in both anaerobic and aerobic environments in soils that received no exogenous C substrates or other resources (with the exception of the root exudates) throughout the duration of our experiment.
We also observed a small abundance of Methylobacterium spp., a group of facultative methylotrophic bacteria, suggesting that methanol was present in the rhizosphere, possibly from the P. resinosa roots themselves [54].Methylobacterium have also been shown to form symbiotic relationships with plants and to promote plant growth via urease production or enzymes that promote disease resistance and have been used with other plant growth-promoting bacteria as inoculants in agricultural cropping systems [55].Thus, it appears that the microbial community located in the rhizosphere used various forms of C and a variety of electron acceptors.These findings were consistent with our FTIC-MS data that revealed a wide diversity of carbon compounds in the rhizosphere samples compared to the bulk soil.
Although pine naturally exist in a relationship with ectomycorrhizal and saprotrophic fungi which significantly contribute to C cycling by SOM decomposition [56,57], here we focused only on bacteria from the perspective of microbial ecology, to obtain information on the structure and function of the microbial community members that were able to sustain in these humics-free basic conditions.

Molecular Characterization of Rhizospheric Biomass
Molecular characterization of C compounds has become a necessity in studies focused on SOM stabilization mechanisms.Conventional thermal destructive procedures [58], and wet oxidative treatment methods for isolation of SOM in soils, including hydrogen peroxide (H 2 O 2 ) [59], potassium permanganate (KMnO 4 ) [60], disodium peroxodisulphate (Na 2 S 2 O 8 ) [61], or sodium hypochlorite NaOCl [19,62], all result in some alteration of mineral soil constituents.Thus, we used a SOM extraction method employing solvents with distinct polarities.Although the result was not a complete extraction of the SOM material adhered to the mineral surfaces, the portion of the extracted material was representative of the overall labile SOM [34].By this method, we were able to identify a wide variety of bio(macro)molecular signatures, ranging from low molecular weight root exudates (such as small protein-fragments and amino sugars) to high molecular weight substances (such as lignin-like and aromatic compounds).These bio(macro)molecular compounds were not present at the start of the experiment and were only produced in situ.These compounds have varying turnover times, ranging from labile to intermediate and refractory SOM pool components (Figure 1).These results were consistent with our microbial data that revealed greater abundance of microbial communities in the rhizosphere samples compared to the bulk material, where it was difficult to find and amplify a microbial DNA signal.Given that our sterilized soil initially did not contain any SOM material, we believe the majority of these compounds were produced by microbial activity.Additional evidence of this is the presence of a pool of refractory OM that was observed in the rhizospheric sample and was thus assumed to have originated from the microbial communities (i.e., lipid-like and aromatic compounds) (Figure 1).However, in the overall C budget, we must also consider another C source.Plant material produced by rhizodeposition, such as root mucilage, sloughed-off epidermal roots cells, and root hair and cortical dead root cells [63] also provide a variety of low molecular weight solutes [64] readily processed and respired by microorganisms as CO 2 .
The small amount of recalcitrant lipid-like and aromatic compounds in the bulk soil indicated a small signature of necromass that was transported away from the root by bulk flow into the surrounding area.Additionally, microorganisms in this nutrient-limited system needed to be located within the rhizosphere for survival, which explains the lack of microbial signal and the lack of SOM diversity in bulk soil samples.
The FTICR results indicated that a large variety of newly-formed organic C compounds were formed in this system, with their chemical attributes giving a perspective to the projected turnover time of the newly formed SOM (Figure 1).These data suggest that rhizosphere supported a relatively active and diverse ecosystem enriched with a wide range of organic labile and recalcitrant molecules compared to the bulk substrate samples that were mainly enriched with aromatic and lipid-like compounds, a part of the refractory SOM pool.However, given that SOM stabilization and decomposition is a multifaceted dynamic process, this can only provide fundamental information on SOM turnover modeling to help understanding the transport and fate of the SOM.

Correlated Imaging of Extracted Material
As previously mentioned, two solvents were used to extract a wide range of organic molecules from the soil because no single solvent would be able to extract SOM compounds that are known to have different polarities and different chemistries.The process of FTICR analyses using differences in SOM composition between water-and methanol-extracted organic compounds initiated the idea of their correlated visualization by electron microscopy, since to our knowledge, SOM has not been imaged by high resolution electron microscopy before.Comparative TEM imaging of soil organic matter extracted with water and methanol validated our perception of hydrophilic and hydrophobic SOM components, by showing their clear morphological differences (Figure 8).The water-extracted material was indicative of water-soluble low molecular weight organic substances, likely forming colloids suspended in water, consistent with the operational definition of dissolved organic matter (DOM) adopted: "the organic carbon smaller than 0.45 µm in diameter" [65].As revealed by FTICR, the major class of water-soluble dissolved organic carbon compounds contained carboxylic, hydroxyl, and phenolic functional groups.In contrast, the SOM material extracted with methanol contained micro-and nano-spherical structures resembling a lipid emulsion that corresponds with the concept of colloidal OM that contained components larger in molecular weight, and contain fewer carboxylic and hydroxyl functional groups than the water-soluble phase.The functional hydrophilic (C=O) and hydrophobic (C-H) groups of the SOM formed in the rhizosphere determine the wettable (hydrophilic) and water-repellent (hydrophobic) soil features, as they directly influence properties such as soil wettability, water repellency, and microaggregation [66].This in turn influences soil moisture retention, soil fertility, and water movement in soils.
We presume the presence of stabilized yet mobile particulate organic C compounds released in the aqueous phase had originated from microbial decomposition.This process is carried out by the oxidative depolymerization of organic C molecules, and it occurs due to the incorporation of oxygen-containing functional groups to SOM [67].This reduces the C compounds molecular size and ultimately improves water solubility.
Overall, the immense importance of carboxyl groups for both the aqueous solubility and the chemical reactivity of OM in soils has been well recognized.These chemical and physical properties have implications on chemical reactivity, stability, interaction with mineral surfaces, and OM mobility in soils, all of which ultimately define how persistent SOM is in soils [8,17].

SOM Association with Mineral Surfaces
Another microscopy imaging technique, Helium ion microscopy, revealed possible mechanisms of SOM stabilization in this system.The majority of organic compounds were intimately associated with mineral phases, as anticipated based on standard texts [68,69], and our previous observation of a similar system [32].For instance, frequently seen occlusion and intercalation may cause spatial inaccessibility against abiotic and microbial decomposition (Figure 8).Krull at al. 2003 [70] reviewed mechanisms and processes of SOM stabilization, and concluded that adsorption on mineral surfaces and particles aggregation can impede decomposition processes but 'molecular recalcitrance' appears to be the main mechanism by which soil OM can be stabilized for longer periods of time.Subsequently, Von Lützov et al. in their review [29] established that "The chemical structure of organic molecules by itself is not sufficient to account for the extreme variation of soil OM in terms of age and turnover times".Recently, Lehmann and Kleber proposed the soil continuum model that reclassifies the concept of soil mineral protection of organic matter from microbial decomposition [17].They recognized the SOM spatial architecture within the minerals matrix as a key factor how the SOM is protected from microbial access, rather than merely for its 'stable' chemical nature.Thus, the role of spatial inaccessibility is yet to be further discussed.
Given the absence of the initial OM substances in the designed mineral matrix, we attribute all observed and analyzed C compounds to de novo synthesis by rhizodeposition and microbial activity.As such, the robust presence of organic material associated with the mineral substrates even after the solvent treatments (Figure 9) suggests the stabilization of the organic matter adhered on the mineral edges and surfaces and intercalated into the phyllosilicate layers.We observed that mineral surfaces were used as templates for microbial attachment with EPS which suggests that the rhizospheric microorganisms acquired metal cations necessary for metabolic processes via microbial specific adhesion mechanisms [67].The subsequent biofilm formation supports creating the micro-habitat to access mineral-associated organic substrates, for nutrient uptake and mobilizing metal cations from minerals by microbial electron transfer [71,72].
Overall, our results confirmed the incorporation of the organic material into mineral matrix that resulted in formation of stabilized amorphous phase.Based on our multimodal analytical approach, this supports our projection that microbial biomass will be significantly contributing to the refractory C pool in soils.Further work needs be done to quantify the production and incorporation of organic material into the mineral matrix into soil organic C microbial models [73][74][75].

XRD of Minerals with Amorphous SOM
Different soil minerals possess different capacities to adsorb and protect organic matter as a result of mineral-specific properties such as surface area, charge characteristics, and porosity [70].Soil mineral matrices such as clays and poorly crystalline minerals [76] can influence C retention, stability and turnover due to their ability to physically protect organic C against decomposition.Our results confirmed this, as we were able to correlate the SEM images of minerals, especially biotite, with substantial organic C deposits, and with X-ray diffraction patterns characteristic for inclusion of amorphous matter.Although the X-ray diffraction and micro-XRD did not provide clear evidence of the initial minerals' crystallinity alteration or secondary mineral formation due to the mineral-organic associations (MOAs), this could be because of the qualitative and quantitative limitations of these methods.Thus, we do not discount the possibility of the initial stages of mineral weathering that were just not detectable at this time point.However, we noticed higher intensity narrower peaks in bulk soil, as compared to reduced intensity and a higher background "shoulder" in rhizosphere (Figure 10, arrows).These low-intensity, broad peaks are indicative of X-rays scattering distributed in a wide range (2 Theta), characteristic of an amorphous impurity, attributed to an increased amorphous content in the RS sample.This partial amorphization was caused by biogenic organic compounds associated with the mineral surfaces and incorporated within the mineral edges.The increased percentage of the amorphous material indicating the incorporation of the material produced by rhizospheric activity on and within the minerals (the above-mentioned intercalation) remained even in the material that was washed in water and methanol, indicating refractory signatures of newly formed biogenic material firmly attached to the mineral substrates (Figure 11).
Additionally, after the treatment with solvents that removed different amounts of organic matter, as evidenced by SEM imaging, the XRD patterns still showed pronounced signs of incorporation of amorphous material into their crystalline mineral content, indicating a certain degree of protection exerted by the mineral matrix.We assume the solvents removed the labile organic C from the samples, with the majority of intermediate and refractory OM remaining attached to the mineral surfaces.With water as the most common solvent present in natural systems, we assume the concentration of hydrophobic compounds on mineral surface would increase in time.

Conclusions
Our study showed that soil microorganisms in the initially OM-free column experiment were able to transform red pine root exudates into a variety of C compounds with various degrees of inherent recalcitrance, and a substantial amount of biomass became incorporated in the system.We also saw evidence that residual microbial biomass persisted in the system even after extractions with different solvents, which suggested that microbial necromass may contribute to the intermediate and refractory SOM pool in soils.Our short-term results support the idea that SOM binding to soil minerals [77] is a major pathway for organic carbon sequestration in soil.This study contributed to advancing our understanding of the formation and fate of root carbon, and the persistence of SOM requires disentangling the complex interactions between root exudates and debris, microbial communities, and soil minerals.

Figure 2 .
Figure 2. A schematic illustration of the experimental design.

Figure 2 .
Figure 2. A schematic illustration of the experimental design.

Figure 3 .
Figure 3. Rhizosphere soil microbial community composition identified by DNA sequencing.(A) A bar plot of the most abundant microbial families grouped by order.(B) A phylogenetic tree shows the relatedness of these most common microbial families.

Figure 4 .
Figure 4. Relative abundance of operational taxonomic units from the order Rhizobiales separated by genus (x-axis) and family (colored legend).

Figure 3 .
Figure 3. Rhizosphere soil microbial community composition identified by DNA sequencing.(A) A bar plot of the most abundant microbial families grouped by order.(B) A phylogenetic tree shows the relatedness of these most common microbial families.

Figure 3 .
Figure 3. Rhizosphere soil microbial community composition identified by DNA sequencing.(A) A bar plot of the most abundant microbial families grouped by order.(B) A phylogenetic tree shows the relatedness of these most common microbial families.

Figure 4 .
Figure 4. Relative abundance of operational taxonomic units from the order Rhizobiales separated by genus (x-axis) and family (colored legend).

Figure 4 .
Figure 4. Relative abundance of operational taxonomic units from the order Rhizobiales separated by genus (x-axis) and family (colored legend).

Figure 5 .
Figure 5. Van Krevelen diagrams of SOM compounds extracted by water and methanol, from rhizosphere (RS) and bulk substrate (BS).

Figure 6 .
Figure 6.Relative abundance of different compound classes uniquely extracted by water and methanol, from rhizosphere (RS) and bulk substrate (BS).

Figure 5 .
Figure 5. Van Krevelen diagrams of SOM compounds extracted by water and methanol, from rhizosphere (RS) and bulk substrate (BS).

Figure 5 .
Figure 5. Van Krevelen diagrams of SOM compounds extracted by water and methanol, from rhizosphere (RS) and bulk substrate (BS).

Figure 6 .
Figure 6.Relative abundance of different compound classes uniquely extracted by water and methanol, from rhizosphere (RS) and bulk substrate (BS).

Figure 6 .
Figure 6.Relative abundance of different compound classes uniquely extracted by water and methanol, from rhizosphere (RS) and bulk substrate (BS).

Figure 9 .
Figure 9. HeIM images of mineral surfaces after extraction with water and methanol showing a significant quantity of organic matter residue, as compared to the initial clean minerals.Original clean anorthite (A); anorthite after methanol sonicated washes (B,C).Original biotite material (D); biotite after water and methanol extraction, resp.(E,F).

Figure 10 .
Figure 10.Micro-XRD patterns of the bulk substrate (BS) and rhizosphere soil (RS).The broadened "shoulder" peak (arrow), as well as lower intensity of signature peaks of biotite (B), quartz (Q), anorthite (A), and corundum (C) is characteristic for increased amorphous content.

Figure 9 .
Figure 9. HeIM images of mineral surfaces after extraction with water and methanol showing a significant quantity of organic matter residue, as compared to the initial clean minerals.Original clean anorthite (A); anorthite after methanol sonicated washes (B,C).Original biotite material (D); biotite after water and methanol extraction, resp.(E,F).

Figure 9 .
Figure 9. HeIM images of mineral surfaces after extraction with water and methanol showing a significant quantity of organic matter residue, as compared to the initial clean minerals.Original clean anorthite (A); anorthite after methanol sonicated washes (B,C).Original biotite material (D); biotite after water and methanol extraction, resp.(E,F).

Figure 10 .
Figure 10.Micro-XRD patterns of the bulk substrate (BS) and rhizosphere soil (RS).The broadened "shoulder" peak (arrow), as well as lower intensity of signature peaks of biotite (B), quartz (Q), anorthite (A), and corundum (C) is characteristic for increased amorphous content.

Figure 10 .
Figure 10.Micro-XRD patterns of the bulk substrate (BS) and rhizosphere soil (RS).The broadened "shoulder" peak (arrow), as well as lower intensity of signature peaks of biotite (B), quartz (Q), anorthite (A), and corundum (C) is characteristic for increased amorphous content.

Figure 11 .
Figure 11.Results of micro-XRD quantitative analysis of the amorphous material in the rhizosphere and bulk substrate from the columns.Material was extracted with water (H2O), methanol (MeOH), and non-extracted (NE).The initial starting sonicated substrate was used as a control.Even after the extractions, amorphous material remained associated with the minerals.The error bars are within 1.0%.

Figure 11 .
Figure 11.Results of micro-XRD quantitative analysis of the amorphous material in the rhizosphere and bulk substrate from the columns.Material was extracted with water (H 2 O), methanol (MeOH), and non-extracted (NE).The starting sonicated substrate was used as a control.Even after the extractions, amorphous material remained associated with the minerals.The error bars are within 1.0%.