Introduction

Pigs were first domesticated in the Near East around 8500 BC and subsequently brought into Europe by agriculturalists1. Ancient mitochondrial DNA (mtDNA) studies further indicate that, by 4500 BC, domesticated pigs bearing Near Eastern haplotypes appeared in northern Europe2. Soon thereafter, agriculturalists are thought to have incorporated local wild boars into their domesticated swine herds1,2. On the basis of a relatively small and geographically diffuse set of mtDNA data, it has been hypothesized that the Near Eastern mtDNA lineages were rapidly replaced by European haplotypes, apparently within 500 years or less, after the first introduction of Near Eastern domesticates into the region1.

Here, we set out to derive a more precise geographic and temporal framework for the early management of indigenous wild boar populations in northern Europe, a process which ultimately led to domesticated European lineages. We examined ancient mtDNA haplotype data from a large sample of 87 domestic and 29 wild Neolithic Sus specimens (Fig. 1). The study targeted an 80-bp highly informative fragment of the mitochondrial control region that provides ample information about maternal ancestry1,3,4,5. The presence of Near Eastern haplotypes Y1 and Y2 in pigs located in Europe has been previously interpreted as a marker of their domesticated status1, demonstrating the role of humans in facilitating the translocation of suids from the Near East into Europe. Reflecting the difficulty in resolving the domesticated status of animals by lineage alone in regions inhabited by both domesticates and their wild progenitors, Sus with the European haplotype A or C could represent either wild or domestic animals2. In order to better resolve the source and timing of pig management processes, we developed a stochastic model to trace the mtDNA haplotypes from nascent domestic pig populations of northern Europe of 5000 BC to 3000 BC back to their likely origin. Our study is the first to use such a quantitative model to ascertain the origin of prehistoric mtDNA haplotypes of northern domestic pigs.

Figure 1: Location of archaeological Sus samples.
figure 1

The large circle represents nine geographically close archaeological sites. The dashed line separates the northern and the southern group of animals.

Results

Study design and aim

In our sample of 116 Neolithic Sus specimens, we observed the four known mtDNA haplotypes Y1, Y2, A and C1 (Tables 1 and 2). Interestingly, Near Eastern haplotype Y1 appeared as late as 4000 to 3000 BC in northern domestic pigs.

Table 1 mtDNA haplotype distribution in ancient domestic pigs.
Table 2 mtDNA haplotype distribution in ancient wild boars.

For the present study, pig remains recovered from sites located in the northern parts of Germany and the Netherlands were classified as belonging to the β€˜northern group’ (n = 73), irrespective of haplotype (Fig. 1, Table S1). The northern region was home to late Mesolithic hunter-gatherers who eventually turned towards the use of domesticated plants and animals during the Neolithic (ca 4100 BC). Haplotype frequencies obtained from northern group pigs were compared to those among specimens recovered from southern Poland and central/southern Germany (β€˜southern group’, n = 43; Fig. 1), a region where Neolithic farmers had been established since ~5500 BC. The wild or domestic status of each specimen was determined based upon (i) its mtDNA haplotype (with Y1 or Y2 indicating domestic status) and (ii) standard metrical and non-metrical analyses (cp. Methods).

Animals were grouped further according to the time period (t1: 5500–5000 BC, t2: 5000–4000 BC, t3: 4000–3000 BC) and region (north vs. south) they had been living in (Tables 1 and 2). The three time periods and the two geographical regions were characterized by distinct subsistence strategies, landscape interactions and cultural material used. The first time period (t1) was largely defined by the presence of Mesolithic hunter-gatherer-fishers in the north and the first farmers in the south (5500–5000 BC). The second period (t2) (5000–4000 cal BC) was described in the north by the transition to early farmers and in the south by the super-regional Linienbandkeramik culture. The third period (t3) (4000–3000 cal BC) was defined by the emergence of more regional agriculturalist Funnel Beaker and Globular Amphora culture. The mtDNA haplotypes of the t2 sample of wild boars from the north (n = 29) were used as a reference for the haplotype distribution of wild boars in periods t1 and t2 (Table 2). In our stochastic model, we investigated time periods t2 and t3 separately. We determined the probability that the mtDNA of a northern domestic pig derived from an animal of certain status (domestic or wild) living in a certain region (north or south) during the previous time period (i.e., t1 or t2). Based on the occurrence of European haplotypes A and C in southern domestic pigs, it is clear that wild boars had been continuously integrated into the domestic livestock in the south1,2. Therefore, the group of southern domestic pigs consists not only of domestic pigs originally descended from domestic pigs from the Near East, but also of previously domesticated southern wild boars.

Origin of mtDNA haplotypes of northern domestic pigs living 5000–4000 BC (t2)

For time period t1 (5500–5000 BC), no domestic pigs have so far been identified in northern European contexts6,7. Therefore, the mtDNA of a northern t2 domestic pig can only originate from either a southern t1 domestic pig or a northern t1 wild boar (Fig. 2A). A single parameter is sufficient to describe this model, namely the probability pNG that the mtDNA of a northern t2 domestic pig derived from a northern t1 wild boar. Our data yielded a maximum likelihood estimate of 0.20 for pNG for two haplotypes (combining A and C on the one hand, and Y1 and Y2 on the other) and of 0 when all four haplotypes were considered. The 95% confidence interval for pNG was [0, 0.83] for two haplotypes and [0, 0.28] for four haplotypes (Table 3, Figures S1A and B). The two-haplotype confidence interval is considerably larger because the European haplotypes A and C were combined, which masks the fact that in southern t1 domestic pigs, haplotype A is more frequent than haplotype C (A: 64% vs. C: 36% among the European haplotypes) whereas in northern t1 wild boars haplotype C is almost five times more frequent than A (A: 17% vs. C:83%). The four haplotype confidence interval takes the separate information about A and C into account and therefore yields a more precise estimate.

Figure 2: Stochastic models of the origin of mtDNA haplotypes in northern domestic pigs.
figure 2

NG, northern group (green); SG, southern group (purple); red square, domestic pig; blue square, wild boar; pNG, probability that a northern domestic pig descended from a northern pig of the previous time period; pDP, probability that a northern predecessor was domestic; (A): period 5000–4000 BC (t2) (B): period 4000–3000 BC (t3).

Table 3 Parameter estimates for the origin of northern domestic pigs.

Origin of mtDNA haplotypes of northern domestic pigs living 4000–3000 BC (t3)

Since domestic pigs were already present in northern Europe in time period t2, mtDNA haplotypes of northern t3 domestic pigs have three possible origins, namely southern t2 domestic pigs, northern t2 domestic pigs or northern t2 wild boars (Fig. 2B). In addition to parameter pNG, i.e. the probability that the haplotype was derived from a northern rather than a southern t2 animal, a second parameter pDP is required to distinguish between northern domestic and northern wild predecessors (pDP equals the probability that a northern predecessor was domestic rather than wild). For two haplotypes, maximum likelihood estimation did not yield a unique solution in our data. Instead, possible values of pNG and pDP were found to lie on a straight line (Fig. 3A). Consequently, pNG was estimated to lie between 0.78 and 1 whereas pDP was estimated to lie between 0 and 0.57 (Table 3). When all four haplotypes were considered, unique maximum likelihood estimates were obtained as 1 for pNG and 0.57 for pDP (Table 3, Fig. 3B).

Figure 3: Likelihood surfaces for parameters pNG and pDP in time period t3 (4000–3000 BC).
figure 3

pNG, probability that a northern t3 domestic pig descended from a northern t2 pig; pDP, probability that a northern predecessor (from time period t2) was domestic; for the definition of the likelihood function, see Methods, section Parameter estimation, Equation (3); (A): likelihood surface for two haplotypes (A and C vs Y1 andY2) (B): likelihood surface for four haplotypes (A, C, Y1, Y2)

Discussion

Our results ensued from 116 Neolithic specimens, a relatively large number for an ancient DNA study, thus enabling us to perform sensible statistical modelling. To the best of our knowledge, our study has been the first to develop a quantitative model to assess the prehistoric origin of mtDNA haplotypes of northern domestic pigs. However, for precise estimation of the model parameters, the study sample size may have been only moderately sufficient, as is evidenced by the relatively large confidence intervals obtained. Nevertheless, a clear trend is apparent. For early period t2, our analysis suggests none or only little influx of wild boar from the north. Thus, whereas predominant maternal descent from southern domestic pigs is likely, it cannot be excluded that a minor proportion of pigs received their mtDNA from animals from the northern region. The confidence interval for pNG (the probability that the mtDNA of a northern t2 domestic pig derived from a northern t1 wild boar) was smaller in the four-haplotype model, not least owing to the larger amount of information that was taken into account, so that the parameter estimates of the four-haplotype model can be assumed to be more accurate than those of the two-haplotype model. Qualitatively, however, the two- and four-haplotype models agreed well for both time periods, thereby increasing confidence in our results.

Previous zooarchaeological research suggests that domesticated pigs were first exploited in northern Europe at around 4100 BC8. According to a recent aDNA study, however, domesticated Near Eastern lineage animals may have been present in the region several centuries earlier, a result derived from a directly radiocarbon-dated specimen with haplotype Y1 (two sigma range: BC 4720–4582)2. The mtDNA results presented here further suggest that, during the early time period t2 (5000–4000 BC), female domestic pigs were introduced almost exclusively from the south. In the following time period t3 (4000–3000 BC), by contrast, imports from the south contributed little to the northern mtDNA genetic pool and, consequently, the mtDNA of northern European lineage domestic pigs should have originated mainly from local animals, a substantial proportion (between 43% and 100%) of which were wild boars. Our model highlights a clear shift in suid management strategies after 4000 BC, a pattern in agreement with the model generated by Larson et al.1, which drew from a much smaller and geographically more widespread sample. Based on the occurrence of Near Eastern haplotype Y1 as late as t3 (between 4000 and 3000 BC; Tables 1 and 2) it seems that, at least in the north, the replacement of Near Eastern mtDNA haplotypes by European ones could have taken longer (>500 years) than previously suggested1. Equally intriguing, our study revealed a high level of interbreeding between the early domestic northern populations and local wild boars. Notably, significant post-domestication gene flow from wild animals into managed swineherds is further corroborated by a recent genomic study on modern pigs9.

Why did early agriculturalists in northern Europe change their suid management strategy after 4000 BC? It may be that the availability of domesticated Near Eastern pigs in the earliest Neolithic settlements located in northern Europe was initially limited, reflecting their location at the end of the Neolithic β€˜supply chain’. The reliability of Near Eastern-derived livestock was perhaps also unstable, particularly if pig husbandry strategies applied to Near Eastern domesticates initially failed to adjust for northern European environments and dietary conditions. The overall number of domestic pigs carried to the north by the first Neolithic agriculturalists along with their livestock may have been relatively small. As a case in point, domesticated sheep and goats, core constituents of the Near Eastern Neolithic livestock package, were also exploited at low intensities only in the earliest Neolithic settlements of the southwestern Baltic region, with their remains generally comprising no more than 10% of the faunal assemblage10,11 [Makarewicz, unpublished data]. A small domestic source population would have been susceptible to isolation, disease and over-harvesting which would have made it difficult to maintain a stable and viable swine herd. In addition, the material cultural record indicates that the contact between northern and southern European cultures was only sporadic at the time8,12. If this was the case, then the strategy of interbreeding locally abundant female wild boars with domestic males would have facilitated an increase in swine herd size without having to rely upon a limited and uncertain supply of Near Eastern domesticated pigs from the south. Moreover, local stock are usually well adapted to their specific environmental conditions, for instance, by showing greater resistance to endemic pathogens13. Out-crossing therefore may have produced more resilient, fertile and larger offspring12,13. Another explanation for the dominance of local European mtDNA genomes in later pig specimens could be that the corresponding genomes conferred a selective advantage in the north, for example, in terms of energy metabolism that would eventually lead to their fixation in European domesticates.

It has been suggested before that interbreeding between pigs and wild boars was mainly unintentional and resulted from chance encounters due to escaped feral domesticates, loose swine management systems and mobile swine herding9,12. Loose pig husbandry is thought to have been characteristic of European agricultural systems for millennia14. However, this scenario would not explain the introgression of European wild sow mtDNA haplotypes (i.e. of female animals) into expanding Neolithic swineherds; mating between domesticated female sows and European wild boars would not have changed the mtDNA pool of domestic herds. The latter instead would have required the incorporation of additional wild females, which could, for instance, be achieved by capture of wild female piglets as part of an active herd-building strategy. Therefore, a combination of loose management and intentional integration of female wild boars seems to have occurred.

The results of our quantitative analyses strongly suggest that pig management was an ever-evolving process which depended heavily upon interbreeding domesticated animals with local wild boar, in particular wild sows. This insight demonstrates how a change in animal handling ~6000 years ago may have influenced livestock composition up to the present day. What is more, our mathematical model presented here can also be adapted to other proxies (i.e. animals) and haploid DNA markers as well as scenarios and periods.

Methods

Sample specification and determination of the wild or domestic status of samples

The domestication status of the sampled Sus remains was determined in a two-stage process. Since Near Eastern mtDNA haplotypes Y1 and Y2 classify a pig as domestic with very high probability, samples with one of these two haplotypes were assigned domestic status (n = 25). Since pigs with European haplotype A or C can be either domestic or wild, we had to rely upon standard metrical and non-metrical morphological analyses for the classification of Sus with these haplotypes (n = 91).

Metrical data retrieved from the appendicular skeleton and from teeth, in particular the third mandibular molar, are commonly used to evaluate change in animal body size associated with intensive husbandry. Large bodied animals are generally associated with the wild condition. Size diminution over time is thought to reflect the introduction of anthropogenic selective pressures and the onset of domesticated phenotypes with small-bodied animals reflecting their domesticated status15. However, especially in the early stages of the domestication process when progenitor species may have been only loosely managed, the relationship between body size and animal β€˜domestication’ status is less clear-cut6,7. Status determination taking into account body size alone can therefore be error-prone. For this study, pig bone specimens were classified as belonging to β€˜domesticated, or β€˜wild’ according to criteria established independently by numerous zooarchaeologists over the course of several decades. For the majority of pig bone specimens obtained from the northern region, domesticated status was established using pig bone and tooth specimens from the medieval site of Hedeby located in Schleswig-Holstein, northern Germany. For each skeletal element, all Mesolithic or Neolithic bone specimens smaller than Hedeby specimens were classified as β€˜domesticated’ (SchmΓΆlcke, personal communication). Bone and tooth specimens within the range of modern European wild boar were classified as β€˜wild’, while bone and tooth specimens larger than those from Hedeby but smaller than modern wild boar were classified as β€˜indeterminate’. The latter were not considered in the present study. Bone specimens recovered from sites located in regions outside of northern Europe were identified as belonging to wild or domesticated animals according to criteria established by individual analysts (Table S1). In addition, this study included mtDNA haplotype information to assist in classification (e.g. specimens yielding Near Eastern Y1 or Y2 lineages were classified as domesticated). In addition, the genetic analysis of a relatively large number of ancient specimens largely allows for zooarchaeological misclassification errors. The type of skeletal element used for aDNA analysis and the archaeological context from which each specimen is derived are presented in Table S1.

aDNA quality control procedures and mtDNA sequence analysis

All analyses were performed in the aDNA laboratory established at Kiel University, following the strict guidelines and procedures for work with minute amounts of degraded DNA2,16,17. Using previously described methods, mtDNA haplotypes were determined for an 80 bp diagnostic fragment from the control region of the Sus mitochondrial genome2. All sequencing results are shown in Table S1.

Stochastic models for the origin of mtDNA haplotypes in northern domestic pigs

In the following, pig subspecies (ps) are denoted by DP (domestic pig) and WB (wild boar), and regional affiliation (r) are denoted by NG (northern group) and SG (southern group). Then, P(h|ti, r, ps) is the conditional probability that an animal from subspecies ps living in region r during time period ti carried haplotype h. For the sake of brevity, we define pNG(ti) as the probability that a northern ti domestic pig descended from a northern tiβˆ’1 pig (either DP or WB), and pDP(ti) is the probability that a northern ti domestic pig descended from a domestic tiβˆ’1 pig, given that this predecessor came from the northern group.

In modelling the origin of mtDNA haplotypes of domestic pigs from the northern group, we will assume that no domestic pigs were present in that group during time period t1. A northern t2 domestic pig therefore descended either from a northern t1 wild boar or from a southern t1 domestic pig. This results in the following model equation:

Northern t3 domestic pigs are assumed to have descended from (i) northern t2 domestic pigs, (ii) northern t2 wild boars or (iii) from southern t2 domestic pigs. The resulting model equation thus reads

The mathematical models underlying the statistical analyses are relatively simple in several regards. Only one or two parameters are considered and the time and geography are treated as discreet (i.e. three time periods, two geographic regions). These simplifications take into account the limitation of (i) the sample size, which prohibits employment of multi-parameter models, and (ii) the level of accuracy achievable by way of archaeological dating. Therefore, the models appear adequate for the aims and scope of the present study acknowledging that more refined temporal or geographic structures could not be resolved.

Parameter estimation

The two model parameters pNG(ti), the probability that a northern ti domestic pig descended from a northern tiβˆ’1 pig, and pDP(ti), the probability that a northern tiβˆ’1 predecessor was domestic rather than wild, were estimated by maximum likelihood using absolute frequencies nh(ti), for short nh, of haplotypes A, C, Y1 and Y2 in the sampled domestic pigs from the northern group in time period ti. Assuming that the sampled pigs were stochastically independent, these data follow a multinomial distribution:

Parameters pNG(ti) and pDP(ti) enter into the likelihood model via Equations (1) or (2), depending upon the respective time period. Conditional probabilities P(h|t1;SG;DP), P(h|t2;NG;DP), P(h|t2;NG;WB) and P(h|t2;SG;DP) can be estimated from the available mtDNA data. For P(h|t1;NG;WB), however, no such data were available so that we had to assume P(h|t1;NG;WB)β€‰β‰ˆβ€‰P(h|t2;NG;WB).

Calculations were performed with statistics software R18. Numerical maximization of the likelihood was carried out with function optim allowing for box constraints on probabilities19. To deal with multiple maxima, we chose starting values from a grid of 12 equally spaced values between 0 and 1 and inspected the likelihood surface by means of a 2D or 3D plot. For the estimation of pNG for time period t2, confidence intervals were calculated by bootstrapping (n = 100,000 bootstrap samples). For the estimation of pNG and pDP for time period t3, bootstrapping and therefore the calculation of a confidence interval was not feasible because the three-dimensional likelihood surface was too complex.

Additional Information

How to cite this article: Caliebe, A. et al. Insights into early pig domestication provided by ancient DNA analysis. Sci. Rep. 7, 44550; doi: 10.1038/srep44550 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.