Population structure and phylogeography of Toda buffalo in Nilgiris throw light on possible origin of aboriginal Toda tribe of South India

Coincidence of low genetic diversity and increasing population size in wild gaur populations in the Khao Phaeng Ma Non-Hunting Area, Thailand: A challenge for conservation management under human-wildlife conflict

The gaur (Bos gaurus) is found throughout mainland South and Southeast Asia but is listed as an endangered species in Thailand with a decreasing population size and a reduction in suitable habitat. While gaur have shown a population recovery from 35 to 300 individuals within 30 years in the Khao Phaeng Ma (KPM) Non-Hunting Area, this has caused conflict with villagers along the border of the protected area. At the same time, the ecotourism potential of watching gaurs has boosted the local economy. In this study, 13 mitochondrial displacement-loop sequence samples taken from gaur with GPS collars were analyzed. Three haplotypes identified in the population were defined by only two parsimony informative sites (from 9 mutational steps of nucleotide difference). One haplotype was shared among eleven individuals located in different subpopulations/herds, suggesting very low genetic diversity with few maternal lineages in the founder population. Based on the current small number of sequences, neutrality and demographic expansion test results also showed that the population was likely to contract in the near future. These findings provide insight into the genetic diversity and demography of the wild gaur population in the KPM protected area that can inform long-term sustainable management action plans.

High-Level Gene Flow Restricts Genetic Differentiation in Dairy Cattle Populations in Thailand: Insights from Large-Scale Mt D-Loop Sequencing

Domestication and artificial selection lead to the development of genetically divergent cattle breeds or hybrids that exhibit specific patterns of genetic diversity and population structure. Recently developed mitochondrial markers have allowed investigation of cattle diversity worldwide; however, an extensive study on the population-level genetic diversity and demography of dairy cattle in Thailand is still needed. Mitochondrial D-loop sequences were obtained from 179 individuals (hybrids of Bos taurus and B. indicus) sampled from nine different provinces. Fifty-one haplotypes, of which most were classified in haplogroup "I", were found across all nine populations. All sampled populations showed severely reduced degrees of genetic differentiation, and low nucleotide diversity was observed in populations from central Thailand. Populations that originated from adjacent geographical areas tended to show high gene flow, as revealed by patterns of weak network structuring. Mismatch distribution analysis was suggestive of a stable population, with the recent occurrence of a slight expansion event. The results provide insights into the origins and the genetic relationships among local Thai cattle breeds and will be useful for guiding management of cattle breeding in Thailand.

Population demographic history and population structure for Pakistani Nili-Ravi breeding bulls based on SNP genotyping to identify genomic regions associated with male effects for milk yield and body weight

The domestic Nili-Ravi water buffalo (Bubalus bubalis) is the best dairy animal contributing 68% to total milk production in Pakistan. In this study, we identified genome-wide single nucleotide polymorphisms (SNPs) to estimate various population genetic parameters such as diversity, pairwise population differentiation, linkage disequilibrium (LD) distribution and for genome-wide association study for milk yield and body weight traits in the Nili-Ravi dairy bulls that they may pass on to their daughters who are retained for milking purposes. The genotyping by sequencing approach revealed 13,039 reference genome-anchored SNPs with minor allele frequency of 0.05 among 167 buffalos. Population structure analysis revealed that the bulls were grouped into two clusters (K = 2), which indicates the presence of two different lineages in the Pakistani Nili-Ravi water buffalo population, and we showed the extent of admixture of these two lineages in our bull collection. LD analysis revealed 4169 significant SNP associations, with an average LD decay of 90 kb for these buffalo genome. Genome-wide association study involved a multi-locus mixed linear model for milk yield and body weight to identify genome-wide male effects. Our study further illustrates the utility of the genotyping by sequencing approach for identifying genomic regions to uncover additional demographic complexity and to improve the complex dairy traits of the Pakistani Nili-Ravi water buffalo population that would provide the lot of economic benefits to dairy industry.

Genome-wide analysis reveals genetic diversity, linkage disequilibrium, and selection for milk production traits in Chinese buffalo breeds

The water buffalo is an important dual-purpose livestock that is widespread throughout central and southern China. However, there has been no characterization of the population genetics of Chinese buffalo. Using an Axiom buffalo genotyping array (Thermo Fisher Scientific, Wilmington, DE), we analyzed the genetic diversity, linkage disequilibrium pattern, and signature of selection in 176 Chinese buffaloes from 13 breeds. A total of 35,547 SNP passed quality control and were used for further analyses. Population genetic analysis revealed a clear separation between swamp and river types. Ten Chinese indigenous breeds were clustered into the swamp group, the Murrah and Nili-Ravi breeds were clustered into the river group, and the crossbred breed was closer to the river group. Genetic diversity analysis showed that the swamp group had a lower average expected heterozygosity. Linkage disequilibrium decay distance was much shorter in the swamp group compared with the river group, with an average square of correlation coefficient value of 0.2 of approximately 50 kb. Analysis of runs of homozygosity indicated extensive remote and recent inbreeding within swamp and river groups, respectively. Moreover, one genomic region under selection was detected between the river and swamp groups. Our findings contribute to our understanding of the characterization of population genetics in Chinese buffaloes, which in turn may be used in buffalo breeding programs.

Genetic analysis of river, swamp and hybrid buffaloes of north-east India throw new light on phylogeography of water buffalo (Bubalus bubalis)

This study analysed buffaloes from north-east India and compared their nuclear and mitochondrial DNA variations with buffaloes of mainland India, China, Mediterranean and South-East Asia. Microsatellite genotypes of 338 buffaloes including 210 from six north-east Indian buffalo populations and three mainland Indian breeds were analysed to evaluate their genetic structure and evolutionary relationships. Phylogenetic analysis and multidimensional scaling plot of pairwise FST revealed the clustering of all swamp-type buffaloes of north-east India with Lower Assamese (significantly hybrid type) buffaloes in one plane and all the mainland river buffaloes in another plane while the upper Assamese buffaloes being distinct from both these clusters. Analysis of mtDNA D-loop region of 530-bp length was performed on 345 sequences belonging to 23 buffalo populations from various geographical regions to establish the phylogeography of Indian water buffalo. The swamp buffaloes of north-east India clustered with both the lineages of Chinese swamp buffalo. Multidimensional scaling display of pairwise FST derived from mitochondrial DNA data showed clustering of upper Assamese, Chilika and Mediterranean buffaloes distinctly from all the other Indian buffalo populations. Median-joining network analysis further confirmed the distinctness and ancestral nature of these buffaloes. The study revealed north-east region of India forming part of the wider hybrid zone of water buffalo that may probably extend from north-east India to South-East Asia.

Power of exclusion of 19 microsatellite markers for parentage testing in river buffalo (Bubalus bubalis)

In the present study, 19 microsatellite markers were assessed for their power of exclusion to test parentage in river buffalo. Microsatellite genotypes of 216 unrelated buffaloes belonging to five different breeds were utilized for the study. The probabilities of exclusion were calculated for three hypothetical situations viz. paternity testing (PE1), one parental genotype unavailable (PE2) and exclusion of both parents i.e. substituted offspring (PE3). The mean probability of exclusion across 19 investigated markers in buffalo was 0.578 (PE1), 0.405 (PE2) and 0.764 (PE3) respectively. The probability of exclusion for paternity (PE1) ranged between 0.297 and 0.814 across different markers. The exclusion probability for the cases one parent unavailable (PE2) and substituted offspring (PE3) varied from 0.143 to 0.688 and 0.465 to 0.946 respectively. Polymorphism information content and expected heterozygosity were found to have significantly high correlation with probability of exclusion of microsatellite markers. The cumulative PE1 of nine marker loci was estimated to be 0.9999 while in case of absence of one of the parental genotypes, a minimum of 11 markers were required to achieve a cumulative PE2 of 0.999. In conclusion, the present study proposes two multiplex sets with four and five markers respectively for routine parentage testing in buffalo and an additional set of four markers for doubtful cases of paternity.

Transcript profiling of pattern recognition receptors in a semi domesticated breed of buffalo, Toda, of India

The primary objective of this study was to assess the expression profile and levels of toll-like receptor (TLR) mRNAs in the spleen, lung, mediastinal lymph node (MLN), jejunum, rectum, skin and peripheral blood mononuclear cells (PBMC) of Toda and Murrah buffalos. Spleen and PBMC had increased expression of TLR mRNAs 2, 4, 5, 6, 8, 9 and 10; lung had increased expression of TLR mRNAs 2, 4, 5, 6 and 8, MLN TLR mRNA 6, 9, 10 and decrease in TLR 3 and 7 mRNAs in skin. No significant differences were observed in the expression levels of any of the TLR mRNA in jejunum and rectum. Toda buffaloes showed significantly higher expression levels of TLR 9 mRNA in MLN, TLR mRNAs 1, 5, 6, 9 and 10 in skin and TLR mRNAs 2, 4, 7 and 9 in PBMC than Murrah buffaloes living in the vicinity. Toda and Murrah buffaloes were inoculated with TLR5 (flagellin) and TLR9 (CpG ODN) ligands in vivo and expression levels of the respective TLRs analyzed 12h later. Following CpG inoculation, Toda buffaloes had significantly higher levels of TLR 9 mRNA expression but not in Murrah. However, flagellin induction did not increase TLR 5 mRNA expression in both these breeds. Histological sections of the skin were made and infiltrating cell clusters were graded and quantified. Following CpG inoculation, Toda buffaloes showed higher numbers of infiltrating grade 1 and grade 3 cell clusters while Murrah showed lower numbers of infiltrating grade 1 cells as compared to mock-inoculated skin sections. Flagellin treatment revealed no significant differences in infiltrating cell clusters in both the breeds. The results have shown differential expression of TLR mRNAs in various tissues between two divergent buffalo breeds with the highest difference in TLR expression profile seen in the skin, the largest portal of entry of pathogens, of Toda.