Ancient DNA reveals the lost domestication history of South American camelids in Northern Chile and across the Andes
Abstract
The study of South American camelids and their domestication is a highly debated topic in zooarchaeology. Identifying the domestic species (alpaca and llama) in archaeological sites based solely on morphological data is challenging due to their similarity with respect to their wild ancestors. Using genetic methods also presents challenges due to the hybridization history of the domestic species, which are thought to have extensively hybridized following the Spanish conquest of South America that resulted in camelids slaughtered en-masse. In this study we generated mitochondrial genomes for 61 ancient South American camelids dated between 3,500 - 2,400 years before the present (Early Formative period) from two archaeological sites in Northern Chile (Tulán-54 and Tulán-85), as well as 66 modern camelid mitogenomes and 815 modern mitochondrial control region sequences from across South America. In addition, we performed osteometric analyses to differentiate big and small body size camelids. A comparative analysis of these data suggests that a substantial proportion of the ancient vicuña genetic variation has been lost since the Early Formative period as it is not present in modern specimens. Moreover, we propose a domestication hypothesis that includes an ancient guanaco population that no longer exists. Finally, we find evidence that interbreeding practices were widespread during the domestication process by the early camelid herders in the Atacama during the Early Formative period and predating the Spanish conquest.
Data availability
Sequencing data have been deposited in GenBank under accession codes:BankIt2325854 C10mito MT254135BankIt2325854 C11mito MT254136BankIt2325854 C12mito MT254137BankIt2325854 c13mito MT254138BankIt2325854 C14mito MT254139BankIt2325854 C15mito MT254140BankIt2325854 C16mito MT254141BankIt2325854 C18mito MT254142BankIt2325854 c1mito MT254143BankIt2325854 C20mito MT254144BankIt2325854 C21mito MT254145BankIt2325854 C22mito MT254146BankIt2325854 C24mito MT254147BankIt2325854 C25mito MT254148BankIt2325854 C27mito MT254149BankIt2325854 C28mito MT254150BankIt2325854 C29mito MT254151BankIt2325854 C30mito MT254152BankIt2325854 C31mito MT254153BankIt2325854 C32mito MT254154BankIt2325854 C33mito MT254155BankIt2325854 C34mito MT254156BankIt2325854 C35mito MT254157BankIt2325854 C36mito MT254158BankIt2325854 C37mito MT254159BankIt2325854 C38mito MT254160BankIt2325854 C39mito MT254161BankIt2325854 C41mito MT254162BankIt2325854 C42mito MT254163BankIt2325854 C43mito MT254164BankIt2325854 C44mito MT254165BankIt2325854 C46mito MT254166BankIt2325854 C47mito MT254167BankIt2325854 C48mito MT254168BankIt2325854 C49mito MT254169BankIt2325854 C50mito MT254170BankIt2325854 C52mito MT254171BankIt2325854 C54mito MT254172BankIt2325854 C55mito MT254173BankIt2325854 C57mito MT254174BankIt2325854 C59mito MT254175BankIt2325854 C5mito MT254176BankIt2325854 C60mito MT254177BankIt2325854 C61mito MT254178BankIt2325854 C62mito MT254179BankIt2325854 C63mito MT254180BankIt2325854 C64mito MT254181BankIt2325854 C65mito MT254182BankIt2325854 C66mito MT254183BankIt2325854 C67mito MT254184BankIt2325854 C68mito MT254185BankIt2325854 C69mito MT254186BankIt2325854 C6mito MT254187BankIt2325854 C70mito MT254188BankIt2325854 C71mito MT254189BankIt2325854 C72mito MT254190BankIt2325854 C73mito MT254191BankIt2325854 C74mito MT254192BankIt2325854 C75mito MT254193BankIt2325854 C76mito MT254194BankIt2325854 C77mito MT254195
Article and author information
Author details
Funding
Fondo Nacional de Desarrollo Científico, Tecnológico y de Innovación Tecnológica (1070040)
- Paloma Fernández Diaz-Maroto
Fondo Nacional de Desarrollo Científico, Tecnológico y de Innovación Tecnológica (1020316)
- Paloma Fernández Diaz-Maroto
Fondo Nacional de Desarrollo Científico y Tecnológico (1130917)
- Paloma Fernández Diaz-Maroto
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Reviewing Editor
- Jessica C Thompson, Yale University, United States
Publication history
- Received: September 23, 2020
- Accepted: March 15, 2021
- Accepted Manuscript published: March 16, 2021 (version 1)
- Version of Record published: April 8, 2021 (version 2)
Copyright
© 2021, Diaz-Maroto et al.
This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,217
- Page views
-
- 404
- Downloads
-
- 10
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
- Microbiology and Infectious Disease
Increasing numbers of small proteins with diverse physiological roles are being identified and characterized in both prokaryotic and eukaryotic systems, but the origins and evolution of these proteins remain unclear. Recent genomic sequence analyses in several organisms suggest that new functions encoded by small open reading frames (sORFs) may emerge de novo from noncoding sequences. However, experimental data demonstrating if and how randomly generated sORFs can confer beneficial effects to cells are limited. Here, we show that by upregulating hisB expression, de novo small proteins (≤50 amino acids in length) selected from random sequence libraries can rescue Escherichia coli cells that lack the conditionally essential SerB enzyme. The recovered small proteins are hydrophobic and confer their rescue effect by binding to the 5′ end regulatory region of the his operon mRNA, suggesting that protein binding promotes structural rearrangements of the RNA that allow increased hisB expression. This study adds RNA regulatory elements as another interacting partner for de novo proteins isolated from random sequence libraries and provides further experimental evidence that small proteins with selective benefits can originate from the expression of nonfunctional sequences.
-
- Evolutionary Biology
- Microbiology and Infectious Disease
Many bacterial genomes carry prophages whose induction can eliminate competitors. In response, bacteria may become resistant by modifying surface receptors, by lysogenization, or by other poorly known processes. All these mechanisms affect bacterial fitness and population dynamics. To understand the evolution of phage resistance, we co-cultivated a phage-sensitive strain (BJ1) and a poly-lysogenic Klebsiella pneumoniae strain (ST14) under different phage pressures. The population yield remained stable after 30 days. Surprisingly, the initially sensitive strain remained in all populations and its frequency was highest when phage pressure was strongest. Resistance to phages in these populations emerged initially through mutations preventing capsule biosynthesis. Protection through lysogeny was rarely observed because the lysogens have increased death rates due to prophage induction. Unexpectedly, the adaptation process changed at longer time scales the frequency of capsulated cells in BJ1 populations increased again, because the production of capsule was fine-tuned, reducing the ability of phage to absorb. Contrary to the lysogens, these capsulated resistant clones are pan-resistant to a large panel of phages. Intriguingly, some clones exhibited transient non-genetic resistance to phages, suggesting an important role of phenotypic resistance in coevolving populations. Our results show that interactions between lysogens and sensitive strains are shaped by antagonistic co-evolution between phages and bacteria. These processes may involve key physiological traits, such as the capsule, and depend on the time frame of the evolutionary process. At short time scales, simple and costly inactivating mutations are adaptive, but in the long-term, changes drawing more favorable trade-offs between resistance to phages and cell fitness become prevalent.