Root-Derived Short-Chain Suberin Diacids from Rice and Rape Seed in a Paddy Soil under Rice Cultivar Treatments

Suberin-derived substituted fatty acids have been shown to be potential biomarkers for plant-derived carbon (C) in soils across ecosystems. Analyzing root derived suberin compounds bound in soil could help to understand the root input into a soil organic carbon pool. In this study, bound lipids were extracted and identified in root and topsoil samples. Short-chain suberin diacids were quantified under rice (Oryza sativa L.) and rape (Brassica campestris) rotations with different cultivar combinations in a Chinese rice paddy. After removal of free lipids with sequential extraction, the residual bound lipids were obtained with saponification and derivatization before analysis using gas chromatography–mass spectrometry (GC-MS). Diacids C16 and C18 in bound lipids were detected both in rice and rape root samples, while diacids C20 and C22 were detected only in rape root samples. Accordingly, diacids were quantified in both rhizosphere and bulk soil (0–15 cm). The amount of total root-derived diacids in bulk soil varied in a range of 5.6–9.6 mg/kg across growth stages and crop seasons. After one year-round rice-rape rotation, root-derived suberin diacids were maintained at a level of 7–9 mg/kg in bulk soil; this was higher under a super rice cultivar LY than under a hybrid cultivar IIY. While concentrations of the analyzed diacids were generally higher in rhizosphere than in bulk soil, the total diacid (DA) concentration was higher at the time of rape harvest than at rice harvest, suggesting that rape roots made a major contribution to the preservation of diacids in the paddy. Moreover, the net change in the concentration and the ratios of C16:0 DA to C18:1 DA, and of C16:0 DA to C18:0 DA, over a whole growing season, were greater under LY than under IIY, though there was no difference between cultivars within a single growth stage. Overall, total concentration of root-derived suberin diacids was found to be positively correlated to soil organic carbon concentration both for bulk soil and rhizosphere. However, the turnover and preservation of the root suberin biomolecules with soil property and field conditions deserve further field studies.


Introduction
Lipids from plant sources are a significant source of carbon (C) input to soil [1,2]. Characterizing free or bound fractions of lipids can help elucidate soil organic matter (SOM) dynamics with land use and management changes [3][4][5]. Lipid compounds are a group of apoplastic biopolymers with tissue-specific composition in the cell walls of endodermis and exodermis and of periderms such as wound periderm [6]. However, suberin is an analogous bio-polyester, containing polyphenolic and polyaliphatic compounds, that is found mostly in root tissues [7]. Furthermore, suberin-derived substituted fatty acids (SFA) could be preserved after incorporation into soil organic matter in soil aggregates surrounding roots [8], owing to their relatively long carbon chains and molecular hydrophobicity. Thus, suberin compounds have been considered a potential tracer to track C input from plant roots and root-C turnover in soil across various ecosystems [8][9][10].
The aliphatic domain of suberin has been shown to be composed mostly of carboxylic acids, ω-hydroxy carboxylic acids, α, ω-dicarboxylic acids (diacids, DAs) and homologous midchain di-hydroxy or epoxy derivatives [11]. Among these monomers, ω-hydroxy fatty acids are not used as root biomarkers, though they are generally considered to be the most prominent substance class [12], for example, for wheat and maize roots [13]. With the development of extraction and quantification procedures [14], diacids of suberin fatty acids, as root derived biomolecules in bound lipids, have been characterized for exploring SOM dynamics across a wide range of soils [10,13]. Above-ground tissues of grain crops, often removed after harvest, are less suberized than that of woody plants [13] and thus contribute less to SOM in agricultural soils. Instead, diacids extracted from soil may better represent crop root-derived biomolecules for assessing root input to SOM. Furthermore, the aliphatic molecules of suberin diacids could exert higher stability in soil than lignin compounds [8,15]. Therefore, analysis of these diacids in soil could provide important information about crop root-derived C to the stable pool of SOM in agricultural soils.
Rice (Oryza sativa L.) is cultivated on 155 million hectares of world croplands, occupying approximately 14% of all global arable land [16]. Although rice production on paddy soils has been considered an important source of biogenic greenhouse gases [17], paddies, with improved management [18,19], could have higher SOM storage and greater C sequestration potential than they currently have. As one of the major cultivated soil types in China [20], rice paddy soils have high soil organic carbon storage via carbon sequestration from root input carbon [21] through enhanced physical protection [22] and chemical binding with oxyhydrates of Fe/Al [23]. Recently, molecular recalcitrance [24] and biological stabilization [25] have also been suggested for SOM accumulation and stabilization in rice paddies. Over the last few decades, many studies have shown higher increases in SOM in paddy soils than in dry croplands of China [19,26] and fast accumulation of SOM in soils recently converted to rice cultivation [27]. The newly accumulated SOM observed in bulk soil has often been credited for an increase in physically protected SOM in micro-aggregates [27]. However, there have been few studies reporting the abundance of root-derived biomolecules as related to root-C input from rice and rotated crops, and their relative contributions to SOM in rice paddies [28].
Super hybrid rice (super rice in short) is a varietal type of rice that combines different plant varieties with heterosis and achieves super high yields through hybridization between indica and japonica [29,30]. Super rice cultivars are characterized by higher root biomass and higher leaf-area index [31] with enhanced photosynthesis and biomass production. This could, in turn, potentially increase root-C input via rhizodeposition of photosynthesized carbon hydrates. While changes in SOM after continuous cultivation of super rice have been discussed [32], there is no field study addressing root-C contribution to rice soil via root activity in field conditions. This study's objective is to demonstrate the abundance and composition of suberin diacids as root-derived biomolecules of rice and the rape crops, and to characterize the changes in crop-growth stages in a rice paddy under continuous super rice breeding rotated with rape. The root suberin diacids were analyzed using a modified extraction and GC-MS detection assay. Concentrations of the suberin diacids were compared both between bulk and rhizosphere soil samples and between rice cultivars in order to infer potential difference in root-derived C input into paddy soil under rice-rape rotation.

Ethics Statement
In this study, no specific permission was required either for the soil/plant sampling at the site or for the biomolecular analysis at the laboratory. No endangered or protected species were involved in the field work.

Experimental site and crop cultivation
Field work of soil and plant sampling was conducted in a demonstration farm for high yield rice breeding, located in Jianhua village (117°47'E, 31°39'N), Zhonghan Town, Anhui Province, China. The farm was under management by the Rice Research Institute, Anhui Academy of Agricultural Sciences, China. The paddy soil was a Hydroagric Stagnic Anthrosol [33] derived from lacustrine deposits. The experiment with super rice breeding was established in 2005 with continuous cultivation of super rice cultivars in comparison with hybrid rice. The rice paddy had been cultivated under summer rice (Oryza sativa L.) and winter rape seed (Brassica campestris) rotation in recent decades. For summer rice, a super rice cultivar belonging to the LY series and a conventional hybrid rice belonging to the ⅡY series had been cultivated, while for winter rape, two parallel cultivars, both belonging to the QY series, were cultivated in a year round rotation. The experimental design of this study is shown in Table 1. Each cultivar treatment was conducted in triplicate and the plots were arranged in a randomized complete block design.
For the samples used in this study, rice was sown on May 10, 2011 and transplanted with a spacing of 30 cm between rows and 13.3 cm in rows on June 11, 2011. Hybrid and super rice cultivars were fertilized with 180 and 210 kg N/ha of urea, 33 and 39 kg P/ha of superphosphate, and 62 and 124 kg K/ha of potassium chloride, respectively. While phosphorus fertilizer of superphosphate was applied as base fertilizer only, total N was split by 50%, 20% and 30%, respectively, before sowing, tillering and spiking, for both cultivars. Meanwhile, 60% of potassium chloride was applied as basal fertilizer and the other 40% was applied at the spiking stage. The soil moisture regime was managed with a scheme of flooding before spiking and waterlogged intermittently thereafter until harvest. Plant protection was consistent across the cultivars. Rape seed was transplanted after rice harvest with a spacing of 33.3 cm between rows and 25.0 cm in rows, totaling 120,000 plants/ha. Rape seed cultivars were fertilized with the same amount of 210 kg N/ha, with a ratio of 1:0.43:0.71 (N:P:K) across the two cultivars.

Collection and treatment of soil and root samples
In this study, soil samples were collected before, during and after crop growth for analysis of root suberin diacids as well as for soil physico-chemical analysis. Specifically, a sampling of initial bulk topsoil (0-15cm) was done after rape harvest, on the 25 th of May 2011, before a new rotation year of crop production in 2011-2012. Over the course of the rice growing season, both bulk soil and rhizosphere samples were collected in each cultivar plot at the end of tillering (8 th July of 2011), at grain filling stage (5 th September of 2011), and at harvest (13 th October of 2011). Meanwhile, both bulk soil and rhizosphere samples under rape cropping were also collected at harvest (26 th May of 2012). Rhizosphere samples were obtained from the crop root zone by shaking and washing the soil attached to crop roots, following a procedure described in detail by Butler et al. [34] and by Wiesenberg et al. [35]. The remaining root was preserved as root samples. Bulk samples of topsoil, however, were collected randomly in a spacing area between the plants. For all these samples, sampling was done in random triplicate in order to form a composite sample for a single treatment plot. After being shipped to the laboratory, any plant detritus and gravel was removed from the collected samples. After air-drying in room temperature, the soil samples were ground and passed through a 2 mm sieve. The root samples were washed and then oven-dried at 60°C. A portion of soil and root samples was further ground with a stainless steel grinder to pass through a 0.15 mm sieve for root suberin diacids assay and SOM-C analysis.

Analysis of soil properties
Soil basic properties were measured following the standard protocols described by Lu [36]. To summarize, soil pH was measured with a glass electrode using a 1:2.5 soil/water ratio with a precision pH meter (Mettler Toledo Seveneasy, Switzerland). SOM concentration was measured using wet digestion and oxidation with potassium dichromate. Total nitrogen was analyzed using the micro Kjeldahl method. Total phosphorus was digested with perchloric acid and sulfuric acid, and analyzed with colorimetry (TU-1810) at 700 nm. Particulate organic matter (POM) was measured according to the procedure described by Cambardella and Christensen [37] as follows: Dispersing 20 g dried soil (<2 mm) in a solution of 100 mL 0.5 mol/L NaOH, and filtering the suspension through a 53 μm sieve. The material remaining on the sieve was analyzed for carbon with a CNS analyzer (Variomax CNS Analyzer, German Elementar Company, 2003). Dissolved organic matter (DOM) was extracted with CO 2-free distilled water at 25°C, and then measured using a TOC analyzer (Jena Multi N/C 2100) [32]. Microbial biomass carbon (MBC) was determined using the fumigation-extraction method [38]. The cation exchange capacity (CEC) was measured with the ammonium acetate (1 mol/L, pH 7.0) leaching method. Concentration of soil particles was determined using a hydrometer after dispersion with 0.5 mol/L NaOH as described by Zhou et al. [39]. Basic properties and organic matter pools of topsoil are presented in Tables 2 and 3, respectively.
Root suberin diacids analysis 1. Processing of soil samples. Soil samples were prepared for root suberin diacids analysis as follows: Step 1, Delipidation: Before depolymerisation of suberin, a portion of a soil sample (<0.149 mm) was extracted three times using a 1:10 w/v dichloromethane-methanol (2:1, v/v) solution in room temperature for 10 hours in order to remove free lipids [40,41]. The extract was separated from soil particles by centrifugation (3500 rpm, 15 min), and the supernatant was decanted.
Step 2, Saponification: Soil residues obtained in Step 1 were depolymerized to release suberin compounds. The residues obtained were shaken for 22 hours in an aqueous solution of potassium hydroxide (1.0 M KOH) in methanol under a N 2 atmosphere [13,42].
Step 3, Derivatization: The product of saponification was trans-esterified by a 12% BF 3 -MeOH solution under a N 2 atmosphere at 75°C for 0.5 h, and the residual BF 3 solution was quenched with water [43]. The solution was filtered before the residue was washed using a methanol solution (methanol/water 9:1, v/v) [13].
Step 4, Extraction: The filtrate from Step 3 was combined with 4 mL saturated NaCl solution, and subsequently extracted with dichloromethane three times to obtain hydrophobic monomers [40,44]. The extraction was then washed with 0.5 mol/L NaCl solution three times, and the pH was adjusted to 7.0 with KOH solution. After the supernatant was decanted, the residue was combined with anhydrous Na 2 SO 4 and dried in room temperature [45] prior to GC analysis.
2. Processing of root samples. The ground root samples were soaked for 10 hours and then washed thoroughly with a mixture of chloroform and methanol (2:1, v/v). The residues were extracted with chloroform in a Soxhlet extractor for 48 hours. After air-drying in room temperature, they were treated for 16 hours with a mixture of cellulase (5 g/L) and pectinase (1 g/L) in a 0.05 M acetate buffer (pH 4) at 30°C with rotatory shaking at 220 rpm in a Thermostatic Oscillator (HZ-9211KB, Shanghai, China). Subsequently, they were washed thoroughly with distilled water and then dried in an oven at 60°C. The above-mentioned procedure was repeated twice. The final residues were washed with a mixture of chloroform and methanol (2:1, v/v) and then dried in room temperature [46]. Finally, the saponification, derivatization, extraction and drying of the root samples for suberin diacids analysis were conducted using the same procedure as for the soil samples.

GC-MS and GC analysis
1. Qualitative analysis. The dried extracts were dissolved in a solution of pyridine or nheptane and filtered through a 0.45 μm filtration membrane, before being analyzed by gas chromatography-mass spectrometry (GC-MS).
The GC-MS analysis involved an HP-5 capillary column (HP-5 MS 5% Phenyl Methyl Silox: 939.53782, 325°C: 30 m × 250 μm × 0.25 μm film thickness) with Helium carrier gas at 34 mL/min constant flow and in an oven that was kept at 60°C for 2 min, then programmed to shift from 60°C to 310°C at 10°C/min, and then held for 2 minutes at 310°C. Samples were injected in split mode (30:1 ratio at 280°C injector temperature). Being more soluble than n-heptane in preparation experiments, pyridine was used as a supplementary solvent for extraction of the suberin fatty acids [45]. Finally, a mixture of pyridine and n-heptane was used as solvent for quantitative analysis with GC to improve the extraction in this study.
2. Quantitative analysis. Trans-esterified monomers dissolved in the solution of n-heptane mixed with pyridine were separated with an Agilent 7890 gas chromatograph (GC) which was equipped with an Agilent 19091J-413 HP-5 5% Phenyl Methyl Siloxan column (J&W Scientific, CA, USA; 30 m × 320 μm × 0.25 μm). The GC oven was kept at 60°C for 1 min, then programmed to increase from 60°C to 220°C at 40°C /min, subsequently from 220°C to 280°C at 3°C /min, and finally from 280°C to 325°C at 15°C/min. Trans-esterified monomers were quantified with a flame ionization detector (FID) coupled to the GC, using Nonadecanoic methyl ester (NU-CHEK PREP, INC) dissolved in pyridine as an internal standard [13,47].

Data processing and statistics
Data were expressed as means plus/minus standard deviation (mean ± SD). The difference between cultivars was examined using one-way analysis of variance (ANOVA) while the difference between rhizosphere and bulk soil samples was examined using a paired t-test. Bivariate correlation analysis was performed to find the relationship between diacids and SOM-C. A significant difference or correlation was defined at p<0.05. All statistical analyses were carried out using SPSS 16.0.

Basic soil properties and organic matter pools under cultivars
As shown in Table 2, there were similar physico-chemical properties of the paddy soil after continuous rice cultivation across the cultivar treatments. However, there were significantly higher concentrations of POM (by 15%), DOM (by 10%) and MBC (by 25%) under cultivars LY than under II Y. In contrast, no significant difference in total SOM was seen between cultivars LY and II Y ( Table 3).

Monomers of root suberin diacids detected in root and rhizosphere samples
Four major classes of organic compounds; alkanes, fatty acids, diacids and aromatic compounds, were identified in the bound lipids both of rhizosphere and root samples (Fig 1 and S1  Table). Here, the significant difference in aromatic compound abundance between rape and rice samples could be due to the difference in the solvent used. Dissolution was conducted for rape in n-heptane only and for rice root in the mixture of n-heptane and pyridine. Fatty acids were dominant in soil samples, alkanes dominated in rice roots, while fatty acids and diacids were predominant in rape roots. With this method, ω-hydroxy fatty acids were not able to be detected. Dissolution of root suberin compounds in the mixture allowed successful detection of the diacids for rice though dissolution in n-heptane, without pyridine, was also successful for rape roots.
Furthermore, the major root biomolecules under the cultivated rice were the short-chain diacids of C16 and C18, similar to those found in study [12]. Being indicative of rice root-input biomolecules, monomers of hexadecanedioic acid (C 16:0 DA), octadecenedioic acid (C 18:1 DA) and octadecanedioic acid (C 18:0 DA) were all detected in the rice roots. Besides these diacids, low levels of longer-chain monomers of eicosanedioic acid (C 20:0 DA) and docosanedioic acid (C 22:0 DA), typically from rape plants, were also detected in rape root and soil samples but not in rice roots. Thus, the modified procedure was able to detect the root-derived biomolecules both for rice and rape, though not the ω-hydroxy fatty acids monomers, in the rice paddy under rice and rape rotation (S1 and S2 Figs).

Root suberin diacids under rice cropping
Concentrations of root suberin diacid monomers (SC16-18), calculated by averaging the measurements across the growing stages, ranged from 1.0 mg/kg to 8.0 mg/kg of dry soil and were generally higher in rhizosphere than in bulk soil samples (Fig 2). Notably, concentrations of C 16:0 DA and total DAs were significantly higher in rhizosphere than in bulk soil under super rice (LY) while no significant difference occurred under hybrid cultivars (II Y). More root suberin diacids were preserved in rhizosphere samples under super rice LY than under the hybrid rice.
Concentrations of the diacids and SOM-C, across rice growing stages in bulk soil and rhizosphere samples, are presented respectively in Table 4 and Table 5. While no consistent changes in the concentrations of different monomers were observed across the growing stages, there was a slight decrease in the total concentration of diacids (SC16-C18) in the bulk soils, as the crops approached the ripening stage. Compared to concentrations before transplantation, the diacids monomers generally decreased at the filling stage and increased again at harvest, both for bulk soil and rhizosphere samples. There were no significant differences between rice cultivars in total DAs concentration and in different monomers both for bulk soil and rhizosphere samples except in the monomer of C 16:0 DA in bulk soil. Here, a net change over a whole ricegrowing season was calculated and assessed. This showed a net change in concentration of C 16

Root suberin diacids under rape cropping
Under rape cropping the soil contained a wide range of monomers of root suberin diacids viz. C 16:0 DA, C 18:1 DA, C 18:0 DA, C 20:0 DA and C 22:0 DA. Data of the concentrations of these monomers and SOM in bulk soil before rape transplantation (S3), and both in the rhizosphere and bulk soil at rape harvest, (R2), are illustrated in Fig 3. The concentrations of all monomers detected varied in a range of 1.2-5.7 mg/kg for rhizosphere and of 0.8-2.5 mg/kg for bulk soils. After rape cropping, the diacid concentrations generally increased greatly in the rhizosphere and moderately in bulk soils, compared to the bulk soil after rice cropping. Similar to the overall trend under rice (Fig 2), the root suberin diacid concentrations under rape were generally higher in the rhizosphere than in bulk soil. There were generally significant differences in the concentrations of most monomers detected between cultivar treatments. In particular, the concentration of C 22:0 DA, a potential rape root biomarker, was much greater (by up to 40%) both in bulk soil and rhizosphere, in LY plots than in IIY plots. Therefore, higher concentrations of root suberin diacids were preserved after rape cropping under super rice (LY) than under hybrid cultivar (IIY). Table 4. Concentration of SOM-C (g/kg) and root suberin diacids (mg/kg) and their ratios in bulk soil under rice across growing stages.

Growth stage
Cultivar SOM-C C 16

Discussion
Variation of root suberin diacids with rice growing stages and cropping seasons Using the modified protocol with GC-MS, the concentration of root suberin diacids could be successfully determined both in the bulk soil and in the rhizosphere. The data presented above show the differences in the concentration of different monomers, and total diacids between rhizosphere and bulk samples at the different rice growing stages as well as for rape cropping. The protocol used did not allow the measurement of ω-hydroxy acids, which were generally considered to be shoot-derived plant aliphatic biomolecules [13]. Here, short-chain monomers of C 16:0 DA and C 18:1 DA as the major diacids and of C 18:0 DA as the minor diacid were detected under rice cultivation, in addition to these, longer-chain monomers of C 20:0 DA and C 22:0 DA were detected in small quantities under rape cultivation. This is similar to the findings of the study [12] which addresses the different role of the monomers in the transport of water into the plants. The dominance of diacid monomers of C 16:0 DA and C 18:1 DA over C 18:1 DA was also reported for tree-root-derived biomarkers [9]. While there were slight inconsistencies between the cultivars, the diacid concentrations were generally higher in the rhizosphere than in bulk soil, especially after rape cultivation. As relatively fresh OM, root-derived diacids could decompose over time, more rapidly in aerobic conditions under rape cultivation than under intermittent waterlogging with rice cultivation [48]. This could be explained, in part, by a bigger difference in the diacid concentrations between rhizosphere with prominent root input, and bulk soil without prominent root input, under rape cropping. In addition, much higher diacid concentrations in the soil at the rape harvest compared to at the rice harvest, suggests a major contribution by rape roots to the soil storage of root-derived suberin diacids under ricerape rotation. Data in Table 5 exhibit hardly any difference in rhizosphere concentrations of the root suberin diacids across rice growing stages. Furthermore, no significant differences were observed in bulk soil samples except in the concentration of monomer C 16:0 DA across all stages. This lack of difference was in contrast to the findings by Mendez-Millan et al. [13] that root-derived suberin aliphatic compounds were different both in isotopic and monomers composition between C3 (wheat) and C4 (maize) plants. Wiesenberg et al. [5] found that there was a difference in the composition of carboxylic acids, especially the ratio of n-C24 versus n-C22+26 and their respective C-isotope values, for different types of plant roots. Similarly, when calculating the ratios of different monomers, we found differences in ratios of monomer C 16:0 DA to C 18:1 DA and of C 16:0 DA to C 18:0 DA between different stages, especially between tillering and harvest stages for cultivar LY both in the bulk and rhizosphere samples. This corresponded to the difference in concentration of monomer C 16:0 DA between the stages. This finding suggested that the composition pattern of the diacids could have changed with growing stages. There were certainly differences in field soil properties, nutrient status and other soil conditions between the rice growing stages with different rice production managements. In a study [35] using both C3 and C4 plants, the distribution pattern of plant lipid biomolecules was not affected by short term artificial pulse 14 CO 2 labelling but rather it was sensitive to environmental changes. Therefore, there could be considerable differences in the composition patterns with the rice, though differences between cultivars were insignificant. Root biomarker concentrations and soil organic matter measured in carbon (mean ± SE, n = 3) in the soil under rape seed. S3, before transplantation of rape seed (after harvest of the rice), R2, at harvest of rape seed; RH, the rhizosphere soil; BS, bulk soil; LY, super rice; ⅡY, hybrid rice. LY-R, super rice-rape rotation; ⅡY-R, hybrid rice-rape rotation. Different letters represent different statistical class among the cultivars at p<0.05. doi:10.1371/journal.pone.0127474.g003

Relationship between root suberin diacids and soil organic carbon
In this study, concentrations of root suberin diacids of C 16:0 DA, C 18:1 DA, C 18:0 DA, C 20:0 DA and C 22:0 DA, and SOM concentrations both of bulk soil and rhizosphere across the crop growing stages followed a pseudo-normal distribution as shown by 2-tailed single series Kolmogorov-Smirnov Test. There was a significant positive correlation between total DAs and SOM-C concentration both for bulk soil and rhizosphere throughout a whole rice-rape rotation year (Fig 4). This correlation could indicate SOM-C changes in the rice paddy over the cultivar treatment plots in relation to their root-derived organic carbon input. Organic matter accumulation in soils has often been shown to be associated with the preservation of hydrophobic organic matter, mostly in the form of lignin and lipids [49] including the suberin fatty acids analyzed here. Hydrophobic organic matter is known to indicate SOM sequestration from crop straw input in paddy soils rich in free oxyhydrates of Fe/Al [50]. Nierop et al. [51] argued that ester-linked macromolecules such as cutin and suberin could be used as stable indicators of plant-derived organic matter in soils. Moreover, a greater portion of bound lipids to total SOM shows SOM stability, as bound lipids have been shown to be more stable than bulk SOM [15].
The significant correlation of diacids with SOM-C observed here could suggest that root suberin biomolecules of diacids are an important component of plant-derived OM contributing to C storage in soil. Here, we compared the changes in the concentrations of the different monomers after rice cropping and found that the concentration of C 16:0 DA and the ratio of C 16:0 DA to C 18:0 DA decreased less, corresponding to a smaller decrease in SOM-C, under cultivar LY than under IIY (Tables 4 and 5). Meanwhile, a significant increase in POM-C of bulk soil was observed after rice harvest under LY compared to under IIY (Table 3). An increase in POM could represent an increased input of plant-derived organic matter under super rice cultivation, because POM has been generally considered to be relatively fresh organic matter [32]. Meng et al. [32] reported that continuous cultivation of super rice resulted in a reduction in indigenous SOM. However, in their study, a higher percentage of the coarse sand micro-aggregates (2-0.2 mm in size) were found under continuous super rice cultivation compared to under hybrid rice [32]. Presumably, more root-derived C in the soil under the super rice-rape rotation had been preserved compared with the hybrid rice-rape rotation. The increase in POM under super rice LY here coincided with these findings, indicating higher preservation of fresh organic input from rice roots under LY than under IIY.
Furthermore, DOM-C and MBC concentrations in bulk soil at rice harvest were also higher under super rice cultivar LY than under hybrid rice IIY, suggesting relatively high available C substrate under the super rice cultivar (Table 3). It has been well known that microbial growth could be greatly promoted with fresh organic matter input [52], which often induces increased soil respiration and C turnover in agricultural soils [53]. While total SOM-C was not observed as being significantly reduced under LY relative to IIY, the increase in fresh organic matter under super rice did not appear to compromise the potential for enhanced C mineralization in the short term (5 years of continuous super rice breeding in this site) [32]. It is still not clear if the change in SOM concentration with super rice cultivation could lag behind the change in the labile SOM pools such as POM, DOM as well as MBC. In this study, however, the concentrations of C 16:0 DA and C 22:0 DA diacids as root derived biomolecules were very significantly and positively correlated with SOM-C in bulk soil (p = 0.009, 0.022) and in the rhizosphere (p = 0.043, 0.013). Nevertheless, as there was no net increase in the total concentrations of the root suberin diacids over a whole rice growing season, the mechanism of the preservation of these root-derived suberin diacids in the rice soil and their long term fate needs further study.

Conclusions
Root suberin diacids were quantified both in bulk and rhizosphere samples in paddy soil under different cultivars across a whole rice and rape rotation year. While a difference in the monomer concentrations was hardly found between cultivar treatments in a single stage, the composition pattern of the monomers was altered with rice growing stages under cultivar LY. Higher concentrations of rice root suberin diacids were generally observed in rhizosphere than in bulk soil, and after rice cropping, there was a smaller decrease in them under a super rice cultivar LY than under a hybrid cultivar IIY. This change was in line with the changes in labile organic matter pools in the bulk soil. Significant correlations of these root suberin diacids with soil organic matter both in bulk soil and rhizosphere could suggest the potential for fresh root-derived organic matter accumulation in the rice paddy. However, the turnover and preservation of the root suberin biomolecules with soil property and field conditions are still not clear from short-term cultivation, and this deserves further field studies.
Supporting Information S1 Table. Compounds identified in base hydrolysis products of rice and rape root, and soil.