Applications of Bayesian Skyline Plots and Approximate Bayesian Computation for Human Demography.

Bayesian methods have been adopted by anthropologists for their utility in resolving complex questions about human history based on genetic data. The main advantages of Bayesian methods include simple model comparison, presenting results as a summary of probability distributions, and the explicit inclusion of prior information into analyses. In the field of anthropological genetics, for example, implementing Bayesian skyline plots and approximate Bayesian computation is becoming ubiquitous as means to analyze genetic data for the purpose of demographic or historic inference. Correspondingly, there is a critical need for better understanding of the underlying assumptions, proper applications, and limitations of these two methods by the larger anthropological community. Here we review Bayesian skyline plots and approximate Bayesian computation as applied to human demography and provide examples of the application of these methods to anthropological research questions. We also review the two core components of Bayesian demographic analysis: the coalescent and Bayesian inference. Our goal is to describe their basic mechanics in an attempt to demystify them.

Ancient individuals from the North American Northwest Coast reveal 10,000 years of regional genetic continuity.

Recent genomic studies of both ancient and modern indigenous people of the Americas have shed light on the demographic processes involved during the first peopling. The Pacific Northwest Coast proves an intriguing focus for these studies because of its association with coastal migration models and genetic ancestral patterns that are difficult to reconcile with modern DNA alone. Here, we report the low-coverage genome sequence of an ancient individual known as "Shuká Káa" ("Man Ahead of Us") recovered from the On Your Knees Cave (OYKC) in southeastern Alaska (archaeological site 49-PET-408). The human remains date to ∼10,300 calendar (cal) y B.P. We also analyze low-coverage genomes of three more recent individuals from the nearby coast of British Columbia dating from ∼6,075 to 1,750 cal y B.P. From the resulting time series of genetic data, we show that the Pacific Northwest Coast exhibits genetic continuity for at least the past 10,300 cal y B.P. We also infer that population structure existed in the late Pleistocene of North America with Shuká Káa on a different ancestral line compared with other North American individuals from the late Pleistocene or early Holocene (i.e., Anzick-1 and Kennewick Man). Despite regional shifts in mtDNA haplogroups, we conclude from individuals sampled through time that people of the northern Northwest Coast belong to an early genetic lineage that may stem from a late Pleistocene coastal migration into the Americas.

Genetic Diversity and Relationships of Tlingit Moieties.

The Tlingit from Southeast Alaska belong to the Northwest Coast cultural tradition, which is defined by regionally shared sociocultural practices. A distinctive feature of Tlingit social organization is the matrilineal exogamous marriage system among clans from two opposite moieties: the Raven/Crow and Eagle/Wolf. Clan and moiety membership are determined by matrilineal descent, and previous genetic studies of Northwest Coast populations have shown a relationship between clan membership and genetic variation of matrilines and patrilines. To further understand this association, in this study mitochondrial DNA sequences from the Tlingit (n = 154) were examined. By comparing mitochondrial DNA with moiety membership information, the authors explore the impact of marriage traditions among the Tlingit with their observable genetic variation. At the genetic level, the results support cultural persistence of Tlingit maternal moiety identity despite the negative impacts of European colonization. This study additionally illustrates the relevance of data derived from Tlingit oral traditions to test hypotheses about population history on the Northwest Coast.

Early human use of anadromous salmon in North America at 11,500 y ago.

Salmon represented a critical resource for prehistoric foragers along the North Pacific Rim, and continue to be economically and culturally important; however, the origins of salmon exploitation remain unresolved. Here we report 11,500-y-old salmon associated with a cooking hearth and human burials from the Upward Sun River Site, near the modern extreme edge of salmon habitat in central Alaska. This represents the earliest known human use of salmon in North America. Ancient DNA analyses establish the species as Oncorhynchus keta (chum salmon), and stable isotope analyses indicate anadromy, suggesting that salmon runs were established by at least the terminal Pleistocene. The early use of this resource has important implications for Paleoindian land use, economy, and expansions into northwest North America.

Patterns of admixture and population structure in native populations of Northwest North America.

The initial contact of European populations with indigenous populations of the Americas produced diverse admixture processes across North, Central, and South America. Recent studies have examined the genetic structure of indigenous populations of Latin America and the Caribbean and their admixed descendants, reporting on the genomic impact of the history of admixture with colonizing populations of European and African ancestry. However, relatively little genomic research has been conducted on admixture in indigenous North American populations. In this study, we analyze genomic data at 475,109 single-nucleotide polymorphisms sampled in indigenous peoples of the Pacific Northwest in British Columbia and Southeast Alaska, populations with a well-documented history of contact with European and Asian traders, fishermen, and contract laborers. We find that the indigenous populations of the Pacific Northwest have higher gene diversity than Latin American indigenous populations. Among the Pacific Northwest populations, interior groups provide more evidence for East Asian admixture, whereas coastal groups have higher levels of European admixture. In contrast with many Latin American indigenous populations, the variance of admixture is high in each of the Pacific Northwest indigenous populations, as expected for recent and ongoing admixture processes. The results reveal some similarities but notable differences between admixture patterns in the Pacific Northwest and those in Latin America, contributing to a more detailed understanding of the genomic consequences of European colonization events throughout the Americas.

DNA analysis of ancient dogs of the Americas: identifying possible founding haplotypes and reconstructing population histories.

As dogs have traveled with humans to every continent, they can potentially serve as an excellent proxy when studying human migration history. Past genetic studies into the origins of Native American dogs have used portions of the hypervariable region (HVR) of mitochondrial DNA (mtDNA) to indicate that prior to European contact the dogs of Native Americans originated in Eurasia. In this study, we summarize past DNA studies of both humans and dogs to discuss their population histories in the Americas. We then sequenced a portion of the mtDNA HVR of 42 pre-Columbian dogs from three sites located in Illinois, coastal British Columbia, and Colorado, and identify four novel dog mtDNA haplotypes. Next, we analyzed a dataset comprised of all available ancient dog sequences from the Americas to infer the pre-Columbian population history of dogs in the Americas. Interestingly, we found low levels of genetic diversity for some populations consistent with the possibility of deliberate breeding practices. Furthermore, we identified multiple putative founding haplotypes in addition to dog haplotypes that closely resemble those of wolves, suggesting admixture with North American wolves or perhaps a second domestication of canids in the Americas. Notably, initial effective population size estimates suggest at least 1000 female dogs likely existed in the Americas at the time of the first known canid burial, and that population size increased gradually over time before stabilizing roughly 1200 years before present.

Origins of an Unmarked Georgia Cemetery Using Ancient DNA Analysis.

Determining the origins of those buried within undocumented cemeteries is of incredible importance to historical archaeologists and, in many cases, the nearby communities. In the case of Avondale Burial Place, a cemetery in Bibb County, Georgia, in use from 1820 to 1950, all written documentation of those interred within it has been lost. Osteological and archaeological evidence alone could not describe, with confidence, the ancestral origins of the 101 individuals buried there. In the present study, we used ancient DNA extraction methods in well-preserved skeletal fragments from 20 individuals buried in Avondale Burial Place to investigate the origins of the cemetery. Through examination of hypervariable region I (HVR1) in the mitochondrial genome (mtDNA), we determined haplotypes for all 20 of these individuals. Eighteen of these individuals belong to the L or U haplogroups, suggesting that Avondale Burial Place was most likely used primarily as a resting place for African Americans. After the surrounding Bibb County community expressed interest in investigating potential ancestral relationships to those within the cemetery, eight potential descendants provided saliva to obtain mtDNA HVR1 information. Three individuals from Avondale Burial Place matched three individuals with oral history ties to the cemetery. Using the online tool EMPOP, we calculated the likelihood of these exact matches occurring by chance alone (< 1%). The present findings exhibit the importance of genetic analysis of cemetery origins when archaeological and osteological data are inconclusive for estimating ancestry of anonymous historical individuals.

Ancient marine hunter-gatherers from Patagonia and Tierra Del Fuego: Diversity and differentiation using uniparentally inherited genetic markers.

OBJECTIVES: The human population history from Patagonia and Tierra del Fuego has been of great interest in the context of the American peopling. Different sources of evidence have contributed to the characterization of the local populations, but some main questions about their history remain unsolved. Among the native populations, two marine hunter-gatherers groups inhabited the Patagonian channels below the 478S: Kawéskar and Yámana. Regardless of their geographical proximity and cultural resemblance, their languages were mutually unintelligible. In this study we aim to evaluate the genetic diversity of uniparental genetic markers in both groups and to test if there is a high genetic differentiation between them, mirroring their linguistic differences. MATERIAL AND METHODS: Ancient DNA was extracted from 37 samples from both populations. We compared their genetic variability of their mitochondrial lineages and Y-STR as well as with other modern native populations from the area and further north. RESULTS AND DISCUSSION: We observed an important differentiation in their maternal lineages: while Kawéskar shows a high frequency of D (80%), Yámana shows a high frequency of C (90%). The analysis of paternal lineages reveals the presence of only Q1a2a1a1 and little variation was found between individuals. Both groups show very low levels of genetic diversity compared with modern populations. We also notice shared and unique mitochondrial DNA variants between modern and ancient samples of Kawéskar and Yámana.

How much DNA is lost? Measuring DNA loss of short-tandem-repeat length fragments targeted by the PowerPlex 16® system using the Qiagen MinElute Purification Kit.

The success in recovering genetic profiles from aged and degraded biological samples is diminished by fundamental aspects of DNA extraction, as well as its long-term preservation, that are not well understood. While numerous studies have been conducted to determine whether one extraction method was superior to others, nearly all of them were initiated with no knowledge of the actual starting DNA quantity in the samples prior to extraction, so they ultimately compared the outcome of all methods relative to the best. Using quantitative PCR to estimate the copy count of synthetic standards before (i.e., "copies in") and after (i.e., "copies out") purification by the Qiagen MinElute PCR Purification Kit, we documented DNA loss within a pool of 16 different-sized fragments ranging from 106 to 409 bp in length, corresponding to those targeted by the PowerPlex 16 System (Promega, Madison, WI). Across all standards from 10(4) to 10(7) copies/µL, loss averaged between 21.75% and 60.56% (mean, 39.03%), which is not congruent with Qiagen's claim that 80% of 70 bp to 4 kb fragments are retained using this product (i.e., 20% loss). Our study also found no clear relationship either between DNA strand length and retention or between starting copy number and retention. This suggests that there is no molecule bias across the MinElute column membrane and highlights the need for manufacturers to clearly and accurately describe on what their claims are based, and should also encourage researchers to document DNA retention efficiencies of their own methods and protocols. Understanding how and where to reduce loss of molecules during extraction and purification will serve to generate clearer and more accurate data, which will enhance the utility of ancient and low-copy-number DNA as a tool for closing forensic cases or in reconstructing the evolutionary history of humans and other organisms.

Exploring prehistory in the North American southwest with mitochondrial DNA diversity exhibited by Yumans and Athapaskans.

A recent study of mitochondrial DNA variation in Native American populations from the American Southwest detected signatures of a population expansion of subhaplogroup B2a, dated to 2,105 years before present (99.5% confidence interval, 1,273-3,773 YBP), following the introduction and intensification of maize agriculture in the region. Only one Yuman group and no Athapaskan speakers were analyzed in previous studies. Here we report mtDNA haplogroup and hypervariable region (HVR I, and II) sequence data from 263 extant Yuman speakers, representing the major branches of the Yuman language family, in addition to the Western Apache (Athapaskan) to further investigate the demographic context and geographic extent of this expansion. Data presented indicate that the expansion of B2a is only slightly older [2,410 YBP (99.5% CI: 1,458-4,320 YBP)] than previously estimated and not significantly. Despite large confidence intervals there are implications for the origin and expansion of the Yuman language family. Cultural transformations due to the inundation and draining of Lake Cahuilla may explain in part the frequencies of this lineage among the Kumeyaay and other Yuman and Takic groups in Southern California. This may have been the result of group fissions and fusions followed by migration and interaction that included expanded trade networks and intermarriage among Yuman speakers. In addition, a series of in-situ genetic bottlenecks is proposed to have occurred among the Western Apache leading to increasing homogeneity within haplogroup A, culminating in an admixture event with the Yavapai.

Brief communication: Evolution of a specific O allele (O1vG542A) supports unique ancestry of Native Americans.

In this study, we explore the geographic and temporal distribution of a unique variant of the O blood group allele called O1v(G542A) , which has been shown to be shared among Native Americans but is rare in other populations. O1v(G542A) was previously reported in Native American populations in Mesoamerica and South America, and has been proposed as an ancestry informative marker. We investigated whether this allele is also found in the Tlingit and Haida, two contemporary indigenous populations from Alaska, and a pre-Columbian population from California. If O1v(G542A) is present in Na-Dene speakers (i.e., Tlingits), it would indicate that Na-Dene speaking groups share close ancestry with other Native American groups and support a Beringian origin of the allele, consistent with the Beringian Incubation Model. If O1v(G542A) is found in pre-Columbian populations, it would further support a Beringian origin of the allele, rather than a more recent introduction of the allele into the Americas via gene flow from one or more populations which have admixed with Native Americans over the past five centuries. We identified this allele in one Na-Dene population at a frequency of 0.11, and one ancient California population at a frequency of 0.20. Our results support a Beringian origin of O1v(G542A) , which is distributed today among all Native American groups that have been genotyped in appreciable numbers at this locus. This result is consistent with the hypothesis that Na-Dene and other Native American populations primarily derive their ancestry from a single source population.

Ancient mitochondrial DNA analysis reveals complexity of indigenous North American turkey domestication.

Although the cultural and nutritive importance of the turkey (Meleagris gallopavo) to precontact Native Americans and contemporary people worldwide is clear, little is known about the domestication of this bird compared to other domesticates. Mitochondrial DNA analysis of 149 turkey bones and 29 coprolites from 38 archaeological sites (200 BC-AD 1800) reveals a unique domesticated breed in the precontact Southwestern United States. Phylogeographic analyses indicate that this domestic breed originated from outside the region, but rules out the South Mexican domestic turkey (Meleagris gallopavo gallopavo) as a progenitor. A strong genetic bottleneck within the Southwest turkeys also reflects intensive human selection and breeding. This study points to at least two occurrences of turkey domestication in precontact North America and illuminates the intensity and sophistication of New World animal breeding practices.

Evaluating the Farming/Language Dispersal Hypothesis with genetic variation exhibited by populations in the Southwest and Mesoamerica.

The Farming/Language Dispersal Hypothesis posits that prehistoric population expansions, precipitated by the innovation or early adoption of agriculture, played an important role in the uneven distribution of language families recorded across the world. In this case, the most widely spread language families today came to be distributed at the expense of those that have more restricted distributions. In the Americas, Uto-Aztecan is one such language family that may have been spread across Mesoamerica and the American Southwest by ancient farmers. We evaluated this hypothesis with a large-scale study of mitochondrial DNA (mtDNA) and Y-chromosomal DNA variation in indigenous populations from these regions. Partial correlation coefficients, determined with Mantel tests, show that Y-chromosome variation in indigenous populations from the American Southwest and Mesoamerica correlates significantly with linguistic distances (r = 0.33-0.384; P < 0.02), whereas mtDNA diversity correlates significantly with only geographic distance (r = 0.619; P = 0.002). The lack of correlation between mtDNA and Y-chromosome diversity is consistent with differing population histories of males and females in these regions. Although unlikely, if groups of Uto-Aztecan speakers were responsible for the northward spread of agriculture and their languages from Mesoamerica to the Southwest, this migration was possibly biased to males. However, a recent in situ population expansion within the American Southwest (2,105 years before present; 99.5% confidence interval = 1,273-3,773 YBP), one that probably followed the introduction and intensification of maize agriculture in the region, may have blurred ancient mtDNA patterns, which might otherwise have revealed a closer genetic relationship between females in the Southwest and Mesoamerica.

Freshwater and anadromous fishing in Ice Age Beringia.

While freshwater and anadromous fish have been critical economic resources for late prehistoric and modern Native Americans, the origin and development of fishing is not well understood. We document the earliest known human use of freshwater and anadromous fish in North America by 13,000 and 11,800 years ago, respectively, from primary anthropogenic contexts in central Alaska (eastern Beringia). Fish use appears conditioned by broad climatic factors, as all occurrences but one are within the Younger Dryas chronozone. Earlier Bølling-Allerød and later early Holocene components, while exhibiting similar organic preservation, did not yield evidence of fishing, suggesting that this was a response to changing environmental factors, perhaps reductions in higher ranked resources such as large terrestrial mammals. Late Pleistocene and recent Indigenous peoples harvested similar fish taxa in the region (salmon, burbot, whitefish, and pike). We characterize late Pleistocene fishing in interior Beringia as an important element of broad-spectrum foraging rather than the intensive communal fishing and storage common among recent peoples.

The genetic impact of Aztec imperialism: ancient mitochondrial DNA evidence from Xaltocan, Mexico.

In AD 1428, the city-states of Tenochtitlan, Texcoco, and Tlacopan formed the Triple Alliance, laying the foundations of the Aztec empire. Although it is well documented that the Aztecs annexed numerous polities in the Basin of Mexico over the following years, the demographic consequences of this expansion remain unclear. At the city-state capital of Xaltocan, 16th century documents suggest that the site's conquest and subsequent incorporation into the Aztec empire led to a replacement of the original Otomí population, whereas archaeological evidence suggests that some of the original population may have remained at the town under Aztec rule. To help address questions about Xaltocan's demographic history during this period, we analyzed ancient DNA from 25 individuals recovered from three houses rebuilt over time and occupied between AD 1240 and 1521. These individuals were divided into two temporal groups that predate and postdate the site's conquest. We determined the mitochondrial DNA haplogroup of each individual and identified haplotypes based on 372 base pair sequences of first hypervariable region. Our results indicate that the residents of these houses before and after the Aztec conquest have distinct haplotypes that are not closely related, and the mitochondrial compositions of the temporal groups are statistically different. Altogether, these results suggest that the matrilines present in the households were replaced following the Aztec conquest. This study therefore indicates that the Aztec expansion may have been associated with significant demographic and genetic changes within Xaltocan.

Nondestructive sampling of human skeletal remains yields ancient nuclear and mitochondrial DNA.

Museum curators and living communities are sometimes reluctant to permit ancient DNA (aDNA) studies of human skeletal remains because the extraction of aDNA usually requires the destruction of at least some skeletal material. Whether these views stem from a desire to conserve precious materials or an objection to destroying ancestral remains, they limit the potential of aDNA research. To help address concerns about destructive analysis and to minimize damage to valuable specimens, we describe a nondestructive method for extracting DNA from ancient human remains. This method can be used with both teeth and bone, but it preserves the structural integrity of teeth much more effectively than that of bone. Using this method, we demonstrate that it is possible to extract both mitochondrial and nuclear DNA from human remains dating between 300 BC and 1600 AD. Importantly, the method does not expose the remains to hazardous chemicals, allowing them to be safely returned to curators, custodians, and/or owners of the samples. We successfully amplified mitochondrial DNA from 90% of the individuals tested, and we were able to analyze 1-9 nuclear loci in 70% of individuals. We also show that repeated nondestructive extractions from the same tooth can yield amplifiable mitochondrial and nuclear DNA. The high success rate of this method and its ability to yield DNA from samples spanning a wide geographic and temporal range without destroying the structural integrity of the sampled material may make possible the genetic study of skeletal collections that are not available for destructive analysis.

Ancient DNA analysis of archaeological specimens extends Chinook salmon's known historic range to San Francisco Bay's tributaries and southernmost watershed.

Understanding a species' historic range guides contemporary management and habitat restoration. Chinook salmon (Oncorhynchus tshawytscha) are an important commercial and recreational gamefish, but nine Chinook subspecies are federally threatened or endangered due to anthropogenic impacts. Several San Francisco Bay Area streams and rivers currently host spawning Chinook populations, but government agencies consider these non-native hatchery strays. Through the morphology-based analysis of 17,288 fish specimens excavated from Native American middens at Mission Santa Clara (CA-SCL-30H), Santa Clara County, circa 1781-1834 CE, 88 salmonid vertebrae were identified. Ancient DNA sequencing identified three separate individuals as Chinook salmon and the remainder as steelhead/rainbow trout (Oncorhynchus mykiss). These findings comprise the first physical evidence of the nativity of salmon to the Guadalupe River in San Jose, California, extending their documented historic range to include San Francisco Bay's southernmost tributary watershed.

Mitochondrial DNA of protohistoric remains of an Arikara population from South Dakota: implications for the macro-Siouan language hypothesis.

Mitochondrial DNA (mtDNA) was extracted from skeletal remains excavated from three Arikara sites in South Dakota occupied between AD 1600 and 1832. The diagnostic markers of four mtDNA haplogroups to which most Native Americans belong (A, B, C, and D) were successfully identified in the extracts of 55 (87%) of the 63 samples studied. The frequencies of the four haplogroups were 42%, 29%, 22%, and 7%, respectively, and principal coordinates analysis and Fisher's exact tests were conducted to compare these haplogroup frequencies with those from other populations. Both analyses showed closer similarity among the Mohawk, Arikara, and Sioux populations than between any of these three and any other of the comparison populations. Portions of the first hypervariable segment (HVSI) of the mitochondrial genome were successfully amplified and sequenced for 42 of these 55 samples, and haplotype networks were constructed for two of the four haplogroups. The sharing of highly derived lineages suggests that some recent admixture of the Arikara with Algonquian-speaking and Siouan-speaking groups has occurred. The Arikara shared more ancient lineages with both Siouan and Cherokee populations than with any other population, consistent with the Macro-Siouan language hypothesis that Iroquoian, Siouan, and Caddoan languages share a relatively recent common ancestry.

The earliest farmers of northwest China exploited grain-fed pheasants not chickens.

Though chickens (Gallus gallus domesticus) are globally ubiquitous today, the timing, location, and manner of their domestication is contentious. Until recently, archaeologists placed the origin of the domestic chicken in northern China, perhaps as early as 8,000 years ago. Such evidence however complicates our understanding of how the chicken was domesticated because its wild progenitor - the red jungle fowl (Gallus gallus) - lives in tropical ecosystems and does not exist in northern China today or in the recent past. Increasingly, multiple lines of evidence suggest that many of the archaeological bird remains underlying this northern origins hypothesis have been misidentified. Here we analyze the mitochondrial DNA of some of the earliest purported chickens from the Dadiwan site in northern China and conclude that they are pheasants (Phasianus colchicus). Curiously, stable isotope values from the same birds reveal that their diet was heavy in agricultural products (namely millet), meaning that they lived adjacent to or among some of the earliest farming communities in East Asia. We suggest that the exploitation of these baited birds was an important adaptation for early farmers in China's arid north, and that management practices like these likely played a role in the domestication of animals - including the chicken - in similar contexts throughout the region.

Subduing the influence of PCR inhibitors on amplifying aged, degraded, and low copy number DNA: PCR enhancer cocktail-p and rescue PCR.

PCR inhibitors are a formidable problem to the study of aged, degraded, and/or low copy number DNA. As a result, there is a need to find alternate methods that ameliorate the efficacy of PCR. In this study, we attempted to use genetic methods to identify the species of salmonid (Oncorhynchus spp.) remains recovered from archaeological sites along the Feather River located in northern California, United States. In the process of doing so, we compared the efficacy of a PCR enhancer cocktail called "PEC-P" and a reagent rich PCR recipe called "rescue PCR" over standard PCR. Across all treatments (full concentration and 1:10 dilute eluates subjected to standard PCR, PEC-P, and rescue PCR) species identification was possible for 74 of 93 archaeological fish specimens (79.6%). Overall, six of the 93 samples (6.5%) consistently yielded species identification across all treatments. The species of ten specimens (10.8%) were uniquely identified from amplicons produced with either PEC-P or rescue PCR or both. Notably, the species of seven samples (7.5%) were uniquely identified with standard PCR over the alternative treatments. Considering both full concentration and 1:10 dilute eluates (N = 186), standard PCR performed as well as PEC-P (p = 0.1451) and rescue (p = 0.6753). Yet, considering results from full concentration eluates alone (N = 93), PEC-P (60.2%) outperformed both standard PCR (44.1%; p = 0.0277) and rescue PCR (40.9%; p = 0.0046). Stochasticity observed in our study cautions us against choosing a "best" performing method of those explored here and suggests their respective potentials to improve success may be sample dependent. When working with samples compromised by PCR inhibitors, it is useful to have alternative methodologies for subduing the problem. Both PEC-P and rescue PCR represent useful alternative methods for the study of aged, degraded, and/or low copy number DNA samples compromised by PCR inhibitors.