Comprehensive selection signature analyses in dairy cattle exploiting purebred and crossbred genomic data

The main objective of the current research was to locate, annotate, and highlight specific areas of the bovine genome that are undergoing intense positive selection. Here, we are analyzing selection signatures in crossbred (Bos taurus X Bos indicus), taurine (Bos taurus), and indicine (Bos indicus) cattle breeds. Indicine cattle breeds found throughout India are known for their higher heat tolerance and disease resilience. More breeds and more methods can provide a better understanding of the selection signature. So, we have worked on nine distinct cattle breeds utilizing seven different summary statistics, which is a fairly extensive approach. In this study, we carried out a thorough genome-wide investigation of selection signatures using bovine 50K SNP data. We have included the genotyped data of two taurine, two crossbreds, and five indicine cattle breeds, for a total of 320 animals. During the 1950s, these indicine (cebuine) cattle breeds were exported with the aim of enhancing the resilience of taurine breeds in Western countries. For this study, we employed seven summary statistics, including intra-population, i.e., Tajima's D, CLR, iHS, and ROH and inter-population statistics, i.e., FST, XP-EHH, and Rsb. The NCBI database, PANTHER 17.0, and CattleQTL database were used for annotation after finding the important areas under selection. Some genes, including EPHA6, CTNNA2, NPFFR2, HS6ST3, NPR3, KCNIP4, LIPK, SDCBP, CYP7A1, NSMAF, UBXN2B, UGDH, UBE2K, and DAB1, were shown to be shared by three or more different approaches. Therefore, it gives evidence of the most intense selection in these areas. These genes are mostly linked to milk production and adaptability traits. This study also reveals selection regions that contain genes which are crucial to numerous biological functions, including those associated with milk production, coat color, glucose metabolism, oxidative stress response, immunity and circadian rhythms.

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.

Uncovering genes underlying coat color variation in indigenous cattle breeds through genome-wide positive selection

The identification of candidate genes related to pigmentation and under selective sweep provides insights into the genetic basis of pigmentation and the evolutionary forces that have shaped this variation. The selective sweep events in the genes responsible for normal coat color in Indian cattle groups are still unknown. To find coat color genes displaying signs of selective sweeps in the indigenous cattle, we compiled a list of candidate genes previously investigated for their association with coat color and pigmentation. After that, we performed a genome-wide scan of positive selection signatures using the BovineSNP50K Bead Chip in 187 individuals of seven indigenous breeds. We applied a wide range of methods to find evidence of selection, such as Tajima's D, CLR, iHS, varLD, ROH, and FST. We found a total of sixteen genes under selective sweep, that were involved in coat color and pigmentation physiology. These genes are CRIM1 in Gir, MC1R in Sahiwal, MYO5A, PMEL and POMC in Tharparkar, TYRP1, ERBB2, and ASIP in Red Sindhi, MITF, LOC789175, PAX3 and TYR in Ongole, and IRF2, SDR165 and, KIT in Nelore, ADAMTS19 in Hariana. These genes are related to melanin synthesis, the biology of melanocytes and melanosomes, and the migration and survival of melanocytes during development.

Properties and physiological effects of dietary fiber-enriched meat products: a review

Meat is a rich source of high biological proteins, vitamins, and minerals, but it is devoid of dietary fiber, an essential non-digestible carbohydrate component such as cellulose, hemicellulose, pectin, lignin, polysaccharides, and oligosaccharides. Dietary fibers are basically obtained from various cereals, legumes, fruits, vegetables, and their by-products and have numerous nutritional, functional, and health-benefiting properties. So, these fibers can be added to meat products to enhance their physicochemical properties, chemical composition, textural properties, and organoleptic qualities, as well as biological activities in controlling various lifestyle ailments such as obesity, certain cancers, type-II diabetes, cardiovascular diseases, and bowel disorders. These dietary fibers can also be used in meat products as an efficient extender/binder/filler to reduce the cost of production by increasing the cooking yield as well as by reducing the lean meat content and also as a fat replacer to minimize unhealthy fat content in the developed meat products. So, growing interest has been observed among meat processors, researchers, and scientists in exploring various new sources of dietary fibers for developing dietary fiber-enriched meat products in recent years. In the present review, various novel sources of dietary fibers, their physiological effects, their use in meat products, and their impact on various physicochemical, functional, and sensory attributes have been focused.

Traditional halal meat production without stunning versus commercial slaughter with electrical stunning of slow-growing broiler chicken: impact on meat quality and proteome changes

Impact of traditional halal meat production without stunning (NST) and commercial slaughter with electrical stunning (ST) of 100 slow-growing broiler chicken on blood plasma and different biochemical, enzymatic, hormonal, meat quality, and proteomic changes was evaluated. The results revealed lower (P < 0.05) postmortem pH values and higher redness (a*) scores for ST samples relative to NST group. Myofibrillar fragmentation index and bleeding efficiency (%) were lower (P < 0.05) in ST compared to NST samples. The ST group had higher (P < 0.05) creatinine, total protein, alanine aminotransferase (ALT), and triiodothyronine (T3) than NST group, however, no difference (P > 0.05) in blood glucose, lactate dehydrogenase (LDH), creatine kinase (CK), thyroxine (T4), cortisol, and aspartate aminotransferase (AST) was observed relative to NST samples. The 2-dimensional gel electrophoresis (2-DE) coupled to MALDI-TOF MS of meat samples has identified 14 differentially abundant proteins between 2 groups. Proteins demonstrating positive correlation with stress namely adenylate kinase isoenzyme-1, Rho guanine nucleotide exchange factor (NST), and apolipoprotein A-I (ST) were overabundant. From the current study, it is concluded that electrical stunning of broilers prior to slaughter or traditional halal slaughter without stunning does not adversely affect the meat quality.

Biomarker discovery and authentication of cold-slaughtered chicken through classical analytical procedures and mass spectrometry based proteomic approaches

1. This study evaluated the suitability of routine analytical procedures and used mass spectrometry-based proteomic approaches to distinguish meat from dead chicken/ cold-slaughtered birds (CS), electrically stunned and slaughtered birds, as per standard protocols (ES), and birds slaughtered according to halal guidelines (HS).2. Meat from CS birds had lower (P < 0.05) pH, water-holding capacity and higher (P < 0.05) lipid oxidation, haem iron content, residual blood and total viable counts relative to ES and HS meat indicating poor quality.3. The results demonstrated the presence of unique protein bands on SDS-PAGE only in CS meat that can be used for routine screening.4. Protein analysis using MALDI-TOF mass spectrometry identified haemoglobin subunit alpha-A and alpha-D; Adenylate kinase isoenzyme 1 as reliable and stable marker proteins for authentication of dead chicken meat under raw and cooked conditions and halal slaughtered chicken, respectively.5. The methods used may be employed by the food safety and regulatory agencies for regular screening of meat quality and to authenticate CS or HS chicken.

Temporal dysregulation of genes in Lamb testis cell during sheeppox virus infection

The present study was aimed to elucidate the host-virus interactions using RNA-Seq analysis at 1h and 8h of post-infection of SPPV in LTC. The differentially expressed genes (DEGs) and the underlying mechanisms linked to the host immune responses were obtained. The protein-protein interaction (PPI) network analysis and Ingenuity pathway analysis (IPA) illustrated the interaction between the DEGs and their involvement in cell signalling responses. Highly connected hubs viz. AURKA, CHEK1, CCNB2, CDC6, and MAPK14 were identified through PPI network analysis. IPA analysis showed that IL-6 and ERK5 mediated signalling pathways were highly enriched at both time points. The TP53 gene was identified to be the leading upstream regulator that directly responded to SPPV infection, resulting in downregulation at both time points. The study provides an overview of how the lamb testis genes and their underlying mechanisms link to growth and immune response during SPPV infection.

Selection and validation of suitable reference gene for qPCR gene expression analysis in lamb testis cells under Sheep pox virus infection

Virus infection alters host gene expression, therefore ideal and stable reference housekeeping genes are required to normalise the expression of other expressed host genes in quantitative real-time PCR (qRT-PCR). The suitable reference gene may vary in response to different viral infections in different hosts or cells. In the present study, we cultured primary lamb testis cells (LTC) and assessed the expression stability of seven widely used housekeeping genes (B2M, HMBS, HPRT1, HSP-90, POLR2A, 18s_RNA, GAPDH) as reference genes in Sheeppox virus (SPPV) infected and control (uninfected-0h) LTC at 0.5h, 4.0h, 8.0h, and 12.0h post-infection) using NormFinder, Bestkeeper, geNorm, and the comparative ΔCT method in RefFinder based on their expression levels. Analysis revealed that HSP90, 18s_RNA, HPRT, POLR2A, and B2M were the most stable genes from the panel in the individual analysis group in 0h, 0.5h, 4.0h, 8.0h, and 12.0h, respectively. Furthermore, B2M was shown to be the most stable reference gene in the combined control with the respective and overall infected groups, except the control group of 4.0hpi of SPPV infection. In this study, we selected the most suitable reference genes in LTC for particular time points of SPPV infection. The identified most suitable housekeeping gene can be used during normalization of expression of other targeted genes at aspecific time point of SPPV infection.

Identification of important genomic footprints using eight different selection signature statistics in domestic cattle breeds

In the present study, the population genomic data of different cattle breeds were explored to decipher the genomic regions affected due to selective events and reflected in the productive, reproductive, thermo-tolerance, and health-related traits. To find out these genomic deviations due to selective sweeps, we used eight different statistical tools (Tajima's D, Fu & Li's D*, CLR, ROH, iHS, F_ST, FLK, and hapFLK) on seven indigenous and five exotic cattle breeds. We further performed composite analysis by comparing their covariance matrix. Several candidate genes were found to be related to milk production (ADARB, WDR70, and CA8), reproductive (PARN, FAM134B2, and ZBTB20), and health-related traits (SP110, CXCL2, CLXCL3, CXCL5, IRF8, and MYOM1). The outcome of this investigation provides a basis for detecting selective sweeps that explain the genetic variation of traits. They may possess functional importance for multiple cattle breeds in different subcontinents. However, further studies are required to improve the findings using high-density arrays or whole-genome sequencing with higher resolution and greater sample sizes.

Characterisation and comparison of immune response mechanisms in an indigenous and a commercial pig breed after classical swine fever vaccination

The live attenuated classical swine fever (CSF) vaccine has been successfully used to prevent and control CSF outbreaks for 6 decades. However, the immune response mechanisms against the vaccine remain poorly understood. Moreover, very few reports exist regarding the breed differences in the response to CSF vaccine. In this study, we generated the peripheral blood mononuclear cell transcriptomes of indigenous Ghurrah and commercial Landrace pig breeds, before and 7 days after CSF vaccination. Subsequently, between and within-breed differential gene expression analyses were carried out. Results revealed large differences in pre-vaccination peripheral blood mononuclear cell transcriptome profiles of the two breeds, which were homogenised 7 days after vaccination. Before vaccination, gene set enrichment analysis showed that pathways related to antigen sensing and innate immune response were enriched in Ghurrah, while pathways related to adaptive immunity were enriched in Landrace. Ghurrah exhibited greater immunomodulation compared to Landrace following the vaccination. In Ghurrah, cell-cycle processes and T-cell response pathways were upregulated after vaccination. However, no pathways were upregulated in Landrace after vaccination. Pathways related to inflammation were downregulated in both the breeds after vaccination. Key regulators of inflammation such as IL1A, IL1B, NFKBIA and TNF genes were strongly downregulated in both the breeds after vaccination. Overall, our results have elucidated the mechanisms of host immune response against CSF vaccination in two distinct breeds and revealed common key genes instrumental in the global immune response to the vaccine.

Etiopathological and hematobiochemical profiles in goats with gastrointestinal disorders

Background and Aim

Gastrointestinal (GI) disorders in small ruminants limit production efficiency and productivity growth in the livestock sector, thereby directly preventing farmers from augmenting their income. This study aimed to provide detailed insight into the etiology, hematobiochemical parameters, and epidemiological risk factors of GI disorders in goats and to determine the pathology associated with the disorders.

Materials and Methods

Over the period of 2018-2019, 500 goats in and around Bhubaneswar, Odisha, India, were screened for GI disorders based on clinical signs. Blood samples from the control (n=10) and treatment (n=25) groups were collected for both hematological and serum biochemical alterations. Fecal examinations (n=220) were conducted for parasitic, bacterial, and virological assessments. Detailed necropsy and histopathological evaluations were conducted on 27 goats.

Results

The GI disorder prevalence rate and mortality rate among the 500 goats analyzed were 44.4% and 12.27%, respectively. Chi-square analysis showed a significantly higher occurrence of GI disorders among the goats that were between 6 months and 1.5 years old (58.72%), were of the Ganjam breed (45.49%), had a poor body condition (71.11%), and were housed with an earthen floor (55.22%). The most common etiological risk factor observed was parasitic infection (65.45%), followed by bacterial (18.18%) and mixed infection (9.54%). Blood analysis showed neutrophilia and eosinophilia in infected goats, in addition to anemia; significant decreases in total protein, globulin, albumin, and glucose levels; and significant increases in aspartate transaminase and alanine aminotransferase levels. The major histopathological findings were infiltration of mononuclear cells and desquamation of the intestinal and ruminal mucosa.

Conclusion

Stakeholders should focus not only on parasitic infections and other important etiological risk factors for GI disorders in goats but also on proper farming management practices to help enhance the income of farmers. The hematobiochemical alterations and pathomorphological changes reported in this study can be used by field veterinarians as guidelines for clinical evaluation and disease severity assessment.

Genome-wide detection of copy number variations in Tharparkar cattle

Copy number variations (CNVs) are major forms of genetic variation with an increasing importance in animal genomics. This study used the Illumina BovineSNP 50 K BeadChip to detect the genome-wide CNVs in the Tharparkar cattle. With the aid of PennCNV software, we noticed a total of 447 copy number variation regions (CNVRs) across the autosomal genome, occupying nearly 2.17% of the bovine genome. The average size of detected CNVRs was found to be 122.2 kb, the smallest CNVR being 50.02 kb in size, to the largest being 1,232.87 Kb. Enrichment analyses of the genes in these CNVRs gave significant associations with molecular adaptation-related Gene Ontology (GO) terms. Most CNVR genes were significantly enriched for specific biological functions; signaling pathways, sensory responses to stimuli, and various cellular processes. In addition, QTL analysis of CNVRs described them to be linked with economically essential traits in cattle. The findings here provide crucial information for constructing a more comprehensive CNVR map for the indigenous cattle genome.

Identification of genes affecting milk fat and fatty acid composition in Vrindavani crossbred cattle using 50 K SNP-Chip

The aim of this study was to identify candidate genes associated with milk fat per cent and fatty acid (FA) composition in Vrindavani cattle using the Illumina 50 K single-nucleotide polymorphism (SNP) array. After quality control, a total of 41,427 informative and high-quality SNPs were used for a genome-wide association study (GWAS) for milk fat percentage and 16 different types of fatty acids. Lactation stage, parity, test day milk yield, and proportion of exotic inheritance were included as fixed effects in the GWAS model. A total of 67 genome-wide significant (P < 1.20 × 10^-06) SNPs and 176 suggestive significant (P < 2.41 × 10^-05) SNPs were identified. Out of these, 15 SNPs were associated with more than one trait. The strongest associations were found on BTA14 for milk fat percentage and on BTA2 and BTA16 for polyunsaturated fatty acids. Several significant SNPs were identified close to or within the genes ELOVL6, FABP4, PMP2, PLIN1, MFGE8, GHRL2, and LDLRAD3 which are known to be associated with fat percentage and FA composition in dairy cattle breeds. This study is a step forward to better characterize the molecular mechanisms of phenotypic variation in milk fatty acids in a taurine-indicine composite cattle breed reared in tropical environments.

Genomic scans for selection signatures revealed candidate genes for adaptation and production traits in a variety of cattle breeds

Domestication and selection are the major driving forces responsible for the determinative genetic variability in livestock. These selection patterns create unique genetic signatures within the genome. BovineSNP50 chip data from 236 animals (seven indicine and five taurine cattle breeds) were analyzed in the present study. We implemented three complementary approaches viz. iHS (Integrated haplotype score), ROH (Runs of homozygosity), and F_ST, to detect selection signatures. A total of 179, 56, and 231 regions revealed 518, 277, and 267 candidate genes identified by iHS, ROH, and F_ST methods, respectively. We found several candidate genes (e.g., NCR3, ARID5A, HIST1H2BN, DEFB4, DEFB7, HSPA1L, HSPA1B, and DNAJB4) related to production traits and the adaptation of indigenous breeds to local environmental constraints such as heat stress and disease susceptibility. However, further studies are warranted to refine the findings using a larger sample size, whole-genome sequencing, and/or high density genotyping.

SNPs with intermediate minor allele frequencies facilitate accurate breed assignment of Indian Tharparkar cattle

Tharparkar cattle breed is widely known for its superior milch quality and hardiness attributes. This study aimed to develop an ultra-low density breed-specific single nucleotide polymorphism (SNP) genotype panel to accurately quantify Tharparkar populations in biological samples. In this study, we selected and genotyped 72 Tharparkar animals randomly from Cattle & Buffalo Farm of IVRI, India. This Bovine SNP50 BeadChip genotypic datum was merged with the online data from six indigenous cattle breeds and five taurine breeds. Here, we used a combination of pre-selection statistics and the MAF-LD method developed in our laboratory to analyze the genotypic data obtained from 317 individuals of 12 distinct breeds to identify breed-informative SNPs for the selection of Tharparkar cattle. This methodology identified 63 unique Tharparkar-specific SNPs near intermediate gene frequencies. We report several informative SNPs in genes/QTL regions affecting phenotypes or production traits that might differentiate the Tharparkar breed.

Signatures of Selection in Composite Vrindavani Cattle of India

Vrindavani is an Indian composite cattle breed developed by crossbreeding taurine dairy breeds with native indicine cattle. The constituent breeds were selected for higher milk production and adaptation to the tropical climate. However, the selection response for production and adaptation traits in the Vrindavani genome is not explored. In this study, we provide the first overview of the selection signatures in the Vrindavani genome. A total of 96 Vrindavani cattle were genotyped using the BovineSNP50 BeadChip and the SNP genotype data of its constituent breeds were collected from a public database. Within-breed selection signatures in Vrindavani were investigated using the integrated haplotype score (iHS). The Vrindavani breed was also compared to each of its parental breeds to discover between-population signatures of selection using two approaches, cross-population extended haplotype homozygosity (XP-EHH) and fixation index (F_ST). We identified 11 common regions detected by more than one method harboring genes such as LRP1B, TNNI3K, APOB, CACNA2D1, FAM110B, and SPATA17 associated with production and adaptation. Overall, our results suggested stronger selective pressure on regions responsible for adaptation compared to milk yield.

Selection and validation of differentially expressed metabolic and immune genes in weaned Ghurrah versus crossbred piglets

In the present investigation, differentially expressed genes (DEGs) were studied using RNA sequencing (RNA-seq) technique in porcine peripheral blood mononuclear cells (PBMC) of weaned Ghurrah and crossbred piglets at 3-month age. Transcriptomic analysis was done using three different packages, namely, EBSeq, DESeq2, and edgeR, to identify the DEGs between Ghurrah and crossbred piglets. Total 7717 DEGs were commonly identified by all three packages, out of which 4151 genes found to be up-regulated, and 3566 genes were down-regulated. Functional annotation of these DEGs indicated metabolism as the most commonly enriched category followed by the immune response. Genes related to metabolism and growth were up-regulated in crossbred piglets as compared with Ghurrah piglets, whereas immunity-related genes were up-regulated in Ghurrah piglets elucidating the disease resistance nature of this indigenous breed over crossbred counterparts. Further, eight DEGs, namely, LRP-1, ADCY4, ERRFI1, LDHD, ARG1, OASL, MGARP, and S100A8, were validated by qRT-PCR in a separate set of biological samples and found to be in concordance with RNA-seq results. Finding in the present study provides insight into genes and their molecular mechanisms governing difference in growth performance between Ghurrah and crossbred pigs.

Detection of recent infection of Japanese encephalitis virus in swine population using IgM ELISA: A suitable sentinel to predict infection in humans

Japanese encephalitis (JE) is a mosquito-borne flaviviral zoonotic disease and is one of the major causes of encephalitis in children. Swine, being an amplifier host of Japanese encephalitis virus (JEV), play an important role in its epidemiology. Therefore, early detection of either JEV or antibodies against JEV in swine is a feasible alternative for initiating necessary measures to prevent the spread of infection to humans. Since IgM antibodies appear early in swine sera, recombinant NS1 protein based indirect IgM ELISA was developed in the present study with the objective to know the recent infection of swine population with JEV. The relative diagnostic sensitivity and specificity of the developed ELISA was 95.34% and 98.6%, respectively. The developed ELISA was found to have excellent reproducibility on inter-laboratory and inter-institutional validation studies. A total of 3,027 field swine sera samples were screened using the developed ELISA and 488 samples were found positive for IgM against JEV with an overall sero-positivity of 16.12% in swine population of India. The highest sero-positivity was observed in swine population of Eastern zone of India which coincided with the maximum number of human JE cases reported from this zone during the same period. Further, antibody kinetics study revealed that the IgM antibodies against NS1 protein of JEV started appearing in swine sera at day 5 and disappeared completely by day 40. The IgG antibodies started appearing at day 7, and remained for more than 365 days indicating the suitability of IgM ELISA to know the recent infection of JEV. The developed IgM ELISA can be readily incorporated into surveillance programs for detection of JEV activity in swine population so that outbreaks in humans can be prevented by taking suitable preventive measures.

Genome-wide assessment of genetic diversity, linkage disequilibrium and haplotype block structure in Tharparkar cattle breed of India

Knowledge about genetic diversity is very essential for the management and sustainable utilization of livestock genetic resources. In this study, we presented a comprehensive genome-wide analysis of genetic diversity, ROH, inbreeding, linkage disequilibrium, effective population size and haplotype block structure in Tharparkar cattle of India. A total of 24 Tharparkar animals used in this study were genotyped with Illumina BovineSNP50 array. After quality control, 22,825 biallelic SNPs were retained, which were in HWE, MAF > 0.05 and genotyping rate >90%. The overall mean observed (H_O) and expected heterozygosity (H_E) were 0.339 ± 0.156 and 0.325 ± 0.129, respectively. The average minor allele frequency was 0.234 with a standard deviation of ± 0.131. We identified a total of 1832 ROH segments and the highest autosomal coverage of 13.87% was observed on chromosome 23. The genomic inbreeding coefficients estimates by F_ROH, F_HOM, F_GRM and F_UNI were 0.0589, 0.0215, 0.0532 and 0.0160 respectively. The overall mean linkage disequilibrium (LD) for a total of 133,532 pairwise SNPs measured by D' and r² was 0.6452 and 0.1339, respectively. In addition, we observed a gradual decline in effective population size over the past generations.

Signature of genome wide gene expression in classical swine fever virus infected macrophages and PBMCs of indigenous vis-a-vis crossbred pigs

The genetic basis of differential host immune response vis-à-vis transcriptome profile was explored in PBMCs of indigenous (Ghurrah) and crossbred pigs after classical swine fever vaccination and in monocyte derived macrophages (MDMs) challenged with virulent classical swine fever (CSF) virus. The humoral immune response (E2 antibody) was higher (74.87%) in crossbred than indigenous pigs (58.20%) at 21st days post vaccination (21dpv). The rate of reduction of ratio of CD4⁺/CD8⁺ was higher in crossbred pigs than indigenous pigs at 7th days post vaccination (7dpv). The immune genes IFIT1, IFIT5, RELA, NFKB2, TNF and LAT2 were up regulated at 7dpv in RNA seq data set and was in concordance during qRT-PCR validation. The Laminin Subunit Beta 1 (LAMB1) was significantly (p ≤ 0.05) down-regulated in MDMs of indigenous pigs and consequently a significantly (p ≤ 0.01) higher copy number of virulent CSF virus was evidenced in macrophages of crossbred pigs than indigenous pigs. Activation of LXR:RXR pathway at 60 h post infection (60hpi) in MDMs of indigenous versus crossbred pigs inhibited nuclear translocation of NF-κB, resulted into transrepression of proinflammatory genes. But it helped in maintenance of HDL level by lowering down cholesterol/LDL level in MDMs of indigenous pigs. The key immune genes (TLR2, TLR4, IL10, IL8, CD86, CD54, CASP1) of TREM1 signaling pathway were upregulated at 7dpv in PBMCs but those genes were downregulated at 60hpi in MDMs indigenous pigs. Using qRT-PCR, the validation of differentially expressed, immunologically important genes (LAMB1, OAS1, TLR 4, TLR8 and CD86) in MDMs revealed that expression of these genes were in concordance with RNA-seq data.

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.

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.

Comparative analysis of five different methods to design a breed-specific SNP panel for cattle

Single nucleotide polymorphisms (SNPs) have now replaced microsatellite markers in several species for various genetic investigations like parentage assignment, genetic breed composition, assessment for individuality and, most popularly, as a useful tool in genomic selection. However, such a resource, which can offer to assist breed identification in a cost-effective manner is still not explored in cattle breeding programs. In our study, we have tried to describe methods for reducing the number of SNPs to develop a breed-specific panel. We have used SNP data from Dryad open public access repository. Starting from a global dataset of 178 animals belonging to 10 different breeds, we selected five panels each comprising of similar number of SNPs using different methods i.e., Delta, Pairwise Wright's F_ST, informativeness for assignment, frequent item feature selection (FIFS) and minor allele frequency-linkage disequilibrium (MAF-LD) based method. MAF-LD based method has been recently developed by us for construction of breed-specific SNP panels. The STRUCTURE software analysis of MAF-LD based method showed appropriate clustering in comparison to other panels. Later, the panel of 591 breed-specific SNPs was called to their respective breeds using Venny 2.1.0 and UGent web tools software. Breed-specific SNPs were later annotated by using various Bioinformatics softwares.

Development of recombinant nonstructural 1 protein based indirect enzyme linked immunosorbent assay for sero-surveillance of Japanese encephalitis in swine

Japanese encephalitis virus (JEV) causes severe neurological disease in humans, especially among children. The disease is endemic in several South Asian countries including India. Swine play a major role as amplifier host for JEV and act as a source of infection to humans through mosquito bite. Early detection of either virus or antibodies in swine will aid to undertake control measures to prevent virus spread to humans. Swine seldom show symptoms of JEV infection and the viraemic phase lasts for a short period of 3 to 4 days indicating the potential of detection of antibodies, which remain for relatively longer period, as a suitable alternative. Cost effective and sensitive assays for the detection of JEV antibodies in swine are not available indigenously. Hence, we have developed a recombinant nonstructural protein 1 (rNS1) based enzyme linked immunosorbent assay for the detection of IgG antibodies against JEV in swine. The test is robust, highly sensitive (91%), specific (97%), reproducible and affordable. Field validation of the assay was done by screening 3628 swine Serum samples collected from different parts of India. The overall sero-positivity was found to be 32.22%. The developed ELISA can be readily incorporated into surveillance programs for detection of Japanese encephalitis virus activity in swine population thereby aiding in prediction of outbreaks in humans.

Genome-wide integrated analysis of miRNA and mRNA expression profiles to identify differentially expressed miR-22-5p and miR-27b-5p in response to classical swine fever vaccine virus

The present study was conducted to identify the differentially expressed miRNAs (DE miRNAs) in the peripheral blood mononuclear cells of crossbred pigs in response to CSF vaccination on 7 and 21 days of post vaccination as compared to unvaccinated control (0 dpv). Simultaneously, set of miRNA was predicted using mRNA seq data at same time point. The proportion of CD4^-CD8⁺ and CD4⁺CD8⁺ increased after vaccination, and the mean percentage inhibition was 86.89% at 21 dpv. It was observed that 22 miRNAs were commonly expressed on both the time points. Out of predicted DE miRNAs, it was found that 40 and 35 DE miRNAs were common, obtained from miRNA seq analysis and predicted using mRNA seq data on 7 dpv versus 0 dpv and 21 dpv versus 0 dpv respectively. Two DE miRNAs, ssc-miR-22-5p and ssc-miR-27b-5p, were selected based on their log₂ fold change and functions of their target genes in immune process/pathway of viral infections. The validations of DE miRNAs using qRT-PCR were in concordance with miRNA seq analysis. Two set of target genes, CD40 and SWAP70 (target gene of ssc-miR-22-5p) and TLR4 and Lyn (target gene of ssc-miR-27b-5p), were validated and were in concordance with results of RNA seq analysis at a particular time point (except TLR4). The first report of genome-wide identification of differentially expressed miRNA in response to live attenuated vaccine virus of classical swine fever revealed miR-22-5p and miR-27b-5p were differentially expressed at 7 dpv and 21 dpv.

PCR-SSCP analysis of MDGI gene and its association with milk production traits in river buffalo (Bubalus bubalis)

In this study, we investigated the genetic variation within 3'UTR of Mammary-Derived Growth Inhibitor (MDGI) gene of buffalo using PCR-SSCP and sequencing; and also analyzed association of polymorphism with the milk production traits. The study revealed two conformational patterns, 'A' and 'B' among 234 Mehsana buffaloes maintained with their records in the field and at farm. The frequency of SSCP variant 'A' was found to be invariably high in the buffalo population under study. Further, association analysis of SSCP variants with various milk production and milk quality traits indicated no significant effect on any of the traits investigated. Sequencing of SSCP variant 'A' showed homozygous G/G and A/A and 'B' had heterozygous G/C and A/G at positions +124 and +140 respectively, in the 3'UTR of buffalo MDGI. The preliminary results showed the substantial variations in the distribution of SSCP variants' frequencies within Mehsana buffaloes, however these variants had non-significant association with milk yield, fat yield and fat percentage in Mehsana buffaloes.

Hemato-biochemical alterations and urinalysis in dogs suffering from benign prostatic hyperplasia

Aim:

The study was designed to evaluate the hemato-biochemical alterations, urinalysis along with histomorphological and histological changes of prostate glands in dogs affected with benign prostatic hyperplasia (BPH) in and around Bhubaneswar, Odisha, India.

Materials and Methods:

In toto, 445 dogs presented to the Teaching Veterinary Clinical Complex of the College of Veterinary Sciences and Animal Husbandry, one Government Veterinary Hospital and two pet clinics in and around Bhubaneswar screened for the presence of BPH. Most of the 57 dogs were 6 years and above as reported by the owners. Only 57 dogs found positive for BPH basing on the presence of typical clinical signs subjected for a detailed hemato-biochemical study. Most of the 57 dogs were 6 years and above as reported by the owners. Routine and microscopic urinalyses were done as per the routine procedure. Histomorphological evaluations of prostate glands were done through manual rectal palpation. Histological examinations of prostate tissue sections of two dead dogs were conducted with routine hematoxylin and eosin stain.

Results:

The study revealed about 12.8% (57/445) of dogs was suffering from BPH. Typical clinical signs - such as passing small thin tape-shaped feces, holding tail away from backward, tenesmus, and straining during urination and defecation - were seen in most of the cases. Urine samples of affected dogs were positive for glucose, occult blood, and protein. A significant decrease in lymphocytes and increase in eosinophil counts in dogs with BPH was recorded. Serum biochemical analysis showed a nonsignificant increase in creatinine and blood urea nitrogen with a significant decrease in total protein, albumin, globulin, A:G ratio. Histology of prostate glands collected during postmortem was characterized by fibrosis of prostate gland, and hyperplasia of the acinar epithelium.

Conclusions:

High rate of the prevalence of BPH in dogs poses an alarming condition which if diagnosed at an early stage can certainly prolong the longevity of the dogs.

Identification of genetic variation in NOD-like receptor 2 gene and influence of polymorphism on gene structure and function in buffalo (Bubalus bubalis)

Nucleotide-binding oligomerization domain (NOD)-like receptor 2 is one of the important mediators of innate as well as adaptive immune response to microbial infections. In this study, NOD-like receptor-2 was characterized by determining the full gene sequence and analyzing genetic diversity in Indian buffaloes. Sequence analysis of buffalo NOD2 revealed 3042 nucleotides long ORF, encoding 1013 amino acids from 12 exons. Domain structure analysis indicated existence of 8 leucine-rich repeat (LRR) domains in buffalo, cattle, sheep and mouse, along with central NACHT/NOD domain and two N-terminal CARD domains. Comparative sequence analysis among different buffalo breeds identified 46 polymorphic sites in NOD2 gene. Among coding region SNPs, 10 were non-synonymous, 7 synonymous and 3 were present in 5'UTR. Genotyping of two nsSNPs, revealed significant differences in the allele frequencies, distinguishing swamp and riverine buffaloes, having different utilities. Association analysis with mastitis in dairy buffaloes indicated significant variation in allelic frequencies at G1135A locus, between mastitis affected and non-affected animals. Further, NOD2 gene expression was quantified in different riverine buffalo tissues, using real-time PCR and lymph node displayed highest expression, compared to others organs included in the study. Overall, the study revealed buffalo NOD2 gene attributes, important to understand species specific immune response in ruminants.

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.

EvaGreen-based Multiplex Real-time PCR Assay for Rapid Differentiation of Wild-Type and Glycoprotein E-Deleted Bovine Herpesvirus-1 Strains

Bovine herpesvirus-1 (BoHV-1) is an important viral pathogen causing significant economic losses to the cattle industry. Glycoprotein E-deleted marker vaccines form the basis for BoHV-1 control programs widely, wherein detection and differentiation of wild-type and gE-deleted vaccine strains is of crucial importance for proper disease management. In the present study, we report an EvaGreen-based multiplex real-time polymerase chain reaction (EGRT-PCR) assay for rapid differentiation of wild-type and glycoprotein E-deleted strains of BoHV-1. The EGRT-PCR assay could simultaneously detect two viral genes (glycoprotein B and E) and an internal positive control gene (bovine growth hormone- bGH), in a single-tube reaction. The analytical sensitivity of the EGRT-PCR assay was as little as 10 copies of the BoHV-1 DNA per reaction. The modified real-time PCR assay could successfully differentiate wild-type and gE-deleted BoHV-1 strains based on gene specific melting temperatures (T_m) peaks. Our results have shown that the EGRT-PCR developed in this study might prove to be a promising tool in disease management by enabling rapid differentiation of wild-type and gE-deleted strains of BoHV-1.

A double-stranded probe coupled with isothermal amplification for qualitative and quantitative detection of avian reovirus

We applied a probe-based real-time loop-mediated isothermal amplification (Cy5-RTqLAMP) technique targeting the avian reovirus (ARV) S3 gene to develop a rapid, sensitive, and specific method for virus detection and quantification. This test specifically detected the presence of ARV, but not other viruses or bacteria present in clinical or artificially spiked samples, including Newcastle disease virus, infectious bursal disease virus, fowl adenovirus, Marek's disease virus, Escherichia coli, and Salmonella spp. This test can detect ARV in less than one hour with an analytical sensitivity of 10 viral gene copies and 1 fg of total cDNA. The Cy5-RTqLAMP does not yield false positive results and is 100 times more sensitive than conventional PCR. This test was shown to be able to detect the presence of ARV in clinical samples. A similar strategy may be used for detection of other important human and animal viral pathogens.

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.

Development and quality evaluation of dehydrated chicken meat rings using spent hen meat and different extenders

It is recommended that for effective utilization of spent hen meat, it should be converted into value added or shelf stable meat products. Since we are lacking in cold chain facilities, therefore there is imperative need to develop shelf stable meat products. The present study was envisaged with the objective to develop dehydrated chicken meat rings utilizing spent hen meat with different extenders. A basic formulation and processing conditions were standardized for dehydrated chicken meat rings. Extenders such as rice flour, barnyard millet flour and texturized soy granule powder at 5, 10 and 15 % levels were incorporated separately replacing the lean meat in pre standardized dehydrated chicken meat ring formulation. On the basis of physico-chemical properties and sensory scores optimum level of incorporation was adjudged as 10 %, 10 % and 5 % for rice flour, barnyard millet flour and texturized soy granule powder respectively. Products with optimum level of extenders were analysed for physico-chemical and sensory attributes. It was found that a good quality dehydrated chicken meat rings can be prepared by utilizing spent hen meat at 90 % level, potato starch 3 % and refined wheat flour 7 % along with spices, condiments, common salt and STPP. Addition of an optimum level of different extenders such as rice flour (10 %), barnyard millet flour (10 %) and TSGP (5 %) separately replacing lean meat in the formulation can give acceptable quality of the product. Rice flour was found to be the best among the three extenders studied as per the sensory evaluation.

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.

Chemical composition of solar dried blood and the ruminal content and its effect on performance of Japanese quails

Aim:

The aim was to determine the chemical composition of solar dried blood and rumen content (DBRC) and further ascertain the concentration at which DBRC could be included in Japanese quail diets without any adverse effect on its performance.

Materials and Methods:

Feeding trial on the effect of DBRC on performance of Japanese quails was studied up to 5 weeks. 252 numbers of day old (Nandanam Type III breed) Japanese quails were purchased from Poultry Research Station, Madhavaram and divided into 7 batches (control+ six treatments) each consisting of 36 birds. The DBRC was included at 0%, 5%, 10%, 15%, 20%, 25% and 30% in diets as control, treatment-1 (T1), treatment-2 (T2), treatment-3 (T3), treatment-4 (T4), treatment-5 (T5) and treatment-6 (T6) respectively in a completely randomized design to replace soybean meal in Japanese quail feed. The birds were provided with ad-labidum feed and drinking water ad-libitum during the entire experimental period.

Results:

The crude protein (CP), crude fiber (CF), ether extract (EE) and ash contents of DBRC were 35.87%, 17.40%, 3.6% and 12.6%, respectively. The amount of essential amino acids and non-essential amino acid content were found to be 12.98 and 4.87 (g/100 g of feed) respectively in DBRC feed. Result showed that all birds fed DBRC diets performed better than the control group. Mortality was unaffected by dietary treatments. There was a significant difference (p<0.01) observed in weight gain in treatment groups compared to the control.

Conclusion:

Up to 30% DBRC could be incorporated in the diets of Japanese quails without any adverse effects on its performance.

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.

Development of a SYBR Green I based duplex real-time PCR for detection of bovine herpesvirus-1 in semen

Bovine herpesvirus-1 (BoHV-1) is a viral pathogen found in infected bull semen, which is transmitted to inseminated cows by artificial insemination. BoHV-1 infection can cause reproductive disorders leading to significant economic loss to cattle industry. To detect BoHV-1 in semen, in this study, a SYBR Green I based duplex real-time PCR was developed. The assay included primers from BoHV-1 glycoprotein C (gC) and bovine growth hormone (bGH) genes for simultaneous detection in single tube. The result was interpreted by analysing melting temperature (Tm) peaks obtained after melt curve analysis of the amplified products at the end of reaction. The Tm peaks for BoHV-1-gC indicated presence of BoHV-1 while the bGH peak indicated reaction without inhibition. The sensitivity of the assay was to detect ten BoHV-1 genome copies per reaction. The analytical sensitivity was to detect 0.21 TCID50 infectious BoHV-1 in spiked semen. The assay was validated with 80 semen samples collected from breeding bulls. The diagnostic sensitivity and specificity of the assay was 100% with OIE recommended TaqMan probe based real-time PCR.

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.

Sequence analysis and identification of new variations in the 5'-flanking region of αS2-casein gene in Indian zebu cattle

Regulatory region of milk protein alpha S2-casein (αS2-CN) gene sequence was characterized and analyzed for nucleotide variations in animals representing 13 Indian zebu cattle (Bos indicus) breeds. A total of 15 variations; 11 in promoter region (1.56 Kb): -1481 (C>T), -1412 (C>T), -1342 (C>T), -1084 (G>A), -979 (A>G), -657 (A>T), -508 (A>G), -186 (T>C), -184 (T>C), -151 (T>C) and -135 (C>T); 1 in 5'-UTR (44 bp): 7 (C>T) while, 3 in intron-I region (73 bp): 186 (C>T), 194 (A>C) and 301 (A>T) were identified. Additionally, single deletion was observed at -975 (A>-) but not involve any known potential transcription factor binding sites (TFBS). Comparison with Bos taurus sequence revealed two additional variations -1085 (T>C) and -739 (A>G). Out of the total 18 variations observed between indicine and taurine αS2-CN regulatory region sequence, 15 were novel to B. indicus and are reported for the first time. Among these, four variations were located within the potential TFBSs; -1342 (C>T) within HNF-3beta, -739 (A>G) within C/EBP-alpha while -657 (A>T) and -508 (A>G) were found within glucocorticoid receptor TFBSs. Variations located within or in proximity to putative TFBSs could possibly influence the binding affinity of nuclear factors towards DNA binding domains, thus affecting transcriptional rate of αS2-CN gene. Phylogenetically, as expected, Indian zebu cattle were grouped close to B. taurus and were most distantly placed in comparison to human. The study indicated possible genetic variations in the regulatory regions of αS2-CN gene within Indian native cattle (B. indicus) and also its comparison with evolutionary different B. taurus breeds.

Single nucleotide polymorphisms in caprine calpastatin gene

The calpains and calpastatin (CAST) make up a major cytosolic proteolytic system, the calpain-calpastatin system, found in mammalian tissues. The relative levels of the components of the calpain-calpastatin system determine the extent of meat tenderization during postmortem storage. Calpastatin (CAST) is a protein inhibitor of the ubiquitous calcium-dependent proteases-micro-calpain and m-calpain. Polymorphisms in the bovine, ovine and pig CAST gene have been associated with meat tenderness but little is known about how caprine CAST gene may affect goat meat quality traits. In this study we selected different parts of the CAST gene: 1) that have been previously reported to be polymorphic, intron 5 and 12 and 3'UTR; 2) first time explored (exon 3, 7 and 8 and part of intron 7 and 8) to investigate polymorphic status of caprine CAST gene. Using comparative sequencing ten novel SN Ps located in exon 3 and intron 5, 7 and 8 were identified. Previously reported SNPs in intron 5, 3'UTR and intron 12 were absent. Sequence analysis revealed a non synonymous amino acid variation in exon 3, which would result in Lys/Arg substitution in the corresponding protein sequence. Considerable variation was detected in intronic regions. Twenty-four InDel were also recognized in intronic regions (15) and 3'UTR (9). All the sequences shared high homology with published bovine and ovine sequences. Three PCR-RFLP loci have been established for further analyzing genetic polymorphism in indigenous goats.

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.

Selection of suitable reference genes for quantitative gene expression studies in milk somatic cells of lactating cows (Bos indicus)

We assessed the suitability of 9 internal control genes (ICG) in milk somatic cells of lactating cows to find suitable reference genes for use in quantitative PCR (qPCR). Eighteen multiparous lactating Sahiwal cows were used, 6 in each of 3 lactation stages: early (25 ± 5 d in milk), mid (160 ± 15 d in milk), and late (275 ± 25 d in milk) lactation. Nine candidate reference genes [glyceraldehyde 3-phosphate dehydrogenase (GAPDH), protein phosphatase 1 regulatory subunit 11 (PPP1R11), β-actin (ACTB), β-2 microglobulin (B2M), 40S ribosomal protein S15a (RPS15A), ubiquitously expressed transcript (UXT), mitochondrial GTPase 1 (MTG1), 18S rRNA (RN18S1), and ubiquitin (UBC)] were evaluated. Three genes, β-casein (CSN2), lactoferrin (LTF), and cathelicidin (CAMP) were chosen as target genes. Very high amplification was observed in 7 ICG and very low level amplification was observed in 2 ICG (UXT and MTG1). Thus, UXT and MTG1 were excluded from further analysis. The qPCR data were analyzed by 2 software packages, geNorm and NormFinder, to determine suitable reference genes, based on their stability and expression. Overall, PPP1R11, ACTB, UBC, and GAPDH were stably expressed among all candidate reference genes. Therefore, these genes could be used as ICG for normalization of qPCR data in milk somatic cells through lactation.

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.

Differentiation of four Indian aboriginal cattle populations revealed by STR markers

Cattle are the most important livestock in India and play a pivotal role in agrarian economy. There are 34 recognized breeds of cattle and number of unexplored lesser known populations. The present study is a contribution towards determining genetic variation and understanding the relationship among four lesser known populations. A total of 194 unrelated DNA samples from three cattle populations of Orissa (Binjharpuri, Ghumsuri, Motu) and Hill cattle of Kumaun (Kumauni) were collected from respective breeding tracts. Genotyping was done with 23 bovine microsatellite markers as suggested by International Society for Animal Genetics (ISAG) and FAO (DAD-IS) on automated sequencer. The average observed heterozygosity in the four populations lie within the narrow range of 0.623 +/- 0.04 in Binjharpuri to 0.664 +/- 0.03 in Kumauni. Mean estimates of observed and expected heterozygosity over all loci and breeds were 0.651 +/- 0.02 and 0.720 +/- 0.01, respectively. In the overall population, the homozygote excess (F(IT)) of 0.132 +/- 0.03, was partly due to the genetic differentiation among breeds (F(ST) = 0.044 +/- 0.01) and, to a larger extent, to a significant homozygote excess within breeds (F(IS) = 0.094 +/- 0.03). The phylogenetic reconstruction from a UPGMA clustering based on Nei's Standard genetic distance yielded a tree with Binjharpuri and Ghumsuri on a single node and Motu and Kumauni on separate nodes. The most probable clustering detected by STRUCTURE in population was three. Binjharpuri and Ghumsuri animals were assigned to one cluster with high proportion of membership.

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.

Microsatellite analysis of genetic population structure of zebu cattle (Bos indicus) breeds from north-western region of India

The present study aims to understand the existing genetic diversity and structure of six native cattle breeds (Rathi, Tharparkar, Nagori, Mewati, Gir, and Kankrej) adapted to the north-western arid and semi-arid region of India based on microsatellite loci. Various diversity estimates, mean number of alleles (12.84); effective number of alleles (5.02); gene diversity (0.769), and observed heterozygosity (0.667) reflected the existence of substantial within-breed diversity in all the investigated cattle breeds. Mean estimates of F-statistics: F(IT) = 0.144 ± 0.023, F(IS) = 0.071 ± 0.021, and F(ST) = 0.078 ± 0.014 were significantly different from zero (P < 0.05). The interbreed relationships indicated moderate level of breed differentiation between the six cattle breeds with least differentiation between Kankrej-Mewati pair. The phylogeny structuring further supported close grouping of Kankrej and Mewati breeds. Correspondence analysis plotted Rathi, Tharparkar, and Gir individuals into three separate areas of multivariate space; whereas, Kankrej, Mewati, and Nagori cattle showed low breed specific clustering. This reflected the existence of discrete genetic structure for Tharparkar, Rathi, and Gir, the prominent dairy breeds of the region; whereas, admixture was observed for Kankrej, Mewati, and Nagori individuals.