Integrated analysis of long noncoding RNA and mRNA expression profiles reveals the potential role of long noncoding RNA in different bovine lactation stages

Long noncoding RNA (lncRNA) play a critical role in mammary development and breast cancer biology. Despite their important role in the mammary gland, little is known of the roles of lncRNA in bovine lactation, particularly regarding the molecular processes underlying it. To characterize the role of lncRNA in bovine lactation, 4 samples of Holstein cow mammary gland tissue at peak and late lactation stages were examined after biopsy. We then profiled the transcriptome of the mammary gland using RNA sequencing technology. Further, functional lncRNA-mRNA coexpression pairs were constructed to infer the function of lncRNA using a generalized linear model, followed by gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. More than 1,000 putative lncRNA were identified, 117 of which were differentially expressed between peak and late lactation stages. Bovine lncRNA were shorter, with fewer exon numbers, and expressed at significantly lower levels than protein-coding genes. Seventy-two differentially expressed (DE) lncRNA were coexpressed with 340 different protein-coding genes. The KEGG pathway analysis showed that target mRNA for DE lncRNA were mainly related to lipid and glucose metabolism, including the peroxisome proliferator-activated receptors and 5′ adenosine monophosphate-activated protein kinase signaling pathways. Further bioinformatics and integrative analyses revealed that 12 DE lncRNA potentially played important roles in bovine lactation. Our findings provide a valuable resource for future bovine transcriptome studies, facilitate the understanding of bovine lactation biology, and offer functional information for cattle lactation.