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.     

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.     

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

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.     

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.     

Genetic diversity of mitochondrial cytochrome b gene in Chinese native buffalo

The origins of the domestic water buffalo remain contentious. To better understand the origins of Chinese water buffalo, we sequenced the complete mitochondrial cytochrome b (MT-CYB) gene from 270 individuals representing 13 Chinese domestic swamp buffalo populations. We found genetic evidence of introgression of river buffalo into Chinese swamp buffalo herds. Swamp buffalo haplotypes can be divided into two highly divergent lineages (A and B), suggesting that Chinese native swamp buffalo have two maternal origins. We found that the A→G transition in the buffalo MT-CYB gene stop codon resulted in buffalo haplotypes being terminated by one of two stop codons: AGA or AGG. AGA is common to river buffalo and lineage A of swamp buffalo, while AGG is specific to lineage B of swamp buffalo. Lineage A appears to have been domesticated in China. Further genetic evidence is required to clarify the origins of lineage B.     

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.     

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.     

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.     

Molecular phylogeography and genetic diversity of East Asian goats

The domestic goat is one of the most important livestock species, but its origins and genetic diversity still remain uncertain. Multiple highly divergent maternal lineages of goat have been reported in previous studies. Although one of the mitochondrial DNA lineages, lineage B, was detected only in eastern and southern Asia, the geographic distribution of these lineages was previously unclear. Here, we examine the genetic diversity and phylogeographic structure of Asian goats by mitochondrial DNA sequences and morphological characteristics. The analyses of a total of 1661 Asian goats from 12 countries revealed a high frequency of lineage B in Southeast Asia. The frequency of this lineage tended to be higher in mountain areas than in plain areas in Southeast Asian countries, and there was a significant correlation between its frequency and morphological traits. The results suggest an original predominance of lineage B in Southeast Asia and the recent infiltration of lineage A into Southeast Asian goats.     

Mitochondrial DNA diversity patterns in Pakistani buffalo

The Indian subcontinent is considered to be the likely centre of river buffalo domestication, based on population dynamics, archaeological evidence and genetic diversity. Recent studies on mitochondrial DNA diversity have drawn useful conclusions about the domestication history of Bubalus bubalis. The conclusions of these studies are, however, incomplete, unless samples can also be analysed from Pakistan, which contains the second largest buffalo population of the world. Here, we report the results of the first study on mitochondrial D-loop sequence diversity in five breeds of Pakistani buffalo. Analysis of sequence variations in 503-bp of the D-loop region of 123 animals revealed 52 haplotypes, including 40 singletons. Multidimensional display of breed pairwise F(ST) values revealed no strong clustering of breeds. Bayesian, maximum parsimony, neighbour joining and UPGMA trees revealed a topology consistent with domestication as well as subsequent introgression of multiple maternal lineages from the wild stocks. Reduced median network analysis provided evidence of population expansion from more than one set of haplotypes. The study also confirmed that Pakistani buffalo are of the river type. The observed mitochondrial D-loop sequence diversity suggests that Pakistani areas bordering India might have contributed to the initiation of domestication of the present-day river buffalo.     

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.     

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.     

Phylogeography of the Coastal Mosquito Aedes togoi across Climatic Zones: Testing an Anthropogenic Dispersal Hypothesis

The coastal mosquito Aedes togoi occurs more or less continuously from subarctic to subtropic zones along the coasts of the Japanese islands and the East Asian mainland. It occurs also in tropical Southeast Asia and the North American Pacific coast, and the populations there are thought to have been introduced from Japan by ship. To test this hypothesis, the genetic divergence among geographic populations of A. togoi was studied using one mitochondrial and three nuclear gene sequences. We detected 71 mitochondrial haplotypes forming four lineages, with high nucleotide diversity around temperate Japan and declining towards peripheral ranges. The major lineage (L1) comprised 57 haplotypes from temperate and subarctic zones in Japan and Southeast Asia including southern China and Taiwan. Two other lineages were found from subtropical islands (L3) and a subarctic area (L4) of Japan. The Canadian population showed one unique haplotype (L2) diverged from the other lineages. In the combined nuclear gene tree, individuals with mitochondrial L4 haplotypes diverged from those with the other mitochondrial haplotypes L1—L3; although individuals with L1—L3 haplotypes showed shallow divergences in the nuclear gene sequences, individuals from Southeast Asia and Canada each formed a monophyletic group. Overall, the genetic composition of the Southeast Asian populations was closely related to that of temperate Japanese populations, suggesting recent gene flow between these regions. The Canadian population might have originated from anthropogenic introduction from somewhere in Asia, but the possibility that it could have spread across the Beringian land bridge cannot be ruled out.     

Y-chromosome and mitochondrial DNA studies on the population structure of the Christmas Island community

Christmas Island is a remote Australian territory located close to the main Indonesian island of Java. Y-chromosome and mitochondrial DNA (mtDNA) markers were used to investigate the genetic structure of the population, which comprises communities of mixed ethnic origin. Analysis of 12 Y-chromosome biallelic polymorphisms revealed a high level of gene diversity and haplotype frequencies that were consistent with source populations in southern China and Southeast Asia. mtDNA hypervariable segment I (HVS-I) sequences displayed high levels of haplotype diversity and nucleotide diversity that were comparable to various Asian populations. Genetic distances revealed extremely low mtDNA differentiation among Christmas Islanders and Asian populations. This was supported by the relatively high proportion of sequence types shared among these populations. The most common mtDNA haplogroups were M* and B, followed by D and F, which are prevalent in East/Southeast Asia. Christmas Islanders of European descent were characterized by the Eurasian haplogroup R*, and a limited degree of admixture was observed. In general, analysis of the genetic data indicated population affinities to southern Chinese (in particular from the Yunnan Province) and Southeast Asia (Thailand, Malaysia, and Cambodia), which was consistent with historical records of settlement. The combined use of these different marker systems provides a useful and appropriate model for the study of contemporary populations derived from different ethnic origins.     

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.     

Mitochondrial DNA variation in Nicobarese Islanders.

The aboriginal populations living in the Nicobar Islands are hypothesized to be descendants of people who were part of early human dispersals into Southeast Asia. However, analyses of ethnographic histories, languages, morphometric data, and protein polymorphisms have not yet resolved which worldwide populations are most closely related to the Nicobarese. Thus, to explore the origins and affinities of the Nicobar Islanders, we analyzed mitochondrial DNA (mtDNA) hypervariable region 1 sequence data from 33 Nicobarese Islanders and compared their mtDNA haplotypes to those of neighboring East Asians, mainland and island Southeast Asians, Indians, Australian aborigines, Pacific Islanders, and Africans. Unique Nicobarese mtDNA haplotypes, including five Nicobarese mtDNA haplotypes linked to the COII/tRNA(Lys) 9-bp deletion, are most closely related to mtDNA haplotypes from mainland Southeast Asian Mon-Kmer-speaking populations (e.g., Cambodians). Thus, the dispersal of southern Chinese into mainland Southeast Asia may have included a westward expansion and colonization of the islands of the Andaman Sea.     

Molecular diversity and phylogeography of the Asian leopard cat, Felis bengalensis, inferred from mitochondrial and Y-chromosomal DNA sequences

To investigate genetic diversity and phylogeography of the Asian leopard cat (Felis bengalensis), mitochondrial DNA (mtDNA) sequences were determined for 39 individuals from various areas. Sequences combining the complete cytochrome b gene (1,140 bp) with the partial control region (646-810 bp) were classified into 24 haplotypes: 21 types from 21 animals, one from eight animals from Tsushima Islands, one from eight animals from Iriomote Island, and one from two animals from Southeast Asia. Phylogenetic trees of the 24 haplotypes clearly showed three clades: a Northern Lineage and Southern Lineages 1 and 2. The Northern Lineage consisted of animals from Tsushima Islands, the Korean Peninsula, the continental Far East, Taiwan, and Iriomote Island. Within the Northern Lineage, genetic contacts could have occurred between geographically neighboring populations before isolation by straits. Southern Lineage 1, comprising Southeast Asian animals, showed higher genetic diversity. Southern Lineage 2 had large genetic distances from other lineages. Within the control region, the Asian leopard cats shared two to four repetitive motifs, and the number of motifs and their constitution were highly variable among individuals. The motifs were polymorphic even within individuals and could be classified into 31 types. Finally, males of mtDNA Southern Lineage 1 had either of two types of the Y-chromosomal gene ZFY, whereas all males of Northern Lineage shared only one type. Our results indicate that the diversity of southern populations is higher and that genetic differentiation among northern local populations reflects past geographical isolation.     

Genetic evidence for the proto-Austronesian homeland in Asia: mtDNA and nuclear DNA variation in Taiwanese aboriginal tribes.

Previous studies of mtDNA variation in indigenous Taiwanese populations have suggested that they held an ancestral position in the spread of mtDNAs throughout Southeast Asia and Oceania (Melton et al. 1995; Sykes et al. 1995), but the question of an absolute proto-Austronesian homeland remains. To search for Asian roots for indigenous Taiwanese populations, 28 mtDNAs representative of variation in four tribal groups (Ami, Atayal, Bunun, and Paiwan) were sequenced and were compared with each other and with mtDNAs from 25 other populations from Asia and Oceania. In addition, eight polymorphic Alu insertion loci were analyzed, to determine if the pattern of mtDNA variation is concordant with nuclear DNA variation. Tribal groups shared considerable mtDNA sequence identity (P>.90), where gene flow is believed to have been low, arguing for a common source or sources for the tribes. mtDNAs with a 9-bp deletion have considerable mainland-Asian diversity and have spread to Southeast Asia and Oceania through a Taiwanese bottleneck. Only four Taiwanese mtDNA haplotypes without the 9-bp deletion were shared with any other populations, but these shared types were widely dispersed geographically throughout mainland Asia. Phylogenetic and principal-component analyses of Alu loci were concordant with conclusions from the mtDNA analyses; overall, the results suggest that the Taiwanese have temporally deep roots, probably in central or south China, and have been isolated from other Asian populations in recent history.