R. S. Kataria

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.

Identification of suitable housekeeping genes for normalization of quantitative real‐time PCR data during different physiological stages of mammary gland in riverine buffaloes (Bubalus bubalis)

Gene expression analysis unravels the complex changes or relations at transcriptomic level. To nullify all type of errors that can be incorporated during any stage of RNA extraction into cDNA synthesis and for reliable results, the data obtained from qPCR have to be normalized using the appropriate/suitable housekeeping genes (HKGs). Unfortunately, till date, no such HKG has been reported for bubaline mammary gland. The objective of the present study was thus to identify and validate the potential HKGs for the gene expression studies in buffalo mammary gland. Mammary tissues from twelve buffaloes during different physiological stages: pre-pubertal (heifer), lactation and involution were obtained for the present study. A total of 16 potential HKGs (GAPDH, β-actin, UXT, β2M, A2M, RPl4, RPS9, RPS15A, RPS18, RPS23, HMBS, HPRT1, GTP, EEF1A1, UB1 and RPL22) from different functional classes were evaluated. The analysis revealed that the expression of EEF1A1, RPl4, β2M and RPS15A was most consistent across different physiological stages of buffalo mammary gland. On the other hand, β-actin, A2M, RPL22 and GAPDH were the least stable genes making them unsuitable as HKGs. Based on our analysis, we recommend the use of EEF1A1, RPl4, β2M and RPS15A genes as suitable HKGs for accurate normalization of gene expression data in bubaline mammary gland.

Characterization of Osteopontin gene of Bubalus bubalis

Osteopontin, a glycoprotein, is expressed in several tissues including the mammary gland. The gene has been reported to be associated with milk and its constituents in various livestock species. This gene was sequenced in buffalo and it coded for the protein of 280 amino acids with the conserved GRGDS domain. The sequence was confirmed from the cDNA sequence derived from the mammary gland of buffalo. The earlier-reported 9T/10T variation in the upstream region of the gene was investigated for its effect on milk in buffalo and it was found to be non-significant.

Identification of internal control genes in milk-derived mammary epithelial cells during lactation cycle of Indian zebu cow

The present study aims to evaluate the suitability of 10 candidate genes, namely GAPDH, ACTB, RPS15A, RPL4, RPS9, RPS23, HMBS, HPRT1, EEF1A1 and UBI as internal control genes (ICG) to normalize the transcriptional data of mammary epithelial cells (MEC) in Indian cows. A total of 52 MEC samples were isolated from milk of Sahiwal cows (major indigenous dairy breed of India) across different stages of lactation: Early (5-15 days), Peak (30-60 days), Mid (100-140 days) and Late (> 240 days). Three different statistical algorithms: geNorm, Normfinder and BestKeeper were used to assess the suitability of these genes. In geNorm analysis, all the genes exhibited expression stability (M) values below 0.5 with EEF1A1 and RPL4 showing the maximum expression stability. Similar to geNorm, Normfinder also identified EEF1A1 and RPL4 as two of the most stable genes. In Bestkeeper algorithm as well, all the 10 genes showed consistent expression levels. The analysis showed that four genes, that is, EEF1A1, RPL4, GAPDH and ACTB exhibited higher coefficient of correlation to the Bestkeeper index, lower coefficient of variance and standard deviation, indicating their superiority to be used as ICG. The present analysis has provided evidence that RPL4, EEF1A1, GAPDH and ACTB could probably act as most suitable genes for normalizing the transcriptional data of milk-derived mammary epithelial cells of Indian cows.

Kinetics of lipogenic genes expression in milk purified mammary epithelial cells (MEC) across lactation and their correlation with milk and fat yield in buffalo

Expression patterns of lipogenic genes (LPL, ABCG2, ACSS2, ACACA, SCD, BDH, LIPIN1, SREBF1, PPARα and PPARγ) were studied in milk purified MEC across different stages of lactation (15, 30, 45, 60, 90, 120 and 240 days relative to parturition) in buffalo. PPARα was the most abundant gene while ABCG2 and ACSS2 had moderate level of expression; whereas expression of SREBF and PPARγ was very low. The expression patterns of some genes (BDH1, ACSS2, and LIPIN1) across lactation were positively correlated with milk yield while negatively correlated with fat yield. SCD also showed weak correlation with milk yield (p, 0.53) and fat yield (p, -0.47). On the other hand, expression pattern of ACACA was negatively correlated with milk yield (p, -0.88) and positively correlated with fat yield (p, 0.62). Strong correlation was observed between genes involved in de novo milk fat synthesis (BDH1, ACSS2, LIPIN2 and SCD) and milk yield.

Expression Analysis of Genes Associated with Prolificacy in FecB Carrier and Noncarrier Indian Sheep

ABSTRACT 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.

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.

Molecular Characterization of Buffalo Haptoglobin: Sequence Based Structural Comparison Indicates Convergent Evolution Between Ruminants and Human

ABSTRACT Haptoglobin (Hp) protein has high affinity for hemoglobin (Hb) binding during intravascular hemolysis and scavenges the hemoglobin induced free radicals. Earlier reports indicate about uniqueness of Hp molecule in human and cattle, but in other animals, it is not much studied. In this paper, we characterized buffalo Hp molecule and determined its molecular structure, evolutionary importance, and tissue expression. Comparative analysis and predicted domain structure indicated that the buffalo Hp has an internal duplicated region in α-chain only similar to an alternate Hp2 allele in human. This duplicated part encoded for an extra complement control protein CCP domain. Phylogenetic analysis revealed that buffalo and other ruminants were found to group together separated from all other non-ruminants, including human. The key amino acid residues involved in Hp and Hb as well as Hp and macrophage scavenger receptor, CD163 interactions in buffalo, depicted a significant variation in comparison to other non-ruminant species. Constitutive expression of Hp was also confirmed across all the vital tissues of buffalo, for the first time. Results revealed that buffalo Hp is both structurally and functionally conserved, having internal duplication in α-chain similar to human Hp2 and other ruminant species, which might have evolved separately as a convergent evolutionary process. Furthermore, the presence of extra Hp CCP domain possibly in all ruminants may have an effect during dimerization of molecule in these species.

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.

Transcriptome Analysis of Circulating PBMCs to Understand Mechanism of High Altitude Adaptation in Native Cattle of Ladakh Region

Ladakhi cattle is native population of Leh and Ladakh region and constantly exposed to hypobaric hypoxia over many generations. In present study, transcriptome signatures of cattle from Ladakh region (~5500 m) and Sahiwal cattle from tropical regions were evaluated using Agilent 44 K microarray chip. The top up-regulated genes in Ladakhi cows were INHBC, ITPRI, HECA, ABI3, GPR171, and HIF-1α involved in hypoxia and stress response. In Sahiwal cows, the top up-regulated genes eEF1A1, GRO1, CXCL2, DEFB3 and BOLA-DQA3 were associated with immune function and inflammatory response indicating their strong immune potential to combat the pathogens prevalent in the tropical conditions. The molecular pathways highly impacted were MAPK signaling, ETC, apoptosis, TLR signaling and NF- kB signaling pathway indicating signatures of adaptive evolution of these two cattle types in response to diverse environments. Further, qPCR analysis revealed increased expression of DEGs such as HIF-1, EPAS-1, VEGFA, NOS2, and GLUT-1/SLC2A1 in cattle types from high altitude suggesting their pivotal role in association with high altitude adaptation. Based on data generated, native cattle of Ladakh region was found to be genetically distinct from native cattle adapted to the tropical region of India.