Microsatellite based Genetic Variation among the Buffalo Breed Populations in Pakistan

Genetic characterization and diversity assessment in 'Bhangor' indigenous swamp buffalo population using heterologous microsatellite markers

'Bhangor' newly identified swamp buffalo population from North East Indian, was characterized using microsatellite markers. Genomic DNA was isolated from blood samples of 76 unrelated animals, 15 microsatellite markers (CSSM33, BM1818, CSRM60, HEL13, ILSTS019, ILSTS025, ILSTS028, ILSTS029, ILSTS033, ILSTS036, ILSTS056, ILSTS058, ILSTS061, ILSTS089 and ETH003) were found to be highly polymorphic in the population of the selected markers. A total of 114 alleles were observed, which ranged from 3 in CSRM60 and ILSTS025 locus to 12 in ILSTS056 and ILSTS061. The mean effective number of alleles across all polymorphic loci was found to be 3.76. The overall mean expected heterozygosity and unbiased expected heterozygosity values were 0.67 and 0.68, ranging from 0.067 (ILSTS025) to 0.85 (ILSTS058) and 0.068 (ILSTS025) to 0.86 (ILSTS058), respectively. Within the population, the inbreeding estimates (FIS) ranged between -0.4352 and 0.804, with an average FIS of 0.114 ± 0.033. The outcome for infinite allele model (IAM), two-phase model (TPM) and test for mode shift revealed the absence of any recent bottleneck in the investigated buffalo population. The population was found to be in optimum diversity based on polymorphic microsatellite markers. With fast changing agro-climatic conditions; there is an urgent need to characterize the nondescript livestock populations.

A Novel Insight into the Identification of Potential SNP Markers for the Genomic Characterization of Buffalo Breeds in Pakistan

Domestic buffaloes (Bubalus bubalis), known as water buffaloes, play a key role as versatile multipurpose agricultural animals in the Asiatic region. Pakistan, with the second-largest buffalo population in the world, holds a rich domestication history of buffaloes. The overall trends in buffalo production demand the genomic characterization of Pakistani buffalo breeds. To this end, the resequencing data of Pakistani breeds, along with buffalo breeds from 13 other countries, were retrieved from our previous study. This dataset, which contained 34,671,886 single-nucleotide polymorphisms (SNPs), was analyzed through a pipeline that was developed to compare possible allele differences among breeds at each SNP position. In contrast, other available tools only check for positional SNP differences for breed-specific markers. In total, 1918, 1549, 404, and 341 breed-specific markers were identified to characterize the Nili, Nili-Ravi, Azakheli, and Kundi breeds of Pakistani buffalo, respectively. Sufficient evidence in the form of phenotypic data, principal component analysis, admixture analysis, and linkage analysis showed that the Nili breed has maintained its distinct breed status despite sharing a close evolutionary relationship with the Nili-Ravi breed of buffalo. In this era of genome science, the conservation of these breeds and the further validation of the given selection markers in larger populations is a pressing need.

Genetic Diversity Relationship in Azakheli Buffalo Inferred from mtDNA and MC1R Sequences Comparison

Azakheli is relatively smaller riverine breed with a very peculiar characteristics kept under unique traditional husbandry practices in comparison with rest of the Pakistani buffalo breeds; however, milk production is comparable. The present study was conducted to evaluate the genetic diversity of Azakheli breed. A total of sixty-six blood samples were collected for the amplification of mtDNA D-loop region and MC1R gene sequencing analysis. Median-joining network analysis of 191 mtDNA D-loop sequences of Azakheli and eight Indian riverine buffalo breeds clustered into three clades. Ancient Azakheli Region 1 clade was the oldest with the highest mutation steps and was present close to the root of UPGMA phylogenetic tree. There was 5 mutated lines distance between Pakistan buffalo and Indian riverine buffaloes. The populations of neighboring countries did not share any haplotypes with Azakheli buffalo of Pakistan. Possibly, residing for so long in the cold atmosphere and high elevation regions caused the mutation in mtDNA D-loop, though these conditions did not affect the overall performance of Azakheli as milch buffalo breed of Pakistan. MC1R analyses showed high mutations in Azakheli of Albino phenotype and all the black phenotype individuals of Azakheli buffalo share haplotypes with dominant Chinese and Indian black phenotypes buffaloes in MC1R median-joining network, indicating the reason of black coat color is due to MC1R gene. The haplotype diversity and nucleotide diversity was (H. 0.923, Pi: 0.00895) in Azakheli. Current results illustrated Asian ancestry for Azakheli buffalo, and mtDNA and MC1R analyses provided further evidence. Additional genetic analyses and archeological studies may provide further insight into the domestication period and history of Azakheli buffalo breed. The further studies are required on different coat colors with different genes on Azakheli buffalo to understand the phenotype variation.

Evaluation of Genetic Diversity and Structure of Turkish Water Buffalo Population by Using 20 Microsatellite Markers

The present study was aimed to investigate the genetic diversity among 17 Turkish water buffalo populations. A total of 837 individuals from 17 provincial populations were genotyped, using 20 microsatellites markers. The microsatellite markers analyzed were highly polymorphic with a mean number of alleles of (7.28) ranging from 6 (ILSTS005) to 17 (ETH003). The mean observed and expected heterozygosity values across all polymorphic loci in all studied buffalo populations were 0.61 and 0.70, respectively. Observed heterozygosity varied from 0.55 (Bursa (BUR)) to 0.70 (Muş (MUS)). It was lower than expected heterozygosity in most of the populations indicating a deviation from Hardy-Weinberg equilibrium. The overall value for the polymorphic information content of noted microsatellite loci was 0.655, indicating their suitability for genetic diversity analysis in buffalo. The mean F_IS value was 0.091 and all loci were observed significantly deviated from Hardy-Weinberg Equilibrium (HWE), most likely based on non-random breeding. The 17 buffalo populations were genetically less diverse as indicated by a small mean F_ST value (0.032 ± 0.018). The analysis of molecular variance (AMOVA) analysis indicated that about 2% of the total genetic diversity was clarified by population distinctions and 88 percent corresponded to differences among individuals. The information produced by this study can be used to establish a base of national conservation and breeding strategy of water buffalo population in Turkey.

Determination and phylogenetic analysis of the complete mitochondrial genome of Bubalus bubalis Linnaeus, 1758 breed Nili-Ravi (Artiodactyla: Bovidae)

Nili-Ravi buffalo (Bubalus bubalis Linnaeus, 1758 breed Nili-Ravi, NRB) is a famous water buffalo breed in the world. It is the first time that the complete mitochondrial genome sequence of the NRB was reported. The total length of the mtDNA is 16,356 bp, It contains the typical structure, including 22 transfer RNA genes, 2 ribosomal RNA genes, 13 protein-coding genes, and 1 non-coding control region (D-loop region). The overall composition of the mtDNA was estimated to be 33.11% for A, 26.45% for T, 26.55% for C, and 13.89% for G. Phylogenetic analyses using neighbor-joining (N-J) computational algorithms showed that the analyzed 18 Ruminantia species are divided into four major clades: Bovidae, Cervidae, Giraffidae, and Atilocapridae. In addition, our work confirmed that NRB has a close genetic relationship with B. bubalis isolate India 4.

Genetic characterization and assessment of diversity in Pandharpuri buffalo breed of India using heterologous microsatellite markers

The present study was aimed at genetic characterization and diversity analysis of Pandharpuri buffalo population using the Food and Agriculture Organization (FAO) recommended bovine microsatellite markers. Genomic DNA was isolated from blood samples of 50 unrelated animals and a total of 23 microsatellite loci were amplified using polymerase chain reaction. Among 23 recommended microsatellite markers, 14 markers (BM2113, BM1818, CSSM66, HEL13, INRA037, ILSTS05, HAUT27, INRA023, INRA035, HEL5, ETH3, NRA063, MM12 and ETH10) were found to be highly polymorphic in Pandharpuri population. The amplified products were run on polyacrylamide gel electrophoresis (PAGE) apparatus along with a ladder; subsequently, the allele typing was done based on silver staining of the gel. The effective and observed number of alleles, heterozygosity (expected and unbiased) estimates and polymorphic information content (PIC) levels were estimated for each locus. A total of 87 alleles were secured for 14 polymorphic loci studied giving an overall average of 6.21 alleles per locus. The number of alleles ranged from 2 (INRA063) to 9 (BM1818 and HEL13). The mean effective number of alleles across all polymorphic loci was found to be 4.28. The overall mean expected heterozygosity and unbiased expected heterozygosity values were 0.77 and 0.76, ranging from 0.50 (INRA063) to 0.88 (BM1818) and 0.50 (INRA063) to 0.88 (BM1818), respectively. The average PIC estimate across all polymorphic loci was 0.73 ranging from 0.373 (INRA063) to 0.864 (BM1818). In the present study, the characterization and diversity estimates are reported for Pandharpuri population. It was found to maintain optimum diversity based on 14 polymorphic microsatellite markers.