Hematological and serum biochemical parameters and profiling of cytokine genes in lumpy skin disease in Vrindavani cattle

Lumpy skin disease (LSD) is a notifiable re-emerging transboundary viral disease of bovines that inflicts heavy losses in affected livestock farms. Genetic variations contribute substantially to the inter-individual differences in the immune response against disease agents. The present study aimed to evaluate the genetic basis of differential immune response in Vrindavani cattle by comparing the hematological, biochemical and cytokine genes' expression profiles of LSD-affected and unaffected animals. After 21 days of the outbreak at the farm, the animals were grouped as affected (those who developed symptoms) and unaffected/healthy (those who did not). Standard hematological and biochemical parameters were evaluated in both the groups. Expression profiling of important Th1 (IL2, INFG and GMCSF) and Th2 (IL4, IL6 and IL10) cytokines was also performed via a relative quantification approach using real-time PCR. Erythrogram and leucogram analyses revealed significant differences in total leucocyte count (TLC: 14.18 ± 0.74 versus 11.38 ± 0.68 x10³/µL), hemoglobin (Hb: 8.66 ± 0.42 versus 10.84 ± 0.17 g%) and percentage of neutrophils (46.40 ± 1.98 versus 35.40 ± 2.11%), lymphocytes (49.40 ± 1.99 versus 62.40 ± 1.86) and monocytes (4.20 ± 0.37 versus 2.40 ± 0.40) between the affected and healthy animals, respectively. The production of liver enzymes (SGOT and SGPT) was significantly higher in affected animals (74.18 ± 4.76 and 59.51 ± 2.75) when compared to the healthy counterparts (65.95 ± 9.18 and 39.21 ± 3.31). The expression profiling of Th1 and Th2 cytokines revealed significant differences between the two groups, except IL10. The expression of IL2, GMCSF and IL6 were upregulated in healthy animals while that of INFG, IL4 and IL10 were upregulated in LSD-affected animals. The highest abundance was observed for IL2 transcripts in healthy animals among all assessed cytokines with log₂fold change of 1.61 as compared to affected counterparts. Overall, the immune response in healthy animals (after exposure to LSD virus) was predominated by the expression of Th1 cell proliferation and there was an increased production of pro-inflammatory cytokines as compared to the affected counterparts. The results revealed that the effective immune response to LSD in cattle consists of changes in hematological and biochemical parameters and altered expression profile of cytokines with enhanced phagocytosis and lymphocyte recruitment. Furthermore, optimal expression of Th1 cytokines is required for maintaining optimal health against infectious insult with LSD virus in cattle.

Systems biology under heat stress in Indian cattle

Transcriptome profiling of Vrindavani and Tharparkar cattle (n = 5 each) revealed that more numbers of genes were dysregulated in Vrindavani than in Tharparkar. A contrast in gene expression was observed with 18.9 % of upregulated genes in Vrindavani downregulated in Tharparkar and 17.8% upregulated genes in Tharparkar downregulated in Vrindavani. Functional annotation of genes differentially expressed in Tharparkar and Vrindavani revealed that the systems biology in Tharparkar is moving towards counteracting the effects due to heat stress. Unlike Vrindavani, Tharparkar is not only endowed with higher expression of the scavengers (UBE2G1, UBE2S, and UBE2H) of misfolded proteins but also with protectors (VCP, Serp1, and CALR) of naïve unfolded proteins. Further, higher expression of the antioxidants in Tharparkar enables it to cope up with higher levels of free radicals generated as a result of heat stress. In this study, we found relevant genes dysregulated in Tharparkar in the direction that can counter heat stress.

Ancestry informative markers derived from discriminant analysis of principal components provide important insights into the composition of crossbred cattle

The cost of SNP genotyping to screen different breeds and to estimate the exact proportion of ancestry level is quite high, which can be compensated through deriving a small panel of ancestry informative markers (AIMs). Hence, we carried out the present study to provide an insight into ancestry level inferred from a panel of informative markers in the crossbred Vrindavani population developed at ICAR-IVRI, India. We have performed a new method i.e., discriminant analysis of principal components (DAPC) for the first time on the dataset of Vrindavani cattle. To confirm our method, we had performed DAPC on two other well-known crossbred cattle, i.e., Frieswal and Beefmaster. Three sets of panels (500, 1000 and 2000 markers) were tested for clustering of individuals. Among all the panels, we found the panel (1000 markers) with DAPC based contribution method was of the smallest size and comparatively of the highest accuracy.

Distribution analysis of single nucleotide polymorphisms related to feet and legs and their association with lameness in Vrindavani cattle

Lameness is one of the most serious economic problems in dairy animals. The present study was undertaken to identify the SNPs associated with lameness in Vrindavani, a synthetic crossbred cattle developed in India. A total of 187 cows were screened for lameness at Cattle and Buffalo Farm, Indian Veterinary Research Institute, Izatnagar, India and 16.04% prevalence was noticed. The effect of age, body weight, parity, and peak yield was not significant on lameness score. Heritability estimate of lameness score was low (0.07 ± 0.01). A total of 90 Vrindavani cows were genotyped by PCR-RFLP for 10 SNPs that were reported to have association with feet and legs in cattle. Five SNPs displayed significant departure from Hardy Weinberg equilibrium. Varying levels of heterozygosity, polymorphic information content, and allelic diversity were noticed at different SNPs. Odds ratio and least squares analysis revealed a significant association (P < 0.05) of three SNPs with lameness in the herd. Cows with AG, TT, and GG genotypes respectively at rs41632254, rs41603160, and rs41636945 loci were less prone to lameness. These genetic markers would certainly aid in development of an early detection system for lameness.

Revelation of genomic breed composition in a crossbred cattle of India with the help of Bovine50K BeadChip

The aim of the present study was to assess the population structure and admixture levels in the Vrindavani composite population in India by using Bovine50KSNP BeadChip data. Genotypic data were generated for randomly selected animals (n = 72) of Vrindavani population and the data for parental breeds i.e., Hariana (n = 10), Holstein-Friesian (n = 63), Jersey (n = 28) and Brown Swiss (n = 22) were retrieved from a public repository. The indices of population structure were calculated using PLINK software and R-program. The merged dataset was analysed for assessing admixture levels and population stratification using three different approaches i.e., principal component analysis (PCA), multi-dimensional scaling (MDS) approach and the model-based approach in STRUCTURE software. The average minor allele frequency (MAF) value for Vrindavani population was estimated to be 0.235. Vrindavani population was found to possess an average ancestry of 39.5, 22.9, 26.9, and 10.7% inheritance levels from Holstein Friesian, Jersey, Hariana and Brown Swiss cattle breeds, respectively.

Association of ATP1A1 gene polymorphism with thermotolerance in Tharparkar and Vrindavani cattle

Aim:

One of the major biochemical aspects of thermoregulation is equilibrium of ion gradient across biological membranes. Na⁺/K⁺-ATPase, a member of P type-ATPase family, is a major contributor to the mechanism that actively controls cross-membrane ion gradient. Thus, we examined ATP1A1 gene that encodes alpha-1 chain of Na⁺/K⁺-ATPase, for genetic polymorphisms.

Materials and Methods:

A total of 100 Vrindavani (composite cross strain of Hariana x Holstein-Friesian/Brown Swiss/Jersey) and 64 Tharparkar (indigenous) cattle were screened for genetic polymorphism in ATP1A1 gene, using polymerase chain reaction single-strand conformation polymorphism and DNA sequencing. For association studies, rectal temperature (RT) and respiration rate (RR) of all animals were recorded twice daily for 3 seasons.

Results:

A SNP (C2789A) was identified in exon 17 of ATP1A1 gene. Three genotypes namely CC, CA, and AA were observed in both, Vrindavani and Tharparkar cattle. The gene frequencies in Tharparkar and Vrindavani for allele A were 0.51 and 0.48, and for allele C were 0.49 and 0.52, respectively, which remained at intermediate range. Association study of genotypes with RT and RR in both cattle population revealed that the animals with genotype CC exhibited significantly lower RT and higher heat tolerance coefficient than CA and AA genotypes.

Conclusion:

Differential thermoregulation between different genotypes of ATP1A1 gene indicate that the ATP1A1 gene could be potentially contributing to thermotolerance in both, Tharparkar, an indigenous breed and Vrindavani, a composite crossbred cattle.