Research article | Open | Published:
Functional development of the adult ovine mammary gland—insights from gene expression profiling
BMC Genomicsvolume 16, Article number: 748 (2015)
The mammary gland is a dynamic organ that undergoes dramatic physiological adaptations during the transition from late pregnancy to lactation. Investigation of the molecular basis of mammary development and function will provide fundamental insights into tissue remodelling as well as a better understanding of milk production and mammary disease. This is important to livestock production systems and human health.
Here we use RNA-seq to identify differences in gene expression in the ovine mammary gland between late pregnancy and lactation.
Between late pregnancy (135 days of gestation ± 2.4 SD) and lactation (15 days post partum ± 1.27 SD) 13 % of genes in the sheep genome were differentially expressed in the ovine mammary gland. In late pregnancy, cell proliferation, beta-oxidation of fatty acids and translation were identified as key biological processes. During lactation, high levels of milk fat synthesis were mirrored by enrichment of genes associated with fatty acid biosynthesis, transport and lipogenesis. Protein processing in the endoplasmic reticulum was enriched during lactation, likely in support of active milk protein synthesis. Hormone and growth factor signalling and activation of signal transduction pathways, including the JAK-STAT and PPAR pathways, were also differently regulated, indicating key roles for these pathways in functional development of the ovine mammary gland. Changes in the expression of epigenetic regulators, particularly chromatin remodellers, indicate a possible role in coordinating the large-scale transcriptional changes that appear to be required to switch mammary processes from growth and development during late pregnancy to synthesis and secretion of milk during lactation.
Coordinated transcriptional regulation of large numbers of genes is required to switch between mammary tissue establishment during late pregnancy, and activation and maintenance of milk production during lactation. Our findings indicate the remarkable plasticity of the mammary gland, and the coordinated regulation of multiple genes and pathways to begin milk production. Genes and pathways identified by the present study may be important for managing milk production and mammary development, and may inform studies of diseases affecting the mammary gland.
The mammary gland, and the physiological control of lactation, evolved as a vital part of the mammalian reproduction strategy [1, 2]. Milk provides an essential source of nutrients to newborn mammals, as well as immune factors, including anti-microbial, anti-inflammatory and immune-modulatory agents, that offer protection against infections, and have beneficial effects on intestinal flora and gut health [2, 3]. Human epidemiological studies provide evidence that breast milk plays a role in protection against gastrointestinal and respiratory tract infections in infants, and in programming metabolism and disease later in life [4–7]. Additionally, humans have long exploited the production of milk by domestic ruminants for the manufacture of dairy products, making milk an important part of human nutrition. Milk production is also important in pastoral livestock systems for the production of meat and fibre. In such systems, milk produced by the dam is the sole source of nutrients for newborn offspring and may influence their survival, growth to weaning, and future productive performance [8–10]. The ability to manipulate lactation performance is an area of increasing interest, and knowledge of the biological pathways and mechanisms that govern mammary gland development and lactation is commercially important.
The mammary gland displays a high level of developmental plasticity, able to undergo repeated cycles of growth, differentiation, and regression, coordinated by the reproductive state . Mammary development follows the reproductive cycle and is regulated by endocrine hormones [11–13]. The dynamic nature of the mammary gland makes it an ideal model for studying molecular regulation of development and cellular differentiation. Mammary glands are comprised of two main tissue components: the parenchyma and the stroma . The parenchyma contains the functional secretory and ductal tissue. The mammary stroma, or fat pad, contains the supportive tissue, including connective tissue, fibroblasts, adipose tissue, nerve tissue, and endothelial cells associated with blood vessels and lymph vessels [13, 14].
While mammogenesis is initiated during embryonic life, the majority of mammary development occurs post-natally, particularly during pregnancy when there is marked expansion of the lobulo-alveolar network (clusters of alveoli, which are spherical structures, comprised of secretory epithelial cells surrounding a central lumen, and basal myoepithelial cells which contract to allow milk let-down) . The number of mammary epithelial cells is correlated to milk yield [16–19], thus establishment of secretory tissue during this developmental stage is critical for subsequent lactational performance. Prior to parturition, functional differentiation of epithelial cells is initiated, referred to as lactogenesis stage 1 . This stage is characterised by cytological changes as well as increasing expression of genes encoding milk proteins such as caseins and beta-lactoglobulin [20–22]. The second stage of lactogenesis (lactogenesis stage two or secretory activation phase) is initiated at parturition and involves the closing of tight junctions between alveolar epithelial cells and further increases in expression of milk proteins, including LALBA (alpha-lactalbumin) which promotes lactose synthesis and associated increases in milk volume [21, 23].
In adult sheep, mammary development is essentially complete by parturition and there is only limited mammary growth during early lactation . This is in contrast to the more extensive continued development observed in the mammary glands of litter bearing species, including rodents , pigs , and some dairy cows and goats [16, 27]. During early lactation, cellular metabolic activity and nutrient transport increases to provide substrates for synthesis of milk components .
Such coordinated regulation of developmental and functional events is likely mediated by large-scale changes in gene expression [29, 30]. While the structural changes, and roles of hormones, in events from late pregnancy to lactation are well known [29, 31–33], little is known about the underlying molecular mechanisms that regulate differentiation of the mammary epithelium to a secretory phenotype.
Here we use RNA-seq to examine global gene expression in sheep mammary glands during the transition from late pregnancy to lactation providing insights into the physiological and metabolic adaptations that occur during this transition. Understanding the regulation of this transition is essential to our ability to intervene in mammary development and disease, and to manipulate lactation outcomes important to animal production and human nutrition.
Animals and sample collection and processing
The animal study was conducted at the Massey University Keeble Sheep and Beef farm, 5 km south of Palmerston North, New Zealand. The study was approved by the Massey University Animal Ethics Committee, Palmerston North, New Zealand.
Mammary tissue was sampled from 2-year-old, primiparous, Romney ewes during late pregnancy (n = 27, ewe age: 733.9 days ± 1.66 SD and 135 days of gestation ± 2.4 SD) and again during early lactation (n = 18, ewe age: 761.0 ± 2.11 SD days and 15 ± 1.27 SD days post partum). Mammary parenchymal tissue (30–50 mg) was collected with a needle biopsy (Bard® Magnum® reusable core biopsy gun and 12G, 10 cm core biopsy needles, Bard Biopsy Systems, AZ, USA) as described by Norgaard et al. . Tissue samples were immediately frozen in liquid nitrogen, then stored at − 80 °C until RNA extraction.
Total RNA was isolated from mammary tissue samples using Trizol (Invitrogen, CA, USA) and purified using RNeasy mini kit (Qiagen, Netherlands). Genomic DNA contamination was eliminated via on-column digestion with DNase (Qiagen), as per the manufacturer’s protocol, and verified by RT-PCR. The concentration and quality of RNA was measured using a Nanodrop ND-1000 spectrophotometer (Thermo Scientific, MA, USA) and integrity was assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, CA, USA). Only RNA with RNA integrity numbers (RINs) above seven was used in this study. One μg of total RNA was used as template to perform cDNA synthesis using the SuperScript VILO cDNA Synthesis Kit (Invitrogen) as per the manufacturer’s protocol.
RNA sequencing and data analysis
For RNA sequencing, we attempted to minimise the effects of individual variation between animals by pooling RNA, as described by Paten et al. . Briefly, each pool contained a total of 2 μg of RNA, sampled from three randomly selected animals. Nine pools were generated for late pregnancy samples (n = 27) and six for lactation samples (n = 18).
RNA-seq data was generated from pooled RNA (as detailed above) using an Illumina Hi-Seq 2000 (GEO accession number GSE71424). Mapping of the 100 bp paired-end sequence reads to the Ovis aries genome (version 3.2) and analysis of expression data was carried out using CLC Genomics Workbench (CLC Bio). Read mapping statistics are supplied in Table S1 (Additional file 1: Table S1). Data were subjected to quantile normalization and RPKM (reads per kilobase (kb) per million mapped reads statistic (RPKM = total exon reads mapped/ mapped reads in millions × exon length in kb)) for each annotated gene (10,175 genes) were calculated. Differentially expressed genes were identified using the Baggerly test . The Baggerly test is similar to a two sample t-test but the test statistic is weighted according to the number of reads in each sample . P-values were corrected for multiple testing using the method of Benjamini and Hochberg . Genes that had a corrected P value of less than 0.05 were considered to be differentially expressed (Additional file 2: Tables S2 and S3). Principal components and heat map analyses were performed to evaluate general patterns of variation in gene expression between treatment groups (Additional file 3: Figure S1). The findings of our RNA-seq data were validated by reverse transcription quantitative (RT-qPCR) analysis of a sub-set of 18 genes, expression was normalized against established reference genes for the ovine mammary gland . The details of which are reported in Additional file 4.
Gene ontology and pathways analyses
Analysis of the RNA-seq data provided a list of differentially expressed genes. In order to contextualise this with physiological and metabolic events occurring in the mammary gland during late pregnancy and lactation, DAVID [38, 39] was used to identify enriched gene ontology (GO) and perform KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways analyses (Fisher’s exact P-value < 0.01, no correction for multiple testing was applied). One limitation to these analyses is that, at the time of this analysis, many sheep genes did not have assigned functional categories assigned to them. To overcome this, mouse orthologs of sheep genes were identified based on homolog gene annotation in the NCBI gene database (http://www.ncbi.nlm.nih.gov/gene). Gene lists consisting of genes more highly expressed during late pregnancy (Additional file 2: Table S2) and lactation (Additional file 2: Table S3) were submitted to DAVID (http://david.abcc.ncifcrf.gov/tools.jsp) and a background list was constructed consisting of all of the genes expressed in the ovine mammary gland for which a mouse gene ortholog could be identified (n = 6733). Functional annotation clustering was used to cluster similar GO terms together and results were ranked according to the Group Enrichment Score (the geometric mean (in-log scale) of member’s P-values in a corresponding annotation cluster). Functional annotation clusters were given an overall term which summarised the general theme of each cluster and any clusters with enrichment scores below 2.0 (P = 0.01) were discarded (Additional file 5: Tables S5 and S6).
Results and discussion
We detected a total of 10,132 genes expressed in ovine mammary gland tissue during late pregnancy and 10,096 genes expressed during lactation (9622 genes expressed in both late pregnancy and lactation), which accounts for approximately 49 % of genes in the sheep genome. Unfortunately a number of these genes have not been functionally annotated, including 4717 genes that were either predicted or uncharacterised. Validation of a subset of these genes by RT-qPCR (Additional file 4: Figure S2) indicates that our RNA-seq data is high quality and biologically relevant.
Highly expressed genes
A very small number of genes accounts for 60 % of all RNA-seq reads (30 and 24 genes in late pregnancy and lactation, respectively). The six most highly expressed genes during lactation (BLG (β-lactoglobulin), CSN2 (β-casein), CSN1S1 (α-S1-casein), LALBA (α-lactalbumin), CSN3 (κ-casein), GLYCAM1 (glycosylation dependent cell adhesion molecule-1) and CSN1S2 (casein-α-S2)) have been reported to be highly expressed in cells derived from milk in other species (e.g. cattle  and human ) using RNA-sequencing. High-expression of some of these genes has been detected using other experimental approaches, including RT-qPCR (e.g. [42, 43]), supporting the conclusion that the ultra-high expression that we report for these genes reflects the biology of the mammary gland in late pregancy and lactation and is not an experimental artifact.
This sub-set of ultra-highly expressed genes included the major milk protein genes (caseins, α-lactalbumin and β-lactoglobulin), genes encoding ribosomal proteins, genes with products involved in energy metabolism and genes with products involved in immunity and inflammation (Fig. 1).
During late pregnancy, 18 ribosomal protein genes were ultra-highly expressed (Fig. 1) and may reflect increased levels of ribosome biogenesis in preparation for high levels of milk protein synthesis. During lactation, several metabolism-related genes were also ultra-highly expressed, which reflects the increased metabolic activity of the mammary gland . Highly expressed metabolism-related genes included COX3, ND2 and ATP6, which encode enzymes of the oxidative phosphorylation pathway, and FASN and FABP3, which are involved in milk-fat synthesis.
Differentially expressed genes
The high expression of a small number of genes effectively dilutes the expression of other genes in the mammary gland  and may limit the ability to detect differences in their expression as well as introduce bias into normalization strategies. To identify differentially expressed genes analysis was performed both including and excluding these highly expressed transcripts. Although there were subtle differences to the fold-changes and RPKMs reported in the two analyses, the inclusion of the highly expressed genes did not alter the identity or number of genes identified as differentially expressed and here focus on the analysis with the highly expressed genes included. Despite this potential limitation, we were able to detect a large proportion (27 %, 2750 genes) of genes in the ovine mammary gland that were differentially expressed between late pregnancy and early lactation. A greater proportion of these genes decreased expression between late pregnancy and early lactation (1510 genes, 55 % of differentially expressed genes) while 1240 genes (45 % of differentially expressed genes) increased in expression between late pregnancy and early lactation. This trend is consistent with that of Finucane et al.  who used microarray analysis to examine the bovine mammary gland during lactation. Late pregnancy is a critical period for mammary development, in which the functional secretory tissue of the gland is established and differentiates prior to parturition. In sheep , there is limited mammary development post-partum, which may explain why the number of genes expressed in the mammary gland appears to be greater during pregnancy than in lactation. The study of Finucane et al.  drew parallels to small-animal models (mouse [47, 48], rat , and rabbit ), in which milk protein expression is initiated during mid-pregnancy and increases throughout pregnancy, plateauing during lactation. Thus, similarly, in sheep, the genes required for lactation may already be adequately expressed during late pregnancy.
Gene ontology analysis using DAVID revealed that nine functional annotation clusters (enrichment score (ES) > 2, equivalent to P < 0.01 ) and 11 KEGG pathways (P < 0.01) were enriched in mammary glands during late pregnancy (Figs. 2a–c and Additional files 4 and 5: Tables S5 and S7). Enriched genes were associated with gene ontology categories involved in energy metabolism (mitochondria, glycolysis, and fatty acid metabolism; in particular beta-oxidation), protein synthesis (translation initiation), cell proliferation (cell cycle, p53 signalling, chromosome and cytoskeleton) and response to hormone stimulus.
A second set of nine functional annotation clusters (ES > 2) and five KEGG pathways (P < 0.01) were enriched in lactating mammary glands (Figs. 2d–f and Additional files 4 and 5: Tables S6 and S8). The most enriched gene ontology categories processes were associated with lipid synthesis and secretion (endoplasmic reticulum, fatty acid biosynthesis, lipid synthesis, synthesis, and vesicle). Additionally, there was enrichment for genes associated with synthesis and post-translational processing of proteins (unfolded protein response, heat shock protein, amino acid biosynthesis).
Expression of genes involved in cellular proliferation, survival and apoptosis in the mammary gland
It is well-known that milk production is a function of the number and activity of secretory epithelial cells [16–19], thus the mammary tissue establishment phase that occurs during pregnancy is crucial to subsequent lactation performance. Physiologically, late pregnancy is characterised by extensive structural remodelling, involving expansion of the lobulo-alveolar network, and functional differentiation of alveolar cells to secretory cells in preparation for producing milk in lactation [21, 24]. On a molecular level this is reflected in the present study through enrichment of genes associated with: regulation of cell cycle, including cyclins and cyclin-dependent kinases which form complexes to regulate cell cycle progression ; DNA synthesis including genes encoding the components of the MCM complex, which form part of the pre-replication complex , and; genes encoding histone fold proteins (POLE3 and CHRAC1) (Fig. 3 and full list of genes in Additional file 5: Table S5). Also enriched were mitosis related genes, specifically those associated with the spindle assembly checkpoint, which ensures proper chromosome alignment and segregation during mitotic cell division [53, 54] (Fig. 3 and Additional file 5: Table S5). These genes are regulated by E2F family and Myc transcription factors . Expression of MYC was higher during late pregnancy, while none of the E2F transcription factors differed in expression. However, E2F is regulated by the Retinoblastoma tumour suppressor protein (pRB), such that when pRB is hypophosphorylated it binds to E2F, inhibiting transcription of S phase target genes . Genes encoding mitogenic factors known to promote hyperphosphorylation of Rb and E2F activity were more highly expressed during late pregnancy (Fig. 3), confirming E2F activity.
Given the increased DNA replication, cell division and proliferation occurring during late pregnancy, mechanisms need to be in place to account for DNA damage repair and ensure proliferation occurs in a controlled and regulated manner. The p53 signalling pathway was enriched during late pregnancy, such that expression of TP53 was higher during lactation, as were p53-responsive DNA repair genes such as GADD45A, GADD45G, CCNG1 and CCNG2 (Additional file 6: Table S7). GADD45A may also use its DNA repair function to remove DNA methylation marks, potentially playing a role in epigenetic regulation of gene expression .
Our RNA-seq data implies that during lactation, cellular growth and proliferation processes are down-regulated, indicating that continued growth of mammary tissue is limited after parturition in sheep. These findings are consistent with transcriptomic studies in other species [46, 58, 59] and observations of cellular proliferation rates (using Ki-67) in cows, which have demonstrated a reduction in epithelial cellular expansion during lactation compared with late pregnancy . In the present study we identified cell-survival markers, such as BCL2L15, BCAR1, BAG1, BAG3, BAG5 and BEX2, as more highly expressed in mammary tissue during lactation compared with late pregnancy. The products of these genes may be involved in maintaining the population of milk secretory cells and milk production (galactopoiesis) during lactation, and may therefore serve as useful biomarkers or target genes for manipulation of lactation persistence.
Expression of genes involved in energy and fat metabolism in the mammary gland
Limited growth of the mammary gland during lactation [24, 46, 60] would imply that greater activity of secretory cells is required to increase milk yield to peak lactation . In support of this, we detected up-regulation of genes linked to metabolism and transport processes during lactation, consistent with studies in bovine  and murine [61, 62] mammary glands. In particular, genes involved in fat metabolism were enriched during lactation (full list of genes in Additional file 5: Table S6), such that there was higher expression of genes encoding enzymes involved in fatty acid biosynthesis, fatty acid activation, glycerol synthesis, triacylglycerol synthesis, cholesterol synthesis, milk fat secretion and genes involved in signal transduction of lipogenic pathways (e.g., SREBF1, PPARD, INSIG1 and ESR1). Genes involved in fatty acid catabolic pathways, particularly beta-oxidation, were down regulated during lactation (Fig. 4 and Additional file 5: Table S5). These findings are consistent with the large amount of milk fat produced by the mammary gland during lactation and with findings in mice , pigs  and cows [46, 58].
Genes involved in energy metabolism, particularly genes encoding components of the mitochondrial respiratory chain, were differentially, yet highly expressed during both late pregnancy and lactation, reflecting the high-energy demands of cellular division and growth during late pregnancy, milk synthesis and secretion during lactation (Fig. 4 and Additional file 5: Tables S5 and S6). Higher expression of genes involved in amino acid biosynthesis and amino acid catabolism during lactation (Additional file 5: Table S6) is consistent with sparing of fatty acids by switching to amino acids as substrates for energy generation, as seen in bovine mammary glands .
Expression of genes involved in protein synthesis in the mammary gland
During both late pregnancy and lactation, we detected increased expression of genes associated with different aspects of protein synthesis. During late pregnancy, one of the most enriched categories of genes (functional annotation clustering enrichment score of 4.02) included genes involved in translation, translation elongation and translation initiation (Fig. 5 and Additional file 4: Table S5). Whilst the increased expression of mRNA transcripts does not necessarily reflect an increase in the levels of active protein, the large number of changes we observe in expression of genes involved in protein synthesis may reflect an increased translational requirement associated with milk proteins and enzymes being synthesised in preparation for lactation. Previous studies have also detected a decrease in expression of protein synthesis components (including genes involved in translational elongation) during lactation [43, 63]. It has been hypothesised that this down-regulation may be a mechanism for the mammary gland to prioritise translation of milk specific proteins  this may account for the significant changes seen in milk protein composition but only slight increases in milk protein synthesis in certain experimental regimes .
In rodents, the JAK2-STAT5 pathway plays an essential role in mediating expression of milk-protein genes and expression of STAT5 genes increases during lactation in response to stimulus from lactogenic hormones. In cows, however, JAK2-STAT5 signalling appears to have only a limited role in milk-protein synthesis . In the present study, we detected higher expression of STAT5 genes (STAT5A and STAT5B: 1.5 and 1.4 fold, respectively, FDR P < 0.01) during late pregnancy, compared with lactation in sheep, which may be related to the role of STAT5 in mammary growth [29, 65, 66] as well as induction of milk-protein gene transcription. The lower expression of STAT5 during lactation, despite the dramatic increase in expression of major milk-protein genes, may indicate that, similarly to cows , STAT5 may not play a central role in regulating milk-protein synthesis during lactation in sheep. However, it is important to consider that regulation of STAT5 activity is predominantly through protein activation [67, 68]. Bionaz and Loor  suggested a possible role of STAT5 in regulating milk-protein synthesis through the E74-like factor 5 transcription factor, encoded by the ELF5 gene. Consistent with this, expression of ELF5 was increased during lactation (2.2 fold, FDR P < 0.001) in the present study, demonstrating a role of ELF5 and possibly STAT5 in transcriptional regulation of milk protein synthesis in the sheep.
Our data showed that during lactation there was also enrichment of post-translational processing of proteins, such that there was higher expression of genes in the unfolded protein response (UPR) pathway in the endoplasmic reticulum (ER) (e.g., ATF4 transcription factor and EIF2AK2, also known as PERK) and the ER-associated degradation (ERAD) pathway, which are activated as a quality control during high levels of protein synthesis (Fig. 5). Genes encoding molecular chaperones (e.g., HSPA8, HSPA5, DNAJB2 and DNAJC3) were more highly expressed in the lactating mammary gland along with genes encoding members of the HSP70 family (genes listed in Additional file 5: Table S6), which are activated at the transcriptional level by endoplasmic stress. In addition to their more traditional molecular chaperone role regulating protein folding and processing, heat shock proteins (HSPs) are also thought to have functions in immunity, inflammation and suppression of apoptosis . HSPA5 and HYOU1 may function in angiogenesis in the mammary gland through promoting VEGF processing and signalling (via MAPK) inducing endothelial cell proliferation . These genes may be useful targets for enhancing blood flow, and thus, nutrient supply to the mammary gland for milk production.
Genes associated with biosynthesis and transport of amino acids were more highly expressed in the ovine mammary gland during lactation, compared with late pregnancy (Additional file 5: Table S6). Supplementation of amino acids can improve milk production [70, 71] and bioavailability of amino acids, particularly lysine and methionine, and their rate of transport into MECs can be a major limiting factor for milk protein synthesis. Thus, in addition to amino acid supplementation, targeted expression and translation of amino acid transporter genes may be a useful intervention to improve milk protein synthesis and potentially promote synthesis of specific milk proteins.
Mechanisms of transcriptional regulation: hormones and epigenetics
It is apparent that coordinated shifts in gene expression govern the switch from proliferation and differentiation to secretion in the mammary gland. These changes must be regulated by mechanisms acting at a higher level, such as endocrine and/or epigenetic regulation. The role of hormones and growth factors in the mammary gland is well known [29, 31–33], and there is growing interest in the role of epigenetics in mammary development and lactation [12, 72–74].
During both late pregnancy and lactation there was enrichment of genes associated with hormone-signalling pathways (genes listed in Additional file 4: Tables S5 and S6). Genes associated with oestrogen and corticosteroids in particular, were more highly expressed during late pregnancy compared with lactation; however, expression of genes encoding the hormones and hormone receptors themselves did not generally differ. Putative oestrogen responsive genes have been identified in the mammary gland of cattle , but we do not see a clear pattern of expression of these genes in our dataset; with seven out of 70 oestrogen responsive genes being more highly expressed in late pregnancy and ten out of 70 more highly expressed during lactation. This may reflect species-specific differences in genes that are responsive to oestrogen, or may reflect the complexity of oestrogen signalling to the mammary gland . Similarly we did not see a clear pattern of expression associated with glucocorticoid exposure .
Genes associated with growth factor activity, such as IGF binding proteins two, four and seven, and V-Erb-B2 Avian Erythroblastic Leukemia Viral Oncogene Homolog 2 (ERBB2), which encodes a member of the epidermal growth factor (EGF) family also known as HER2, were more highly expressed during late pregnancy. Both the IGFs and EGFs are potent mitogens and have been implicated in mediating hormone effects on mammary growth [65, 77].
A number of genes encoding hormone receptors were found to be more highly expressed during lactation, e.g., oestrogen type 1 receptor (ESR1), thyroid hormone receptor (ERBA BETA1), growth hormone receptor (GHR), insulin receptor substrate 1 (IRS1), glucocorticoid receptor subfamily 3, group c, members 1 and 2 (NR3C1 and NR3C2), transforming growth factor beta receptor type 1 (TGFBR1), as well as parathyroid hormone-related protein (PTH-RP), indicating there may be increased sensitivity of the mammary gland to these hormones during lactation. Many of these encode receptors that mediate the effects of lactogenic and galactopoietic hormones (e.g., insulin, GH and glucocorticoids) and are important for the initiation and maintenance of milk production [33, 65].
In addition to regulation of transcription by hormones, the transcriptional potential of cells in the mammary gland may be influenced by epigenetics. A sub-set of differentially expressed genes, more highly expressed during late pregnancy, were associated with chromosome and chromatin organisation (genes listed in Additional file 5: Table S5). While many of these genes may simply be related to cell division (chromosome replication and segregation) during the mammary growth phase, some may be involved in epigenetic regulation of gene expression. Alterations to the structure of chromatin can influence the recruitment of transcription factors to DNA, which may impact upon the expression of genes. Differentially expressed genes which were associated with chromatin remodelling included: CENPA, H3F3A and H3F3B which encode variants of the histone H3 protein, HDAC2, a histone deacetylase associated with transcriptional repression of genes, and CBX3 (chromobox homolog 3), which has recently been shown to bind to gene bodies and play a role regulating genes, such as cell cycle associated genes, through transcriptional regulation, RNA processing and alternative splicing . Additionally, another set of genes more highly expressed during late pregnancy were associated with RNA splicing (e.g., HNRNPA3, HNRNPM, HNRNPF, HNRNPU, HNRNPH3, HNRNPD, HNPNPA0 and CD2BP2). We also identified a small, but significant (p = 0.03 for late pregnancy, p = 7.6 × 10−3 for lactation), number of gene clusters or “gene neighbourhoods” in our sets of differentially expressed genes (Additional file 7: Table S9). These clusters may function as chromatin domains or gene neighbourhoods that may be epigenetically regulated . Collectively, these results may indicate a coordinated role for genes involved in chromatin remodelling and RNA splicing in the transcriptional regulation of the mammary gland during the transition period from late pregnancy to lactation.
This is the first study on the global expression profile of the ovine mammary gland during late pregnancy and lactation. We demonstrated that 27 % of genes expressed in the ovine mammary gland are differentially expressed between late pregnancy and lactation. Our findings indicated a strong transcriptional regulation of cell proliferation, lipid metabolism and protein translation and processing, such that expression of genes involved in cell cycle, translation and fatty acid catabolism were down-regulated, while expression of genes involved in fatty acid and amino acid biosynthesis and transport, lipogenesis, and protein processing, were up-regulated during lactation. Furthermore, hormones and growth factors, signalling pathways, e.g. JAK-STAT, SREBF and PPARD, and epigenetic regulation were highlighted as having a potential key role in mediating the adaptive transcriptional changes undertaken by the ovine mammary gland to support lactation. The identification of enriched genes and pathways in the present study will provide a platform for future research into the management of mammary development, function and disease.
Capuco AV, Akers RM. The origin and evolution of lactation. J Biol. 2009;8(4):37.
Goldman AS, Chheda S, Garofalo R. Evolution of immunologic functions of the mammary gland and the postnatal development of immunity. Pediatr Res. 1998;43(2):155–62.
Goldman AS. Modulation of the gastrointestinal tract of infants by human milk. interfaces and interactions. an evolutionary perspective. J Nutr. 2000;130(2):426.
Goldman AS. The immune system of human milk: antimicrobial, antiinflammatory and immunomodulating properties. Pediatr Infect Dis J. 1993;12(8):664–72.
Hanson LA. Breastfeeding provides passive and likely long-lasting active immunity. Ann Allergy Asthma Immunol. 1998;81(6):523–37.
Jackson KM, Nazar AM. Breastfeeding, the immune response, and long-term health. JAOA. 2006;106(4):203–7.
Oddy WH. Breastfeeding protects against illness and infection in infants and children: a review of the evidence. Breastfeed Rev. 2001;9(2):11.
Degen AA, Benjamin R. Milk and herbage intakes and growth rate of lambs from 32 to 130 days of age raised on natural pasture in the semi-arid Negev. Small Rumin Res. 2005;58(1):39–45.
Jordan D, Mayer D. Effects of udder damage and nutritional plane on milk yield, lamb survival and lamb growth of Merinos. Anim Prod Sci. 1989;29(3):315–20.
Mellor DJ. Nutritional and placental determinants of fetal growth-rate in sheep and consequences for the newborn lamb. Brit Vet J. 1983;139(4):307–24.
Knight CH, Peaker M. Development of the mammary-gland. J Reprod Fertil. 1982;65(2):521–36.
Capuco AV, Akers RM: Management and environmental influences on mammary gland development and milk production. In: Managing the Prenatal Environment to Enhance Livestock Productivity. Netherlands: Springer; 2010: 259–292.
Hurley WL, Loor JJ. Mammary gland | Growth, development and involution. In: Fuquay JW, editor. Encyclopedia of Dairy Sciences. Secondth ed. San Diego: Academic Press; 2011. p. 338–45.
Nickerson SC, Akers RM. Mammary Gland| Anatomy. In: Fuquay JW, editor. Encyclopedia of Dairy Sciences. Secondth ed. San Diego: Academic Press; 2011. p. 328–37.
McManaman JL, Neville MC. Mammary physiology and milk secretion. Adv Drug Deliv Rev. 2003;55(5):629–41.
Boutinaud M, Guinard-Flament J, Jammes H. The number and activity of mammary epithelial cells, determining factors for milk production. Reprod Nutr Dev. 2004;44(5):499–508.
Capuco A, Wood D, Baldwin R, McLeod K, Paape M. Mammary cell number, proliferation, and apoptosis during a bovine lactation: relation to milk production and effect of bST. J Dairy Sci. 2001;84(10):2177–87.
Keys JE, Capuco AV, Akers RM, Djiane J. Comparative study of mammary gland development and differentiation between beef and dairy heifers. Domest Anim Endocrinol. 1989;6(4):311–9.
Knight CH, Peaker M. Mammary development and regression during lactation in goats in relation to milk secretion. Q J EXP PHYSIOL CMS. 1984;69(2):331–8. Cambridge, England.
Tucker HA. Physiological control of mammary growth, lactogenesis, and lactation. J Dairy Sci. 1981;64(6):1403–21.
Akers RM, Capuco AV. Lactation. In: Hubert R, editor. Encyclopedia of Dairy Sciences. Oxford: Elsevier; 2002. p. 1442–6.
Brisken C, Rajaram RD. Alveolar and lactogenic differentiation. J Mammary Gland Biol Neoplasia. 2006;11(3–4):239–48.
Neville MC, McFadden TB, Forsyth I. Hormonal regulation of mammary differentiation and milk secretion. J Mammary Gland Biol Neoplasia. 2002;7(1):49–66.
Anderson RR. Mammary-gland growth in sheep. J Anim Sci. 1975;41(1):118–23.
Tucker H. Factors affecting mammary gland cell numbers. J Dairy Sci. 1969;52(5):720–9.
Kim SW, Hurley WL, Han IK, Easter RA. Changes in tissue composition associated with mammary gland growth during lactation in sows. J Anim Sci. 1999;77(9):2510–6.
Knight C. The importance of cell division in udder development and lactation. Livest Prod Sci. 2000;66(2):169–76.
Capuco AV, Ellis SE, Hale SA, Long E, Erdman RA, Zhao X, et al. Lactation persistency: Insights from mammary cell proliferation studies. J Anim Sci. 2003;81 suppl 3:18–31.
Hennighausen L, Robinson GW. Signaling pathways in mammary gland development. Dev Cell. 2001;1(4):467–75.
Hennighausen L, Robinson GW. Information networks in the mammary gland. Nat Rev Mol Cell Biol. 2005;6(9):715–25.
Collier RJ, McNamara J, Wallace CR, Dehoff MH. A review of endocrine regulation ofmetabolism during lactation. J Anim Sci. 1984;59(2):498–510.
Hovey RC, Trott JF, Vonderhaar BK. Establishing a framework for the functional mammary gland: From endocrinology to morphology. J Mammary Gland Biol Neoplasia. 2002;7(1):17–38.
Tucker H. Hormones, mammary growth, and lactation: a 41-year perspective. J Dairy Sci. 2000;83(4):874–84.
Norgaard JV, Nielsen MO, Theil PK, Sorensen MT, Safayi S, Sejrsen K. Development of mammary glands of fat sheep submitted to restricted feeding during late pregnancy. Small Rumin Res. 2008;76(3):155–65.
Paten AM, Pain SJ, Peterson SW, Blair HT, Kenyon PR, Dearden PK, et al. Identification of reference genes for RT-qPCR in ovine mammary tissue during late-pregnancy, lactation and in response to maternal nutritional programming. Physiol Genomics. 2014;46(15):560–70. doi:10.1152/physiolgenomics.00030.2014.
Baggerly KA, Deng L, Morris JS, Aldaz CM. Differential expression in SAGE: accounting for normal between-library variation. Bioinformatics. 2003;19(12):1477–83.
Benjamini Y, Hochberg Y. Controlling the false discovery rate - a practical and powerful approach to multiple testing. J Roy Stat Soc B Met. 1995;57(1):289–300.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4(1):44–57.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.
Wickramasinghe S, Rincon G, Islas-Trejo A, Medrano JF. Transcriptional profiling of bovine milk using RNA sequencing. BMC Genomics. 2012;13:45.
Lemay DG, Ballard OA, Hughes MA, Morrow AL, Horseman ND, Nommsen-Rivers LA. RNA sequencing of the human milk fat layer transcriptome reveals distinct gene expression profiles at three stages of lactation. Plos One. 2013;8(7):e67531.
Bionaz M, Loor JJ. Gene networks driving bovine milk fat synthesis during the lactation cycle. BMC Genomics. 2008;9:366.
Bionaz M, Loor JJ. Gene networks driving bovine mammary protein synthesis during the lactation cycle. Bioinf Biol Insights. 2011;5:83–98.
Bauman DE, Bruce Currie W. Partitioning of nutrients during pregnancy and lactation: a review of mechanisms involving homeostasis and homeorhesis. J Dairy Sci. 1980;63(9):1514–29.
Bionaz M, Loor JJ. Identification of reference genes for quantitative real-time PCR in the bovine mammary gland during the lactation cycle. Physiol Genomics. 2007;29(3):312–9.
Finucane KA, McFadden TB, Bond JP, Kennelly JJ, Zhao F-Q. Onset of lactation in the bovine mammary gland: gene expression profiling indicates a strong inhibition of gene expression in cell proliferation. Funct Integr Genomics. 2008;8(3):251–64.
Rijnkels M, Wheeler D, De Boer H, Pieper F. Structure and expression of the mouse casein gene locus. Mamm Genome. 1997;8(1):9–15.
Robinson GW, McKnight RA, Smith GH, Hennighausen L. Mammary epithelial cells undergo secretory differentiation in cycling virgins but require pregnancy for the establishment of terminal differentiation. Development. 1995;121(7):2079–90.
Rosen J, Woo S, Comstock J. Regulation of casein messenger RNA during the development of the rat mammary gland. Biochemistry. 1975;14(13):2895–903.
Shuster RC, Houdebine LM, Gaye P. Studies on the synthesis of casein messenger RNA during pregnancy in the rabbit. Eur J Biochem. 1976;71(1):193–9.
Sherr CJ, Roberts JM. Living with or without cyclins and cyclin-dependent kinases. Genes Dev. 2004;18(22):2699–711.
Tye BK. MCM proteins in DNA replication. Annu Rev Biochem. 1999;68(1):649–86.
Musacchio A, Hardwick KG. The spindle checkpoint: structural insights into dynamic signalling. Nat Rev Mol Cell Biol. 2002;3(10):731–41.
Taylor SS, Scott MI, Holland AJ. The spindle checkpoint: a quality control mechanism which ensures accurate chromosome segregation. Chromosom Res. 2004;12(6):599–616.
Bracken AP, Ciro M, Cocito A, Helin K. E2F target genes: unraveling the biology. Trends Biochem Sci. 2004;29(8):409–17.
Slansky J, Farnham P. Introduction to the E2F family: protein structure and gene regulation. In: Transcriptional Control of Cell Growth. Springer; 1996. p. 1–30.
Barreto G, Schafer A, Marhold J, Stach D, Swaminathan SK, Handa V, et al. Gadd45a promotes epigenetic gene activation by repair-mediated DNA demethylation. Nature. 2007;445(7128):671–5.
Bionaz M, Periasamy K, Rodriguez-Zas SL, Everts RE, Lewin HA, Hurley WL, et al. Old and New Stories: Revelations from Functional Analysis of the Bovine Mammary Transcriptome during the Lactation Cycle. Plos One. 2012;7(3):e33268.
Shu D-P, Chen B-L, Hong J, Liu P-P, Hou D-X, Huang X, et al. Global transcriptional profiling in porcine mammary glands from late pregnancy to peak lactation. OMICS. 2012;16(3):123–37.
Sorensen MT, Norgaard JV, Theil PK, Vestergaard M, Sejrsen K. Cell turnover and activity in mammary tissue during lactation and the dry period in dairy cows. J Dairy Sci. 2006;89(12):4632–9.
Rudolph MC, McManaman JL, Hunter L, Phang T, Neville MC. Functional development of the mammary gland: Use of expression profiling and trajectory clustering to reveal changes in gene expression during pregnancy, lactation, and involution. J Mammary Gland Biol Neoplasia. 2003;8(3):287–307.
Rudolph MC, Neville MC, Anderson SM. Lipid synthesis in lactation: diet and the fatty acid switch. J Mammary Gland Biol Neoplasia. 2007;12(4):269–81.
Lemay DG, Neville MC, Rudolph MC, Pollard KS, German JB. Gene regulatory networks in lactation: identification of global principles using bioinformatics. BMC Syst Biol. 2007;1:56.
Brophy B, Smolenski G, Wheeler T, Wells D, L’Huillier P, Laible G. Cloned transgenic cattle produce milk with higher levels of beta-casein and kappa-casein. Nat Biotechnol. 2003;21(2):157–62.
Akers RM. Major advances associated with hormone and growth factor regulation of mammary growth and lactation in dairy cows. J Dairy Sci. 2006;89(4):1222–34.
Naylor MJ, Oakes SR, Gardiner-Garden M, Harris J, Blazek K, Ho TWC, et al. Transcriptional changes underlying the secretory activation phase of mammary gland development. Mol Endocrinol. 2005;19(7):1868–83.
Wheeler TT, Broadhurst MK, Sadowski HB, Farr VC, Prosser CG. Stat5 phosphorylation status and DNA-binding activity in the bovine and murine mammary glands. Mol Cell Endocrinol. 2001;176(1):39–48.
Yang J, Kennelly J, Baracos V. The activity of transcription factor Stat5 responds to prolactin, growth hormone, and IGF-I in rat and bovine mammary explant culture. J Anim Sci. 2000;78(12):3114–25.
Lee AS. Glucose-regulated proteins in cancer: molecular mechanisms and therapeutic potential. Nat Rev Cancer. 2014;14(4):263–76.
Nichols JR, Schingoethe DJ, Maiga HA, Brouk MJ, Piepenbrink MS. Evaluation of Corn Distillers Grains and Ruminally Protected Lysine and Methionine for Lactating Dairy Cows. J Dairy Sci. 1998;81(2):482–91.
Schwab CG, Satter LD, Clay AB. Response of lactating dairy cows to abomasal infusion of amino acids. J Dairy Sci. 1976;59(7):1254–70.
Park CS. Role of compensatory mammary growth in epigenetic control of gene expression. Faseb J. 2005;19(12):1586–91.
Rijnkels M, Kabotyanski E, Montazer-Torbati MB, Beauvais CH, Vassetzky Y, Rosen JM, et al. The epigenetic landscape of mammary gland development and functional differentiation. J Mammary Gland Biol Neoplasia. 2010;15(1):85–100.
Singh K, Erdman RA, Swanson KM, Molenaar AJ, Maqbool NJ, Wheeler TT, et al. Epigenetic Regulation of Milk Production in Dairy Cows. J Mammary Gland Biol Neoplasia. 2010;15(1):101–12.
Li RW, Meyer MJ, Van Tassell CP, Sonstegard TS, Connor EE, Van Amburgh ME, et al. Identification of estrogen-responsive genes in the parenchyma and fat pad of the bovine mammary gland by microarray analysis. Physiol Genomics. 2006;27(1):42–53.
Hoijman E, Rocha-Viegas L, Kalko SG, Rubinstein N, Morales-Ruiz M, Joffe EB, et al. Glucocorticoid alternative effects on proliferating and differentiated mammary epithelium are associated to opposite regulation of cell-cycle inhibitor expression. J Cell Physiol. 2012;227(4):1721–30.
Forsyth IA. The insulin-like growth factor and epidermal growth factor families in mammary cell growth in ruminants: Action and interaction with hormones. J Dairy Sci. 1996;79(6):1085–96.
Smallwood A, Hon GC, Jin F, Henry RE, Espinosa JM, Ren B. CBX3 regulates efficient RNA processing genome-wide. Genome Res. 2012;22(8):1426–36.
Lemay DG, Pollard KS, Martin WF, Freeman Zadrowski C, Hernandez J, Korf I, et al. From genes to milk: genomic organization and epigenetic regulation of the mammary transcriptome. Plos One. 2013;8(9):e75030.
The authors wish to thank Dr Anne Ridler for her advice and assistance with the mammary biopsies and Meaghan O’Neill for her technical assistance with some of the RT-qPCR experiments. This work was funded by Massey University, Genetics Otago, and Gravida; National Centre for Growth and Development. AP was funded by a Gravida; National Centre for Growth and Development PhD scholarship.
This project was funded through grants from Gravida; the National Centre for Growth and Development to HTB and PKD.
The authors declare that they have no competing interests.
AMP performed the RNA extractions, cDNA synthesis and RT-qPCR experiments with assistance from EJD. AMP, PKD and EJD designed the molecular study, analysed the data and interpreted the results. SJP, SWP, HTP and PRK managed animals. AMP performed tissue collection with assistance from SJP, SWP, HTB, and PRK. SJP, HTP, and PRK sourced funding for these experiments. AMP, PKD and EJD drafted the manuscript. All authors edited and approved the final version of the manuscript.
Availability of data and materials
Summary of RNA-seq reads and mapping statistics for each sample. (DOCX 59 kb)
RNA-seq data for genes more highly expressed during late pregnancy. Table S3: RNA-seq data for genes more highly expressed during lactation. (XLSX 1006 kb)
And supplementary text showing similarity between the RNA-seq samples using principle components analysis (PCA) and heatmaps. (DOCX 243 kb)
Supplementary text detailing RT-qPCR validation of the RNA-seq data. Table S4. Candidate genes measured by RT-qPCR in late pregnant and lactating ovine mammary tissue. Figure S2. Correlation of RNA-seq and RT-qPCR gene expression data in the ovine mammary gland. (DOCX 213 kb)
Gene ontology functional annotation clusters enriched in the ovine mammary gland during late pregnancy. Table S6. Gene ontology functional annotation clusters enriched in the ovine mammary gland during lactation. (XLSX 64 kb)
Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways enriched in the ovine mammary gland during late pregnancy. Table S8: Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways enriched in the ovine mammary gland during lactation. (XLSX 38 kb)
And supplementary text detailing the analysis used to identify clusters of differentially expressed genes in the ovine genome. (DOCX 114 kb)
Comparison of two approaches to identify differentially expressed genes during the transition from late pregnancy to lactation; CLC genomics and EdgeR. (DOCX 137 kb)