Analysis of mitochondrial D-loop region casts new light on domestic water buffalo (Bubalus bubalis) phylogeny

The phylogeny of water buffaloes (Bubalus bubalis) is still a matter of discussion, especially if the two types of domestic water buffalo (swamp and river) derived from different domestication events or if they are products of human selection. To obtain more insight, we analyzed the entire mitochondrial D-loop region of 80 water buffaloes of four different breeds, i.e., 19 swamp buffaloes (Carabao) and 61 river buffaloes (Murrah, Jafarabadi, and Mediterranean), sampled in Brazil and Italy. We detected 36 mitochondrial haplotypes with 128 polymorphic sites. Pooled with published data of South-East Asian and Australian water buffaloes and based on comprehensive median-joining network and population demography analyses we show evidence that both river and swamp buffaloes decent from one domestication event, probably in the Indian subcontinent. However, the today swamp buffaloes have an unravelled mitochondrial history, which can be explained by introgression of wild water buffalo mtDNA into domestic stocks. We are also discussing indications for an independent domestication of buffaloes in China.     

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.     

Mitochondrial DNA diversity of Pakistani goats

Pakistan contains numerous domestic goat breeds, but until now there has been no comprehensive study on genetic diversity or a phylogenetic analysis of Pakistani goats. In this study, we analysed the complete mitochondrial DNA D-loop and the cytochrome b gene of 13 Pakistani domestic goat breeds (Capra hircus) and one wild goat, the Sindh Ibex (Capra aegagrus blythi). The phylogenetic analyses and sequence divergence (SD) established four distinct mt-lineages termed as A, B and C (previously reported) and a new lineage D. The Sindh Ibex appeared as an outgroup of domestic goats. The estimated divergence times between the most recently evolved mt-lineages A and D were from 260,483 to 371,052 YA. This suggested that at least four different strains of wild Capra might have been the source of the modern domestic goats. The new mt-lineage D revealed high SD from mt-lineage A and may be the oldest branch under domestication, while mt-lineages B and C showed lower SD and might have been domesticated during an advanced stage of the domestication process.     

Mitochondrial sequence-based evolutionary analysis of riverine–swamp hybrid buffaloes of India indicates novel maternal differentiation and domestication patterns

In this study, mitochondrial D-loop sequence data on riverine, swamp and hybrid buffaloes from India have been generated and compared with other reported Indian riverine, Chinese and Bangladeshi swamp buffalo populations. Sequence analysis revealed the presence of 132 haplotypes, with a haplotype diversity of 0.9611 ± 0.0045 and a nucleotide diversity of 0.04801 ± 0.00126. For the first time, the existence of riverine–swamp hybrids among the Indian Chilika buffalo population has been recorded, having 49 chromosomes, which was also confirmed by mitochondrial haplotype sharing between Chilika and Indian swamp as well as Chinese swamp buffalo populations in the network analysis. Phylogenetic analysis documents the sharing of reported pre-domestication haplogroups ‘SA1’, ‘SA2’, ‘SA3’ and ‘SB1’ between the Chilika and swamp buffalo populations of India, China and Bangladesh, an indication of the migration of swamp buffaloes towards Bangladesh and adjoining lower parts of India and north towards Chinese domestication sites. The results have also been supplemented by multidimension scaling, grouping Indian and Chinese swamp buffaloes more closely together with Bangladeshi buffaloes, but into a separate quadrant, whereas Chilika grouped away from other riverine as well as swamp buffaloes. These findings thus confirm the previous reports that the northeast region of India, close to the Indo-China border, is the point of evolution of swamp buffaloes with multiple sites of domestication.     

Phylogeography and domestication of Indian river buffalo

Background  

The water buffalo- Bu balus bubalis holds tremendous potential in livestock sector in many Asian countries, particularly India. The origin, domestication and genetic structure of the Indian river buffalo are poorly understood. Therefore, to understand the relationship among the maternal lineages of Indian river buffalo breeds and their domestication process, we analysed mitochondrial D-loop region of 217 animals representing eight breeds from eight different locations in India along with published sequences of Mediterranean buffalo.     

Genetic diversity of Asian water buffalo (Bubalus bubalis): mitochondrial DNA D-loop and cytochrome b sequence variation

Swamp and river buffalo mitochondrial DNA (mtDNA) was sequenced for 303 bp of the cytochrome b gene for 54 animals from 14 populations, and for 158 bp of the D-loop region for 80 animals from 11 populations. Only one cytochrome b haplotype was found in river buffalo. Of the four haplotypes identified in swamp buffalo, one found in all populations is apparently ancestral both to the other swamp haplotypes and to the river haplotype. The phylogenetic relationships among the 33 D-loop haplotypes, with a cluster of 11 found in swamp buffalo only, also support the evolution of domesticated swamp and river buffalo from an ancestral swamp-like animal, most likely represented today by the wild Asian buffalo ( Bubalus arnee ). The time of divergence of the swamp and river types, estimated from the D-loop data, is 28 000 to 87 000 years ago. We hypothesise that the species originated in mainland south-east Asia, and that it spread north to China and west to the Indian subcontinent, where the rive type evolved and was domesticated. Following domestication in China, the domesticated swamp buffalo spread through two separate routes, through Taiwan and the Philippines to the eastern islands of Borneo and Sulawesi, and south through mainland southeast Asia and then to the western islands of Indonesia.     

贵州黄牛mtDNA D-loop 遗传多样性研究

对贵州4个地方黄牛品种共计82个个体的线粒体DNA D-loop区全序列910 bp进行分析,检测到31种单倍型,其核苷酸多态位点65个,约占所测核苷酸总长的7.14%,其中有62个转换,2个颠换,1个转换/颠换共存。贵州4个黄牛品种mtDNA D-loop区核苷酸多样度（π值）为2.16%～2.61%,单倍型多样度（H）为0.695～0.909,表明贵州黄牛mtDNA遗传多样性比较丰富。根据单倍型构建了贵州4个黄牛品种的NJ分子系统树。聚类表明,贵州黄牛有普通牛和瘤牛2大母系起源,其影响较为均一。并探讨了用核苷酸多样度π值的大小来衡量黄牛群体遗传分化程度的可行性。      

Independent maternal origin of Chinese swamp buffalo (Bubalus bubalis)

To obtain more knowledge on the origin and genetic diversity of the swamp buffalo (Bubalus bubalis) in China, the complete mitochondrial D-loop sequences of 119 samples representing seven native types were compared. Two mitochondrial DNA (mtDNA) lineages (lineages A and B) were determined for the Chinese swamp buffalo. Examination of the diversity patterns suggest that lineage A has undergone a population expansion event. Divergence of lineages A and B was estimated at 18,000 years ago. Combined analyses of mtDNA sequences from Chinese, Indian, Brazilian/Italian and Southeast Asian/Australian buffalo samples showed independent domestication events in the swamp buffalo from China and the river buffalo from the India subcontinent. The spread of swamp and river buffalo from China and India respectively to mainland Southeast Asia suggests that Southeast Asia is a hybrid zone for buffalo. Our data support the hypothesis of the evolution of domesticated swamp and river buffalo from ancestral swamp-like animals. These ancestral animals were extensively distributed across mainland Asia and most likely are represented today by the wild Asian buffalo (Bubalus arnee).     

DNA fingerprinting of Pakistani buffalo breeds (Nili-Ravi, Kundi) using microsatellite and cytochrome b gene markers

Cytochrome b gene markers have been proved as an efficient and powerful tool for breed characterization and species identification of buffaloes. This study represents the substantial analysis of mitochondrial DNA variation in Pakistani buffalo breeds and provides information about their genetic diversity. In this study partial amplification of cytochrome b gene of 1,061 bp was done and sequencing results showed ten haplotypes. Comparing all fifty samples from two buffalo breeds of Pakistan, fifteen polymorphic sites were observed out of which, twelve codons 42, 71, 118, 120, 199, 235, 269, 297, 318, 327, 350, 355 of mitochondrial cytochrome b gene are monomorphic which translate same amino acids as in the reference protein sequence due to silent mutation while different in DNA sequence. Similarly three codons 163, 246, 337 of mitochondrial cytochrome b are polymorphic and different from the reference sequence with respect to DNA as well as protein sequence. For the further confirmation a panel of nine microsatellite markers was used with high polymorphism information content (PIC). The frequency distribution of these alleles varies from three to eight allele at locus CSSM66 and ILST029 respectively. The results obtained from this study may contribute to the establishment of routine genotyping service of buffalo breeds for buffalo farmers for animal forensic application in case of any dispute. Additionally this study may help for breed characterization and phylogeny of aforementioned breeds of buffalo.     

Near East mtDNA haplotype variants in Roman cattle from Augusta Raurica, Switzerland, and in the Swiss Evolène breed

Typical Near East mitochondrial haplotypes of the T2 lineage were found in one cattle metacarpus sample from the Roman period and in two present-day Evolène cattle in Switzerland. Sequences from eight additional Evolène and four Raetian Grey aligned to the European haplotype T3. Analysis of nucleotide diversity within the mitochondrial D-loop of both studied Swiss cattle breeds revealed high haplotype diversity and similar diversity to a European cattle reference group. Mitochondrial T3 haplotypes radiated star-like from two similarly frequent haplotypes, possibly indicating two different expansion routes. The breed structure of Evolène cattle can be explained either by an introduction of diverse female lineages from the domestication centre or by later admixture. The introduction of the Near East lineage to Switzerland must have happened during the Roman time or earlier.     

Phylogeography and Domestication of Chinese Swamp Buffalo

To further probe into whether swamp buffaloes were domesticated once or multiple times in China, this survey examined the mitochondrial DNA (mtDNA) Control Region (D-loop) diversity of 471 individuals representing 22 populations of 455 Chinese swamp buffaloes and 16 river buffaloes. Phylogenetic analysis revealed that Chinese swamp buffaloes could be divided into two distinct lineages, A and B, which were defined previously. Of the two lineages, lineage A was predominant across all populations. For predominant lineage A, Southwestern buffalo populations possess the highest genetic diversity among the three hypothesized domestication centers (Southeastern, Central, and Southwestern China), suggesting Southwestern China as the most likely location for the domestication of lineage A. However, a complex pattern of diversity is detected for the lineage B, preventing the unambiguous pinpointing of the exact place of domestication center and suggesting the presence of a long-term, strong gene flow among swamp buffalo populations caused by extensive migrations of buffaloes and frequent human movements along the Yangtze River throughout history. Our current study suggests that Southwestern China is the most likely domestication center for lineage A, and may have been a primary center of swamp buffalo domestication. More archaeological and genetic evidence is needed to show the process of domestication.     

Mitochondrial DNA analyses of Indian water buffalo support a distinct genetic origin of river and swamp buffalo

Water buffalo (Bubalus bubalis) is broadly classified into river and swamp categories, but it remains disputed whether these two types were independently domesticated, or if they are the result of a single domestication event. In this study, we sequenced the mitochondrial D-loop region and cytochrome b gene of 217 and 80 buffalo respectively from eight breeds/locations in northern, north-western, central and southern India and compared our results with published Mediterranean and swamp buffalo sequences. Using these data, river and swamp buffalo were distinguished into two distinct clades. Based upon the existing knowledge of cytogenetic, ecological and phenotypic parameters, molecular data and present-day distribution of the river and swamp buffalo, we suggest that these two types were domesticated independently, and that classification of the river and swamp buffalo as two related subspecies is more appropriate.     

我国部分黄牛品种线粒体D-loop区遗传多样性与起源分化

为了解我国地方黄牛品种线粒体DNA的遗传变异情况, 文章测定了16个地方黄牛品种206个个体线粒体D-loop区的全序列, 共检测到101个变异位点; 99种单倍型, 其中73种是普通牛单倍型, 26种是瘤牛单倍型; 平均核苷酸差异为22.6920, 单倍型多样度为0.9320, 核苷酸多样度为0.0227, 表明我国黄牛品种遗传多样性非常丰富。构建的NJ进化树显示16个品种来源于两大母系: 普通牛和瘤牛; 构建的Network图表明73种普通牛单倍型可以分为3大单倍型群; 26种瘤牛单倍型分为5种单倍型群, 推测我国瘤牛在迁移过程中, 至少经历了4次群体扩张事件。通过分析比较地方黄牛品种与内罗门牛共有的 H3单倍型, 发现其中只有16%的序列与内罗门牛的H3单倍型非常相似, 其余84%的序列均发生了鸟嘌呤变异, 推测这些变异很可能是我国瘤牛固有的变异。     

Genetic diversity of Guangxi chicken breeds assessed with microsatellites and the mitochondrial DNA D-loop region

The domestic chicken (Gallus gallus domesticus) is an excellent model for genetic studies of phenotypic diversity. The Guangxi Region of China possesses several native chicken breeds displaying a broad range of phenotypes well adapted to the extreme hot-and-wet environments in the region. We thus evaluated the genetic diversity and relationships among six native chicken populations of the Guangxi region and also evaluated two commercial breeds (Arbor Acres and Roman chickens). We analyzed the sequences of the D-loop region of the mitochondrial DNA (mtDNA) and 18 microsatellite loci of 280 blood samples from six Guangxi native chicken breeds and from Arbor Acres and Roman chickens, and used the neighbor-joining method to construct the phylogenetic tree of these eight breeds. Our results showed that the genetic diversity of Guangxi native breeds was relatively rich. The phylogenetic tree using the unweighed pair-group method with arithmetic means (UPGAM) on microsatellite marks revealed two main clusters. Arbor Acres chicken and Roman chicken were in one cluster, while the Guangxi breeds were in the other cluster. Moreover, the UPGAM tree of Guangxi native breeds based on microsatellite loci was more consistent with the genesis, breeding history, differentiation and location than the mtDNA D-loop region. STRUCTURE analysis further confirmed the genetic structure of Guangxi native breeds in the Neighbor-Net dendrogram. The nomenclature of mtDNA sequence polymorphisms suggests that the Guangxi native chickens are distributed across four clades, but most of them are clustered in two main clades (B and E), with the other haplotypes within the clades A and C. The Guangxi native breeds revealed abundant genetic diversity not only on microsatellite loci but also on mtDNA D-loop region, and contained multiple maternal lineages, including one from China and another from Europe or the Middle East.     

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.     

Molecular phylogeny of Pakistani riverine buffalo based on genetic variability of mitochondrial cytochrome b gene

Mitochondrial cytochrome b gene is considered to be one of the best markers for breed characterization as well as studying the ancestry in the vertebrates due to its exclusive maternal inheritance. DNA fingerprinting by single nucleotide polymorphism is most reliable and widely used molecular technique in modern forensics and is being considered in this study. Partial sequencing of 1,061?bp of aforementioned gene from 14580 to 15643 was conducted in two famous Pakistani buffalo breeds named Nili-Ravi and Kundi. In which we explore seven haplotypes within earlier and none in the latter breed. Nili-Ravi is polymorphic at four codons of this gene, and the protein translation is also different from the reference sample while monomorphic at three codons with no amino acid replacement. Haplotypes frequency distribution of these four haplotypes named NR3, NR4, NR5, NR7 revealed that the prevalence of each haplotype is 0.04?% in the Pakistani buffalo population of this Nili-Ravi breed while complete homoplasmy was observed in the Kundi breed population. Nili-Ravi breed of buffalo is genetically more variable than the Kundi breed as far as the gene in subject is concerned. It means later breed has spent more time to propagate its wild type haplotype which make this breed more ancestral as compare to Nili-Ravi. Secondly both breeds share their common ancestors with regional water buffalo rather than the swamp one.     

Introgression and isolation contributed to the development of Hungarian Mangalica pigs from a particular European ancient bloodline

Background

Mangalica breeds are indigenous to Hungary and their breeding history dates back to about 200–250 years ago. They are fat-type pigs and have a rare curly hair phenotype. The aim of our study was to establish the relationships between these unique breeds and other European breeds.

Results

Based on a core sequence of 382 bp present in 2713 mitochondrial D-loop sequences from pigs belonging to 38 local breeds from nine countries, five cosmopolitan breeds and wild boars from 14 countries, we identified 164 haplotypes. More than half of the 2713 sequences belonged to either four haplotypes characteristic of continental European breeds or two haplotypes characteristic of British/cosmopolitan breeds; each haplotype is present in more than 100 individuals. Most Mangalica individuals belonged either to one of these common continental European haplotypes or to two Mangalica-specific haplotypes that were absent in all other breeds. In addition, we identified the ancestral mitochondrial D-loop signature present in these 2713 sequences and found that ~ 80% carried the European ancient signatures, ANC-Aside and ANC-Cside or their closely related signatures, while most of the remaining sequences carried a modern Asian signature, ANC-Easia. Mangalica individuals carried the ANC-Aside signature, but not the ANC-Cside or ANC-Easia signatures.

Conclusions

In all the Mangalica individuals, a unique ancient European signature was found in the mitochondrial DNA D-loop region, but they belonged almost exclusively to either certain very abundant European or two Mangalica-specific D-loop haplotypes. This indicates that the present-day Mangalica population in Hungary evolved either by introgression of other European breeds and wild boars or via total isolation after the divergence of European ancient porcine bloodlines.     

Five ovine mitochondrial lineages identified from sheep breeds of the near East

Archaeozoological evidence indicates that sheep were first domesticated in the Fertile Crescent. To search for DNA sequence diversity arising from previously undetected domestication events, this survey examined nine breeds of sheep from modern-day Turkey and Israel. A total of 2027 bp of mitochondrial DNA (mtDNA) sequence from 197 sheep revealed a total of 85 haplotypes and a high level of genetic diversity. Six individuals carried three haplotypes, which clustered separately from the known ovine mtDNA lineages A, B, and C. Analysis of genetic distance, mismatch distribution, and comparisons with wild sheep confirmed that these represent two additional mtDNA lineages denoted D and E. The two haplogroup E sequences were found to link the previously identified groups A and C. The single haplogroup D sequence branched with the eastern mouflon (Ovis orientalis), urial (O. vignei), and argali (O. ammon) sheep. High sequence diversity (K = 1.86%, haplogroup D and O. orientalis) indicates that the wild progenitor of this domestic lineage remains unresolved. The identification in this study of evidence for additional domestication events adds to the emerging view that sheep were recruited from wild populations multiple times in the same way as for other livestock species such as goat, cattle, and pig.     

Mitochondrial D-loop sequence variation among the 16 maternal lines of the Lipizzan horse breed

Mitochondrial DNA from 49 Lipizzan horses representing 16 maternal lines from the original stud at Lipica was used for SSCP analysis and DNA sequencing. The SSCP analysis of the 444 bp long fragment of the D-loop region extending from the tRNA(Pro) gene to the central conserved sequence block revealed three distinct groups of SSCP patterns. Both ends of the D-loop region (378 bp and 310 bp), which are considered as the most variable regions within the mammalian mitochondrial DNA, were sequenced. According to 49 polymorphic sites identified within the both parts of the D-loop region, the 16 maternal lines were grouped into 13 distinct mitochondrial haplotypes. The minimal difference between two different haplotype DNA sequences was one nucleotide and the maximal 24 nucleotides. The inheritance of mitochondrial haplotypes was stable and no sequence variation potentially attributable to mutation within maternal line was observed. Considerable DNA sequence similarity of Lipizzan mitochondrial haplotypes with the haplotypes from other breeds was observed. Phylogenetic analysis of the sequence data revealed a dendrogram with three separated branches, supporting the historical data about the multiple origin of the Lipizzan breed.     

中国山羊mtDNA D-loop遗传多样性及其起源研究

采用DNA测序技术分析了中国9个山羊品种（板角山羊、成都麻羊、贵州黑山羊、贵州白山羊、黔北麻羊、马头山羊、陕南白山羊、黄淮山羊和雷州山羊）共计128个个体的mtDNA D-loop全序列。结果表明：山羊mtDNA D-loop全序列长度为1212-1213bp,检测到102个变异位点,约占分析位点总数的8．42％,可变位点中转换占99个,颠换2个,1个转换/颠换共存;界定了92种单倍型,有78种为各品种独享单倍型,另外14种为群体内或群体间共享单倍型。9个山羊品种单倍型多样度为0．9333-1．0000,核苷酸多样度为0．7062％-1．8265％,表明中国山羊品种遗传多样性丰富。根据92种mtDNA单倍型构建了中国山羊的NJ分子系统树,聚类表明,中国山羊mtDNA D-loop序列单倍型分为支系A和支系B两大类。支系A包括75种单倍型,代表95个样本,占总数的74．22％;支系B包括17种单倍型,代表33个样本,占总数的25．78％,说明中国山羊存在支系A和支系B两大母系起源。对中国山羊mtDNA D-loop的支系A和支系B进行核苷酸不配对分布曲线分析和Fu的Fs中性检验,分析表明,支系A的分布曲线呈单峰形,Fs值为-24．6491,P值为0．0000,显著偏离中性,表明山羊支系A曾经历群体扩张;支系B呈近似双峰分布,Fs值为-3．3947,P值为0．0980,中性检验差异不显著,表明山羊支系B没有经历群体扩张,群体大小保持相对稳定。山羊支系B可能起源于中国。