MARKER ASSISTED EVALUATION OF MORPHOLOGICAL AND GENETIC ATTRIBUTES OF SUB-POPULATIONS OF NILI-RAVI BUFFALO: A VULNERABLE DAIRY TYPE RIVERINE BREED OF INDIA

In the present study, we report the distribution of true to type and atypical Nili-Ravi buffalo, a vulnerable dairy type riverine breed of North India and its underlying genetic structure. Out of total investigated buffaloes 73.5% had bilateral wall eyes while 5.4% had unilateral wall eyes and 21.1% had no wall eyes. 41.15% of Nili-Ravi buffaloes maintained in the breeding farm were having typical true to the type characteristics (both eyes walled, white markings in forehead, muzzle/chin, all the four legs and tail) while only 28.5% of Nili-Ravi buffaloes were true to the type under field conditions. Genotypic data were generated in four groups of Nili-Ravi buffalo (FMTNR--Typical Nili-Ravi from farm; FMANR--Atypical Nili-Ravi from farm; FDTNR--Typical Nili-Ravi from field; FDANR--Atypical Nili-Ravi from field) at 16 microsatellite loci. Comparative genetic analysis of various groups of Nili-Ravi buffaloes with Murrah revealed significant between group differences with an estimated global F(ST) of 0.063. Pair-wise F(ST) values ranged from 0.003 (between FDTNR and FDANR) to 0.112 (between FMTNR and FDTNR). Phylogenetic analysis and multi-dimensional scaling revealed clustering of FDTNR and FDANR together while FMTNR and FMANR clustered separately with Murrah in between farm and field Nili-Ravi buffaloes. Based on the results, the paper also proposes three pronged strategy for conservation and sustainable genetic improvement of Nili-Ravi buffalo in India.

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.

Expression profiling of glucose transporter 1 (GLUT1) and apoptotic genes (BAX and BCL2) in milk enriched mammary epithelial cells (MEC) in riverine buffalo during lactation

Lactation is an important physiological process in dairy animals. During lactation, up to 85% of the body glucose is directed toward the mammary glands for milk synthesis. Studies related to lactation physiology are generally carried out on mammary biopsies, which may adversely affect animal health. In the present study, milk enriched MEC were used to study the expression pattern of GLUT1 and apoptotic genes (BAX and BCL2) across different stages of lactation in riverine buffalo in relation to milk yield. MEC were enriched from milk using cytokeratin-8 antibodies coated magnetic beads. Total RNA isolated from enriched MEC showed significant correlation (r(2) = 0.92 ± 0.02) with the milk yield at different stages of lactation. GLUT1 expression pattern correlated with the milk yield as highest GLUT1 expression (4.68 ± 0.79) was observed during peak-lactation (90 days post-parturition), whereas low GLUT1 expression (1.01 ± 0.1, 15 d; 0.71 ± 0.03, 30 d) was observed during early lactation. The BAX/BCL2 ratio was high (1.02 ± 0.2, 15 d; 0.94 ± 0.06, 30 d) during the early phase of lactation, indicating high rate of apoptosis, whereas low BAX/BCL2 ratio (0.25 ± 0.03, 60 d; 0.42 ± 0.04, 90 d) was observed during mid-lactation coinciding with the increase in RNA concentration and milk yield. Highest BAX/BCL2 ratio (1.41 ± 0.3, 120 d; 4.02 ± 0.6, 240 d) was observed during late lactation i.e., 240 days, which was also reflected as decline in milk yield and RNA concentration. Also, BAX/BCL2 ratio in milk enriched MEC was in accordance with RNA concentration in MEC and milk yield at different phases of lactation. Our study showed that expression pattern of genes under study (GLUT1, BAX, and BCL2) in milk enriched MEC correlated well with important physiological properties such as milk yield in buffalo.

Analysis of genetic variations across regulatory and coding regions of kappa-casein gene of Indian native cattle (Bos indicus) and buffalo (Bubalus bubalis)

The promoter region of kappa-casein (κ-CN) gene in Indian native cattle and buffalo breeds was sequenced and analyzed for nucleotide variations. Sequence comparison across breeds of Indian cattle revealed a total of 7 variations in the promoter region, of which − 515 G/T, − 427 C/T, − 385 C/T, − 283 A/G and − 251 C/T were located within consensus binding sites for octamer-binding protein (OCT1)/pregnancy specific mammary nuclear factor (PMF), activator protein-2 (AP2), hepatocyte nuclear factor (HNF-1) and GAL4 transcription factors (TFs), respectively. These variations might be involved in gain or loss of potential transcription factor binding sites (TFBSs). Unlike the other 4 variants, the − 283 (A/G) variant located within HNF-1 TFBS was specific to Indian cattle as this change has not been observed in the Bos taurus sequence. Other TFBSs viz., MGF, TBP, NF-1, milk box and C/EBP were conserved across species. For the Indian native buffalo breeds, only 3 changes were identified in the promoter region; − 305 (A/C), − 160 (T/C) and − 141 (A/G) and most of the TFBSs were found to be conserved. However, deletion of two adjacent nucleotides located in and around binding site for C/EBP TF was identified in buffalo when compared with promoter sequence of bovine κ-CN. For κ-CN of Indian native cattle, a strong linkage disequilibrium (LD) was observed for variations 515 G/T, − 427 C/T and − 385 C/T in the promoter region; and for variations at codons 136 and 148 of exon-IV. Further, among intragenic haplotypes, variation − 427 C/T was found to be in LD with variations at codons 136 and 148. The information generated in the present work provides comprehensive characterization of κ-CN gene promoter and coding regions in Indian cattle and buffaloes and reported variations could become important candidates for carrying out further research in dairy traits.

•

The κ-casein gene (κ-CN) and its promoter region was sequence characterized in 15 Indian native cattle and 8 buffalo breeds.

•

Among the identified variations, four were located within TFBSs in Indian cattle, while all TFBSs were conserved in buffaloes.

•

The κ-CN CDS was highly conserved with only 3 and 2 non-synonymous changes for Indian cattle and buffaloes, respectively.

•

Across promoter and CDS, − 427 C/T was found to be in LD with variation at codons 136 and 148.

Peripheral blood mononuclear cells: a potential cellular system to understand differential heat shock response across native cattle (Bos indicus), exotic cattle (Bos taurus), and riverine buffaloes (Bubalus bubalis) of India

Circulating leukocytes can be used as an effective model to understand the heat stress response of different cattle types and buffaloes. This investigation aimed to determine the temporal profile of HSPs (HSP40, HSP60, HSP70, and HSP90) expression in circulating peripheral blood mononuclear cells (PBMCs) of Murrah buffaloes, Holstein-Friesian (HF), and Sahiwal cows in response to sublethal heat shock at 42 °C. The viability data indicated HF PBMCs to be the most affected to the heat shock, whereas Sahiwal PBMCs were least affected, indicating its better survivability during the heat stress condition. The qRT-PCR expression data showed significant increase in mRNA expression of the analyzed HSPs genes after heat stimuli to the PBMCs under in vitro condition. In each case, the HSPs were most upregulated at 2 h after the heat stress. Among the HSPs, HSP70 was relatively more expressed followed by HSP60 indicating the action of molecular chaperones to stabilize the native conformation of proteins. However, PBMCs from different cattle types and buffaloes showed difference in the extent of transcriptional response. The level of expression of HSPs throughout the time period of heat stress was highest in buffaloes, followed by HF and Sahiwal cows. The higher abundance of HSP70 mRNA at each time point after heat stress showed prolonged effect of heat stress in HF PBMCs. The data presented here provided initial evidence of transcriptional differences in PBMCs of different cattle types and buffaloes and warrant further research.

Classical swine fever in the pygmy hog

The pygmy hog is a rare, small and highly endangered mammal belonging to the Suidae family, and it is presently found only in the Assam state of India. While investigating the cause of death of pygmy hogs housed at a conservation centre for captive breeding and research at Basistha, Assam, it was confirmed that they were susceptible to and died as a result of contracting classical swine fever (CSF), caused by CSF virus (CSFV), which is a highly infectious endemic disease of domestic pigs in India. The post-mortem findings and serum CSFV-specific antibody titres, along with the isolation of CSFV from two pygmy hogs, and further confirmation by CSFV genomic E2 and 5' untranslated region (UTR) gene amplification in PCR (polymerase chain reaction), clearly established the cause of death of the pygmy hogs. Further, on phylogenetic analysis, the pygmy hog CSFV 5' UTR sequences were grouped in the genotype 1.1 cluster of Indian CSFVs, and hence the strains causing infection were closely related to CSFV isolates circulating in domestic pigs. Therefore, the occurrence of CSF in this endangered species may pose a potent threat to their existence unless properly controlled, and thus it needs urgent attention. To the authors' knowledge this is the first report on CSF in pygmy hogs.

Genetic polymorphisms in the bovine toll-like receptor 4 (TLR4) and monocyte chemo attractant protein-1(CCL2) genes: SNPs distribution analysis in Bos indicus Sahiwal cattle breed

Toll-like receptor 4 gene (TLR4) that recognizes the Gram negative bacterial ligand LPS was sequenced in the Bos indicus Sahiwal cattle breed. Ninety four single nucleotide polymorphisms (SNPs) were detected within 10.8 kb gene region. Seventeen of the SNPs were in the coding regions and the one at position 9589(A > G) in exon3 resulted in an amino acid change from Valine to Isoleucine. These SNPs led to generation of 27 TLR4 gene haplotypes. All the Sahiwal animals studied presently showed the occurrence of the genotype CC at gene position 9662, which codes for the amino acid threonine at position 674 of the TLR4 protein, and which had been reported to be associated with lower somatic cell score and, therefore, a lower susceptibility to mastitis, in Taurus cattle. This nucleotide configuration of the Toll-like receptor 4 gene of the Bos indicus Sahiwal cattle breed could possibly indicate toward a lower susceptibility to mastitis in the Sahiwal animals. Monocyte chemo-attractant protein-1 (CCL2) gene encoding for small inducible cytokine A2 that belongs to the CC chemokine family was also sequence characterized in these Sahiwal animals. The CCL2 gene was observed to have 12 polymorphic sites in 3.3 kb region of which one SNP at position 2500 (A > G) in exon 3 resulted in amino acid change from Valine to Isoleucine at position 46 of the mature CCL2 peptide. Seventeen haplotypes of the CCL2 gene were predicted corresponding to 12 genotypes detected.

Genetic diversity within 5'upstream region of Toll-like receptor 8 gene reveals differentiation of riverine and swamp buffaloes

In this study the nucleotide diversity in the 5'untranslated region (UTR) of TLR8 gene in riverine as well as swamp buffaloes has been described. Analysis of the 5'UTR of TLR8 gene showed presence of two SNPs in this region, g.-139G>T and g.-128A>G. A PCR-RFLP assay designed for genotyping of g.-139G>T SNP across 667 samples from 2 buffalo populations revealed a striking difference in allele frequency distribution across the swamp and riverine buffaloes. The frequency of T allele was higher in swamp buffalo as compared to riverine buffalo, ranging from 0.71 to 1. The G allele on the other hand exhibited a higher frequency across all the Indian riverine buffalo breeds/populations. The principal component analysis revealed separate clusters for the riverine and swamp buffaloes, as expected; however, the riverine type Assamese buffalo population of eastern India formed a distinct cluster. Since most of the buffalo populations in the eastern region are swamp type, this demarcation may be related to the difference in immune response in riverine and swamp buffaloes. These preliminary results indicate that the genetic variation observed in 5'upstream region of TLR8 gene, which differentiates swamp and riverine buffalo, needs to be further explored for association with disease susceptibility in buffalo, an important dairy and meat animal of Southeast Asia.

Caprine Toll-like receptor 8 gene sequence characterization reveals close relationships among ruminant species

TLR8 mediates antiviral immunity by recognizing ssRNA viruses and triggers potent antiviral and antitumor immune responses. In this study, approximately 3.5 Kb nucleotide sequence data of caprine TLR8 gene were generated from one sample each of twelve different Indian goat breeds belonging to different geographical regions. Cloning and characterization of cDNA synthesized from RNA purified from goat spleen revealed TLR8 ORF to be of 3102 nucleotides long coding for 1033 amino acids similar to other ruminant species, that is sheep, buffalo and cattle. The sequence analysis at nucleotide level revealed goat TLR8 to be closer to buffalo sharing 99.6% homology, followed by cattle and sheep. Simple Modular Architecture Research Tool (SMART) used for the structural analysis of goat TLR8 showed the presence of 16 leucine-rich repeats (LRRs) along with single Toll/interleukin-1 receptor (TIR) domain. TIR domain when compared with other livestock species was found to be conserved in ruminant species goat, sheep, cattle and buffalo. The phylogenetic analysis also revealed grouping of all ruminant species together, goat being closer to buffalo followed by cattle and sheep. Total 4 polymorphic sites were observed in TLR8 gene of one specimen goat representing each of 12 different breeds studied, all of which were synonymous and present within the coding region. Of these 4 SNPs, two were in ectodomains, one in TIR domain and one was found to be present in transmembrane domain. PCR-RFLP genotyping of two of the SNPs indicated variations in allele frequencies among different goat breeds. The expression profiling in 13 tissues of goat showed maximum expression of TLR8 gene in kidney followed by spleen, lung and lymph node. Overall, our results indicate conservation of TLR8 gene among the ruminant species and low variation within Indian goat breeds.

Novel polymorphisms in UTR and coding region of inducible heat shock protein 70.1 gene in tropically adapted Indian zebu cattle (Bos indicus) and riverine buffalo (Bubalus bubalis)

Due to evolutionary divergence, cattle (taurine, and indicine) and buffalo are speculated to have different responses to heat stress condition. Variation in candidate genes associated with a heat-shock response may provide an insight into the dissimilarity and suggest targets for intervention. The present work was undertaken to characterize one of the inducible heat shock protein genes promoter and coding regions in diverse breeds of Indian zebu cattle and buffaloes. The genomic DNA from a panel of 117 unrelated animals representing 14 diversified native cattle breeds and 6 buffalo breeds were utilized to determine the complete sequence and gene diversity of HSP70.1 gene. The coding region of HSP70.1 gene in Indian zebu cattle, Bos taurus and buffalo was similar in length (1,926 bp) encoding a HSP70 protein of 641 amino acids with a calculated molecular weight (Mw) of 70.26 kDa. However buffalo had a longer 5' and 3' untranslated region (UTR) of 204 and 293 nucleotides respectively, in comparison to Indian zebu cattle and Bos taurus wherein length of 5' and 3'-UTR was 172 and 286 nucleotides, respectively. The increased length of buffalo HSP70.1 gene compared to indicine and taurine gene was due to two insertions each in 5' and 3'-UTR. Comparative sequence analysis of cattle (taurine and indicine) and buffalo HSP70.1 gene revealed a total of 54 gene variations (50 SNPs and 4 INDELs) among the three species in the HSP70.1 gene. The minor allele frequencies of these nucleotide variations varied from 0.03 to 0.5 with an average of 0.26. Among the 14 B. indicus cattle breeds studied, a total of 19 polymorphic sites were identified: 4 in the 5'-UTR and 15 in the coding region (of these 2 were non-synonymous). Analysis among buffalo breeds revealed 15 SNPs throughout the gene: 6 at the 5' flanking region and 9 in the coding region. In bubaline 5'-UTR, 2 additional putative transcription factor binding sites (Elk-1 and C-Re1) were identified, other than three common sites (CP2, HSE and Pax-4) observed across all the analyzed animals. No polymorphism was found within the 3'-UTR of Indian cattle or buffalo as it was found to be monomorphic. The promoter sequences generated in 117 individuals showed a rich array of sequence elements known to be involved in transcription regulation. A total of 11 nucleotide changes were observed in the promoter sequence across the analyzed species, 3 of these changes were located within the potential transcription factor binding domains. We also identified 4 microsatellite markers within the buffalo HSP70.1 gene and 3 microsatellites within bovine HSP70.1. The present study identified several distinct changes across indicine, taurine and bubaline HSP70.1 genes that could further be evaluated as molecular markers for thermotolerance.

Detection of polymorphism and sequence characterization of Toll-like receptor 7 gene of Indian goat revealing close relationship between ruminant species

In this study, approximately 3.4 kb nucleotide sequence of caprine TLR7 (Toll-like receptor 7) gene was generated from twelve different Indian goat breeds belonging to different geographical regions. Goat TLR7 gene ORF (Open Reading Frame) was found to be 3141 nucleotides long coding for 1046 amino acids similar to sheep. The sequence analysis at nucleotide level revealed goat TLR7 having 99.5% homology with sheep, followed by other livestock species. Simple Modular Architecture Research Tool (SMART) was used for the structural analysis of goat TLR7 that showed the presence of 22 leucine rich repeats (LRRs) along with single Toll/interleukin-1 receptor (TIR) domains. TIR domain, when compared, was found to be similar in ruminant species, goat, sheep, cattle, and buffalo. The phylogenetic analysis also revealed grouping of all ruminant species together, goat being closer to sheep followed by cattle and buffalo. A total of 22 polymorphic sites were observed in TLR7 gene of 24 goats representing 12 different breeds, out of which 19 were present within the coding region and three in 3'UTR. Out of the seven nonsynonymous SNPs, two were in ectodomains and one in TIR domain. Overall our results indicate substantial variation within goat TLR7 gene, which could be exploited for association with disease susceptibility.

Sequence characterization of river buffalo Toll-like receptor genes 1-10 reveals distinct relationship with cattle and sheep

The present study was undertaken to characterize the full-length transcripts of Toll-like receptor (TLR) genes 1-10 of river buffalo. The conceptualized amino acid identity of bubaline TLRs ranged between 86% to 100% with ruminants, while it ranged between 45% to 91% with other vertebrate species. Simple modular architecture tool (SMART) analysis revealed the presence of TIR domains and varying numbers of leucine-rich repeat motifs in all the buffalo TLRs. With respect to TIR domains, TLRs 1, 2 and 3 of river buffalo were found to have 99.3% identity with cattle and 100% identity of TLRs 4, 6 and 10 with sheep. Phylogenetic analysis of TLRs of buffalo and different vertebrate species revealed the clustering of major TLR gene subfamilies with high bootstrap values. The evolutionary relationship between buffalo and other ruminant species was found to vary among different TLRs. In order to understand the relationship between TLRs of different ruminant species, multidimensional scaling (MDS) analysis of pairwise amino acid differences between different species within each TLR was performed. Buffalo and cattle were found to be closely related only with respect to TLRs 1, 2 and 7, while buffalo and sheep were found to be clustering together with respect to TLRs 3, 6, 8 and 10. The distinct relationship of bubaline TLRs with cattle and sheep revealed the possible differences in the pathogen recognition receptor systems in these animals and consequently the differences in their susceptibility/resistance to various invading organisms.

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.

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

We report the genetic structure and evolutionary relationship of the endangered Toda buffalo of Nilgiris in South India with Kanarese and two other riverine buffalo breeds. The upgma phylogeny drawn using Nei's distance grouped South Kanara and Toda buffaloes at a single node while Marathwada and Murrah together formed a separate node. Principal component analysis was performed with pairwise interindividual chord distances which revealed clustering of Murrah and Marathwada buffaloes distinctly, while individuals of Toda and South Kanara breeds completely intermingled with each other. Furthermore, there were highly significant group variances (p < 0.01) when the breeds were grouped based on phylogeny, thus revealing the existence of cryptic genetic structure within these buffalo breeds. To know the evolutionary relationship among these breeds, 537-bp D-loop region of mitochondrial DNA was analysed. The phylogenetic analysis of mtDNA haplotypes following NJ algorithm with Chinese swamp buffalo as outgroup revealed a major cluster that included haplotypes from all the four investigated breeds and two minor clusters formed by South Kanara and Toda haplotypes. Reduced median network analysis revealed haplotypes of South Kanara and Toda to be quite distinct from the commonly found haplotypes indicating that these might have been ancestral to all the present-day haplotypes. Few mutations in two of the haplotypes of South Kanara buffalo were found to have contributed to ancestral haplotypes of Toda buffalo suggesting the possible migration of buffaloes from Kanarese region towards Nilgiris along the Western Ghats. Considering the close social, economic and cultural association of Todas with their buffaloes, the present study supports the theory of migration of Toda tribe from Kanarese/Mysore region along with their buffaloes.

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

5'-UTR-based phylogenetic analysis of Classical swine fever virus isolates from India

Classical swine fever (CSF) caused by Classical swine fever virus (CSFV) is a globally significant disease of pigs. Genetic typing of CSFV isolates can help in understanding the epidemiology of disease and trace down the source of outbreak. 5'-UTR sequence analysis and subsequent genetic classification of nine CSFV field isolates from India indicated that 3 isolates belonged to genotype 2.1 and were closely related to European CSFV strains. The remaining 6 isolates belonged to genotype 1 that contained old and new strains. However, the genotype 2.1 group consisted of recent field isolates only. The study showed circulation of both genotypes 1 and 2.1 in north-eastern part of India.