Subscriber access provided by University of South Dakota
Agricultural and Environmental Chemistry
Soil Perturbation in Mediterranean Ecosystems Reflected by Differences in Free Lipid Biomarker Assemblages Pilar Tinoco, Gonzalo Almendros, and Jesús Sanz J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.8b01483 • Publication Date (Web): 04 Sep 2018 Downloaded from http://pubs.acs.org on September 5, 2018
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 34
Journal of Agricultural and Food Chemistry
== AGRICULTURAL AND ENVIRONMENTAL CHEMISTRY ==
RESEARCH ARTICLE:
RUNNING TITLE: Soils lipid signature
Soil Perturbation in Mediterranean Ecosystems Reflected by Differences in Free Lipid Biomarker Assemblages
†
‡*
§
PILAR TINOCO , GONZALO ALMENDROS , JESÚS SANZ
†
Universidad Alfonso X el Sabio, Campus de Villanueva de la Cañada. Av. Universidad 1, E-28691 Madrid ‡
§
MNCN, CSIC, Serrano 115B, 28006-Madrid, Spain
Instituto de Química Orgánica General, CSIC, Juan de la Cierva 3, E-28006 Madrid, Spain
*
Corresponding author. E-mail:
[email protected] (G. Almendros)
1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
ACS Paragon Plus Environment
Page 2 of 34
Page 3 of 34
Journal of Agricultural and Food Chemistry
1
Environmental information provided by free lipids in soil samples collected from control and
2
disturbed plots (Madrid, Spain), was assessed by comparing molecular assemblages of terpenoids
3
and distribution patterns of alkanes and fatty acids (FAs), analyzed by gas chromatography-mass
4
spectrometry (GC-MS). Wildfires in pine forests led to increased proportions of retene,
5
dehydroabietin and simonellite. Friedo-oleananes were characteristic in soils under angiosperms,
6
and norambreinolide-type diterpenes in soils encroached by Cistus bushes. Steroids were major
7
compounds in pastured site. Enhanced Shannon’s lipid biodiversity index in disturbed soils
8
compared to control soils suggested patterns of recent lipids overlapped with a preserved original
9
lipid signature. The extent of the environmental impacts was illustrated as Euclidean distances
10
between paired control and disturbed sites calculated using as descriptors the compounds in alkyl
11
homologous series. As expected, reforestation, bush encroachment, wildfires and cultivation were
12
reflected by changes in the molecular record of lipids in soils.
13 14 15 16
KEYWORDS: Alkane, Biomarker, Terpene, Fatty acid, Molecular tracer, Signature lipid
17
1
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
18
INTRODUCTION
19
The lipid fraction of soil is often considered a valuable source of environmental information,
20
representing a continuous molecular record shedding light on climate change and the intensity of
21
organic matter (OM) turnover1–4. In fact, whereas the soil lipid fraction originally consists of a
22
heterogeneous molecular assemblage inherited mainly from plants and microorganisms5,6, further
23
abiotic or microbial transformation of these biogenic lipids provides additional compounds in the
24
soil7–11. The dynamics of soil lipids is complex since, apart from the above processes, lipid mixtures
25
are subjected to continuous biodegradation of their comparatively labile molecules. Nevertheless,
26
lipid molecules can also be included into organo-mineral structures in progressively transformed
27
soil organic matter pools12–14 with an improved preservation of compounds associated with small
28
size aggregate fractions15. Other fractions of lipids present in soil, which require chemical
29
treatments for their release in the form of free molecules, may correspond both to cellular
30
constituents in not still degraded biomass, and to condensation products incorporated to humic
31
substances. Therefore, the balance between biodegradation and humification processes could be
32
monitored by the minor but diagnostic fraction of lipids temporarily free in the soil, which are the
33
subject of this study. Hence, the surviving lipids which can be directly isolated from soil in the form
34
of free compounds could represent a molecular signature for reconstructing recent and past soil
35
processes, since their occurrence could in most cases depend on the environmental impact on
36
terrestrial ecosystems16–19.
37
Finally, specific soil lipids are also important due to their role on soil processes, acting as
38
antimicrobial agents20–23 or in the alellopathic interaction between higher plants1,24,25; or in insect-
39
host plant relationships 26 as well as through the effect on soil physical properties, mainly
40
aggregation and soil water repellence27–30.
41
In the present study, the aim was the analytical comparison between free soil lipids from relict
42
forests and those from soils in adjacent sites under the same climatic conditions and the same
43
original geological substrate, but affected by environmental perturbation processes typical of 2
ACS Paragon Plus Environment
Page 4 of 34
Page 5 of 34
Journal of Agricultural and Food Chemistry
44
continental Mediterranean ecosystems in central Spain. Both individual signature compounds, as
45
well as the distribution patterns of the major series of alkyl compounds, were analysed in 16 well
46
characterised ecosystems in central Spain, consisting of control forests or sites affected by (i)
47
clearing and bush encroachment, (ii) cultivation, (iii) wildfire and (iv) reforestation with pine. The
48
specific objectives would be to assess the extent of the biogeochemical changes undergone by the
49
soils due to the environmental perturbation, to test the response of the biomarker assemblages to
50
environmental changes and to compare hypothetical new proxies for these perturbations.
51 52
MATERIAL AND METHODS
53 54 55
Sampling Soil samples from 8 contrasting continental Mediterranean forest ecosystems in Madrid (central
56
Spain) representative for sclerophyllic (oak), mesophylic (chestnut and ash) or pine forests were
57
collected. In addition to the control sites, another set consisting of 8 altered neighbor soils with the
58
same climatic, geologic and topographic features were sampled31. A two character labelling code
59
was used to refer to the sampling sites: the first letter (see below) was the code for the soil
60
series, followed by an odd number in the case of the relict ecosystems and an even number for
61
the perturbed ones. The distance between sampling points in the paired soils always was < 100 m.
62
Some general features of the soil samples are listed in Table 1.
63
Sampling was carried out in duplicate between April and June: samples were collected after
64
removing the litter layer, when existing, and the soil material (the whole O horizon) was collected
65
with a spade. In order to obtain representative samples averaging the possible spatial variability in
66
the plots, composite samples were taken from each plot. Each individual sample was prepared by
67
mixing three soil subsamples of ca. 1 L from the points of a virtual triangle of ca. 100 m side. The
68
samples were air-dried, litter and root fragments were hand-picked and the resulting material was
69
sieved to 2 mm and used to determine routine analytical characteristics, calculating the least 3
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
70
significant difference illustrating spatial variability in soil taxonomic characteristics. A combined
71
subsample for the molecular characterization of the lipid fraction was prepared by mixing the
72
subsamples.
73
The characteristics of the 16 ecosystems (8 control sites, coded with a label including an odd
74
numeral + 8 disturbed sites, coded with even numerals) are shown in Table 2. Series R corresponds
75
to soil areas covered by the original oak forest (Quercus ilex subsp. ballota: samples R1 and R3). At
76
some places, the oak had been removed and the site reforested with Pinus pinea or Pinus pinaster
77
(R2 and R4, respectively). The series of cleared forests corresponded to oak (C1) or ash (Fraxinus
78
angustifolia, C3) forest dedicated to cereal cultivation (C2) or grazing pasture (C4). In other cases,
79
oak forests were cleared in a historical period for wood extraction (B1, B3), then occupied by
80
Mediterranean bush with Cistus (B2) or Cytisus (B4). The effect of wildfire was examined by
81
comparing forests of Pinus halepensis (F1) or Pinus sylvestris (F3) with the adjacent sites affected
82
by high intensity (F2) or medium intensity (F4) forest fires.
83 84
Lipid analysis
85
Soil samples (50-g) were Soxhlet extracted with petroleum ether (40–60 °C) for 24 h. This
86
solvent was chosen in order to prevent the removal of dark-colored nonvolatile macromolecular
87
material, probably oligomers of humic substances or lignins, that is co-extracted when using more
88
polar solvents3. Preliminary experiments using polar solvents such as CH2Cl2-MeOH produced, in
89
the case of forest soils, very dark brown extracts not suitable for direct GC analysis. The extract was
90
filtered and concentrated in a rotary evaporator to approximately 50 mL, dried under a stream of
91
N2 at room temperature (20–25 °C) and the residue weighed and methylated with CH2N2/Et2O.
92
The lipids were separated and identified using GC-mass spectrometry (GC-MS) with an HP 5890
93
gas-chromatograph connected to an HP 5971 mass detector (EI, 70 eV). A cross-linked OV-1
94
column (25 m × 0.25 mm i.d. × 0.25 µm film thickness) with and He flow of 1 mL min-1 was used.
95
Constant pressure, and split injection were used. The oven temperature was programmed from 70 4
ACS Paragon Plus Environment
Page 6 of 34
Page 7 of 34
Journal of Agricultural and Food Chemistry
-1
96
°C to 320 °C (held 20 min) at 4 °C min . Injection temperature was 300 °C. Some splitless
97
injections were also carried out for trace compound identifications; no significant differences in
98
quantitative results were found between both injection techniques. Total ion current (TIC) areas
99
(acquisition range 39–530 amu) were used for quantitative measurements. Compounds were
100
assigned from their electron impact mass spectra and confirmed, when possible, by comparison
101
with those in the spectral databases or with standards of the authentic compounds.
102
Identification was confirmed by using retention indices calculated for all the compounds from
103
their retention times and those of the n-alkane series, in programmed-temperature conditions.
104
Coelution caused overlapped peaks within FA methyl ester and alkanol series, which were
105
determined from their single ion characteristic traces using the ions at m/z 74 and 69,
106
respectively. Integration values of the peaks in the traces for these single ions were transformed
107
into total ion counts by using suitable correction factors that were calculated from the
108
fragmentation patterns of pure compounds32.
109 110 111
RESULTS and DISCUSSION The total lipid concentration ranged between 0.05 and 9.50 g kg-1 (Table 2). Fig. 1 illustrates the
112
composition of the petroleum ether extract from some representative samples. About 100
113
compounds in each sample, mainly n-alkanes, n-alkanols and FAs were tentatively assigned (Table
114
3). Depending on the soil, variable proportions of cyclic compounds mainly monoterpenes,
115
sesquiterpenes, diterpenes, triterpenoids and steroids, were present (Figs. 2 and 3).
116
Many of the identified compounds could be considered as typical biomarkers, to the extent
117
that they have definitive chemical structures, which can be related directly or indirectly through a
118
set of structural alterations to biogenic sources, and they cannot be synthesized by abiogenic
119
processes7. Typically, monoterpenes and cyclic diterpenes (the latter as hydrocarbons, acids,
120
alcohols, or other derivatives) are major constituents of gymnosperms resins33–40. Other
121
compounds as sterols may derive mainly from the cell walls of animals and plants whereas 5
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
122
Page 8 of 34
37
triterpenoids are widely distributed in the epicuticular constituents of higher plants .
123 124
Cyclic compounds
125
Table 4 and Fig. 2 show the compounds with a terpenoid structure. They amounted to between
126
0.9 and 51.9% as proportions of the total ion chromatographic area in the different samples.
127
Cyclic diterpenes (as hydrocarbons, acids, alcohols, or other derivatives) were the most
128
abundant cyclic compounds (0.9–45.2%, Table 5) and belonged to three families: abietane,
129
pimarane and labdane (Fig. 2). The abietane-type diterpenoids were the most abundant, mainly
130
in samples derived from Pinus pinea (R2: 22.3%) and Pinus halepensis (F1: 25.3%).
131
Dehydroabietic acid (compound 95, Table 4), was the most frequent diterpenoid acid in all
132
samples. In fact, this compound is the resin acid most common in the geosphere where it
133
originates from the rapid conversion of abietic acid41–42.
134
There was also an appreciable abundance of abietic acid (compound 101), occurring in
135
samples R2 and F1, and levopimaric acid (99) in samples R2, F1 and F4. Seco-dehydroabietic
136
acid33 which appeared as two isomers (1α– and 2β–) was frequent in the soil lipids under pine
137
forest. In addition, samples R2, R4, F1 and F4 had appreciable proportions of alteration products
138
of dehydroabietic acid, such as 7-oxodehydroabietic (110), 7-hydroxydehydroabietic (111) and
139
15-hydroxydehydroabietic acids (112).
140
A series of abietane hydrocarbons, also presumably derived from the above acids, was also
141
present. The highest proportion of these hydrocarbons corresponded to lipids from soils under
142
burned forest, F2. Dehydroabietin was present in all the soils under pine43 but other typical
143
hydrocarbons, formed by transformation of resin acids8, were also present, such as
144
dehydroabietane (58), 19-norabieta-4,8,11,13-tetraene (40), simonellite (37) and retene (31).
145 146
The pimarane-type diterpenoids (Fig. 2) are traditionally considered to be less stable than the abietanes and, consequently, were present in lower proportions, between 0.1 and 18.5% (Table
6
ACS Paragon Plus Environment
Page 9 of 34
Journal of Agricultural and Food Chemistry
147
5). Pimaric acid (96), was frequent in all samples under pine forest, but occurred only as traces in
148
samples C3, C4 and B4.
149
Sandaracopimaric (98), 8,15-isopimaradien-18-oic acid (97) and isopimaric (101) acids were
150
identified. The latter was in high proportion (5.9%) in R2, under Pinus pinea. The labdane
151
diterpenes (Fig. 2) were not abundant (0.1–4.8%) appearing only in R4, B4, F1, F3 and F4. The
152
highest proportion of labdanes occurred in B4, under Cistus bush. This sample also contained
153
eperuic acid (104), cativic acid (105), labdanolic acid (113) and norambreinolid (36). The major
154
compound was labdanolic acid (1.4%). Other labdane-type diterpenoids present in soils under
155
pine were manoyl oxyde (73), pinifolic acid (126) and anticopalic acid (103).
156
As expected, the highest proportion of the above resin acids occurred in the soils under pine
157
(R2, R4, F1, F2, F3 and F4). The presence of small relative amounts of these compounds in C3
158
(ash forest), C4 (pasture) and B2 (Spanish broom) could be explained by wind transport of
159
aerosols from adjacent pine forests, which would contribute to the fact that terpenoids are likely
160
over represented in the sedimentary record43. Nevertheless, "classical" conifer biomarkers of the
161
abietane series have also been reported to be widespread in several cyanobacterial strains44. In
162
any case, in environmental conditions abietanes readily turn into dehydroabietic acid, and have
163
no tendency to accumulate in the advanced transformation stages of soil lipids3.
164
Decarboxylation and aromatization are considered to be trends during ageing of resin acids8
165
leading to diterpene hydrocarbons at the final transformation stages. These hydrocarbons, with
166
the highest proportion in F2, would reflect the effect of fires in producing decarboxylation and
167
aromatization reactions45,46.
168
The fire-affected soils showed decreased proportions of monoterpenes, sesquiterpenes and
169
diterpenoids and concomitantly increased proportions of diterpene hydrocarbons. In fact,
170
wildfires in pine forests would lead to thermal dehydrogenation and decarboxylation of
171
diterpene resin acids47 resulting into aromatic hydrocarbons such as retene, dehydroabietin and
172
simonellite. Sample F4 showed an increased proportion of diterpenes (Table 5), as could 7
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 10 of 34
45
173
correspond to thermal distillation of plant waste incorporated into soil . In particular
174
sesquiterpene patterns were quite responsive to the intensity of the fire (F2 > F4).
175
Triterpenoids (2.1–3.9% of the total chromatographic area) were present in R4, C4, B1, B4, F2
176
and F4 (Table 5 and Fig. 2). This group included typical sterols (27, 28 or 29 C atoms). Sample C4
177
displayed the greatest number of sterols, the main ones being cholesterol, campesterol and β-
178
stigmasterol. In R4 the major sterol was stigmast-5-en-3β-ol and in F4 a steroidal ketone
179
(stigmast-4-en-3-one) was found. Other triterpenoids had a pentacyclic skeleton34,35. Friedelan-
180
3-one (friedelin) was detected (Fig. 2) in B1 and B4, whereas friedo-olean-14-en-3-one and
181
ursolic acid were found in F2.
182
Mainly F1 (Pinus halepensis) contained monoterpenes (4.7%); α-pinene and β-pinene, being
183
the most abundant (2% each). Myrcene, γ-terpinene, terpinolene, α-terpinene and p-cymenol
184
(Table 4, Figs. 1 and 2) were also identified. Several sesquiterpenes were also present in F1
185
(Table 4, Fig. 2), mainly caryophyllene (14, 0.4%) and caryophyllene oxide (23, 1.1%). Other
186
sesquiterpenes present in lower proportions were cadalene, calamenene, α-humulene, α-
187
cubebene, α-copaene, α-muurolene, δ-cadinene, β-eudesmol and cembrene (9, 13, 15–19, 24
188
and 59, respectively). No sesquiterpenes were found in other samples, although the lipid
189
fraction of R4 (Pinus pinaster) included a trace of cadalene (0.1%). The fact that monoterpenes
190
were not major compounds in the samples may be explained by the fact that monoterpenes
191
have a relatively low boiling point and probably they are not strongly retained in Mediterranean
192
soils48. In addition, the presence of oxygen-containing functional groups as well as unsaturated
193
bonds in lipid molecules makes lipids suitable for strong retention in the soil organo-mineral
194
matrix, including physical occlusion and hydrophobic interactions3. In active terrestrial
195
ecosystems the occurrence and intensity of such processes should also be related to the
196
selective preservation of free lipid assemblages.
197
Concerning sterols, the dominance of cholesterol, of animal origin, in the pastured site C4 agree
198
with the current use of this soil. 8
ACS Paragon Plus Environment
Page 11 of 34
199
Journal of Agricultural and Food Chemistry
Finally, in two samples (R4 and B1) there were detectable proportions of DDE (1,1,1-trichloro-
200
2,2-bis(p-chlorophenyl)ethylene), derived from the degradation of
201
dichlorodiphenyltrichloroethane (DDT).
202 203
Alkyl series
204
Figure 1 and Table 3 show that the major lipid compounds corresponded to alkyl compounds:
205
alkanes (0.5–68.1%), FAs (4.9–55.9%) and alkanols (5.0–53.5%). Fig. 3 illustrates some distribution
206
patterns (Cno vs. relative abundance) of the n-chain homologues of the major series.
207
Homologous series of n-alkanes, n-alkanols and n-alkanoic acids also show differential
208
characteristics depending on their origin. In general the homologues > C20 with strong even-to-
209
odd C-number preference are often considered as characteristic of wax from vascular plants,
210
whereas those < C20 are believed to derive mainly from microorganisms8. In particular the α,ω-
211
alkanedioic acids and ω-hydroxyacids (C12, C14 and C16) are constituents of plant polyesters such
212
are cutin and suberin37. Other compounds as the isoprenoid ketone (6,10,14-
213
trimethylpentadecan-2-one) probably derive from phytol49. The n-alkanes (Fig. 3) had in general
214
a clear dominance of homologues > C20 (Table 3), with frequent maxima at C29 and C31, and an
215
odd/even carbon preference index (CPI) typical of epicuticular wax in vascular plants36.
216
It is often assumed that soil alkanes < C20 may have a microbial origin8. They had the highest
217
relative abundance in R2 and R4 from conifer vegetation. Alkane abundance in these samples
218
showed a trend towards a bimodal distribution (maxima at ca. C23 and C29 or C31) that could be
219
interpreted as caused by a dual origin for the homologues in the series.
220
In some distribution patterns the difference between soil pairs was small; this could indicate a
221
comparatively low environmental impact, i.e. a low intensity of the perturbation, or high soil
222
resilience, or both. Visual inspection of the histograms (Fig. 2) showed major changes in soils
223
affected by fire, which presumably produce alkane fragmentation as indicated by independent
224
studies in laboratory-controlled conditions45. The extent of this fragmentation could be reflected 9
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 12 of 34
50
225
by the lack of an odd/even preference in the low molecular weight (MW) (< C20) hydrocarbons .
226
This production of a variety of short chain compounds at pyrolysis temperatures could be favored
227
by their fixation into mineral matrices, which hinder their thermo-evaporation51.
228
In any case, the proportions of alkanes and FAs did not show the same response to fire (Table
229
3), which is expected from the very different intensity of wildfires affecting soils F3 and F4 (Table
230
1).
231
The alkane patterns of the soils under deciduous forest (oak, chestnut) showed comparatively
232
small changes after clearing and cultivation, or after bush encroachment (series C or B), but
233
reforestation with pine was associated with accumulation of homologues of comparatively lower
234
Cno (Table 3).
235
In the studied soils, the classical indices based in the average chain length or the carbon
236
preference (Table 3), reflected mainly local differences between neighbor ecosystems rather
237
than systematic environmental changes. General differences were examined taking advantage of
238
additional indices such are the diversity indices of the whole alkyl series, as a possible proxy of
239
the complexity of the whole ecosystem52 or the Euclidean distances between paired control and
240
perturbed soils, to provide an objective index of the extent of the corresponding change in the
241
whole alkyl series3.
242
In fact, the extent of the changes in the alkane assemblages in samples from neighbor (relict
243
and perturbed) plots could be quantified by way of a statistical dissimilarity calculated using the
244
proportions of the different homologues in the series from the adjacent soil samples (e.g. R1 vs.
245
R2, R3 vs. R4, etc). For this purpose3, we selected the Euclidean distance ED:
246
EDi , j =
n
∑ (c − c ni
nj
)2
i =1
247
where Cn corresponds to the total abundance (normalized by the total peak area of all
248
homologues) of each compound in the series, and i and j are the different soils compared. The
249
resulting distances EDi,j which are shown in Table 3 provide a relative measure of the extent of the 10
ACS Paragon Plus Environment
Page 13 of 34
250
Journal of Agricultural and Food Chemistry
change in the whole alkyl series.
251
The interpretation of the differences between neighbor paired sites (Table 3) confirms the
252
substantial alteration produced by intense fire (F1 F2) compared with the case in which fire
253
affected more the forest canopy than the epipedon (F3 F4). On the other hand, the lowest
254
change was observed in the C1 C2 pair, which suggested as a high stability of the OM in the
255
centennial oak forest. This was also probably the case with the cleared oak forest after
256
encroachment by ameliorant bush (B1 B2) which could be interpreted as a low capacity of the
257
weak density secondary Mediterranean vegetation to introduce change in the historical period, in
258
the biogeochemical behavior of the original system31.
259
While differences between paired samples can be estimated from the Euclidean distances, the
260
complexity of the composition of a given sample could be quantitatively estimated using diversity
261
indices. The Shannon-Wiener index takes into account the number and size of subgroups in a given
262
population. It was calculated for each alkyl series using the relative abundances of the different
263
homologues; its value should increase with the complexity of the molecular population. Table 4
264
and Figs. 2 and 3 clearly show that the diversity index (H’) in the perturbed sites tended to be
265
higher than in the original ones. This trend to increased diversity was observed mainly for the
266
alkanes and, to lesser extent, FAs. This behaviour was observed in all types of perturbations,
267
including series B and C, often considered as leading to desertification or to a simplification of the
268
ecosystem structure. The increase in H’ values occurred both in situations in which the total
269
abundance of lipids in soil decreased after the perturbation (most cases) or when it increased
270
(series B, under Mediterranean bush with essential oils, adapted to xerophytic conditions).
271
There are two probable reasons for the enhanced molecular diversity after perturbation. First,
272
selectively preserved molecules from the original ecosystem would coexist with those recently
273
incorporated from the new vegetation. Other explanation could be the abiotic transformations and
274
microbial reworking of the original lipid assemblage in perturbed soils51. Since increasing molecular
275
complexity has been considered to parallel the extent to which the composition of soil lipid in 11
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
3
276
forest ecosystems changes as regards that of the original plant lipid , this index could inform on
277
the transformation stages of the soil system.
278
The FA patterns showed a clear even/odd preference, with C16, C22 and C24 as major members.
279
Comparison of the FAs from soils under oak forest (R1, R3) with those from soils reforested with
280
Pinus pinea and Pinus pinaster (R2, R4) showed that the former series could be differentiated by
281
maxima at C22 and C24 and a higher chain length ratio (> C20 /< C20) (Table 3). Typical C15, C17 and C19
282
iso- and anteiso-branched FAs often considered as markers of microbial activity 53 were found in
283
all samples, except B1.
284
Like the alkanes, the differences between FA series from paired plots were quantified via
285
Euclidean distances. Large transformations, in both cases accompanied by the relative
286
accumulation of low MW FAs were found after reforestation and after wildfire. There were small
287
changes in the sites cleared for cultivation or for pasture.
288
The presence of unsaturated acids was also frequent in the samples, except B1, C1 and C2, the
289
most frequent being palmitoleic (C16:1) and oleic (C18:1). In general unsaturated FAs present low
290
abundances and they could be considered as indicative of recent biogenesis, since such
291
compounds are comparatively reactive and are rapidly degraded or condensed with the active sites
292
of the soil matrix54.
293
In general terms, the information supplied by the study of FAs was similar to that suggested
294
by the alkane series. To a large extent, soil FA composition is considered to reflect the original
295
vegetation type55. All samples (Table 3) showed dominance of compounds with an even number
296
of carbons (even/odd >3) typical of a recent biogenic origin8,36. As in the case of the alkanes, for
297
FAs from series F the samples affected by fire showed a higher proportion of homologues < C20,
298
as would correspond to the thermal breakdown of higher chain homologues45,50.
299
When the impact of the disturbation was assessed in terms of the Euclidean distance between
300
FA patterns (Table 3), it is observed small changes in the case C1 C2 , which could be interpreted
301
as a substantial resilience of the relict oak forest, with a high clay content and thick humus horizon 12
ACS Paragon Plus Environment
Page 14 of 34
Page 15 of 34
Journal of Agricultural and Food Chemistry
302
(10 cm). As in the case of alkanes, and when comparing plots after bush encroachment, the most
303
intense changes were produced in the soil occupied by Cistus (B4). Finally, the most marked
304
perturbation of the whole pattern of FAs corresponds, as in the case of alkanes, to the soil affected
305
by severe wildfire (F2). In these soils subjected to heating, the changes (mainly related to a
306
selective accumulation of short chain FAs, as in the case of alkanes), were proportional to the
307
intensity of the fire (F2 > F4).
308
As for the FAs, the n-alkanols (C13–C28) showed a clear odd/even preference, with a maximum at
309
C24 or C26. Besides the n-alkanols, there was a less abundant series of branched alcohols (mainly
310
C24, C26 and C28). The C14 and C16 alkanediols, frequent constituents of epicuticular lipids of
311
higher plants37 also occurred in R2. Nevertheless,the alkanol patterns did not show systematic
312
changes depending on the perturbation in the relative chain length, or in the carbon preference
313
index, but provided complementary local information on features that were not straightforwardly
314
reflected by the above alkyl compounds. In particular, the effects of bush encroachment (B2, B4)
315
could be distinguished by the increased abundance of the low MW homologues. In series R, and
316
after reforestation with pine, a dominance of the C24–C28 alkanols was observed, but in F2 and
317
F4 the chain length was shorter as a probable effect of fire.
318 319
Environmental changes reflected by the lipid signature
320
As a whole, the results suggest that the lipid molecular record in soils from Mediterranean
321
ecosystems is highly responsive to the environmental perturbations. The finding that diversity
322
indices of the alkane and FA series increase in disturbed soils vs. control soils could be
323
interpreted in terms of supposing that preexisting compounds from the original lipid signature
324
coexist with new lipids generated after the perturbation. In general, these alkyl compounds
325
were found to be an important source of environmental proxies when described via suitable
326
indices calculated from the whole homologue series.
327
Of the four pine forests, the biomarkers investigated remain efficient for plant source 13
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
328
reconstruction in spite of the environmental perturbations: the most stabilized ecosystems
329
showed comparatively higher proportions of dehydroabietic acid and lower of abietane- and
330
pimarane-type diterpenoids. Soils affected by forest fires showed increased proportions of
331
shorter chain alkanes and FAs as well as of natural polycyclic aromatic hydrocarbons such as
332
retene, dehydroabietin and simonellite.
333
On the other hand, some individual signature compounds were specifically found in a
334
reduced number of samples: friedooleananes were found in three soils under angiosperm
335
vegetation, whereas norambreinolid-type diterpenes occurred in soil under Cistus bush. The
336
steroids were found to prevail in a pastured site. In conclusion, the comparison of free lipid compounds in soils plots whose characteristics
337 338
have changed as a result of disturbances seems to indicate as if there were a "memory of soil
339
lipids" that revealed environmental impact in soil. In all cases, the signature of lipid assemblages
340
reflected changes in land use with alkyl homologous series being especially responsive to the
341
impact of fires. From the chemometric viewpoint, the environmental impact was associated with
342
enhanced diversity of alkane homologues, and Euclidean distance between alkane abundances
343
in paired control and disturbed sites acts as a proxy for the extent of soil disturbation.
344 345
ACKNOWLEDGMENTS
346
Financial support by Spanish MINECO project CGL2013-43845-P is acknowledged. The authors
347
would like to express their sincere appreciation to Dr. J.R. Maxwell by his helpful comments and
348
suggestions to a previous version of the document.
349 350
REFERENCES
351
(1)
Wiley: New York, 1982, 172–194.
352 353
Stevenson, F.J. Soil Lipids. In Humus Chemistry: Genesis, Composition, Reactions, Chap. 7,
(2)
Eglinton, T.I.; Eglinton, G. Molecular proxies for paleoclimatology. Earth Planet. Sci. Lett., 14
ACS Paragon Plus Environment
Page 16 of 34
Page 17 of 34
Journal of Agricultural and Food Chemistry
2008, 275 (1-2) 1–16.
354 355
(3)
Mediterranean forests. Eur. J. Soil Sci. 1996, 47, 183–196.
356 357
Almendros, G.; Sanz J.; Velasco, F. Signature of lipid assemblages in soils under continental
(4)
Almendros, G.; Tinoco, P.; Rodríguez-Sánchez, S.; Jiménez-González, M.A.; Sanz, J.
358
Comparing chromatographic data for soil lipid compounds as density surfaces in the space
359
defined by their atomic ratios. In Advances in Chromatography and Related Techniques,
360
(Book of Abstracts of the XVI Scientific Meeting), González-Pérez, J.A.; Almendros, G.;
361
González-Vila, F.J.; De la Rosa, J.M. Eds.; Spanish Society of Chromatography and Related
362
Techniques: Seville, 2016, 163.
363
(5)
pp. 45–53.
364 365
Hegnauer, R. Chemotaxonomie der Pflanzen, Birkhauser: Basel und Stuttgart, 1966. Vol. 4.
(6)
Zocatelli, R.; Jacob, J.; Gogo, S.; Le Milbeau, C.; Rousseau, J.; Laggoun-Défarge, F. Spatial
366
variability of soil lipids reflects vegetation cover in a French peatland. Org. Geochem.
367
2014, 76, 173–183.
368
(7)
Simoneit, B.R.T.; Pisani, O.; Ekpo, B.O.; Fubara, E.P.; Nna, P.J.; Ekpa, O.D. Lipid biomarker
369
analysis of suspended particulate matter from the Great Kwa River, SE Nigeria: origins and
370
environmental implications of biogenic and anthropogenic organic compounds. Aquat.
371
Geochem. 2017, 23, 89–108.
372
(8)
Simoneit, B.R.T.; Mazurek, M.A. Organic matter of the troposphere–II. Natural background
373
of biogenic lipid matter in aerosols over the rural Western United States. Atmos. Environ.
374
1982, 16; 2139–2159.
375
(9)
Amblès, A.; Jacquesy, R.; Tejada, A.; Jambu, P.; Fustec, E. Origin and dynamics of lipids in a
376
podzol soil. In Studies about Humus. Transactions of the VIIth International Symposium
377
(Novak, B., Ed.; Institute of Crop Protection: Prague, 1983, Vol. I, 5–8.
378 379
(10) de Leeuw, J.W. Molecular fossil lipids as indicators for sources of input and short–term diagenesis. Stuifmail 1988, 5, 18–29. 15
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
380 381 382
(11) Eglinton, G.; Parkes, R.J.; Zhao, M. Lipid biomarkers in biogeochemistry: Future roles?. Mar. Geol. 1993, 113 (1–2) 141–145. DOI: 10.1016/0025-3227(93)90155-O. (12) González-Vila, F.J.; Almendros, G.; Martin, F. An evaluation of the differences in the
383
composition of humic acids in soils under oak and pine forests by GC–MS after mild
384
degradation. Plant Soil 1987, 103, 83–88.
385 386 387 388 389 390
(13) Amblès, A.; Jacquesy, J.C.; Jambu, P.; Joffre, J.; Maggi-Churin, R.; Polar lipid fraction in soil: a kerogen–like matter. Org. Geochem. 1991. 17, 341–349. (14) Spaccini, R.; Baiano, S.; Gigliotti, G.; Piccolo, A. Molecular characterization of a compost and its water-soluble fractions. J. Agric. Food Chem. 2008, 56, 1017–1024. (15) Jansen, B.; Wiesenberg, G.L.B. Opportunities and limitations related to the application of plant-derived lipid molecular proxies in soil science. SOIL 2017, 3, 211–234.
391
(16) Stevenson, F.J. Lipids in soil. Amer. Oil Chemists' Soc., 1966, 43, 203–210.
392
(17) Bailly, J.R. Quelques aspects de la biochimie des lipides du sol, In I Lipidi Nelle Piante e nel
393
Terreno, XIII Simposio Internazionale di Agrochimica. Supplementary issue 1, Agrochimica,
394
1985; 64–92.
395
(18) Dinel, H.; Schnitzer, M.; Mehuys, G.R. Soil lipids: origin, nature, content, decomposition,
396
and effect on soil physical properties. In Soil Biochemistry; Bollag, J.-M.; Stotzky, G. Eds.;
397
Dekker, New York, 1990; Vol. 6, 397–429.
398
(19) Guillemot, T.; Bichet, V.; Gauthier, E.; Zocatelli, R.; Massa, C.; Richard, H.; Environmental
399
responses of past and recent agropastoral activities on south Greenlandic ecosystems
400
through molecular biomarkers. Holocene 2017, 27, 783–795.
401
(20) Dhar, P.; Chan, P.; Cohen, D.T.; Khawam, F.; Gibbons, S.; Snyder-Leiby, T.; Dickstein, E.; Rai,
402
P.K.; Watal, G. Synthesis, antimicrobial evaluation, and structure-activity relationship of
403
α-pinene derivatives. J. Agric. Food Chem. 2014, 62, 3548–3552.
404 405
(21) Cristani, M.; D'Arrigo, M.; Mandalari, G.; Castelli, F.; Sarpietro, M.G.; Micieli, D.; Venuti, V.; Bisignano, G.; Saija, A.; Trombetta, D. Interaction of four monoterpenes contained in 16
ACS Paragon Plus Environment
Page 18 of 34
Page 19 of 34
Journal of Agricultural and Food Chemistry
406
essential oils with model membranes: Implications for their antibacterial activity. J. Agric.
407
Food Chem. 2007, 55, 6300–6308.
408
(22) Xu, J.; Kang, J. ; Cao, X. ; Sun, X. ; Yu, S. ; Zhang, X. ; Sun, H. ; Guo, Y. Characterization of
409
diterpenes from Euphorbia prolifera and their antifungal activities against
410
phytopathogenic fungi. J. Agric. Food Chem. 2015, 63 (25), 5902–5910.
411
(23) Heumann, S.; Schlichting, A.; Böttcher, J.; Leinweber, P. Sterols in soil organic matter in
412
relation to nitrogen mineralization in sandy arable soils. J Plant Nutr. Soil Sci., 2011, 174,
413
576–586.
414
(24) Lynch, J.M. Products of soil microorganisms in relation to plant growth. Crit. Rev. Microbiol. 1976, 5, 67–107.
415 416
(25) Torres-Romero, D.; King-Díaz, B.; Strasser, R.J.; Jiménez, I.A.; Lotina-Hennsen, B.; Bazzocchi,
417
I.L. Friedelane triterpenes from Celastrus vulcanicola as photosynthetic inhibitors. J. Agric.
418
Food Chem. 2010, 5, 10847–10854.
419
(26)
Fraga, B.M.; Díaz, C.E.; Guadaño, A.; González-Coloma, A. Diterpenes from Salvia
420
broussonetii transformed roots and their insecticidal activity. J. Agric. Food Chem. 2005,
421
53, 5200–5206.
422 423 424
(27) De Blas, E.; Almendros, G.; Sanz, J.; Molecular characterization of lipid fractions from extremely water-repellent pine and eucalyptus forest soils. Geoderma 2013, 206, 75–84. (28) Atanassova, I.; Hristov, B.; Shishkov, T.; Doerr, S. Lipid biomarkers and their environmental
425
significance in mine soils from Eastern Europe. Arch. Agron. Soil Sci., 2017, 63, 1–14.
426
(29) Sarker, T.C.; Incerti, G.; Spaccini, R.; Piccolo, A.; Mazzoleni, S.; Bonanomi, G. Linking organic
427
matter chemistry with soil aggregate stability: Insight from 13C NMR spectroscopy. Soil
428
Biol. Biochem. 2018, 117, 175–184.
429
(30) Siteur, K.; Mao, J., Nierop, K.G.J.; Rietkerk, M.; Dekker, S.C.; Eppinga, M.B. Soil water
430
repellency: A potential driver of vegetation dynamics in coastal dunes. Ecosystems, 2016,
431
19, 1210–1224. 17
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
432
(31) Tinoco, P., Caracterización Molecular de la Materia Orgánica de Suelos Afectados por
433
Distintos Tipos de Degradación en la Comunidad de Madrid. PhD Thesis, Universidad
434
Autónoma de Madrid, 2000.
435
(32) Almendros, G.; Hernández, Z.; Sanz, J.; Jiménez-González, M.A.; Rodríguez-Sánchez, S.;
436
González-Pérez, J. A. Graphical statistical approach to soil organic matter resilience using
437
analytical pyrolysis data. J. Chromatogr. A, 2018, 1533, 164–173.
438
(33) Zinkel, D.; Zank, L.C.; Wesolowski, M.R. Diterpene Resin Acids. A Compilation of Infra Red,
439
Mass, Nuclear Magnetic Resonance, Ultra Violet Spectra and Gas Chromatographic
440
Retention Data (of the methyl esters), USDA, Forest Service, Forest Products Laboratory,
441
Madison, Wisconsin, 1971.
442
(34) Nordby, H.E.; McDonald, R.E. Friedelin, the major component of grapefruit epicuticular wax. J. Agric. Food Chem. 1994, 42, 708–713.
443 444
(35)
Szakiel, A., Pa̧ czkowski, C., Huttunen, S. Triterpenoid content of berries and leaves of
445
bilberry Vaccinium myrtillus from Finland and Poland J. Agric. Food Chem., 2012, 60,
446
11839–11849.
447
(36) Stránský, K.; Streibl, M.; Herout V. On natural waxes VI. Distribution of wax hydrocarbons
448
in plants at different evolutionary levels. Coll. Czech. Chem. Commun. 1967, 32, 3213–
449
3219.
450 451
(37) Walton T.J. Waxes, cutin and suberin. In Methods in Plant Biochemistry; Harwood, J. L; Boyer, J. R., Eds.; Academic Press, New York, 1990; Vol. 4; 105–158.
452
(38) Gough, L.J.; Conifer resin constituents. Chem. Industry, 1964, 2059–2060.
453
(39) Barjaktarović, B.; Sovilj, M.; Knez, Ž. Chemical composition of Juniperus communis L. fruits
454
supercritical CO2 extracts: Dependence on pressure and extraction time. J. Agric. Food
455
Chem. 2005, 53, 2630–2636.
456
(40) Zavarin, E.; Snajberk, K. Oleoresins of pinions. J. Agric. Food Chem., 1980, 28, 829–834.
457
(41) Chang, T. L.; Mead, T. E. Mass spectra of diterpene resin acids methyl esters. J. Amer. Oil 18
ACS Paragon Plus Environment
Page 20 of 34
Page 21 of 34
458 459 460 461
Journal of Agricultural and Food Chemistry
Chemists’ Soc. 1971, 48, 455–461. (42) Ekman, R.; Modified resin acids in the reaction zone of Fomes annosus–affected sapwood of Norway spruce. Acta Acad. Aboensis Ser. B. 1979, 39, 1–7. (43) Giri, S.J.; Diefendorf, A.F.; Lowell, T.V. Origin and sedimentary fate of plant-derived
462
terpenoids in a small river catchment and implications for terpenoids as quantitative
463
paleovegetation proxies. Org. Geochem. 2015, 82, 22–32.
464
(44) Costa, M.S.; Rego, A.; Ramos, V.; Afonso, T. B.; Freitas, S.; Preto, M.; Lopes, V.;
465
Vasconcelos, V.; Magalhães, C.; Leao, P.N. The conifer biomarkers dehydroabietic and
466
abietic acids are widespread in Cyanobacteria, 2016, Sci. Rep. 2016, 6, 23436.
467
(45) Almendros, G.; González-Vila, F.J. Wildfires, soil carbon balance and resilient organic
468
matter in Mediterranean ecosystems. A review. Spanish J. Soil Sci. 2012, 2, 8–33.
469
(46) González-Vila, F.J.; Tinoco, P.; Almendros, G.; Martin, F. Pyrolysis-GC-MS analysis of the
470
formation and degradation stages of charred residues from lignocellulosic biomass. J.
471
Agric. Food Chem. 2001, 49, 1128–1131.
472
(47) Robinson, N.; Evershed, R.P.; James Higgs, W.; Jerman, K.; Eglinton, G. Proof of a pine
473
wood origin for pitch from Tudor (Mary Rose) and Etruscan shipwrecks: Application of
474
analytical organic chemistry in archaeology. Analyst 1987, 112, 637–644.
475 476 477
(48) Thorin, J.; Nommik, H. Monoterpene composition of cortical oleoresin from different clones of Pinus sylvestris. Phytochemistry 1974, 13, 1879–1881. (49) Brooks P.W.; Maxwell J.R; Patience R.L. Stereochemical relationships between phytol and
478
phytanic acid dihydrophytol and C18 ketone in recent sediments. Geochim. Cosmochim.
479
Acta. 1978, 42, 1175–1180.
480
(50) Almendros, G.; González-Vila, F. J.; Martin, F. Fire–induced transformation of soil organic
481
matter from an oak forest. An experimental approach to the effects of fire on humic
482
substances. Soil Sci. 1990, 149, 158–168.
483
(51) De la Rosa, J.M.; Merino, A.; Jiménez Morillo N.T.; Jiménez-González, M.A.; González19
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
484
Pérez. J.A.; González-Vila, F.J.; Knicker, H.; Almendros, G. Unveiling the effects of fire on
485
soil organic matter by spectroscopic and thermal degradation methods. In: Fire Effects in
486
Soil Properties. Current Knowledge and Methods Used. CSIRO Publishing: Clayton South
487
VIC 3169, Australia, 2018 (in press).
488 489 490
(52) Jiménez-González, M.A.; Álvarez, A.M.; Hernández, Z.; Almendros, G. Soil carbon storage predicted from the diversity of pyrolytic alkanes. Biol. Fertil. Soils. 2018, 54, 617–629. (53) Kaneda, T.; Iso- and anteiso-fatty acids in bacteria: biosynthesis, function, and taxonomic
491
significance. Microbiol. Reviews 1991, 55, 288–302. https://doi.org/10.1007/s00374-018-
492
1285-6.
493 494 495
(54) Moucawi, J.; Fustec, E.; Jambu, P.; Amblès, A.; Jacquesy, R. Biooxidation of added and natural hydrocarbons in soils: Effect of iron. Soil Biol. Biochem. 1981, 13, 335–342. (55) Zhang, J.; Wang, M.; Liang, P.; Cao, Y.; Cao, Z.; Wong, M.H.; Christie, P.; Wu, S., Effects of
496
land use change on soil organic carbon sources and molecular distributions: 6280 years of
497
paddy rice cropping revealed by lipid biomarkers. J. Soils Sediments, 2017, 18, 1–12.
498
20
ACS Paragon Plus Environment
Page 22 of 34
Page 23 of 34
499
Journal of Agricultural and Food Chemistry
Figure captions
500 501
Fig. 1. TIC trace of methylated soil lipids from Quercus ilex subsp. ballota (R1), Pinus pinea (R2),
502
Pinus halepensis (F1) and Pinus halepensis affected by forest fire (F2). FA, fatty acids; AL, alkanes;
503
OH, alkanols; numbers on the peaks refer to Table 4.
504 505
Fig. 2. Monoterpenes, diterpenes and sesquiterpenes (as methyl esters), triterpenoids and
506
sterols identified in soil lipids.
507 508
Fig. 3. Distribution patterns (Cno vs. relative abundance) of selected series of alkyl compounds
509
(for information on the remaining series see Table 3). The n-alkanes indicate change as results of
510
wildfires (F1F2 and F3F4). The FAs indicate changes due to clearing and cultivation [oak
511
forest to cereal fields (C1C2) or ash forest to grazing pasture (C3C4)]. The alkanol series
512
reflect changes attributable to bush encroachment by species of Fabaceae (broom) or Cistaceae
513
(B2 and B4 respectively). Sample labels refer to Table 1. H’, diversity index; EDS1,S2, Euclidean
514
distance between S1 and S2, calculated to quantify the extent of the change after perturbation;
515
the total abundances (normalized to the total peak area of all homologues) of alkanes or FAs were
516
used as descriptors.
517 518 519
21
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 24 of 34
Table 1 Soil classification and general analytical characteristics of samples soils in Central Spain Ref. Ecosystem R1 R2 R3 R4
Relict evergreen oak forest Reforested pine forest Relict evergreen oak forest Reforested pine forest
C1
Relict oak forest
C2
Cultivated wheat field
C3
Relict ash forest
C4
Pastured site
B1
Relict chesnut forest
B2
Secondary bush
B3
Control evergreen oak forest
B4
secondary bush
F1
Pine forest
F2
Burned pine forest (high intensity fire)
F3
Pine forest
F4
Burned pine forest (mainly affecting forest canopy)
Soil type (IUSS Working Group WRB, 2015) Eutric Cambisol (Loamic, Humic ) Eutric Cambisol (Loamic, Humic) Cambic Folic Phaeozem (Loamic, Humic) Eutric Cambisol (Loamic, Humic) Eutric Folic Cambisol (Loamic, Humic) (Rendzic Leptosol (Loamic) Eutric Cambisol (Loamic, Humic) Eutric Cambisol (Loamic, Humic) Eutric Folic Cambisol (Loamic, Humic) Dystric Cambisol (Loamic, Humic) Eutric Cambisol (Loamic, Humic Eutric Cambisol (Loamic, Humic Calcaric Cambisol (Loamic, Humic) Calcaric Cambisol (Loamic, Humic) Cambic Umbrisol (Loamic, Hyperdystric Hyperdystric Cambisol (Loamic, Humic)
Universal Transverse Parent rock Vegetation Mercator (UTM) Quercus ilex subsp. ballota in addition to Rosmarinus officinalis, Daphne gnidium, 4488–411 Granite Cistus ladanifer Pinus pinea in addition to Cistus ladanifer 4488–411 Granite 4490–403
Granite
Quercus pyrenaica, in addition to Prunus spinosa, Rosa canina
4490–403
Granite
Pinus pinaster, in addition to Cistus ladanifer
4472–476
Limestone
Quercus ilex subsp. ballota
4472–476
Limestone
Triticum aestivum
4526–451
Granite
Fraxinus angustifolia, in addition to Rosa canina, Paeonia coriacea
4526–451
Granite
Poa bulbosa, Trifolium dubium, Trifolium campestre in addition to Fraxinus angustifolia, Rosa canina, Paeonia coriacea and Micropyrum tenellum
4465–371
Granite
Castanea sativa, Genista sp., Rosa canina, Vicia sp. and Poaceae
4465–371
Granite
Cytisus scoparius, Genista sp., Retama sphaerocarpa and Lavandula stoechas
4493–401
Granite
Quercus ilex subsp. ballota, Cistus ladanifer, Rosmarinus officinalis and Daphne gnidium
4493–401
Granite
Cistus ladanifer, Daphne gnidium, Retama sphaerocarpa, and residual forest of Quercus ilex subsp. ballota,.
4459–450
Limestone
Pinus halepensis, Eryngium campestre and Reseda sp.
4459–450
Limestone
Pinus halepensis, Eryngium campestre and Reseda sp.
4553–52
Gneiss
Pinus sylvestris, Erica sp., Pteridium aquilinum and Cytisus scoparius
4553–452
Gneiss
Pinus sylvestris, Erica sp., Pteridium aquilinum and Cytisus scoparius
22
ACS Paragon Plus Environment
Page 25 of 34
Journal of Agricultural and Food Chemistry
Table 2
Analytical characteristics of soil samples (0–10 cm depth)
Sample
Altitude
Slope
pH
Clay
m asl
(%)
(H2O)
g kg
Ca
Textural type (USDA)
-1
g kg
Lipidb
C/N -1
CECc
-1
g kg
Sd
cmolc kg
-1
-1
cmolc kg
R1
880
5
7.5
101
Sandy loam
170
16.5
1.45
47.1
33.9
R2
855
0
6.0
30
Loamy sand
67
24.7
0.59
12.2
7.6
R3
1150
15
5.9
158
Sandy loam
39
11.3
0.17
17.2
10.4
R4
1240
15
6.2
55
Sandy loam
65
29.6
0.54
19.0
13.3
C1
870
0
7.9
27
Silt loam
95
15.3
0.29
41.2
41.2
C2
870
0
8.4
228
Silt loam
17
12.7
0.10
21.2
21.2
C3
950
2
6.8
113
Sandy loam
92
13.3
0.49
21.8
20.8
C4
950
2
6.0
126
Sandy loam
49
13.1
0.18
12.6
10.2
B1
840
20
6.2
134
Loam
61
15.5
0.08
17.1
11.1
B2
825
5
6.2
45
Loamy sand
57
13.0
0.13
10.6
4.7
B3
1015
15
7.1
51
Sandy loam
44
14.2
0.13
15.0
7.5
B4
990
20
6.5
86
Sandy loam
88
17.1
0.37
25.0
16.8
F1
630
0
6.9
120
Silt loam
213
14.5
9.50
70.6
49.7
F2
624
8
8.7
125
Sandy loam
39
12.3
0.05
13.5
13.5
F3
1580
20
4.6
77
Sandy loam
69
14.7
0.42
23.5
2.1
1615
15
5.7
64
Sandy loam
64
8.2
0.14
41.2
4.0
0.1
30
5
1.4
1.8
5.11
F4 LSD a
e b
c
d
0.02 +
+
2+
2+
e
Total oxidizable soil C; Soxhlet extraction with petroleum ether; cation exchange capacity; sum of exchangeable bases (Na + K + Ca + Mg ); least significant difference based
on adjacent spatial replicates.
23
ACS Paragon Plus Environment
Journal of Agricultural and Food Chemistry
Page 26 of 34
Table 3 Relative concentrationa, ratios and statistical indices calculated from alkanes, FAs and alkanols in the soils Compound type, ratios
R1
R2
Total alkanes
27.0
R3
3.0
R4
22.2
3.1
B3
B4
68.1
49.6
45.9
22.4
4
283
197
509
202
C2
C3
C4
B1
47.3
33.8
38.8
21.9
Linear/branched
187
2n+1/2n
B2
C1
F1
9
2
4
4
7
6
12
4
6
8
9
8
> C20/C20/