Gene regulatory networks in lactation: identification of global principles using bioinformatics
© Lemay et al; licensee BioMed Central Ltd. 2007
Received: 22 August 2007
Accepted: 27 November 2007
Published: 27 November 2007
The molecular events underlying mammary development during pregnancy, lactation, and involution are incompletely understood.
Mammary gland microarray data, cellular localization data, protein-protein interactions, and literature-mined genes were integrated and analyzed using statistics, principal component analysis, gene ontology analysis, pathway analysis, and network analysis to identify global biological principles that govern molecular events during pregnancy, lactation, and involution.
Several key principles were derived: (1) nearly a third of the transcriptome fluctuates to build, run, and disassemble the lactation apparatus; (2) genes encoding the secretory machinery are transcribed prior to lactation; (3) the diversity of the endogenous portion of the milk proteome is derived from fewer than 100 transcripts; (4) while some genes are differentially transcribed near the onset of lactation, the lactation switch is primarily post-transcriptionally mediated; (5) the secretion of materials during lactation occurs not by up-regulation of novel genomic functions, but by widespread transcriptional suppression of functions such as protein degradation and cell-environment communication; (6) the involution switch is primarily transcriptionally mediated; and (7) during early involution, the transcriptional state is partially reverted to the pre-lactation state. A new hypothesis for secretory diminution is suggested – milk production gradually declines because the secretory machinery is not transcriptionally replenished. A comprehensive network of protein interactions during lactation is assembled and new regulatory gene targets are identified. Less than one fifth of the transcriptionally regulated nodes in this lactation network have been previously explored in the context of lactation. Implications for future research in mammary and cancer biology are discussed.
Lactation is one of the most remarkable products of evolution. The signature feature and basis of the competitive emergence of mammals, including humans, is the production of complete early nourishment of neonates by the mother. The processes of lactation include the development of mammary tissue, as well as the synthesis and secretion of milk. At weaning, the mammary gland morphologically returns to a near pre-pregnant state. Thus, in addition to the important nutritional implications, lactation provides a model for basic biological processes such as the proliferation, differentiation, survival and death of cells.
Although lactation is believed to be a product of Darwinian selective pressure, little is known of its molecular origins or its regulation. Current knowledge of the molecular regulation of mammary development and lactation has largely been derived from dissection of signaling networks in cell culture systems and phenotypic characterization of genetically altered mice. Some proteins modulated during pregnancy and lactation have been identified and characterized in the context of hormonal and metabolic pathways (reviewed in [1, 2]). Beyond these signaling pathways, the regulation of mammary gland development and lactation is incompletely understood. Of particular interest are the major molecular events that govern macroscopic and histological changes in the mammary gland during secretory differentiation, secretory activation (the lactation switch), and the onset of involution (the involution switch). Unbiased genome-wide approaches are likely to identify novel genes and gene products involved in the regulation of lactation, particularly when incorporated into a larger picture of mammary development and function.
In this study, bioinformatic techniques are applied to transcriptomic and proteomic data to enhance understanding of how the mammary gland is regulated through pregnancy, lactation, and involution. Using non-hypothesis-driven analyses, transcriptional and post-transcriptional trends are described and putative key regulatory targets are identified. Gene products and their interactions unexplored in the current literature are visualized as a network, providing a framework on which to base future research. Such exploratory methods can be applied to other areas of biological inquiry to establish a quantitative representation of current knowledge and to facilitate the generation of new hypotheses.
Global transcriptional trends during mammary development
Because few genes are strongly correlated with the second and third principal components, other gene lists associated with these trends may be of greater biological interest. For the second component, the 2,064 genes that are statistically up-regulated during involution day 2 (I2) compared to lactation day 9 (L9) are provided as part of Additional data file 5. With respect to the third principal component, the 634 genes that are down-regulated by lactation day 2 (L2) relative to pregnancy day 17 (P17) and up-regulated by I2 relative to L9 are provided in Additional data file 2.
How do these principal components in Figure 1 relate to the known histology (reviewed in [4, 5]) of the mouse mammary gland? In the first part of pregnancy, the epithelial compartment of the gland undergoes a remarkable expansion, from less than 10% of the total volume to about 90% by birth. This proliferative activity is greatest prior to time point pregnancy day 12 (P12). The rise between pregnancy day 1 (P1) and P12 in the third principal component may reflect an increase in expression of epithelial genes due to the increase in this cell population. Secretory differentiation, characterized by an increase in lipid droplets in the cells, begins about mid-pregnancy and continues through P17, at which point the droplets mostly fill the alveolar cells until parturition, pregnancy day 19 (P19), when the lipid droplets move into the lumen representing the first onset of milk secretion. Approximately half of the variance in the microarray data can be described by the first principal component (see Figure 1A) in which gene expression increases from P12 to P19. This suggests that the transcriptional changes necessary to prepare the gland for secretion during this time period are the strongest trend across the entire developmental cycle. Soon after parturition, P19, the lipid droplets disappear from the lumen as milk is removed by the suckling offspring and the gland begins copious milk production. Lactation reaches full maturation by about lactation day 7 and certainly by time point L9. The third principal component implies that the transcription of a subset of genes may be suppressed during lactation.
At weaning, the mammary gland returns to the pre-pregnant state within four to ten days . In the microarray study, the I2 time point may be interpreted as early involution because it occurred two days after the pups were removed. The second principal component shows a sharp upward trend in gene expression between L9 and I2. Very few genes are strongly correlated with this principal component, despite the fact that over 2,000 genes are statistically up-regulated between L9 and I2. This suggests that these genes that are markedly induced during involution are not coordinately regulated during pregnancy or lactation. In other words, involution represents a totally new phase of mammary development that follows a pathway distinct from the pathways that lead to milk secretion.
First principal component, significant biological processes
Positive regulation of interleukin-6 biosynthesis
Positive regulation of calcium-mediated signaling
Positive regulation of cytokine biosynthesis
Positive regulation of lymphocyte proliferation
Positive regulation of transferase activity
Metal ion transport
The genes correlated with the first principal component were also analyzed for enrichment in metabolic and signaling pathways. None of the 152 canonical pathways were statistically significant (p < 0.05) or even marginally significant (unadjusted p < 0.05). It is likely that these genes that are coordinately regulated during mammary development participate in pathways that are not yet annotated.
To assess cumulative transcriptional changes over the full lactation cycle, the number of genes whose expression differs significantly (p < 0.001) from the first time point, P1, is plotted in Figure 2B. The late pregnancy transcriptional state is quite divergent from the initial pregnancy state and the greatest extent of divergence has occurred by full lactation. This graph hints at the possibility of a convergence to the pre-lactation state at the time of early involution, but it does not directly test the number of genes different between the late pregnancy time points and involution. The number of genes whose expression differs significantly (p < 0.001) from the involution time point, I2, is plotted in Figure 2C. The time point with the closest transcriptional state to I2 is the late pregnancy time point, P17. Of the 2,468 genes with statistically different expression in the transition from L9 to I2 (Figure 2A), the expression of 1,688 of these genes is not different between P17 and I2. In other words, 68% of the transcriptional changes between L9 and I2 can be attributed to a partial return to the pre-lactation state. Conceptually, the MECs may be exploring new frontiers of the transcriptional state space to satisfy the unusual requirements of lactation, frontiers from which they must retreat before they can proceed with the involution program.
Transcriptionally regulated biological functions during mammary development
Pregnancy Gene Set, significant biological processes
Golgi vesicle transport
Transition metal ion transport
Protein targeting to ER
ER to Golgi vesicle-mediated transport
Establishment of localization
Inorganic anion transport
Fatty acid beta-oxidation
Fatty acid oxidation
Fatty acid metabolism
Blood vessel morphogenesis
Blood vessel development
Organic acid metabolism
Carboxylic acid metabolism
Lactation Gene Set, significant biological processes
Intracellular protein transport
Establishment of cellular localization
Fatty acid beta-oxidation
Copper ion homeostasis
Cell organization and biogenesis
Establishment of protein localization
Nuclear mRNA splicing, via spliceosome
RNA splicing, via transesterification reactions with bulged adenosine as nucleophile
RNA splicing, via transesterification reactions
Modification-dependent macromolecule catabolism
Modification-dependent protein catabolism
Fatty acid oxidation
Ubiquitin-dependent protein catabolism
Cellular protein metabolism
Cellular macromolecule metabolism
Cellular protein catabolism
Proteolysis during cellular protein catabolism
Antigen presentation, endogenous antigen
Actin cytoskeleton organization and biogenesis
Protein import into nucleus
Cytoskeleton organization and biogenesis
Involution Gene Set, significant biological processes
Cell organization and biogenesis
Establishment of cellular localization
Modification-dependent macromolecule catabolism
Modification-dependent protein catabolism
Intracellular protein transport
Copper ion homeostasis
Ubiquitin-dependent protein catabolism
Establishment of protein localization
Response to oxidative stress
Insulin-like growth factor receptor signaling pathway
Nitrogen compound biosynthesis
Amino acid biosynthesis
Pyrimidine nucleotide metabolism
Pregnancy Gene Set: early pregnancy (P1) to parturition (P19)
By the end of pregnancy, P19, 985 genes are up-regulated and 1,265 genes are down-regulated relative to the start of pregnancy, P1. To screen for transcriptionally modulated biological functions during pregnancy, this Pregnancy Gene Set (see Methods) was analyzed for enriched GO terms. Statistically significant biological process GO terms associated with genes up-regulated during pregnancy are listed in Table 2. Biological processes associated with genes up-regulated during this period are dominated by the synthesis, transport, localization, and secretion of proteins and by transport of ions. MECs appear to be involved in a massive operation of synthesizing and packaging materials before the onset of lactation. Eight of the nine enriched GO terms are related to protein transport, localization, and secretion. This analysis suggests that many of the genes for milk synthesis and secretion are transcribed prior to the onset of lactation.
One of the nutritional properties of milk is the co-delivery of mineral ions in excess of their normal poor bioavailability in aqueous matrices, implying that the processes of lactation must include substantial ion transport capabilities. GO analysis (Table 2) implies that the increased transcription of genes associated with ion transport begins during pregnancy. With respect to ion transport, both zinc and iron ion transport systems are up-regulated. However, some genes associated with iron transport are also down-regulated. Ceruloplasmin (Cp, [GenBank:U49430]), ferritin heavy chain 1 (Fth1, [GenBank:X52561]), lactotransferrin (Ltf, [GenBank:J03298]), sideroflexin 2 (Sfxn2, [GenBank:AA189555]), and transferrin receptor (Tfrc, [GenBank:X57349]) are up-regulated, while ferritin light chain 1 (Ftl1, [GenBank:L39879]), hephaestin (Heph, [GenBank:AF082567]), hemochromatosis (Hfe, [GenBank:Y12650]), sideroflexin 1 (Sfxn1, [GenBank:NM_027324]), and solute carrier family 11 member 1 (Slc11a1, [GenBank:L13732]) are down-regulated. As zinc, copper, and iron are preferentially channeled through the MECs into the milk , this transcriptional modulation of ion transport systems probably reflects this process.
Table 2 also lists statistically significant biological process GO terms associated with genes down-regulated during pregnancy. These GO terms point to three major functions – anion transport, fatty acid oxidation, and blood vessel development. 'Phosphate transport' is the most highly enriched GO term associated with genes down-regulated by late pregnancy. Genes annotated with this GO term are almost all pro-collagen genes, and it is this collagen-forming function that more accurately describes genes up-regulated during pregnancy than 'phosphate transport'. The 'phosphate transport' annotation appears to be out of context when applied to the mammary gland. It was previously suggested that down-regulation of collagen genes is either due to a reduction of the stromal compartment or the expansion of the epithelial compartment in the mammary gland . Table 2 also lists 'fatty oxidation'. It is expected that oxidation of fatty acids would be down-regulated during lactation, and these results show that, consistent with prior study , this function is already down-regulated by late pregnancy.
Surprisingly, genes usually associated with vasculature development appear to be down-regulated, a finding that is not consistent with what is known to be true at the macroscopic level. The gland experiences significant angiogenesis during pregnancy as this process is a necessary prerequisite for lactation. In the rat, the vasculature doubles by mid-pregnancy through angiogenesis by sprouting and intersucception . The 'angiogenesis' GO annotations are likely dictated by cancer research, so there may be genes in the genetic circuit controlling vascularization of normal mammary development that are not yet annotated as such. To facilitate further study, genes associated with these GO terms are provided in Additional files 2, 3, 4. Although some of the biological functions listed in Table 2 are expected given prior study , screening with GO terms indicates the following: some GO terms may generally be out of context when applied to the mammary gland, elements of the secretory pathway are up-regulated prior to parturition, and the vascularization of the normal mammary gland may be incompletely understood at the level of transcription.
Lactation Gene Set: parturition (P19) to mature lactation (L9)
Relative to parturition (P19), 122 genes are up-regulated and 1,704 are down-regulated by the time of mature lactation, L9. This Lactation Gene Set (see Methods) was analyzed for enriched GO terms. Surprisingly, no GO terms are statistically enriched among genes up-regulated during this time period. This leads one to question whether genes up-regulated during lactation are not annotated. However, both up-regulated and down-regulated probes are annotated to approximately the same level, 87 and 89%, respectively. It is possible that the gene set is too small to find statistically enriched terms, or that the genes up-regulated during lactation have disparate functions, or that the common functions have not yet been annotated for their lactation-specific functions.
Enriched GO terms associated with genes down-regulated during lactation relative to parturition are given in Table 3. As during pregnancy, beta-oxidation of fatty acids is suppressed. Because one of the major functions of the mammary gland during lactation is the synthesis and packaging of lipids, it would be expected that oxidative consumption of such lipids for energy by the same cells would be counterproductive. Surprisingly, the transcription of genes involved in protein transport is suppressed during lactation. Thus, the main functions of packaging and exporting proteins by the traditional secretory pathway during the latter part of pregnancy may be transcriptionally attenuated during lactation. As shown in Table 3, 'protein metabolism' – both biosynthesis and catabolism – has also been transcriptionally turned down during lactation. In particular, ubiquitin-dependent protein catabolism is suppressed. Likewise, genes annotated for 'spliceosome assembly' – a ribonucleoprotein apparatus that catalyzes nuclear mRNA splicing – are also down-regulated, implying that both mRNA and protein processing are essentially halted or dramatically slowed during lactation. However, synthesis of the major milk proteins is known to be highly up-regulated . Thus, it is likely that a small subset of genes is specifically utilized by the lactating gland to process proteins for secretion.
Genes associated with 'antigen processing' or 'antigen presentation' are also statistically down-regulated during lactation relative to late pregnancy. Although the complete immune function of colostrum may not be known, it is well-accepted that colostrum, which would be expected to be functionally tied to the genes transcribed in late pregnancy, has substantial immune components. Perhaps these genes, or a subset of them, are unique to colostrum and thus, their transcription is switched off by the time of full lactation. In summary, the lactation period appears to be characterized by widespread suppression of transcription across many functional classes.
Involution Gene Set: mature lactation (L9) to involution (I2)
By involution, time point I2, 2,064 genes are up-regulated and 404 genes are down-regulated relative to mature lactation, L9. GO analysis of this Involution Gene Set (see Methods) yields enriched GO terms for both up-regulated and down-regulated genes (Table 4). Seventeen of the twenty enriched GO terms (85%) associated with up-regulation during involution exactly match GO terms of the genes down-regulated during lactation. Transcriptionally, the switch to involution largely reverses the state of lactation. Genes involved in protein transport are transcribed again, re-enabling the traditional secretory pathway, for example. Protein catabolism is up-regulated, an expected major function during involution of the mammary gland. The 'response to oxidative stress' GO term is probably associated with the massive protein catabolism that is taking place, either to protect the MECs or the surrounding cells. Genes involved in metal ion homeostasis, especially copper, are also transcribed again.
Proton transport is up-regulated during involution. Because both ATP synthesis-coupled and ATP hydrolysis-coupled proton transport are up-regulated, the up-regulation of these genes is probably not reflective of an energetic transition of the mammary cells. It is possible that this set of genes is up-regulated to facilitate the degradation of milk components and the transfer of the small molecules to the interstitial compartment.
Far fewer genes, only 404, are down-regulated during involution. Nevertheless, a few clustered biological functions are represented (Table 4). In general, biosynthetic and metabolic processes are down-regulated, as would be expected of involuting cells. Overall, early involution appears to be characterized by a suppression of specific lactation functions and up-regulation of genes associated with novel functional categories.
Transcriptionally regulated pathways during mammary development
To identify which of the known metabolic and signaling pathways are transcriptionally regulated, pathway analysis was applied to the Pregnancy, Lactation, and Involution Gene Sets (see Methods). Of the 152 canonical pathways evaluated, 16 are marginally significant during pregnancy, 22 during lactation, and 40 during involution (unadjusted p < 0.05). When a multiple testing correction is applied (see Methods), no pathways are significant during pregnancy, 8 are significant during lactation, and 24 are significant during involution (p < 0.05). The complete list of pathways, their p-values, and the molecules regulated in those pathways are provided (Additional data file 19), as well as color-coded diagrams for the significant pathways (Additional data file 20, 21, 22). An interpretation of all of these pathways is well beyond the scope of this paper. Therefore, the three pathways that most significant during lactation are briefly explored in the following sections.
The PI3K-AKT pathway is highly significant during lactation (p = 0.0101) and involution (p < 0.0001) (Figure 3B–C). Thus, the pathway analysis suggests that Akt is a central regulator throughout the lactation cycle. Also, Akt1 appears to be the sole down-regulated protein in the pathway during involution, implying a critical regulatory role during this developmental stage. Experimental data supports these findings. Akt expression increases dramatically at the onset of lactation but is relatively low during involution . In transgenic mice with constitutively active Akt1, involution is delayed, suggesting that Akt1 is an essential regulator of involution [13–15]. Interestingly, the overexpression of Akt1 in these mice during the pregnancy and lactation stages pathologically enhanced lipogenesis . The fat content of the milk was 65–70%, compared to 25–30% in wild-type mice. Thus, experimental evidence supports Akt1 as a central regulator during all phases of lactation. In Figure 3, genes that are transcribed to form products upstream and downstream of Akt1 are also significantly transcriptionally regulated. The enrichment of the entire PI3K-AKT pathway for differentially regulated genes in the three time periods analyzed in this paper suggests that the entire pathway, and not solely Akt1, is coordinately regulated throughout the lactation cycle.
Protein ubiquitination pathway
Biological principles of the lactation proteome
Interestingly, more secreted proteins (transcripts with signal peptides) are down-regulated during pregnancy and lactation than are up-regulated. In fact, fewer than 100 transcripts with signal peptides are up-regulated during pregnancy or lactation. This suggests that the secreted proteins are limited to those important to milk composition and that other secretory functions are suppressed as the gland single mindedly pursues its lactational course. Thus, it is quite possible that the secretory pathway components serve primarily to synthesize and secrete large volumes of only a few proteins.
Protein interactions of genes expressed in the mammary gland
Networks with few highly-connected nodes of high degree (i.e., with many links) and with most nodes of low degree (i.e., with few links) are considered to be scale-free. The networks in Figure 8A, C, and 8D–F all have degree distributions with similar shapes that appear to follow a power law, as would be expected of scale-free biological networks (data not shown). The degree distribution of the network in Figure 8B is, trivially, a single point.
From the network in Figure 8C, it is clear that genes up-regulated during involution do directly interact with each other at the protein level and are highly connected. Genes up-regulated by late pregnancy interact with each other to a limited extent (Figure 8A), either because there are fewer protein-protein interactions or because the protein interactions in which they are involved in mammary cells have not been observed in the high-throughput biology experiments used to derive protein interaction data. Few genes are up-regulated during lactation, and here (Figure 8B) there is no evidence that they directly interact with each other at the protein level. Again, some protein interactions may occur that are unique to mammary cells and would therefore not be present in these networks derived from OPHID.
The star patterns in Figure 8E suggest that a couple of the genes up-regulated during lactation directly interact with many other proteins. Network analysis indicates that these four lactation proteins are the most highly connected: KH domain-containing, RNA-binding, signal transduction-associated protein 1 (SAM68, [PIR:Q07666]), guanine nucleotide-binding protein subunit beta 4 (GBB4, [PIR:Q9HAV0]), C-jun-amino-terminal kinase-interacting protein 1 (JIP1, [PIR:Q9UQF2]), and insulin-like growth factor IB precursor (IGF1B, [PIR:P05019]). SAM68, GBB4, JIP1, and IGF1B interact with 35, 16, 12, and 11 proteins, respectively. Annotations in the UniProt database  provide possible functions for these proteins. SAM68 is an RNA binding protein that is involved in mRNA processing. GBB4 is a G protein likely involved as a modulator or transducer in various trans-membrane signaling systems. JIP1 may function as a regulator of vesicle transport, through interactions with the JNK-signaling components and motor proteins or as an anti-apoptotic protein. IGF1B is an insulin-like growth factor. In summary, the networks in Figure 8 suggest that the extent of interactions between proteins that are transcriptionally regulated may be a function of reproductive stage, with very few interactions occurring during lactation.
Development of a lactation protein interaction network
To develop a more comprehensive lactation protein-protein interaction network, the scientific literature was mined for genes involved in lactation. Three gene sets – the Literature Gene Set (see Methods), up-regulated genes from the Pregnancy Gene Set, and up-regulated genes from the Lactation Gene Set – were combined and used to query the protein interaction database. Restricting the network only to interactions between proteins in the query list, this network contains 313 nodes and 438 linkages. Approximately 71% of the nodes belong to a single major component (subnetwork). The remaining nodes are associated with subnetworks consisting of 6 or fewer nodes each.
To understand what known functions are associated with this comprehensive lactation network, GO analysis was applied to the proteins in this network relative to all human UniProt/SwissProt Accessions available through Ensembl . The top three non-synonymous GO terms enriched in the lactation network (p < 0.001) are 'JAK-STAT cascade,' 'anti-apoptosis,' and 'COPI coated vesicle membrane'. All might have been predicted from previous understanding of genes regulated during lactation.
Pathway analysis was also applied to the proteins in this lactation network. The pathways most significantly enriched (p < 0.005) with these proteins are the 'neuregulin', 'JAK/Stat', 'NF-kappaB', 'PI3K/AKT', 'IL-4', 'PTEN', and 'IGF-1' signaling pathways. In total, 38 pathways are marginally statistically enriched (unadjusted p < 0.05), 22 of which are significantly enriched (p < 0.05). However, no more than 20% of any one of these significant pathways is represented by nodes from the lactation network. These facts together suggest that proteins in this network are interacting in a way that is not well-characterized by any individual canonical signaling or metabolic pathways of this set of 152 curated pathways.
Major hubs of the comprehensive lactation network
No. of Links
UniProt Accession/HGNC Symbol
Epidermal growth factor receptor precursor
Ornithine aminotransferase, mitochondrial
Tyrosine-protein kinase JAK2
Receptor tyrosine-protein kinase erbB-2 precursor
Signal transducer and activator of transcription 3
Vacuolar ATP synthase subunit E
Breast cancer type 1 susceptibility protein
Signal transducer and activator of transcription 5A
6-phosphogluconate dehydrogenase, decarboxylating
Coatomer subunit beta
Growth hormone receptor precursor
Tyrosine-protein phosphatase non-receptor type 6
Eukaryotic translation initiation factor 4E
Seryl-tRNA synthetase, cytoplasmic
Protein transport protein Sec23B
Histone acetyltransferase type B catalytic subunit
Suppressor of cytokine signaling 3
Insulin receptor substrate 1
Suppressor of cytokine signaling 1
Mothers against decapentaplegic homolog 2
Ubiquitin-fold modifier 1 precursor
Insulin receptor precursor
Eukaryotic translation initiation factor 3 subunit 4
CCAAT/enhancer-binding protein beta
The identification of all genes responsible for lactation was previously proposed . However, the lactome is spatially and temporally complex. For successful lactation to occur, molecular events outside of the mammary gland would be expected to support maternal-offspring behaviors, offspring sucking and rooting reflexes, offspring digestive response to the milk, and so forth. If one limits the spatial scope to the mammary gland, this study demonstrates temporal complexity as well. Even with a very strict Bonferonni multiple testing correction (p < 0.001), nearly a third of the mouse transcriptome is in flux from pregnancy through involution. It is possible that some of these changes are due to the expansion and contraction of the epithelial compartment and concurrent proportional representation of the stromal cells. Ultimately, housekeeping genes of these changing cell populations may need to be separated from genes that are uniquely regulated in the context of the lactation apparatus. Despite these obvious limitations, an initial list of genes that comprise the lactome – those differentially regulated by late pregnancy relative to early pregnancy and by full lactation relative to late pregnancy – are provided with this study (see Pregnancy and Lactation Gene Sets in Methods).
Annotation of the lactome
Gene ontology annotations are not necessarily applicable to the lactome as these annotations have arisen from studies in culture cells of other tissues. Here, the use of GO annotations provides a method to screen for unexpected biological functions (i.e., down-regulation of known angiogenesis-related genes during pregnancy) and to identify trends (i.e., reciprocal modulation of gene classes between lactation and involution). Annotation with respect to the specific biological context of lactation will ultimately be required. The comprehensive lactation protein interaction network derived in this paper provides a working model for future work on molecular interactions in mammary cells and highlights gene targets that have not yet been annotated for lactation in the literature. Only 17% of the regulated gene products in the network (Figure 9) co-occur in the literature with the keyword, 'lactation'. Even among the most central regulators – the hubs of the network in Table 5 – eleven of the 27 (40.7%) do not yet have a described role with respect to lactation or mammary development.
The lactation switch
The transcriptional profiles of genes by functional class (i.e., milk proteins, lipid synthesis, etc.) have been well-described by a previous study , and the molecular events surrounding secretory activation have been recently reviewed . This study illustrates, on a broader perspective, very gradual transcriptional changes with progression of reproductive stage until involution (Figure 2A). In other words, aside from small subsets of proteins such as the major milk proteins, there is no sudden transcriptional switch around the time of parturition. Preparations of the mammary gland for lactation include modifications to the transcriptional program, but the onset of lactation appears to be primarily controlled by post-transcriptional mechanisms.
Regulation of proteins upstream of mRNA transcription during lactation is not a new concept. In the presence of lactogenic hormones, the mRNA of casein genes accumulates rapidly due to increased RNA stability . More recently, enhancement of translation of the mRNA of the beta-casein protein was found to occur by lengthening the poly(A) tract via cytoplasmic polyadenylation synergized by the activity of both prolactin and insulin signaling . The rate of translation has also been found to be reduced by amino acid deprivation .
This study provides further targets for exploration of post-transcriptional regulation in the mammary gland. Of the mere 82 genes up-regulated during lactation relative to late pregnancy, the gene with the most known interactions at the protein level – the hub of the network with the highest number of links – is SAM68, an RNA binding signal transduction protein. SAM68, also known as KHDRBS1, is a putative regulator of mRNA splicing, translation, and nuclear export [27–29] and has recently been shown to play a role in regulation of apoptotic genes . SAM68, therefore, has the potential to be an important post-transcriptional regulator of both milk secretion and mammary cell survival during lactation. Other potential post-transcriptional mediators of lactation exist among the hubs of the comprehensive lactation network: the eukaryotic translation initiation factors 4E (IF4E) and 3 subunit 4 (IF34). Long et al. hypothesized that IF4E may effect translation of certain mRNA or the acceleration of overall protein synthesis . The existence of IF4E as a hub in the lactation network certainly supports a central regulatory role.
Network analysis in this study demonstrates that the proteins encoded by genes transcriptionally up-regulated in the mammary gland during lactation do not interact with each other. Thus, regulation in MECs may substantially occur through protein-protein interactions of gene products that pre-exist in the cell. To some extent, the differences seen in network connectivity, density, and shape may be an artifact of gene set size. There are roughly twice as many genes up-regulated during involution as there are during pregnancy and more than ten times as many compared to lactation. Also, because protein interaction data is derived by high-throughput screens of protein interactions in other cell systems, it is possible that unique protein interactions exist in the mammary gland have not yet been observed.
Mammary epithelial cells during lactation
During lactation, MECs are essentially biofactories of lipids, proteins, and carbohydrates for milk. The transcriptional trends highlighted in this study imply that these cells become biofactories not by gain of function, but by a broad suppression of function to effectively push all of the cell's resources towards a very few important tasks such as the massive synthesis of lipids and of a large quantity of a minority of proteins. During lactation, MECs are highly specialized, inhabiting an extreme transcriptional state space that is the most divergent from the early pregnant state. The near-universal attenuation of integrin signaling suggests that MECs isolate themselves from their environment while they carry out their singular agenda. The striking suppression of the proteolytic machinery in both the GO and pathway analysis, particularly the ubiquitination pathway, suggests that an early step in mammary evolution was the acquisition of tolerance for such quantities of protein structures beyond that normally acceptable within cells. Such a suppression of proteolysis may have been necessary to allow the emergence of nutritionally functional proteins, particularly those with little stabilizing secondary and tertiary structure. MECs during lactation are truly unique and could serve as useful model for biological questions that transcend mammary biology.
Secretory diminution is the gradual decline in milk production after peak lactation has been attained. The fact that the machinery for milk synthesis appears to be in place by the time of late pregnancy lends evidence towards a new theory of secretory diminution. With the genes contributing to the machinery down-regulated during lactation, the machinery is not replenished as it degrades over time. Thus, secretory diminution is neither programmed nor pathological as recently reviewed ; rather, it is an inevitable feature of this biological design.
The involution switch
Using microarray technology, transcriptional profiles during involution have been studied by numerous groups [7, 8, 32, 33]. Here, analysis of the transcriptome more clearly demonstrates that a massive transcriptional switch occurs between lactation and early involution (Figure 2A). Five times as many transcripts are differentially regulated at weaning relative to other times in the development of the gland from the onset of the pregnancy. Conceptually, a genome-wide transcriptionally-driven switch at involution may be biologically necessary because the timing of weaning is determined by the offspring and therefore, is unpredictable. Thus, the mammary gland cannot make gradual preparations for its onset. The signals responsible for this event are likely to be of interest for many biological processes particularly related to tissue remodeling.
The analysis in this study implies that transcriptional changes during early involution may function, in part, to revert the mammary gland to the pre-lactation state. First, the state space of the transcriptome gradually diverges from the onset of pregnancy through lactation and then, at involution, the transcriptome partially converges to the pre-lactation state (Figure 2B). More than half of the transcriptional changes during the transition from lactation to early involution bring the gland back to the transcriptional state of late-pregnancy. Second, GO analysis demonstrated that functionally, the transcriptional state is largely a reversion to the pre-lactation state – 85% of the enriched GO terms associated with up-regulation during involution exactly match GO terms of the genes down-regulated during lactation. To study involution, prior syntheses of microarray data have been applied exclusively to involution time points . However, comparison of these involution time points to transcriptional changes during pregnancy should improve the ability to sort the many signals that are multiplexed during this very complex time period.
The transcriptional state space of a full lactation cycle was previously explored by Master et al using a three-dimensional visualization of the dimensionality-reduced pairwise distances between developmental time points . The observation in this study that the state space of the transcriptome gradually diverges from the onset of pregnancy through lactation is consistent with their model. However, the observed partial convergence of the early involution transcriptional state towards the pre-lactation state is not suggested by their visualization. There are a couple of possibilities for this discrepancy. It is possible that the early involution time point is much closer to the near-parturition time point in their model than appears on the page in two dimensions. Second, the shape of the entire trajectory may be overly constrained in low dimensions such that the goodness-of-fit between any two of the developmental time points is compromised. In this study, for the particular question of whether the early involution transcriptional state is a partial reconvergence to the pre-lactation state, the pairwise differences are directly measured so the observed distances are unequivocal.
The network of interactions between proteins that are encoded from genes up-regulated during early involution is well-connected. This suggests that the process of early involution is coordinately regulated. That is, genes up-regulated together produce proteins that interact with each other. Furthermore, the cellular processes are unlikely to be unique to involution of the mammary gland because they are described by protein interaction data derived from other cell cultures. The protein interaction network at the time of involution (Additional data file 16) provides a roadmap for future molecular studies of involution.
The milk proteome
GO analysis in this study suggests that much of the machinery for the secretory pathway is transcribed prior to lactation. Because very few genes are up-regulated during lactation relative to late pregnancy, this secretory pathway is likely devoted to the secretion of large amounts of a few proteins. Highly expressed milk protein genes, such as the caseins, account for as much as 30% of total RNA (personal communication, MC Neville). As a group, genes transcribed to form these major milk proteins are up-regulated sharply around the time of parturition.
Assessments of the number of unique proteins in the milk proteome differ due to technical limitations in available proteomic technologies. In this study, fewer than 100 genes whose proteins are destined for secretion were up-regulated during either pregnancy or lactation. This suggests that the endogenous portion of the milk proteome that is synthesized through the secretory pathway is derived from fewer than 100 transcripts. Some proteins such as immunoglobins, transferrin, and albumin enter the milk by the transcytosis pathway [35–37]. Thus, if the diversity of the milk proteome proves to be substantially greater than 100 proteins, that diversity is due to either to post-translational modifications of these 100 transcripts or to the import by the gland of many more exogeneous proteins than previously realized.
The small number of up-regulated transcripts with signal peptides also suggests that secretory functions not related to lactation are turned off while only those important to milk composition are retained. This biological principle is consistent with the view that the MECs have reached a very specialized transcriptional state that is single mindedly devoted to milk secretion.
Implications for cancer research
During pregnancy, lactation, and involution, nearly a third of the mouse transcriptome is in flux within the mammary gland during which the proliferation, differentiation, and death of cells are exquisitely controlled. The survival of mammalian species depends on the ability of these cells to traverse this vast transcriptional state space repeatedly without malignancy. As such, normal mammary development provides a rich model system for the study of cellular development and the effects of perturbation.
Intriguingly, GO analysis revealed that genes known to be involved in angiogenesis in other tissues are, as a class, down-regulated during pregnancy. On a macroscopic level, the vascularization of the gland clearly increases. In the rat, the vasculature doubles by mid-pregnancy . The 'angiogenesis' GO annotations are likely dictated by cancer research, so there may be genes that are important to angiogenesis in normal mammary development that are not yet annotated as such. Angiogenesis within the normal developing mammary gland may fundamentally differ from that of the neoplastic transformation. Other evidence of this hypothesis is provided by HIF1alpha. In human tumors, HIF1alpha was up-regulated, particularly in breast tumors that exhibit high rates of proliferation[38, 39]. However, in the HIF1alpha null mouse mammary gland, vasculature development was unchanged . Taken together, these studies imply that the angiogenesis of breast cancer is not merely the untimely enabling of the normal mammary vascularization circuit.
The techniques and data used in this study have several limitations worth noting. First, principal component analysis reduces the dimensionality of a data set by retaining components of the data that contribute most to its variance. These lower order components are generally the most important; however, it is possible that these components are unimportant or that the discarded higher order components are biologically relevant. Second, the application of GO analysis is unlikely to completely describe the function of genes transcribed during mammary development as many genes have not yet been annotated in this biologically unusual context. Third, the statistical tests commonly available to assess the significance of enrichment in the GO and pathway analyses are anti-conservative and based on flawed assumptions . A multiple testing correction was used to mitigate the former, but solutions for the latter are not readily available. Also, enrichment of small gene sets, such as the set of genes that form products of the JAK-Stat pathway, are overly sensitive to small perturbations in the gene set under test, a property that will be inherent to any statistical test of enrichment. Fourth, the network analysis in this study is limited by the integrity of the underlying data. Protein interaction data, derived largely from large-scale yeast two-hybrid experiments, is known to contain both false positives and missed interactions . It is also highly likely that protein interactions important to mammary development are incompletely represented by experiments with other cell types. Lastly, the microarray chips utilized were limited to 12,488 probes so it is possible that some differentially regulated transcripts were not measured.
Bioinformatic techniques were applied to various data sets to derive global biological principles governing the molecular events in the mammary gland during pregnancy, lactation, and involution. The key biological principles identified include the following: (1) nearly a third of the transcriptome fluctuates to build, run, and disassemble the lactation apparatus; (2) genes encoding the secretory machinery are transcribed prior to lactation; (3) the diversity of the endogenous portion of the milk proteome is derived from fewer than 100 transcripts; (4) while some genes are differentially transcribed near the onset of lactation, the lactation switch is primarily post-transcriptionally mediated; (5) the secretion of materials during lactation occurs not by up-regulation of novel genomic functions, but by widespread transcriptional suppression of functions such as protein degradation and cell-environment communication; (6) the involution switch is primarily transcriptionally mediated; and (7) during early involution, the transcriptional state is partially reverted to the pre-lactation state. These guiding principles begin to establish a road map for future study of mammary development during pregnancy, lactation, and involution.
The seeds for new hypotheses are also provided in this study. First, secretory diminution, usually attributed to programmed cell death or oxidative damage, may actually be due to lack of replenishment of the secretory machinery at the mRNA level. GO analysis in this analysis suggests that once lactation has commenced, transcription of RNA to maintain the secretory machinery is not sustained. It is, therefore, inevitable that milk production will gradually decrease over time. Second, angiogenesis-related genes are down-regulated during pregnancy, while vasculature in the gland clearly increases. This suggests that the angiogenesis of breast cancer is fundamentally pathological, rather than an inappropriate launch of normal vascularization processes.
In addition to the biological principles derived in this paper, extensive supplementary materials are provided as a mineable resource for mammary biologists. The 22 additional data files include sets of differentially expressed genes, associated GO terms, libraries of annotated pathway diagrams, and developmental stage-specific protein interaction networks. The comprehensive lactation network of protein interactions assembled in this study prioritizes regulatory gene targets for future study and provides a framework from which to base future research. The network enables the visualization of connectivity between those gene products known to be involved in lactation and those previously unexplored in the lactation literature. These integrative methods can be applied to other areas of biological inquiry to establish a snapshot of current knowledge and to facilitate the generation of new hypotheses.
Microarray study design
This microarray study was previously described in Rudolph et al. . Briefly, mammary gland RNA samples were collected from FVB mice, isolated, and hybridized onto Affymetrix MG_U74Av2 chips. In total, ten time points in mammary development with four biological replicates for each time point were produced: P1, the day a vaginal plug was observed; P3, pregnancy day 3; P7, pregnancy day 7; P12, pregnancy day 12; P17, pregnancy day 17; P19, pregnancy day 19; L1, early lactation, the first day pups are observed in the cage, L2, lactation day 2; L9, lactation day 9; and I2, two days after pup removal on lactation day 9. Histologically, the mammary proliferative stage is represented by P1, P3, and P7, the secretory differentiation stage by P12, P17, and P19, early lactation by L1 and L2, full lactation by L9, and involution by I2. These data have been deposited in NCBI's Gene Expression Omnibus  and are accessible through GEO Series accession number GSE8191.
Microarray data analysis
GeneSpring GX 7.3.1 was used to analyze the data. First, GC-RMA preprocessing was applied to all CEL files. Signal intensity values were normalized as follows. Values below 0.01 were set to 0.01. Each measurement was divided by the 50.0th percentile of all measurements in that sample. Each gene was divided by the median of its measurements in all samples. If the median of the raw values was below 10, then each measurement for that gene was divided by 10 if the numerator was above 10, otherwise the measurement was discarded.
The 12,488 genes were then pre-filtered to remove those genes with unreliable or undetectable signals. Those that did not have an Affymetrix call of 'Present' in at least four of the forty samples were discarded. Additionally, genes with consistently low expressing signal across all ten developmental stages were removed using a minimum threshold of 0.712, empirically decided using GeneSpring's Cross Gene Error Model. Genes that did not exceed this threshold in at least one of the ten developmental stages were removed. The remaining 7,534 genes were used for subsequent analysis.
To find genes that were significantly differentially expressed during different developmental stages of the mammary gland, from initial pregnancy to involution, a one-way ANOVA was applied with a Bonferroni multiple testing correction, p-value cutoff of 0.001. No genes would be expected to pass this restriction by chance. 4,832 genes passed. To find genes significantly differentially expressed between individual development stages, all pairwise comparisons were examined with a Tukey post hoc test (p < 0.001).
Pregnancy, lactation, and involution gene sets
The Pregnancy Gene Set, listed in Additional data file 3, is the cluster of genes that are significantly differentially regulated between the start of pregnancy (time point P1) and late pregnancy (time point P19). The Lactation Gene Set, Additional data file 4, is the group of genes significantly differentially regulated between late pregnancy (time point P19) and full lactation (time point L9). The Involution Gene Set, Additional data file 5, refers to those genes significantly differentially regulated between full lactation (time point L9) and involution (time point I2). The Tukey post hoc test was applied to all pairwise comparisons to determine significance (p < 0.001).
Literature Gene Set
The Single Gene Biological Term Mapper feature of CoPub Mapper  was used to extract all of the genes co-occuring in PubMed publications with the biological process keyword 'lactation'. In total, 685 genes have been co-published with 'lactation' two or more times. Of these genes, 456 co-occur with a relative score greater than zero. The relative score, previously described , is a measure of the frequency of co-occurrence adjusted for the frequency of occurrence of each item individually. This set of 456 genes – those genes known to be involved in lactation – comprise the Literature Gene Set, provided in Additional data file 6.
Gene ontology analysis
MAPPFinder 2.0 within GenMAPP 2.0  was used to identify gene ontology (GO) terms that are over-represented by significantly regulated genes compared to the full genome-wide set of Affymetrix ProbeIDs given the most recent mouse gene database, Mm-Std_20060628.gdb. Custom Perl and shell scripts were written to format the data for analysis. To determine statistical significance, a Benjamini and Hochberg multiple testing correction was applied to GenMAPP's permuted p values with a custom R script using the multtest library [46, 47]. In each analysis, statistically enriched GO terms are those terms with 3 or more associated genes, a positive z-score, and a Benjamini and Hochberg adjusted p-value less than or equal to 0.05. Thus, a False Discovery Rate of 5% or less would be expected for these enriched GO terms. These same statistical tests were applied to determine enriched GO terms among genes correlated with the principal components and among gene products represented in the lactation network.
Ingenuity Pathways Analysis  was used to identify metabolic and signaling pathways that are over-represented by significantly regulated genes compared to all of the genes that are part of the Ingenuity Pathways Knowledge Base. The Ingenuity Pathways Analysis library of canonical pathways includes 80 metabolic and 72 signaling pathways that have been incorporated from various resources and hand-curated. Genes that were statistically up-regulated or down-regulated between the time points studied (i.e. P1 to P19, P19 to Lac9, Lac9 to Inv2) and were associated with a canonical pathway were considered for the analysis. A Fischer's exact test was used to calculate a p-value to determine the probability that the enrichment of the canonical pathway with these genes is explained by chance alone. In this paper, this p-value is referred to as the unadjusted p-value. To improve the stringency of the test, a Benjamini and Hochberg multiple testing correction was applied to the unadjusted p values with a custom R script using the multtest library [46, 47]. Unless otherwise stated, pathways reported to be statistically significant are those with a Benjamini and Hochberg adjusted p-value less than or equal to 0.05. These same statistical tests were applied to determine enriched pathways among genes correlated with the principal components and among gene products represented in the lactation network.
Protein interaction network analysis
Gene lists were converted to Human Swissprot IDs using tables from the Ensembl database, release 41 . For each list of Human Swissprot IDs, interactions between those gene products were obtained from OPHID  and post-processed using custom scripts to determine all of the linkages in the network and to generate a network file. This network file was then explored using Pajek 1.16 , a program for large network analysis.
List of abbreviations
involution day 2
lactation day 2
lactation day 9
mammary epithelial cells
Online Predicted Human Interaction Database
pregnancy day 1
pregnancy day 12
pregnancy day 17
pregnancy day 19.
The authors acknowledge CJ Dillard for editorial assistance with the manuscript, AS Hinrichs for technical assistance with gnuplot, and Agilent Technologies for providing the GeneSpring GX software. This research was supported in part by the California Dairy Research Foundation, UC Discovery, NIEHS Superfund grant P42 ES04699, and by the CHARGE study (grant P01 ES11269). The acquisition of the microarray data set was supported by NIH grant P01-HD38129 to MCN.
- Hennighausen L, Robinson GW: Information networks in the mammary gland. Nat Rev Mol Cell Biol. 2005, 6 (9): 715-725. 10.1038/nrm1714PubMedView ArticleGoogle Scholar
- Anderson SM, Rudolph MC, McManaman JL, Neville MC: Key stages in mammary gland development. Secretory activation in the mammary gland: it's not just about milk protein synthesis!. Breast Cancer Res. 2007, 9 (1): 204- 10.1186/bcr1653PubMed CentralPubMedView ArticleGoogle Scholar
- 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. 10.1023/B:JOMG.0000010030.73983.57PubMedView ArticleGoogle Scholar
- Palmer CA, Neville MC, Anderson SM, McManaman JL: Analysis of lactation defects in transgenic mice. J Mammary Gland Biol Neoplasia. 2006, 11 (3-4): 269-282. 10.1007/s10911-006-9023-3PubMedView ArticleGoogle Scholar
- Richert MM, Schwertfeger KL, Ryder JW, Anderson SM: An atlas of mouse mammary gland development. J Mammary Gland Biol Neoplasia. 2000, 5 (2): 227-241. 10.1023/A:1026499523505PubMedView ArticleGoogle Scholar
- Atabai K, Sheppard D, Werb Z: Roles of the innate immune system in mammary gland remodeling during involution. J Mammary Gland Biol Neoplasia. 2007, 12 (1): 37-45. 10.1007/s10911-007-9036-6PubMed CentralPubMedView ArticleGoogle Scholar
- Master SR, Hartman JL, D'Cruz CM, Moody SE, Keiper EA, Ha SI, Cox JD, Belka GK, Chodosh LA: Functional microarray analysis of mammary organogenesis reveals a developmental role in adaptive thermogenesis. Mol Endocrinol. 2002, 16 (6): 1185-1203. 10.1210/me.16.6.1185PubMedView ArticleGoogle Scholar
- Clarkson RW, Wayland MT, Lee J, Freeman T, Watson CJ: Gene expression profiling of mammary gland development reveals putative roles for death receptors and immune mediators in post-lactational regression. Breast Cancer Res. 2004, 6 (2): R92-109. 10.1186/bcr754PubMed CentralPubMedView ArticleGoogle Scholar
- Kelleher SL, Lonnerdal B: Molecular regulation of milk trace mineral homeostasis. Mol Aspects Med. 2005, 26 (4-5): 328-339. 10.1016/j.mam.2005.07.005PubMedView ArticleGoogle Scholar
- Djonov V, Andres AC, Ziemiecki A: Vascular remodelling during the normal and malignant life cycle of the mammary gland. Microsc Res Tech. 2001, 52 (2): 182-189. 10.1002/1097-0029(20010115)52:2<182::AID-JEMT1004>3.0.CO;2-MPubMedView ArticleGoogle Scholar
- Yang ZZ, Tschopp O, Baudry A, Dummler B, Hynx D, Hemmings BA: Physiological functions of protein kinase B/Akt. Biochem Soc Trans. 2004, 32 (Pt 2): 350-354. 10.1042/BST0320350PubMedView ArticleGoogle Scholar
- Schwertfeger KL, McManaman JL, Palmer CA, Neville MC, Anderson SM: Expression of constitutively activated Akt in the mammary gland leads to excess lipid synthesis during pregnancy and lactation. J Lipid Res. 2003, 44 (6): 1100-1112. 10.1194/jlr.M300045-JLR200PubMedView ArticleGoogle Scholar
- Schwertfeger KL, Richert MM, Anderson SM: Mammary gland involution is delayed by activated Akt in transgenic mice. Mol Endocrinol. 2001, 15 (6): 867-881. 10.1210/me.15.6.867PubMedView ArticleGoogle Scholar
- Hutchinson J, Jin J, Cardiff RD, Woodgett JR, Muller WJ: Activation of Akt (protein kinase B) in mammary epithelium provides a critical cell survival signal required for tumor progression. Mol Cell Biol. 2001, 21 (6): 2203-2212. 10.1128/MCB.21.6.2203-2212.2001PubMed CentralPubMedView ArticleGoogle Scholar
- Ackler S, Ahmad S, Tobias C, Johnson MD, Glazer RI: Delayed mammary gland involution in MMTV-AKT1 transgenic mice. Oncogene. 2002, 21 (2): 198-206. 10.1038/sj.onc.1205052PubMedView ArticleGoogle Scholar
- Kass L, Erler JT, Dembo M, Weaver VM: Mammary epithelial cell: Influence of extracellular matrix composition and organization during development and tumorigenesis. Int J Biochem Cell Biol. 2007, 39 (11): 1987-1994. 10.1016/j.biocel.2007.06.025PubMed CentralPubMedView ArticleGoogle Scholar
- McMahon CD, Farr VC, Singh K, Wheeler TT, Davis SR: Decreased expression of beta1-integrin and focal adhesion kinase in epithelial cells may initiate involution of mammary glands. J Cell Physiol. 2004, 200 (2): 318-325. 10.1002/jcp.20011PubMedView ArticleGoogle Scholar
- Fink JL, Aturaliya RN, Davis MJ, Zhang F, Hanson K, Teasdale MS, Kai C, Kawai J, Carninci P, Hayashizaki Y, Teasdale RD: LOCATE: a mouse protein subcellular localization database. Nucleic Acids Res. 2006, 34 (Database issue): D213-7. 10.1093/nar/gkj069PubMed CentralPubMedView ArticleGoogle Scholar
- Brown KR, Jurisica I: Online predicted human interaction database. Bioinformatics. 2005, 21 (9): 2076-2082. 10.1093/bioinformatics/bti273PubMedView ArticleGoogle Scholar
- The Universal Protein Resource (UniProt). Nucleic Acids Res. 2007, 35 (Database issue): D193-7.
- Hubbard TJ, Aken BL, Beal K, Ballester B, Caccamo M, Chen Y, Clarke L, Coates G, Cunningham F, Cutts T, Down T, Dyer SC, Fitzgerald S, Fernandez-Banet J, Graf S, Haider S, Hammond M, Herrero J, Holland R, Howe K, Howe K, Johnson N, Kahari A, Keefe D, Kokocinski F, Kulesha E, Lawson D, Longden I, Melsopp C, Megy K, Meidl P, Ouverdin B, Parker A, Prlic A, Rice S, Rios D, Schuster M, Sealy I, Severin J, Slater G, Smedley D, Spudich G, Trevanion S, Vilella A, Vogel J, White S, Wood M, Cox T, Curwen V, Durbin R, Fernandez-Suarez XM, Flicek P, Kasprzyk A, Proctor G, Searle S, Smith J, Ureta-Vidal A, Birney E: Ensembl 2007. Nucleic Acids Res. 2007, 35 (Database issue): D610-7. 10.1093/nar/gkl996PubMed CentralPubMedView ArticleGoogle Scholar
- Grigorov MG: Global properties of biological networks. Drug Discov Today. 2005, 10 (5): 365-372. 10.1016/S1359-6446(05)03369-6PubMedView ArticleGoogle Scholar
- Ward RE, German JB: Understanding milk's bioactive components: a goal for the genomics toolbox. J Nutr. 2004, 134 (4): 962S-7S.PubMedGoogle Scholar
- Guyette WA, Matusik RJ, Rosen JM: Prolactin-mediated transcriptional and post-transcriptional control of casein gene expression. Cell. 1979, 17 (4): 1013-1023. 10.1016/0092-8674(79)90340-4PubMedView ArticleGoogle Scholar
- Choi KM, Barash I, Rhoads RE: Insulin and prolactin synergistically stimulate beta-casein messenger ribonucleic acid translation by cytoplasmic polyadenylation. Mol Endocrinol. 2004, 18 (7): 1670-1686. 10.1210/me.2003-0483PubMedView ArticleGoogle Scholar
- Moshel Y, Rhoads RE, Barash I: Role of amino acids in translational mechanisms governing milk protein synthesis in murine and ruminant mammary epithelial cells. J Cell Biochem. 2006, 98 (3): 685-700. 10.1002/jcb.20825PubMedView ArticleGoogle Scholar
- Reddy TR, Xu W, Mau JK, Goodwin CD, Suhasini M, Tang H, Frimpong K, Rose DW, Wong-Staal F: Inhibition of HIV replication by dominant negative mutants of Sam68, a functional homolog of HIV-1 Rev. Nat Med. 1999, 5 (6): 635-642. 10.1038/9479PubMedView ArticleGoogle Scholar
- McLaren M, Asai K, Cochrane A: A novel function for Sam68: enhancement of HIV-1 RNA 3' end processing. Rna. 2004, 10 (7): 1119-1129. 10.1261/rna.5263904PubMed CentralPubMedView ArticleGoogle Scholar
- Paronetto MP, Achsel T, Massiello A, Chalfant CE, Sette C: The RNA-binding protein Sam68 modulates the alternative splicing of Bcl-x. J Cell Biol. 2007, 176 (7): 929-939. 10.1083/jcb.200701005PubMed CentralPubMedView ArticleGoogle Scholar
- Long E, Capuco AV, Zhao X: Cloning of bovine eukaryotic translation initiation factor 4E (eIF-4E) and its expression in the bovine mammary gland at different physiological stages. DNA Seq. 2001, 12 (5-6): 319-329.PubMedGoogle Scholar
- Hadsell D, George J, Torres D: The Declining Phase of Lactation: Peripheral or Central, Programmed or Pathological?. J Mammary Gland Biol Neoplasia. 2007, 12 (1): 59-70. 10.1007/s10911-007-9038-4PubMedView ArticleGoogle Scholar
- Stein T, Morris JS, Davies CR, Weber-Hall SJ, Duffy MA, Heath VJ, Bell AK, Ferrier RK, Sandilands GP, Gusterson BA: Involution of the mouse mammary gland is associated with an immune cascade and an acute-phase response, involving LBP, CD14 and STAT3. Breast Cancer Res. 2004, 6 (2): R75-91. 10.1186/bcr753PubMed CentralPubMedView ArticleGoogle Scholar
- Clarkson RW, Watson CJ: Microarray analysis of the involution switch. J Mammary Gland Biol Neoplasia. 2003, 8 (3): 309-319. 10.1023/B:JOMG.0000010031.53310.92PubMedView ArticleGoogle Scholar
- Stein T, Salomonis N, Gusterson BA: Mammary Gland Involution as a Multi-step Process. J Mammary Gland Biol Neoplasia. 2007, 12 (1): 25-35. 10.1007/s10911-007-9035-7PubMedView ArticleGoogle Scholar
- Hunziker W, Kraehenbuhl JP: Epithelial transcytosis of immunoglobulins. J Mammary Gland Biol Neoplasia. 1998, 3 (3): 287-302. 10.1023/A:1018715511178PubMedView ArticleGoogle Scholar
- Ollivier-Bousquet M: Transferrin and prolactin transcytosis in the lactating mammary epithelial cell. J Mammary Gland Biol Neoplasia. 1998, 3 (3): 303-313. 10.1023/A:1018767528017PubMedView ArticleGoogle Scholar
- Monks J, Neville MC: Albumin transcytosis across the epithelium of the lactating mouse mammary gland. J Physiol. 2004, 560 (Pt 1): 267-280. 10.1113/jphysiol.2004.068403PubMed CentralPubMedView ArticleGoogle Scholar
- Bos R, Zhong H, Hanrahan CF, Mommers EC, Semenza GL, Pinedo HM, Abeloff MD, Simons JW, van Diest PJ, van der Wall E: Levels of hypoxia-inducible factor-1 alpha during breast carcinogenesis. J Natl Cancer Inst. 2001, 93 (4): 309-314. 10.1093/jnci/93.4.309PubMedView ArticleGoogle Scholar
- Zhong H, De Marzo AM, Laughner E, Lim M, Hilton DA, Zagzag D, Buechler P, Isaacs WB, Semenza GL, Simons JW: Overexpression of hypoxia-inducible factor 1alpha in common human cancers and their metastases. Cancer Res. 1999, 59 (22): 5830-5835.PubMedGoogle Scholar
- Seagroves TN, Hadsell D, McManaman J, Palmer C, Liao D, McNulty W, Welm B, Wagner KU, Neville M, Johnson RS: HIF1alpha is a critical regulator of secretory differentiation and activation, but not vascular expansion, in the mouse mammary gland. Development. 2003, 130 (8): 1713-1724. 10.1242/dev.00403PubMedView ArticleGoogle Scholar
- Goeman JJ, Buhlmann P: Analyzing gene expression data in terms of gene sets: methodological issues. Bioinformatics. 2007, 23 (8): 980-987. 10.1093/bioinformatics/btm051PubMedView ArticleGoogle Scholar
- Parrish JR, Gulyas KD, Finley RL: Yeast two-hybrid contributions to interactome mapping. Curr Opin Biotechnol. 2006, 17 (4): 387-393. 10.1016/j.copbio.2006.06.006PubMedView ArticleGoogle Scholar
- NCBIs Gene Expression Omnibus . http://www.ncbi.nlm.nih.gov/geo/
- Alako BT, Veldhoven A, van Baal S, Jelier R, Verhoeven S, Rullmann T, Polman J, Jenster G: CoPub Mapper: mining MEDLINE based on search term co-publication. BMC Bioinformatics. 2005, 6: 51- 10.1186/1471-2105-6-51PubMed CentralPubMedView ArticleGoogle Scholar
- Doniger SW, Salomonis N, Dahlquist KD, Vranizan K, Lawlor SC, Conklin BR: MAPPFinder: using Gene Ontology and GenMAPP to create a global gene-expression profile from microarray data. Genome Biol. 2003, 4 (1): R7- 10.1186/gb-2003-4-1-r7PubMed CentralPubMedView ArticleGoogle Scholar
- Multtest R library. http://bioconductor.org/
- Pollard KS, Dudoit S, van der Laan MJ: Bioinformatics and Computational Biology Solutions Using R and Bioconductor. Edited by: Gentleman R, Carey V, Huber W, Irizarry R, Dudoit S. 2005, 251-272.Multiple Testing Procedures: R multtest Package and Applications to Genomics, Springer,Statistics for Biology and HealthGoogle Scholar
- Ingenuity Pathways Analysis. http://www.ingenuity.com
- Batagelj V, Mrvar A: Pajek -- Analysis and Visualization of Large Networks. Graph Drawing Software. Edited by: Jünger M, Mutzel P. 2003, 77-103. Berlin , SpringerGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.