Use of Thermally Assisted Hydrolysis and Methylation (THM-GC-MS) to Unravel Inﬂuence of Pottery Production and Post-Depositional Processes on the Molecular Composition of Organic Matter in Sherds from a Complex Coastal Settlement

: Ceramic fragments from the Islet of Guidoiro Areoso (NW Spain), covering a wide range of cultural periods (Neolithic to Late Bronze Age), have been studied by color analysis, elemental analysis of carbon (C) and nitrogen (N), and molecular analysis (thermally assisted hydrolysis and methylation, THM-GC-MS), in order to identify the organic matter (OM) in the prehistoric pottery and reveal information on ceramic production techniques, food remains and post-depositional effects. Results showed that the strong marine inﬂuence (sherds recovered from coastal deposits) and microbial activity (recovery from waste deposits, “cuncheiros”) had a profound effect on C/N ratio and molecular composition (N-rich protein and chitin structures). Other organic ingredients originated from the material used for creating the ware (detected as pyrogenic OM) and possibly food remains (fatty acid ﬁngerprints). Dark-colored ware was enriched in both pyrogenic OM from incomplete combustion and non-bacterial fatty acids. Fatty acid patterns could not be related to possible vessel use, and markers of aquatic resources were scarce, or absent. It is argued that THM-GC-MS of pottery fragments is useful for understanding how an archaeological deposit developed in time, what kinds of OM are present, and possibly to make a pre-selection of samples with high potential for more cost-demanding dietary molecular assessments. Beaker Imitation type wares [30] were recovered. This kind of pots are related to Late Neolithic chronologies. The Early Bronze Age is mainly represented by different types of Bell Beaker pottery (undecorated Bell Beakers and some of the regional variations of the Bell Beaker pottery for Galicia [31]). The heterogeneity of Early Bronze Age pottery suggests several moments of activity along this period. During the Middle Bronze Age the midden continues to grow in all directions, containing many in situ potsherds. Then, at the beginning of Late Bronze Age (Late Bronze Age I), ceremonial activities related to the midden halted, and the midden was buried by sediment. On this new soil surface, large ware had been placed on top. Finally, during the Late Bronze Age II, different structures were built on the mounds. Stones, cists, millstones and typical Late Bronze Age pottery were found.


Introduction
In spite of the advances in organic chemistry applications for archaeometry, there is still potential to be explored. Molecular research on organic constituents in ceramic materials have traditionally focused on total lipid extracts, aiming to characterize the preserved food residues, for instance by biomarker analysis (liquid and gas chromatography (GC) coupled to mass spectrometry (MS), such as GC-MS) [1][2][3][4] or compound-specific isotope ratio MS of lipid extracts from ceramic fragments [5,6], to obtain invaluable dietary information. However, recovered ceramic fragments contain organic matter (OM) from a wide range of possible sources, namely (1) the raw material clay fabric (incomplete burn-off during firing), (2) organic tempers and binders (incomplete burn-off), (3) soot deposited into or onto the ceramic during firing or during use (cooking), (4) decoration and finishing materials (varnishes, pigments, etc.), (5) use (food) residues, (6) post-use accumulations from the burial environment (soil OM, dissolved OM), (7) contamination from handling during excavation or laboratory treatments [7][8][9]. Identification and characterization of Beaker Imitation type wares [30] were recovered. This kind of pots are related to Late Neolithic chronologies. The Early Bronze Age is mainly represented by different types of Bell Beaker pottery (undecorated Bell Beakers and some of the regional variations of the Bell Beaker pottery for Galicia [31]). The heterogeneity of Early Bronze Age pottery suggests several moments of activity along this period. During the Middle Bronze Age the midden continues to grow in all directions, containing many in situ potsherds. Then, at the beginning of Late Bronze Age (Late Bronze Age I), ceremonial activities related to the midden halted, and the midden was buried by sediment. On this new soil surface, large ware had been placed on top. Finally, during the Late Bronze Age II, different structures were built on the mounds. Stones, cists, millstones and typical Late Bronze Age pottery were found.
Separations 2021, 8, x FOR PEER REVIEW a ceremonial role (including banquettes/feasts), from the final Late Neolithic stage t Middle Bronze Age. At the first moment of activity in the midden, Penha and Bell B Imitation type wares [30] were recovered. This kind of pots are related to Late Neo chronologies. The Early Bronze Age is mainly represented by different types of Bell B pottery (undecorated Bell Beakers and some of the regional variations of the Bell B pottery for Galicia [31]). The heterogeneity of Early Bronze Age pottery suggests se moments of activity along this period. During the Middle Bronze Age the midden co ues to grow in all directions, containing many in situ potsherds. Then, at the beginni Late Bronze Age (Late Bronze Age I), ceremonial activities related to the midden h and the midden was buried by sediment. On this new soil surface, large ware had placed on top. Finally, during the Late Bronze Age II, different structures were bu the mounds. Stones, cists, millstones and typical Late Bronze Age pottery were foun Figure 1. Location of the Guidoiro Areoso islet in the Ría de Arousa (modified after [32]).

Sample Selection and Preparation
More than 3600 fragments were recovered [33], 145 of which were studied in d of which a final selection of 50 samples was made for color, elemental and molecular yses ( Table 1). The selection aimed to obtain a representative number of fragments each cultural period. Therefore, we selected fragments that allowed for a detailed ch terization and well-defined position within the deposit stratigraphy.
The surfaces of the pottery fragments were carefully and gently cleaned from e nally bound sediment using distilled water, and the samples were ground to very powder (<50 µm).

Sample Selection and Preparation
More than 3600 fragments were recovered [33], 145 of which were studied in detail, of which a final selection of 50 samples was made for color, elemental and molecular analyses ( Table 1). The selection aimed to obtain a representative number of fragments from each cultural period. Therefore, we selected fragments that allowed for a detailed characterization and well-defined position within the deposit stratigraphy.
The surfaces of the pottery fragments were carefully and gently cleaned from externally bound sediment using distilled water, and the samples were ground to very fine powder (<50 µm).

Color and Elemental Analyses
Quantitative color analysis in the CIE L*a*b*C*h (CIELAB) color space was done using a Konica Minolta CR-5 colorimeter [34]. Samples were analyzed in the powder form. L* represents luminosity (L), primary color components (a*: green-red axis, and b*: blue-yellow axis) chromaticity (C*) and hue (h). Carbon and nitrogen contents were obtained by combustion of 100 mg of finely milled samples using a Truspec CHNS analyzer (Leco) hosted at the RIAIDT services of the University of Santiago de Compostela.

THM-GC-MS
An aliquot of 10 µL of tetramethylammonium hydroxide (TMAH, in water, 25%, Sigma-Aldrich, St. Louis, MO, USA) was added to the fire-polished quartz tubes in which 1-2 mg of sample was embedded. The setpoint temperature of the THM reaction was 650 • C, maintained for 20 s, using a Pyroprobe pyrolyzer (CDS Analytical, Oxford, PA, USA). The pyrolyzer was connected online to an 8860 GC (Agilent, Santa Clara, USA). Analytical-grade helium was the carrier gas (1 mL/min constant flow). The pyrolyzer's interface oven, the GC inlet and the GC-MS interface were held isothermally at 325 • C. The GC was equipped with a HP-5MS non-polar column, and was in split mode (1:25). The GC was heated from 60 to 325 • C at a rate of 20 • C/min. The initial and final isothermal periods were 5 and 2.5 min, respectively (total runtime 19.75 min). The MS operated in electron ionization (70 eV) mode, scanning in the range of 50 to 500 amu, using a solvent delay period of 5 min to allow the residues of the TMAH and aqueous solvent to elute before the MS was activated. Relative abundances between the THM products (identified by literature and NIST libraries) were based on peak areas of their dominant m/z fragments (Appendix A) and expressed as percentage of total quantified peak area (% TQPA). For background on THM-GC-MS, readers may refer to Challinor [24], and for recent applications in diverse fields, to He et al. [35].

Data Evaluation
We applied principal component analysis (PCA) on the relative proportions data from THM-GC-MS, after applying a centred logratio transformation (clr) to avoid the close data nature of compositional data. We used one-way analysis of variance (ANOVA) to explore differences between samples of different periods and different depositional features (alleged influence of external factors such as marine influence). All statistical tests were performed in SPSS 20.0 software (IBM, New York, NY, USA).

Color and Elemental Analysis
The ceramic fragments analyzed generally had low luminosity (L*, n = 50; 44.1 ± 7.6; Table 1) and chromaticity was low as well (C*, 17.6 ± 7.4). The color components a* and b* had positive values, indicating that the basal color is a mixture of red and yellow, with predominance of yellow over red as the hue exceeds 45 • in all measured samples.
The elemental analysis showed that carbon (C) contents ranged from 0.65 to 4.99% (1.68 ± 0.87%; Table 1), while nitrogen (N) content ranged from 0.05 to 0.21% (0.10 ± 0.03%). The C/N molar ratio ranged from 11.2 to 52.9 (19.9 ± 7.8). Figure 2 shows some relationships between luminosity and elemental composition (%C and C/N). Clearly, dark-colored (i.e., low luminosity value) samples have higher C content and higher C/N than light-colored ceramic fragments. Hence, OM is the main component that decreases luminosity. The negative correlation between L* and C/N indicates that the darkening agent is enriched in C relative to N.

Thermally Assisted Hydrolysis and Methylation (THM-GC-MS)
The THM-GC-MS chromatograms ( Figure 3; Supplementary Material) allowed to identify and semi-quantify 104 compounds (Appendix A). Most of the THM products have a molecular structure based on polymethylene chains (32.5 ± 17.9% of TQPA; Figure  4). The dominant subgroup of these compounds are fatty acid methyl esters (FAMEs), accompanied with fatty diacids (DAMEs), linear alkanes and isoprenoid alkanes. The carbon number of the FAMEs ranges between 7 and 18. These compounds can originate from oils and fats (food remains) but also from fatty acid moieties from plant remains and microbial sources (plankton, fungi, bacteria, etc.). The presence of branched (iso/anteiso) C15 FAME probably indicates that at least some of the FAMEs originate from bacterial sources of OM. We identified traces of only one FAME with mid-chain methoxy substitution, probably from cutin (plant cuticles) in selected samples. We searched for some of the markers reported by Shoda et al. [25] and identified a minor trace of a compound that was tentatively ascribed to an isoprenoid pristane-based FAME with m/z 88 and 312 at the expected retention time of pristanic acid methyl ester (ME), and possibly m/z 101 and 326 of phytanic acid ME (too small for peak area calculation). It is concluded that the abundance (relative to total OM that was accounted for in the THM-GC-MS data evaluation) in the studied sherds is much lower than that of the samples studied by these authors [25], or not present at all, and biomarker evaluation for palaeodiet in the present data is not sustained.

Thermally Assisted Hydrolysis and Methylation (THM-GC-MS)
The THM-GC-MS chromatograms ( Figure 3; Supplementary Material) allowed to identify and semi-quantify 104 compounds (Appendix A). Most of the THM products have a molecular structure based on polymethylene chains (32.5 ± 17.9% of TQPA; Figure 4). The dominant subgroup of these compounds are fatty acid methyl esters (FAMEs), accompanied with fatty diacids (DAMEs), linear alkanes and isoprenoid alkanes. The carbon number of the FAMEs ranges between 7 and 18. These compounds can originate from oils and fats (food remains) but also from fatty acid moieties from plant remains and microbial sources (plankton, fungi, bacteria, etc.). The presence of branched (iso/anteiso) C 15 FAME probably indicates that at least some of the FAMEs originate from bacterial sources of OM. We identified traces of only one FAME with mid-chain methoxy substitution, probably from cutin (plant cuticles) in selected samples. We searched for some of the markers reported by Shoda et al. [25] and identified a minor trace of a compound that was tentatively ascribed to an isoprenoid pristane-based FAME with m/z 88 and 312 at the expected retention time of pristanic acid methyl ester (ME), and possibly m/z 101 and 326 of phytanic acid ME (too small for peak area calculation). It is concluded that the abundance (relative to total OM that was accounted for in the THM-GC-MS data evaluation) in the studied sherds is much lower than that of the samples studied by these authors [25], or not present at all, and biomarker evaluation for palaeodiet in the present data is not sustained.

Thermally Assisted Hydrolysis and Methylation (THM-GC-MS)
The THM-GC-MS chromatograms ( Figure 3; Supplementary Material) allowed to identify and semi-quantify 104 compounds (Appendix A). Most of the THM products have a molecular structure based on polymethylene chains (32.5 ± 17.9% of TQPA; Figure  4). The dominant subgroup of these compounds are fatty acid methyl esters (FAMEs), accompanied with fatty diacids (DAMEs), linear alkanes and isoprenoid alkanes. The carbon number of the FAMEs ranges between 7 and 18. These compounds can originate from oils and fats (food remains) but also from fatty acid moieties from plant remains and microbial sources (plankton, fungi, bacteria, etc.). The presence of branched (iso/anteiso) C15 FAME probably indicates that at least some of the FAMEs originate from bacterial sources of OM. We identified traces of only one FAME with mid-chain methoxy substitution, probably from cutin (plant cuticles) in selected samples. We searched for some of the markers reported by Shoda et al. [25] and identified a minor trace of a compound that was tentatively ascribed to an isoprenoid pristane-based FAME with m/z 88 and 312 at the expected retention time of pristanic acid methyl ester (ME), and possibly m/z 101 and 326 of phytanic acid ME (too small for peak area calculation). It is concluded that the abundance (relative to total OM that was accounted for in the THM-GC-MS data evaluation) in the studied sherds is much lower than that of the samples studied by these authors [25], or not present at all, and biomarker evaluation for palaeodiet in the present data is not sustained.     Many compounds with at least one atom of nitrogen (N-compounds) were identified (18.2 ± 11.2%; Appendix A), further indicating the contribution of microbial OM to the molecular fingerprints. These compounds include alkylpyrroles, alkylindoles, anilines and succinic acid ME. In addition to microbial N, these compounds (except for succinic acid ME) can also originate from plant protein, but as discussed below, on the basis of statistical data evaluation, this is an unlikely source of N. Finally, the N-compound benzonitrile was also identified in the chromatograms of most samples. The most probable source of benzonitrile is pyrogenic OM (PyOM) [20]. Many compounds with at least one atom of nitrogen (N-compounds) were identified (18.2 ± 11.2%; Appendix A), further indicating the contribution of microbial OM to the molecular fingerprints. These compounds include alkylpyrroles, alkylindoles, anilines and succinic acid ME. In addition to microbial N, these compounds (except for succinic acid ME) can also originate from plant protein, but as discussed below, on the basis of statistical data evaluation, this is an unlikely source of N. Finally, the N-compound benzonitrile was also identified in the chromatograms of most samples. The most probable source of benzonitrile is pyrogenic OM (PyOM) [20].
This PyOM is also the likely source of polycyclic aromatic hydrocarbons (PAHs; 20.7 ± 9.8%), which are more abundant than is usual for THM fingerprints of biomass. In all likelihood, this can be explained by the substantial portion of charred remains in the OM in the sherds, formed during firing (incomplete burn-off) or cooking (soot, charred food). The first source of PyOM signals is usually dominant for ceramics produced (partially) at low temperature or with truncated firing (reflected by a dark-colored core in contrast to light-colored surfaces, also referred to as "sandwich" morphologies, with a heat-isolated dark-colored nucleus; [8,36,37]).
Monocyclic aromatic products include penta-and hexamethylbenzene (14.5 ± 11.4%), which may originate from marine sources [38] and which had not been observed in a study of fragments from non-coastal settlements in NW Spain [20]. This may indicate that some of N-compounds also originate from marine sources such as plankton, as will be discussed later.
Finally, we identified several compounds containing sulfur (pronounced in two samples from the Middle Bronze Age) or phosphorus (most abundant in two samples from the Neolithic period; Figure 4), other compounds from contamination such as tert-butylphenol derivatives (plastic additive, surfactant), and compounds which may be formed during the THM reactions as secondary rearrangements such as branched alkanes. We are aware that the THM reaction may generate artifacts, which can be catalyzed by the presence of Fe oxides [39].
In summary, the THM-GC-MS analyses indicate the presence of OM that has a prefiring origin (temper-derived PyOM, perhaps some of the FAMEs), vessel use origin (soot-derived PyOM, perhaps some food-derived FAMEs) and from after burial (microbial OM from soil and marine fluxes, lignin and FAMEs).
The PCA ( Figure 5) provided 3 PCs that explain 60% of total variance. PC1 has elevated negative loadings for the short-chain FAMEs including bacterial FAMEs and positive ones for polyalkyl aromatic compounds. PC1 is not correlated to other variables and is not of major interest for ceramic use identification, as it indicates the balance between different post-depositional inputs and, possibly, also analytical artifacts.
PC2 has high positive loadings for N-containing and lignin products (G6, S6) and negative ones for C 18 FAME. A possible explanation is that PC2 marks the input of postdepositional "exogenous" inputs (from the deposit), such as microbial and plant-derived soil OM constituents.
Regarding differences between samples, and based on the argued main likely sources of the THM products (Table 2), samples with the highest relative proportion of PyOM are CAM4A-037, -040, -080, -085 and -100, corresponding to various periods. We believe that the most likely source of the PyOM is the incomplete combustion of soil OM from collected raw material and organic binder/temper, either way "pre-firing". If so, this means that these samples were subjected to a firing regime that facilitated preservation of OM, e.g., due to low temperature or short duration; or other parameters such as use of thick material and/or relatively high loads of temper. Dealing with relative proportions implies that it can also be simply that these fragments are relatively "clean" (few food or post-depositional organic inclusions). A source in surface-bound soot is unlikely due to associations between PC3 and color analysis (indicating bulk OM contributing to low L*, not effects of a thin layer of surface-adhered soot).  Part of the variation can indicate the balance between fatty acids and compounds of microbial origin (FAMEs, poly-alkyl MAH and PAH, PC1). If part of the non-microbial FAMEs detected indeed originate from food remains, this would imply a great value for selecting samples for future analysis. In this regard, samples CAM4A-018 and -081 have signals that may be the most likely associated with food remains. Samples CAM4A-007, -  Part of the variation can indicate the balance between fatty acids and compounds of microbial origin (FAMEs, poly-alkyl MAH and PAH, PC1). If part of the non-microbial FAMEs detected indeed originate from food remains, this would imply a great value for selecting samples for future analysis. In this regard, samples CAM4A-018 and -081 have signals that may be the most likely associated with food remains. Samples CAM4A-007, -082, -097, -101 and -102 would have relatively strong marine influence. For fragment CAM4A-101 this coincides with a sherd that is colonized by mollusks.
Samples with the highest proportions of lignin and microbial N are CAM4A-014 and -081. The most likely explanation is a relative strong signal from terrestrial post-depositional inclusions (fungi and plant remains). But an effect of decay intensity cannot be discarded as chitin and lignin are probably among the labile sources of OM in the present environment. Bacterial FAMEs are most abundant in CAM4A-010 and -081, and least abundant in -097 and -101. Finally, sample CAM4A-082 has a relatively strong signal of leaf cuticle-derived OM (cutin).

Comparing THM Data with Color and Elemental Composition
The strongest correlation between the different datasets is the negative link between luminosity (L*) and the PyOM-related PAHs such as naphthalene, fluorene, biphenyls and phenanthrene/anthracene (Figure 2), benzonitrile and PC3. This link between darkness and PyOM is a strong indication of incomplete burn-off as the cause. It was also the main link between color and Py-GC-MS analyses of pottery fragments from non-coastal settlements in NW Spain [20]. It is remarkable that the THM-GC-MS data also highlight such an effect as the main products of PyOM, benzene and toluene, are not detected by THM-GC-MS. The total sum of PAH is not correlated to L* probably because of contributions from non-PyOM, such as plankton, to the poly-alkyl PAHs (mainly polyalkylnaphthalenes and -indenes). Indeed, trimethylindene and hexamethylbenzenes are correlated positively with L*. There is no apparent effect of cultural period on L* or molecular proxies of PyOM (low values for L* and high values for PyOM for samples CAM4A-009, -037, -040, -080, -085 and -105). These data reflect low firing temperature and/or ceramic thickness (preservation of binder/temper in inner sections), but we have not elaborated this further as Py-GC-MS is more suitable for these purposes and because of the likely catalytic effects of Fe oxides [18,40].
The negative correlations between L*, %C, C/N and PyOM indicates that the PyOM has an effect on L* that the OM from N-rich marine and terrestrial microbial sources does not have. The relationships between other color parameters and the molecular fingerprints are relatively weak and may be due to spurious co-variation effects. Suffice it to mention that the PyOM tends to mitigate the expression of a* (redness) component, probably related to Fe oxides, to a larger extent than b*, C* or h. This negative effect of PyOM on redness may also be related to the association of PyOM with low temperatures (caused by incomplete burn-off) that would also reduce the efficiency of Fe (oxy)hydroxide dehydration to strongly pigmenting red oxides such as hematite. Note that this negative effect makes it unlikely that Fe oxide-catalyzed secondary reactions are the main formation mechanism of PyOM products, lending support to our approach and interpretation of PC3. The %C is negatively correlated with polyalkyl MAH and PAH, probably due to the fact that the THM-GC-MS data is a closed dataset: molecular fingerprints of samples with low OM content will be more strongly affected by marine-derived OM than samples with high OM content. Indeed, we did not find consistent trends between the intensity of the marine signal from THM and position in the deposits or age/cultural period. If this argument is correct, the variations in the fingerprints indicate that not only the PyOM but also the mentioned FAMEs are mainly pre-depositional. Then, the marine, and probably also terrestrial signals (microbial N-rich moieties and plant-derived lignin and perhaps cutin) developed mainly post-depositionally.
The correlations between %N and THM products are different to those for %C. For %N, positive correlations (p < 0.001) were observed for iso-C 15 FAME (bacterial lipids), C 9 -diacid (degraded C 18 fatty acids), N compounds (alkylindoles, N-methylsuccinimide) and the phenolic products P6 y G6. This highlights the influence of soil microbes (bacteria and fungi) on OM composition. The association between microbial and lignin products further supports a post-depositional source of lignin, instead of incomplete burn-off of organic temper or soil OM that was present in the raw material before firing. Perhaps the inclusion of plant signals is stimulated by microbial activity on sherd surfaces, or because enhanced microbial alteration provokes stronger interactions between the ceramic fabric and the soil environment in general. Relationships between THM fingerprint and C/N are consistent, with correlations between C/N with FAMEs and PyOM markers.
The ANOVA analysis using cultural period (Neolithic, Early Bronze Age, Middle Bronze Age, Late Bronze Age I and Late Bronze Age II; 7-15 samples/period) revealed few significant differences in molecular composition of the OM. The tentatively identified isoprenoid FAME (pristanic acid ME) is enriched in samples from Late Bronze Age I, which could reflect an effect of different vessel use between periods (marine oils), but in the light of planktonic aquatic sources of OM and the uncertainty in identification of this compound, this observation cannot be considered more than a reason to pursue further analyses using alternative methods of the samples from the Late Bronze Age I in particular. Furthermore, an aliphatic ether is enriched in samples from the Early Bronze Age, and phosphoric acid ME in samples from the Neolithic. Color data are more strongly differentiated on the basis of period, with low values of a*, b* and C* in samples from the Late Bronze Age II and Neolithic, and high values for the other samples, in particular Middle Bronze Age and Late Bronze Age I. For hue, the tendency is reverse. This probably indicates differences in raw material selection and firing conditions, but because these results are not related to OM chemistry, they were not explored further.
An additional ANOVA aimed to identify differences between recovery locations. The midden was expected to coincide with samples with higher likelihood of post-depositional contamination (VAR1; high-intermediate-low) and contamination type (VAR2; waste deposit-post waste deposit-marine influence) than samples from the superimposed layers (Table 1). It appeared that samples with high general contamination likeliness (VAR1) were enriched in N-containing THM products (alkylindoles and -pyrroles), confirming a stronger microbial influence in the midden (including possibly mollusk-derived OM), whereas C 16 FAME is enriched in samples with low contamination. There is no color effect on VAR1. For VAR2 (contamination type), samples from the waste deposit have positive PC1 scores from THM (low fatty acid proportions and high microbial influence) whereas the other two categories have negative PC1. Indeed, C 16 and C 18 FAME are depleted in the samples from the waste deposit. These results might indicate that samples that are prolific in FAMEs are good candidates for dietary analyses of the fatty acid profiles using compound-specific IRMS of lipid extracts.

Conclusions
The THM-GC-MS analyses have shown that the OM present in the ceramic fragments from Guidoiro Areoso Islet originated mainly from: (1) marine and terrestrial microorganisms (protein, chitin, bacterial fatty acids), (2) charred residues (PyOM, formed during firing, from soil OM or temper/binder), (3) uncharred plant remains (traces of lignin and cutin from post-depositional inclusion of plant-derived soil OM) and (4) fatty acids of unknown source (possibly from food remains). It is clear that the burial environment, with high loads of organic waste and marine influence, caused a relatively complex mixture of OM sources that hampered the unambiguous identification of food remains, let alone identifying precise ingredients. Nevertheless, the study revealed some interesting links between color, elemental and molecular composition of ancient ceramic fragments that allowed for a better understanding of the changes in deposition conditions and will be used to select a subset of samples of which relatively "clean" fatty acid patterns are expected, in particular dark-colored ones, with high C/N and with high non-bacterial FAME proportions. We believe that light-colored fragments, with low fatty acid contents and low C/N, would reveal fatty acid profiles that are too severely affected by microbial inputs. Shinya Shoda (University of York, UK) for sharing THM-GC-MS chromatograms of Neolithic pottery, which helped to discard the biomarker potential of our obtained data. We acknowledge the valuable comments from two anonymous reviewers.

Conflicts of Interest:
The authors declare no conflict of interest.