New insights into the use of Neolithic pottery in Guangxi of South China: organic residue analysis of experimental and archaeological pottery

As one of the world’s earliest pottery origins in Eurasia, the role of pottery during Neolithic period in South China remains unclear. Here, we present a new understanding on the utilization pattern of early pottery from the Neolithic sites of Dayan, Dingsishan and Zengpiyan in Guangxi province, China. Compound-specific stable isotope analysis of pottery and bone lipids show that a persistent use of terrestrial C3 and C4 resources throughout the pottery utilization history highlight the culinary practice and terrestrial exploitation during Neolithic period. Contrary to the traditional view that the primary use of these pottery was closely linked with aquatic resource processing, no aquatic biomarkers were found in the potsherds, yet these compounds were detected in the cooking experiments in this study, indicating that aquatic resources in these sites were probably processed or consumed in other ways. Moreover, a remarkable diachronic change from extensive use of C3 and C4 resources to the intensive use of certain kinds of C3 resources in different Neolithic stages, further indicating that an ecological tradition rooted in the subtropical and tropical environments as well as the Neolithic social evolution exerted great influence upon the initial production and the subsequent utilization of early pottery in prehistoric South China.


Introduction
The phenomenon of pottery emergence towards the end of the Late Pleistocene made East Eurasia the recognized area of pottery origin.A series of early pottery vessel traditions that appeared simultaneously in South China, the Japanese Islands and the Russian Far East between approximately 20000-10000 BP [1][2][3][4], launched the ongoing discussion on the pottery use and the trajectories of prehistoric social evolution associated with them [5,6].In contrast to the traditional view that pottery use was linked with agricultural communities [7][8][9], new advances in archaeological research have been successful in linking pottery use with the adaptation to changing ecological and social structures even before farming within Eurasia [10,11] and demonstrating that the use of early pottery and resource processing had a strong association with the newly established ecological conditions in the Japanese archipelago and the Russian Far East during the transition to Holocene [12,13].
As one of the oldest regions of pottery origin in Eurasia, the invention and utilization of pottery in South China have always attracted much attention [1,5,10,11].The earliest pottery appeared in several sites in South China at the end of the terminal Pleistocene, such as pottery containers excavated from Yuchanyan cave in Hunan province, Xianrendong cave in Jiangxi province and Dayan and Zengpiyan caves in Guangxi province [1,4,14], and gradually adopted in the subsequent millennia with further social evolution.It has been suggested that the initial function of these early pottery had more connection with resource processing, particularly for aquatic resources [14][15][16][17].Significantly, the most prominent feature concerning the pottery-bearing sites across the area is the apparent increase in the intensity of shellfish, which was not seen in quantities before the early Neolithic periods [18,19].Accordingly, extensive aquatic resource exploitation during the local Neolithic has been demonstrated through stable carbon and nitrogen isotope analysis of human and animal bone collagen [20,21].
However, recent molecular analysis of lipid residues in early pottery from Neolithic shell mounds in South China did not detect reliable evidence of aquatic biomarkers, such as isoprenoid fatty acids (e.g., 4,8,12-TMTD, phytanic and pristanic acids) and ⍵-(o-alkylphenyl) alkanoic acids (APAAs), and their presence is indicative of cooking aquatic foods [22].Conversely, a series of terrestrial biomarkers, such as phytosterols, terpenes, alkanes and animal cholesterols, predominate in these early potsherds [23,24], suggesting that the vessels were primarily used for processing terrestrial plants and animals.Indeed, it has long been argued that the use of early pottery vessels was affected by the ecological and social traditions across different regions [5,10,25,26].This is particularly true in South China, where abundant subtropical and tropical resources remain relatively stable through millions of years [27].Similarly, the archaeological records of lithic assemblages and quantities of the terrestrial remains coexisted with the pottery sherds in the study sites may suggest a deep relationship between pottery use and the adaption to the ecological and social traditions in South China [5,10].Although archaeological evidence indicates that resource processing was a crucial factor for pottery use in prehistoric South China, detailed pottery residue analysis remains limited.
This study aims to investigate the function of early pottery and diachronic change of resource utilization patterns throughout the Neolithic times in South China by analyzing lipid residues in pottery and contemporaneous animal bones from the Neolithic sites (ca.12,000-6000 cal BP) of Dayan, Dingsishan and Zengpiyan in Guangxi province, China.While it has long been assumed that these early pottery vessels were closely linked with the processing of aquatic resources, little evidence of aquatic exploitation was found in previous studies [23,24].Cooking experiments were also conducted in this study to determine the formation of the aquatic biomarkers through different cooking conditions and further improve our understanding of the characteristic of aquatic lipids in pottery vessels after cooking.

Principles of pottery residue analysis
Organic residues in archaeology include DNA, proteins, lipids and other biomolecular from archaeological deposits.Of which, lipids are defined by their small atomic mass units (> 200 but < 1000 amu) and hydrophobicity, this property forms the basis for their preservation and extraction from archaeological materials [28].Biomarkers, as "chemical fingerprints" of lipids, are the compounds that originated from specific biological and biochemical components of natural products [29].These compounds are usually preserved in archaeological pottery, via surface adsorption or trapping within molecularsized pores of the inside wall.Based on the concept of archaeological biomarkers, lipid residues extracted from pottery could provide molecular evidence for the foodstuffs (such as aquatic resources, terrestrial animals, edible plants), and various products (such as resins, waxes, tars), which were likely to have been exploited in the past [29,30].On the other hand, plants have distinct δ 13 C values depending on their photosynthetic pathways, such as C 3 , C 4 or CAM plants, the δ 13 C values of animal and plant tissues generally reflect the isotopic level of local environments [31].Moreover, the δ 13 C values of food and fluid ingested by animals also have a strong impact on the isotopic compositions of the tissues they synthesize [32].Thus, the application of compound-specific stable isotope analysis of the fatty acids that composed of lipid could contribute to further distinguishing the animal and plant type, such as ruminants, non-ruminants and dairy fats [30].To date, the use of organic residue analysis has shown great potential in revealing the diversified resource utilization in ancient pottery among different regions, such as the earliest use of aquatic resources in the Japanese archipelago and North Africa [33,34]; terrestrial plant processing in China, Southeast Asia and North America [35][36][37][38]; and animal domestication and beeswax utilization in Europe [39][40][41].

Archaeological background
The archaeological sites of Dayan, Zengpiyan and Dingsishan are located in Guangxi province of China with Dayan and Zengpiyan sites lying in the karst plain along the Li Jiang River in Northern Guangxi, while the Dingsishan site is situated on the banks of Yong Jiang River in Southern Guangxi (Fig. 1).These sites were Archaeological deposits in these sites can be divided into four major phases: terminal Paleolithic (15000-12000 BP), Early Neolithic (12000-9000 BP), Middle Neolithic (9000-7000 BP) and Late Neolithic (7000-6000 BP), through radiocarbon dating and material cultural associations [42,43].The chronology and cultural frameworks based on archaeological assemblages in these sites suggested a typical local culture technocomplex characterized by flexed burials, flake lithic assemblages and other stone, bone and shell tools, as well as an abundance of aquatic and terrestrial animal and plant remains [14,19,45].Notably, early pottery vessels (Fig. 2) were found in these sites, with the earliest pottery container (mainly round-bottomed 'fu') dated to approximately 12000 BP at Dayan and 10000 BP at Dingsishan, and are comparable to contemporary pottery from the sites of Zengpiyan cave in Guangxi province and Yuchanyan cave in Hunan province, China [42].
There are several shared features between the early pottery (dated to approximately 12000-6000 BP) in these sites, such as small quantity but large size, simple cord-marked and thick walls (about 2.9 cm thick for earliest potsherd from Dayan cave, Fig. 2c, d).These pottery vessels are tempered with large particles of calcite or quartzite, indicating low firing temperature and pristine technology in Early Neolithic [43].From the middle Neolithic, an increase in the quantity, quality and pottery types, particularly, a significant shift towards smaller pottery sizes and thin walls, suggests subsequent development in pottery technology [42,43].Besides the manufacturing technology, the functional need for resource preparation is also determined by the reconstruction of life histories of pot vessels through absorbed lipids [23,24], but the role of these early pottery within a changing social context in prehistoric South China remains unclear.

Materials
A total of 62 archaeological samples, including 52 potsherds (examples of the potsherds are represented in Fig. 2) and 10 animal bones which include Water buffalo (n = 3), Deer (n = 2), Small carnivore (n = 1), Wild boar (n = 1) and Freshwater fish (n = 3), from Neolithic sites of Dayan and Dingsishan and Zengpiyan in Guangxi, South China; combined with 36 modern replica potsherds and foodstuffs from cooking experiments (Additional file 1: Figure S1) are selected for lipid residue analysis through Gas Chromatography-Mass Spectrometry (GC-MS) and Gas Chromatography-Combustion-Isotope Ratio Mass Spectrometry (GC-C-IRMS).

Cooking experiments
Modern replica potsherds (n = 30) and food products (n = 6) in cooking experiments were used to generate comparative data for further understand the characteristics of aquatic lipids in pottery through different time and temperature.
For the cooking experiments, modern replica pottery vessels used in this study were wheel-thrown and manufactured in Jingdezhen, Jiangxi Province, China.The mineral composition of these pots mainly includes local Pot Clay and sand, and the firing temperatures were relatively high (around 1000 ℃) to avoid fire-crack.The food products (Additional file 1: Table S1), including shellfish (Viviparus), freshwater fish (Cyprinus carpio), meat (lamb), taro, Chinese yam, and leafy vegetables, were placed into these replica pots over the charcoal fire stoves.
To determine the temperature and time required for the formation of aquatic lipids in pottery vessels, particularly for aquatic biomarkers (such as, ω-(o-alkylphenyl)alkanoic acids (APAAs), phytanic, pristanic acids and 4,8,12-TMTD), series of temperatures on 230 °C, 250 °C and 280 °C adjustment from previous studies [46,47], were set as 5 h and 10 h (1 × 10 times) for three groups respectively, each group used the same portions of food ingredients (Additional file 1: Table S1).In addition, a small amount of roasting and raw meat from shellfish were also prepared for comparison (Additional file 1: Table S1).The experiments were conducted at the field workstation of the Institute of Archaeology, Chinese Academy of Social Science.The food products used in this study were bought from supermarket in Beijing and Guangxi province, China.

Lipid extraction
1-2 mm of the inside surface of the potsherd were first mechanically removed using a hand-drill in order to eliminate exogenous residue and possible contamination.About 1 g of potsherd was then ground to powder and placed into a glass tube.4 mL methanol was added and ultrasonicated for 15 min.800 μL of concentrated sulfuric acid was added and heated for 4 h at 70 ℃, after centrifuged for 5 min, the supernatant was further extracted with hexane (4 × 2 × 2 mL) and collected via a prepared Pasteur pipette three times.Finally, the extracts were dried under a gentle stream of nitrogen.An internal standard of 10 μL n-tetratriacontane (C 34:0 ) or n-hexatriacontane (C 36:0 ) was added to each sample before analysis by GC-MS.
The archaeological animal bones were used as isotopic baseline for pottery lipids in this study.These bone samples were taken from the same one where bone collagen and hydroxyapatite were well preserved [21,48].While given the preservation of the bone lipids, a twostep extraction [49] was conducted to obtain the endogenous lipid in bones.About 1 g of bone fragments (limb) were mechanically cleaned and ground to powder using a pestle and mortar.The bone powder was first washed with solvent extracted with dichloromethane/methanol (2:1 vol/vol, 4 × 2 × 2 mL); after being ultrasonicated for 15 min and centrifuged for 5 min, the supernatant was removed and repeated three times and then dried completely under a gentle stream of nitrogen.The second step was using acidified methanol procedure to the remaining bone powder exactly as above.

Analytical methods
All the lipid extracts were then analysed by GC-MS and GC-C-IRMS.Since the lipid composition of the archaeological potsherds (n = 52) has already been investigated in our previous work [23], GC-MS analysis were conducted on the archaeological bones (n = 10) and modern materials (n = 36) from experimental cooking.
GC-MS analysis of the methylated extracts was carried out using an Agilent 7890A gas chromatograph attached to a 5975C mass detector (Agilent Technologies) in electron ionization mode (70 eV) at the Institute of Vertebrate Paleontology and Paleoanthropology, Chinese Academy of Sciences.A capillary column DB-5HT (15 m × 0.32 mm × 0.1 μm) was used to separate samples. 1 μL of the sample was injected in a split-less mode.The temperature program was as follows: initial temperature 50 °C for 1 min, increased to 100 °C at 15 °C/min, then up to 375 °C at 10 °C/min, with 10 min isothermal.Helium was used as the carrier gas in a constant flow rate of 2.0 mL/min.The MS ion source temperature was 230 °C; the interface temperature was 300 °C; the MS quadrupole temperature was 150 °C.The spectrometer is run in fullscan mode and in selected ion monitoring (SIM) mode in order to monitor aquatic biomarkers.The group of ions at m/z 262, 290, 318, 346 were selected to assess APAAs of carbon chain lengths C 16 -C 22 , and the fragment ion of the base peak m/z 105; ions of m/z 74, 87, 213, 270 corresponding 4,8,12-trimethyltridecanoic acid (TMTD) fragmentation; m/z 74, 88, 101, 312 corresponding to pristanic acid; m/z 74, 101, 171, 326 corresponding to phytanic acid, respectively [22,50].Identification of common lipids was performed by comparison with the NIST mass spectral database or mass spectra from the published literature.
Compound-specific carbon stable isotope analysis of archaeological samples was carried out on a Thermo Trace Gas Chromatography Ultra coupled with an Isotope Ratio Mass Spectrometer (Thermo Fisher Scientific MAT253) via a combustion reactor held at 1000 °C, at the Institute of Atmospheric Physics, Chinese Academy of Sciences.All samples were diluted with hexane and 1 μL of each sample was injected into a DB-5 MS fused silica capillary column (30 m × 0.25 mm, 0.25 μm film thickness).The GC oven temperature was set for 50 °C (2 min) to 120 °C at 15 °C/min, then from 120 to 300 °C at 5 °C/ min, with a final isotherm at 300 °C for 16 min.Highpurity grade helium was used as a carrier gas maintained at a constant flow rate of 2.0 mL/min.Compounds were combusted to CO 2 and ionised in a mass spectrometer.The ion intensities of m/z 44, 45 and 46 were monitored so as to automatically compute the 13 C/ 12 C ratio for each sample peak by comparison with a calibrated CO 2 reference gas of known isotopic composition through Isodat 3.0 Gas Isotope Ratio MS Software (version 3.0; Thermo Fisher).The analysis results are reported as per mil (‰) relative to internationally defined standards for carbon (Vienna Pee Dee Belemnite, VPDB).Isotopic values of each sample were monitored by use of an International standard mixture of fatty acid methyl esters with known δ 13 C values (Indiana F8-3 standard; obtained from Biogeochemical Laboratories, Indiana University), the measurement error of δ 13 C values of each compound based on duplicate analysis was less than ± 0.4‰.Moreover, the δ 13 C values were corrected for the potential carbon atoms that occurred during methylation using the measured carbon isotopic value of the methanol used in pretreatments.
Notably, our results show that while the heating conditions at 250 ℃ for 5 to 10 h can produce partial APAAs with carbon length C 18 and C 20 , abundant long-chain APAAs (carbon length of C 16 -C 22 ) are generated at a higher temperature of 280℃ for 5 h (Additional file 1: Table S1).Moreover, contrary to a previous study [47], APAAs were not detected at a lower temperature of 230 ℃ for either 5 or 10 h (Fig. 3 and Additional file 1: Table S1), suggesting an increased temperature at ≥ 250 ℃ should be the prerequisite for the formation of APAAs.In addition, the sealed condition is seemingly not necessary for the production of APAAs, as abundant APAAs were detected in fired meat of shellfish under direct open fire at 280 ℃ for only 30 min, which further suggests that high temperature is needed to form these compounds (Additional file 1: Table S1).
APAAs with different carbon lengths and low concentrations of phytanic acid were detected either in the separate vessels containing aquatic and plant foods or in a vessel containing both of them (Additional file 1: Table S1).While 4,8,12-TMTD were only detected in the vessels containing aquatic foods (Fig. 4 and Additional file 1: Table S1), which aligns with previous studies that such a compound occur stably in aquatic environments owing to the fat catabolism of aquatic organisms [22,51].Moreover, as 4,8,12-TMTD was also detected in the raw shellfish as well as vessels containing aquatic foods cooked at a lower temperature of 230 ℃ for 5 h (Additional file 1: Table S1) in this study, it seems to suggest that the presence of this compound may not require high temperature or protracted heating, but rather reflects a typical characteristic of aquatic lipids.
Thus, based on the cooking experiments in this study, we suggest that despite the formation of APAAs in pottery needs high heating conditions, 4,8,12-TMTD may not.Instead, it suggests that if the aquatic foods were ever processed in pottery even at a lower heating condition, certain amount of 4,8,12-TMTD could be detected.However, the post-depositional survival of this compound under different environmental conditions needs further investigation.

Lipid molecular signature in archaeological bones
The preservation of bone lipids is an important prerequisite for lipid residue analysis.For the archaeological bones in this study, the total lipid concentrations decreased from 1128.8 to 750.5 μg g −1 (Additional file 1: Table S2) by using the modified procedure of twostep extraction [49], probably due to the loss of both endogenous lipids and exogenous contaminants, such as phthalates from plasticizers, during the repeat of solvent washes in the first step.An obvious removal of several phthalates in deer bone confirmed the efficiency of this method (Additional file 1: Figure S2).Furthermore, the lipid distributions in animal bones after the solvent washed were similar between species, mainly dominated by saturated fatty acids of C 12:0 -C 30:0 , branched fatty acids of C 14:0 and C 16:0 and monounsaturated fatty acids of C 18:1 and C 22:1 .Cholesterol derivatives were also present in trace amounts (Additional file 1: Figure S2).These characteristics indicate the endogenous lipids had been preserved in bone fats and provide an essential assurance for the compound-specific stable isotope analysis.

Stable isotope composition of archaeological bones and pottery lipids
Ten archaeological animal bones, including Water buffalo (n = 3), Deer (n = 2), Small carnivore (n = 1), Wild boar (n = 1) and Freshwater fish (n = 3), from Neolithic sites of Zengpiyan and Dayan were selected as the isotopic baseline for compound-specific δ 13 C values of pottery residues.The δ 13 C values of palmitic (C 16:0 ) and stearic acids (C 18:0 ) extracted from bone lipids (presented in Additional file 1: Table S2) ranged from −30.0 to −26.6‰ for C 16:0 and −30.0 to −28.3‰ for C 18:0 (n = 10), indicating a widely C 3 based lipid origins.However, a more complicated isotopic fractionation related to lipid synthesis and different dietary metabolism among animal species should be considered when compared to isotopic data of lipid, collagen and apatite from the same bone [21,48].
The δ 13 C values of bone components from different species are shown in Fig. 5.The general depletion patterns of 13 C in the C 16:0 and C 18:0 fatty acids relative to collagen and apatite in terrestrial mammals are consistent with previous studies due to discrimination against 13 C during the metabolic pathway of lipid synthesis which resulted in more depleted δ 13 C values than other tissues, e.g., protein, carbohydrate and cholesterol [52,53].For most of the terrestrial mammalian bones in Fig. 5, including deer, wild boar and small carnivores, the δ 13 C values of fatty acid, collagen and apatite indicate a C 3 lipids origin with dietary digestion and biosynthetic process.Notably, subtle variation of δ 13 C in apatite (−11.3‰)relative to fatty acids (−29.8‰ for C 16:0 and −28.3‰ for C 18:0 ) and collagen (−21.9 ± 0.1‰) from wild boar may relate to a small intake of C 4 foods.This difference is particularly clear in the bone composition of water buffalo (Fig. 5 and Additional file 1: Table S2), in which collagen (−12.8 ± 3.7‰) and apatite (−8.3‰) indicated a more C 4 diet (enriched in 13 C), while the fatty acids (−27.0 ± 0.5‰ for C 16:0 and −28.9 ± 0.3‰ for C 18:0 ) show a depleted in 13 C by 14.2‰ (Δδ 13 C = C 16:0FA-COLL ) and 18.7‰ (Δδ 13 C = C 16:0FA-APAT ).This may be mainly reflected by variations of isotopic fractionation during carbon metabolic and biosynthesis in the organism [54,55].
Controlled animal feeding experiments have found significant differences (−3.3‰ to + 12.2‰) in δ 13 C values of tissues from animals fed on a diet of constant isotopic composition [53,56] and pointed out that the fractionation of 13 C from diet to tissue is not identical on different diets, such as the proportion of C 3 vs C 4 , or aquatic vs terrestrial in diet.Considering the habitat environments and foraging behavior of water buffalo, its food resource may contain terrestrial and aquatic C 3 and C 4 plants.This mixed diet could lead to a larger isotopic difference in these bone components, although sex, age and physical condition may also have secondary effects [53].
Carbon sources are complex in the freshwater ecosystem and can give rise to a wide range of isotopic ratios in freshwater environments.The δ 13 C values of bone components plotted in Fig. 5 shows a pattern that fatty acids are more depleted than their corresponding collagen and apatite values for freshwater fish, suggesting a typical lipid origin similar to the terrestrial C 3 pathway.As observed for most of the terrestrial mammals in Fig. 5, Samples without plant biomarkers have the molecular composition characterized by higher abundance of C 18:0 and lower abundance of C 16:0 [23].These samples showing a relatively narrow distribution with δ 13 C values ranging from −31.6 to −26.4‰ for C 16:0 and −31.9 to −28.0‰ for C 18:0 (n = 15), and therefore reflect that lipids were derived entirely from C 3 animals.Except for a few extreme data pointing to an increased C 4 resource, the δ 13 C values of most samples fall into the isotopic ranges between terrestrial C 3 and aquatic resources in Fig. 6a.However, aquatic lipid origin cannot be verified solely from isotopic data, due to some overlapping in the ẟ 13 C values of aquatic lipids and terrestrial C 3 lipids, especially aquatic biomarkers were not detected in these pottery samples [23].Based on the results of cooking experiments shown in Fig. 4 and Additional file 1: Table S1, potsherds that have been processed for aquatic stuffs under high heating conditions yielded aquatic biomarkers including isomers of APAAs (C 16 -C 22 ) and at least one isoprenoid fatty acids.Other potsherds that underwent low temperature conditions as well as raw aquatic materials had partial aquatic biomarkers, especially 4,8,12-TMTD.Thus, we suggest that if pottery had been used to process aquatic resources, partial aquatic biomarkers at least should be detected.Therefore, a predominantly terrestrial resource origin is convincing based on the current results.
This result is further supported in Fig. 6b by plotting Δ 13 C (δ 13 C 18:0 -δ 13   in biosynthetic process between non-ruminant and ruminant fat types [57].From the data shown in Fig. 6b, the classification of ruminant dairy fats from adipose fats and non-ruminant adipose fats can be separated through Δ 13 C values of animal bone fatty acids by displaying the offsets lower than 0.0‰ and −3.3‰ [58], and provide the reference proportion that approximately 40% (22/52) of the ruminant resources and 60% (30/52) of non-ruminant resources, including terrestrial plants, were being processed in these early pots, which is consistent with archaeobotanical and zooarchaeological studies in the site area [19,44,59].
One sample fell within the reference range of ruminant dairy fats with a clear C 4 signal which was beyond our expectation.However, no clear lipid components of dairy fats, such as triacylglycerols, were identified.Low abundance of short-chain (C 8:0 , C 12:0 -C 14:0 ) fatty acids detected in this sample [23], were probably from the degradation of free fatty acids through thermal cracking, or by bacterial action [60].Given that only one sample falls into this range and no archaeological evidence related to dairy products or activities have been found during the Neolithic times in the site area, we suggest that dairying was unlikely to take place during the early stage of the Neolithic in this region, although occasional utilization cannot be ruled out.
Overall, isotope results combined with molecular information demonstrate that a mixture of terrestrial C 3 and C 4 resources, including plants, ruminants and nonruminants, dominate the utilization history of the early pottery vessels in both sites.However, the evidence of frequent processing of aquatic products within these pots remain less clear.The findings suggest that despite the abundant edible aquatic resources in the study area, the Neolithic people had a clear preference for processing terrestrial resources via pottery vessels, and indicated that the culinary practice was homogeneous in this region and was persistent throughout Neolithic times.

A new understanding on the resource utilization pattern through pottery during Neolithic times in South China
The emergence and utilization of pottery in South China have been widely debated, mainly focusing on the function behind its invention [5,10,61].A persistent use of terrestrial C 3 and C 4 resources through our systematically stable isotope analysis of pottery lipids together with plenty of terrestrial biomarkers (such as a series of diterpenes, triterpenes, sterols, cholesterols and longchain ketones) in previous studies, suggests these Neolithic vessels have mainly been used to process terrestrial plants and animals, rather than aquatic resource, which is in contrast with the traditional view.Moreover, our new results also different from that in Japan and the Eastern Russia, where most of the early pots yielded aquatic signals, suggesting the extensive exploitation of aquatic resources [12,13].The pottery use pattern in this study represents a unique phenomenon among prehistoric shell mounds, providing new insight into the local traits of ecological tradition and social evolution during Neolithic times in South China.
In the ecological context of South China, the subtropical and tropical environments, compared to other regions in East Eurasia, promote the rich and unique endemic flora and fauna locally [62].The forests landscape reshaped from the onset of the Pleistocene and the hunting tradition of the large mammals which documented by the paleolithic phases in the study sites can be viewed as the ecological tradition for prehistoric people on one hand [59,62], these settings in turn provide the environmental background for pottery production during the late Pleistocene in south China on the other hand.From this point of view, though aquatic resources were included in the dietary spectrum of the Neolithic people, the continuous use of the terrestrial resource might be explained by the adaptation to the ecological traditions that inherited by generations of the local Paleolithic people.Besides, plant resources, particularly those with higher sugar contents, have been reported as the potential factor resulting in human oral diseases during Paleolithic times [63].This also raises another issue that related to the need for cooked food.Compared to animals, some plants can be processed to remove toxicity and became easily digestible [64].Almost 70% (37 of 52) of the pottery samples that yielded plant biomarkers in the study sites, particular in Dingsishan site, together with the population size which reflected by total number of the human skeleton remains in the Dingsishan site (no less than 300 individuals [65]), suggests some kinds of terrestrial resources probably became part of the staple foods with social development and population growth during Neolithic times.
Furthermore, pottery itself, as a symbol of cultural evolution, represents a materialized form of public activities embedded in prehistoric social life and constantly linked with the transformation of prehistoric social structures, such as farming [5].Although the transition to agriculture in South China is less clear than in North China and West Asia, where archaeological evidence made the hunter-gather to farmer connection much easier to trace [66,67].We suggest that such a case can also be reflected by pottery residue analysis.As shown in Fig. 7, our isotopic data combined with the published data from the Zengpiyan site (Additional file 1: Table S4), provide a remarkable diachronic change from extensive to the intensive tendency of resources used during different Neolithic stages, particularly within the early and middle Neolithic stages.Broader distribution of δ 13 C values demonstrates the resources derived from a wide range of C 3 and C 4 ecotone.This was further supported by the evidence of zooarchaeology and stable isotope analysis of human and animal bone apatite [19,68], which showed a variety of niches around the study areas, including dense tropical forests, open C 3 woodland with patches of C 4 grassland [68], and such ecosystem may expand the hunting scope and point to relatively high or seasonal mobility for hunting and gathering food resources during these two stages.
By contrast, the intensive use of C 3 resources (Fig. 7) during the late Neolithic may closely map onto the plant and animal management or a more sedentary lifestyle.The scale of human selection on plant species in this stage can be seen in the number of which are common in the subtropical and tropical regions, including root plants (e.g.Dioscorea, Araceae, Taccaceae, Alismataceae), Coix and Triticeae Dumort [69,70].Despite domesticated animals being fewer in species [19], the feeding strategy of other animals may also occur, such as deer, water buffalo and wild boar which were found in large quantities in the study sites [19,59].Furthermore, an increased number of processing tools, such as smaller pottery vessels and grinding and digging tools, on a regional scale during the late Neolithic [71], also indicates that intentional food production or cultivation could arise in this stage [71,72].As illustrated in Fig. 7, compared to a mixture of C 3 and C 4 resources in the previously two stages, the intensive use of certain C 3 plants and animals in late Neolithic, implied that the early mobile groups were possibly replaced by a new social structure of sedentary groups that build on food production during this time, laying the foundation for the subsequent agricultural development and social evolution in prehistoric South China.
In summary, resource processing is undoubtedly a major driving force for pottery invention and usage in South China.The ecological tradition and the prehistoric social evolution path hidden implied by the resource utilization pattern are worthy of further investigation.Perhaps it is owing to the unique ecological environments in South China, where the resource utilization patterns of early pottery and the subsequent social development behind pottery production are different from those of other regions in East Eurasia during Neolithic times.
The absence of aquatic biomarkers in these pottery samples was somewhat surprising, since evidence from faunal and stable isotope of human and animal bone collagen shows clear consumption of aquatic foods.A similar phenomenon was also found in the previous work [24].Both findings seemingly contradict the traditional assumption that the pottery in these sites was primarily used for processing aquatic resources [16,18].For this issue, we suggest the following possibilities.
Firstly, although the samples we selected in this study represent the typical pottery container for cooking [18,43], the small sample size or type may have clouded the detection of aquatic lipids.Secondly, lipids were extracted in every potsherd, and their concentrations were all above the minimum amount required for interpretation (> 5 μg g −1 ), with the highest concentration up to 6666.7 μg g −1 , indicative of a relatively good lipid preservation [23].Even considering the lipid degradation during the long-term burial, it is still difficult to explain the presence of a series of terrestrial lipids but only absence of the aquatic ones.In fact, these compounds remain relative stable after burial, and were commonly found in a wide range of archaeological pottery [13,22,50,[73][74][75][76][77], though sometimes in trace amounts [22].Thus, another possible explanation may be related to the cooking conditions.The cooking experiments in this study suggested that certain kinds of aquatic biomarkers could be detected in pottery at lower firing conditions.However, if these products were cooked in a shorter time, such as gently boiled in a few minutes, the firing conditions may not enough to form various aquatic biomarkers in pottery.Indeed, a large collection of well-preserved mollusk shell remains in the study sites may indicate that they have not undergone intense burning.
The last possibility may be related to dietary habits or cooking tools for aquatic resources during Neolithic period.Eating raw food has been viewed as a dietary habit in ancient and modern societies in South China, thus the possibility of eating certain kinds of aquatic food without processing, such as fish, should not be completely excluded.In addition, prehistoric people may choose other vessels to cook aquatic foods.Forest resources, such as bamboo and wood, are widespread in South China and adjacent Southeast Asia, and they are still common tools for the daily cooking of aquatic food in this vast area today [78,79].Although it is presently unclear how aquatic resources were processed and cooked in South China during the Neolithic times, the lack of evidence for aquatic products in our pottery samples suggests that they were probably prepared in other ways.In future, organic residue analysis on larger archaeological assemblages including pottery vessels, stone tools and other possible artifacts may help understand the culinary practices then.

Conclusion
This study conducted a systematic analysis of archaeological pottery residues and bone lipids from Dayan, Dingsishan and Zengpiyan sites in Guangxi province of China, revealing the extensive processing of terrestrial C 3 and C 4 resources in pottery vessels throughout the Neolithic times in this area.The persistent use of terrestrial resources and the remarkable diachronic changes from extensive use of C 3 and C 4 resources to the intensive use of certain kinds of C 3 resources in different Neolithic stages further suggest that the resource utilization patterns in early pottery may be closely related to the ecological tradition and the subsequent social evolution in prehistoric South China.
Notably, the divergence between the findings of zooarchaeology as well as stable isotope analysis in these sites, which suggest an extensive exploitation of aquatic resources, and the lipid residue evidence from the potsherds in our study, which strongly indicates the processing of terrestrial resources, offers a new insight into the pottery use and culinary practice during Neolithic times in South China.However, a larger investigation is needed to further understand this discrepancy.
In fact, the commonality of pottery utilization built on cumulative archaeological evidence not only reflects the dietary innovation at the onset of the Neolithic period across different regions, but also implies the potential difference of culinary traditions under diverse environmental contexts.A combined research of ecological, economic and social factors based on organic residue analysis is of significant importance for exploring the pottery invention and utilization in ancient China and beyond.

Fig. 1
Fig. 1 Map showing the location of the study sites within South China

Fig. 2
Fig. 2 Early pottery samples from Dayan and Dingsishan sites: a-d the earliest pottery sherd from Dayan site (a pottery sample in this study, b, c general view of the earliest pottery sherd, d restoration image of the earliest pottery with the sample position, pictured from Xinrong Jiang), e-l pottery sherds from Dingsishan site (e, f early Neolithic, g-j middle Neolithic, k, l late Neolithic)

Fig. 3
Fig. 3 Total ion current (TIC) and m/z 105 of the cooking shellfish.Bottom: fatty acids distributions without APAAs heated at 230 ℃ for 5 h.Top: Mass chromatograms of base peak (m/z 105) with molecular ions of C 18 -C 22 APAAs heated at 280 ℃ for 5 h

Fig. 4
Fig. 4 TIC chromatograms from cooking shellfish (top) and plants (bottom).Top: APAAs based on ions m/z 105, 262, 290, 318 and 346 with carbon length C 16 -C 22 .4,8,12-TMTD based on ions m/z 87 and 270.Phytanic acid was monitored with ions based on m/z 101.Bottom: APAAs with carbon length C 18 and phytanic acid were also detected after cooking plants C 16:0 ) values against δ 13 C 16:0 values.The offset (Δ 13 C) reflected the different routing of C 18:0 and C 16:0 during lipid synthesis in ruminant dairy fats and adipose fats, emphasizing the difference

Fig. 5
Fig. 5 The offsets of δ 13 C values between bone components of different animals

Fig. 6
Fig. 6 Plot of δ 13 C values (a) and the offset values (Δ 13 C, b) of C 16:0 and C 18:0 in archaeological bone and pottery lipids from Zengpiyan, Dayan and Dingsishan sites

Fig. 7
Fig. 7 Diachronic change of δ 13 C values of C 16:0 and C 18:0 in Neolithic South China