Bioarchaeology aids the cultural understanding of six characters in search of their agency (Tarquinia, ninth-seventh century BC, central Italy).

Etruria contained one of the great early urban civilisations in the Italian peninsula during the first millennium BC, much studied from a cultural, humanities-based, perspective, but relatively little with scientific data, and rarely in combination. We have addressed the unusual location of twenty inhumations found in the sacred heart of the Etruscan city of Tarquinia, focusing on six of these as illustrative, contrasting with the typical contemporary cremations found in cemeteries on the edge of the city. The cultural evidence suggests that the six skeletons were also distinctive in their ritualization and memorialisation. Focusing on the six, as a representative sample, the scientific evidence of osteoarchaeology, isotopic compositions, and ancient DNA has established that these appear to show mobility, diversity and violence through an integrated bioarchaeological approach. The combination of multiple lines of evidence makes major strides towards a deeper understanding of the role of these extraordinary individuals in the life of the early city of Etruria.

Imputation of Ancient Whole Genome Sus scrofa DNA Introduces Biases Toward Main Population Components in the Reference Panel.

Sequencing ancient DNA to high coverage is often limited by sample quality and cost. Imputing missing genotypes can potentially increase information content and quality of ancient data, but requires different computational approaches than modern DNA imputation. Ancient imputation beyond humans has not been investigated. In this study we report results of a systematic evaluation of imputation of three whole genome ancient Sus scrofa samples from the Early and Late Neolithic (∼7,100-4,500 BP), to test the utility of imputation. We show how issues like genetic architecture and, reference panel divergence, composition and size affect imputation accuracy. We evaluate a variety of imputation methods, including Beagle5, GLIMPSE, and Impute5 with varying filters, pipelines, and variant calling methods. We achieved genotype concordance in most cases reaching above 90%; with the highest being 98% with ∼2,000,000 variants recovered using GLIMPSE. Despite this high concordance the sources of diversity present in the genotypes called in the original high coverage genomes were not equally imputed leading to biases in downstream analyses; a trend toward genotypes most common in the reference panel is observed. This demonstrates that the current reference panel does not possess the full diversity needed for accurate imputation of ancient Sus, due to missing variations from Near Eastern and Mesolithic wild boar. Imputation of ancient Sus scrofa holds potential but should be approached with caution due to these biases, and suggests that there is no universal approach for imputation of non-human ancient species.

Variance components for bovine tuberculosis infection and multi-breed genome-wide association analysis using imputed whole genome sequence data.

Bovine tuberculosis (bTB) is an infectious disease of cattle generally caused by Mycobacterium bovis, a bacterium that can elicit disease humans. Since the 1950s, the objective of the national bTB eradication program in Republic of Ireland was the biological extinction of bTB; that purpose has yet to be achieved. Objectives of the present study were to develop the statistical methodology and variance components to undertake routine genetic evaluations for resistance to bTB; also of interest was the detection of regions of the bovine genome putatively associated with bTB infection in dairy and beef breeds. The novelty of the present study, in terms of research on bTB infection, was the use of beef breeds in the genome-wide association and the utilization of imputed whole genome sequence data. Phenotypic bTB data on 781,270 animals together with imputed whole genome sequence data on 7,346 of these animals' sires were available. Linear mixed models were used to quantify variance components for bTB and EBVs were validated. Within-breed and multi-breed genome-wide associations were undertaken using a single-SNP regression approach. The estimated genetic standard deviation (0.09), heritability (0.12), and repeatability (0.30) substantiate that genetic selection help to eradicate bTB. The multi-breed genome-wide association analysis identified 38 SNPs and 64 QTL regions associated with bTB infection; two QTL regions (both on BTA23) identified in the multi-breed analysis overlapped with the within-breed analyses of Charolais, Limousin, and Holstein-Friesian. Results from the association analysis, coupled with previous studies, suggest bTB is controlled by an infinitely large number of loci, each having a small effect. The methodology and results from the present study will be used to develop national genetic evaluations for bTB in the Republic of Ireland. In addition, results can also be used to help uncover the biological architecture underlying resistance to bTB infection in cattle.

Ancient Ethiopian genome reveals extensive Eurasian admixture throughout the African continent.

Characterizing genetic diversity in Africa is a crucial step for most analyses reconstructing the evolutionary history of anatomically modern humans. However, historic migrations from Eurasia into Africa have affected many contemporary populations, confounding inferences. Here, we present a 12.5× coverage ancient genome of an Ethiopian male ("Mota") who lived approximately 4500 years ago. We use this genome to demonstrate that the Eurasian backflow into Africa came from a population closely related to Early Neolithic farmers, who had colonized Europe 4000 years earlier. The extent of this backflow was much greater than previously reported, reaching all the way to Central, West, and Southern Africa, affecting even populations such as Yoruba and Mbuti, previously thought to be relatively unadmixed, who harbor 6 to 7% Eurasian ancestry.

Paging through history: parchment as a reservoir of ancient DNA for next generation sequencing.

Parchment represents an invaluable cultural reservoir. Retrieving an additional layer of information from these abundant, dated livestock-skins via the use of ancient DNA (aDNA) sequencing has been mooted by a number of researchers. However, prior PCR-based work has indicated that this may be challenged by cross-individual and cross-species contamination, perhaps from the bulk parchment preparation process. Here we apply next generation sequencing to two parchments of seventeenth and eighteenth century northern English provenance. Following alignment to the published sheep, goat, cow and human genomes, it is clear that the only genome displaying substantial unique homology is sheep and this species identification is confirmed by collagen peptide mass spectrometry. Only 4% of sequence reads align preferentially to a different species indicating low contamination across species. Moreover, mitochondrial DNA sequences suggest an upper bound of contamination at 5%. Over 45% of reads aligned to the sheep genome, and even this limited sequencing exercise yield 9 and 7% of each sampled sheep genome post filtering, allowing the mapping of genetic affinity to modern British sheep breeds. We conclude that parchment represents an excellent substrate for genomic analyses of historical livestock.

Genome-wide association study for calving traits in Holstein-Friesian dairy cattle.

Dystocia and perinatal mortality are quantitative traits that significantly impact animal productivity and welfare. Their economic importance is reflected by their inclusion in the national breeding goals of many cattle populations. The genetic architecture that influences these traits, however, has still yet to be thoroughly defined. Regions of the bovine genome associated with calving difficulty (direct and maternal) and perinatal mortality were detected in this study using a Bayesian approach with 43 204 single nucleotide polymorphisms (SNPs) on up to 1970 Holstein-Friesian bulls. Several SNPs on chromosomes 5, 6, 11, 12, 17,18 and 28 were detected to be strongly associated with these calving performance traits. Novel genomic regions with previously reported associations with growth, stature, birth weight and bone morphology were identified in the present study as being associated with the three calving performance traits. Morphological abnormalities are a known contributor to perinatal mortality and the most significantly associated SNP for perinatal mortality in the present study was located in a region in linkage disequilibrium with the gene SLC26A7. This gene, SLC26A7, has similarities and colocalises with SLC4A2, which has previously been associated with osteoporosis and mortality in cattle populations. The HHIP gene that is known to be associated with stature in humans was strongly associated with direct calving difficulty in the present study; large calves are known to, on average, have a greater likelihood of dystocia. A stemloop microRNA, bta-mir-1256, on chromosome 12, involved in post-transcriptional regulation of gene expression was associated with maternal calving difficulty. Previously reported quantitative trait loci associated with calving performance traits in other populations were again identified in this study; with one genomic region on chromosome 18 supporting very strong evidence of an underlying causative mutation and accounting for 2.1% of the genetic variation in direct calving difficulty. Overlapping genomic regions associated with one or more of the calving traits were also detected substantiating the known genetic covariances existing between these traits. Moreover, some genomic regions were only associated with one of the calving traits implying the selective genomic breeding programs exploiting these regions could help resolve genetic antagonisms.

High concordance of bovine single nucleotide polymorphism genotypes generated using two independent genotyping strategies.

Single nucleotide polymorphisms (SNPs) represent the most common form of DNA sequence variation in mammalian livestock genomes. While the past decade has witnessed major advances in SNP genotyping technologies, genotyping errors caused, in part, by the biochemistry underlying the genotyping platform used, can occur. These errors can distort project results and conclusions and can result in incorrect decisions in animal management and breeding programs; hence, SNP genotype calls must be accurate and reliable. In this study, 263 Bos spp. samples were genotyped commercially for a total of 16 SNPs. Of the total possible 4,208 SNP genotypes, 4,179 SNP genotypes were generated, yielding a genotype call rate of 99.31% (standard deviation ± 0.93%). Between 110 and 263 samples were subsequently re-genotyped by us for all 16 markers using a custom-designed SNP genotyping platform, and of the possible 3,819 genotypes a total of 3,768 genotypes were generated (98.70% genotype call rate, SD ± 1.89%). A total of 3,744 duplicate genotypes were generated for both genotyping platforms, and comparison of the genotype calls for both methods revealed 3,741 concordant SNP genotype call rates (99.92% SNP genotype concordance rate). These data indicate that both genotyping methods used can provide livestock geneticists with reliable, reproducible SNP genotypic data for in-depth statistical analysis.

Variation in chicken populations may affect the enzymatic activity of lysozyme.

The chicken lysozyme gene encodes a hydrolase that has a key role in defence, especially in ovo. This gene was resequenced in global chicken populations [red, grey, Ceylon and green jungle fowl (JF)] and related bird species. Networks, summary statistics and tests of neutrality indicate that although there is extensive variation at the gene, little is present at coding sites, with the exception of one non-synonymous site. This segregating site and a further fixed non-synonymous change between red JF and domestic chicken populations are spatially close to the catalytic sites of the enzyme and so might affect its activity.

Maternal and paternal genealogy of Eurasian taurine cattle (Bos taurus).

Maternally inherited mitochondrial DNA (mtDNA) has been used extensively to determine origin and diversity of taurine cattle (Bos taurus) but global surveys of paternally inherited Y-chromosome diversity are lacking. Here, we provide mtDNA information on previously uncharacterised Eurasian breeds and present the most comprehensive Y-chromosomal microsatellite data on domestic cattle to date. The mitochondrial haplogroup T3 was the most frequent, whereas T4 was detected only in the Yakutian cattle from Siberia. The mtDNA data indicates that the Ukrainian and Central Asian regions are zones where hybrids between taurine and zebu (B. indicus) cattle have existed. This zebu influence appears to have subsequently spread into southern and southeastern European breeds. The most common Y-chromosomal microsatellite haplotype, termed here as H11, showed an elevated frequency in the Eurasian sample set compared with that detected in Near Eastern and Anatolian breeds. The taurine Y-chromosomal microsatellite haplotypes were found to be structured in a network according to the Y-haplogroups Y1 and Y2. These data do not support the recent hypothesis on the origin of Y1 from the local European hybridization of cattle with male aurochsen. Compared with mtDNA, the intensive culling of breeding males and male-mediated crossbreeding of locally raised native breeds has accelerated loss of Y-chromosomal variation in domestic cattle, and affected the contribution of genetic drift to diversity. In conclusion, to maintain diversity, breeds showing rare Y-haplotypes should be prioritised in the conservation of cattle genetic resources.

Cytochrome b sequences of ancient cattle and wild ox support phylogenetic complexity in the ancient and modern bovine populations.

Mitochondrial DNA has been the traditional marker for the study of animal domestication, as its high mutation rate allows for the accumulation of molecular diversity within the time frame of domestic history. Additionally, it is exclusively maternally inherited and haplotypes become part of the domestic gene pool via actual capture of a female animal rather than by interbreeding with wild populations. Initial studies of British aurochs identified a haplogroup, designated P, which was found to be highly divergent from all known domestic haplotypes over the most variable portion of the D-loop. Additional analysis of a large and geographically representative sample of aurochs from northern and central Europe found an additional, separate aurochs haplotype, E. Until recently, the European aurochs appeared to have no matrilinear descendants among the publicly available modern cattle control regions sequenced; if aurochs mtDNA was incorporated into the domestic population, aurochs either formed a very small proportion of modern diversity or had been subsequently lost. However, a haplogroup P sequence has recently been found in a modern sample, along with a new divergent haplogroup called Q. Here we confirm the outlying status of the novel Q and E haplogroups and the modern P haplogroup sequence as a descendent of European aurochs, by retrieval and analysis of cytochrome b sequence data from twenty ancient wild and domesticated cattle archaeological samples.

Bayesian inference of population expansions in domestic bovines.

The past population dynamics of four domestic and one wild species of bovine were estimated using Bayesian skyline plots, a coalescent Markov chain Monte Carlo method that does not require an assumed parametric model of demographic history. Four domestic species share a recent rapid population expansion not visible in the wild African buffalo (Syncerus caffer). The estimated timings of the expansions are consistent with the archaeological records of domestication.

Differentiation of European cattle by AFLP fingerprinting.

The Neolithic introduction of domestic cattle into Europe was followed by differential adaptation, selection, migration and genetic isolation, leading ultimately to the emergence of specialized breeds. We have studied the differentiation of European cattle by amplified fragment length polymorphism (AFLP) fingerprinting. Combining AFLP data sets from two laboratories yielded 81 biallelic polymorphic markers scored in 19-22 individual animals from 51 breeds. Model-based clustering differentiated Podolian cattle as well as French and Alpine breeds from other European cattle. AFLP genetic distances correlated well with microsatellite-based genetic distances calculated for the same breeds. However, the AFLP data emphasized the divergence of taurine and indicine cattle relative to the variation among European breeds and indicated an Eastern influence on Italian and Hungarian Podolian breeds. This probably reflects import from the East after the original introduction of domestic cattle into Europe. Our data suggest that Italian cattle breeds are relatively diverse at the DNA sequence level.

Early history of European domestic cattle as revealed by ancient DNA.

We present an extensive ancient DNA analysis of mainly Neolithic cattle bones sampled from archaeological sites along the route of Neolithic expansion, from Turkey to North-Central Europe and Britain. We place this first reasonable population sample of Neolithic cattle mitochondrial DNA sequence diversity in context to illustrate the continuity of haplotype variation patterns from the first European domestic cattle to the present. Interestingly, the dominant Central European pattern, a starburst phylogeny around the modal sequence, T3, has a Neolithic origin, and the reduced diversity within this cluster in the ancient samples accords with their shorter history of post-domestic accumulation of mutation.

Mitochondrial DNA sequence diversity in extant Irish horse populations and in ancient horses.

Equine mitochondrial DNA sequence variation was investigated in three indigenous Irish horse populations (Irish Draught Horse, Kerry Bog Pony and Connemara Pony) and, for context, in 69 other horse populations. There was no evidence of Irish Draught Horse or Connemara Pony sequence clustering, although the majority of Irish Draught Horse sequences (47%) were assigned to haplogroup D. Conversely, 31% of the Kerry Bog Pony sequences were assigned to the rare haplogroup E. In addition to the extant population analyses, ancient DNA sequences were generated from three out of four Irish archaeological specimens, all of which were assigned to haplogroup A.

Combination of multiple microsatellite data sets to investigate genetic diversity and admixture of domestic cattle.

Microsatellite markers are commonly used for population genetic analyses of livestock. However, up to now, combinations of microsatellite data sets or comparison of population genetic parameters from different studies and breeds has proven difficult. Often different genotyping methods have been employed, preventing standardization of microsatellite allele calling. In other cases different sets of markers have been genotyped, providing differing estimates of population genetic parameters. Here, we address these issues and illustrate a general two-step regression approach in cattle using three different sets of microsatellite data, to combine population genetics estimates of diversity and admixture. This regression-based method is independent of the loci genotyped but requires common breeds in the data sets. We show that combining microsatellite data sets can provide new insights on the origin and geographical distribution of genetic diversity and admixture in cattle, which will facilitate global management of this livestock species.

Levels of interpopulation differentiation among different functional classes of immunologically important genes.

It has been postulated that gene function may influence the degree to which allele frequencies differ among populations. In order to evaluate this effect, genotypic data from resequencing studies of genes classified as cytokines, cytokine receptors, cell adhesion molecules, Toll-like receptors and coagulation proteins were analysed for genetic differentiation (FST) between population samples of European and African descent. FST values did not differ statistically among functional groups when all polymorphic sites were included in the analyses. However, analysis based on nonsynonymous SNPs alone suggested weak heterogeneity among functional classes (P=0.0424). Particularly high levels of differentiation were shown by individual nonsynonymous SNPs at some genes, most notably ICAM1 and some Toll-like receptors. These genes interact directly with pathogens, and may therefore have been subject to geographically localised natural selection. Such loci warrant particular attention in studies of genetic disease risk and local adaptation to environmental conditions.

The phylogenetic position of the 'giant deer' Megaloceros giganteus.

The giant deer, or 'Irish elk', has featured extensively in debates on adaptation, sexual selection, and extinction. Its huge antlers--the largest of any deer species, living or extinct--formed a focus of much past work. Yet the phylogenetic position of the giant deer has remained an enigma. On the basis of its flattened antlers, the species was previously regarded as closely related to the living fallow deer. Recent morphological studies, however, have challenged that view and placed the giant deer closer to the living red deer or wapiti. Here we present a new phylogenetic analysis encompassing morphological and DNA sequence evidence, and find that both sets of data independently support a sister-group relationship of giant and fallow deer. Our results include the successful extraction and sequencing of DNA from this extinct species, and highlight the value of a joint molecular and morphological approach.

The mtDNA catalogue of all Portuguese autochthonous goat (Capra hircus) breeds: high diversity of female lineages at the western fringe of European distribution.

Goat is believed to be the first true livestock domesticated and, apart from its historical importance, keeps playing an essential economic role in very diverse human societies. We have analysed the female gene pool of all Portuguese autochthonous breeds (Bravia, Serrana, Charnequeira, Serpentina and Algarvia) through the mtDNA HVI sequencing of 288 unrelated animals sampled throughout the country. All breeds proved to be extremely diverse (average haplotype diversity of 0.977), in contrast with the Portuguese peripheral geographic situation in the distribution range of the species. Moreover, observed genetic distances between breeds do not correlate with microgeography inside Portugal. These observations are consistent with recurrent refreshment of the breeding stock through the introduction of exotic animals. Fitting the new data into the still loosely defined female genetic pool landscape of goats, all Portuguese animals, one sample excepted (belonging to Bravia), are classified into haplogroup A and haplotype sharing is geographically very sparse, including a Far East match. Our results confirm that goats stand out among most of domesticates as exceptionally diverse and showing an unparalleled degree of mobility of animals (at least females) used for reproductive purposes.

Evidence for the hypothesis of strain formation in Pomphorhynchus laevis (Acanthocephala): an investigation using mitochondrial DNA sequences.

A hypothesis has been erected stating that in the British Isles the acanthocephalan, Pomphorhynchus laevis can be separated into 3 strains, an English, Irish and marine strain. Ecological and morphological evidence exists in support of this hypothesis. An investigation at the molecular level was conducted in order to test the validity of the existing evidence. A mitochondrial gene, subunit I of cytochrome c oxidase was partially sequenced from 3 Irish populations of P. laevis, 1 Scottish population and 3 English populations. P. laevis sequences from brown trout from Ireland, England and Scotland were very similar, showing a mean sequence divergence of 0.7%. Sequences from two populations of P. laevis from English chub and bullhead were also similar to each other (0.35% divergence). These two groups of sequences, the brown trout group and the chub/bullhead group were 2.2% different. These data confirm the existence of at least 2 strains in Ireland and Britain, although there is evidence to suggest that these strains are defined by their host species rather than by their geographical distributions.