Genetic Diversity of mtDNA D-loop and Maternal Origin of Three Chinese Native Horse Breeds

Article information

Asian-Australas J Anim Sci. 2012;25(7):921-926
*Corresponding Author: Tao Zhang. Tel: +86-916-2641661, Fax: +86-916-2641564, E-mail: zl780823@163.com
School of Bioscience and Engineering, Shaanxi University of Technology, Hanzhong 723000, Shaanxi, China
Received 2011 December 13; Accepted 2012 February 16; Revised 2012 May 04.

Abstract

In order to protect the genetic resource of native horse breeds, the genetic diversity of mitochondrial DNA (mtDNA) D-loop of three native horse breeds in western China were investigated. Forty-three 600 bp mtDNA D-loop sequences were analyzed by PCR and sequencing techniques, 33 unique haplotypes with 70 polymorphic sites were detected in these horses, which account for 11.67% of 600 bp sequence analyzed, showing the abundant genetic diversity of the three native horse breeds in western China. The Neighbour-Joining (NJ) phylogenetic tree based on 247 bp of 43 D-loop sequences demonstrated the presence of seven major lineages (A to G), indicating that the three native horse breeds in western China originated from multiple maternal origins. Consistent with the front, the NJ phylogenetic tree based on 600 bp of mtDNA D-loop sequences of 43 Chinese western native horses and 81 sequences of six horse breeds from GenBank indicated that the three horse breeds had distributed into the seven major lineages (A to G). The structure of the phylogenic tree is often blurred because the variation in a short segment of the mitochondrial genome is often accompanied by high levels of recurrent mutations. Consequently, longer D-loop sequences are helpful in achieving a higher level of molecular resolution in horses.

INTRODUCTION

Horses have played a pivotal role in the development of human society. As an important domestic animal, the origin and history of modern domestic horses have been studied for decades in the biology and archaeology fields. At present, it is widely believed that horses were domesticated from several wild populations 4,000 to 6,000 years ago (Vila et al., 2001; Outram et al., 2009). Subsequently, they differentiated into many types after long-term domestication (Jansen et al., 2002). In history, horses differed from other domestic animals that were kept as food sources, and attained a prominent role as animals of transport and warfare, changing societies on a continent-wide scale. Literature reports the importation of horses into China from different parts of the world throughout history (Xie, 1987). As a strategic military asset, governments paid much attention to horse breeding throughout Chinese history. Some archaeological evidence and genetic data supports China as one of the earliest areas of horse domestication (Lei et al., 2009). China has rich horse genetic resources, the number of domestic horses in China was 7.639 million in 2004 (Statistic report in China, 2004). Chinese indigenous horse breeds are classified into five types based on geographical distribution and characteristics. They are Mongolian, Kazakhstan, Hequ, Southwest and Tibetan types (Chang, 2009). Four horse types are distributed in western China and still act as an important source of animal power because regions such as Guizhou, Shaanxi and Xinjiang are mountainous and pastoral. Ningqiang horse is an ancient breed. Individuals of this breed with a height at the withers less than or equal to 106 cm are called Ningqiang pony, while the other with a height greater than 106 cm are called Ningqiang horse (Hou, 1989). With the development of mechanization, the role of horse as a provider of transport and draft power is decreasing. The amount of horses declined quickly in recent years. So it is urgent to study and protect horse genetic resources.

The mtDNA is a small extra-chromosomal genome that is maternally inherited. The equine mtDNA has an approximate length of 16,660 bp and consists of a coding region which contains 37 genes and a non-coding region called the D-loop region of approximately 1,192 bp (Wolstenholme, 1992; Boore, 1999; Bowling et al., 2000). The mtDNA sequences have been widely used to study taxonomy, phylogeny, genetic structure and biological identity (Peng et al., 2007; Ji et al., 2008; Ki et al., 2010). D-loop mtDNA is an important tool to study equine diversity, phylogenetic relationships and origins, such as the Mongolian horse (Lei et al., 2009), Cheju horse (Yang et al., 2002), Debao pony (Jiang et al., 2011), Japanese horse (Kakoi et al., 2007), Tibetan horse (Xu et al., 2007), Przewalskii’s wild horse (Ishida et al., 1995) and donkey (Lei et al., 2005). The diversity and origin of the Ningqiang pony, Guizhou horse and Kazakh horse as indicated by microsatellite markers and mtDNA D-loop have been reported (Du et al., 2009; Lei et al., 2009). However, the mtDNA D-loop sequence (247 bp) studied in the literature was short and might not fully show the diversity and origin of the three horse breeds in western China. Consequently, a 600 bp mtDNA control region of 43 animals from Ningqiang pony, Guizhou and Kazakh breeds and 81 mtDNA D-loop sequences of 6 breeds from GenBank (Table 1) were used in this present study to determine the diversity and origin of the three native horse breeds in western China.

Cited equine D- loop sequences from GenBank

MATERIALS AND METHODS

Blood sample collection and DNA isolation

Blood samples were obtained from 43 horses of three western China native horse breeds including Ningqiang pony, Guizhou and Kazakh horse breeds and stored at −70°C. Ningqiang pony samples were collected from Ningqiang county, Hanzhong city, Shaanxi province; Guizhou horse samples came from Bijie region, Guizhou province; Kazakh horse samples were collected from Shihezi city, Xinjiang Autonomous Region. In addition, 81 mtDNA D-loop sequences from six domestic horse breeds and one donkey sequence as an outgroup were selected from GenBank (Table 1). The DNA used in PCR amplification was isolated from blood samples using the phenol-chloroform procedure.

PCR amplification and sequencing

The mtDNA D-loop sequences was amplified by the following primers: F: 5′-CGACAACAATTCACCCTCAT-3′, R: 5′-GAAGAAGGGTTGACAGATTTA-3′. Amplification was performed in a 25 μl reaction containing 50 ng of genomic DNA, 0.25 μM of each primer, 0.25 mM dNTPs, 2.5 mM MgCl2, 1.5 U Taq DNA polymerase and 10×buffer. Thermal cycling was carried out on an MJ Research 200 Thermal cycler. The procedure included predenaturation at 94°C (3 min), followed by 34 cycles of 94°C (30 s), 55°C (30 s), 72°C (90 s), and then final extension at 72°C for 10 min. The 1,100 bp PCR production was purified and two-way sequenced by ABI 3130 automated sequencer.

Statistical analysis

In this study, 43 mtDNA D-loop sequences of three horse breeds were sequenced, and aligned with the DNAMAN4.0 software (Tamura et al., 2007). Because of the number of repeats was variable within individuals (Kim et al., 1999), the repeat region was excluded and these sequences were truncated to 600 bp. Polymorphic sites, haplotype numbers, haplotype diversity (h) and nucleotide diversity (p) were calculated by DNASP 5.0 software (Rozas et al., 2003). The NJ tree of different horse breeds were constructed based on the Kimum 2-parameter model.

RESULTS

Genetic polymorphism of mtDNA D-loop region of three Chinese native horse breeds

Equine 600 bp of mtDNA D-loop sequence was examined in 43 Chinese domestic horses. The average ratio of the four nucleotides A, T, C, G was 28.3%, 27.5%, 27.0% and 17.2%, respectively.The contents of A+T exceeded G+C, showing that A+T was richer at the mtDNA D-loop region. This was in accordance with the order of nucleotide composition in the vertebrate mitochondrial genome, which is A>C>T>G with a skewed base composition and with more A+T than G+C base pairs (Asakawa et al., 1991; Ji et al., 2008). Thirty three mtDNA haplotypes were identified. 70 polymorphic sites (Figure 1) were detected, which represented 11.67% of the total mtDNA sequence analyzed (600 bp). Only one of the 70 polymorphisms represented insertion/deletion of single base pairs. The remaining were 67 transitions and two transversions, indicating a strong bias towards transitions that was common in mammalian mitochondrial evolution (Kim et al., 1999).

Figure 1

The polymorphic sites of 33 mtDNA D-loop haplotypes of three Chinese horse breeds.

Thirty-three mtDNA haplotypes were detected in 43 Chinese western native horse samples (Table 2). The number of haplotypes identified in the three western China native horse breeds ranged from nine to 15. The haplotype diversity and nucleotide diversity of the mtDNA D-loop region in three Chinese horse breeds are shown in Table 2. The haplotype diversity values ranged from 0.977±0.023 in the Kazakh breed to 0.989±0.031 in the Guizhou breed. The comparing result of nucleotide diversity was consistent with that of haplotype diversity. The data indicated an abundant genetic diversity in Chinese western domestic horse breeds.

Analysis of haplotype and nucleotide diversity of mtDNA D-loop region within different horse breeds

Phylogenetic analysis of different horse breeds based on mtDNA D-loop sequence

A Neighbor-Joining (NJ) phylogenetic tree was constructed based on 247 bp of mtDNA D-loop hypervariable region sequences of 43 Chinese western native horses (Figure 2). Based on the previously published seven major mtDNA lineages (A to G) in horses (Vila et al., 2001; Jansen et al., 2002; McGahern et al., 2006), the mtDNA haplotypes in Chinese three western horse breeds demonstrated the presence of these seven major lineages (A to G) (Figure 2). Ningqiang pony breed was distributed into two lineages (A, F), Guizhou breed had four lineages (A, B, D, G), and Kazakh breed had five lineages (A, C, D, E, F). Consistent with the front, when the NJ phylogenetic tree was constructed based on 600 bp of equine mtDNA D-loop sequences of 43 Chinese western native horses and 81 sequences of different horse breeds from GenBank (Table 1) it indicated that the three Chinese western horse breeds had distributed into the seven major lineages (A to G) (Figure 3). In detail, according to the NJ tree based on 600 bp of mtDNA D-loop sequences, Ningqiang pony breed was distributed into two lineages (A, F), Guizhou breed was distributed into three lineages (A, B, D, G), and Kazakh breed was distributed into four lineages (A, C, E, F). In addition, Baise breed was grouped into four lineages (A, C, D, G), Mongolian breed was grouped into four lineages (A, B, D, F), Debao pony breed was grouped into six lineages (A, B, C, D, F, G), Arabian breed was grouped into one lineage (A), Cheju breed was grouped into four lineages (A, B, C, F), Yunnan breed was grouped into five lineages (A, C, D, F, G). Among the three horse breeds and six breeds available from GenBank, lineage A was predominant. The results indicated that China western native horse breeds derived from multi-matrilineal origins and interchanged genetic material with each other widely during their evolutionary history.

Figure 2

Neighbour-Joining tree based on 247 bp of mtDNA D-loop region in three Chinese western horse breeds. Bootstrap values greater than 50 in 1,000 replications are shown at the respective nodes.

Figure 3

Neighbour-Joining tree based on 600 bp of mtDNA D-loop region in seven horse breeds in China (four breeds from GenBank) and two foreign breeds available from GenBank. Bootstrap values greater than 50 in 1,000 replications are shown at the respective nodes.

DISCUSSION

In recent years, many researches have indicated the diversity and origin of different horse populations, and suggested that domestic horses were derived from several wild populations (Lister, 1998; Vila., 2001; Jansen, 2002). Wallner (2003) studied the unique nucleotide sequence of the Y chromosome, and found that Przewalskii haplotype was the ancestor. Lei et al. (2009) showed that the domestication of horses was not a single discrete event, multiple horse domestication events may have occurred later in other regions, and that China was one of these domestication regions.

We found 33 haplotypes in 43 horses from three western China native horse breeds based on 600 bp of mtDNA control region sequence, and all these haplotypes grouped into seven lineages. These three horse breeds shared A lineage, and also had their respective lineages. Nine haplotypes were identified in ten Ningqiang ponies, and were distributed into two lineages (A, F). There were 13 haplotypes in 14 Guizhou horses representing three lineages (A, B, D, G). In 19 Kazakh horses, 15 haplotypes were distributed into four lineages (A, C, E, F). These data showed that the three western China native horses possessed abundant mtDNA diversity and indicated multiple maternal origins. The results were in accordance with other results from Yunnan, Mongolian, Debao, Guanzhong and other horse breeds in China (Wang et al., 1994; McGahern et al., 2006; Li et al., 2009). Some Ningqiang ponies, Guizhou horses and Kazakh horses were clustered into A lineage, indicating lineage A is predominant. Ningqiang pony breed demonstrated only A and F lineages, indicating close genetic relationship within this breed, which was consistent with the results of Lei et al. (2009). However, in the study of Lei et al. (2009), the Ningqiang breed indicated A, B, D and F lineages, Kazakh breed indicated A, B, C, D, F and G lineages, and Guizhou breed indicated A to G lineages. The possible reason for this difference with the present study is that more samples were analyzed in the study of Lei et al. (2009). It is necessary to point out that the lineages in the study of Lei et al were based on 247 bp while in the present study 600 bp D-loop sequences were used in these three horse breeds. In many other mammalian species, the variation in this short segment of the mitochondrial genome is often accompanied by high levels of recurrent mutations, thus blurring the structure of the tree and rendering the distinction between some important ancient branches within the tree virtually impossible (Achilli et al., 2009). So, to improve mtDNA phylogeny, the longer D-loop sequences in mitochondrial genomes will be helpful to achieve the highest level of molecular resolution in horses (Achilli et al., 2012).

In this study, three typical horse breeds were selected to clarify the origin, evolutionary relationship and exchange background of western horses in China. Kazakh horse is an ancient breed, and distributed in the Chinese northwest region, while the Guizhou horse is a southwest horse breed. Ningqiang pony is a breed from the border land between northwest and southwest and possesses some characteristics of northern and southern horse populations. Many researchers claim that Ningqiang pony belongs to the Southwest type. The recent results of archaeological studies indicated that there was a chama gudao (Tea-Horse road in ancient times) in southern Shaanxi and that Ningqiang county was an important center of the Tea-Horse trade system. Combined with the geography, history and the data from genetics, we believe that Ningqiang ponies were affected by northern and southern horse populations, such as the Mongolian, Guizhou and Kazakh breeds. Our data demonstrated that Chinese native horse breeds had extensive genetic flow from breeds all over the world. The results were in accordance with the historical records that horses were continuously imported for military forces and trade during Chinese history (Xie, 1987).

ACKNOWLEDGEMENTS

The work was supported by Scientific Research Fund of Education Department of Shaanxi Province (No. 06JK328 and 07JK212).

References

Achilli A, Bonfiglio S, Olivieri A, Malusa A, Pala M, Hooshiar Kashani B, Perego UA, Ajmone-Marsan P, Liotta L, Semino O, Bandelt HJ, Ferretti L, Torroni A. 2009;The multifaceted origin of taurine cattle reflected by the mitochondrial genome. PLoS One 4:e5753.
Achilli A, Olivieri A, Soares P, Lancioni H, Hooshiar Kashani B, Perego UA, Nergadze SG, Carossa V, Santagostino M, Capomaccio S, Felicetti M, Al-Achkar W, Penedo MC, Verini-Supplizi A, Houshmand M, Woodward SR, Semino O, Silvestrelli M, Giulotto E, Pereira L, Bandelt HJ, Torroni A. 2012;Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication. Proc Natl Acad Sci USA 109:2449–2454.
Asakawa S, Kumazawa Y, Araki T, Himeno H, Miura K, Watanabe K. 1991;Strand-specific nucleotide composition bias in echinoderm and vertebrate mitochondrial genomes. J Mol Evol 32:511–20.
Boore JL. 1999;Animal mitochondrial genomes. Nucleic Acids Res 27:1767–1780.
Bowling AT, Del Valle A, Bowling M. 2000;A pedigree-based study of mitochondrial D-loop DNA sequence variation among Arabian horses. Anim Genet 31:1–7.
Chang H. 2009;Animal genetic resources. Science Press
Du D, Deng L, Zhao CJ, Han GC, Wu CX. 2009;Study on the genetic diversity of Ningqiang pony and Mongolian horse using microsatellite markers. Chinese J Anim Sci 45:10–13.
Hou WT, Li GM, Wolstenholme DR. 1992;Genetic novelties in mitochondrial genomes of multicellular animals. Curr Opin Genet Dev 2:918–925.
Ishida N, Oyunsuren T, Mashima S, Mukoyama H, Saitou N. 1995;Mitochondrial DNA sequences of various species of the Genus Equus with special reference to the phylogenetic relationship between Przewalskii’s wild horse and domestic horse. J Mol Evol 41:180–188.
Jansen T, Forster P, Levine MA, Oelke H, Hurles M, Renfrew C, Weber J, Olek K. 2002;Mitochondrial DNA and the origins of the domestic horse. Proc Natl Acad Sci USA 99:10905–10910.
Ji X, Wu X, Yan P, Amato G. 2008;Complete sequence and gene organization of the mitochondrial genome of Siamensis Crocodile (Crocodylus siamensis). Mol Biol Rep 35:133–138.
Jiang Q, Wei Y, Huang Y, Jiang H, Guo Y, Lan G, Liao J. 2011;The complete mitochondrial genome and phylogenetic analysis of the Debao pony (Equus caballus). Mol Biol Rep 38:593–599.
Kakoi H, Tozaki T, Gawahara H. 2007;Molecular analysis using mitochondrial DNA and microsatellites to infer the formation process of Japanese native horse populations. Biochem Genet 45:375–395.
Ki JS, Hwang DS, Park TJ, Han SH, Lee JS. 2010;A comparative analysis of the complete mitochondrial genome of the Eurasian otter Lutra lutra (Carnivora; Mustelidae). Mol Biol Rep 37:1943–1955.
Kim KI, Yang YH, Lee SS, Park C, Ma R, Bouzat JL, Lewin HA. 1999;Phylogenetic relationships of Cheju horses to other horse breeds as determined by mtDNA D-loop sequence polymorphism. Anim Genet 30:102–108.
Lei CZ, Chen H, Yang GS, Sun WB, Lei XQ, Ge QL, Wang ZF, Lu N, Gao X, Hou WT. 2005;Study on mitochondrial DNA D-loop polymorphism in Chinese donkeys. Yi Chuan Xue Bao 32:481–486.
Lei CZ, Su R, Bower MA, Edwards CJ, Wang XB, Weining S, Liu L, Xie WM, Li F, Liu RY, Zhang YS, Zhang CM, Chen H. 2009;Multiple maternal origins of native modern and ancient horse populations in China. Anim Genet 40:933–944.
Li JL, Shi Y, Fan C, Manglai DJ. 2008;mtDNA diversity and original of Chinese Mongolian horse. Asian-Aust J Anim Sci 21:1696–1702.
Lister AM, Kadwell M, Kaagan LM, Jordan WC, Richards MB, Stanley HF. 1998;Ancient and modern DNA in study of horse domestication. Anc Biomol 2:267–280.
McGahern A, Bower MA, Edwards CJ, Brophy PO, Sulimova G, Zakharov I, Vizuete-Forster M, Levine M, Li S, MacHugh DE, Hill EW. 2006;Evidence for biogeographic patterning of mitochondrial DNA sequences in Eastern horse populations. Anim Genet 37:494–497.
Outram AK, Stear NA, Bendrey R, Olsen S, Kasparov A, Zaibert V, Thorpe N, Evershed RP. 2009;The earliest horse harnessing and milking. Science 323:1332–1335.
Peng R, Zeng B, Meng X, Yue B, Zhang Z, Zou F. 2007;The complete mitochondrial genome and phylogenetic analysis of the giant panda (Ailuropoda melanoleuca). Genes 397:76–83.
Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R. 2003;DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19:2496–2497.
Statistic Report in China Editing Commission. 2004. Statistic Report in China Statistic Report in China Editing Commission. Beijing, China:
Tamura K, Dudley J, Nei M, Kumar S. 2007;MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol 24:1596–1599.
Vila C, Leonard JA, Gotherstom A, Marklund S, Sandberg K, Linden K, Wayne RK, Ellegren H. 2001;Wide spread origin of domestic horse lineages. Science 291:474–477.
Wallner B, Brem G, Müller M, Achmann R. 2003;Fixed nucleotide differences on the Y chromosome indicate clear divergence between Equus przewalskii and Equus caballus. Anim Genet 34:453–456.
Wang W, Liu AH, Lin SY, Lan H, Su B, Xie DW, Shi LM. 1994;Multiple genotypes of mitochondrial DNA within a horse population from a small region in Yunnan Province of China. Biochem Genet 32:371–378.
Wolstenholme DR. 1992;Genetic novelties in mitochondrial genomes of multicellular animals. Curr Opin Genet Dev 2:918–925.
Xie CX. 1987. Horse and ass breeds in China Shanghai Scientific and Technical Publishing House. Shanghai, China:
Xu S, Luosang J, Hua S, He J, Ciren A, Wang W, Tong X, Liang Y, Wang J, Zheng X. 2007;High altitude adaptation and phylogenetic analysis of Tibetan horse based on the mitochondrial genome. J Genet Genomics 34:720–729.
Yang YH, Kim KI, Cothran EG, Flannery AR. 2002;Genetic diversity of Cheju horses (Equus caballus) determined by using mitochondrial DNA D-loop polymorphism. Biochem Genet 40:175–186.

Article information Continued

Figure 1

The polymorphic sites of 33 mtDNA D-loop haplotypes of three Chinese horse breeds.

Figure 2

Neighbour-Joining tree based on 247 bp of mtDNA D-loop region in three Chinese western horse breeds. Bootstrap values greater than 50 in 1,000 replications are shown at the respective nodes.

Figure 3

Neighbour-Joining tree based on 600 bp of mtDNA D-loop region in seven horse breeds in China (four breeds from GenBank) and two foreign breeds available from GenBank. Bootstrap values greater than 50 in 1,000 replications are shown at the respective nodes.

Table 1

Cited equine D- loop sequences from GenBank

Breed Abbreviation Number GenBank accession number
Arabian ARA 5 AY246180-AY 246183, AY246185
Baise BS 19 GQ203128-GQ203144, GQ222059-GQ222060
Cheju CJ 16 AF014405-AF014408, AF014411-AF014412, AF056071, AF354425, AF354427-AF354430, AY246204-AY246205, AY246207-AY246208
Debao DB 24 EU826536, EU939445, FJ392562-FJ392580, GQ203125-GQ203127
Mongolian MG 10 AF014413-014415, DQ297635, DQ297637, GU561995-GU561999
Yunnan YN 7 AF014416-AF014417, FJ718896-FJ718898, EF597512, EF597514
Donkey DK 1 NC001788

Table 2

Analysis of haplotype and nucleotide diversity of mtDNA D-loop region within different horse breeds

Breeds Samples Haplotypes Polymorphic sites Haplotype diversity Nucleotide diversity
Ningqiang pony 10 9 32 0.978±0.054 0.0199±0.0020
Guizhou 14 13 43 0.989±0.031 0.0226±0.0022
Kazakh 19 15 49 0.977±0.023 0.0244±0.0012