Front. Ecol. Evol. Frontiers in Ecology and Evolution Front. Ecol. Evol. 2296-701X Frontiers Media S.A. 10.3389/fevo.2022.840199 Ecology and Evolution Original Research Integrating Lipid and Starch Grain Analyses From Pottery Vessels to Explore Prehistoric Foodways in Northern Gujarat, India García-Granero Juan José 1 2 * † Suryanarayan Akshyeta 3 4 5 † Cubas Miriam 6 7 Craig Oliver E. 7 Cárdenas Marc 8 Ajithprasad P. 9 Madella Marco 5 10 11 * 1HUMANE–Human Ecology and Archaeology Research Group, Department of Archaeology and Anthropology, Institució Milà i Fontanals, Spanish National Research Council (IMF-CSIC), Barcelona, Spain 2School of Archaeology, University of Oxford, Oxford, United Kingdom 3Department of Archaeology, University of Cambridge, Cambridge, United Kingdom 4Université Côte d’Azur, CNRS, CEPAM, Nice, France 5CaSEs Research Group, Department of Humanities, Universitat Pompeu Fabra, Barcelona, Spain 6Departamento de Historia y Filosofía, Universidad de Alcalá, Facultad de Filosofía y Letras, Área Prehistoria, Alcalá de Henares, Spain 7Department of Archaeology, BioArCh, University of York, York, United Kingdom 8Faculty of Archaeology, Universiteit Leiden, Leiden, Netherlands 9Department of Archaeology and Ancient History, Maharaja Sayajirao University of Baroda, Vadodara, India 10ICREA—Institució Catalana de Recerca i Estudis Avançats, Barcelona, Spain 11School of Geography, Archaeology and Environmental Studies, The University of the Witwatersrand, Johannesburg, South Africa

Edited by: Xiaoyan Yang, Institute of Tibetan Plateau Research (CAS), China

Reviewed by: Yuichi Naito, Nagoya University, Japan; Subir Bera, University of Calcutta, India

 *Correspondence: Juan José García-Granero, jjgarciagranero@imf.csic.es Marco Madella, marco.madella@upf.edu

†These authors have contributed equally to this work

This article was submitted to Paleontology, a section of the journal Frontiers in Ecology and Evolution

 16 03 2022 2022 10 840199 20 12 2021 22 02 2022 Copyright © 2022 García-Granero, Suryanarayan, Cubas, Craig, Cárdenas, Ajithprasad and Madella. 2022 García-Granero, Suryanarayan, Cubas, Craig, Cárdenas, Ajithprasad and Madella

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

This study attempts a holistic approach to past foodways in prehistoric northern Gujarat, India, by considering evidence of food production, distribution, preparation and consumption. We present here the results of a pilot residue study, integrating lipid and starch grain analyses, conducted on 28 ceramic vessels from three Chalcolithic/Harappan settlements (c. 3300–2000 cal. BC) in northern Gujarat, which are discussed in the light of previous evidence of plant and animal acquisition and preparation strategies in this region. We aim to explore how the prehistoric inhabitants of northern Gujarat transformed ingredients into meals, focusing on how different foodstuffs were processed. When assessed on their own, the lipid and compound-specific isotopic data suggest that animal fats were primarily processed in ceramic vessels, specifically non-ruminant fats. However, lipid residue analysis favors the detection of fat-rich animal products and is often unable to disentangle signatures resulting from the mixing of plant and animal products. The incorporation of starch grain analyses provides evidence for the processing of a range of plants in the vessels, such as cereals, pulses and underground storage organs. Together, the results provide a holistic perspective on foodways and a way forward in overcoming preservational and interpretational limitations.

 food culinary practices archeology lipids starch grains South Asia Ministerio de Educación, Cultura y Deporte 10.13039/501100003176 Agència de Gestió d'Ajuts Universitaris i de Recerca 10.13039/501100003030 Ministerio de Economía y Competitividad 10.13039/501100003329 Ministerio de Ciencia e Innovación 10.13039/501100004837 H2020 Marie Skłodowska-Curie Actions 10.13039/100010665 Ministerio de Ciencia, Innovación y Universidades 10.13039/100014440 Introduction

Archeological research has long focused on the cultural, social and economic practices concerning food production and consumption (often referred to as “foodways”; Staller and Carrasco, 2009; Peres, 2017), especially since the systematic recovery of macrobotanical and microfaunal remains through flotation techniques became widespread in the 1970s. Through the analysis of plant and animal remains, often combined with the analysis of food-related artifacts such as ceramics and grinding stones, archeologists have explored issues such as the social division of labor during food production (e.g., Bolger, 2010), distribution (e.g., Welch and Scarry, 1995) and preparation (e.g., VanDerwarker and Detwiler, 2002), intra- and inter-cultural culinary preferences (e.g., Kirch and O’Day, 2003) and the role of food in social aggregation and cooperation (e.g., Bray, 2003), among other topics.

Despite the diversity of methodological and theoretical approaches used to explore prehistoric foodways and the abundant literature discussing this topic, few studies have considered both plant and animal resources holistically (e.g., Spielmann, 2002; Bogaard et al., 2009; Kansa et al., 2009; Twiss et al., 2009; Ivanova et al., 2018; Gaastra et al., 2019; McClatchie et al., 2019; Dunne et al., 2021). However, the acquisition, preparation and consumption of both plants and animals are firmly tied together in integrated agro-pastoral systems and, therefore, both resources need to be considered together in order to explore prehistoric foodways and pursue a better understanding of how food systems operated in the past.

In the last few decades, past food preparation and consumption activities have been explored through the analysis of residues from archeological artifacts and human dental calculus. Chemical (lipids and proteins; e.g., Craig et al., 2015; Hendy et al., 2018) and microbotanical (starch grains and phytoliths; e.g., Lu et al., 2005) analyses have greatly contributed to our understanding of the use of plant and animal ingredients in prehistoric cuisines, often highlighting the presence of foodstuffs not detected through conventional archeobotanical and zooarcheological methods (e.g., Högberg et al., 2009; Salque et al., 2012; Saul et al., 2013). Phytoliths are often found on non-edible plant parts but can also represent the use of taxa often underrepresented in conventional archeobotanical analyses, such as spices (Saul et al., 2013); whereas starch grains are mainly produced in the edible plants parts (seeds, fruits, underground storage organs, etc.) and are usually regarded as a direct evidence of the consumption of starchy plants, including, among others, cereals, pulses and tubers (Torrence and Barton, 2006). In parallel, the analysis of lipids from pottery vessels has thrown light on our understanding of the exploitation of animal fats and the origin of dairying practices (e.g., Craig et al., 2005; Evershed et al., 2008), as well as on the use of leafy vegetables, plant oils (Charters et al., 1993, Charters et al., 1997; Copley et al., 2001, 2005; Dunne et al., 2016) and apicultural products (Roffet-Salque et al., 2015).

Individually, phytoliths, starch grains and lipids provide valuable information on the consumption of plant and animal resources. However, their independent analysis might underestimate the importance of certain resources. For example, fat-rich animal products are more easily identified in lipid residue analysis, while the contribution of lipid content of plant-based products, generally at least tenfold lower, is likely masked by animal fats (Charters et al., 1995; Hammann and Cramp, 2018; Miller et al., 2020). Plant waxes, sterols and seed oils have been identified in pottery vessels, but generally in archeological settings with good organic preservation (e.g., Copley et al., 2001, 2005; Dunne et al., 2016). Additionally, there are only a few established biomarkers for specific plant products (brassica leaf waxes, maize and some millet species; Charters et al., 1997; Reber and Evershed, 2004; Heron et al., 2016). These preservational and interpretational biases can potentially skew our interpretation of past resource- and vessel-use in respect to plants, rendering the use of plant products in pottery vessels “invisible”.

At present, microbotanical remains and lipid residues are seldom analyzed as part of the same study, and when this happens they often come from different vessels, thus impeding an effective integration of the results (see for example, microbotanical and lipid residue analyses conducted at the Neolithic site of Stavroupoli, in northern Greece; García-Granero et al., 2018; Whelton et al., 2018). The integrated analysis of lipid residue and microbotanical proxies has the potential to widen the spectrum of identifiable food resources and examine differential pathways of the processing and consumption of food (Kooiman et al., 2021). Moreover, an integrated analysis can help overcome the interpretative limitations of individual proxies.

We present here the results of a pilot study aiming at assessing the potential of an integrated lipid residue and microbotanical analysis to explore prehistoric foodways in Chalcolithic/Harappan northern Gujarat, western India (Figure 1). The materials analyzed in this study comprise 28 ceramic vessels potentially used for food preparation and consumption from three settlements in northern Gujarat: Chalcolithic deposits from Datrana (c. 3300–3000 cal. BC) and Loteshwar (c. 2700–2300 cal. BC), and late Urban Harappan deposits from Shikarpur (c. 2200–2000 cal. BC)—see Table 1 for a summary of the archeological contexts and Supplementary Material for a more detailed description of the full stratigraphy of the sites. The results are discussed in the light of previous evidence of plant and animal acquisition and preparation strategies in this region to explore how the Chalcolithic/Harappan inhabitants of northern Gujarat transformed ingredients into meals, thinking through the potential culinary pathways of different foodstuffs.

Map of Gujarat showing the case studies (red triangles) and other archeological sites mentioned in the text (green circles). Background map prepared by Francesc C. Conesa.

Summary of the available archeological information for the deposits where the samples analyzed in this study were recovered.

Site Chronology of the analyzed samples Site description Ceramics Lithics Other material culture References
Datrana c. 3300–3000 BC Dune occupied by agro-pastoral groups focused on the production of lithic blades Mostly pre-prabhas, very little Early Harappan Sindh and Anarta Tools and debitage, mostly made of chalcedony, some Rohri chert blades, grinding tools Carnelian beads Ajithprasad, 2002, 2011; Gadekar et al., 2013; Rajesh et al., 2013b; García-Granero et al., 2017a
Loteshwar c. 2700–2300 BC Dune seasonally occupied by semi-nomadic agro-pastoral groups Mostly Anarta, very little Harappan Tools and debitage, mostly made of chalcedony and chert, grinding tools Terracotta, shell and copper objects, steatite and semi-precious stone beads Rajesh et al., 2013a; Gadekar et al., 2014a; García-Granero et al., 2016
Shikarpur c. 2200–2000 BC Fortified rural settlement Mostly Classical and Sorath Harappan, some Anarta Chert blades and cores (including Rohri chert), grinding tools Terracotta, shell and copper objects, steatite and semi-precious stone beads Bhan and Ajithprasad, 2009; Gadekar et al., 2014b; Chase et al., 2020
Area of Study

For the purposes of this study, we define northern Gujarat broadly, encompassing the northern part of mainland Gujarat (often referred to as North Gujarat), the island of Kachchh and the northern part of the peninsula of Saurashtra, a region characterized by a semiarid climate (400–600 mm of average annual precipitation), with most of the rainfall occurring during the Indian summer monsoon (June-September). Several ceramic traditions co-exist in northern Gujarat during the Chalcolithic/Harappan period, including the Anarta (c. 3700–2250 cal. BC), the Pre-Prabhas (c. 3300–2900 cal. BC) and the Classical and Sorath Harappan (c. 2500–1900 cal. BC). Anarta pottery has been recovered in different proportions from over 60 prehistoric sites, but it is most common in seasonal camps occupied by semi-nomadic agro-pastoralists, such as Loteshwar (Ajithprasad and Sonawane, 2011; Rajesh et al., 2013a), whereas Pre-Prabhas pottery has only been recovered at Datrana, a lithic blade workshop (Ajithprasad, 2011; Gadekar et al., 2013; Rajesh et al., 2013b), and a few other sites (Rajesh et al., 2018). Classical Harappan pottery of the Indus Civilization is found mostly in walled settlements, with the characteristic Indus city plan and associated material culture, ranging from villages such as Shikarpur (Bhan and Ajithprasad, 2009) to major urban centers such as Dholavira (Bisht, 2015). Finally, Sorath Harappan pottery is the dominant regional ceramic tradition associated with the Harappan sites (mostly) of Saurashtra (see e.g., Farooqui et al., 2013, p. 2632 and references therein).

Recent archeobotanical and zooarcheological research in northern Gujarat has provided extensive evidence for the production and distribution of plant and animal resources during the Chalcolithic and Harappan period. The analysis of macroscopic plant remains (charred and mineralized fruits and seeds) from Anarta and Harappan sites shows that prehistoric populations relied on the cultivation of monsoon-adapted crops native to South Asia, particularly small millets such as little millet (Panicum sumatrense), browntop millet (Brachiaria ramosa), bristly foxtail (Setaria verticillata), green foxtail (S. viridis), Kodo millet (Paspalum scrobiculatum) and barnyard millet (Echinochloa spp.), as well as tropical pulses such as horsegram (Macrotyloma uniflorum), mung bean (Vigna radiata) and black gram (V. mungo) (Fuller, 2011; Pokharia et al., 2011, 2017; García-Granero et al., 2015, 2016). Sesame (Sesamum indicum) was also recovered from Loteshwar (García-Granero et al., 2016) and a few Harappan sites in the region (Pokharia et al., 2011, 2017). Small amounts of barley (Hordeum vulgare) and, to a lesser degree, free-threshing wheat (Triticum aestivum/turgidum) are also recovered from all studied sites; these were probably not cultivated locally but traded in from other Indus regions more suitable for the cultivation of winter crops (García-Granero, 2015). Wild plants are also normally recovered from Chalcolithic/Harappan sites in northern Gujarat but there exists no clear evidence of their use for human consumption, with the possible exception of Egyptian crowfoot grass (Dactyloctenium aegyptium) at Loteshwar (García-Granero et al., 2016).

There is abundant evidence for animal husbandry in Chalcolithic and, especially, Harappan northern Gujarat, where hunting appears to have played a minor role, particularly during the Harappan period. Cattle (Bos sp.) herding seems to have been the main animal-related activity during the Anarta period (Patel, 2009), later complemented with sheep (Ovis aries) and goats (Capra hircus) in Harappan settlements—and, possibly, water buffaloes (Bubalus bubalis) and pigs (Sus sp.), the domestic status of which has not been clarified to date (Thomas et al., 1996, 1997; Chase, 2010, 2014; Joglekar and Goyal, 2011; Goyal, 2013; Joglekar et al., 2013). The domestic animal slaughter patterns at Harappan Bagasra and Shikarpur suggest that cattle and buffalo were generally kept for secondary products and/or animal traction prior to consumption, whereas sheep and goats were raised primarily for meat (Chase, 2010, 2014). Fishing, particularly of marine habitats, also seems to have been an important activity in Harappan Gujarat (Abhayan, 2016). Judging by the abundance of otoliths and other fish bones at archeological sites in northern Gujarat, fish must also have been an important resource for prehistoric populations, although its role in their subsistence is currently unknown due to the scant attention fish remains have received in South Asian zooarcheology (Abhayan et al., 2016).

Husbandry practices were further explored at Harappan Bagasra, Shikarpur, Jaidak and Kotada Bhadli through stable isotope analyses (87Sr/86Sr, δ13C, and δ18O) of animal tooth enamel. At Bagasra, strontium ratios suggest that most sheep and goats were raised locally, whereas around half of the analyzed cattle were raised further afield, in several locations throughout central Saurashtra (Chase et al., 2014, 2018). Carbon and oxygen values further suggest that cattle were predominantly fed agricultural fodder (millet cultivation by-products), as shown by the presence of C4 plants in their diet, whereas sheep and goats consumed a mixed diet made of both agricultural fodder and wild plants (mixed C3–C4 diet) (Chase et al., 2014). Results from nearby Shikarpur and Jaidak (Chase et al., 2020) seem to confirm the patterns observed at Bagasra. Δ13C-values from Kotada Bhadli, however, showed that both cattle/buffaloes and sheep/goats consumed a mixed C3–C4 diet (Chakraborty et al., 2018), thus suggesting that livestock management practices were not uniform throughout Harappan northern Gujarat.

Despite the relatively rich literature discussing food production and distribution in Chalcolithic northern Gujarat, less attention has been devoted to food preparation activities, and the remains of food consumption have been largely unexplored. Evidence for the pre-consumption preparation of plant resources comes from the analysis of microbotanical remains (starch grains and phytoliths) from grinding stones from Loteshwar, Datrana and Shikarpur, which show that these were mostly used to grind small millets and pulses, with a minor presence of other resources such as wheat/barley and underground storage organs (García-Granero et al., 2015, 2016, 2017a,2017b). Evidence for the preparation of animal resources comes from the analysis of cut-marks in cattle/buffalo and sheep/goat from Bagasra, which showed different butchery practices between those residing within and without the walled settlement (Chase, 2012). Using a different approach, Goyal (2017) suggested that preparation via roasting was more common in wild animals, particularly deer, than in domestic animals at Kanmer. Presence of burnt marks and the scarcity of cut marks in the fish assemblages from Bagasra, Kanmer and Shikarpur suggest that at least some fish resources were also prepared via roasting (Abhayan, 2016: 298).

 Materials and Methods Sample Details

A total of 28 pottery vessels were analyzed in this study: 11 from Datrana (excavated in 2010), six from Loteshwar (excavated in 2009) and 11 from Shikarpur (excavated in 2012). All samples from Datrana come from bases of medium or large pots/dishes characterized as Pre-Prabhas. Samples from Loteshwar come from body sherds of medium pots, five characterized as Anarta and one characterized as probably Sorath Harappan (but definitely not Anarta). Finally, samples from Shikarpur come from bases of a variety of Classical Harappan vessels, including four pots/vases, four goblets, two pots/jars and one flat platter (a detailed description of the vessels can be found in Supplementary Table 1).

After retrieval from the archeological matrix potsherds were wrapped in aluminum foil and stored at the Department of Archeology and Ancient History, M.S. University of Baroda, Vadodara, India, where they were sampled for lipid residue and microbotanical analyses in November 2013. Sampling took place in a controlled environment—a closed room with no airstream. A small portion of each sherd (c. 1 cm2) was detached from the main sherd, wrapped in aluminum foil and sent to the BioArch Laboratory, University of York, United Kingdom, for lipid residue analyses.

Microbotanical residue recovery consisted of a two-step process in which the outer layer of sediment was first dry brushed from the inner surface of the vessel (dry sample), and then the inner layer of sediment was brushed with deionized water (wet sample). By removing the outer layer of sediment the likelihood of contamination from the burial environment decreases considerably (Hart, 2011), so our analysis focused on the wet samples. Microbotanical samples were transferred to the BioGeoPal Laboratory, IMF-CSIC, Barcelona, Spain, where wet samples were immediately dried at 40°C and stored. Gloves were not used during the sampling of potsherds or the extraction and analysis of starch grains to prevent starch contamination (Crowther et al., 2014).

 Lipid Residue Analysis: Extraction and Analytical Protocol

Lipid extracts were obtained and prepared from the pottery sherds (n = 28) at the BioArch Laboratory using previously reported protocols of extraction and methylation (Craig et al., 2013). Each sherd was first cleaned with a high-speed drill to eliminate any surface contamination. Ceramic was then drilled from the interior surface (1 g). The ceramic powder was weighed and sealed in glass vials prior to all analyses. Methanol (4 ml) was added to the powdered pottery and the mixture was sonicated for 15 min. Then, sulfuric acid (800 μl) was added and heated at 70°C for 4 h. Lipids were extracted from centrifuged pottery powder with n-hexane (3 × 2 ml) and dried under a gentle stream of N2. Internal standard (10 μl of hexatriacontane C36:0) was added in all the samples prior analysis to quantify the relative abundance of lipids. All samples were analyzed by Gas Chromatography-Flame Ionisation Detection (GC-FID) (n = 28), Gas Chromatography-Mass Spectrometry (GC-MS) (n = 28) and Gas Chromatography-Combustion-Isotope Ratio Mass Spectrometry (GC-C-IRMS) (n = 27).

GC-FID was carried out using an Agilent 7890S gas chromatograph (Agilent Technologies, Cheadle, Cheshire, United Kingdom). The sample (1 μl) was injected into the GC at 300°C with a splitless injector, using helium as carrier gas (2 ml min–1). The GC column was a polymide coated fused-silica DB-1 (15 m × 320 μm × 0.1 μm; J&W Scientific, Folsom, CA, United States). The GC oven was set at 100°C for 2 min, then increased by 20°C min–1 until 325°C, where it was held for 3 min.

GC-MS was carried out using an Agilent 7890 B Series Gas Chromatograph attached to an Agilent 5977 B Mass Spectrometer with a quadrupole mass analyzer (Agilent technologies, Cheadle, Cheshire, United Kingdom). All samples were initially screened using a split/splitless injector in splitless mode, which was maintained at 300°C. The GC carrier gas was helium, configured at a constant flow rate of 1 ml min–1. The column (HP-5MS) was coated with 5% phenyl-methylpolysiloxane (30 m × 0.25 mm × 0.25 μm; Agilent technologies, Cheadle, Cheshire, United Kingdom). The oven temperature was set at 50°C for 2 min, then raised by 10°C min–1 until 325°C was reached, where it was held for 15 min until the end of the run. The ionization energy of the mass spectrometer was 70 eV and spectra were obtained in scanning mode between m/z 50 and 800.

In order to assess the presence of miliacin, each Total Ion Chromatogram (TIC) was scanned for m/z 189, m/z 204, m/z 231, m/z 425, m/z 440, corresponding to miliacin fragmentation. Additionally, ω-(o-alkylphenyl) alkanoic acids (APAAs) were screened by searching the TIC for the molecular ions (M +) for APAAs of C16 –C22 at m/z 262, 290, 318, and 346 and the fragment ion of the base peak m/z 105.

The stable carbon isotope values of palmitic and stearic FAMES were measured by GC-c-IRMS, using an Agilent 7890B series GC (Agilent Technologies, Santa Clara, CA, United States), linked by an Isoprime GC5 interface (Isoprime Cheadle, United Kingdom) to an Isoprime 100 (Isoprime, Cheadle, United Kingdom) and to an Agilent 5975C inert mass spectrometer detector (MSD). Samples were re-dissolved in hexane and 1 μl was injected into DB-5MS ultra-inert fused-silica column (60 m × 0.250 mm × 0.25 μm, J&W Scientific, Folsom, CA, United States). The temperature program was 50°C for 0.5 min, 25°C min–1 to 175°C, 8°C min–1 to 325°C, isothermal hold for 20 min. The carrier gas used was ultra-high purity grade helium (3 ml min–1). The gas flow eluting from the column was split into one stream that was directed to the MSD for compound identification, and another stream that was directed through the CuO furnace tube at 850°C to convert all the carbon species to CO2. Ion intensities (44, 45, and 46 m/z) of eluted products were recorded and the corresponding 13C/12C ratios were computed. Data was analyzed using IonVantage and IonOS software (Isoprime, Cheadle, United Kingdom) and the samples were compared with a standard reference gas (CO2) of known isotopic composition. The results are expressed in per mill (‰) relative to an international standard (VPDB). Within each batch, a mixture of n-alkanoic acid ester standards of known isotopic composition (Indiana standard F8-3) was used to check instrument accuracy (± 0.3‰) and precision on repeated measurements (± 0.5‰). Each sample was measured in duplicate. The resulting data were corrected to account for methylation of the carboxyl group through comparisons with a C16:0 and C18:0 fatty acid standard of known isotopic composition that were processed with each batch under identical conditions. As only a couple modern reference dairy fats from South Asia are available (Craig et al., 2005), data obtained from other published modern ruminant, dairy, non-ruminant fats, marine fats and plant oils (Copley et al., 2003; Craig et al., 2005; Spangenberg et al., 2006; Gregg et al., 2009; Outram et al., 2009; Steele et al., 2010; Dunne et al., 2012) were compared with the obtained results.

 Microbotanical Remains

Processing and analysis of the wet samples from potsherds took place at the BioGeoPal Laboratory, IMF-CSIC, Barcelona, Spain. Samples were chemically processed for the extraction of starch grains and phytoliths following the protocols described in García-Granero et al. (2017b), which involve the use of 5% sodium hexametaphosphate [(NaPO3)6] to deflocculate clays, sodium polytungstate [Na6(H2W12O40)] at a specific gravity of 1.8 g/cm3 to isolate starch grains, 5% hydrochloric acid (HCl) to dissolve calcium carbonate, 33% oxygen peroxide (H2O2) to oxidize organic matter and sodium polytungstate at a specific gravity of 2.35 g/cm3 to isolate phytoliths. Samples were observed under plain and cross-polarized transmitted light in a Leica DM2500 optical microscope with a Leica DFC490 camera for microphotography. 10% of the starch residue from each sample was mounted in 50% glycerine and fully scanned at × 200 magnifications, whereas c. 1 mg of the phytolith residue from each sample was mounted in Entellan® and scanned at × 630 magnifications until 250 identifiable single-cell phytoliths were observed. Samples where less than 100 identified phytoliths had been encountered after scanning 10% of the slide were considered sterile and the analysis did not proceed any further. Microbotanical remains encountered in archeological samples were compared with the modern plant reference collection housed at the BioGeoPal Laboratory (Supplementary Figure 1).

 Results Lipid Preservation and Molecular Evidence (Gas Chromatography-Flame Ionisation Detection and Gas Chromatography-Mass Spectrometry)

Lipid yields ranged from 4.6 to 339.1 μg/g (x¯ = 74.5 μg/g) (Table 2). All but one sample (DTR116) had lipid yields above 5 μg/g (n = 27, 96%). The lipid yields are relatively higher than those reported from other lipid residue studies in South Asia, such as at Kotada Badli, a Sorath Harappan settlement in Gujarat (Chakraborty et al., 2020, using a slightly different extraction protocol; Correa-Ascencio and Evershed, 2014), and Harappan sites in northwest India (Suryanarayan et al., 2021, using the same extraction protocol as in this study). Comparisons of lipid yields suggest there are no differences in lipid yields across sites {Kruskal-Wallis test of effect of site [χ2(2) = 2.6, p = 0.27]}.

Results of the chemical and starch grain analyses from pottery vessels from Datrana, Loteshwar and Shikarpur.

Molecular evidence GC-c-IRMS Starch grains
Lipid conc. (μ g/g) P/S ratio δ 13C C16:0 δ 13C C18:0 Δ 13C (C18:0-C16:0) Small millets Job’s tears Wheat/barley Pulses Ginger USO (other) UNID starch
Datrana
DTR101 130.4 0.8 –28.3 –25.8 2.5 0 0 0 0 0 0 0
DTR104 107.5 0.9 –29.6 –29.0 0.6 0 0 0 0 0 0 0
DTR105 43.0 1.1 –25.4 –28.1 –2.7 0 0 2 0 0 0 1
DTR106 67.0 0.7 –25.0 –23.5 1.5 0 0 0 0 0 0 0
DTR107 41.7 0.9 –29.7 –28.5 1.2 0 0 0 0 0 0 0
DTR109 17.4 0.8 –30.2 –29.8 0.4 0 0 0 0 0 0 0
DTR113 30.0 1.1 –30.0 –29.6 0.3 0 0 1 0 0 0 0
DTR115 9.2 0.9 –29.8 –29.6 0.2 0 0 0 0 0 0 0
DTR116 4.6 – – – – 0 0 256 0 0 0 0
DTR117 94.8 0.9 –30.0 –29.8 0.2 0 0 0 0 0 0 0
DTR120 27.9 1.2 –24.5 –26.0 –1.5 1 0 0 0 1 0 1
Total Datrana −⁣− – – – – 1 0 259 0 1 0 2
Ubiquity Datrana −⁣− – – – – 9% 0% 27% 0% 9% 0% 18%
Loteshwar
011109/4 34.6 0.9 –28.2 –27.5 0.7 0 1 0 0 0 0 0
021109/2 52.4 1.1 –28.5 –29.0 –0.5 0 0 0 4 0 0 2
021109/8 12.8 1.0 –30.3 –29.6 0.7 0 0 0 0 0 0 0
031109/11 130.6 0.9 –30.1 –29.5 0.6 0 1 0 0 0 0 0
071109/19 60.3 0.9 –30.2 –29.7 0.5 0 0 0 0 0 0 0
071109/20 42.8 1.0 –28.4 –29.1 –0.7 0 0 0 4 0 0 0
Total Loteshwar −⁣− – – – – 0 2 0 8 0 0 2
Ubiquity Loteshwar −⁣− – – – – 0% 33% 0% 33% 0% 0% 17%
Shikarpur
240112/23 99.2 0.9 –29.9 –29.5 0.4 0 0 1 0 0 0 0
240112/29 97.7 0.9 –26.9 –28.5 –1.6 0 0 0 2 0 0 0
250112/8 12.8 1.0 –30.5 –29.8 0.7 0 0 0 0 0 1 0
270112/8 86.1 0.9 –30.2 –29.6 0.6 0 0 0 0 0 0 0
010212/16 339.1 1.0 –26.3 –28.0 –1.8 0 0 0 0 0 0 0
010212/19 79.8 0.9 –26.6 –27.8 –1.2 0 0 0 1 0 0 0
020212/9 52.2 0.9 –29.9 –29.5 0.3 0 0 0 0 0 0 0
050212/12 44.9 0.9 –29.4 –29.3 0.1 0 0 0 1 0 0 0
050212/8 25.4 1.0 –27.2 –29.9 –2.7 0 0 0 7 1 0 0
070212/12 139.1 0.9 –30.3 –29.6 0.7 0 0 0 0 0 0 1
090212/7 135.1 0.7 –26.6 –23.5 3.2 0 0 0 0 0 0 0
Total Shikarpur −⁣− – – – – 0 0 1 11 1 1 1
Ubiquity Shikarpur −⁣− – – – – 0% 0% 9% 36% 9% 9% 9%

The description of the features used to taxonomically identify starch grains can be found in Table 3.

All lipid profiles were dominated by mid-chain (C14:0, C16:0 and C18:0) fatty acids and small peaks of long-chain fatty acids (C20:0–C24:0). All the profiles also contained odd-chain fatty acids such as C15:0 and C17:0, as well as branched-chain fatty acids such as C15 and C17. Such profiles are characteristic of degraded animal fats (Dudd and Evershed, 1998; Dudd et al., 1999; Halmemies-Beauchet-Filleau et al., 2013, 2014; Supplementary Figure 2). The palmitic/stearic (P/S) ratio of the extracts ranged between 0.7 and 1.2, which is also indicative of animal products.

No peaks of n-alkanes, which can be indicative of plant waxes (Kolattukudy, 1970), were detected in the extracts. Similarly, miliacin, the chemical “biomarker” for Panicum spp. and some other small millets (Lu et al., 2009; Bossard et al., 2013), could not be identified in the extracts. Other chemical biomarkers, such as compounds formed during exposure to high temperatures or those indicating heated fish products—mid-chain ketones and long-chain (> C18) ω-(o-alkylphenyl) alkanoic acids and isoprenoid fatty acids, respectively—could also not be identified in the lipid extracts using the methods described.

 Compound-Specific Isotopic Evidence (Gas Chromatography-Combustion-Isotope Ratio Mass Spectrometry)

The fatty acid δ13C-values from the analyzed extracts ranged between –30.5 and –24.5‰ (C16:0) and between –29.9 and –23.5‰ (C18:0), reflecting a relatively narrow spread of values. The Δ13C-values ranged from 3.2 to –2.7‰, which are consistent with modern reference fats from non-ruminant and ruminant adipose fats (Figure 2, Table 2, and Supplementary Figure 3). Most of the extracts (n = 21, 78%) had Δ13C-values which fall within the range of reference non-ruminant fats, while the remaining (n = 6, 22%) fall within the range for ruminant adipose fats, suggesting they were predominantly used to process the carcass fats of cattle/buffalo, sheep/goat or wild ruminant animals such as deer or nilgai (Boselaphus tragocamelus). None of the extracts had values consistent with reference dairy-based products. Some extracts also have fatty acid δ13C-values that are consistent with those reported for plant oils, especially C3 oilseed plants such as sesame (Steele et al., 2010). However, plants have a much higher palmitic/stearic ratio than animal fats, and they produce significant deviations in Δ13C-values depending on the absolute δ13C-values of the end-members (Steele et al., 2010; Hendy et al., 2018). Clear evidence for plant products was not detected in the lipid extracts in this study.

Compound-specific isotopic values by site. The difference in δ13C-values for C16:0 and C18:0 fatty acids (Δ13C) are plotted against δ13C-values for C16:0 fatty acids against a global reference range for modern porcine, ruminant carcass, dairy fats and plant oils (Copley et al., 2003; Spangenberg et al., 2006; Gregg et al., 2009; Outram et al., 2009; Steele et al., 2010; Dunne et al., 2012; Hendy et al., 2018).

The δ13C16:0 values of the extracts indicate the input of both C3 and C4 plants, likely routed through animal diet (Halmemies-Beauchet-Filleau et al., 2013, 2014). A moderate negative correlation coefficient between the δ13C16:0 and Δ13C values was observed (–0.35, p = 0.012), suggesting that samples with more negative Δ13C-values produced fatty acids enriched in 13C, likely from tissues of ruminant animals that were consuming C4 plants. However, there was no effect of site or vessel form on δ13C16:0 values [χ2(2) = 1.3, p = 0.53 and χ2(6) = 4.6, p = 0.59, respectively] or on δ13C18:0 values [χ2(2) = 0.4, p = 0.82 and χ2(6) = 2.3, p = 0.89, respectively]. This indicates that no significant differences in the usage of vessels across sites and vessel forms can be detected (see Supplementary Figure 4).

 Microbotanical Evidence

Starch grains were generally scarce in most samples (0–8 grains per sample) but very abundant in sample DTR116 (n = 256) (Table 2). The presence and relative abundance of starch types varies greatly among sites (Figure 3 and Table 3). Thus, whereas at Datrana most starch grains (99% of the assemblage) belong to the Hordeeae tribe (Poaceae, grasses), these are a minor component of the assemblage at Shikarpur (7%) and completely absent at Loteshwar. Plants from the Hordeeae tribe are not native to Gujarat, suggesting these starch grains most likely represent wheat (Triticum sp.) and/or barley (Hordeum vulgare) traded in from other areas of the Indus Civilization, such as the Indus plain in present-day Sindh (Pakistan). Conversely, starch grains from the Faboideae subfamily (Fabaceae, pulses) predominate at both Shikarpur (73%) and Loteshwar (67%) but are completely absent from Datrana. The starch assemblage further includes one small millet (Paniceae tribe, Poaceae) grain from Datrana, two Job’s tears (Coix lacryma-jobi, Andropogoneae tribe, Poaceae) grains from Loteshwar, two ginger (Zingiber sp., Zingiberaceae) grains (one from Datrana and one from Shikarpur), one grain from an unidentified underground storage organ from Shikarpur and a few grains that could not be identified due to severe damage (n = 1), the lack of comparable modern reference material (n = 1) or the lack of diagnostic features in the starch grains (n = 3). Phytoliths were virtually absent from all but three vessels (Supplementary Table 3) and are therefore not considered in the discussion of the results.

Starch grains observed in samples from pottery vessels from Datrana, Loteshwar and Shikarpur: (A) Paniceae (small millet) grain form sample DTR120; (B) Andropogoneae (Job’s tears) grain from sample 031109/11; (C) cluster of Hordeeae (wheat/barley) grains from sample DTR116; (D) Faboideae (pulses) grains from sample 071109/20; (E) Zingiber sp. (ginger) grain from sample 050212/8; and (F) undetermined underground storage organ from sample 250112/8. Scale bar: 20 μm (for comparative reference material, see Supplementary Figure 1).

Shape and size of the starch grain morphotypes found in pottery vessels from Datrana, Loteshwar and Shikapur.

Taxonomy No Shape Size (μ m)
Mean (SD) Range Nm
Paniceae (small millets) 1 Polyhedral 7,851 – 1
Andropogoneae (Job’s tears) 2 Polyhedral 18,734 (1,699) 17,533–19,936 2
Hordeeae Type A (wheat/barley) 24 Rounded/ovate, discoidal 23,552 (3,193) 17,942–26,666 7
Hordeeae Type B (wheat/barley) 239* Rounded, spherical 6,302 (1,743) 3,779–7,785 4
Faboideae (pulses) 19 Ovate, ovoid 22,797 (8,049) 10,907–36,994 17
Zingiber sp. (ginger) 2 Ovate, discoidal 23,653 (6,339) 19,171–28,136 2
Unidentified underground storage organ 1 Globose Broken, not measured – –
Unknown 1 Ovate, ovoid 14,305 – 1
Damaged unidentified 1 – Broken, not measured – –

*This morphotype occurs in numerous plants and therefore was only assigned to wheat/barley when found in association with Hordeeae Type A starch grains (e.g., Figure 3C). No = number of observed starch grains. Nm = number of measured starch grains.

 Discussion

Direct archeological evidence for food preparation and consumption is often limited to exceptional archeological finds (e.g., Lu et al., 2005). In most cases, the reconstruction of past foodways is hampered by the poor survival of organic remains, taphonomic pathways and depositional and post-depositional processes. In this study, the low preservation of starch grains and survival of only free fatty acids within the lipid extracts limits the extent to which past foodways can be reconstructed. The interpretation of the results from northern Gujarat is further hindered by few previous biomolecular studies and lack of modern reference fats from the region; in fact, this is one of the few lipid residue studies on South Asian archeological material (see also Chakraborty et al., 2020; Suryanarayan et al., 2021) and the only study in South Asian archeology integrating chemical and microbotanical analysis from the same vessels. Nevertheless, a holistic approach to past foodways has the potential to overcome biases by integrating evidence, and opens up new questions to examine the production, distribution, preparation and consumption of both animal and plant resources.

Plant and Animal Ingredients in Vessels?

The obtained lipid residue results suggest the presence of degraded animal fats in the vessels. When compared with available modern reference fats from other parts of the world, the fatty acid isotopic values suggest the dominance of the processing of non-ruminant fats in the vessels, possibly porcine carcass fats or fats of other monogastric animals, such as birds. However, the presence of non-ruminant fats in vessels does not correlate with the faunal record from Chalcolithic/Harappan northern Gujarat, as there is limited evidence for the exploitation of pigs (domestic or wild) or other omnivorous taxa. Pig remains are relatively scarce (0–15% NISP) in Chalcolithic/Harappan settlements in Gujarat compared to the remains of ruminants, particularly sheep/goat (15–40% NISP) and cattle/buffalo (50–70% NISP), which dominate the faunal assemblages (Thomas et al., 1996, 1997; Patel, 2009; Chase, 2010, 2014; Joglekar and Goyal, 2011; Goyal, 2013; Joglekar et al., 2013). Similar results were reported from Harappan sites in northwest India, where porcine remains are similarly scarce in the zooarcheological assemblage (Suryanarayan et al., 2021). While it is possible that non-ruminant animal resources were selectively processed in pots, it is also worth exploring alternative explanations for their predominance in the pottery vessels from northern Gujarat.

Mixing models to investigate how mixtures of plant and animal products processed in ancient vessels can affect compound-specific isotope results have demonstrated that it can be challenging to interpret fatty acids isotopic results when plant and animal products are mixed in vessels, particularly from environments where both C3 and C4 plants are available (Hendy et al., 2018; Bondetti et al., 2020; Suryanarayan, 2020; Suryanarayan et al., 2021). Hendy et al. (2018) showed that Δ13C-values similar to non-ruminant fats can be created when vessels are used to process mixtures of ruminant adipose products and C3 plants. Considering the available zooarcheological evidence and the predominance of C3 plants in the starch assemblage from pottery vessels from Datrana, Loteshwar and Shikarpur, it is possible that the Δ13C-values of the extracts resulted from mixtures of plant and animal products rather than non-ruminant animal fats. Although the low lipid contribution of plants may not be enough to influence the δ13C-value of the animal-derived lipids (Miller et al., 2020), a likely possibility could be the use of C3 plant oils, such as sesame, which has been found at Loteshwar (García-Granero et al., 2016), which have higher lipid content. These hypotheses, however, need to be tested through additional mixing models and experimental research.

A small number of vessels from Datrana and Shikarpur have extracts with isotopic values that fall within the range of ruminant carcass fats and have δ13C16:0 values between –27 and –24‰. These values suggest that fats processed in these vessels were from animals with a mixed C3–C4 plant diet (e.g., Roffet-Salque et al., 2016). The zooarcheological record in prehistoric northern Gujarat includes both domestic (cattle, sheep and goats) and wild ruminants (e.g., nilgai). The stable carbon isotope signature of animals grazing on wild vegetation would reflect the consumption of mostly C3 plants, which form the majority of the local flora (Chase et al., 2014: 4). On the other hand, the stable carbon isotope signature of animals grazing on agricultural fodder would reflect the consumption of both C4 (e.g., small millets) and C3 plants (e.g., tropical pulses), as shown by stable carbon isotope analyses of tooth enamel of domestic ruminants at Bagasra (Chase et al., 2014, 2018), Shikarpur (Chase et al., 2020), Jaidak (Chase et al., 2020) and Kotada Bhadli (Chakraborty et al., 2018). Therefore, the presence of C4 plants in the diet of the ruminants processed at Datrana and Shikarpur seems to indicate that they were partly fed agricultural by-products, and thus the molecular evidence likely represents the processing of domestic animals in vessels. Moreover, the study of burnt animal bones from Kanmer suggested that, unlike domestic taxa, wild mammals such as deer were regularly prepared via roasting (Goyal, 2017), further suggesting that the ruminant fats from Datrana and Shikarpur derive from domestic animals.

 Prehistoric Foodways in Northern Gujarat

Archeobotanical and zooarcheological research has shown that the inhabitants of Chalcolithic/Harappan northern Gujarat consumed a wide array of plant and animal resources. Lipid residue and starch grain analyses, however, only identified a few of these resources in the pottery vessels from Datrana, Loteshwar and Shikarpur. Although we acknowledge that the small set of samples analyzed may not be generalizable, the results open up questions about the culinary pathways of certain foodstuffs.

Among the food ingredients potentially consumed by prehistoric populations in northern Gujarat, two were not detected in the extracts: fish and dairy products. Molecular markers for fish products were not detected in the analyzed extracts, and the compound-specific isotopic results do not suggest aquatic input. In contrast, aquatic products were tentatively identified from a single vessel at Kotada Bhadli (Chakraborty et al., 2020). Although “absence of evidence is not evidence of absence”—and the techniques required to detect fish products in lipid extracts may need to be more sensitive in nature—fish might have had a specific culinary pathway at the sites examined in this study. Fish (both marine and riverine) seem to have been an important component of the Harappan diet, but how this resource was part of the Harappan cuisine has rarely been explored. Ichthyoarcheological studies on the materials from several Indus Civilization sites highlighted the paucity of cut marks on the recovered bones (< 2% of the assemblages in all studied sites; Abhayan, 2016; Belcher, 1998). At Bagasra, Kanmer and Shikarpur, most cut marks were noticed on the vertebral elements, which probably indicates standardization in chopping (Abhayan, 2016; Abhayan et al., 2018) and suggests that fish might have been cooked mostly whole and not dismembered. Visible signs of heat alteration on some of the fish bones (Abhayan, 2016, p. 298) is also commensurate with roasting or frying rather than cooking in vessels, potentially explaining why fish products were not detected in the vessels analyzed in this study.

None of the fatty acid δ13C-values from the lipid extracts were consistent with modern references of dairy products, including an ethnographic dairy vessel from Gujarat (Craig et al., 2005). Our results are similar to those obtained from northwest India, where only a small percentage of vessels had evidence of dairy products (Suryanarayan et al., 2021), but contrast with lipid residue analysis of vessels from Kotada Bhadli, where 6 out 21 vessels (mostly bowls) had evidence of dairy use (Chakraborty et al., 2020). Although the origins of dairying in prehistoric South Asia are still unknown, it is assumed that secondary products such as dairy were an important part of the economy, especially during the Harappan period (Gouin, 1990; Bourgeois and Gouin, 1995; Miller, 2004; Wright, 2010; Chase et al., 2014). Once again, “absence of evidence is not evidence of absence”; however, a number of possibilities must be explored. It is possible that the Chalcolithic semi-nomadic agro-pastoralists at Loteshwar and inhabitants of Datrana were not engaging in dairying; however, this is harder to explain for Shikarpur, where zooarcheological analysis suggests that cattle and buffalo were probably kept for secondary products and/or animal traction (Chase, 2010, 2014). It may also be possible that at Datrana, Loteshwar and Shikarpur milk and dairy products were processed in containers, such as bowls, that were not part of this study, or in vessels that have not survived in the archeological record (e.g., wooden vessels). In present-day traditional northern Gujarat milk is collected in small or medium pots with a slightly wide, open mouth and consumed in bowls or cups (P. Ajithprasad personal observation), and milk collection vessels are seldom used for any other purpose. Thus, it is possible that their representation in archeological contexts might also be limited. Finally, it is also possible that the mixing of resources in vessels has made it harder to detect the direct processing of dairy products in vessels.

By comparing the starch evidence from the ceramic vessels with the evidence for the use of starch-rich plant parts from previously published macro and microbotanical analyses at Datrana, Loteshwar and Shikarpur we can attempt to disentangle the different culinary pathways potentially followed by each plant resource (Figure 4 and Table 4). The starch evidence from potsherds was generally scarce and therefore any interpretation using this proxy must be considered with caution. Only 10% of the starch residue was analyzed from each sample to allow comparability with previous starch analyses on grinding stones, which may have contributed to the low number of starch grains recovered in this study—which suggests that grinding stones might provide a more suitable environment for the preservation of starch grains, although such a discussion is beyond the scope of the present study. In any case, the marked taxonomic differences observed between the assemblages from grinding stones and pottery vessels in all studied sites suggests culinary choices and the different pathways in which plant ingredients were processed and prepared. At Datrana, macrobotanical remains and phytoliths were extremely scarce. However, the starch evidence from grinding stones suggested the inhabitants of this lithic blade workshop consumed mainly small millets and Job’s tears, complemented by small amounts of pulses, wheat/barley and ginger (García-Granero et al., 2017a). In striking contrast, small millet and Job’s tears starch grains were virtually absent from pottery vessels, which mostly included wheat/barley and a single ginger starch grain. At Loteshwar and Shikarpur, both the macrobotanical assemblage and the microbotanical evidence from grinding stones suggested that locally cultivated small millets and, to a lesser degree, Job’s tears and pulses formed the basis of the diet of the inhabitants of these settlements, complemented by wheat and barley traded in from other regions within the greater Indus Valley more prone to the cultivation of winter cereals (García-Granero et al., 2015, 2016). The starch evidence from pottery vessels, on the other hand, was mostly composed of pulses, with a minor presence of Job’s tears at Loteshwar and underground storage organs (including ginger) and wheat/barley at Shikarpur.

Plant foodways at Datrana, Loteshwar and Shikarpur. Comparison of the main groups of starch-rich plants identified in archeological sediments (acquired resources), grinding stones (ground resources) and pottery vessels (cooked resources). Very few macrobotanical remains were found at Datrana (García-Granero et al., 2017a) and therefore have not been included in this figure. Raw data can be found in Supplementary Table 2. Panicoideae: several small millet species and Job’s tears; Hordeeae: wheat and barley; Faboideae: several pulses, including horsegram and Vigna sp.; Underground storage organs: ginger and other unidentified resources.

Plant foodways at Datrana, Loteshwar and Shikarpur.

Acquired resources Ground resources Cooked resources
Datrana
Cultivated: – Gathered: underground storage organs (USOs) Traded: small millets, Job’s tears, barley (and wheat?), pulses and ginger Small millets, Job’s tears, wheat/barley, pulses and USOs (incl. ginger) Small millets, wheat/barley and USOs (incl. ginger)
Loteshwar
Cultivated: browntop millet, Setaria spp., other small millets, Job’s tears, horsegram and sesame Gathered: Egyptian crowfoot grass and USOs Traded: wheat (and barley?) Small millets, Job’s tears, wheat/barley, Egyptian crowfoot grass, pulses and USOs Job’s tears and pulses
Shikarpur
Cultivated: browntop millet, Setaria spp., Job’s tears and Vigna sp. Gathered: USOs Traded: barley (and wheat?), rice and ginger Small millets, Job’s tears, wheat/barley, pulses and USOs Wheat/barley, pulses and USOs (incl. ginger)

Comparison of macroscopic remains of cultivated, gathered and traded edible plants most commonly recovered from archeological sediments (acquired resources), starch grains recovered from grinding stones (ground resources) and starch grains recovered from pottery vessels (cooked resources). Most common plant resources in each category are highlighted in bold.

The evidence from the analyzed pottery vessels thus suggests that plant resources were not only acquired in different ways (cultivated, gathered or traded in) but also followed different culinary pathways. Small millets were probably ground to prepare flour-based products. Job’s tears seems to have been consumed in a similar way, possibly mixed with small millet flour. Miliacin was not detected in the lipid extracts (although perhaps more sensitive techniques are required), which might suggest millets were not cooked in vessels. Starch grains of pulses, on the other hand, are detected in both grinding stones and from pottery vessels, which might indicate they were incorporated in a wider range of meals using different processing techniques. Grinding stones may have been used both for grinding pulses into flour or used for preparing other dishes—and for splitting the seeds to form the basis of stews or soups. Starches of wheat and barley are not common in the studied sites but appear in abundance in vessel DTR116 from Datrana. Other types of starch grains were not found in this vessel, and it had a very low lipid concentration (below 5 μg/g). The detection of cereals via lipid residue analysis remains challenging, since biomarkers such as plant sterols and alkylresorcinols are highly susceptible to degradation, and their uptake into the ceramic matrix is very limited (Hammann and Cramp, 2018). All the available evidence, thus, seems to suggest that vessel DTR116 was used exclusively to prepare wheat/barley-based foods. Finally, ginger may have been ground on grinding stones and also incorporated during cooking. Based on a qualitative assessment of our reference collection, taxa within the ginger family produce notably less starch grains than cereals and pulses. Therefore, the presence of ginger at Datrana and Shikarpur, though minor, attests to its culinary use at these sites. Ginger starch grains were also documented in human dental calculus from Harappan Farmana, in Haryana (northern India) (Weber and Kashyap, 2010; Kashyap and Weber, 2013), highlighting that their use as food condiments may have been widespread across the greater Indus Valley during the Chalcolithic/Harappan period.

 Conclusion

The combination of methods used in this study provides a unique means to explore the culinary use of both plant and animal ingredients that may not be detectable via conventional techniques in archeobotany and zooarcheology. They also provide a way to overcome interpretational limitations posed by individual methods, especially for the detection of plant remains within vessels. This preliminary study suggests the dominance of animal fats in vessels, although the interpretation of the compound-specific results is challenging, and pottery vessels may have been used to process C3 plants/oils and fats from ruminants, which is supported by microbotanical analyses. The presence of Δ13C-values similar to non-ruminant fats or other omnivorous taxa might have resulted from the combination of these ingredients, either as part of a single dish combining plant and animal foodstuffs or as a result of multiple cooking events. No differences in culinary practices can be detected across the studied sites, which had markedly different ceramic traditions and relied on different subsistence strategies. In particular, the predominance of non-ruminant fats and/or admixtures of plants and ruminant fats in pottery vessels and the use of small millets for producing flour have been observed in all studied sites. Overall, the interpretations offered are, of course, tentative given the sample size and issues of differential preservation of the molecular and microbotanical evidence but nonetheless illustrative of how finer grained interpretations may be gleaned from multi-proxy approaches of this nature. Future research in this and neighboring areas will assess whether the observed culinary continuity (spanning over a thousand years) is representative of deeper cultural practices common to prehistoric populations in northern Gujarat and other Indus regions.

 Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

 Author Contributions

PA and MM acquired funding, provided samples, and designed the study. JG-G, AS, MCu, OC, and MCá analyzed the data. JG-G, AS, MCu, OC, and MM interpreted the data. JG-G and AS drafted the main manuscript and prepared the figures. All authors reviewed the manuscript.

 Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

 Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Funding

Fieldwork at Datrana and Loteshwar was funded by the Spanish Ministry of Education, Culture and Sport (Programa de Ayudas para Proyectos Arqueológicos en el Exterior 2009 and 2010) and the Catalan Government (Programa EXCAVA 2009), and fieldwork at Shikarpur was funded by the Department of Archaeology and Ancient History of The Maharaja Sayajirao University of Baroda (Gujarat, India), the Archaeological Survey of India and the Gujarat State Department of Archaeology and Museums. Laboratory work was funded by the Spanish Ministry of Economy and Competitiveness (Programa I + D, HAR2010-16052) and the Spanish Ministry of Science and Innovation (Programa CONSOLIDER INGENIO 2010, CSD2010–00034). JG-G acknowledges funding from the European Commission (H2020 MSCA 2015, Grant No. 704867) and the Spanish Ministry of Science, Innovation and Universities (IJC2018-035161-I). AS thanks the Cambridge Trust, Nehru Trust for Cambridge University and Sidney Sussex College, Cambridge, for providing funding for research and travel. MCu was funded by a Ramón y Cajal fellowship (RYC2019-026697-I) of the Spanish Ministry of Science and Innovation.

We are grateful to the members of the NoGAP and the Department of Archaeology and Ancient History of The Maharaja Sayajirao University of Baroda, India for their help during the field and laboratory work in India, and to the Archaeological Survey of India and the Gujarat State Department of Archaeology and Museums for granting excavation permissions.

 Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2022.840199/full#supplementary-material

 References Abhayan G. S. (2016). Preparation and Application of Fish Bone Identification Keys for Harappan Sites in Gujarat with Special Reference to Bagasra, Kanmer and Shikarpur. Ph.D. thesis. Gurgaon: Deemed University. Abhayan G. S. Joglekar P. P. Ajithprasad P. Krishnan K. Bhanand K. K. Rajesh S. V. (2018). “Exploitation during the Harappan Period at Bagasra in Gujarat, India in An Ichthyoarchaeological Approach,” in Walking with the Unicorn: Social Organization and Material Culture in Ancient South Asia. Jonathan Mark Kenoyer Felicitation, eds Frenez D. Jamison G. M. Law R. W. (Oxford: Archaeopress Archaeology), 1–18. 10.2307/j.ctv19vbgkc.6 Abhayan G. S. Joglekar P. P. Rajesh S. V. Gopan A. S. Chase B. (2016). Fish Otoliths from Navinal, Kachchh, Gujarat: identification of Taxa and Its Implications. Herit. J. Multidiscip. Stud. Archaeol. 4 218–227. Ajithprasad P. (2002). “The pre-Harappan cultures of Gujarat,” in Indian archaeology in retrospect, vol II. Protohistory: archaeology of the Harappan Civilization, eds Settar S. Korisettar R. (New Delhi: Manohar), 129–158. Ajithprasad P. (2011). “Chalcolithic Cultural Patterns and the Early Harappan Interaction in North Gujarat,” in Cultural Relations between the Indus and the Iranian Plateau during the Third Millennium BCE, eds Osada T. Witzel M. (Cambridge: Harvard University), 11–40. Ajithprasad P. Sonawane V. H. (2011). “The Harappa culture in North Gujarat: A regional paradigm,” in Linguistics, Archaeology and Human Past, eds Osada T. Endo H. (Kyoto: Research Institute for Humanity and Nature), 223–269. Belcher W. B. (1998). Fish exploitation of the Baluchistan and Indus Valley traditions: an ethnoarchaeological approach to the study of fish remains. Ph.D. thesis. Madison: University of Wisconsin-Madison. Bhan K. Ajithprasad P. (2009). Excavations at Shikarpur, Gujarat 2008-2009. Available online at: https://www.harappa.com/sites/default/files/pdf/Shikarpur-2008-2009.pdf (accessed February 28, 2022). Bisht R. S. (2015). Excavations at Dholavira (1989-90 to 2004-05). Available online at: https://www.scribd.com/document/262316120/Excavations-at-Dholavifra-1989-2005-RS-Bisht-2015 (accessed February 28, 2022). Bogaard A. Charles M. Twiss K. C. Fairbairn A. Yalman N. (2009). Private pantries and celebrated surplus: storing and sharing food at Neolithic Çatalhöyük. Central Anatolia. Antiquity 83 649–668. 10.1017/S0003598X00098896 Bolger D. (2010). The dynamics of gender in early agricultural societies of the Near East. Signs J. Women Cult. Soc. 35 503–531. 10.1086/605512 Bondetti M. Lucquin A. Savel’ev N. A. Weber A. W. Craig O. E. (2020). Resource processing, early pottery and the emergence of Kitoi culture in CisBaikal: insights from lipid residue analysis of an Early Neolithic ceramic assemblage from the Gorelyi Les habitation site, Eastern Siberia. Archaeol. Res. Asia 24:100225. 10.1016/j.ara.2020.100225 Bossard N. Milbeau C. L. Sauze J. Terwilliger V. (2013). Distribution of miliacin (olean-18-en-3β-ol methyl ether) and related compounds in broomcorn millet (Panicum miliaceum) and other reputed sources: implications for the use of sedimentary miliacin as a tracer of millet. Org. Geochem. 63 48–55. 10.1016/j.orggeochem.2013.07.012 Bourgeois G. Gouin P. R. (1995). ésultats d’une analyse de traces organiques fossiles dans une «faisselle» harappéenne. Paléorient 21 125–128. 10.3406/paleo.1995.4613 Bray T. L. (ed.) (2003). The Archaeology and Politics of Food and Feasting in Early States and Empires. London: Kluwer Academic/Plenum Publishers. 10.1007/b100538 Chakraborty K. S. Chakraborty S. Le Roux P. Miller H. M. L. Shirvalkar P. Rawat Y. (2018). Enamel isotopic data from the domesticated animals at Kotada Bhadli, Gujarat, reveals specialized animal husbandry during the Indus Civilization. J. Archaeol. Sci. Rep. 21 183–199. 10.1016/j.jasrep.2018.06.031 Chakraborty K. S. Slater G. F. Miller H. M. L. Shirvalkar P. Rawat Y. (2020). Compound specific isotope analysis of lipid residues provides the earliest direct evidence of dairy product processing in South Asia. Sci. Rep. 10:16095. 10.1038/s41598-020-72963-y 32999318 Charters S. Evershed R. P. Quye A. Blinkhorn P. W. Reeves V. (1997). Simulation Experiments for Determining the Use of Ancient Pottery Vessels: the Behaviour of Epicuticular Leaf Wax During Boiling of a Leafy Vegetable. J. Archaeol. Sci. 24 1–7. 10.1006/jasc.1995.0091 Charters S. Evershed R. P. Blinkhorn P. W. Denham V. (1995). Evidence for the Mixing of Fats and Waxes in Archaeological Ceramics. Archaeometry 37 113–127. Charters S. Evershed R. P. Goad L. J. Leyden A. Blinkhorn P. W. Denham V. (1993). Quantification and distribution of lipid in archaeological ceramics: implications for sampling potsherds for organic residue analysis and the classification of vessel use. Archaeometry 35 211–223. Chase B. (2010). Social change at the Harappan settlement of Gola Dhoro: a reading from animal bones. Antiquity 84 528–543. 10.1017/s0003598x00066758 Chase B. (2012). “Crafting Harappan Cuisine on the Saurashtran Frontier of the Indus Civilization,” in The Menial Art of Cooking: Archaeological Studies of Cooking and Food Preparation, eds Graff S. Rodriguez-Alegria E. (Boulder: The University Press of Colorado), 145-172. Chase B. (2014). On the pastoral economies of Harappan Gujarat: faunal analyses at Shikarpur in context. Herit. J. Multidiscip. Stud. Archaeol. 2 1–22. Chase B. Meiggs D. Ajithprasad P. (2020). Pastoralism, climate change, and the transformation of the Indus Civilization in Gujarat: faunal analyses and biogenic isotopes. J. Anthropol. Archaeol. 59:101173. 10.1016/j.jaa.2020.101173 Chase B. Meiggs D. Ajithprasad P. Slater P. A. (2014). Pastoral land-use of the Indus Civilization in Gujarat: faunal analyses and biogenic isotopes at Bagasra. J. Archaeol. Sci. 50 1–15. 10.1016/j.jas.2014.06.013 Chase B. Meiggs D. Ajithprasad P. Slater P. A. (2018). What is left behind: advancing interpretation of pastoral land-use in Harappan Gujarat using herbivore dung to examine biosphere strontium isotope (87Sr/86Sr) variation. J. Archaeol. Sci. 92 1–12. 10.1016/j.jas.2018.01.007 Copley M. S. Berstan R. Dudd S. N. Docherty G. Mukherjee A. J. Straker V. (2003). Direct chemical evidence for widespread dairying in prehistoric Britain. Proc. Natl. Acad. Sci. U. S. A. 100 1524–1529. 10.1073/pnas.0335955100 12574520 Copley M. S. Bland H. A. Rose P. Horton M. Evershed R. P. (2005). Gas chromatographic, mass spectrometric and stable carbon isotopic investigations of organic residues of plant oils and animal fats employed as illuminants in archaeological lamps from Egypt. Analyst 130 860–871. 10.1039/b500403a 15912234 Copley M. S. Rose P. J. Clapham A. Edwards D. N. Horton M. C. Evershed R. P. (2001). Detection of Palm Fruit Lipids in Archaeological Pottery from Qasr Ibrim, Egyptian Nubia. Proc. R. Soc. Lond. B Biol. Sci. 268 593–597. 10.1098/rspb.2000.1394 11297176 Correa-Ascencio M. Evershed R. P. (2014). High throughput screening of organic residues in archaeological potsherds using direct acidifed methanol extraction. Anal. Methods 6 1330–1340. 10.1039/c3ay41678j Craig O. E. Chapman J. Heron C. Willis L. Bartosiewicz L. Taylor G. (2005). Did the first farmers of central and eastern Europe produce dairy foods?. Antiquity 79 882–894. 10.1017/S0003598X00115017 Craig O. E. Saul H. Lucquin A. Nishida Y. Taché K. Clarke L. (2013). Earliest evidence for the use of pottery. Nature 496 351–354. 10.1038/nature12109 23575637 Craig O. E. Shillito L. M. Albarella U. Viner-Daniels S. Chan B. Cleal R. (2015). Feeding Stonehenge: cuisine and consumption at the Late Neolithic site of Durrington Walls. Antiquity 89 1096–1109. 10.15184/aqy.2015.110 Crowther A. Haslam M. Oakden N. Walde D. Mercader J. (2014). Documenting contamination in ancient starch laboratories. J. Archaeol. Sci. 49 90–104. 10.1016/j.jas.2014.04.023 Dudd S. N. Evershed R. P. (1998). Direct demonstration of milk as an element of archaeological economies. Science 282 1478–1481. 10.1126/science.282.5393.1478 9822376 Dudd S. N. Evershed R. P. Gibson A. M. (1999). Evidence for varying patterns of exploitation of animal products in different prehistoric pottery traditions based on lipids preserved in surface and absorbed residues. J. Archaeol. Sci. 26 1473–1482. 10.1006/jasc.1998.0434 Dunne J. Evershed R. P. Salque M. Cramp L. Bruni S. Ryan K. (2012). First dairying in green Saharan Africa in the fifth millennium BC. Nat. Lett. 486 390–394. 10.1038/nature11186 22722200 Dunne J. Jórdeczka M. Chłodnicki M. Hardy K. Kubiak-Martens L. Hoyo M.-M.-D. (2021). Holocene resource exploitation along the Nile: diet and subsistence strategies of Mesolithic and Neolithic societies at Khor Shambat 1, Sudan. Antiquity 95 1426–1445. 10.15184/aqy.2021.141 Dunne J. Mercuri A. Evershed R. Bruni S. di Lernia S. (2016). Earliest direct evidence of plant processing in prehistoric Saharan pottery. Nat. Plants 3:16194. 10.1038/nplants.2016.194 27991880 Evershed R. P. Payne S. Sherratt A. Coolidge J. Urem-Kotsu D. (2008). Earliest date for milk use in the Near East and southeastern Europe linked to cattle herding. Nature 455 528–531. 10.1038/nature07180 18690215 Farooqui A. Gaur A. S. Prasad V. (2013). Climate, vegetation and ecology during Harappan period: excavations at Kanjetar and Kaj, mid-Saurashtra coast, Gujarat. J. Archaeol. Sci. 40 2631–2647. 10.1016/j.jas.2013.02.005 Fuller D. Q. (2011). Finding plant domestication in the Indian subcontinent. Curr. Anthropol. 52 S347–S362. 10.1086/658900 Gaastra J. S. de Vareilles A. Vander Linden M. (2019). Bones and seeds: an integrated approach to understanding the spread of farming across the Western Balkans. Environ. Archaeol. 27 44–60. 10.1080/14614103.2019.1578016 Gadekar C. Ajithprasad P. Madella M. (2013). Crested ridge technique and lithic assemblage from Datrana, Gujarat. Herit. J. Multidiscip. Stud. Archeol. 1 16–28. Gadekar C. Ajithprasad P. Madella M. Balbo A. Rajesh S. V. Rondelli B. (2014a). Continuation of a tradition over five thousand years: lithic assemblage from Loteshwar, North Gujarat, western India. Herit. J. Multidiscip. Stud. Archaeol. 2 283–304. Gadekar C. Rajesh S. V. Ajithprasad P. (2014b). Shikarpur lithic assemblage: new questions regarding Rohri chert blade production. J. Lithic Stud. 1, 137–149. García-Granero J. J. (2015). From gathering to farming in semi-arid northern Gujarat (India): a multi-proxy approach. Ph.D. thesis. Barcelona: University of Barcelona. Available online at: https://www.tdx.cat/handle/10803/351960#page=1 García-Granero J. J. Gadekar C. Esteban I. Lancelotti C. Madella M. (2017a). What is on the craftsmen’s menu? Plant consumption at Datrana, a 5000-year-old lithic blade workshop in North Gujarat, India. Archaeol. Anthropol. Sci. 9 251–263. 10.1007/s12520-015-0281-0 García-Granero J. J. Lancelotti C. Madella M. (2017b). A methodological approach to the study of microbotanical remains from grinding stones: a case study in northern Gujarat (India). Veget. Hist. Archaeobot. 26 43–57. 10.1007/s00334-016-0557-z García-Granero J. J. Lancelotti C. Madella M. (2015). A tale of multi-proxies: integrating macro-and microbotanical remains to understand subsistence strategies. Veget. Hist. Archaeobot. 24 121–133. 10.1007/s00334-014-0486-7 García-Granero J. J. Lancelotti C. Madella M. Ajithprasad P. (2016). Millets and herders: the origins of plant cultivation in semiarid north Gujarat (India). Curr. Anthropol. 57 149–163. 10.1086/685775 García-Granero J. J. Urem-Kotsou D. Bogaard A. Kotsos S. (2018). Cooking plant foods in the northern Aegean: microbotanical evidence from Neolithic Stavroupoli (Thessaloniki, Greece). Quat. Int. 496 140–151. 10.1016/j.quaint.2017.04.007 Gouin P. (1990). Rapes, jarres et faisselles: la production et l’exportation des produits laitiers dans l’Indus du 3e millénaire. Paléorient 16 37–54. 10.3406/paleo.1990.4531 Goyal P. (2013). Multiple Role of Cattle in the Harappan Economy at Kanmer Gujarat. Herit. J. Multidiscip. Stud. Archeol. 1 63–77. Goyal P. (2017). The preparation and discard of animal products at Kanmer, Gujarat. Shodha Patrika 68 52–61. Gregg M. W. Banning E. B. Gibbs K. Slater G. F. (2009). Subsistence practices and pottery use in Neolithic Jordan: molecular and isotopic evidence. J. Archaeol. Sci. 36 937–946. 10.1016/j.jas.2008.09.009 Halmemies-Beauchet-Filleau A. Vanhatalo A. Toivonen V. Heikkilä T. Lee M. R. Shingfield K. J. (2013). Effect of replacing grass silage with red clover silage on ruminal lipid metabolism in lactating cows fed diets containing a 60:40 forage-to-concentrate ratio. J. Dairy Sci. 96 5882–5900. 10.3168/jds.2013-6872 23849641 Halmemies-Beauchet-Filleau A. Vanhatalo A. Toivonen V. Heikkilä T. Lee M. R. Shingfield K. J. (2014). Effect of replacing grass silage with red clover silage on nutrient digestion, nitrogen metabolism, and milk fat composition in lactating cows fed diets containing a 60:40 forage-to-concentrate ratio. J. Dairy Sci. 97 3761–3776. 10.3168/jds.2013-7358 24679932 Hammann S. Cramp L. J. E. (2018). Towards the detection of dietary cereal processing through absorbed lipid biomarkers in archaeological pottery. J. Archaeol. Sci. 93 74–81. 10.1016/j.jas.2018.02.017 Hart T. C. (2011). Evaluating the usefulness of phytoliths and starch grains found on survey artifacts. J. Archaeol. Sci. 38 3244–3253. 10.1016/j.jas.2011.06.034 Hendy J. Colonese A. C. Franz I. Fernandes R. Fischer R. Orton D. (2018). Ancient proteins from ceramic vessels at Çatalhöyük West reveal the hidden cuisine of early farmers. Nat. Commun. 9:4064. 10.1038/s41467-018-06335-6 30283003 Heron C. Shoda S. Breu Barcons A. Czebreszuk J. Eley Y. (2016). First molecular and isotopic evidence of millet processing in prehistoric pottery vessels. Sci. Rep. 6:38767. 10.1038/srep38767 28004742 Högberg A. Puseman K. Yost C. (2009). Integration of use-wear with protein residue analysis – a study of tool use and function in the south Scandinavian Early Neolithic. J. Archaeol. Sci. 36 1725–1737. 10.1016/j.jas.2009.03.030 Ivanova M. De Cupere B. Ethier J. Marinova E. (2018). Pioneer farming in southeast Europe during the early sixth millennium BC: climate-related adaptations in the exploitation of plants and animals. PLoS One 13:e0197225. 10.1371/journal.pone.0197225 29775469 Joglekar P. Goyal P. (2011). Faunal remains from Shikarpur, a Harappan site in Gujarat, India. Iran. J. Archaeol. Stud. 1 15–25. Joglekar P. P. Abhayan G. S. Mungur-Medhi J. Nath J. A. (2013). Glimpse of Animal Utilization Pattern at Khirsara, Kachchh, Gujarat, India. Herit. J. Multidiscip. Stud. Archeol. 1 1–15. Kansa S. W. Kennedy A. Campbell S. Carter E. (2009). Resource exploitation at late Neolithic Domuztepe: faunal and botanical evidence. Curr. Anthropol. 50 897–914. 10.1086/605910 Kashyap A. Weber S. (2013). “Starch Grain Analysis and Experiments Provide Insights into Harappan Cooking Practices,” in Connections and Complexity: New Approaches to the Archaeology of South Asia, eds Abraham S. A. Gullapalli P. Raczek T. P. Rizvi U. Z. (California: Left Coast Press), 177–194. Kirch P. O’Day S. J. (2003). New archaeological insights into food and status: a case study from pre-contact Hawaii. World Archaeol. 34 484–497. 10.1080/0043824021000026468 Kolattukudy P. (1970). Plant Waxes. Lipids 5 9–275. 10.1007/BF02532477 Kooiman S. M. Albert R. K. Malainey M. E. (2021). Multiproxy Analysis of Adhered and Absorbed Food Residues Associated with Pottery. J. Archaeol. Method Theory 10.1007/s10816-021-09537-3 [Epub ahead of print]. Lu H. Yang X. Ye M. Liu K. B. Xia Z. Ren X. (2005). Culinary archaeology: millet noodles in late Neolithic China. Nature 437:967. 10.1038/437967a 16222289 Lu H. Zhang J. Liu K. B. Wu N. Li Y. Zhou K. (2009). Earliest domestication of common millet (Panicum miliaceum) in East Asia extended to 10,000 years ago. Proc. Natl. Acad. Sci. U. S. A. 106 7367–7372. 10.1073/pnas.0900158106 19383791 McClatchie M. Schulting R. McLaughlin R. Colledge S. Boggard A. Barratt P. (2019). Food Production, Processing and Foodways in Neolithic Ireland. Environ. Archaeol. 27 80–92. 10.1080/14614103.2019.1615215 Miller L. (2004). Urban economies in urban states: the secondary products revolution in the Indus Civilisation. Ph.D. thesis. New York: New York University. Miller M. Whelton H. L. Swift J. A. Maline S. Hammann S. Cramp L. J. E. (2020). Interpreting ancient food practices: stable isotope and molecular analyses of visible and absorbed residues from a year-long cooking experiment. Sci. Rep. 10:13704. 10.1038/s41598-020-70109-8 32855436 Outram A. K. Stear N. A. Bendrey R. Olsen S. Kasparov A. Zaibert V. (2009). The earliest horse harnessing and milking. Science 323 1332–1335. 10.1126/science.1168594 19265018 Patel A. K. (2009). Occupational histories, settlements, and subsistence in western India: what bones and genes can tell us about the origins and spread of pastoralism. Anthropozoologica 44 173–188. 10.5252/az2009n1a8 33311142 Peres T. M. (2017). Foodways archaeology: a decade of research from the southeastern United States. J. Archaeol. Res. 25 421–460. 10.1007/s10814-017-9104-4 Pokharia A. K. Agnihotri R. Sharma S. Bajpai S. Nath J. Kumaran R. N. (2017). Altered cropping pattern and cultural continuation with declined prosperity following abrupt and extreme arid event at~ 4,200 yrs BP: evidence from an Indus archaeological site Khirsara, Gujarat, western India. PLoS One 12:e0185684. 10.1371/journal.pone.0185684 28985232 Pokharia A. K. Kharakwal J. S. Rawat R. S. Osada T. Nautiyal C. M. (2011). Archaeobotany and archaeology at Kanmer, a Harappan site in Kachchh, Gujarat: evidence for adaptation in response to climatic variability. Curr. Sci. 100 1833–1846. Rajesh S. V. Gadekar C. Abhayan G. S. Ajithprasad P. Sharma B. P. (2018). Janan: a Pre-Urban Harappan site on Khadir Island, Kachchh District, Gujarat. Puratattva 48 173–183. Rajesh S. V. Krishnan K. Ajithprasad P. Sonawane V. H. (2013a). Evaluating the Anarta tradition in the light of material culture from Loteshwar and other sites in Gujarat. Man Environ. 38 10–45. Rajesh S. V. Krishnan K. Ajithprasad P. Madella M. (2013b). Pre-Prabhas assemblage from Gujarat, western India. Herit. J. Multidiscip. Stud. Archaeol. 1 181–209. Reber E. A. Evershed R. P. (2004). How did Mississippians prepare maize? The application of compound-specific carbon isotope analysis to absorbed pottery residues from several Mississippi Valley sites. Archaeometry 46 19–33. 10.1111/j.1475-4754.2004.00141.x Roffet-Salque M. Lee M. R. F. Timpson A. Evershed R. P. (2016). Impact of modern cattle feeding practices on milk fatty acid stable carbon isotope compositions emphasise the need for caution in selecting reference animal tissues and products for archaeological investigations. Archaeol. Anthropol. Sci. 9 1343–1348. 10.1007/s12520-016-0357-5 Roffet-Salque M. Regert M. Evershed R. P. Outram A. K. Cramp L. J. Decavallas O. (2015). Widespread exploitation of the honeybee by early Neolithic farmers. Nature 527 226–230. 10.1038/nature15757 26560301 Salque M. Radi G. Tagliacozzo A. Uria B. P. Wolfram S. Hohle I. (2012). New insights into the Early Neolithic economy and management of animals in Southern and Central Europe revealed using lipid residue analyses of pottery vessels. Anthropozoologica 47 45–63. 10.5252/az2012n2a4 33311142 Saul H. Madella M. Fischer A. Glykou A. Hartz S. Craig O. E. (2013). Phytoliths in pottery reveal the use of spice in European prehistoric cuisine. PLoS One 8:e70583. 10.1371/journal.pone.0070583 23990910 Spangenberg J. E. Jacomet S. Schibler J. (2006). Chemical analyses of organic residues in archaeological pottery from Arbon Bleiche 3, Switzerland – evidence for dairying in the late Neolithic. J. Archaeol. Sci. 33 1–13. 10.1016/j.jas.2005.05.013 Spielmann K. A. (2002). Feasting, craft specialisation, and the ritual mode of production in small-scale societies. Am. Anthropol. 104 195–207. 10.1525/aa.2002.104.1.195 Staller J. Carrasco M. (eds) (2009). Pre-Columbian Foodways: Interdisciplinary Approaches to Food, Culture, and Markets in Ancient Mesoamerica. Berlin: Springer Science & Business Media. 10.1007/978-1-4419-0471-3 Steele V. Stern B. Stott A. W. (2010). Olive oil or lard? Distinguishing plant oils from animal fats in the archeological record of the eastern Mediterranean using gas chromatography/combustion/isotope ratio mass spectrometry. Rapid Commun. Mass Spectrom. 24 3478–3484. 10.1002/rcm.4790 21072805 Suryanarayan A. (2020). What’s cooking in the Indus Civilisation? Investigating Indus food through ceramic lipid residue analysis. Ph.D. thesis. Cambridge: University of Cambridge. 10.17863/CAM.50249 28635591 Suryanarayan A. Cubas M. Craig O. E. Heron C. P. Shinde V. S. (2021). Organic residues in pottery from the Indus Civilisation in northwest India. J. Archaeol. Sci. 125:105291. 10.1016/j.jas.2020.105291 33519031 Thomas P. K. Joglekar P. P. Deshpande A. Pawankar S. J. Matsushima Y. (1997). Subsistence based on animals in the Harappan culture of Gujarat, India. Anthropozoologica 2 767–776. Thomas P. K. Matsushima Y. Deshpande A. (1996). “Faunal remains,” in Kuntasi: A Harappan Emporium on West Coast, eds Dhavalikar M. K. Raval M. H. Chitalwala Y. M. (Maharashtra: Deccan College), 297-330. Torrence R. Barton H. (eds) (2006). Ancient Starch Research. Walnut Creek: Left Coast Press. Twiss K. C. Bogaard A. Charles M. Henecke J. (2009). Plants and animals together: interpreting organic remains from building 52 at Çatalhöyük. Curr. Anthropol. 50 885–895. 10.1086/644767 VanDerwarker A. M. Detwiler K. R. (2002). Gendered practice in Cherokee foodways: a spatial analysis of plant remains from the Coweeta Creek Site. Southeast. Archaeol. 21 21–28. Weber S. Kashyap A. (2010). Harappan plant use revealed by starch grains from Farmana, India. Antiquity 84:326. Welch P. D. Scarry C. M. (1995). Status-related variation in foodways in the Moundville chiefdom. Am. Antiq. 60 397–419. 10.2307/282257 Whelton H. L. Roffet-Salque M. Kotsakis K. Urem-Kotsou D. Evershed R. P. (2018). Strong bias towards carcass product processing at Neolithic settlements in Northern Greece revealed through absorbed lipid residues of archaeological pottery. Quat. Int. 496 127–139. 10.1016/j.quaint.2017.12.018 Wright R. (2010). The Ancient Indus: Urbanism, Economy and Society. Cambridge: Cambridge University Press.