Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Retroviral DNA Sequences as a Means for Determining Ancient Diets

  • Jessica I. Rivera-Perez ,

    jessica.i.rivera88@gmail.com

    Affiliation Department of Biology, University of Puerto Rico, Rio Piedras Campus, San Juan, Puerto Rico

  • Raul J. Cano,

    Affiliation Center for Applications in Biotechnology, Biological Sciences Department, California Polytechnic State University, San Luis Obispo, California, United States of America

  • Yvonne Narganes-Storde,

    Affiliation Center for Archaeological Research, University of Puerto Rico, Rio Piedras Campus, San Juan, Puerto Rico

  • Luis Chanlatte-Baik,

    Affiliation Center for Archaeological Research, University of Puerto Rico, Rio Piedras Campus, San Juan, Puerto Rico

  • Gary A. Toranzos

    Affiliation Department of Biology, University of Puerto Rico, Rio Piedras Campus, San Juan, Puerto Rico

Abstract

For ages, specialists from varying fields have studied the diets of the primeval inhabitants of our planet, detecting diet remains in archaeological specimens using a range of morphological and biochemical methods. As of recent, metagenomic ancient DNA studies have allowed for the comparison of the fecal and gut microbiomes associated to archaeological specimens from various regions of the world; however the complex dynamics represented in those microbial communities still remain unclear. Theoretically, similar to eukaryote DNA the presence of genes from key microbes or enzymes, as well as the presence of DNA from viruses specific to key organisms, may suggest the ingestion of specific diet components. In this study we demonstrate that ancient virus DNA obtained from coprolites also provides information reconstructing the host’s diet, as inferred from sequences obtained from pre-Columbian coprolites. This depicts a novel and reliable approach to determine new components as well as validate the previously suggested diets of extinct cultures and animals. Furthermore, to our knowledge this represents the first description of the eukaryotic viral diversity found in paleofaeces belonging to pre-Columbian cultures.

Significance Statement

In the midst of an increasing awareness towards the importance of our microbiome, i.e. the microorganisms associated to our bodies, this study is the first molecular characterization of the endogenous retroviruses present in fecal samples of pre-Columbian cultures. Thus, this study gives us a glimpse of the virome and diets of ethnic groups in the Caribbean before the arrival of Europeans. Results obtained from this study complement previous archaeological, botanical and anthropological analyses with molecular data. Furthermore, viral sequences present in coprolites allow for the detection of diet components that cannot be discovered by conventional archaeological studies. Overall, this study has allowed us to further understand the everyday lives of these ancient Caribbean cultures.

Introduction

Coprolites (mummified or fossilized feces) have been shown to document indispensable clues concerning paleodiets, diseases and even ancient cooking habits [13]. In previous studies, the origins and diets of ancient species were unraveled with the detection of key genes in their fecal DNA [4,5]. Notably, recent advances in sequencing technology have allowed for a high-throughput molecular approach to studying paleofeces, and thus exponentially enhance the amount of information obtained from these rare specimens [6]. For example, microbial DNA isolated from coprolites not only provides valuable information on the fecal and gut microbiome of the extinct host, but can also suggest key diet components of these ancient specimens [710].

In contrast, the viral communities present in ancient feces remain largely unknown, possibly because most viral sequences identified in feces belong to modern-day enteric pathogenic viruses; pathogens whose effects on the host’s intestinal system are likely to prevent the formation of coprolites. Despite previous skepticism, recent studies show that DNA from paleofeces can yield further insight in the composition of ancient viral communities (other than pathogenic enteric viruses) within the human gut [11]. However, the information that could be derived from these ancient viromes concerning specific details of the lives of these extinct organisms has barely been explored so far. To this end, the purpose of this study was to determine if diet-related information could be obtained from the viral DNA in coprolites from two pre-Columbian indigenous cultures (i.e. Saladoid and Huecoid) that inhabited Puerto Rico. Previous carbon-dating indicated the sample ages ranged from 1,400 to approximately 540 uncalibrated years before present.

Materials and Methods

Sample description

Coprolites previously obtained from three archaeological deposits in Sorcé Estate, located on the Caribbean island of Vieques were used (Table 1) [8]. Carbon-dating analysis was conducted in two separate laboratories (Teledyne Isotopes in Westwood, NJ; BETA Analytic, Inc. in Miami, FL), using standard methods. Excavations in La Hueca site, Sorcé Estate, were lead by archaeologists Luis Chanlatte Blaik and Yvonne Narganes Storde of the Archaeological Research Center at the University of Puerto Rico, Rio Piedras Campus. The University of Puerto Rico is an academic institution and as such is not required to obtain special permits for academic research and archaeological excavations. Excavations on private property in Sorcé, Vieques were done with the owner’s consent.

thumbnail
Table 1. Description of coprolites analyzed in this study.

https://doi.org/10.1371/journal.pone.0144951.t001

Preventing contamination

Upon excavation, coprolites were individually stored in sterile sample bags until further use. Samples were individually processed inside a class II biosafety cabinet (BSC) used exclusively for ancient DNA analysis. Strategies used for the obtainment of reliable data were implemented [12]. The BSC was routinely cleaned with 70% ethanol and irradiated with UV for 30 min before and after use. Solutions used for DNA extractions were dispensed into single-use aliquots. Previously unused, sterile micropipettes were used; non-disposable equipment were sterilized by autoclave and baked overnight at >100°C. Baked and autoclaved stainless steel utensils were used to separate the inner and outer layers of the coprolites and were sterilized between samples using 70% ethanol. Previously published data comparing the microbial profiles detected in the inner and outer regions of the coprolites was used as an additional control [8].

DNA extraction

First, the exterior layers of the coprolites were aseptically removed in a class II biosafety hood used strictly for aDNA, and subjected to controls to avoid extant DNA contamination. DNA was isolated using the commercial PowerSoil DNA extraction kit (MO BIO Laboratories, Inc.) Coprolite cores were ground and 5g of each sample were hydrated overnight in sterile C1 buffer at 4°C. In the case of the Huecoid culture, an additional coprolite was needed to obtain sufficient material from the interior of the coprolites. DNA yield was then concentrated to 10ul using standard glycogen precipitation and pooled to one composite for each ethnic group (MixS1 and MixH1).

DNA sequencing

The two DNA composites were sequenced in a separate laboratory (MR DNA Research Lab, Shallowater, TX) using a non-targeted metagenomic approach. The library was prepared using Nextera DNA Sample preparation kit (Illumina) following the manufacturer's user guide. The initial concentration of DNA was evaluated using the Qubit® dsDNA HS Assay Kit (Life Technologies). Because of too low DNA concentration for both the samples, whole genome amplification was carried out by using REPLI-g Midi kit (Qiagen) followed by Nextera DNA Sample preparation. The linear amplified DNA was purified using PowerClean DNAClean-Up Kit (MO BIO Laboratories) and the concentration was again evaluated using the Qubit® dsDNA HS Assay Kit (Life Technologies) (S1 Table). Samples were then diluted accordingly to achieve the recommended DNA input of 50ng at a concentration of 2.5ng/uL. Subsequently, the samples underwent simultaneous fragmentation and addition of adapter sequences. These adapters were utilized during a limited-cycle (5 cycles) PCR in which a unique index was added to the sample. Following the library preparation, the final concentration of the library (S1 Table) was measured using the Qubit® dsDNA HS Assay Kit (Life Technologies), and the average library size was determined using Experion (Bio-Rad). The libraries were then pooled in equimolar ratios of 4nM, and 13.5pM of the library pool was sequenced paired end for 500 cycles using the MiSeq system (Illumina).

Data processing and identification

DNA reads were filtered using MG-RAST [13] default parameters for quality scores and read length. Overlapping pair-end reads were merged using default parameters, however non-overlapping reads were also retained to maximize the amount of information obtained despite possible fragmentation due to taphonomic processes. Screened reads were downloaded and submitted to diagrid.org hub server (Purdue Research Foundation, West Lafayette, IN, July 2014). Translated nucleotide BLASTX and taxonomical BLASTN analyses were done using the BLASTer tool [14] against the non-redundant NCBI database (National Center for Biotechnology Information). Cut off values for functional and taxonomical identification included 85% minimum similiarity and E value ≤1e-15. All identified sequences with the terms ‘virus’ and/or ‘phage’ in the description were then pooled using the grepi command in the shell command prompt (Apple OS X) for further evaluation.

Project data

Sequence data discussed in this study is included in S1 Dataset.

Results and Discussion

To determine if eukaryote viral DNA in coprolites could suggest the dietary components ingested by ancient specimens, we isolated and sequenced DNA of coprolites from two Caribbean pre-Columbian cultures (n = 4 Saladoid; n = 5 Huecoid, see Table 1 for more details). Although they are believed to be from South American origins, the Saladoid and Huecoid cultures co-inhabited Sorcé, Vieques for over one thousand years. Unfortunately, the extent to which these ethnic groups influenced modern-day Puerto Rican genetic heritage is unkown.

A brief description of the raw sequence data obtained from metagenomic analysis of pooled samples of each culture is included in Table 2. Similar to previous studies, while sequences belonging to bacteriophages infecting the intestinal and fecal microbiota were the most abundant in our data, a small number of DNA fragments from eukaryote viruses were also detected (S2 Table) [15,16]. Although the bacteriophages detected in this study will be further discussed in another communication, the presence of key bacteriophages already hinted at possible components of these cultures’ diets. For example, Vibrio phages were detected in both cultures; these bacteria can be harbored by a variety of marine organisms, such as crustaceans, that have been proposed by archaeologists as diet components of these cultures. In addition, archaean and protozoan viruses were also detected in lesser proportions. To our surprise, several eukaryote retroviral sequences presented strong evidence of the diets for both cultures (Table 3). For instance, avian pox, as well as grouper iridovirus and frog virus sequences were present in Saladoid feces, confirming their consumption of birds, fish and amphibians. This was further supported by the detection of DNA sequences similar to those described in Xenopus sp. (S3 Table). Avian pox viruses (Poxviridae), are chordopoxviruses that infect domestic and wild game birds including hawks, seagulls, parrots, canaries, pigeons, hawks and fowl, causing cutaneous, respiratory and gastrointestinal lesions in their host [17]. Of these, the fowlpox virus is one of the most studied, its genome fully sequenced and annotated [18]. On the other hand, Iridoviruses (Iridoviridae) only infect invertebrates and poikilothermic vertebrates such as fish, amphibians and reptiles, often causing serious systemic infections and even death [19]. Our results were further supported by the detection of sequences from protozoan parasites infecting similar viral hosts, such as Histomonas meleagridis, which infects a wide range of birds, and Perkinsus marinus, an oyster pathogen (this will be further discussed in another communication). Caribbean indigenous cultures are known to have indulged in various bird and fish species as part of their diets. However, the consumption of amphibians was, to the best of our knowledge, previously unreported.

thumbnail
Table 2. Nucleotide fragments detected via metagenomic analysis of coprolite DNA.

https://doi.org/10.1371/journal.pone.0144951.t002

thumbnail
Table 3. Viral sequences detected as components of the Saladoid and Huecoid diets.

https://doi.org/10.1371/journal.pone.0144951.t003

Overall, our results overlap synergistically with those previously suggested by osteological and botanical findings associated with these archaeological deposits and other pre-Columbian cultures in the Caribbean [27,28]. Notably, the majority of the eukaryote viruses identified in this study were known retroviruses, strongly suggesting the presence of proviruses and/or the lateral acquisition of genes during the host-parasite interaction. This was no surprise as viral “junk” DNA is known to occupy approximately 4.8% of the human genome; in contrast, eukaryote protein-coding genes comprise only 3% of the same genome [29]. Furthermore, the presence and expression of horizontally-acquired genes are now known to be common in host-parasite interactions involving vertebrates and invertebrates [30,31]. Endogenous retroviruses (ERV) in particular are viral gene remnants that seem to have integrated into the host genome after infection [32]. ERVs are passed onto the host’s progeny and may persist for over a hundred million years in the germline; therefore, they serve as an excellent snapshot of a species’ history of infection [33,34]. In addition, within the eukaryote genome ERVs are subject to lower mutation rates as opposed to non-retroviral RNA viruses [35].

We observed DNA sequences from various endogenous eukaryote retroviruses, such as those infecting humans (HERVs), as well as those infecting rodents, flatworms, nematodes, bats and fungi, in coprolites from both cultures (Table 3). Pre-Columbian indigenous cultures in the Caribbean and other regions of the world have previously been reported to eat cooked and possibly lightly cooked meat [20], as well as roasted, desiccated and even raw vegetables [36]; such eating habits could explain the presence of detectable remnant eukaryotic DNA in their feces. Osteological studies suggest that rodents such as Isolobodon portoricensis (hutía) and Heteropsomys insulans (spiny rat) were normal diet components of the Saladoid and Huecoid, among other Caribbean cultures [37]. Furthermore, evolution and colonization studies of the Antilles suggest the presence of a much higher diversity of rodents (among other animals) in the Caribbean islands than what was previously believed [38]. It was therefore no surprise that we detected sequences similar to those of viruses infecting the Murinidae family, some of these rodents’ few relatives sequenced to date (Table 3). This was further supported by the detection of rodent DNA in these samples, particularly similar to Octodon degus (S3 Table). Although bats endemic to the area at that time were often symbolized in their religious artifacts, until now it remained uncertain if these animals, such as Phyllonycteris major for example, were a part of the Saladoid and Huecoid diets. The detection of bat ERVs in feces from both cultures strongly suggests their consumption. Finding oncogenic retroviral sequences in these coprolites opens the possibility of looking at viral pathogens as additional factors that could have aided in the sudden decline in the populations of certain rodents in this island, such as the hutía and the spiny rat, which to date is commonly thought to have become endangered as a result of weather patterns or over-predation. Although natives of most Caribbean islands fed on the “hutía”, it is peculiar that these rodents only went extinct in Vieques and Puerto Rico, but still thrive in nearby islands. Similarly, bones of the Puerto Rican spiny rat are rarely found in archaeological excavations on the island, indicating that these organisms were not commonly hunted by these cultures, however they are also believed extinct.

Ancient diets have been reconstructed using a variety of methods. For example, morphological and chemical analyses of the bacterial residues found in a dinosaur’s coprolite which clearly suggested an herbivorous diet [39]. Similarly, transient DNA in the gut (and feces) has also been studied in order to detect key genes associated to the host diet, such as chloroplast-specific sequences, for example [40,41]. However, these analyses have yet to include the gut virome, an indispensable and highly populated community of the gut biome. This study clearly demonstrates the novel information that may be obtained from the fecal virome and provirome of ancient organisms through the use of a shotgun sequencing analysis as opposed to selective gene-screening methods. For instance, the detection of over forty DNA fragments corresponding to proteins encoded by Capsaspora owczarzaki, a snail symbiont, strongly suggests a previously unknown yet common diet component of the Saladoid culture (data not shown). This will be further discussed in another communication. Also, the presence of retroviral DNA from marmoset New World monkeys may support the hypothesis of ancient organic trade between Caribbean and South American cultures [21]. In addition, as shown here the presence of host-specific viral and proviral DNA can help validate hypotheses generated by other methods and disciplines studying the diets and customs of ancient specimens. It is of high interest that sequences belonging to a European plant virus were found (Morus notabilis); and this finding opens up more questions than it answers. Although highly hypothetical, and perhaps speculative, these sequences may belong to a Caribbean, pre-Columbian virus that was later introduced to Europe capable of infecting the mulberry tree. This demonstrates one limitation of relying solely on DNA sequences to answer complex questions. It remains to be seen if this hypothesis holds for other unexpected sequences found. Notably, no DNA sequences belonging to either Cassava spp., Ipomea spp. and Zea spp. or their viruses were found in these coprolites when using metagenomic sequencing, in spite of them being staples for these cultures. Curiously, our group had previously detected Zea and Ipomea sequences when using 18S rRNA sequencing. This could be due to the DNA amplification steps used in microbial profiling methods, which could presumably target low abundance DNA. This observation also points to the limitations intrinsic to each method, particularly those common in metagenomic analyses without the pre-amplification of DNA. A possible conclusion may be that these foods were thoroughly cooked prior to consumption, thus although consumed in large amounts, the target DNA may be degraded and thus limiting the possibility for detection by metagenomic analyses. Interestingly, DNA of the pea aphid Acyrthosiphon pisum was found in Huecoid feces, suggesting the consumption of these legumes (S3 Table). Some legumes are endemic to the Caribbean islands, others are often found in South America, where the Saladoid and Huecoid cultures are thought to have originated [23,35].

Conclusions

Advances in metagenomics have increased our perception of the highly complex host-microbiome interactions ocurring within a holobiont. Similarly, the effects certain diets may induce on the dynamics of the gut microbiota are still being determined. This study demonstrates the novel approach of demonstrating the usefulness of the virome and provirome in ancient specimens for the reconstruction of the diets of extinct cultures. In fact, the apparent potential and abundance of proviruses in coprolites point to their importance in future molecular paleovirological studies. Our results complement the currently limited molecular data available on the diets of indigenous cultures prior to the arrival and colonization by Europeans [42]. Fan-Ng, et al have also shown that viral sequences found in paleofeces of extinct animals can be used as a means of taking a peek at the past and an opportunity to reconstruct their diet [43]. Our study and previously published results are also slowly changing the apparent dogmata on the half-life of DNA and its resilience over time [44,45]; however, DNA resiliency studies have been done using purified DNA, thus shedding little light on the survival of intracellular DNA during taphonomic processes over millennia [46]. It is clear from our studies, though, that intracellular DNA is in fact more resilient than previously expected.

Supporting Information

S1 Table. DNA and library concentrations.

1Indicates whole genome amplified and purified DNA. Reported size includes sequencing adaptors, measuring approximately 120bp each.

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

(DOCX)

S2 Table. Virus genes detected by translated-nucleotide query of coprolite DNA from both cultures.

1gi| NCBI gene identification number. 2Identification code of the corresponding DNA sequence detected in Huecoid (H-A4LNU…) and Saladoid (S-A4LNU…) coprolite samples. (See S1 Dataset for complete DNA sequences).

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

(DOCX)

S3 Table. Description of eukaryote diet-associated genes detected after conducting a translated-nucleotide query of coprolite DNA from both cultures.

1gi| NCBI gene identification number. 2Identification code of the corresponding DNA sequence detected in Huecoid (H-A4LNU…) and Saladoid (S-A4LNU…) coprolite samples. (See S1 Dataset for complete DNA sequences).

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

(DOCX)

S1 Dataset. Fasta format file containing the trimmed, post-QC ancient DNA sequences described in this study.

These sequences correspond to the genes identified through blast query search.

https://doi.org/10.1371/journal.pone.0144951.s004

(DOCX)

Acknowledgments

We would like to thank MR DNA Research lab (Shallowater, TX) for their reliable sequencing services and excellent troubleshooting skills. Also, many thanks to MoBio for donating the Power-Soil DNA Extraction Kit used in this study.

Author Contributions

Conceived and designed the experiments: JIRP RJC YNS LCB GAT. Performed the experiments: JIRP RJC YNS LCB. Analyzed the data: JIRP RJC. Contributed reagents/materials/analysis tools: RJC YNS LCB GAT. Wrote the paper: JIRP RJC YNS GAT.

References

  1. 1. Reinhard KJ, Confalonieri UE, Herrmann B, Ferreira LF, Araújo A. Recovery of Parasite Remains From Coprolites and Latrines : Aspects of Paleoparasitological Technique. Homo. 1986;37:217–39.
  2. 2. Horne PD. A review of the evidence of human endoparasitism in the pre-Columbian new world through the study of coprolites. J. Archaeol. Sci. 1985;12:299–310.
  3. 3. De Araujo AJ., Ferreira LF, Confalonieri UE, Nunez L, Ribeiro Filho BM. The Finding of Enterobius vermicularis Eggs in Pre-Columbian Human Coprolites. Mem Inst Oswaldo Cruz. 1985;80:141–3. pmid:3915336
  4. 4. Hofreiter M, Poinar HN, Spaulding WG, Bauer K, Martin PS, Possnert G, et al. A molecular analysis of ground sloth diet through the last glaciation. Mol. Ecol. 2000;9:1975–84. pmid:11123610
  5. 5. Poinar HN, Kuch M, Sobolik KD, Barnes I, Stankiewicz AB, Kuder T, et al. A molecular analysis of dietary diversity for three archaic Native Americans. 2001;98:4317–22. pmid:11296282
  6. 6. Tito RY, Knights D, Metcalf J, Obregon-Tito AJ, Cleeland L, Najar F, et al. Insights from characterizing extinct human gut microbiomes. PLoS One [Internet]. 2012;7:e51146. Available: www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3521025&tool=pmcentrez&rendertype=abstract. Accessed 8 October 2013. pmid:23251439
  7. 7. Santiago-Rodriguez TM, Narganes-Storde YM, Chanlatte L, Crespo-Torres E, Toranzos GA, Jimenez-Flores R, et al. Microbial communities in pre-columbian coprolites. PLoS One [Internet]. 2013; 8:e65191. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3673975&tool=pmcentrez&rendertype=abstract. Accessed 7 November 2013. pmid:23755194
  8. 8. Cano RJ, Rivera-Perez JI, Toranzos GA, Santiago-Rodriguez TM, Narganes-Storde YM, Chanlatte-Baik L, et al. Paleomicrobiology: Revealing Fecal Microbiomes of Ancient Indigenous Cultures. White BA, editor. PLoS One [Internet]. 2014; 9:e106833. Available: http://dx.plos.org/10.1371/journal.pone.0106833. Accessed 11 September 2014. pmid:25207979
  9. 9. Appelt S, Armougom F, Bailly M, Robert C, Drancourt M. Polyphasic Analysis of a Middle Ages Coprolite. PLoS One. 2014;9:e88376. pmid:24586319
  10. 10. Cano RJ, Tiefenbrunner F, Ubaldi M, Del Cueto C, Luciani S, Cox T, et al. Sequence analysis of bacterial DNA in the colon and stomach of the Tyrolean Iceman. Am. J. Phys. Anthropol. 2000;112:297–309. pmid:10861348
  11. 11. Appelt S, Fancello L, Le Bailly M, Raoult D, Drancourt M, Desnues C. Viruses in a 14th-century coprolite. Appl. Environ. Microbiol. [Internet]. 2014; 80:2648–55. Available: http://www.ncbi.nlm.nih.gov/pubmed/24509925. Accessed 22 March 2014. pmid:24509925
  12. 12. Drancourt M, Raoult D. Molecular Detection of Past Pathogens. In: Raoult D, Drancourt M, editors. Paleomicrobiology Past Hum. Infect. Berlin Heidelberg: Springer-Verlag; 2008. p. 55–68.
  13. 13. Meyer F, Paarmann D, D’Souza M, Olson R, Glass EM, Kubal M, et al. The metagenomics RAST server—a public resource for the automatic phylogenetic and functional analysis of metagenomes. BMC Bioinformatics [Internet]. 2008; 9:1–8. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2563014&tool=pmcentrez&rendertype=abstract. Accessed 11 July 2014.
  14. 14. Cotton B, Thompson C, Raub B. BLASTer: A hub-based tool for bioinformatics. Las Vegas: WORLDCOMP2013; 2013. p. 1–5. Available: www.diagrid.org.
  15. 15. Minot S, Sinha R, Chen J, Li H, Keilbaugh S a, Wu GD, et al. The human gut virome: inter-individual variation and dynamic response to diet. Genome Res. [Internet]. 2011;21:1616–25. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3202279&tool=pmcentrez&rendertype=abstract. Accessed 14 July 2014. pmid:21880779
  16. 16. Minot S, Bryson A, Chehoud C, Wu GD, Lewis JD, Bushman FD. Rapid evolution of the human gut virome. PNAS. 2013;110:12450–5. pmid:23836644
  17. 17. Bolte AL, Meurer J, Kaleta EF. Avian host spectrum of avipoxviruses. Avian Pathol. 1999;28:415–32.
  18. 18. Afonso CL, Tulman ER, Lu Z, Zsak L, Kutish GF, Rock DL. The genome of fowlpox virus. J. Virol. 2000;74:3815–31. pmid:10729156
  19. 19. Song WJ, Qin QW, Qiu J, Huang CH, Hew CL, Wang F. Functional Genomics Analysis of Singapore Grouper Iridovirus : Complete Sequence Determination and Proteomic Analysis. J. Virol. 2004;78:12576–90. pmid:15507645
  20. 20. de las Casas FB. Apologética Historia Sumaria. Tomo II. Jalisco: Universidad Nacional Autónoma de México; 1967.
  21. 21. Canter Ronald. Yucatan Channel and Trade. FAMSI J. Anc. Am. 2006.
  22. 22. Chanlatte-Baik L, Narganes-Sorde Y. Vieques-Puerto Rico: Asiento de una nueva cultura aborigen antillana. San Juan, Puerto Rico: Arqueologia Puertorriquena; 1983.
  23. 23. Narganes-Sorde Y, Chanlatte-Baik L, Garcia-Padilla A, Hillyer G. La Cultura Saladoide en Puerto Rico. Museo de Historia, Antropologia y Arte Universidad de Puerto Rico; 2002.
  24. 24. Chanlatte Baik L, Narganes-Sorde Y. Cultura La Hueca. 1st ed. San Juan, Puerto Rico: Museo de Historia, Antropologia y Arte: Universidad de Puerto Rico, Recinto de Rio Piedras; 2005.
  25. 25. Laffoon JE, Rodríguez Ramos R, Chanlatte Baik L, Narganes Storde Y, Rodríguez Lopez M, Davies GR, et al. Long-distance exchange in the precolonial Circum-Caribbean: A multi-isotope study of animal tooth pendants from Puerto Rico. J. Anthropol. Archaeol. 2014;35:220–33.
  26. 26. Rodríguez-Ramos R, Pagán-Jiménez JR, Santiago-Blay J, Lambert JB, Craig PR. Some indigenous uses of plants in pre-Columbian Puerto Rico. Life Excit. Biol. [Internet]. 2013;1:83–90. Available: http://essential.metapress.com/openurl.asp?genre=article&id=doi:10.9784/LEB1(1)Rodriguez.09
  27. 27. Pagan-Jimenez J. De antiguos pueblos y culturas botanicas en el Puerto Rico indigena: El archipielago borincano y la llegada de los primeros pobladores agroceramistas. BAR S1687. England: Archaeopress Publishers of British Archaeological Reports; 2007.
  28. 28. Oviedo F de. Historia General y Natural de las Indias. Biblioteca. Madrid: Biblioteca de Autores EspaÑoles Madrid;
  29. 29. Nelson PN, Hooley P, Roden D, Davari Ejtehadi H, Rylance P, Warren P, et al. Human endogenous retroviruses: transposable elements with potential? Clin. Exp. Immunol. [Internet]. 2004;138:1–9. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1809191&tool=pmcentrez&rendertype=abstract. Accessed 17 December 2014. pmid:15373898
  30. 30. Gilbert C, Schaack S, , Brindley PJ, Feschotte C. A role for host-parasite interactions int he horizontal transfer of transposons across phyla. Nature. 2010;464:1347–50. pmid:20428170
  31. 31. Crisp A, Boschetti C, Perry M, Tunnacliffe A, Micklem G. Expression of multiple horizontally acquired genes is a hallmark of both vertebrate and invertebrate genomes. Genome Biol. [Internet]. 2015;16:1–13. Available: http://genomebiology.com/2015/16/1/50.
  32. 32. Weiss RA. On the concept and elucidation of endogenous retroviruses. Phil. Trans. R. Soc. B. 2013;368:1–6.
  33. 33. Boeke J, Stoye J. Retrotransposons, Endogenous Retroviruses, and the Evolution of Retroelements. In: , , , editors. Retroviruses. Cold Spring Harbor (NY): Cold Spring Harbor Laboratory Press.; 1997. p. 1–3.
  34. 34. Lee A, Nolan A, Watson J, Tristem M. Identification of an ancient endogenous retrovirus, predating the divergence of the placental mammals. Phil. Trans. R. Soc. B. 2013;368:1–12.
  35. 35. Nelson PN, Carnegie PR, Martin J, Ejtehadi HD, Hooley P, Roden D, et al. Demystified … Human endogenous retroviruses. J Clin Pathol Mol Pathol. 2003;56:11–8.
  36. 36. Hard RJ, Mauldin RP, Raymond GR. Mano Size, Stable Carbon Isotope Ratios, and Macrobotanical Remains as Multiple Lines of Evidence of Maize Dependence in the American Southwest. J. Archaeol. Method Theory. 1996;3:253–315.
  37. 37. Turvey ST, Oliver JR, Narganes Storde YM, Rye P. Late Holocene extinction of Puerto Rican native land mammals. Biol. Lett. [Internet]. 2007; 3:193–6. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2375922&tool=pmcentrez&rendertype=abstract. Accessed 6 January 2015. pmid:17251123
  38. 38. Brace S, Turvey ST, Weksler M, Hoogland MLP, Barnes I, Turvey ST. Unexpected evolutionary diversity in a recently extinct Caribbean mammal radiation. Proc. R. Soc. B. 2015;282:20142371. pmid:25904660
  39. 39. Hollocher TC, Chin K, Hollocher KT, Kruge MA. Bacterial Residues in Coprolite of Herbivorous Dinasours: Role of Bacteria in Mineralization of Feces. Soc. Sediment. Geol. 2014;16:547–65.
  40. 40. Leach JD, Sobolik KD. High Dietary Intake of Prebiotic Inulin-Type Fructans from Prehistoric Chihuahuan Desert. Br. J. Nutr. 2010;103:1558–61. pmid:20416127
  41. 41. Alcover JA, Perez-Obiol R, Yll E-I, Bover P. The diet of Myotragus balearicus Bate 1909 (Artiodactyla: Caprinae), an extinct bovid from the Balearic Islands: evidence from coprolites. Biol. J. Linn. Soc. [Internet]. 1999;66:57–74. Available: http://doi.wiley.com/10.1111/j.1095-8312.1999.tb01917.x
  42. 42. Sistiaga A, Mallol C, Galván B, Summons RE. The Neanderthal meal: a new perspective using faecal biomarkers. PLoS One [Internet]. 2014;9:e101045. Available: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=4071062&tool=pmcentrez&rendertype=abstract. Accessed 28 November 2014. pmid:24963925
  43. 43. Ng TFF, Chen L-F, Zhou Y, Shapiro B, Stiller M, Heintzman PD, et al. Preservation of viral genomes in 700-y-old caribou feces from a subarctic ice patch. Proc. Natl. Acad. Sci. [Internet]. 2014;111:16842–7. Available: pmid:25349412
  44. 44. Fulton T, Barnett R, Larson G, Rohland N, Wagner S, Shapiro B, et al. Ancient DNA Methods and Protocols. Shapiro B, Hofreiter M, editors. New York: Humana Press; 2014.
  45. 45. Alvarez AJ, Yumet GM, Santiago CL, Toranzos G a. Stability of manipulated plasmid DNA in aquatic environments. Environ. Toxicol. Water Qual. 1996;11:129–35.
  46. 46. Willerslev E, Cooper A. Review Paper: Ancient DNA. Proc. R. Soc. B. 2005;272:3–16. pmid:15875564