Transcriptomic profiles of the bovine mammary gland during lactation and the dry period

The initiation and maintenance of lactation are complex phenomena governed by biochemical and endocrine processes in the mammary gland (MG). Although DNA-based approaches have been used to study the onset of lactation, more comprehensive RNA-based techniques may be critical in furthering our understanding of gene alterations that occur to support lactation in the bovine MG. To further determine how gene profiles vary during lactation compared with the dry period, RNA-seq transcriptomic analysis was used to identify differentially expressed genes (DEG) in bovine MG tissues from animals that were lactating and not lactating. A total of 881 DEG (605 upregulated and 276 downregulated) were identified in MG of 3 lactating Chinese Holstein dairy cows versus the 3 dry cows. The subcellular analysis showed that the upregulated genes were most abundantly located in “integral to membrane” and “mitochondrion,” and the top number of downregulated genes existed in “nucleus” and “cytoplasm.” The functional analysis indicated that the DEG were primarily associated with the support of lactation processes. The genes in higher abundance were most related to “metabolic process,” “oxidation-reduction process,” “transport” and “signal transduction,” protein synthesis-related processes (transcription, translation, protein modifications), and some MG growth-associated processes (cell proliferation/cycle/apoptosis). The downregulated genes were mainly involved in immune-related processes (inflammatory/immune/defense responses). The KEGG analysis suggested that protein synthesis-related pathways (such as protein digestion and absorption; protein processing in endoplasmic reticulum; and glycine, serine, and threonine metabolism) were highly and significantly enriched in the bovine MG of lactating cows compared to dry cows. The results suggested that the dry cows had decreased capacity for protein synthesis, energy generation, and cell growth but enhanced immune response. Collectively, this reduced capacity in dry cows supports the physiological demands of the next lactation and the coordinated metabolic changes that occur to support these demands. A total of 51 identified DEG were validated by RT-PCR, and consistent results were found between RT-PCR and the transcriptomic analysis. This work provides a profile of gene-associated changes that occur during lactation and can be used to facilitate further investigation of the mechanisms underlying lactation in dairy cows.