Hypoxia related genes modulate in similar fashion in skin fibroblast cells of yak (Bos grunniens) adapted to high altitude and native cows (Bos indicus) adapted to tropical climate during hypoxia stress.

The present study was conducted to understand transcriptional response of skin fibroblast of yak (Bos grunniens) and cows of Bos indicus origin to hypoxia stress. Six primary fibroblast cell lines derived from three individuals each of Ladakhi yak (Bos grunniens) and Sahiwal cows (Bos indicus) were exposed to low oxygen concentration for a period of 24 h, 48 h and 72 h. The expression of 10 important genes known to regulate hypoxia response such as HIF1A, VEGFA, EPAS1, ATP1A1, GLUT1, HMOX1, ECE1, TNF-A, GPx and SOD were evaluated in fibroblast cells of Ladakhi yak (LAY-Fb) and Sahiwal cows (SAC-Fb) during pre- and post-hypoxia stress. A panel of 10 reference genes (GAPDH, RPL4, EEF1A1, RPS9, HPRT1, UXT, RPS23, B2M, RPS15, ACTB) were also evaluated for their expression stability to perform accurate normalization. The expression of HIF1A was significantly (p < 0.05) induced in both LAY-Fb (2.29-fold) and SAC-Fb (2.07-fold) after 24 h of hypoxia stress. The angiogenic (VEGFA), metabolic (GLUT1) and antioxidant genes (SOD and GPx) were also induced after 24 h of hypoxia stress. However, EPAS1 and ATP1A1 induced significantly (p < 0.05) after 48 h whereas, ECE1 expression induced significantly (p < 0.05) at 72 h after exposure to hypoxia. The TNF-alpha which is a pro-inflammatory gene induced significantly (p < 0.05) at 24 h in SAC-Fb and at 72 h in LAY-Fb. The induction of hypoxia associated genes indicated the utility of skin derived fibroblast as cellular model to evaluate transcriptome signatures post hypoxia stress in populations adapted to diverse altitudes.

Editing of HSF-1 and Na/K-ATPase α1 subunit by CRISPR/Cas9 reduces thermal tolerance of bovine skin fibroblasts to heat shock in vitro.

A follow-up to our previous findings, the present study was planned to evaluate the role of Na/K-ATPase alpha1-subunit (ATP1A1) gene in heat shock tolerance. The primary fibroblast culture was established using ear pinna tissue samples of Sahiwal cattle (Bos indicus). The knockout cell lines of Na/K-ATP1A1 and HSF-1 (heat shock factor-1, as a positive control) genes were developed by CRISPR/Cas9 method and the gene-editing was confirmed by the genomic cleavage detection assay. The two knockout cell lines (ATP1A1 and HSF-1) and wild-type fibroblasts were exposed to heat shock at 42 °C in vitro and different cellular parameters viz., apoptosis, proliferation, mitochondrial membrane potential (ΔΨm), oxidative stress, along with expression pattern of heat-responsive genes were studied. The results showed that in vitro heat shock given to knockout fibroblast cells of both ATP1A1 and HSF-1 genes resulted in decreased cell viability, while increasing the apoptosis rate, membrane depolarization, and ROS levels. However, the overall impact was more in HSF-1 knockout cells as compared to ATP1A1 knockout cells. Taken together, these results indicated that the ATP1A1 gene plays a critical role as HSF-1 under heat stress and helps cells to cope with heat shock.

Allelic diversity at BoLA DRB3 locus and association with predisposition to clinical mastitis in indicus and crossbred cattle.

Bovine lymphocyte antigen (BoLA) DRB3 locus in healthy and mastitis affected cattle has been genotyped by a polymerase chain reaction and restriction fragment length polymorphisms (PCR-RLFP) using RsaI restriction enzyme, followed by sequencing. In 130 farm animals, 25 BoLA DRB3 alleles have been detected by PCR-RFLP. Three distinct allelic patterns significantly associated with mastitis in Karan Fries crossbred and Sahiwal indicus cattle have been identified, whereas, four other allelic patterns were significantly high in frequency among healthy animals. Sequencing of RFLP genotypes revealed 25 and 47 alleles among healthy Sahiwal and Karan Fries, respectively, while 17 and 38 patterns observed in mastitis affected Sahiwal and Karan Fries animals, respectively. From Tajima's D-test of neutrality, it was concluded that alleles associated with mastitis were expanding in the population, whereas those of healthy were under contraction. Phylogenetic analysis carried out to delineate the evolutionary relationship of the farm and field animals at DRB3 locus, differentiating allelic patterns into six different clusters. Among the phylogenetic lineages, five patterns DRB3*028:01, DRB3*011:03, DRB3*031:01, DRB3*001:01 and DRB3*043:01, were previously reported, whereas one novel allelic variant was observed in indicus and crossbred cattle. This information will help in further exploring the association between BoLA-DRB3 genetic diversity and disease resistance in distinct cattle breeds, important in designing breeding strategies for increasing the distribution of favorable alleles.

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

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

Genetic admixture and population structure analysis of Indian water buffaloes (Bubalus bubalis) using STR markers.

BACKGROUND: India has a vast riverine and swamp buffalo diversity adapted to various agro-ecological conditions. In the present study, genetic diversity data for 10 different buffalo populations of India, using 20 highly polymorphic microsatellite markers has been generated for the genetic diversity analysis. The buffalo populations of Eastern Odisha state, were the primary focus. METHODS AND RESULTS: The minimal spanning network based on Bruvo's distance, PCA (Principal Component Analysis) based on the Fst (Fixation Index) values, and genetic admixture analysis using both the STRUCTURE and 'snapclust' were performed. The analysis could identify the Manda population as distinct from other Odisha buffalo breeds as well as adjoining Chhattisgarhi buffalo breeds. The total observed number of alleles ranged between 143 (Manda) and 301 (Paralakhemundi) with an average of 204 alleles per breed. The Sambhalpuri buffalo population also clustered into two separate subpopulations, half of the unique sub-population located geographically south-wards, displayed no admixture with any of the adjacent buffalo populations. The Manda buffalo population has shown sufficient allelic richness and heterozygosity under random mating being practiced in the field conditions. CONCLUSIONS: The study has led to the identification of the Manda as a distinct buffalo population, and the germplasm has been registered as a new Indian buffalo breed. Whereas, the Sambhalpuri population requires elaborate analysis to confirm the existence of two distinct sub-populations.

Genetic characterization and diversity assessment in 'Bhangor' indigenous swamp buffalo population using heterologous microsatellite markers.

'Bhangor' newly identified swamp buffalo population from North East Indian, was characterized using microsatellite markers. Genomic DNA was isolated from blood samples of 76 unrelated animals, 15 microsatellite markers (CSSM33, BM1818, CSRM60, HEL13, ILSTS019, ILSTS025, ILSTS028, ILSTS029, ILSTS033, ILSTS036, ILSTS056, ILSTS058, ILSTS061, ILSTS089 and ETH003) were found to be highly polymorphic in the population of the selected markers. A total of 114 alleles were observed, which ranged from 3 in CSRM60 and ILSTS025 locus to 12 in ILSTS056 and ILSTS061. The mean effective number of alleles across all polymorphic loci was found to be 3.76. The overall mean expected heterozygosity and unbiased expected heterozygosity values were 0.67 and 0.68, ranging from 0.067 (ILSTS025) to 0.85 (ILSTS058) and 0.068 (ILSTS025) to 0.86 (ILSTS058), respectively. Within the population, the inbreeding estimates (FIS) ranged between -0.4352 and 0.804, with an average FIS of 0.114 ± 0.033. The outcome for infinite allele model (IAM), two-phase model (TPM) and test for mode shift revealed the absence of any recent bottleneck in the investigated buffalo population. The population was found to be in optimum diversity based on polymorphic microsatellite markers. With fast changing agro-climatic conditions; there is an urgent need to characterize the nondescript livestock populations.

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

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

Characterization of thermo-physiological, hematological, and molecular changes in response to seasonal variations in two tropically adapted native cattle breeds of Bos indicus lineage in hot arid ambience of Thar Desert.

The selection of climate resilient animal is necessary to secure the future of sustainable animal production. The present investigation therefore was an effort to unravel answers to the adaptation at physiological, hematological, and molecular levels in cows of hot arid region that helps them to survive harsh environment, to continue production and reproduction. This investigation was carried out in indicine cows over a period of one year, encompassing four seasons, wherein physiological data of 50 animals, hematological data of 15 animals, and gene expression profile of 5 animals from each of Sahiwal and Kankrej breeds per season was generated. In total, 5600 physiological observations, 1344 hematological observations, and 480 molecular samples were processed. The meteorological data revealed a high diurnal variation of temperature across seasons, with THI exceeding 80 during the months of summer and hot-humid seasons, indicating significant heat stress (HS). The physiological parameters showed an increasing trend with the incremental THI, with significantly (p < 0.05) higher values of rectal temperature (RT), respiration rate (RR), pulse rate (PR), and body surface temperature (BST) at ventral (VT), lateral (LT), dorsal (DT), and frontal (FT), in both breeds recorded during HS. The hematological pictures also revealed significant (p < 0.05) seasonal perturbations in erythrocytic and leucocytic parameters. Moreover, the molecular response was driven by a significant (p < 0.05) upregulation of all the key HSPs, HSP70, HSP90, HSP60, and HSP40, except HSP27 during the hotter months of summer and hot-humid seasons. The expression of HSF1, an important transcriptional regulator of  HSP70 was also significantly (p < 0.05) upregulated during summer season in both breeds. All the molecular chaperones revealed a significant upregulation during the summer season, followed by a decreasing trend by hot-humid season. The study indicated a well-developed thermotolerance mechanism in animals of both breeds, with Kankrej cows exhibiting better thermotolerance compared to Sahiwal cows.

In Silico Analysis of HSP70 Gene Family in Bovine Genome.

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

Casein (CSN) gene variants and parity affect the milk protein traits in crossbred (Bos taurus x Bos indicus) cows in sub-tropical climate.

Genotypes at four casein (CSN) loci-A26181G of CSN1S1, C6227T of CSN1S2, A8101C of CSN2, and A13104C of CSN3-along with non-genetic factors were studied for their effects on various milk protein traits in 100 crossbred cows with major inheritance of Holstein Frisian (Bos taurus) and Tharparkar (Bos indicus). Results demonstrated the presence of all CSN genotypes with a predominance of heterozygotes. At CSN2 (A8101C; His67Pro) locus, the A2 allele, desirable for human health, was present in 62% as heterozygous and 29% in homozygous condition. Among non-genetic factors, parity of the cows had a significant influence on the milk protein traits in these crossbreds. The genotypes at CSN1S1, CSN2, and CSN3 loci were found to influence (p<0.05 to 0.01) the casein and whey protein yields and composition traits. The casein index and total milk yield were most influenced by the CSN1S2 locus. The AA (A1 milk) genotype of CSN2 had significantly higher yields and percentages of casein and whey proteins. Positive influence of CC genotype of CSNS3 on milk proteins of was observed similar to Bos taurus cows; however, such influence of AA genotype of CSN2 locus may be distinctive to the crossbred cows maintained in subtropical condition. Overall, the results revealed the diverse effects of CSN genotypes on milk proteins in crossbred cattle.

Genetic characterization and population structure of different coat colour variants of Badri cattle.

The present study aimed to genetically characterize the Badri cattle and its three colour variants and assess their population structure using 24 microsatellite markers. Out of 96 animals analyzed, 32 each were collected from grey (GVBC), brown (BrVBC) and black (BVBC) colour variants of Badri cattle (BC). The genetic diversity parameters including allele frequencies, observed and effective number of alleles, observed and expected heterozygosity, PIC, Shannon's indices and F-statistics were estimated using POPGENE software. Bottleneck analysis was performed using both qualitative and quantitative approaches. A total of 274 alleles (50 private and 224 shared) were scored for BC, GVBC, BrVBC and BVBC with mean number of 11.417, 9.083, 9.125 and 9.083 alleles, respectively. All populations exhibited average heterozygosity estimate > 0.5 indicating existence of substantial genetic variability, concurrent with revelations from Shannon's indices. Observed mean PIC estimates (> 0.74) were indicative of optimum informativeness of used microsatellite markers. The mean inbreeding estimates (F) in GVBC, BrVBC and BVBC were 0.041, - 0.024 and 0.016, respectively. The pair wise genetic (> 0.91) pointed towards similarity between different colour variant populations. STRUCTURE analysis also revealed clear admixture for the three Badri colour variants indicating absence of genetic differentiation. The present study revealed first-hand information that populations of Badri cattle with different phenotypes with respect to coat colour are genetically related and can be considered as a single breed. The comprehensive knowledge generated for Badri cattle will help in designing breeding plan for its genetic improvement and deciding the conservation priorities.

Characterizing binding sites of heat responsive microRNAs and their expression pattern in heat stressed PBMCs of native cattle, exotic cattle and riverine buffaloes.

It is generally believed that due to evolutionary differences and adaptation to tropical conditions, Indian native cattle has superior heat tolerant ability than Bos taurus cattle. In the present study, 3'-UTR of two most important heat responsive genes i.e., heat shock protein 70.1 (HSP70.1) and heat shock factor- 1 (HSF-1) were sequence characterized in different breeds of Indian native cattle to identify the variations and miRNA binding sites. In addition, the impact of heat stress was assessed in a total of 57 PBMCs samples of native Sahiwal cows (Bos indicus), exotic Holstein cows (Bos taurus) and Murrah buffaloes (Bubalus bubalis) using various cellular parameters like cell viability, cytotoxicity and apoptosis. Further, expression profile of 12 heat responsive miRNAs were also evaluated in unstressed and stressed PBMCs to understand post transcriptional changes in native cows, exotic cows and Murrah buffaloes. The sequence data showed 3'-UTR of HSP70.1 gene of Indian cattle to be exactly similar to Bos taurus with no miRNA binding site. Whereas, sequencing of 3'-UTR of HSF-1 gene revealed 3 SNPs at positions G1762T; C1811T and C1983T with 7 well conserved miRNA binding sites. The impact of heat stress on various cellular parameters in terms of cell viability, cytotoxicity and apoptosis was highest in PBMCs of Holstein cows followed by Murrah buffaloes and Sahiwal cows. Further, in contrast to Holstein Frisian cows and Murrah buffaloes, the expression pattern of 12 heat responsive miRNAs, in heat stressed PBMCs of Sahiwal cows were quite distinct. There was a significant (p < 0.05) induction in expression of most of the miRNAs after heat stress in PBMCs of Sahiwal cows followed by a rapid decline. The distinct cellular response and pattern of miRNA expression across cattle types and buffaloes might be influencing their PBMCs tolerance level to heat stress.

Evaluation of Milk Colostrum Derived Lactoferrin of Sahiwal (Bos indicus) and Karan Fries (Cross-Bred) Cows for Its Anti-Cancerous Potential.

Lactoferrin (Lf) is an iron-binding glycoprotein protein known to have immune-modulatory role and recently, its anticancerous effect against different cancer cell types was emphasized. In the present investigation, a comparative evaluation of anticancer potential of colostrum-derived lactoferrin from Indian native zebu cow (Sahiwal, SAC), crossbred (Karan Fries, KFC) and commercially available (C-Lf) lactoferrin from exotic cow using cellular models was made. A protocol was standardized successfully to purify Lf protein from colostrum of both breeds using HPLC and purity was confirmed by LC-MS. A standardized dose of 750 µg/mL Lf was used to treat two cell types MDA-MB-231 and MCF-7 with Lf from three different sources; SAC-Lf, KFC-Lf and C-Lf for 48 h and 72 h. Different cellular parameters including cytotoxicity, viability, apoptosis and cell proliferation were determined. Comparatively, Lf from commercial source (C-Lf) had maximum effect in both cell types followed by SAC-Lf and KFC-Lf. Further, transcriptional changes in genes associated with apoptosis (Bax and Bcl-2), tumor progression (p53, p21, CD44 and NF-κß) and survival (survivin) were evaluated in Lf treatment. The overall results strongly emphasized to the fact that Lf purified from cow colostrum has the capacity to inhibit the in vitro growth of cancerous cell lines albeit to a varied extent.

Overexpression of genes associated with hypoxia in cattle adapted to Trans Himalayan region of Ladakh.

Ladakh is an important part of the Trans-Himalayan region located between the Kunlun mountain range in the north and the main Great Himalayas to the south in the state of Jammu and Kashmir of India. The local cattle from Leh and Ladakh region, known as "Ladakhi cattle" is a unique germplasm having an excellent adaptation potential to high altitude hypobaric stress. In the present study, an effort was made to evaluate the transcriptional pattern of hypoxia inducing factor-1 (HIF-1) and several of its regulated genes in PBMCs of local Ladakhi cattle, Holstein Frisian crosses, Jersey (exotic) maintained at high altitude region and Sahiwal (Bos indicus) and Karan Fries (cross bred) cattle maintained in tropical environment. The combined data set indicated increased expression of HIF-1 and its regulated genes viz., glucose transporter 1 (GLUT1), vascular endothelial growth factor (VEGF), and hexokinase (HK2) in high altitude cattle indicating their importance in maintaining cellular homeostasis during high altitude hypoxia. The data indicated that hypoxia associated genes accumulated under hypoxic conditions are part of an essential adaptive component for adaptation to the high altitude of the trans-Himalayan region. In contrary, higher expression of molecular chaperons' viz., HSP70 and HSP90 in tropically adapted cattle give tolerance to high ambient temperature prevalent in tropical condition. In conclusion, HIF-1 and its regulatory genes could be termed as important candidates for producing homeostatic responses to hypoxia in cattle populations reared in higher altitudes of the Trans-Himalayan region.

Matrix-based three-dimensional culture of buffalo mammary epithelial cells showed higher induction of genes related to milk protein and fatty acid metabolism.

Demanding transcriptomic studies in livestock animal species could be replaced by good in vitro models mimicking the function of mammary gland. Mammary epithelial cells (MEC) are the functional unit of the mammary gland. Extracellular matrix is known to be a key factor providing normal homeostasis in three-dimensional (3D) environment as important signals are lost when cells are cultured in two-dimensional (2D) environment. The aims of this study were to establish a buffalo mammary epithelial cells (BMECs) in 3D culture using extracellular matrix and to determine whether such a 3D culture model has different expression pattern than 2D counterpart. The purified MEC generated after several passages were used to establish 3D culture using Geltrex matrix. The expression of milk casein genes viz., alpha S1-casein (CSN1S1), alpha S2-casein (CSN1S2), beta-casein (CSN2), kappa-casein (CSN3); and fatty acid metabolism genes viz., butyrophilin (BTN1A1), glycerol-3-phosphate acyltransferase (GPAM), fatty acid-binding protein 3 (FABP3), and stearoyl-CoA desaturase (SCD) was assessed in 3D culture in comparison to traditional monolayer culture using qRT-PCR. Notable morphological differences were observed for BMECs grown in 3D culture in comparison to 2D culture. Morphologically, epithelial structures grown in Geltrex matrix (3D) environment showed enhanced functional differentiation in comparison to 2D culture. In 3D culture, lumen and dome-like structures were formed by day 5, whereas polarized acinus-like structure were formed within 15 days of culturing. The expression data showed higher mRNA induction of milk casein and fatty acid metabolism genes in 10-day-old 3D BMECs culture in comparison to 2D monolayer culture. The result suggests that 3D organization of epithelial cells has favorable effect on induction of milk and fatty acid metabolism-related genes. Therefore, matrix-based 3D culture of MEC that recapitulate the structural and functional context of normal tissues could provide a better in vitro model to understand the mammary gland functioning of buffaloes.

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

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

Stage specific expression of ATP-binding cassette and solute carrier superfamily of transporter genes in mammary gland of riverine buffalo (Bubalus bubalis).

In the present study, expression level of various ATP-binding cassette (ABC) viz., ABCA1, ABCA7, ABCG1, ABCG2, and ABCG5; associated transcription factors viz., SREBF1, LXRα (NR1H3), PPARA, and Solute Carriers (SLC); or Glucose transporters (GLUT) viz., SLC2A1(GLUT1), SLC2A4 (GLUT4), SLC2A8 (GLUT8), and SLC2A12 (GLUT12) superfamily of transporters were compared across physiological stages of buffalo mammary gland. The relative expression of ABCA1, and ABCG1 was significantly (p < 0.05) higher in mammary gland of heifer followed by involution and lactation stages. Similarly, ABCA7 gene expression was highest in heifer mammary gland followed by lactation and involution stages. ABCG2 gene expression was significantly (p < 0.05) high in lactating mammary gland in comparison to involution and heifer stages. On the other hand, ABCG5 gene expression was highest in involuting mammary gland followed by lactation and involution stages. Additionally, the expression of LXRα SREBF1, and PPARA which are known to regulate some of the ABC tranporters were also analyzed. The expression of LXRα gene was high in involuting as compared to lactating mammary gland. In contrast, SREBF1 and PPARA expression was significantly (p < 0.05) high in lactating mammary gland. Among the several SLC transporters studied, SLC2A1, SLC2A4, and SLC2A8 showed significant (p < 0.05) higher expression during lactation stage, whereas SLC2A12 expression was greater during heifer stage suggesting SLC2A1, SLC2A4, and SLC2A8 to be the major transporters associated with glucose uptake in buffalo mammary gland. The expression profile of (lactoferrin) LTF, known to be expressed at high level in mammary gland during involution was also studied. As expected, its expression was significantly (p < 0.05) higher during involution in comparison to lactating mammary gland.in buffaloes as well. The inclusion of LTF as a control gene further provided the confidence in the buffalo mammary gland expression data generated in the present study. This study thus helped to provide information about the distinct expression pattern of various transporters and their regulators in buffalo mammary gland during different physiological states.

Selection of stable reference genes in heat stressed peripheral blood mononuclear cells of tropically adapted Indian cattle and buffaloes.

Eleven reference genes from different functional categories were assessed for their stable expression pattern in heat stressed PBMCs across Indian buffalo and cattle using GeNorm, NormFinder and BestKeeper algorithms. As the first report, we suggest B2M, RPS9 and RPS15a as suitable reference genes for accurate normalization of PBMC transcript data while the ACTB gene is not recommended.

ddRAD sequencing based genotyping of six indigenous dairy cattle breeds of India to infer existing genetic diversity and population structure.

The present investigation aimed to identify genome wide SNPs and to carry out diversity and population structure study using ddRAD-seq based genotyping of 58 individuals of six indigenous milch cattle breeds (Bos indicus) such as Sahiwal, Gir, Rathi, Tharparkar, Red Sindhi and Kankrej of India. A high percentage of reads (94.53%) were mapped to the Bos taurus (ARS-UCD1.2) reference genome assembly. Following filtration criteria, a total of 84,027 high quality SNPs were identified across the genome of 6 cattle breeds with the highest number of SNPs observed in Gir (34,743), followed by Red Sindhi (13,092), Kankrej (12,812), Sahiwal (8956), Tharparkar (7356) and Rathi (7068). Most of these SNPs were distributed in the intronic regions (53.87%) followed by intergenic regions (34.94%) while only 1.23% were located in the exonic regions. Together with analysis of nucleotide diversity (π = 0.373), Tajima's D (D value ranging from - 0.295 to 0.214), observed heterozygosity (HO ranging from 0.464 to 0.551), inbreeding coefficient (FIS ranging from - 0.253 to 0.0513) suggested for the presence of sufficient within breed diversity in the 6 major milch breeds of India. The phylogenetic based structuring, principal component and admixture analysis revealed genetic distinctness as well as purity of almost all of the 6 cattle breeds. Overall, our strategy has successfully identified thousands of high-quality genome wide SNPs that will further enrich the Bos indicus representation basic information about genetic diversity and structure of 6 major Indian milch cattle breeds which should have implications for better management and conservation of valuable indicine cattle diversity.

Identification of stably expressed Internal Control Genes (ICGs) for normalization of expression data in liver of C57BL/6 mice injected with beta casomorphins.

In recent years, beta-casomorphin peptides (BCM7/BCM9) derived from the digestion of cow milk have drawn a lot of attention world over because of their proposed impact on human health. In order to evaluate the transcriptional modulation of target genes through RT-qPCR in response to these peptides, availability of appropriate reference or internal control genes (ICGs) will be the key. The present study was planned to identify a panel of stable ICGs in the liver tissue of C57BL/6 mice injected with BCM7/BCM9 cow milk peptides for 3 weeks. A total of ten candidate genes were evaluated as potential ICGs by assessing their expression stability using software suites; geNorm, NormFinder and BestKeeper. The suitability of the identified ICGs was validated by assessing the relative expression levels of target genes, HP and Cu/Zn SOD. Based on geNorm, PPIA and SDHA gene pair was identified to be most stably expressed in liver tissue during the animal trials. Similarly, NormFinder analysis also identified PPIA as the most stable gene. BestKeeper analysis showed crossing point SD value for all the genes in the acceptable range that is closer to 1. Overall, the study identified a panel of stable ICGs for reliable normalization of target genes expression data in mice liver tissues during BCM7/9 peptides trial.