Stable isotope chemistry reveals plant-dominant diet among early foragers on the Andean Altiplano, 9.0–6.5 cal. ka

Abstract

Current models of early human subsistence economies suggest a focus on large mammal hunting. To evaluate this hypothesis, we examine human bone stable isotope chemistry of 24 individuals from the early Holocene sites of Wilamaya Patjxa (9.0–8.7 cal. ka) and Soro Mik’aya Patjxa (8.0–6.5 cal. ka) located at 3800 meters above sea level on the Andean Altiplano, Peru. Contrary to expectation, Bayesian mixing models based on the isotope chemistry reveal that plants dominated the diet, comprising 70–95% of the average diet. Paleoethnobotanical data further show that tubers may have been the most prominent subsistence resource. These findings update our understanding of earliest forager economies and the pathway to agricultural economies in the Andean highlands. The findings furthermore suggest that the initial subsistence economies of early human populations adapting to new landscapes may have been more plant oriented than current models suggest.

Introduction

The extent to which early human subsistence economies relied on meat versus plant foods is debated [1, 2]. Current understanding of the earliest subsistence economies of the Andean highlands suggest that meat was the major subsistence resource. Early Holocene assemblages, 11–5 cal. ka, consistently reveal abundant camelid and deer remains, projectile points, and scrapers suggesting hunting-oriented economies [3–11]. For example, Rick [4] concluded that, “The settlement pattern [of the Junín region] and the faunal collections [from the site of Pachamachay] strongly support the hypothesis that vicuña, or similar camelids, were the major food source for puna [ecosystem] hunter-gatherers.” Such observations are furthermore consistent with diet breadth models, which suggest that early hunter-gatherer populations would tend to target high-ranked large mammals before resorting to plant foods [12–14].

Recent analyses of materials from the Archaic Period sites of Soro Mik’aya Patjxa (8.0–6.5 cal. ka) and Wilamaya Patjxa (ca. 8.9 cal. ka), located on the Andean Altiplano (High Plateau), provide new opportunities to evaluate these economic models for early highland foragers. Similar to previous research, an abundance of projectile points, scrapers, and lithic debitage indicate considerable investment in the hunting of large terrestrial mammals, likely camelid and deer [15–17]. Preliminary zooarchaeological investigations reveal abundant large-mammal bone, consistent with a hypothesis of large mammal hunting [15, 18]. Groundstone artifacts, though informal and infrequent, suggest some degree of investment in plant resources [15, 18]. Distinctive dental wear patterns on the upper incisors, known as lingual surface attrition of the maxillary anterior teeth suggest intensive tuber processing [19]. Collectively, studies of the Soro Mik’aya Patjxa and Wilamaya Patjxa materials indicate diverse diets of large mammals and plants with an emphasis on large-mammal hunting, consistent with previous findings at other Andean highland sites.

Despite general agreement of various lines of evidence, the evidence remains indirect. Preservation biases that favor projectile points and animal bone compared to plant materials could, to some extent, inflate the hunting signal [20–22]. The biases of previous researchers who have generally been males from a culture in which hunting is a distinctly masculine pursuit could furthermore inflate the hunting signal [23]. It was ostensibly for this reason that ethnography famously revealed plant foods to play a prominent role in forager economies—Arctic economies aside—in contrast to earlier models that emphasized hunting [24–26]. Thus, current archaeological models and evidence leave considerable room for interpretive error.

A more direct but previously unexplored measure of early Andean diets is stable isotope chemistry of human bone. A study of bone isotope chemisty of six early Holocene (8.2–8.0 cal. ka) individuals, including four children and two adults, from the Andean highlands of Argentina finds evidence for tuber and herbivore consumption with breastfeeding and environmental aridity enriching the isotopic values [27]. The current analysis examines the diets of 24 Archaic Period foragers at the highland archaeological sites of Soro Mik’aya Patjxa and Wilamaya Patjxa using stable carbon and nitrogen isotopes in conjunction with more traditional zooarchaeological and paleoethnobotanical approaches. These assemblages date to the early Holocene, collectively spanning 9.0–6.5 cal. ka [16, 18]. Given current models of early highland economies, which point to a mixed diet of animals and plants with an emphasis on large mammals, we should expect to observe the human osteological samples from Soro Mik’aya Patjxa and Wilamaya Patjxa to exhibit dietary carbon (δ¹³C_diet) and nitrogen (δ¹⁵N_diet) values between those of local fauna and flora with a bias toward the means of the faunal δ¹³C and δ¹⁵N values.

This expectation follows from well established relationships between the isotopic composition of human bone and the foods that humans consume [28, 29]. Stable nitrogen isotopes (δ¹⁵N) vary with trophic level. Stable carbon isotopes (δ¹³C) vary with photosynthetic pathway. Such isotopic values in human bone chemistry can provide insights into major subsistence resources including C₃ plants, C₄ plants, mammals, and fish. Although stable isotope analysis does not offer taxonomic specificity beyond those broad categories, by coupling isotopic insights with zooarchaeological and paleoethnobotanical insights, it may be possible to move beyond preservation biases to gain more accurate estimates of human diets. Given current understanding of highland Archaic diets, we expect zooarchaeological analysis of the Archaic Period sites to reveal an abundance of faunal remains including vicuña, guanaco, or taruca with few small mammal, fish, or bird remains. Paleoethnobotanical analyses should reveal abundant wild chenopod seeds or tuber remains.

Materials and methods

Recently discovered human burials and other cultural pit features—possibly roasting or storage pits—at the Early-Late Archaic Period archaeological sites of Soro Mik’aya and Wilamaya Patjxa afford an opportunity to evaluate models of early subsistence practices on the Andean Altiplano. A series of radiocarbon dates place Soro Mik’aya Patjxa securely in the Middle to Late Archaic Periods (8.0–6.5 cal. ka) [16]. Radiocarbon dates and artifact typology broadly place Wilamaya Patjxa in the Early to Late Archaic Periods (ca. 11–5 cal. ka), with two direct showing occupation around 9.0 cal. ka [18, 30].

Portions of the two sites were systematically excavated with site matrix and feature fill screened using 6 mm and 1 mm screens, respectively. For each cultural feature, 10-liter bulk soil samples were taken for flotation analysis unless the feature consisted of less than 10L, in which case all feature sediment was collected for flotation. Flotation procedures followed d’Alpoim Guedes et al. [31] and Lennstrom and Hastorf [32] using a modified version of Watson’s [33] flotation machine.

The excavations revealed 18 cultural pit features at Soro Mik’aya Patjxa and 39 at Wilamaya Patjxa. From these, 16 individuals were discovered at Soro Mik’aya Patjxa and 12 at Wilamaya Patjxa. Here, we describe the laboratory methods for the three analytical approaches including isotopic, zooarchaeological and paleoethnobotanical approaches.

Stable isotope analysis of human bone

Human bone samples were excavated and exported under Peruvian Ministry of Culture Permit numbers 064-2013-DGPA-VMPCIC/MC and 138-2015-VMPCIC/MC. Stable isotope chemistry is performed in four different labs including the University of California Davis Stable Isotopes Facility (UCDSIFS), University of California Irvine W.M. Keck Carbon Cycle Accelerator Mass Spectrometer (KCCAMS) Facility, the University of Arizona Accelerator Mass Spectrometry Lab (AMS), and the Penn State University Laboratory for Isotopes and Metals in the Environment (LIME). Collagen extraction for the UCDSIFS and LIME submissions is performed at the UC Davis Archaeometry Lab following the protocol of Eerkens et al. [34]. Collagen extraction for the KCCAMS submissions follows the protocol described by Haas et al. [18].

To assess the extent of diagenetic alterations to bone collagen, we consider atomic C/N ratios with the expectation that reliable readings will exhibit C/N ratios in the range of 3.1–3.6 [35]. As an additional quality control measure, we compare the resultant δ¹³C and δ¹⁵N values for three of the individuals—SMP9, SMP16, and WMP6—to values previously reported in radiocarbon analyses performed by The University of Arizona AMS laboratory [16] and the KCCAMS facility [18].

Baseline δ¹³C and δ¹⁵N values for candidate subsistence resources, including C₃ plants, C₄ plants, camelid, or freshwater fish, are compiled from published sources [36–40]. To the extent possible, control samples are restricted to archaeological samples from the central Andean highlands. For any modern samples included, δ¹³C values are corrected by +1.5‰ for Suess effects. For any lowland samples included, δ¹³C values are offset +2‰ and δ¹⁵N values -1.5‰ based on the regression equations of Szpak et al. [38]. δ¹³C bone collagen samples are adjusted using a -2.4‰ offset if terrestrial [41] and -3.7‰ offset if aquatic [42, 43] to adjust for meat-bone offset.

Bayesian mixing models are used to estimate the dietary composition of the Soro Mik’aya Patjxa and Wilamaya Patjxa individuals using the δ¹³C and δ¹⁵N values of the consumers and potential food resources. Children (n = 4) are excluded from the model to prevent breast-feeding effects from influencing the results. Previous research shows that trophic offsets can vary widely due to a variety of environmental and trophic effects [27, 44, 45]. Given the antiquity of the system under investigation, the mobility of humans, and the volatility of isoscapes over time, it is difficult to know which trophic correction factors apply to a dietary regime under consideration. We therefore consider the range of possible trophic enrichment factors along with variance terms for both δ¹³C and δ¹⁵N. For δ¹³C, we use trophic enrichment factors ranging between 4.5–6.0‰, in 0.5‰ increments, with a standard deviation of 0.63‰ [45]. For δ¹⁵N, we use trophic enrichment factors ranging between from 3–6‰, in 1‰ increments, with a standard deviation of 0.74‰ [45]. Considering all possible combinations of trophic enrichment factors results in 16 models for evaluation.

All models assume uniformed priors given the lack of prior knowledge on the relative dietary contributions of the broad resource classes. Although it might be tempting to draw on zooarchaeological and archaeobotanical data for informative priors, differential preservation of faunal and floral artifacts precludes this possibility.

Model runs assume both residual and process error, a chain length of 1 million, burn-in of 5000, thinning of 500, and three chains. For dietary estimates, we report median values and 95% credible intervals for each subsistence resource and each model. Model convergence is assessed using Gelmen-Rubin and Geweke diagnostics, and the models are compared to one another using leave-one-out cross-validation information criterion (LOOic) and Akaike Information Criterion (AIC) weights [46]. All computation is performed using R statistical computing language [47] with Bayesian mixing modeling performed using MixSIAR package [46]. Although other packages for mixing modeling are available (e.g., FRUITS [48]), we use MixSIAR [49] because of its currency, documentation, integration with R statistical computing language, and accessibility via open-source Linux operating systems [50]. All code is made available in S1 Data.

Zooarchaeological analysis

New faunal data are reported for Soro Mik’aya Patjxa. Wilamaya Patjxa data are derived from a previous investigation by Noe [18]. For newly reported materials, all animal bone is weighed and counted. Although abundant, the animal bone is highly fragmented, likely due to intensive processing, making more precise taxonomic identification difficult. The method presented here serves to broadly distinguish between human, small mammal, large mammal, bird, and fish bone. Animal bone fragments are distinguished from human bone based on bone macrostructure where (a) human bone tends to be more porous than animal bone, (b) cortical bone tends to be thicker relative to bone diameter in animals compared to humans, (c) diaphyseal trabecula tends to be present in human but absent in animal bone, and (d) human cranial vaults tend to have thick dipole while animal cranial vaults tend to be more compact [51]. More detailed faunal analysis is ongoing, but the coarse analytical approach taken here is sufficient to address the broad dietary question at hand.

Paleoethnobotanical analysis

All Soro Mik’aya Patjxa features, which include burial pits and pits of unknown function, are subject to macrobotanical analysis. Samples are sorted using a stereoscopic light microscope with 10 to 40X magnification. Due to environmental conditions in the central Altiplano and the antiquity and exposed nature of the sites, it is highly unlikely that uncarbonized plant remains would have preserved, so analysis is restricted to carbonized remains. Macrobotanical remains are sorted into different tissue categories including seeds, wood, and parenchyma. Parenchyma refers to plant storage tissue. These distinct carbonized tissues with thin-walled cells are believed to be tuber fragments by paleoethnobotanists working in the Altiplano. All specimens are identified to the most specific taxonomic level possible. Paleoethnobotanical analysis is restricted to Soro Mik’aya Patjxa with Wilamaya Patjxa paleoethnobotanical analysis ongoing.

Results

Isotopic control samples compiled from the Andean literature include 96 large mammals, 84 C₃ plants, 29 C₄ plants, and 10 fish samples. The data reveal strong clustering of carbon and nitrogen values by category (Fig 1, Table 1, S1 Table) providing an ideal baseline for comparison with the human bone samples reported here (Table 2). Twenty-four individuals from Soro Mik’aya Patxja (SMP) and Wilamaya Patjxa (WMP) show δ¹⁵N_diet values ranging from 2.0‰ to 8.3‰ with a mean of 3.4‰ (see Fig 1A) and δ¹³C_diet values ranging from -24.3‰ to -22.9‰ with a mean of -23.7‰ (see Fig 1B).

[Figure omitted. See PDF.]

a) δ¹³C_diet values are consistent with those of C₃ plants with slight enrichment from some other resource types. b) δ¹⁵N_diet values are most consistent with those of plants. c) Biplot of δ¹³C and δ¹⁵N values are consistent with a mixed diet principally based on C₃ plants with low levels of enrichment from some other resource. Dietary values assume δ¹³C TEF = 5.0‰, δ¹⁵N TEF = 6.0‰ based on mixing model results (see Table 4). Dots = individual samples, ellipses = 95% variance ranges for each category, and crosshairs = mean values by category (see Tables 1 and 2).

[Figure omitted. See PDF.]

See S1 Table for sample data.

[Figure omitted. See PDF.]

See S2 Table for additional metadata.

Quality control measures indicate that the archaeological stable isotope values reported here are reliable. Atomic C/N ratios for all but one sample—WMP1—fall within the acceptable range of 3.1–3.6 (see Table 2), indicating that diagenetic processes have not significantly altered the collagen [35]. Furthermore, the δ¹³C and δ¹⁵N values from the four previously reported radiocarbon dates [16, 18, 30] are in close agreement with less than 1.3‰ separating values reported among the three labs (Table 3).

[Figure omitted. See PDF.]

UCDSIFS compared to AMS and KCCAMS.

The δ¹³C and δ¹⁵N values for the human individuals fall closest to the mean of C₃ plants with slight enrichment from some other set of resources indicating C₃ dominant diets and excluding the possibilities that mammals, C₄ plants, or fish comprised large portions of the diets (see Fig 1C). All 16 of the Bayesian mixing models indicate that adult diets were dominated by C₃ plants (Table 4). Median C₃ plant contribution estimates range from 60–95% and median mammal estimates range from 3–34% for all 16 models. Gelman and Geweke model diagnostics are consistently zero or near zero, indicating that all models are plausible.

[Figure omitted. See PDF.]

The best-fit Bayesian mixing model indicates that C₃ plants comprised approximately 84% (73–92%) of the average adult diet with meat comprising just 9% (0–24%), fish 4% (0–13%), and C₄ plants 2% (0–6%; Fig 2). This model generated the lowest LOOic value and the highest AIC weight, and is based on a δ¹⁵N trophic enrichment factor of 6.00±0.74‰ and a δ¹³C trophic enrichment factor of 5.00±0.63‰. However, five other models produced nearly equivalent AIC weights greater than 0.10 suggesting virtually equivalent model performance. Among these models, median C₃ plant estimates range from 70–95% and mammal estimates range between 3–23%. These results show that the particular trophic enrichment factors, ranging from 3.0–6.0‰ for δ¹⁵N and 4.5–6.0 for δ¹³C, have little effect on the broad dietary estimates. All models indicate a plant dominant diet with median values for C₃ plants ranging between 70–95% and mammals ranging between 3–23%. Thus all credible subsistence models indicate that plant foods comprised the majority of individual diets and meat played a secondary role. These findings are inconsistent with the working hypothesis of a meat-dominant diet and instead suggest a plant-dominant diet among early forager populations of the Andean Altiplano, 9.0–6.5 cal. ka.

[Figure omitted. See PDF.]

a. model 14, b. model 9, c. model 13, d. model 15, e. model 16, f. model 10 (see Table 4).

Zooarchaeological and paleoethnobotanical observations offer some taxonomic precision beyond the broad food categories used in the stable isotope analysis. Excavations at Soro Mik’aya produced 3193 fragments of animal bone (number of individual specimens or NISP) from across the site. Of the total specimens, the most frequently identified category included 200 large mammal fragments with only trace amounts of small mammal, bird, and fish bones reported [15]. Most of these large mammal bones were burned (n = 150; 75%). A previously published assessment of 341 faunal bone elements from Wilamaya Patjxa faunal assemblage similarly revealed that camelid and deer bone were the most frequently identified taxa with 17 camelid and 5 deer elements observed [18]. Notably absent, again, are small mammals, birds, and fish. These data indicate that the slight enrichment observed in the carbon and nitrogen values in the human bone was likely due to large mammal consumption and not small mammals, fish, or birds.

The most abundant plant-food specimens in the Soro Mik’aya Patjxa paleoethnobotanical assemblage is parenchyma tissue, identified as tuber fragments with 52 specimens found in 9 of 15 features (Fig 3, Table 5). Chenopod seeds are nearly absent with just three wild specimens observed among three features. The most abundant paleoethnobotanical samples are non-food resources including wood fragments (n = 448) and grass (Poaceae) seeds (n = 161) observed in nearly every feature. The wood most likely reflects use as fuel. Similarly the small grass seeds likely reflect fuel use whether from burning grasses or the dung of camelids [57]. The small grass seeds are typical of the region, and none of the taxa are suitable for human consumption but are excellent forage for camelids. This finding is consistent with a model of early wild tuber consumption and other lines of evidence that suggest intensive tuber use in the region during the Archaic Period [19, 58, 59].

[Figure omitted. See PDF.]

Ubiquity is measured as the proportion of archaeological features containing each taxon. Values based on 15 features, 300L of flotation, and 688 paleoethnobotanical artifacts (see Table 5).

[Figure omitted. See PDF.]

Discussion

The stable isotope, faunal, and paleoethnobotanical evidence from the sites of Soro Mik’aya Patjxa and Wilamaya Patjxa converge to indicate that C₃ plants, likely wild tubers, comprised the major component of early forager diets on the Andean Altiplano and that meat, including vicuña and taruca, played a secondary role. C₄ plants, small mammals, fish, and birds appear to have played negligible roles in these early subsistence economies. The findings presented here depart from current thinking about early Andean highland diets and force a reconsideration of existing economic models.

One possible explanation for the unexpected emphasis on plant foods among this early highland population may be that large mammal populations had been severely reduced by 9 cal. ka. Although the Wilamaya Patjxa assemblage includes the earliest archaeological period of the region—the Early Archaic Period, 11.0–9.0 cal. ka—it is restricted to the latter end of the period [18]. It is currently unclear when human populations first arrived in the region [60], but if they arrived as early as 11 cal. ka., then humans would have been hunting the region for 2,000 years prior to the earliest individuals under investigation. This would certainly have been enough time to decimate the region’s animal population in the absence of animal conservation strategies [4]. A second possibility is that the earliest populations simply did not engage in hunting to the extent previously thought. Previous research suggests that prey choice models may over-estimate the dietary values of large mammals whether due to physical risk [61], prey behavior [62], or economic risks associated with long encounter intervals [63]. A third explanatory possibility is that early highland populations relied heavily on hunting large mammals but incorporated animal digesta into their diets [64], which could simultaneously account for both the archaeological signatures of hunting—animal bone and projectile points—and the depleted nitrogen values in human bone chemistry observed here. Evaluating these hypotheses will ultimately require investigation of earlier archaeological assemblages.

These findings furthermore hold implications for our understanding of domestication in the high Andes. Pearsall [65] and Kuznar [66] proposed that plant management commenced after camelids were managed in corrals. In this view, camelids transported plant seeds to corrals where they thrived in soils fertilized by camelid dung, creating a mutually beneficial relationship between camelids, plants, and human communities. Versions of this model suggest that Chenopodium spp. (including the crops quinoa and kañawa) and the tuber maca (Lepidium mevenii Walp.) were domesticated in this way. These coevolutionary processes also likely led to the domestication of potatoes (Solanum tuberosum L.) and up to 15 other species of roots and tubers in the Andes [67]. Consistent with this model, chenopod seeds, maca tubers, and managed camelids appeared after about 4,000 years ago at the rock shelter site of Panalauca where the sizes of chenopod and maca specimens increased through time [65]. Recent research into domestication mutualism supports the Andean version of the camp follower hypothesis showing that chenopods, tubers, and camelids were likely domesticated in tandem as complementary foods [68].

While current models suggest initial co-evolutionary processes involving maca, chenopod, and camelid intensification in the highlands, the paleoethnobotanical evidence presented here fails to find strong evidence of early intensive chenopod use on the Altiplano during the Archaic period. This may reflect preservation biases given the small and delicate nature of wild chenopod seeds. Alternatively, it may be that chenopods did not become economically important in the region until sometime after 6.5 cal. ka. A later incorporation of chenopods into the diet would be consistent with prey choice models given the low post-encounter return rates of small seeds relative to tubers [69].

The finding of tuber fragments at Soro Mik’aya Patjxa is consistent with the role of tubers in the early stages of the co-evolutionary process. However, the tuber fragments are unlikely to be maca, which is a more northerly taxon. Several tuber species could potentially account for the tubers observed at Soro Mik’aya Patjxa, but they are most likely associated with wild potato species, which are concentrated in the region and were likely domesticated there [59, 70–72]. Additional sites on the Altiplano should be examined with particular attention to contemporaneous sites that would assess replicability of the current findings and to non-contemporaneous sites that would afford diachronic comparisons.

In sum, the results presented here are consistent with a model of human-tuber-camelid co-evolutionary dynamics beginning approximately 9,000 years ago on the Andean Altiplano. These findings support an updated model of Archaic Period subsistence practices in the central high Andes in which forager subsistence economies 9.0–6.5 cal. ka. emphasized plant foraging with lesser attention to large mammal hunting and a virtual absence of small animal hunting and fishing. This resource base may have catalyzed potato and camelid domestication in the subsequent Terminal Archaic Period after 5 cal. ka [58].

These findings further highlight the need to re-evaluate anthropological understanding of early forager diets more generally. Current perspectives vary with some models emphasizing the primacy of plants and others of animals [1, 2] with plant foraging becoming increasingly important relatively late in time on the eve of agriculture [73]. This may still be so, but the current analysis suggests that the shift to plant-foraging economies may have happened relatively rapidly, evidently having transpired in less than 2,000 years in the Andean case. This observation resonates with recent archaeological theory and findings that reveal a prominent role for plant foods in early forager diets [31, 74, 75] and ethnographic findings of the 1960s when, contrary to dominant thinking of the time, many subsistence economies once thought to be meat-dominant were shown to be plant-dominant [24, 25]. Stable isotope chemistry gives archaeologists the opportunity to reliably extend such investigations into the deep past. The current study arrives at a similar place to the earlier ethnographic findings—plant foods were central to early human economies.

Supporting information

S1 Data. Bayesian mixing model code in R language.

https://doi.org/10.1371/journal.pone.0296420.s001

(R)

S1 Table. Isotopic control data.

https://doi.org/10.1371/journal.pone.0296420.s002

(CSV)

S2 Table. Isotopic data for Soro Mik’aya Patjxa and Wilamaya Patjxa individuals.

https://doi.org/10.1371/journal.pone.0296420.s003

(CSV)

Acknowledgments

Field support was provided by Collasuyo Archaeological Research Institute, and the communities of Mulla Fasiri and Totorani, Peru. Bryna Hull (UC Davis) provided support for stable isotope analysis. Samples of human remains were excavated and exported under Peruvian Ministry of Culture Permit numbers 064-2013-DGPA-VMPCIC/MC and 138-2015-VMPCIC/MC.

Citation: Chen JC, Aldenderfer MS, Eerkens JW, Langlie BS, Viviano Llave C, Watson JT, et al. (2024) Stable isotope chemistry reveals plant-dominant diet among early foragers on the Andean Altiplano, 9.0–6.5 cal. ka. PLoS ONE 19(1): e0296420. https://doi.org/10.1371/journal.pone.0296420

About the Authors:

Jennifer C. Chen

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Anthropology, The Pennsylvania State University, State College, PA, United States of America

ORICD: https://orcid.org/0000-0001-7402-5946

Mark S. Aldenderfer

Roles: Data curation, Funding acquisition, Investigation, Resources, Writing – review & editing

Affiliation: Department of Anthropology and Heritage Studies, University of California, Merced, CA, United States of America

Jelmer W. Eerkens

Roles: Methodology, Resources, Supervision, Validation, Writing – review & editing

Affiliation: Department of Anthropology, University of California Davis, Davis, CA, United States of America

BrieAnna S. Langlie

Roles: Formal analysis, Investigation, Methodology, Resources, Writing – review & editing

Affiliation: Department of Anthropology, Binghamton University, Binghamton, NY, United States of America

Carlos Viviano Llave

Roles: Data curation, Investigation, Project administration, Writing – review & editing

Affiliation: National Register of Peruvian Archaeologists, Lima, Peru

James T. Watson

Roles: Formal analysis, Investigation, Methodology, Writing – review & editing

Affiliation: Arizona State Museum and School of Anthropology, University of Arizona, Tucson, AZ, United States of America

Randall Haas

Roles: Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

E-mail: [email protected]

Affiliation: Department of Anthropology, University of Wyoming, Laramie, WY, United States of America

ORICD: https://orcid.org/0000-0002-7742-2127

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

[/RAW_REF_TEXT]

References

1. Speth JD. The paleoanthropology and archaeology of big-game hunting. New York: Springer New York; 2010. https://doi.org/10.1007/978-1-4419-6733-6

2. Pontzer H, Wood BM, Raichlen DA. Hunter-gatherers as models in public health. Obesity Reviews 2018;19:24–35. pmid:30511505

3. Lynch TF editor. Guitarrero Cave: early man in the Andes. New York: Academic Press; 1980.

4. Rick JW. Prehistoric hunters of the high Andes. New York: Academic Press; 1980.

5. MacNeish RS, Vierra RK, Nelken-Terner A, Garcia Cook A. Prehistory of the Ayacucho Basin, Peru: the Preceramic way of life(Vol. IV). Ann Arbor: University of Michigan Press; 1983.

6. Lavallée D, Julien M, Wheeler JC, Karlin C. Telarmachay: cazadores y pastores prehistóricos de los Andes(Vol. I). Lima: Instituto Francés de Estudios Andinos; 1995.

7. Aldenderfer MS. Montane foragers. Iowa City: University of Iowa Press; 1998. https://doi.org/10.2307/j.ctt20q1wq5

8. Núñez L, Grosjean M, Cartajena I. Human Occupations and Climate Change in the Puna de Atacama, Chile. Science 2002;298:821–824. pmid:12399589

9. Rick JW, Moore KM. El Precerámico de las punas de Junín: el punto de vista desde Panaulauca. Boletín de Arqueología PUCP 2012;3:263–296.

10. Rademaker K, Hodgins G, Moore K, Zarrillo S, Miller C, Bromley GRM, et al. Paleoindian settlement of the high-altitude Peruvian Andes. Science 2014;346:466–469. pmid:25342802

11. Yacobaccio HD. Peopling of the high Andes of northwestern Argentina. Quaternary International 2017;461:34–40.

12. Alvard MS, Kuznar L. Deferred harvests: the transition from hunting to animal husbandry. American Anthropologist 2001;103:295–311.

13. Aldenderfer MS. Costly signaling, the sexual division of labor, and animal domestication in the Andean highlands. In: Kennett DJ, Winterhalder B editors. Behavioral ecology and the transition to agriculture. Berkeley: University of California Press; 2006. https://doi.org/10.1525/9780520932456-011

14. Bird DW O’Connell JF. Behavioral Ecology and Archaeology. Journal of Archaeological Research 2006;14:143–188.

15. Haas R, Viviano Llave C. Hunter-gatherers on the eve of agriculture: investigations at Soro Mik’aya Patjxa, Lake Titicaca Basin, Peru, 8000–6700 BP. Antiquity 2015;89:1297–1312.

16. Haas R, Stefanescu IC, Garcia-Putnam A, Aldenderfer MS, Clementz MT, Murphy MS, et al. Humans permanently occupied the Andean highlands by at least 7 ka. Royal Society Open Science 2017;4:170331. pmid:28680685

17. Kitchel N, Aldenderfer MS, Haas R. Diet, mobility, technology, and lithics: neolithization on the Andean Altiplano, 7.0–3.5 ka. Journal of Archaeological Method and Theory 2021;29:390–425.

18. Haas R, Watson J, Buonasera T, Southon J, Chen JC, Noe S, et al. Female hunters of the early Americas. Science Advances 2020;6:eabd0310. pmid:33148651

19. Watson JT, Haas R. Dental evidence for wild tuber processing among Titicaca Basin foragers 7000 ybp. American Journal of Physical Anthropology 2017;164:117–130. pmid:28581062

20. Kornfeld M. The big-game focus: reinterpreting the archaeological record of Cantabrian Upper Paleolithic economy. Current Anthropology 1996;37:629–657.

21. Waguespack NM. The organization of male and female labor in foraging societies: implications for early Paleoindian archaeology. American Anthropologist 2005;107:666–676.

22. Adovasio JM, Soffer O, Illingworth JS, Hyland DC. Perishable fiber artifacts and Paleoindians: new implications. North American Archaeologist 2014;35:331–352.

23. Gero JM, Conkey MWeditors. Engendering archaeology. Oxford: Blackwell Publishers; 1991.

24. Lee RB. The! Kung San: men, women and work in a foraging society. Cambridge: Cambridge University Press; 1979.

25. Roosevelt AC. Ecology in human evolution: the origins of humans and their complex societies. In: Scarborough VL editor. A catalyst for ideas: anthropological archaeology and the legacy of Douglas Schwartz. Santa Fe: School for America Research Press; 2005.

26. Lee RB. Hunter-gatherers and human evolution: new light on old debates. Annual Review of Anthropology 2018;47:513–531.

27. Killian Galván V, Martnez J, Cherkinsky A, Mondini M, Panarello H. Stable isotope analysis on human remains from the final early Holocene in the southern Puna of Argentina: the case of Peñas de las Trampas1.1. Environmental Archaeology 2015;21:1–10.

28. Richards MP, Trinkaus E. Isotopic evidence for the diets of European Neanderthals and early modern humans. Proceedings of the National Academy of Sciences 2009;106:16034–16039.

29. Schoeninger MJ. Stable isotope analyses and the evolution of human diets. Annual Review of Anthropology 2014;43:413–430.

30. Smallwood A, Haas R, Jennings T. Lithic usewear confirms the function of Wilamaya Patjxa projectile points. Scientific Reports 2023;13. pmid:37923759

31. d’Alpoim Guedes J, Lu H, Li Y, Spengler RN, Wu X, Aldenderfer MS. Moving agriculture onto the Tibetan Plateau: the archaeobotanical evidence. Archaeological and Anthropological Sciences 2013;6:255–269.

32. Lennstrom HA, Hastorf CA. Testing old wives’ tales in palaeoethnobotany: A comparison of bulk and scatter sampling schemes from Pancán, Peru. Journal of Archaeological Science 1992;19:205–229.

33. Watson PJ. In pursuit of prehistoric subsistence: a comparative account of some contemporary flotation techniques. Midcontinental Journal of Archaeology 1976;1:77–100.

34. Eerkens JW, de Voogt A, Dupras TL, Rose SC, Bartelink EJ, Francigny V. Intra- and inter-individual variation in updelta 13 C and updelta 15 N in human dental calculus and comparison to bone collagen and apatite isotopes. Journal of Archaeological Science 2014;52:64–71.

35. DeNiro MJ. Postmortem preservation and alteration of in vivo bone collagen isotope ratios in relation to palaeodietary reconstruction. Nature 1985;317:806–809.

36. Miller MJ, Capriles JM, Hastorf CA. The fish of Lake Titicaca: implications for archaeology and changing ecology through stable isotope analysis. Journal of Archaeological Science 2010;37:317–327.

37. Barberena R, Méndez C, Mena F, Reyes O. Endangered species, archaeology, and stable isotopes: huemul (Hippocamelus bisulcus) isotopic ecology in central-western Patagonia (South America). Journal of Archaeological Science 2011;38:2313–2323.

38. Szpak P, White CD, Longstaffe FJ, Millaire J-F, Sánchez VFV. Carbon and nitrogen isotopic survey of northern Peruvian plants: baselines for paleodietary and paleoecological studies. PLoS ONE 2013;8:e53763. pmid:23341996

39. Grant J. Of Hunting and Herding: Isotopic evidence in wild and domesticated camelids from the Southern Argentine Puna (2120–420 years BP). Journal of Archaeological Science: Reports 2017;11:29–37.

40. Miller MJ, Kendall I, Capriles JM, Bruno MC, Evershed RP, Hastorf CA. Quinoa, potatoes, and llamas fueled emergent social complexity in the Lake Titicaca Basin of the Andes. Proceedings of the National Academy of Sciences 2021;118. pmid:34845028

41. DeNiro MJ, Epstein S. Influence of diet on the distribution of carbon isotopes in animals. Geochimica et Cosmochimica Acta 1978;42:495–506.

42. Keegan WF, DeNiro MJ. Stable carbon- and nitrogen-isotope ratios of bone collagen used to study coral-reef and terrestrial components of prehistoric Bahamian diet. American Antiquity 1988;53:320–336.

43. Guiry EJ, Hunt BPV. Integrating fish scale and bone isotopic compositions for `deep time’ retrospective studies. Marine Environmental Research 2020;160:104982. pmid:32907720

44. Ugan A, Coltrain J. Variation in collagen stable nitrogen values in black-tailed jackrabbits (Lepus californicus) in relation to small-scale differences in climate, soil, and topography. Journal of Archaeological Science 2011;38:1417–1429.

45. Cheung C, Szpak P. Interpreting past human diets using stable isotope mixing models. Journal of Archaeological Method and Theory 2020;28:1106–1142.

46. Stock BC, Jackson AL, Ward EJ, Parnell AC, Phillips DL, Semmens BX. Analyzing mixing systems using a new generation of Bayesian tracer mixing models. PeerJ 2018;6:e5096. pmid:29942712

47. R Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, https://www.R-project.org/. Vienna 2023.

48. Fernandes R, Millard AR, Brabec M, Nadeau M-J, Grootes P. Food reconstruction using isotopic transferred signals (FRUITS): a bayesian model for diet reconstruction. PLoS ONE 2014;9:e87436. pmid:24551057

49. Stock B, Jackson A, Ward E, Venkiteswaran J. MixSIAR 3.1.9.2018.

50. Cheung C, Szpak P. Interpreting past human diets using stable isotope mixing models—best practices for data acquisition. Journal of Archaeological Method and Theory 2022;29:138–161.

51. Mulhern DM. Differentiating human from nonhuman skeletal remains. In: Ubelaker SB&DH editor. Handbook of forensic anthropology and archaeology. Walnut Creek: Left Coast Press; 2009.

52. Killian Galván V, Grant JL, Escola P, Panarello HO, Olivera DE. Análisis de paleodietas humanas en zonas áridas a través de isótopos estables: el caso de Antofagasta de la Sierra (noroeste argentino). Revista Colombiana de Antropologa 2016;52:199–227.

53. López M. P, Cartajena I, Loyola R, Nú nez L, Carrasco C. The use of hunting and herding spaces: stable isotope analysis of late archaic and early formative camelids in the Tulan transect (Puna de Atacama, Chile). International Journal of Osteoarchaeology 2017;27:1059–1069.

54. Grant J, Mondini M, Panarello HO. Carbon and nitrogen isotopic ecology of Holocene camelids in the southern puna (Antofagasta de la Sierra, Catamarca, Argentina): archaeological and environmental implications. Journal of Archaeological Science: Reports 2018;18:637–647.

55. Buikstra JE, Ubelaker DH. Standards for data collection from human skeletal remains: Proceedings of a seminar at the Field Museum of Natural History. Fayetteville: Arkansas Archeological Survey; 1994.

56. Hogg AG, Heaton TJ, Hua Q, Palmer JG, Turney CSM, Southon J, et al. SHCal20 southern hemisphere calibration, 0–55,000 years cal BP. Radiocarbon 2020;62:759–778.

57. Bruno MC, Hastorf CA. Gifts from the camelids: archaeobotanical insights into camelid pastoralism through the study of dung. In: Capriles Flores JM, Tripcevich Neditors. The archaeology of Andean pastoralism. Albuquerque: University of New Mexico Press; 2016.

58. Rumold CU, Aldenderfer MS. Late Archaic–Early Formative Period microbotanical evidence for potato at Jiskairumoko in the Titicaca Basin of southern Peru. Proceedings of the National Academy of Sciences 2016;113:13672–13677. pmid:27849582

59. Jorgensen K, Garcia OA, Kiyamu M, Brutsaert TD, Bigham AW. Genetic adaptations to potato starch digestion in the Peruvian Andes. American Journal of Biological Anthropology 2022;180:162–172.

60. Haas R. Early settlement in the high Andes. Oxford Research Encyclopedia of Anthropology 2023.

61. Lupo KD, Schmitt DN. When bigger is not better: The economics of hunting megafauna and its implications for Plio-Pleistocene hunter-gatherers. Journal of Anthropological Archaeology 2016;44:185–197.

62. Lupo KD, Schmitt DN, Madsen DB. Size matters only sometimes: the energy-risk trade-offs of Holocene prey acquisition in the Bonneville basin, western USA. Archaeological and Anthropological Sciences 2020;12.

63. Bird DW, Bird RB, Codding BF. In pursuit of mobile prey: Martu hunting strategies and archaeofaunal interpretation. American Antiquity 2009;74:3–29.

64. Garvey R. Human consumption of large herbivore digesta and its implications for foraging theory. Evolutionary Anthropology: Issues, News, and Reviews 2023:1–9. pmid:37066856

65. Pearsall DM. Adaptation of prehistoric hunter-gatherers to the high Andes: the changing role of plant resources. In: Harris DR, Hillman GCeditors. Foraging and farming: the evolution of plant exploitation.: Routledge; 1989.

66. Kuznar LA. Mutualism between chenopodium, herd animals, and herders in the South Central Andes. Mountain Research and Development 1993;13:257–65.

67. Flores HE, Walker TS, Guimaraes RL, Bais HP, Vivanco JM. Andean root and tuber crops: underground rainbows. HortScience 2003;38:161–167.

68. Langlie BS, Capriles JM. Paleoethnobotanical evidence points to agricultural mutualism among early camelid pastoralists of the Andean central Altiplano. Archaeological and Anthropological Sciences 2021;13:107.

69. Joyce K, Louderback LA, Robinson E. Direct evidence for geophyte exploitation in the Wyoming Basin. American Antiquity 2021;87:236–247.

70. Hawkes JG. The domestication of roots and tubers in the American tropics. In: Harris DR, Hillman GCeditors. Foraging and farming: the evolution of plant exploitation.: Routledge; 1989.

71. Hijmans RJ, Spooner DM. Geographic distribution of wild potato species. American Journal of Botany 2001;88:2101–2112. pmid:21669641

72. Spooner DM, McLean K, Ramsay G, Waugh R, Bryan GJ. A single domestication for potato based on multilocus amplified fragment length polymorphism genotyping. Proceedings of the National Academy of Sciences 2005;102:14694–14699. pmid:16203994

73. Stiner MC. Thirty years on the "broad spectrum revolution" and Paleolithic demography. Proceedings of the National Academy of Sciences 2001;98:6993–6996. pmid:11390968

74. Henry AG, Brooks AS, Piperno DR. Microfossils in calculus demonstrate consumption of plants and cooked foods in Neanderthal diets (Shanidar III, Iraq; Spy I and II, Belgium). Proceedings of the National Academy of Sciences 2010;108:486–491. pmid:21187393

75. Hardy K, Bocherens H, Miller JB, Copeland L. Reconstructing Neanderthal diet: The case for carbohydrates. Journal of Human Evolution 2022;162:103105. pmid:34923240