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.     

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).     

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.     

Genetic differentiation of water buffalo (Bubalus bubalis) populations in China, Nepal and south-east Asia: inferences on the region of domestication of the swamp buffalo

Data from three published studies of genetic variation at 18 microsatellite loci in water buffalo populations in China (18 swamp type, two river type), Nepal (one wild, one domestic river, one hybrid) and south-east Asia (eight swamp, three river) were combined so as to gain a broader understanding of genetic relationships among the populations and their demographic history. Mean numbers of alleles and expected heterozygosities were significantly different among populations. Estimates of θ (a measure of population differentiation) were significant among the swamp populations for all loci and among the river populations for most loci. Differentiation among the Chinese swamp populations (which was due primarily to just one population) was much less than among the south-east Asian. The Nepal wild animals, phenotypically swamp type but genetically like river type, are significantly different from all the domestic river populations and presumably represent the ancestral Bubalus arnee (possibly with some river-type introgression). Relationships among the swamp populations (D(A) genetic distances, principal component analysis and structure analyses) show the south-east Asian populations separated into two groups by the Chinese populations. Given these relationships and the patterns of genetic variability, we postulate that the swamp buffalo was domesticated in the region of the far south of China, northern Thailand and Indochina. Following domestication, it spread south through peninsular Malaysia to Sumatra, Java and Sulawesi, and north through China, and then to Taiwan, the Philippines and Borneo.     

Analysis of genetic diversity of domestic water buffaloes and anoas based on variations in the mitochondrial gene for cytochrome b

There are two major groups of domestic water buffaloes in East and Southeast Asia: swamp buffaloes and river buffaloes. Genetic diversity among swamp and river buffaloes was studied by DNA sequence analysis of the mitochondrial gene for cytochrome b. The results showed that each of the two groups has mitochondrial DNA (mtDNA) with a specific cytochrome b haplotype. The pairwise nucleotide sequence divergence was calculated to be 2·67% between swamp and river buffaloes, suggesting that they might have diverged from the ancestral populations of Asiatic domestic water buffaloes, approximately 1 million years ago. In addition, the sequences of the same gene from three subspecies of anoa (lowland, mountain and quarles anoa) were determined and compared with that of a domestic water buffalo. The sequence divergence was 1·2% for mountain anoa vs quarles anoa, 3·6% for mountain anoa vs lowland anoa and 3·3% for quarles anoa vs lowland anoa. Moreover, the sequence divergence between water buffaloes and anoas was found to be approximately 3·33%. Our results provide molecular evidence to support the taxonomic classification, namely, that Asiatic buffaloes may be classified into four lineages, swamp buffalo, river buffalo, lowland anoa and mountain plus quarles anoa. However, the sequence divergence values among these four groups were lower than the sequence divergence values found in the genus and subgenus levels within the subfamily Bovinae. In particular, in contrast to some proposed taxonomic classifications, our results indicated that mtDNA in the water buffaloes and anoas did not diverge at the genus level.     

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.     

Asian water buffalo: domestication,history and genetics

The domestic Asian water buffalo (Bubalus bubalis) is found on all five continents, with a global population of some 202 million. The livelihoods of more people depend on this species than on any other domestic animal. The two distinct types (river and swamp) descended from different wild Asian water buffalo (Bubalus arnee) populations that diverged some 900 kyr BP and then evolved in separate geographical regions. After domestication in the western region of the Indian subcontinent (ca. 6300 years BP), the river buffalo spread west as far as Egypt, the Balkans and Italy. Conversely, after domestication in the China/Indochina border region ca. 3000–7000 years BP, swamp buffaloes dispersed through south-east Asia and China as far as the Yangtze River valley. Molecular and morphological evidence indicates that swamp buffalo populations have strong geographic genetic differentiation and a lack of gene flow, but strong phenotypic uniformity. In contrast, river buffalo populations show a weaker phylogeographic structure, but higher phenotypic diversity (i.e. many breeds). The recent availability of a high-quality reference genome and of a medium-density marker panel for genotyping has triggered a number of genome-wide investigations on diversity, evolutionary history, production traits and functional elements. The growing molecular knowledge combined with breeding programmes should pave the way to improvements in production, environmental adaptation and disease resistance in water buffalo populations worldwide.     

Y‐chromosomal variation confirms independent domestications of swamp and river buffalo

Y‐chromosomal variation in the water buffalo was analysed by sequencing of DBY, ZFY and SRY gene segments. A clear separation of the paternal lineages of the river and swamp types parallels the differences between their maternal lineages and nuclear DNA. Sequence divergence was found to be comparable to the divergence of taurine cattle and zebu, and this divergence predated domestication, confirming that river and swamp buffalo originated from different wild populations. Within a sample of 23 Thai swamp buffaloes, we identified four haplotypes with different geographical distributions, two of which were shared by Thai wild buffaloes.     

Cytochrome b Genetic Diversity and Maternal Origin of Chinese Domestic Donkey

To clarify the origin of Chinese domestic donkeys, we investigated the mitochondrial Cytb gene from 244 animals from 13 native breeds. We found 55 variable sites in the Cytb gene sequence and subsequently defined 58 haplotypes. Analysis of haplotypes in combination with Cytb sequences revealed two mitochondrial origins in Chinese domestic donkeys, phenotypically expressed by the Somalian and Nubian lineages. The Somalian lineage predominated in Chinese domestic donkey breeds. Five specific Cytb gene SNPs diagnostic of each of the lineages were found in this study: 225(T-C), 237(C-T), 915(C-T), 1014(C-T), and 1134(A-G) mutations. They effectively distinguish the Nubian from the Somalian lineage in the mtDNA Cytb gene. Both lineages are from Africa and thus support the African maternal origins of Chinese domestic donkeys. No obvious geographic structure was found in Chinese domestic donkey breeds, but the population showed abundant genetic diversity.     

Strong and stable geographic differentiation of swamp buffalo maternal and paternal lineages indicates domestication in the China/Indochina border region

The swamp type of the Asian water buffalo is assumed to have been domesticated by about 4000 years BP, following the introduction of rice cultivation. Previous localizations of the domestication site were based on mitochondrial DNA (mtDNA) variation within China, accounting only for the maternal lineage. We carried out a comprehensive sampling of China, Taiwan, Vietnam, Laos, Thailand, Nepal and Bangladesh and sequenced the mtDNA Cytochrome b gene and control region and the Y‐chromosomal ZFY, SRY and DBY sequences. Swamp buffalo has a higher diversity of both maternal and paternal lineages than river buffalo, with also a remarkable contrast between a weak phylogeographic structure of river buffalo and a strong geographic differentiation of swamp buffalo. The highest diversity of the swamp buffalo maternal lineages was found in south China and north Indochina on both banks of the Mekong River, while the highest diversity in paternal lineages was in the China/Indochina border region. We propose that domestication in this region was later followed by introgressive capture of wild cows west of the Mekong. Migration to the north followed the Yangtze valley as well as a more eastern route, but also involved translocations of both cows and bulls over large distances with a minor influence of river buffaloes in recent decades. Bayesian analyses of various migration models also supported domestication in the China/Indochina border region. Coalescence analysis yielded consistent estimates for the expansion of the major swamp buffalo haplogroups with a credibility interval of 900 to 3900 years BP. The spatial differentiation of mtDNA and Y‐chromosomal haplotype distributions indicates a lack of gene flow between established populations that is unprecedented in livestock.     

Genetic variants of ribosomal DNA and mitochondrial DNA between swamp and river buffaloes

To clarify the genetic relationship between Swamp and River buffaloes, the restriction fragment length polymorphisms (RFLPs) of nuclear genomic ribosomal DNA (rDNA) and cytoplasmic mitochondrial DNA (mtDNA) were analysed. Blood or liver samples from 73 Swamp and three River buffaloes were collected in East and South-east Asian countries. DNA samples from cattle, goats and sheep were used for comparisons. The analysis of rDNA allowed water buffaloes, cattle, goats and sheep to be characterized by four distinct repeat-types. However, swamp and river buffaloes showed the same repeat-type. Divergence of water buffalo and cattle is considered to have occurred approximately four to six million years ago. The RFLPs for mtDNA divided water buffaloes into three haplotypes, swamp-1, swamp-2 and river types. Swamp-1 accounted for 91% of all swamp buffaloes while swamp-2 was observed only in water buffaloes from Thailand (9%). All river buffaloes were of the same haplotype. No differences were observed between swamp and river buffaloes at the rDNA level. In contrast, a few distinct differences between them were found at the mtDNA level. Therefore, mtDNA polymorphisms provide an adequate means for classifying water buffaloes into either swamp or river buffaloes.     

Mitochondrial DNA analysis provides new insights into the origin of the Chinese domestic goat

Many studies have focused on modern goats, however, few reports focused on origin and genetic structure of Chinese ancient goats. In this study, we analyzed the 289-bp fragment of mitochondrial DNA (mtDNA) control region from 14 Chinese ancient goats excavated from two archaeological sites in Inner Mongolia, China dating back about 2500 years. 10 haplotypes were successfully obtained from the 14 ancient goats. Phylogenetic analysis revealed the multiple maternal origins of Chinese domestic goats, three mtDNA lineages A, B, and D were identified in the Chinese ancient individuals, in which lineage A was predominant (70%), lineages B was moderate (20%), and lineage D was present at low frequency (10%). The network analysis showed that lineage B was subdivided into two subgroups B1 and B2. One of the Chinese ancient goats shared the founder haplotype in the center of subgroup B1, and the shared sequences of the founder haplotypes of subgroups B1 and B2 distributed mainly in China. These results implied that lineage B including subgroups B1 and B2 probably originated from China, and further supported the hypothesis that China may be one of the goat domestication centers. In addition, the analysis of shared sequences indicated that the ancient goats from Inner Mongolia were closely genetically related to Chinese modern goats, suggesting that the ancient goats from Inner Mongolia had the genetic contribution to Chinese modern goats.     

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.     

Origin, genetic diversity, and population structure of Chinese domestic sheep

To characterize the origin, genetic diversity, and phylogeographic structure of Chinese domestic sheep, we here analyzed a 531-bp fragment of mtDNA control region of 449 Chinese autochthonous sheep from 19 breeds/populations from 13 geographic regions, together with previously reported 44 sequences from Chinese indigenous sheep. Phylogenetic analysis showed that all three previously defined lineages A, B, and C were found in all sampled Chinese sheep populations, except for the absence of lineage C in four populations. Network profiles revealed that the lineages B and C displayed a star-like phylogeny with the founder haplotype in the centre, and that two star-like subclades with two founder haplotypes were identified in lineage A. The pattern of genetic variation in lineage A, together with the divergence time between the two central founder haplotypes suggested that two independent domestication events have occurred in sheep lineage A. Considerable mitochondrial diversity was observed in Chinese sheep. Weak structuring was observed either among Chinese indigenous sheep populations or between Asian and European sheep and this can be attributable to long-term strong gene flow induced by historical human movements. The high levels of intra-population diversity in Chinese sheep and the weak phylogeographic structuring indicated three geographically independent domestication events have occurred and the domestication place was not only confined to the Near East, but also occurred in other regions.     

Genetic diversity and differentiation of Chinese domestic buffalo based on 30 microsatellite markers

To determine genetic diversity and evolutionary relationships among Chinese buffalo populations, 18 indigenous swamp buffalo populations and two introduced river buffalo breeds were genotyped for 30 microsatellite loci. The mean number of alleles across the 30 loci was 8.13, and the expected heterozygosity ranged from 0.517 (Yanjin) to 0.609 (Dehong). Although there was abundant genetic variation, genetic differentiation between Chinese buffalo populations was low, with only 2.8% of the total genetic variance among populations. The genetic differentiation pattern and genetic relationships among Chinese buffalo populations were consistent with their geographic distribution. The Dehong population was discerned as a distinct indigenous population, but suffered genetic admixture with river buffalo. The remaining populations were divided into four major clusters, i.e. the Upper and Middle Reaches of Yangtze Valley cluster (Guizhou, Guizhoubai, Yanjin, Fuling, Enshi and Jianghan), the Lower Reaches of Yangtze Valley cluster (Haizi, Shanqu and Dongliu), the South China cluster (Fuan and Xinfeng) and the Southwest China cluster (Xinglong, Xilin, Diandongnan and Dechang).     

Mitochondrial diversity and phylogeographic structure of Chinese domestic goats

China has numerous native domestic goat breeds, but so far there has been no extensive study on genetic diversity, population demographic history, and origin of Chinese goats. Here, we examined the genetic diversity and phylogeographic structure of Chinese domestic goats by determining a 481-bp fragment of the first hypervariable region of mitochondrial DNA (mtDNA) control region from 368 individuals representing 18 indigenous breeds. Phylogenetic analyses revealed that there were four mtDNA lineages (A-D) identified in Chinese goats, in which lineage A was predominant, lineage B was moderate, and lineages C and D were at low frequency. These results further support the multiple maternal origins of domestic goats. The pattern of genetic variation in goat mtDNA sequences indicated that the two larger lineages A and B had undergone population expansion events. In a combined analysis of previously reported sequences and our sequences belonging to lineage B, we detected two subclades, in which one was unique to eastern Asia and another was shared between eastern and southern Asia. A larger genetic variation in eastern Asia than southern Asia and the pattern of phylogeographic variation in lineage B suggest that at least one subclade of lineage B originated from eastern Asia. There was no significant geographical structuring in Chinese goat populations, which suggested that there existed strong gene flow among goat populations caused by extensive transportation of goats in history.     

Genetic variation within and relationships among populations of Asian water buffalo (Bubalus bubalis)

Genetic variation at 53 protein-coding loci (25 polymorphic) was analysed for 17 water buffalo populations – 12 swamp, three Lankan and two of the Murrah breed (river type), to determine the magnitude of genetic differentiation and the genetic relationships among the populations. In accord with previous cytological studies, the Lankan buffalo clearly are river type. Significant deviations from Hardy–Weinberg equilibrium were shown for a number of locus–population combinations, with all populations but one showing significant heterogeneity in these deviations among loci. By contrast, heterogeneity among populations for each locus was much less, indicating locus-specific deviations, which suggest selection affecting allele frequencies at some loci. There was significant genetic differentiation among populations of both the swamp and river types. The differentiation among the swamp populations may reflect the geography of south-east Asia and the presumed spread of the swamp buffalo through this region. Phylogenies derived from pairwise genetic distance estimates show the clear separation of swamp and river types, but the topology of the swamp populations shows rather poor consistency with their geographic locations. For at least one population (Australia), it is clear that bottleneck effects have distorted the phylogenetic topology. Average genetic distances for both the swamp and river types, as compared with previous studies of livestock breeds, show that the genetic differentiation of each of these sets of populations is of the same order of magnitude as that among well-recognized and established breeds of other species.     

Phylogenetic relationship among all living species of the genusBubalus based on DNA sequences of the cytochromeb gene

The cytochromeb genes of all living species ofBubalus, including the river type and the swamp type of domestic buffaloes (Bubalus bubalis), were sequenced to clarify their phylogenetic relationships. These sequences were compared together with the African buffalo (Syncerus caffer) and banteng (Bos javanicus) sequences as an outgroup. Phylogenetic trees ofBubalus species based on the DNA sequences of the cytochromeb gene demonstrated that the tamaraw (Bubalus mindorensis), endemic to the Philippines, could be classified into the subgenusBubalus, not the subgenusAnoa. The divergence time between the lowland anoa (B. depressicornis) and the mountain anoa (B. quarlesi) was estimated at approximately 2.0 million years (Myr), which is almost the same as the coalescence time for theBubalus sequences. This large genetic distance supports the idea that the lowland anoa and the mountain anoa are different species. An unexpectedly large genetic distance between the river and the swamp type of domestic buffaloes suggests a divergence time of about 1.7 Myr, while the swamp type was noticed to have the closest relationship with the tamaraw (1.5 Myr). This result implies that the two types of domestic buffaloes have differentiated at the full species level.     

Genetic diversity of mitochondrial cytochrome b gene in swamp buffalo

The swamp buffalo (Bubalus carabanensis) is mainly bred for meat, transport and rice cultivation in China and Southeast Asian countries. In the current study, we investigated the genetic diversity, maternal origin and phylogenetic relationship of swamp buffalo by analyzing 1,786 mitochondrial cytochrome b (cytb) sequences from China, Vietnam, Laos, Thailand, India and Bangladesh. Our results indicated that the swamp buffalo can be divided into two major lineages (SA and SB with the sublineages) and three rare lineages (SC, SD and SE), which showed strong geographic differentiation. The SA1 lineage represented a major domestication event, which involved population expansion. Regions III and V showed higher haplotype diversity than the other regions, indicating that the regions of Southwest China and IndoChina are potential domestication centers for the swamp buffalo. In addition, the swamp buffalo showed a closer phylogenetic relationship to tamaraw. In conclusion, our findings revealed a high level of genetic diversity and the phylogenetic pattern of the swamp buffalo.