Udder and teat morphometry in relation to clinical mastitis in dairy cows.

The present study was undertaken to explore the association of occurrence of clinical mastitis with udder and teat morphometry in cattle. The relationship of nine udder and eight teat morphometric traits with the occurrence of clinical mastitis was analysed using a binomial logistic regression model. A significant association of the udder and teat morphometric traits with the occurrence of mastitis was observed both in Karan Fries and Sahiwal cattle. The rear udder height, rear udder width, udder width and fore udder attachment were found to be the important udder morphometric traits associated with the occurrence of mastitis. In both groups, mastitic animals had broader rear udder width, lower rear udder height and loose fore udder attachment with lengthy and thicker teats. The inclusion of udder and teat morphometric traits in the selection and breeding program of Karan Fries and Sahiwal cattle may help in selection of animals for mastitis resistance.

Reduced representation approach for identification of genome-wide SNPs and their annotation for economically important traits in Indian Tharparkar cattle.

The present study was carried out in Tharparkar cattle for identification of genome-wide SNPs and microsatellites, and then annotate the identified high-quality SNPs to milk production, fertility, carcass, adaptability and immune response of economically important traits. A total of 146,011 SNPs were identified with respect to Bos taurus reference genome which are indicus specific, out of which 10,519 SNPs were found to be novel. Similarly, a total of 87,047 SNPs were identified with respect to Bos indicus reference genome. After final annotation of SNPs identified with respect to Bos indicus reference genome, 2871 SNPs were found to be associated in 383 candidate genes having to do with milk production, fertility, carcass, immune response and adaptability traits. Following that, 2571 microsatellites were identified. The information mined from the data might be of importance for the future breed improvement programs, conservation efforts and for enhancing the SNPs density of the existing bovine SNP chips.

Genomic divergence reveals unique populations among Indian Yaks.

The present study focused upon identification of genome-wide SNPs through the reduced representation approach and to study the genomic divergence of the Indian yak populations. A total of 80 samples belonging to Arunachali yak (N = 20), Himachali yak (N = 20), Ladakhi yak (N = 20) and Sikkimi yak (N = 20) of India were used in the study. The results of the study revealed a total of 579575 high quality SNPs along with 50319 INDELs in the Indian yaks. The observed heterozygosity was found to be high in Himachali yak, followed by Arunachali yak, Ladakhi yak and Sikkimi yaks. The Sikkimi yaks was found to be genetically distant, followed by Ladakhi yaks which was observed to have some few individuals from Arunachali and Himachali yaks. Arunachali and Himachali yaks are found to get clustered together and are genetically similar. The study provides evidence about the genomic diversity in the Indian yak populations and information generated in the present study may help to formulate a suitable breeding plan for endangered Indian yaks. Moreover, the unique yak populations identified in the study will further help to focus attention for future characterization and prioritization of the animals for conservation purposes through the ddRAD approach.

WITHDRAWN:Principal component analysis of linear udder type traits and their relationship with milk yield and composition in indigenous Sahiwal cattle.

Ahead of Print article withdrawn by publisher.

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.

Genome-wide discovery of SNPs in candidate genes related to production and fertility traits in Sahiwal cattle.

The present study was carried out to identify genome-wide genetic markers and variants in candidate genes for production and reproduction traits in Sahiwal cattle using a cost-effective reduced representation sequencing method. A total of 258,231 genome-wide SNPs were identified in Sahiwal cattle with reference to Bos indicus genome, of which 150,231 were novel SNPs. Among the high-confidence SNPs identified, 91.86% and 27.30% were genotyped in 50% and 100% of the samples. Mapping of the identified SNPs revealed 525 SNPs in candidate genes related to production traits while 333 SNPs were mapped to candidate genes related to reproduction traits. The SNPs identified in this study will facilitate further insights on tropical adaptation, domestication history and population structure of indigenous cattle. The variants in candidate genes identified in this study will serve as useful genetic tools, in the quest for phenotype modifying nucleotide change and help in designing appropriate genetic improvement programs.

Development of PCR based assays for detection of lethal Holstein haplotype 1, 3 and 4 in Holstein Friesian cattle.

Holstein haplotype (HH) 1, 3 and 4 are lethal mutations, responsible for early embryonic losses in Holstein Friesian (HF) cattle, worldwide. Three PCR based assays - tetra Amplification Refractory Mutation System PCR, PCR primer induced restriction analysis and PCR-restriction fragment length polymorphism techniques for screening of HH1, 3 and 4, respectively were developed and validated. During screening, six among 60 HF bulls were found as carrier for either of three mutations. These PCR assays are highly accurate and reproducible and can be used for screening of the haplotypes in HF cattle.

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.

Genomewide identification and annotation of SNPs in Bubalus bubalis.

The present study was carried out to identify and annotate the genome wide SNPs in Murrah buffalo genome. A total of 21.2 million raw reads from 4 pooled female Murrah buffalo samples were obtained using restriction enzyme digestion followed by sequencing with Illumina Hiseq 2000. After quality filtration, the reads were aligned to Murrah buffalo genome (ICAR-NBAGR) and Water buffalo genome (UMD_CASPUR_WB_2.0) which resulted in 99.37% and 99.67% of the reads aligning, respectively. A total of 130,688 high quality SNPs along with 35,110 indels were identified versus the Murrah bufffalo genome. Similarly 219,856 high quality SNPs along with 15,201 indels were identified versus the Water buffalo genome. We report 483 SNPs in 66 genes affecting Milk Production, 436 SNPs in 38 genes affecting fertility and 559 SNPs in 72 genes affecting other major traits. The average genome coverage was 13.4% and 14.8% versus the Murrah and Water buffalo genomes, respectively.

Novel SNP identification in exon 3 of HSP90AA1 gene and their association with heat tolerance traits in Karan Fries (Bos taurus × Bos indicus) cows under tropical climatic condition.

Heat shock proteins (HSPs) act as molecular chaperones those are preferentially transcribed in respose to heat stress and the polymorphism in HSP genes associated with heat tolerance traits in cows. HSP90AA1 gene has been mapped on Bos taurus autosome 21 (BTA-21) and spans nearly 5368 bp comprising of 11 exons out of which the first exon does not translate. The present study was done on Karan Fries (5/8 HF × 3/8 Tharparkar) cows reared in tropical climate with the objectives of identifying single-nucleotide polymorphisms (SNPs) in targeted regions (exon 3) of HSP90AA1 gene and analyzing their association with heat tolerance traits in Karan Fries cows. Respiration rate (RR) and rectal temperature (RT) were recorded once daily for four consecutive days during probable extreme hours in different seasons or temperature humidity index (THI), viz., winter, spring, and summer. For detecting single-nucleotide polymorphisms, sequence data were analyzed using BioEdit software (version 7.2). Comparative sequence analysis of HSP90AA1 gene showed point mutation, viz., g.1209A>G (exon 3) as compared to Bos taurus (NCBI Ref Seq: AC_000178.1). Association analysis indicated that THI was influenced (P < 0.01) by RR, RT, and HTC. Similarly, SNPs at locus g.1209A>G were categorized into three genotypes, i.e., AA, AG, and GG, and the least squares means (LSMEANS) of RR, RT, and HTC for GG (homozygous) genotype were significantly lower (P < 0.01) than AA (homozygous) and AG (heterozygous) genotypes. These findings may partly suggest that cows with GG genotypes were favored for heat tolerance trait, which can be used as an aid to selection for thermo-tolerance Karan Fries cows for better adaptation in subtropical and tropical hot climate.

Novel SNPs in HSPB8 gene and their association with heat tolerance traits in Sahiwal indigenous cattle.

Heat shock proteins (HSPs) are expressed in response to heat stress, and the polymorphism in HSP genes at single-nucleotide level has been reported to be associated with heat tolerance and production performance traits in cattle. HSPB8 gene has been mapped on Bos taurus autosome 17 (BTA-17) spanning nearly 13,252 bp and comprising three exons and two introns. The present study was conducted in Sahiwal cows (n = 108) reared in subtropical climate with the objectives to identify SNPs in all three exons and part of intron 1 of HSPB8 gene and to analyze their association with heat tolerance traits in Sahiwal cows. Respiration rate (RR) and rectal temperature (RT) were recorded once during probable extreme hours in different seasons or Temperature-Humidity Index (THI), i.e., winter, spring, and summer. Heat tolerance coefficient (HTC) was also calculated to check the adaptability of the animals during the period of heat stress. The comparative sequence analysis revealed a total two single-nucleotide polymorphisms (SNPs), i.e., g.507G>A in exon 1 and g.881T>C in intron 1 of HSPB8 gene. Out of these two identified SNPs, only one SNP, i.e., g.507G>A, was found to be significantly associated with heat tolerance indicator traits (RR, RT, and HTC) in Sahiwal cows. The perusal of results across different seasons showed the significant (P < 0.01) difference in RR, RT, and HTC between winter, spring, and summer seasons. RR, RT, and HTC were found to be significantly lower (P < 0.01) in GA as compared to GG genotype of g.507G>A SNP of HSPB8 gene. However, in case of another SNP, i.e., g.881T>C, located on intron 1, the RR, RT, and HTC were having non-significant association with the different genotypes, i.e., TT and TC. These findings may partly suggest that GA genotype of SNP g.507G>A of HSPB8 gene has a probable role in heat tolerance in Sahiwal cattle and can therefore be utilized as a marker in propagation of thermo-tolerance cattle in hot tropical and subtropical climate. Nevertheless, the involvement of other regulatory mechanisms cannot be overruled.

Genetic polymorphisms within exon 3 of heat shock protein 90AA1 gene and its association with heat tolerance traits in Sahiwal cows.

AIM: The present study was undertaken to identify novel single nucleotide polymorphism (SNP) in Exon 3 of HSP90AA1 gene and to analyze their association with respiration rate (RR) and rectal temperature (RT) in Sahiwal cows. MATERIALS AND METHODS: The present study was carried out in Sahiwal cows (n=100) with the objectives to identify novel SNP in exon 3 of HSP90AA1 gene and to explore the association with heat tolerance traits. CLUSTAL-W multiple sequence analysis was used to identify novel SNPs in exon 3 of HSP90AA1 gene in Sahiwal cows. Gene and genotype frequencies of different genotypes were estimated by standard procedure POPGENE version 1.32 (University of Alberta, Canada). The significant effect of SNP variants on physiological parameters, e.g. RR and RT were analyzed using the General Linear model procedure of SAS Version 9.2. RESULTS: The polymerase chain reaction product with the amplicon size of 450 bp was successfully amplified, covering exon 3 region of HSP90AA1 gene in Sahiwal cows. On the basis of comparative sequence analysis of Sahiwal samples (n=100), transitional mutations were detected at locus A1209G as compared to Bos taurus (NCBI GenBank AC_000178.1). After chromatogram analysis, three genotypes AA, AG, and GG with respective frequencies of 0.23, 0.50, and 0.27 ascertained. RR and RT were recorded once during probable extreme hours in winter, spring, and summer seasons. It was revealed that significant difference (p<0.01) among genetic variants of HSP90AA1 gene with heat tolerance trait was found in Sahiwal cattle. The homozygotic animals with AA genotype had lower heat tolerance coefficient (HTC) (1.78±0.04(a)), as compared to both AG and GG genotypes (1.85±0.03(b) and 1.91±0.02(c)), respectively. The gene and genotype frequencies for the locus A1209G were ascertained. CONCLUSIONS: Novel SNP was found at the A1209G position showed all possible three genotypes (homozygous and heterozygous). Temperature humidity index has a highly significant association with RR, RT, and HTC in all the seasons. Perusal of results across different seasons showed the significant (p<0.01) difference in RR, RT, and HTC among winter, spring, and summer seasons. Genetic association with heat tolerance traits reveals their importance as a potential genetic marker for heat tolerance traits in Sahiwal cows.