Multiple maternal origins of native modern and ancient horse populations in China

To obtain more knowledge of the origin and genetic diversity of domestic horses in China, this study provides a comprehensive analysis of mitochondrial DNA (mtDNA) D-loop sequence diversity from nine horse breeds in China in conjunction with ancient DNA data and evidence from archaeological and historical records. A 247-bp mitochondrial D-loop sequence from 182 modern samples revealed a total of 70 haplotypes with a high level of genetic diversity. Seven major mtDNA haplogroups (A–G) and 16 clusters were identified for the 182 Chinese modern horses. In the present study, nine 247-bp mitochondrial D-loop sequences of ancient remains of Bronze Age horse from the Chifeng region of Inner Mongolia in China ( c. 4000–2000a bp ) were used to explore the origin and diversity of Chinese modern horses and the phylogenetic relationship between ancient and modern horses. The nine ancient horses carried seven haplotypes with rich genetic diversity, which were clustered together with modern individuals among haplogroups A, E and F. Modern domestic horse and ancient horse data support the multiple origins of domestic horses in China. This study supports the argument that multiple successful events of horse domestication, including separate introductions of wild mares into the domestic herds, may have occurred in antiquity, and that China cannot be excluded from these events. Indeed, the association of Far Eastern mtDNA types to haplogroup F was highly significant using Fisher's exact test of independence ( P  = 0.00002), lending support for Chinese domestication of this haplogroup. High diversity and all seven mtDNA haplogroups (A–G) with 16 clusters also suggest that further work is necessary to shed more light on horse domestication in China.     

Ancient DNA reveals traces of Iberian Neolithic and Bronze Age lineages in modern Iberian horses

Multiple geographical regions have been proposed for the domestication of Equus caballus . It has been suggested, based on zooarchaeological and genetic analyses that wild horses from the Iberian Peninsula were involved in the process, and the overrepresentation of mitochondrial D1 cluster in modern Iberian horses supports this suggestion. To test this hypothesis, we analysed mitochondrial DNA from 22 ancient Iberian horse remains belonging to the Neolithic, the Bronze Age and the Middle Ages, against previously published sequences. Only the medieval Iberian sequence appeared in the D1 group. Neolithic and Bronze Age sequences grouped in other clusters, one of which (Lusitano group C) is exclusively represented by modern horses of Iberian origin. Moreover, Bronze Age Iberian sequences displayed the lowest nucleotide diversity values when compared with modern horses, ancient wild horses and other ancient domesticates using nonparametric bootstrapping analyses. We conclude that the excessive clustering of Bronze Age horses in the Lusitano group C, the observed nucleotide diversity and the local continuity from wild Neolithic Iberian to modern Iberian horses, could be explained by the use of local wild mares during an early Iberian domestication or restocking event, whereas the D1 group probably was introduced into Iberia in later historical times.     

Mitochondrial D-loop sequence variation among Italian horse breeds

The genetic variability of the mitochondrial D-loop DNA sequence in seven horse breeds bred in Italy (Giara, Haflinger, Italian trotter, Lipizzan, Maremmano, Thoroughbred and Sarcidano) was analysed. Five unrelated horses were chosen in each breed and twenty-two haplotypes were identified. The sequences obtained were aligned and compared with a reference sequence and with 27 mtDNA D-loop sequences selected in the GenBank database, representing Spanish, Portuguese, North African, wild horses and an Equus asinus sequence as the outgroup. Kimura two-parameter distances were calculated and a cluster analysis using the Neighbour-joining method was performed to obtain phylogenetic trees among breeds bred in Italy and among Italian and foreign breeds. The cluster analysis indicates that all the breeds but Giara are divided in the two trees, and no clear relationships were revealed between Italian populations and the other breeds. These results could be interpreted as showing the mixed origin of breeds bred in Italy and probably indicate the presence of many ancient maternal lineages with high diversity in mtDNA sequences.     

Genetic diversity of Cheju horses (Equus caballus) determined by using mitochondrial DNA D-loop polymorphism

We used sequence polymorphism of the mitochondrial DNA D-loop (968 bp excluding the tandem repeat region) to determine genetic diversity of horses inhabiting Cheju (a southern island of Korea). Seventeen haplotypes with frequencies from 1.5 to 21.5% were found among 65 Cheju horse samples. Genetic diversity (h) of the 17 haplotypes was calculated to be 0.91, indicating that the extant Cheju horse population consists of diverse genetic groups in their maternal lineage. Phylogenetic analysis showed that 17 types of Cheju (D-loop sequences determined), 5 Mongolian, 6 Arabian, 3 Belgian, 2 Tsushima, 2 Yunnan, 1 Przewalskii, and 3 Thoroughbred horses (published sequences for the latter seven breeds) showed that Cheju horses were distributed into many different clusters in the tree. Four Mongolian horses clustered with separate Cheju horse groups, showing that some Cheju horses are clearly of Mongolian origin. The analysis of partial sequences (284 bp) of the D-loop of 109 horses showed that Thoroughbred, Mongolian, Lipizzan, and Arabian breeds are as diverse as Cheju horses. Our data together with others' suggest that most horse breeds tested with reasonably sufficient numbers of samples are diverse in their maternal lineages and also are not uniquely different from each other.     

Genetic characterization of horse bone excavated from the Kwakji archaeological site,Jeju, Korea

We determined the nucleotide sequences of the hypervariable D-loop region of mitochondrial DNA (mtDNA) from horse bone (humerus, A.D. 700 to A.D. 800) that was excavated from the Kwakji archaeological site, Jeju, Korea. We compared them with ones from extant horses. We designed three pairs of oligonucleotide primers from the tRNA-Thr and tRNA-Phe gene regions of mtDNA that are highly conserved among many other animal species. We cloned 232, 336, and 644 bp from the horse bone in order to determine the mtDNA D-loop sequence. The sequence was 1,124 bp long; the middle contained 19 tandem repeats of an 8-bp sequence (TGTGCACC) that is specific to equines. The mtDNA D-loop region contained each base (total number, percentage of total) as follows: A (317, 28.20%), C (336, 29.89%), G (169, 15.04%), and T (302, 26.87%). This sequence, like those of other horse populations, was AT rich. Sequence divergence was the lowest (1.71%) between the ancient horse bone and that of the Thoroughbred horse 1. The neighbor-joining and strict consensus tree of three of the most parsimonious trees also suggested that the ancient bone was considerably unrelated to native Jeju horses. The molecular phylogenetic characteristics of the horse bone that was excavated from the Kwakji archaeological site (Jeju, Korea) showed that some horse breeds may have existed on Jeju Island, Korea before Mongolian horses were introduced. The horse bone that was excavated from the Kwakji archaeological site may aid future research on the origin and ancestry of native Jeju horses.     

Phylogenetic relationships of Cheju horses to other horse breeds as determined by mtDNA D-loop sequence polymorphism

Historical records suggest that horses inhabiting the island of Cheju in Korea are descendants of Mongolian horses introduced in 1276. Other studies, however, suggest that horses may have been present on the island prior to the Mongolian introduction. To determine the origin of the Cheju horses we used a phylogenetic analysis of sequences of the mitochondrial DNA (mtDNA) D-loop region, including tRNA Pro and parts of tRNA thr and tRNA Phe sequences (1102-bp excluding the tandem repeat region). Maximum parsimony and neighbor-joining trees were constructed using sequences determined for seven Cheju, four Mongolian, one Przewalskii and two Chinese Yunnan horses, and published sequences for one Swedish and three Thoroughbred horses. Donkey mtDNA was used as an outgroup. We found that the mtDNA D-loop sequence varies considerably within Mongolian, Cheju and Thoroughbred horse breeds, and that Cheju horses clustered with Mongolian horses as well as with horses from other distantly related breeds. On the basis of these findings we propose that horses on Cheju Island are of mixed origin in their maternal lineage, and that horses may have existed and been traded on the island before the Mongolian introduction.     

Polymorphic sequence in the D-loop region of equine mitochondrial DNA

The D-loop regions in equine mitochondrial DNA were cloned from three thoroughbred horses by polymerase chain reaction (PCR). The total number of bases in the D-loop region were 1114bp, 1115bp and 1146bp. The equine D-loop region is A/T rich like many other mammalian D-loops. The large central conserved sequence block and small conserved sequence blocks 1, 2 and 3, that are common to other mammals, were observed. Between conserved sequence blocks 1 and 2 there were tandem repeats of an 8bp equine-specific sequence TGTGCACC, and the number of tandem repeats differed among individual horses. The base composition in the unit of these repeats is G/C rich as are the short repeats in the D-loops of rabbit and pig. Comparing DNA sequences between horse and other mammals, the difference in the D-loop region length is mostly due to the difference in the number of DNA sequences at both extremities. The similarities of the DNA sequences are in the middle part of the D-loop. In comparison of the sequences among three thoroughbred horses, it was determined that the region between tRNA^Pro and the large central conserved sequence block was the richest in variation. PCR primers in the D-loop region were designed and the expected maternal inheritance was confirmed by PCR-RFLP (restriction fragment length polymorphism).     

Molecular phylogenetics and diagnosis of soil and clinical isolates of Halicephalobus gingivalis (Nematoda: Cephalobina: Panagrolaimoidea), an opportunistic pathogen of horses

Phylogenetic relationships among six isolates of Halicephalobus gingivalis (Stefanski, 1954), a species with pathogenic potential in horses and humans, were evaluated using DNA sequences from the nuclear large-subunit ribosomal RNA (LSU rDNA) gene. Sequences from nematodes obtained from in vitro cultures (soil or clinical sources), or isolated from infected horse tissues, were compared. Gene sequences from a fatal equine clinical case from southern California and a free-living isolate recovered from southern California soil showed no fixed differences. Sequences from isolates representing two fatal equine cases from North America, one from Ontario, Canada and another from Tennessee also showed no fixed differences. In contrast, two equine cases from Tennessee had 18 fixed differences for this LSU region, the greatest observed among isolates from horses. Phylogenetic analysis of six Halicephalobus sequences and four outgroup taxa by maximum parsimony yielded one tree with five well-supported clades. This phylogeny did not group isolates of Halicephalobus strictly by region of geographic isolation or source of sample, and depicted one clinical and one soil isolate as sister taxa. These results confirm that free-living environmental isolates are potential sources of infection for horses. The phylogeny also reveals that diverse isolates can cause infections in horses within a relatively limited geographic region, and conversely that genetically similar sister taxa can be recovered from geographically distant localities. PCR primers that selectively amplify Halicephalobus DNA were designed and tested based on comparison of closely related nematodes as inferred from phylogenetic analysis.     

A Statistical Approach to Identify Ancient Template DNA

One of the key problems in the study of ancient DNA is that of authenticating sequences obtained from PCR amplifications of highly degraded samples. Contamination of ancient samples and postmortem damage to endogenous DNA templates are the major obstacles facing researchers in this task. In particular, the authentication of sequences obtained from ancient human remains is thought by many to be rather challenging. We propose a novel approach, based on the c statistic, that can be employed to help identify the sequence motif of an endogenous template, based on a sample of sequences that reflect the nucleotide composition of individual template molecules obtained from ancient tissues (such as cloned products from a PCR amplification). The c statistic exploits as information the most common form of postmortem damage observed among clone sequences in ancient DNA studies, namely, lesion-induced substitutions caused by cytosine deamination events. Analyses of simulated sets of templates with miscoding lesions and real sets of clone sequences from the literature indicate that the c-based approach is highly effective in identifying endogenous sequence motifs, even when they are not present among the sampled clones. The proposed approach is likely to be of general use to researchers working with DNA from ancient tissues, particularly from human remains, where authentication of results has been most challenging. [Reviewing Editor: Dr. Magnus Nordborg]     

Mitochondrial DNA from prehistoric canids highlights relationships between dogs and South-East European wolves

The question of the origins of the dog has been much debated. The dog is descended from the wolf that at the end of the last glaciation (the archaeologically hypothesized period of dog domestication) was one of the most widespread among Holarctic mammals. Scenarios provided by genetic studies range from multiple dog-founding events to a single origin in East Asia. The earliest fossil dogs, dated approximately 17-12,000 radiocarbon ((14)C) years ago (YA), were found in Europe and in the Middle East. Ancient DNA (a-DNA) evidence could contribute to the identification of dog-founder wolf populations. To gain insight into the relationships between ancient European wolves and dogs we analyzed a 262-bp mitochondrial DNA control region fragment retrieved from five prehistoric Italian canids ranging in age from approximately 15,000 to approximately 3,000 (14)C YA. These canids were compared to a worldwide sample of 547 purebred dogs and 341 wolves. The ancient sequences were highly diverse and joined the three major clades of extant dog sequences. Phylogenetic investigations highlighted relationships between the ancient sequences and geographically widespread extant dog matrilines and between the ancient sequences and extant wolf matrilines of mainly East European origin. The results provide a-DNA support for the involvement of European wolves in the origins of the three major dog clades. Genetic data also suggest multiple independent domestication events. East European wolves may still reflect the genetic variation of ancient dog-founder populations.     

Mitochondrial sequence variation in ancient horses from the Carpathian Basin and possible modern relatives

Movements of human populations leave their traces in the genetic makeup of the areas affected; the same applies to the horses that move with their owners This study is concerned with the mitochondrial control region genotypes of 31 archaeological horse remains, excavated from pre-conquest Avar and post-conquest Hungarian burial sites in the Carpathian Basin dating from the sixth to the tenth century. To investigate relationships to other ancient and recent breeds, modern Hucul and Akhal Teke samples were also collected, and mtDNA control region (CR) sequences from 76 breeds representing 921 individual specimens were combined with our sequence data. Phylogenetic relationships among horse mtDNA CR haplotypes were estimated using both genetic distance and the non-dichotomous network method. Both methods indicated a separation between horses of the Avars and the Hungarians. Our results show that the ethnic changes induced by the Hungarian Conquest were accompanied by a corresponding change in the stables of the Carpathian Basin.     

Determination of genetic relatedness from low‐coverage human genome sequences using pedigree simulations

We develop and evaluate methods for inferring relatedness among individuals from low‐coverage DNA sequences of their genomes, with particular emphasis on sequences obtained from fossil remains. We suggest the major factors complicating the determination of relatedness among ancient individuals are sequencing depth, the number of overlapping sites, the sequencing error rate and the presence of contamination from present‐day genetic sources. We develop a theoretical model that facilitates the exploration of these factors and their relative effects, via measurement of pairwise genetic distances, without calling genotypes, and determine the power to infer relatedness under various scenarios of varying sequencing depth, present‐day contamination and sequencing error. The model is validated by a simulation study as well as the analysis of aligned sequences from present‐day human genomes. We then apply the method to the recently published genome sequences of ancient Europeans, developing a statistical treatment to determine confidence in assigned relatedness that is, in some cases, more precise than previously reported. As the majority of ancient specimens are from animals, this method would be applicable to investigate kinship in nonhuman remains. The developed software grups (Genetic Relatedness Using Pedigree Simulations) is implemented in Python and freely available.     

一种实用可靠的古代 DNA 获取方法

古代DNA序列信息能够为物种演化研究提供最直接的分子证据,但获取古代DNA的技术仍存在诸多瓶颈,尤其是扩增中存在受损伤DNA模板的干扰、获取成本高和实验周期长等问题．改进了异丙醇沉淀提取法,并采用了尿嘧啶糖苷酶(UNG)去除受损伤DNA模板后进行扩增的方法,最终可以高效地获取真实的古代DNA序列．实验利用距今4 300～3 900年前的猪牙样本,将改进的古 DNA 获取方法与常规方法进行比较研究,结果表明,改进的异丙醇沉淀法提取结合UNG处理后进行PCR扩增的方法,可以在保证古代DNA获取成功率并提高获得的DNA序列可靠性的前提下,将经费投入和实验周期都各减少至常规方法的50%以下．这可以为开展大规模古代样本检测提供一种切实可行的 DNA 获取方法．     

Phylogeography of lions (Panthera leo ssp.) reveals three distinct taxa and a late Pleistocene reduction in genetic diversity

Lions were the most widespread carnivores in the late Pleistocene, ranging from southern Africa to the southern USA, but little is known about the evolutionary relationships among these Pleistocene populations or the dynamics that led to their extinction. Using ancient DNA techniques, we obtained mitochondrial sequences from 52 individuals sampled across the present and former range of lions. Phylogenetic analysis revealed three distinct clusters: (i) modern lions, Panthera leo ; (ii) extinct Pleistocene cave lions, which formed a homogeneous population extending from Europe across Beringia (Siberia, Alaska and western Canada); and (iii) extinct American lions, which formed a separate population south of the Pleistocene ice sheets. The American lion appears to have become genetically isolated around 340 000 years ago, despite the apparent lack of significant barriers to gene flow with Beringian populations through much of the late Pleistocene. We found potential evidence of a severe population bottleneck in the cave lion during the previous interstadial, sometime after 48 000 years, adding to evidence from bison, mammoths, horses and brown bears that megafaunal populations underwent major genetic alterations throughout the last interstadial, potentially presaging the processes involved in the subsequent end-Pleistocene mass extinctions.     

Chinese Mongolian horses may retain early domestic male genetic lineages yet to be discovered

The Mongolian horse represents one of the most ancient extant horse populations. In this study we determined the male‐specific region of the Y chromosome (MSY) haplotype distribution in 60 Chinese Mongolian horses representing five distinct populations. Cosmopolitan male lineages were predominant in horses from one improved (Sanhe), one Chinese Mongolian subtype (Baicha Iron Hoof) and one indigenous (Abaga Black) population. In contrast, autochthonous Y chromosome diversity was evident among the two landrace populations (Wushen and Wuzhumuqin), as the majority of their MSY haplotypes were situated at root nodes in a network. Our results also suggest gene flow between Chinese Mongolian and Arabian horses, as an appreciable number of Wuzhumuqin horses carried haplotypes that are typically observed in Arabian horses. Although most horses carried modern haplotypes as a direct result of recent breed improvement, authentic Chinese Mongolian horses retain an ancient signature of paternal lineages that has not previously been described in extant horse populations. Therefore, further characterization of MSY variation in these populations will be important for the discovery of lost diversity in modern domestic horses and also for understanding the evolutionary history of equine paternal lineages.     

Tracking down human contamination in ancient human teeth

DNA contamination arising from the manipulation of ancient calcified tissue samples is a poorly understood, yet fundamental, problem that affects the reliability of ancient DNA (aDNA) studies. We have typed the mitochondrial DNA hypervariable region I of the only 6 people involved in the excavation, washing, and subsequent anthropological and genetic study of 23 Neolithic remains excavated from Granollers (Barcelona, Spain) and searched for their presence among the 572 clones generated during the aDNA analyses of teeth from these samples. Of the cloned sequences, 17.13% could be unambiguously identified as contaminants, with those derived from the people involved in the retrieval and washing of the remains present in higher frequencies than those of the anthropologist and genetic researchers. This finding confirms, for the first time, previous hypotheses that teeth samples are most susceptible to contamination at their initial excavation. More worrying, the cloned contaminant sequences exhibit substitutions that can be attributed to DNA damage after the contamination event, and we demonstrate that the level of such damage increases with time: contaminants that are >10 years old have approximately 5 times more damage than those that are recent. Furthermore, we demonstrate that in this data set, the damage rate of the old contaminant sequences is indistinguishable from that of the endogenous DNA sequences. As such, the commonly used argument that miscoding lesions observed among cloned aDNA sequences can be used to support data authenticity is misleading in scenarios where the presence of old contaminant sequences is possible. We argue therefore that the typing of those involved in the manipulation of the ancient human specimens is critical in order to ensure that generated results are accurate.     

Phylogenetic relationships of Argentinean Creole horses and other South American and Spanish breeds inferred from mitochondrial DNA sequences

South American horses constitute a direct remnant of the Iberian horses brought to the New World by the Spanish conquerors. The source of the original horses was Spain, and it is generally assumed that the animals belonged to the Andalusian, Spanish Celtic, Barb or Arabian breeds. In order to establish the relationship between Argentinean and Spanish horses, a portion of the mitochondrial D-loop of 104 animals belonging to nine South American and Spanish breeds was analysed using SSCP and DNA sequencing. The variability found both within and between breeds was very high. There were 61 polymorphic positions, representing 16% of the total sequence obtained. The mean divergence between a pair of sequences was 2.8%. Argentinean Creole horses shared two haplotypes with the Peruvian Paso from Argentina, and the commonest haplotype of the Creole horses is identical to one of the Andalusian horses. Even when there was substantial subdivision between breeds with highly significant Wright's Fixation Index (FST), the parsimony and distance-based phylogenetic analyses failed to show monophyletic groups and there was no clear relationship in the trees between the South American and any of the other horses analysed. Although this result could be interpreted as mixed ancestry of the South American breeds with respect to the Spanish breeds, it is probably indicating the retention of very ancient maternal lineages in the breeds analysed.     

Brief communication: identification of the authentic ancient DNA sequence in a human bone contaminated with modern DNA

We present a method to distinguish authentic ancient DNA from contaminating DNA in a human bone. This is achieved by taking account of the spatial distribution of the various sequence families within the bone and the extent of degradation of the template DNAs, as revealed by the error content of the sequences. To demonstrate the veracity of the method, we handled two ancient human tibiae in order to contaminate them with modern DNA, and then subjected segments of the bones to various decontaminating treatments, including removal of the outer 1-2 mm, before extracting DNA, cloning, and obtaining a total of 107 mitochondrial DNA sequences. Sequences resulting from the deliberate contamination were located exclusively in the outer 1-2 mm of the bones, and only one of these 27 sequences contained an error that could be ascribed to DNA degradation. A second, much smaller set of relatively error-free sequences, which we ascribe to contamination during excavation or curation, was also located exclusively in the outer 1-2 mm. In contrast, a family of 72 sequences, displaying extensive degradation products but identifiable as haplogroup U5a1a, was distributed throughout one of the bones and represents the authentic ancient DNA content of this specimen.     

Evolution, systematics, and phylogeography of pleistocene horses in the new world: a molecular perspective

The rich fossil record of horses has made them a classic example of evolutionary processes. However, while the overall picture of equid evolution is well known, the details are surprisingly poorly understood, especially for the later Pliocene and Pleistocene, c. 3 million to 0.01 million years (Ma) ago, and nowhere more so than in the Americas. There is no consensus on the number of equid species or even the number of lineages that existed in these continents. Likewise, the origin of the endemic South American genus Hippidion is unresolved, as is the phylogenetic position of the "stilt-legged" horses of North America. Using ancient DNA sequences, we show that, in contrast to current models based on morphology and a recent genetic study, Hippidion was phylogenetically close to the caballine (true) horses, with origins considerably more recent than the currently accepted date of c. 10 Ma. Furthermore, we show that stilt-legged horses, commonly regarded as Old World migrants related to the hemionid asses of Asia, were in fact an endemic North American lineage. Finally, our data suggest that there were fewer horse species in late Pleistocene North America than have been named on morphological grounds. Both caballine and stilt-legged lineages may each have comprised a single, wide-ranging species.     

PCR-RFLP analysis of the cytochrome b gene in horse mitochondrial DNA

The mitochondrial DNA sequence of cytochrome b gene in a Thoroughbred horse was determined. By comparing DNA sequences between the Thoroughbred and published sequence data (two horses and one Grevyi zebra), polymerase chain reaction (PCR) primers were designed for amplification of a 590 bp DNA fragment in the cytochrome b gene, and PCR-restriction fragment length polymorphism (RFLP) analysis was studied in 140 horses of six breeds using three restriction enzymes ( AciI, BamHI, RsaI ). Two morphs were found using each of the three enzymes. By combining three enzymes morphs, the 140 horses examined were classified into four types. Type 2 was most frequent in all breeds.