Deciphering the transcriptome of prepubertal buffalo mammary glands using RNA sequencing

Although water buffaloes are the main milk-producing animals in Indian subcontinent, only limited attempts have been made to identify canonical pathways and gene regulatory networks operating within the mammary glands of these animals. Such information is important for identifying unique transcriptome signatures in the mammary glands of diseased animals. In this report, we analyzed the transcription profile of 3 prepubertal buffalo mammary glands and identified common genes (mean FPKM > 0.2 in all samples) operating in the glands. Among 19,994 protein coding genes, 14,678 genes expressed and 5316 unique genes did not express in prepubertal buffalo mammary glands. Of these 14,678 expressed genes, 79% comprised a ubiquitous transcriptome that was dominated by very lowly expressed genes (51%). The percentage of rarely, moderately, and abundantly expressed genes was 25%, 2%, and 1%, respectively. Gene Ontology (GO) terms reflected in the expression of common genes (mean FPKM > 5.0) for molecular function were related to binding and catalytic activity. Products of these genes were involved in metabolic and cellular processes and belong to nucleic acid binding proteins. The canonical pathways for growth of mammary glands included integrin signaling, inflammation, GnRH and Wnt pathways. KEGG enriched pathways revealed many pathways of cancer including ribosome, splisosome, endocytosis, and ubiquitin-mediated proteolysis, pathways for viral infection, and bacterial invasion of epithelial. Highly expressed genes (mean FPKM > 500 included beta-actin (ACTB), beta-2 microglobulin (B2M), caseins (CSN2, CNS3), collagens (COL1A1, COL3A1), translation elongation factors (EEF1A1, EEF1G, EEF2), keratins (KRT15, KRT19), major histocompatibility complex genes (CD74, JSP.1), vimentin (VIM), and osteopontin (SPP1). Interestingly, expression of milk protein genes in prepubertal glands opens possible roles of these genes in development of mammary glands. We report the whole transcriptomic signature of prepubertal buffalo mammary gland and indicated its molecular signature is similar to cancer type.