Diversity assessment of a lesser known buffalo population from Central India and its comparative evaluation reveals presence of sufficient genetic variation and absence of selection.

Planned breeding and conservation strategies for a lesser-known population require an assessment of complete genetic diversity and population structure analysis in addition to its morphometric characteristics. In the present study, a comparative analysis of the genetic structure of a rare buffalo population, namely Chhattisgarhi, was extensively studied using a panel of FAO-recommended microsatellite markers along with well-established breeds namely Murrah, Nili-Ravi, Gojri, Kalahandi, and Nagpuri. Mode shift analysis indicated the absence of genetic bottleneck in the recent past. Assessment of genetic diversity indices across all loci indicated the presence of sufficient genetic variation within and between populations. Analysis of molecular variance between the six different buffalo populations attributed 19.05% of the variations to between-population differentiation. Cluster analyses using DAPC and Bayesian approach along with the phylogenetic tree based on UPGMA grouped six populations into three groups. The Chhattisgarhi population was revealed to be genetically closer to Nagpuri and Kalahandi populations. The study reveals the presence of sufficient genetic diversity within the Chhattisgarhi population and indicates the absence of a systematic selection program. We suggest improvement and conservation programs should be planned for this breed in the near future through short-term selection.

Discriminant analysis reveals differences in breed structure in river buffaloes of Maharashtra (India).

Buffaloes are a major milch species of India, and owning a true-to-breed buffalo herd is a reverie for livestock keepers, as it substantiates their incomes. Multivariate discriminate analysis is an efficient methodology to classify different farm animal populations based on canonical discriminant function. The present study was conducted on 392 adult buffaloes using 18 body biometric traits. Except for body height and body length, all the body characters revealed significant differences amongst the four buffalo populations. Unlike cattle, distance between horns and rump width were more discriminating factors, compared to body height, leg length, chest girth, and body length. We report the most suitable discriminate function which can discriminate with maximum between-group variation in four studied populations. Mahalanobis distances were calculated between populations and they ranged from 0.786 to 1.227. Purnathadi has revealed a distinct phenotypic structure, and these buffaloes have a distinct native breeding region comprising Akola, Amravati, and Buldhana districts of Maharashtra state. On the other hand, Elachpuri animals failed to depict physical and biometric discriminations from Nagpuri, Marathwadi, and Purnathadi buffalo populations, and Elachpuri can be considered as a strain of Nagpuri buffalo.

Expression profile of different classes of proteases in milk derived somatic cells across different lactation stages of indigenous cows (Bos indicus) and riverine buffaloes (Bubalus bubalis).

Proteases play a significant role in milk and its products by affecting flavor, texture and longevity. The expression of endogenous proteases varies across different stages of lactation. The study was conducted to understand the transcriptional pattern of different classes of protease-pathways associated genes (CTSB, CTSD, CTSH, CTSL, CTSK, CTSS, CTSZ, PLAU, PLAT) and potential protease inhibitors (SERPIN E2 and SERPIN F2) in 40 milk somatic cells (MSC) samples isolated during early, peak, mid and late lactation stages of Sahiwal cows and Murrah buffaloes - the two most important dairy breeds of India. In Sahiwal cows, except CTSK and PLAU, the expression of other proteases class was not affected significantly (p > 0.05) across lactation stages. However, in Murrah buffaloes, the expression of different proteases increased as the lactation progressed. Most of the proteases showed lower expression during early and peak lactation stages while their expression tends to increase during mid to late lactation stages. The overall trend was somewhat similar in both the dairy species albeit the level of expression was higher in buffalo MSC as compared to cow MSC. The study has provided valuable information on expression kinetics of different proteases in milk somatic cells of two major dairy breeds of India.

Identification of novel polymorphism in mammary-derived growth inhibitor gene of water buffalo and its expression analysis in the mammary gland.

Mammary-derived growth inhibitor (MDGI), a member of the lipophilic family of fatty acid-binding proteins, plays an important role in the development, regulation, and differentiation of the mammary gland. The aim of the study was to identify polymorphism in the MDGI gene and its expression analysis in the mammary gland at various stages of lactation, in Indian buffalo. Nucleotide sequence analysis of MDGI gene in different breeds of riverine and swamp buffaloes revealed a total of 16 polymorphic sites and one Indel. Different transcription factor binding sites were predicted for buffalo MDGI gene promoter sequence, using online tools and in-silico analysis indicating that the SNPs in this region can impact the gene expression regulation. Phylogenetic analysis exhibited the MDGI of buffalo being closer to other ruminants like cattle, yak, sheep, and goats. Further, the expression analysis revealed that buffalo MDGI being highly expressed in well-developed mammary glands of lactating buffalo as compared to involution/non-lactating and before functional development to start the milk production stage in heifers. Stage-specific variation in expression levels signifies the important functional role of the MDGI gene in mammary gland development and milk production in buffalo, an important dairy species in Southeast Asia.

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.

Identification of novel allelic patterns and evolutionary lineage of BoLA MHC class II DQA locus in indicine and taurine cattle.

Among different cattle types, Bos indicus are known for their ability to better resist the tropical microbial infections comparatively, wherein MHC molecules play a significant role. In this study allelic diversity at MHC locus, DQA of Bos indicus, Bos taurus and crossbred of taurine-indicus has been explored to understand the possible role of MHC region in differential immune response. Thirty nine different DQA alleles were identified, out of which 14 were novel, along with documentation of duplication of DQA alleles. Indicus cattle population presented diverse types of DQA alleles compared to crossbred and exotic. Translated amino acid sequence analysis indicated, codon 64 and 50 of peptide binding sites being highly polymorphic and most of the indicus cattle presented alanine and arginine amino acid at position 64 and 50. Within breed genetic variation found to be higher than between breeds. Because of their ability to bind and subsequently respond to a wide array of antigens, the newly identified DQA alleles with high diversity present in the form of duplicated haplotypes in different combinations in cattle populations provided significant insights into probable role of this MHC locus in better tropical disease combating ability and genetic fitness of indicus cattle.

In Silico Analysis of HSP70 Gene Family in Bovine Genome.

Heat shock proteins (HSPs), members of molecular chaperones families fulfill essential roles under normal conditions and provide protection and adaptation during and after stress. Among different HSPs, HSP70 kDa family of proteins is most abundant and well-studied in human and mouse but has not yet been characterized in bovines. In silico analysis was performed to characterize members of HSP70 gene family in bovine genome and a total of 17 genes of bovine HSP70 gene family were identified. The members of HSP70 family were distributed over 12 chromosomes with gene size ranging from 1911 (HSPA2) to 54,017 bp (HSPA4). Five genes were intronless, while rest of 12 genes were multiexonic. Phylogenetic analysis of HSP70 gene family distinguished them into eight major evolutionary groups wherein members of group 1 were most divergent and quite dissimilar than from rest of the HSP70 sequences. Domain structure of all bovine HSP70 genes was conserved and three signature patterns HSP70_1, HSP70_2, and HSP70_3 were identified. HSPA8, HSP9, and HSPA1A showed comparatively higher expression in majority of tissues. Like humans, bovine HSP70 family was characterized by remarkable evolutionary diversity. The analysis also suggested resemblance of bovine HSP70 family to that of human compared to mouse. Overall, the study indicates the presence of diversity for structure, function, localization, and expression in the bovine HSP70 family chaperons which could form the basis to understand thermotolerance/adaptive changes in the bovines.

Assessment of genetic diversity and bottleneck in Purnathadi buffaloes using short tandem repeat markers.

This study is the maiden attempt for genetic characterization of the Purnathadi, a phenotypically distinct buffalo population of the western Vidarbha region of Maharashtra and to explore genetic diversity using STR markers. A total of 48 unrelated Purnathadi buffaloes from the entire native tract were genotyped using a battery of 25 heterologous microsatellite markers. 5' end of forward primer of each microsatellite marker was labeled with one of the fluorescent dyes, viz., FAM (Blue), VIC (Green), NED (Yellow) or PET (Red) to assess the fragment length of genotyped PCR product with automated DNA sequencer (ABI 3100). 23 microsatellite loci (except ETH003 and ILSTS030) amplified successfully and adequately high allelic diversity (observed: 0.615 ± 0.043 and expected: 0.655 ± 0.037) was reported with 162 distinct microsatellite alleles. Sufficiently high Shannon index and PIC indicated the suitability of markers to evaluate genetic diversity in Purnathadi buffaloes. Within-population inbreeding estimates (FIS) for Purnathadi buffalo ranged between -0.171 and 0.495 with global FIS average of 5.9%. The outcome for IAM, TPM and test for mode shift revealed the absence of any recent bottleneck in Purnathadi buffalo. Pairwise FST (genetic differentiation) and gene flow between Purnathadi, Nagpuri and Marathwadi buffaloes were estimated using genotype data of 19 microsatellite markers. Lowest FST (0.031) was observed between Nagpuri and Purnathadi buffaloes with highest gene flow of 7.91% and highest FST (0.094) was between Purnathadi and Marathwadi populations. Present findings will definitely support in designing breeding plan for genetic improvement, as well as for developing conservation strategies of Purnathadi buffalo population. The comparative molecular study with other breeds of the regions is needed.

Genomic diversity and selection sweeps identified in Indian swamp buffaloes reveals it's uniqueness with riverine buffaloes.

The present investigation was focused to study genomic diversity of Indian swamp buffalo populations through reduced representation approach (ddRAD). The heterozygosity (FST) among the swamp buffaloes was 0.11 between Assam and Manipuri; 0.20 between swamp (Manipuri) and riverine buffaloes; 0.30 between swamp (Manipuri) and cattle. The average observed and expected heterozygosity in swamp buffalo populations was 0.254 and 0.221 respectively. The Inbreeding coefficient (FIS) value was 0.02 among the swamp buffaloes. PCA and structure analysis revealed Manipuri swamp buffalo was genetically distinct and closely related to Nagaland swamp buffalo and least to Assam swamp buffalo. Identification of selective sweeps revealed 1087 regions to have undergone selection related to immune response, adaptation and nervous system. A total of 3451 SSRs were identified in the genome of swamp buffaloes. The study evidenced the genomic diversity in the swamp buffalo populations and its uniqueness in comparison with riverine buffalo and cattle.

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.

Differential expression of cytokines in PBMC of Bos indicus and Bos taurus × Bos indicus cattle due to Brucella abortus S19 antigen.

Brucellosis is the most dreadful disease among bovines, although breed differences have been observed in prevalence of disease, worldwide. In present study, antibody response and relative expression of proinflammatory cytokines was compared in Bos indicus (zebu) and Bos taurus × Bos indicus (crossbred) cattle vaccinated by live attenuated Brucella abortus S19 antigen. Six female calves (4-6 months age) of both groups were vaccinated with B.abortus S19 strain. Blood samples were collected before vaccination (0d) and 7th (7d), 14th (14d) and 28th (28d) days after vaccination. Indirect ELISA showed high (p < .05) anti-Brucella antibody level after vaccination; with no significant difference between the groups. During Real-time expression, IFNγ, TNFα, IL6 and IL10 genes initially showed down regulation followed by upregulation in both the groups; however, the trend was much prominent in crossbreds. The expressions of IFNγ, TNFα and IL6, proinflammatory molecules important for initial containment of the Brucella were significantly (p < .01) higher in crossbred. The study showed that the Sahiwal cattle were less responsive to B.abortus S19 antigen than crossbreds, indicating its lower sensitivity to the Brucella, comparatively. In contrary, higher expression of the proinflammatory molecules in crossbreds could be important for containment of the organism during initial stage of infection.

Evolutionary dynamics of meiotic recombination hotspots regulator PRDM9 in bovids.

Hybrid sterility or reproductive isolation in mammals has been attributed to allelic incompatibilities in a DNA-binding protein PRDM9. Not only is PRDM9 exceptional in being the only known 'speciation gene' in vertebrates, but it is also considered to be the fastest evolving gene in the genome. The terminal zinc finger (ZF) domain of PRDM9 specifies genome-wide meiotic recombination hotspot locations in mammals. Intriguingly, PRDM9 ZF domain is highly variable between as well as within species, possibly activating different recombination hotspots. The present study characterized the full-length coding sequence of PRDM9 in cattle and buffalo and explored the diversity of the ZF array in 514 samples from different bovids (cattle, yak, mithun, and buffalo). Substantial numerical and sequence variability were observed in the ZFs, with the number of repeats ranging from 6 to 9 in different bovines. Sequence analysis revealed the presence of 37 different ZFs in cattle, 3 in mithun, 4 in yak, and 13 in buffaloes producing 41 unique PRDM9 alleles in these species. The posterior mean of dN/dS or omega values calculated using Codeml tool of PAMLX identified sites -5, -1, +2, +3, +4, +5, and +6 in the ZF domain to be evolving positively in the studied species. Concerted evolution which typifies the evolution of this gene was consistently evident in all bovines. Our results demonstrate the extraordinary diversity of PRDM9 ZF array across bovines, reinforcing similar observations in other metazoans. The high variability is suggestive of unique repertoire of meiotic recombination hotspots in each species.

Expression Analysis of Genes Associated with Prolificacy in FecB Carrier and Noncarrier Indian Sheep.

The effect of FecB mutation on the gene expression in FecB carrier and noncarrier estrous synchronized ewes, has been analyzed. For this study the whole ovarian tissues and Graafian follicles were collected from estrus synchronized FecB carrier Garole, and non-carrier Deccani Indian sheep, showing remarkable differences in the numbers of preovulatory follicles among two groups. Eleven potential candidate genes (BMP15, GDF9, BMP4, BMP7, BMPR1B, BMPR1A, SMAD9, LHCGR, FSHR, IGF1R, and STAT5) were selected for their expression analysis by SybrGreen based real-time PCR, across ovaries and Graafian follicles of different fecundity groups, for having better insights into the effect of FecB genotypes on follicular development. Variable expression was observed for almost all the genes included in the present study among high and low fecundity groups that was most significant for the BMP7, BMP4, LHCGR, and FSHR transcripts in the ovarian follicles of high and low fecundity ewes, indicating their importance in governing the fecundity in FecB carrier, Indian Garole sheep. BMP4 expression among the genes studied was significantly higher in FecB carrier Garole sheep. This study confirms the changes in mRNA expression of the genes implicated in follicular development in FecB carrier and noncarrier Indian sheep breeds.

High genetic diversity and distribution of Bubu-DQA alleles in swamp buffaloes (Bubalus bubalis carabanesis): identification of new Bubu-DQA loci and haplotypes.

In this study, genetic diversity analysis of MHC class II-DQA locus helped in identification of 25 new Bubu-DQA nucleotide sequences in swamp buffaloes (Bubalus bubalis carabanesis, Bubu). Phylogenetic analysis revealed the distribution of the buffalo DQA sequences in two major clusters of DQA1 and DQA2 genes, sharing common lineages with corresponding cattle alleles, possibly due to trans-species evolution. However, a highly divergent sequence, Bubu-DQA*2501, homologous to cattle (BoLA) DQA3 allele, was identified, indicating the existence of an additional locus; putative DQA3 in buffalo. PCR-RFLP analysis revealed extensive duplication of DQA locus in swamp buffaloes, sharing DQA1, DQA2, and DQA3 alleles in different combinations in duplicated haplotypes. Higher dN than dS values and Wu-Kabat variability at peptide-binding regions in Bubu-DQA indicated high polymorphism with balancing selection. Levels of genetic diversity within DQA sequences and duplication in a small population of swamp buffalo indicate the genetic richness of the species, important for fitness.

Molecular Characterization of Buffalo Haptoglobin: Sequence Based Structural Comparison Indicates Convergent Evolution Between Ruminants and Human.

Haptoglobin (Hp) protein has high affinity for hemoglobin (Hb) binding during intravascular hemolysis and scavenges the hemoglobin induced free radicals. Earlier reports indicate about uniqueness of Hp molecule in human and cattle, but in other animals, it is not much studied. In this paper, we characterized buffalo Hp molecule and determined its molecular structure, evolutionary importance, and tissue expression. Comparative analysis and predicted domain structure indicated that the buffalo Hp has an internal duplicated region in α-chain only similar to an alternate Hp2 allele in human. This duplicated part encoded for an extra complement control protein CCP domain. Phylogenetic analysis revealed that buffalo and other ruminants were found to group together separated from all other non-ruminants, including human. The key amino acid residues involved in Hp and Hb as well as Hp and macrophage scavenger receptor, CD163 interactions in buffalo, depicted a significant variation in comparison to other non-ruminant species. Constitutive expression of Hp was also confirmed across all the vital tissues of buffalo, for the first time. Results revealed that buffalo Hp is both structurally and functionally conserved, having internal duplication in α-chain similar to human Hp2 and other ruminant species, which might have evolved separately as a convergent evolutionary process. Furthermore, the presence of extra Hp CCP domain possibly in all ruminants may have an effect during dimerization of molecule in these species.

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.

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.

POLYMORPHISM IN THE CODING REGION SEQUENCE OF GDF8 GENE IN INDIAN SHEEP.

The present study was undertaken to identify polymorphism in the coding sequence of GDF8gene across indigenous meat type sheep breeds. A 1647 bp sequence was generated, encompassing 208 bp of the 5'UTR, 1128 bp of coding region (exon1, 2 and 3) as well as 311 bp of 3'UTR. The sheep and goat GDF8 gene sequences were observed to be highly conserved as compared to cattle, buffalo, horse and pig. Several nucleotide variations were observed across coding sequence of GDF8 gene in Indian sheep. Three polymorphic sites were identified in the 5'UTR, one in exon 1 and one in the exon 2 regions. Both SNPs in the exonic region were found to be non-synonymous. The mutations c.539T > G and c.821T > A discovered in this study in the exon 1 and exon 2, respectively, have not been previously reported. The information generated provides preliminary indication of the functional diversity present in Indian sheep at the coding region of GDF8gene. The novel as well as the previously reported SNPs discovered in the Indian sheep warrant further analysis to see whether they affect the phenotype. Future studies will need to establish the affect of reported SNPs in the expression of the GDF8 gene in Indian sheep population.

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.

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.