Milk somatic cell derived transcriptome analysis identifies regulatory genes and pathways during lactation in Indian Sahiwal cattle (Bos indicus).

Abstract

The present study is an effort to understand the genomic drivers of lactation in Sahiwal (Bos indicus), the best milch cattle breed of the tropics. RNA sequencing of four animals from early, mid and late lactation stages was performed using milk somatic cells as source of RNA. The genes encoding the milk casein and whey proteins showed highest expression in early and mid lactation, with a declining trend towards the late stage. The enhanced expression of PLIN2, FABP5 and FABP3 genes in mid lactation suggests enrichment of the PPARα pathway which is linked to fatty acid metabolism. A gradual decline in the percentage of genes involved in metabolism of proteins, mRNA and insulin synthesis from early to late lactation reflected transition from lactogenesis to involution. Major biological pathways maintained throughout lactation were adaptive immune system, FGF signaling, EGFR signaling, activated TLR4 signaling, NFkB and MAP kinases activation mediated by TLR4 signaling repertoire. Differential expression analysis revealed 547, 1010 and 1313 differentially expressed genes (p < 0.05) between early-late, early-mid and mid-late stages, respectively. The topmost regulatory genes identified by network analysis from the differentially expressed genes, were involved in Chemokine receptor, GPCR and EGFR1 pathways. The genes and pathways delineated in this study have regulatory implications in cell morphogenesis, lipid droplet formation and protein synthesis in the course of lactation. The study provides an insight into the expression profile of genes influencing milk properties and lactation in Sahiwal cattle.