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.     

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.     

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.     

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.     

Diversity analysis at MHC class II DQA locus in buffalo (Bubalus bubalis) indicates extensive duplication and trans-species evolution

Variation at MHC Class II-DQA locus in riverine and swamp buffaloes (Bubu) has been explored in this study. Through sequencing of buffalo DQA, 48 nucleotide variants identified from 17 individuals, reporting 42 novel alleles, including one pseudogene. Individual animal displayed two to seven variants, suggesting the presence of more than two Bubu-DQA loci, as an evidence of extensive duplication. dN values were found to be higher than dS values at peptide binding sites, separately for riverine and swamp buffaloes, indicating locus being under positive selection. Evolutionary analysis revealed numerous trans-species polymorphism with alleles from water buffalo assigned to at least three different loci (Bubu-DQA1, DQA2, DQA3). Alleles of both the sub-species intermixed within the cluster, showing convergent evolution of MHC alleles in bovines. The results thus suggest that both riverine and swamp buffaloes share con-current arrangement of DQA region, comparable to cattle in terms of copy number and population polymorphism.     

Identification of Cu/Zn superoxide dismutase in cattle and river buffaloes

All air-living organisms produce superoxide dismutase (SOD) and several antioxidant enzymes that limit oxidative stress by detoxifying reactive oxygen species. SOD1 gene has been investigated in Egyptian native cattle and buffalos at the level of genomic DNA and cDNA that were extracted from leucocytes. An unexpected band at approximately 370 bp was obtained in cattle genomic DNA and cDNA as well as in buffalo cDNA. SOD1 amplified sequence of native cattle genomic DNA and cDNA showed a 93% alignment. Native cattle genomic DNA SOD1 amplicon shares sequence homology with mRNAs of Bos taurus “similar to superoxide dismutase” (SOD1) sequence of the GeneBank database. It also shares sequence homology with “similar to superoxide dismutase” on B. taurus chromosome BTA13. The results indicate that the genomic DNA of Egyptian native cattle contains SOD1 processed pseudo gene. SOD1 primers amplified three fragments in buffalo genomic DNA which indicates that buffalo genome has different copies of SOD1 due to alternative splicing. It failed to produce the 370 bp fragments found in cattle DNA. The protein analysis revealed no differences between Egyptian native cattle and B. taurus SOD1 mRNA. However, one amino acid, aspartic acid (Asp), in Egyptian native cattle and B. taurus SOD1, is substituted with asparagine (Asn) (D26N) in buffaloes. This amino acid substitution may be due to non-synonymous single nucleotide polymorphisms (nsSNPs). The nsSNPs detected in buffaloes may affect the function of the encoded protein. This study is the first investigation reporting that the resistance of the buffalo to diseases and parasites that afflict cattle may not be acquired but may have a genetic basis.     

Evaluation of Genetic Diversity,Population Structure,and Relationship Between Legendary Vechur Cattle and Crossbred Cattle of Kerala State,India

The legendary Vechur cattle of Kerala, described as a very short breed, and the crossbred (CB) Sunandini cattle population exhibited great phenotypic variation; hence, the present study attempted to analyze the genetic diversity existing between them. A set of 14 polymorphic microsatellites were chosen from FAO-ISAG panel and amplified from genomic DNA isolated from blood samples of 30 Vechur and 64 unrelated crossbred cattle, using fluorescent labeled primers. Both populations revealed high genetic diversity as evidenced from high observed number of alleles, Polymorphic Information Content and expected heterozygosity. Observed heterozygosity was lesser (0.699) than expected (0.752) in Vechur population which was further supported by positive F_IS value of 0.1149, indicating slight level of inbreeding in Vechur population. Overall, F_ST value was 0.065, which means genetic differentiation between crossbred and Vechur population was 6.5%, indicating that the crossbred cattle must have differentiated into a definite population that is different from the indigenous Vechur cows. Structure analysis indicated that the two populations showed distinct differences, with two underlying clusters. The present study supports the separation between Taurine and Zebu cattle and throws light onto the genetic diversity and relationship between native Vechur and crossbred cattle populations in Kerala state.     

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.     

Genetic admixture and population structure analysis of Indian water buffaloes (Bubalus bubalis) using STR markers

Molecular Biology Reports - India has a vast riverine and swamp buffalo diversity adapted to various agro-ecological conditions. In the present study, genetic diversity data for 10 different...     

Molecular prevalence of different genotypes of Theileria orientalis detected from cattle and water buffaloes in Thailand

Here we report on an epidemiological study regarding the molecular prevalence of different genotypes of Theileria orientalis present among domestic cattle and water buffalo populations bred in Thailand. A phylogenetic analysis based on the parasitic gene encoding a major piroplasm surface protein revealed the presence of 5 genotypes (Types 1, 3, 5, 7, and N-3) in cattle and 7 genotypes (Types 1, 3, 4, 5, 7, N-2, and N-3) in water buffaloes. Types 4, 7, and N-3 of T. orientalis were reported for the first time in water buffaloes. The previously reported C and Thai types from Thailand clustered as types 7 and 6, respectively, in the present analysis. Great similarities were observed among nucleotide sequences of isolates of the same genotype from cattle and water buffaloes, and, therefore, water buffaloes were considered to serve as a reservoir for these genotypes of T. orientalis in Thailand. In conclusion, T. orientalis parasites circulating in Thailand are more diverse in their genetic characters than previously anticipated.     

沼泽型水牛Y染色体微卫星标记的筛选与多态性检测

为了研究水牛Y染色体的遗传多样性, 文章以滇东南水牛3个地方群体- 红河(HH)、西双版纳(BN)和普洱(PR)共31头公牛为研究对象, 选取14个家牛Y染色体特异性微卫星标记, 以检测这些标记在水牛Y染色体遗传多样性研究中的可行性。结果表明, 3个标记(INRA008、UMN0103和UMN0504)只有1个等位基因, 表现为单态; 3个标记(UMN1113、UMN0304和BC1.2)均为3个等位基因, 但呈单态; 3个标记(UMN0920、UMN0307和UMN3008)呈现无规律的梯状条带, 所以这9个标记都不适用于水牛的Y染色体遗传多样性研究; 只有5个标记(INRA124、INRA189、BM861、PBR1F1和UMN2001)具有多态性, 表明适用于水牛的Y染色体遗传多样性研究。这5个多态性Y染色体特异微卫星标记在滇东南水牛群体中的平均等位基因数(NA)为2.8000, 平均期望杂合度(He)为0.3998, 基因多样性(GD)为0.4144, 多态信息含量(PIC)为0.3245, Shannon信息熵(SI)为0.5849, 表明滇东南水牛群体的Y染色体具有中等遗传多态性。     

Limited genetic differentiation in Labeo rohita (Hamilton 1822) populations as revealed by microsatellite markers

Labeo rohita, popularly known as rohu is a widely cultured species in the whole Indian subcontinent. Knowledge of the genetic diversity of this species is important to support management and conservation programs which will subsequently help in sustainable production of this species. DNA markers, mostly microsatellite markers are excellent tool to evaluate genetic variation of populations. Genetic variation of three wild and one farm population was assessed using eleven microsatellite loci. In analyzing 192 samples, the number of alleles ranged from 4 to 23; observed heterozygosity 0.500 to 0.870 and expected heterozygosity from 0.389 to 0.878. Exact test for Hardy Weinberg disequilibrium revealed that each riverine sample had at least one locus not in equilibrium except one river. Negative inbreeding coefficients (F_IS) were observed across populations indicating very high level of genetic diversity but little genetic differentiation among populations.     

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

Molecular Characterization and Comparison of Phospholipase C zeta (PLCZ1) Gene Between Swamp (Bubalus carabanensis) and Riverine (Bubalus bubalis) Buffaloes: Its Implications and Future Perspectives

Phospholipase C zeta, a novel sperm-specific protein which is widely known to induce oocyte activation following fertilization, had already been characterized in various mammalian species, but not in water buffaloes thus far. The present study was conducted to initially characterize and compare the sequences of PLCZ1 gene of swamp and riverine buffaloes. Semen samples were collected; total RNA was extracted and reverse-transcribed. PLCZ1 cDNA was then amplified, and submitted for sequencing. Buffalo PLCZ1 gene yielded a sequence of 1905 base pair nucleotides translated into 634?bp amino acids. In general, the buffalo PLCZ1 gene was found to have high sequence identity with cattle and other domestic species. Similarly, significant residues and motifs in PLCZ1 gene sequence are found conserved in water buffaloes. However, there are variations in sequences identified between types of water buffaloes that may play a role in species-specific differences in terms of gene and protein expression, physiological mechanisms, and biological functions. The molecular information on buffalo PLCZ1 gene is highly valuable in subsequent works such as correlation studies on the identified gene variations with semen quality and fertility, and the development of biomarkers for bull fertility.     

Short tandem repeat based analysis of genetic variability in Kanarese buffalo of South India

The goal of the present study was assessing genetic diversity within Kanarese buffalo, the dual purpose breed of South India. A total of 48 unrelated animals were genotyped at 23 short tandem repeat (STR markers)loci. The total number of observed alleles was 180 with a mean of 7.83 per locus, which varied from 3 to 12 across different loci. The mean observed and expected heterozygosity in South Kanara buffaloes was estimated to be 0.518 and 0.712 respectively. Within population inbreeding estimate (F _IS) was significantly positive in most of the investigated loci which resulted in significant deviation from Hardy-Weinberg equilibrium at 19 of 23 loci analyzed. Evaluation of South Kanara buffalo population for mutation drift equilibrium revealed no significant heterozygosity excess under three different models of evolution viz. infinite alleles model (IAM), step-wise mutation model (SMM) and two phase model (TPM), thus indicating the absence of any recent genetic bottleneck. The results of the present study will help in formulating rational breeding strategies as well as conservation of this important germplasm.     

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.     

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.     

Evaluating the potential for reintroducing the endangered wild water buffalo (Bubalus arnee) in Kanha National Park,central India

Reintroducing megafauna to their historic range is an effective strategy to halt their extinctions and restore ecosystems. Wild water buffalo (Bubalus arnee) is an endangered megaherbivore that is lost from 95% of its range. About 90% of its global population (less than 4000) resides within India, in two isolated populations: northeast and central India. The central Indian population is on the verge of extinction and warrants urgent conservation interventions. We assess the potential and provide a strategy for reintroducing buffaloes in Kanha National Park, India. Habitat suitability using the global occurrence of buffalo revealed low-lying grasslands with least human pressure found in Kanha (390 km²) to be suitable. Within this suitable range, we evaluated vegetation composition, forage biomass, and potential carrying capacity. Multidimensional ordination classified these suitable sites into moist and dry grassland clusters. Moist grasslands were found to better suit reintroduction due to lower grazing pressure, higher productivity, and availability of perennial water sources. Distinct matriarchal clades of swamp and river buffalo were observed. Within the river buffalo clade, Indian wild buffaloes formed a distinct cluster with close proximity between northeast and central Indian samples, suggesting northeast buffaloes could be sourced for a founding population in Kanha. Following IUCN guidelines, we discuss the reintroduction strategy that could sustain approximately 200 buffaloes in the area and subsequently expand to other suitable habitats in central India. If implemented as proposed, recovery of this lost ecosystem engineer will help to restore grasslands, and swamp habitats, as well as contribute to its global conservation efforts.     

Oestrus synchronization in buffaloes (Bubalus bubalis)

A study was conducted on 34 Surti buffalo cows to determine the feasibility of synchronizing oestrus using prostaglandin F_2α and a 12-day progesterone intravaginal device. Eighteen cycling buffalo cows having palpable corpora lutea were treated with a single intramuscular injection of 30 mg of prostaglandin F_2α. Three cows exhibited oestrus approximately 54 h after treatment and two of these were diagnosed pregnant 90 days after natural breeding. Sixteen randomly selected post partum cows were treated for 12 days with a progesterone intravaginal device. Ten mg of oestradiol 17β in 5 ml of ether was also injected at the time of insertion of the device. Thirteen cows retained the device for 12 days and 10 of them returned to oestrus 4–5 days after its removal. Eight animals were diagnosed pregnant 90 days after natural breeding. The results indicate that short term progesterone intravaginal device treatment is more reliable than prostaglandin for synchronizing oestrus in buffaloes.