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.

Selection of reference genes for normalizing gene expression data across seasons in spermatozoa of water buffalo (Bubalus bubalis).

Selection of the most stably expressed reference genes is key to monitoring accurate target gene expression across any tissue or cell type. The mRNA in spermatozoa stores valuable information related to changes in spermatogenesis due to variations in environmental conditions, especially during heat stress, which affects various sperm functions. Semen quality in buffalo bulls is significantly influenced by the seasons. In the study, a panel of nine genes was evaluated to identify the most stably expressed internal control gene (ICG) for the normalization of real-time gene expression data generated across various seasons for Murrah buffalo bulls' spermatozoa. Sperm cells were purified from the semen samples collected during different seasons, with temperature-humidity index (THI) ranging from 80.80 ± 1.47 (hot summer) to 55.88 ± 1.98 (winter), using the BoviPure™ gradient purification method. The RNA isolated from the purified spermatozoa fraction was quality checked prior to reverse transcription and subjected to qPCR (quantitative real-time PCR) based expression analysis. An automated 'endoGene' pipeline was employed to apply the geNorm, NormFinder, and BestKeeper algorithms for data analysis. The result indicated that GAPDH and PP1A were the most stably expressed among the gene panel, whereas ATPSF1 and ACTB were the two least stable expressed reference genes. Further, the most suitable ICGs identified were validated by normalization of real time expression data of heat stress and sperm quality genes, HSFY2 and AKAP4, respectively. The genes identified would help in generating the most reliable results for the expression profiling of the genes dictating sperm quality and heat stress cope-up mechanism in buffalo spermatozoa, collected during different seasons.

Development and validation of most efficient RNA isolation method from buffalo bull spermatozoa.

BACKGROUND: Being highly fragmented and low in concentration, isolation of good quality RNA from sperm cells is a big challenge. Attempts have been made to evaluate various sperm RNA isolation methods from purified buffalo bull sperm cells. METHODS: Both, non-membrane and membrane-based methods have been evaluated for isolating RNA from Murrah buffalo sperms and compared for their respective efficacies. The traditional TRIzol, TRIzol-heat lysed (H-TRIzol) and cocktail of TCEP-RLT lysis buffer (Qiagen RNeasy mini kit)-TRIzol (C-TRIzol) based isopropanol isolation methods have been evaluated. RESULTS: H-TRIzol yielded best results among conventional methods. The combined T-RLT RNA isolation protocol yielded best quality and quantity compared to other membrane-based methods, due to high lytic property of cocktail of lysis reagents, necessary for complete breakdown of sperm membrane and RNA binding membrane for RNA isolation. Combined lysis performed by treatment with RLT-T and T-RLT differing in order of reagents used were also evaluated. T-RLT combination giving better results compared to RLT-T due to high gDNA contamination and membrane clogging in later protocol steps. CONCLUSION: Overall, in terms of total RNA quantity and quality per million spermatozoa, the heat-lysed TRIzol method (H-TRIzol) performs best among RNA separation techniques employed and is also quite easy to perform. This comparative evaluation of sperm RNA isolation protocols can be useful in deciding the best protocol for isolation of good quality and high concentration sperm RNA from buffalo semen, for transcriptome and other downstream studies.

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.

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.

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.

Current status and unique attributes of Indian Chilika buffalo for adaptation to brackish water ecology.

Chilika buffalo is native to the Eastern coast of India and well adapted to the largest coastal brackish water lagoon of Asia, Chilika Lake. We present here a report on the Chilika buffalo breed emphasizing the conservational urgency based on unique biochemical and molecular evidence related to liver and kidney functions while comparing it with tropically adapted other water buffalo breeds (Bubalus bubalis) of India. It is found that the Chilika buffalo breed has a better ability to withstand a long dehydration period as evident from its better glomerular filtration and higher expression of the ion transport channel. Mitochondrial D-loop sequencing results have shown these buffaloes being closer to swamp-type buffaloes of Bangladesh and northeast India and represent a unique "hybrid zone" on the eastern coast of India. Conservation of such uniquely adapted germplasm is crucial owing to the current global trend, where the introduction of exotic breeds has negatively impact "sui-generis" germplasm and they require higher managerial resource consumption for maintaining higher productivity. Further, the introduction of unconventional fisheries activities has proved detrimental to the lagoon ecosystem, potentially causing more threat to the buffalo's population.

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.

RNA sequencing and transcriptome analysis of buffalo (Bubalus bubalis) blastocysts produced by somatic cell nuclear transfer and in vitro fertilization.

Across farm animal species, the live birth rate obtained with somatic cell nuclear transfer (SCNT) embryos is only <2% compared with >40% obtained with in vitro fertilization (IVF) embryos, primarily due to incomplete nuclear reprogramming which results in aberrant embryonic gene expression. We used RNA sequencing to compare the global transcriptome profile of SCNT and IVF buffalo blastocysts. SCNT blastocysts expressed 17,061 transcripts, of which 941 were unique whereas, IVF blastocysts expressed 17,303 transcripts, of which 1,183 were unique. At ≥2-folds change (p < .05), 331 transcripts were differentially expressed in the two groups among which, 19 were unique, 188 were downregulated and 143 were upregulated in SCNT compared with IVF blastocysts. Many genes affecting pluripotency, trophectoderm development, developmental regulation, and epigenetic modifications were upregulated in SCNT compared with IVF blastocysts. Among the four functional categories analyzed, epigenetic regulators were the most affected. Most of the WNT signaling pathway genes were upregulated whereas, the inhibitors of this pathway, such as DKK1, were downregulated in SCNT blastocysts, suggesting that this pathway is overexpressed in SCNT embryos. Gene Ontology analysis revealed that 25 biological processes, 20 molecular functions, and 24 cellular compartment categories were enriched in SCNT blastocysts. This data can help identify reprogramming errors for improving cloning efficiency.

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.

Buffalo SCNT embryos exhibit abnormal gene expression of ERK/MAPK pathway and DNA methylation.

Inhibition of ERK/MAPK pathway has been shown to decrease DNA methylation via down-regulation of DNA methyltransferases (DNMTs) in several studies suggesting that this pathway plays an important role in regulation of DNA methylation. We examined the relative expression level of seven important genes related to ERK/MAPK pathway and DNMTs (DNMT1, DNMT3a and DNMT3b) by quantitative real-time PCR in buffalo blastocysts produced by Hand-made cloning and compared it with that in blastocyst-stage embryos produced by in vitro fertilization (IVF). The expression level of six of seven genes related to ERK/MAPK pathway examined i.e., p21RAS, RAF1, AKT1, ERK2, PIK3R2 and c-Myc was significantly higher (p < 0.05) in cloned than in IVF embryos. However, the expression level of FOS was lower (p < 0.005) in cloned than in IVF embryos. The relative expression level of DNMT3a and DNMT3b but not that of DNMT1 was significantly higher (p < 0.05) in cloned than in IVF embryos. These results indicate that the cloned embryos exhibit an abnormal expression of several important genes related to ERK/MAPK pathway and DNMTs. Although a direct link between ERK/MAPK pathway and DNMTs was not examined in the present study, it can be speculated that ERK/MAPK pathway may have a role in regulating the expression of DNMTs in embryos, as also observed in other tissues.

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.

Identification of polymorphism in fatty acid binding protein 3 (FABP3) gene and its association with milk fat traits in riverine buffalo (Bubalus bubalis).

The fatty acid binding protein 3 (FABP3) gene, known to be associated with fat percentage of milk and meat in bovines, was screened among swamp and riverine buffaloes for polymorphism detection and further association with milk fat contents. An SNP g.307C > T was identified in the intron 2 (+53 exon 2) region of FABP3 gene of Indian buffaloes. The SNP identified was genotyped in 692 animals belonging to 15 riverine, swamp and hybrid (riverine × swamp) buffalo populations of diverse phenotypes and utilities, by PCR-RFLP. A marked contrast was observed between the C and T allele frequencies in three types of buffaloes. The frequency of C allele ranged from 0.67 to 0.96 in pure swamp buffalo populations, with the highest in Mizoram (0.96). Whereas the frequency of T allele was high across all the Indian riverine buffalo breeds, ranging from 0.57 to 0.96. None of the genotypes at FABP3 g.307C > T locus was found to have significant association with milk fat and other production traits in Mehsana dairy buffalo breed. Our study revealed marked differences in the allele frequencies between riverine and swamp buffaloes at FABP3 g.307C > T locus, without any significant association with different milk traits in riverine buffaloes.

Genetic diversity and relationship of Indian cattle inferred from microsatellite and mitochondrial DNA markers.

BACKGROUND: Indian agriculture is an economic symbiosis of crop and livestock production with cattle as the foundation. Sadly, the population of indigenous cattle (Bos indicus) is declining (8.94% in last decade) and needs immediate scientific management. Genetic characterization is the first step in the development of proper management strategies for preserving genetic diversity and preventing undesirable loss of alleles. Thus, in this study we investigated genetic diversity and relationship among eleven Indian cattle breeds using 21 microsatellite markers and mitochondrial D loop sequence. RESULTS: The analysis of autosomal DNA was performed on 508 cattle which exhibited sufficient genetic diversity across all the breeds. Estimates of mean allele number and observed heterozygosity across all loci and population were 8.784 ± 0.25 and 0.653 ± 0.014, respectively. Differences among breeds accounted for 13.3% of total genetic variability. Despite high genetic diversity, significant inbreeding was also observed within eight populations. Genetic distances and cluster analysis showed a close relationship between breeds according to proximity in geographic distribution. The genetic distance, STRUCTURE and Principal Coordinate Analysis concluded that the Southern Indian Ongole cattle are the most distinct among the investigated cattle populations. Sequencing of hypervariable mitochondrial DNA region on a subset of 170 cattle revealed sixty haplotypes with haplotypic diversity of 0.90240, nucleotide diversity of 0.02688 and average number of nucleotide differences as 6.07407. Two major star clusters for haplotypes indicated population expansion for Indian cattle. CONCLUSIONS: Nuclear and mitochondrial genomes show a similar pattern of genetic variability and genetic differentiation. Various analyses concluded that the Southern breed 'Ongole' was distinct from breeds of Northern/ Central India. Overall these results provide basic information about genetic diversity and structure of Indian cattle which should have implications for management and conservation of indicine cattle diversity.

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.

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.