Functional Microarray Analysis of Mammary Organogenesis Reveals a Developmental Role in Adaptive Thermogenesis
Stephen R. Master,
Jennifer L. Hartman,
Celina M. DCruz,
Susan E. Moody,
Elizabeth A. Keiper,
Seung I. Ha,
James D. Cox,
George K. Belka and
Lewis A. Chodosh
Department of Cancer Biology (S.R.M., J.L.H., C.M.D., S.E.M., E.A.K., S.I.H., J.D.C., G.K.B., L.A.C.); Department of Cell and Developmental Biology (L.A.C.); Department of Medicine (L.A.C.), Division of Endocrinology, Diabetes and Metabolism, and Abramson Family Cancer Research Institute, University of Pennsylvania School of Medicine, Philadelphia, Pennsylvania 19104-6160
Address all correspondence and requests for reprints to: Dr. Lewis Chodosh, Department of Cancer Biology, University of Pennsylvania School of Medicine, 612 Biomedical Research Building II/III, 421 Curie Boulevard, Philadelphia, Pennsylvania 19104-6160. E-mail: chodosh{at}mail.med.upenn.edu.
 |
ABSTRACT
|
---|
The use of DNA microarrays to study vertebrate organogenesis presents unique analytical challenges compared with expression profiling of homogeneous cell populations. We have used a general approach that permits the automated, unbiased identification of biologically relevant patterns of gene expression to study murine mammary gland development. Our studies confirm the utility of this approach by demonstrating the ready identification of cellular processes and pathways of known functional importance in mammary development. Additionally, this approach permitted the identification of genetic pathways with unpredicted patterns of developmental regulation, including those involved in angiogenesis, extracellular matrix synthesis, and the ß-oxidation of fatty acids. Surprisingly, our findings demonstrate that the coordinate regulation of genes involved in the ß-oxidation of fatty acids reflects the presence of an abundant, yet previously unrecognized stromal compartment within the mammary gland that is composed of brown adipose tissue. Our data demonstrate that the amount of brown adipose tissue present in the mammary gland is developmentally regulated; that PPAR
, Ucp1, and genes involved in fatty acid oxidation are spatially and temporally coregulated during development; that the mammary gland plays a functional role in adaptive thermogenesis; and that the transcriptional control of this adaptive response to cold is itself developmentally regulated.
 |
INTRODUCTION
|
---|
THE USE OF microarrays to simultaneously measure the abundance of transcripts from thousands of genes has made it possible to systematically study complex regulatory processes within the cell. To aid in interpreting the large volumes of data generated, clustering algorithms have been used to organize genes into groups of related expression patterns (1, 2). Previously, such clustered expression analyses have been used both to identify groups of genes that function in common pathways (3) and to identify clinically significant molecular differences between phenotypically similar neoplastic cell populations (4, 5). Most microarray studies to date have focused on relatively homogeneous populations of cells, whether cultured in vitro or derived from clonal primary neoplasms. This restricted focus has permitted comparatively straightforward inferences to be drawn regarding genetic regulatory networks within the cell without the analytical complications introduced by mixed cell populations. However, analysis of gene expression patterns in primary tumors has suggested that it may also be possible to identify the presence of multiple cell types (6, 7).
The use of microarray hybridization to study Drosophila metamorphosis (8) suggested the potential applicability of this technique to the analysis of vertebrate development. However, analysis of microarray expression data from an entire organism or organ over its developmental course brings an additional level of complexity because gene expression profiles reflect both the alteration of genetic programs within discrete subsets of cells in the organ as well as dynamic shifts in the relative populations of various cell types. Furthermore, in contrast to studies in which changes in gene expression patterns are monitored over intervals that approximate expected time scales for transcriptional regulation (minutes to hours), studying developmental processes such as mammalian organogenesis requires the ability to reproducibly detect changes in gene expression patterns that occur over weeks to months. As such, whether microarray approaches can be used successfully to develop a coherent picture of the biological programs used during vertebrate development remains an open question.
To address this question, we undertook a gene expression analysis of murine mammary gland development using DNA oligonucleotide microarrays. The mammary gland possesses several unique characteristics that make it an ideal experimental system for this analysis. First, mammary gland organogenesis occurs via many programs of general developmental interest including branching morphogenesis, inductive stromal-epithelial interactions, programmed cell death, extracellular matrix remodeling, and hormone action (9). Second, because much of mammary ontogeny occurs after birth, these developmental events can be conveniently manipulated and analyzed. Finally, the striking extent to which the timing of normal developmental events influences breast cancer risk in humans suggests that a thorough understanding of mammary gland development will be of clinical utility.
Our findings confirm the ability of this nondirected analytical approach to identify cellular processes and pathways of known importance in mammary development as well as shifts in the relative abundance of different cell types within the gland. In addition, this approach permitted the identification of novel aspects of mammary development, particularly with regard to our findings that the murine mammary gland participates in the adaptive response to cold and that this response occurs in a developmentally regulated manner.
 |
RESULTS
|
---|
The murine mammary gland undergoes a series of dramatic phenotypic changes during puberty, pregnancy, lactation, and postlactational involution (Fig. 1A
). At birth, the female gland consists of a rudimentary epithelial tree at one end of a mesenchymally derived fat pad composed of white adipose tissue and a fibroblastic stroma. With the onset of puberty, highly proliferative structures known as terminal end buds appear and rapidly extend the ductal tree through the fat pad via a regulated process of branching morphogenesis (10). During early pregnancy, increases in circulating hormones induce the rapid proliferation of alveolar epithelial cells within alveolar buds. Further proliferation and differentiation within alveoli generates the casein-producing cells required for milk production during lactation. After nursing pups are weaned, widespread programmed cell death and remodeling occur in the mammary gland, resulting in a gland superficially resembling, but distinct from, that of adult nulliparous females.

View larger version (46K):
[in this window]
[in a new window]
|
Figure 1. Mammary Gland Development
A, Whole mounts of the no. 4 murine mammary gland demonstrate the dramatic growth of the epithelial ductal tree during puberty and of the alveolar epithelial compartment during pregnancy. After lactation and postlactational involution, the gland returns to a state superficially resembling the adult nulliparous gland. Glands shown were harvested at 3 wk (onset of puberty), 5 wk (mid-puberty), 10 wk (completion of puberty), 15 wk (adult nulliparous), d 14 pregnancy, d 9 lactation, d 7 involution, and d 28 involution. The nipple (N), epithelial tree (Ep), fat pad (FP), and lymph node (LN) are indicated. B, Multidimensional scaling. The set of pairwise distances between developmental time points was reduced to three dimensions using multidimensional scaling. Lines connecting consecutive time points in pregnancy (red), lactation (white), and involution (black) are shown along with a convex shape bounding all time points in nulliparous development (green). V2: 2-wk G0P0; V5: 5-wk G0P0; V10: 10-wk G0P0; V15: 15-wk G0P0; Pd6: d 6 pregnancy (G1P0); Pd12: d 12 pregnancy (G1P0); Pd18: d 18 pregnancy (G1P0); Lh18: lactation h 18 (G1P1); Ld9: lactation d 9 (G1P1); Id2: involution d 2 (G1P1); Id7: involution d 7 (G1P1); Id28: involution d 28 (G1P1).
|
|
To identify genetic pathways that are regulated during this developmental cycle, total RNA was prepared from pooled mammary glands harvested from cohorts of mice at each of 13 time points during mammary development. One male time point (9.5 wk) is represented, and the remaining 12 female time points span key events in postnatal mammary development including those before puberty (2-wk nulliparous), during puberty (5-wk nulliparous), after puberty (10- and 15-wk nulliparous), and during pregnancy (d 6, 12, and 18), lactation (18 h and 9 d post partum), and postlactational involution (d 2, 7, and 28). Samples were analyzed by hybridization to oligonucleotide microarrays representing approximately 5,300 murine genes and expressed sequence tags, and raw expression levels were scaled to correct for variation among chips.
Global Analysis
To first determine whether global patterns of gene expression during mammary gland development are consistent with previous biological observations, pairwise "distances" between developmental points were calculated using the Pearson correlation coefficient and were visualized in three-dimensional space using multidimensional scaling (Fig. 1B
). Similar methods have previously been used to represent global relationships between individual tumors at the level of gene expression (11). Our analysis of mammary gland gene expression revealed that the four points of nulliparous development occupy a discrete region of gene expression space that is separable from all other surveyed points (Fig. 1B
). In addition, global expression profiles display a monotonic progression through gene expression space as the mammary gland progresses from the nulliparous state through pregnancy toward lactation, with the greatest differences in gene expression existing between the nulliparous and lactating mammary gland. Conversely, as the mammary gland undergoes postlactational involution the pattern of gene expression returns toward a state most similar to that of the nulliparous gland. However, the trajectory in gene expression space that the mammary gland follows from lactation through involution is distinct from that which describes its path from the nulliparous state to that of lactation. Presumably, this reflects the fact that the cellular processes by which the mammary gland progresses from nulliparity to lactation (e.g. proliferation and differentiation) are not simply the inverse of those by which the gland returns to the involuted state (e.g. apoptosis and remodeling). Notably, the global gene expression pattern found after 28 d of involution is similar, but not identical, to the age-matched adult nulliparous gland, consistent with the existence of morphological differences between the nulliparous and parous involuted gland, as well as functional differences with respect to cancer susceptibility. In aggregate, these data are consistent with previous morphological descriptions of postnatal mammary development, suggesting that these gene expression data represent a reasonable measure of underlying developmental processes.
While pairwise correlation coefficients provide a means to estimate similarities between global gene expression profiles, other analytical approaches are required to identify coordinate regulatory patterns among subsets of genes. Accordingly, profiles of genes whose expression exceeded a minimal threshold of regulation during development were normalized and clustered using a two-dimensional self-organizing map (SOM). As expected, these clusters include a variety of physiologically suggestive patterns of developmental expression indicative of regulation during puberty, pregnancy, lactation, and postlactational involution (Fig. 2
). Complete lists of the genes contained within each cluster and their expression values are published as supporting information at www.abramsoninstitute.org/chodoshdata.html.

View larger version (53K):
[in this window]
[in a new window]
|
Figure 2. Clustered Expression Profiles over Mammary Development
An 8 x 5 SOM depicting clusters of normalized mammary gland developmental gene expression patterns. Clusters are numbered (upper left), and the number of genes present in each cluster is indicated (upper right). Developmental points shown are: male (blue; 9.5 wk), female nulliparous (green; 2, 5, 10, and 15 wk), pregnancy (red; d 6, 12, and 18), lactation (white; 18 h, 9 d), postlactational involution (black; d 2, 7, and 28). Gray shading indicates expression ± 1 SD of genes within each cluster. Genes contained within each cluster, along with their expression values, can be found at www.abramsoninstitute.org/chodoshdata.html.
|
|
Associations between Clustered Expression Data and Functional Gene Categories
The detailed interpretation of clustered gene expression profiles is heavily dependent upon the ability to identify functionally relevant groups of genes. However, as the number of sequences that can be monitored using microarray technology is rapidly increasing, this task has become correspondingly difficult. As such, we attempted to utilize existing sources of annotation to highlight statistically significant associations between functional gene categories and temporal patterns of gene expression during mammary development. To generate an annotated data set suitable for this approach in the mouse, the gene sequences represented on the murine oligonucleotide microarrays used in this study were mapped to related SWISS-PROT/TrEMBL entries and keyword annotations were extracted in an automated manner. We then implemented software to search for functional pathways regulated during mammary gland development by identifying SWISS-PROT categories having a statistically significant association with specific clustered developmental expression profiles.
Thirty independent clustering runs were tested for significant associations between expression clusters and gene categories (see Materials and Methods). Overlapping sets of genes within individual categories that exceeded the significance threshold in any individual clustering run were pooled. To confirm the statistical significance of associations identified by this approach, we used data sets in which gene identifications had been randomly permuted with expression data to generate 100 data sets for statistical characterization. This analysis demonstrated that each individual category appearing in the original compiled list has a probability of chance appearance less than 0.03. This approach resulted in the identification of 49 gene categories that are associated in a statistically significant manner with one or more temporal expression profiles during mammary gland development (Table 1
). Conversely, 20 clustered developmental gene expression profiles were statistically associated with one or more functional gene categories.
Automated Identification of Functional Pathways in the Mammary Gland
As a first demonstration of the applicability of an automated, nondirected analytical approach to the study of mammalian organogenesis, we determined whether such an approach could detect discrete cellular processes central to mammary development. For example, the large population of milk-producing alveolar cells required for lactationa cardinal feature of mammary functionis generated by a massive wave of epithelial proliferation during the first half of pregnancy (Fig. 1
; Ref. 9). This wave of proliferation is reflected in the statistically significant association (P < 0.03) between genes in the "Cell cycle" and "DNA replication" categories and clusters manifesting peak expression at d 6 of pregnancy (Fig. 2
, clusters 12 and 17; Table 1
). In fact, the temporal profile of epithelial proliferation rates during pregnancy as revealed by bromodeoxyuridine labeling is essentially identical to that observed for cluster 17 (12). Genes associated with the "Cell cycle" and "DNA replication" categories encode cyclin A2, cdk4, components of the replication machinery such as PCNA and the catalytic subunit of DNA Pol
, and four components of the MCM complex, which is required to license one round of DNA replication (13). Further examination of cluster 17 reveals the presence of additional genes associated with proliferation and cell cycle progression (Fig. 3A
). These include genes whose products transduce proliferative signals (e.g. c-myc) as well as the enzyme that catalyzes the rate-limiting step in guanine nucleotide synthesis, IMP dehydrogenase 2. Cluster 17 also contains several nucleolar genes, including the proliferation-associated argyrophilic nucleolar organizer region protein, nucleophosmin (14), and nucleolin, which is a putative transcriptional target of c-myc (15). Interestingly, while this cluster contains many genes associated with proliferation, multiple members of the cyclin family are preferentially expressed during other stages of mammary development, including puberty (cyclin D2), late pregnancy/lactation (cyclin E), and early involution (cyclin G), suggesting the specific utilization of different cell cycle components during distinct stages of mammary gland development (Fig. 3C
).

View larger version (52K):
[in this window]
[in a new window]
|
Figure 3. Coordinate Gene Expression during Mammary Development
A, Partial set of individual gene expression profiles from cluster 17 generated using TreeView (1 ). Low expression is represented by green, and high expression is represented by red. Note the common peak of expression during early pregnancy. Genes that were previously identified as members of a mammary epithelial proliferation cluster (6 ) are indicated ( ). B, Expression of cathepsin family members is differentially regulated during distinct stages of involution. In addition, cathepsins S and Z have an expression profile highly similar to cathepsins B and K (data not shown). C, Expression profiles of cyclin family members. D, Adipocyte-specific markers show a coordinate decrease from puberty (5 wk G0P0) to pregnancy (d 6, 12, 18). E, Detection of Ob (Lep) mRNA expression by in situ hybridization reveals significant expression in adipocytes (Ad) but not epithelial cells (Ep). Comparison of expression in pubertal (6-wk nulliparous) vs. pregnant (d 14 pregnancy) glands demonstrates that the decrease seen in (A) cannot be accounted for by a drop in per-cell Ob expression. F, Cytokeratin 14increases significantly from juvenile (2 wk G0P0) to adult nulliparous (10 wk G0P0), and from adult nulliparous (15 wk G0P0) through early pregnancy, consistent with the rapid increase in epithelial cell numbers during these stages of development.
|
|
Our results further revealed that clustered expression profiles that peak during lactation were statistically associated (P < 0.03) with several gene categories, including "Milk," "Protein biosynthesis," and "Protein transport" (Tables 1
and 2
; Fig. 2
, e.g. clusters 32, 36, and 38). "Milk" proteins that are significantly associated in this manner during clustering include
, ß,
,
, and
-casein along with whey acidic protein,
-lactalbumin, and milk fat globule-EGF factor 8. "Protein biosynthesis" genes include those encoding eIF-4B, eIF-5, EF-Tu, EF-1
, and EF-2 and a number of aminoacyl-tRNA synthases, consistent with the requirement for large-scale production of milk proteins. "Protein transport" genes include sequences encoding each of the components (
, ß,
) of the heterotrimeric SEC61 complex responsible for protein transport into the endoplasmic reticulum (16) as well as homologs of other genes involved in protein secretion (SEC13, SEC23). Further inspection of clusters with peak expression during lactation reveals genes up-regulated during de novo fatty acid synthesis (ATP citrate-lyase, cluster 36), mammary epithelial differentiation (H-FABP, cluster 36), and the stimulation of adipocyte lipid degradation (zinc alpha-2 glycoprotein, cluster 37) (17, 18, 19).
After lactation, withdrawal of the suckling stimulus leads to mammary gland involution (Fig. 1
), which proceeds via two main stages: an initial stage characterized by the induction of widespread programmed cell death within 48 h, and a subsequent overlapping stage in which extracellular matrix degradation and remodeling occur (20). Consistent with this model, examination of the SOM in Fig. 2
reveals clusters of genes up-regulated exclusively at either d 2 (stage one; cluster 29) or d 7 (stage two; clusters 13 and 14) of involution, as well as clusters of genes expressed at both stages (clusters 19 and 24). Examination of stage one clusters reveals a number of genes, including SGP-2, Stat3, Igfbp5, and C/EBP
, that have been previously associated with apoptosis during early involution (21, 22). Consistent with the widespread degradation and remodeling of the extracellular matrix that occurs during the second phase of mammary gland involution, transcripts encoding "Thiol protease" and "Zymogen" proteins with previously described roles in involution, such as stromelysin 1 (Mmp3), stromelysin 3 (Mmp11), and MT-MMP (Mmp14), are statistically associated (P < 0.03) with this stage of development (cluster 14; see Table 1
) (20). Intriguingly, specific cathepsin isoforms are also observed to be expressed in waves during successive stages of mammary gland involution (Fig. 3B
) consistent with previous suggestions for a general role for cathepsins in apoptosis and extracellular matrix remodeling (23, 24).
Identification of Shifts in Cellular Compartments during Organogenesis
In addition to reflecting changes in transcriptional regulatory patterns within individual cells, changes in gene expression observed during organogenesis may also reflect shifts in the relative abundance of specific cell types. For example, a number of adipocyte markers, including C/EBP
, aP2, Fsp27, ß3-adrenergic receptor, Ob, and PPAR
, were noted to share a similar pattern of expression during puberty and pregnancy (Fig. 3D
) (25, 26). This finding could reflect changes in the level of expression of these genes within each adipocyte or changes in the relative abundance of adipocytes within the gland. To distinguish between these possibilities, the expression of Ob in mammary adipocytes was determined using in situ hybridization. This analysis demonstrated that, in contrast with the marked decrease in Ob expression detected by microarray analysis during pregnancy, Ob expression does not change substantially on a per-cell basis during puberty and pregnancy (Fig. 3E
). This suggests that the progressive reduction in adipocyte-specific mRNAs as a fraction of total mRNA during pregnancy is more likely to be due to dilution by the large increase in epithelial cell abundance and/or levels of epithelial RNA synthesis that occurs during this developmental transition. Indeed, the marked increase in epithelial cells during puberty and pregnancy can be directly visualized using cytokeratin 14 (a marker of myoepithelial cells) as well as other cytokeratins (Fig. 3F
and data not shown). Furthermore, the statistical association (P < 0.01) between members of the "Keratin" gene category and cluster 22 suggests a proportional increase of other epithelial cell subtypes, including subsets of luminal epithelial cells (cytokeratin 19) (27) and putative stem cells (cytokeratin 6) (28). These data suggest that the use of selected cell type-specific genes whose expression across development is relatively constant on a per-cell basis will facilitate the rapid identification of genes expressed in the same cell type and the establishment of canonical cell type-specific expression profiles against which regulatory perturbations may be recognized and measured.
Beyond identifying shifts among abundant cell types, this approach can also readily detect changes in cell types that constitute only a small fraction of the mammary gland. For example, the presence of macrophages in the mammary gland during involution is indicated by the dramatic up-regulation of macrophage-specific genes such as Mac-2, Nramp1, F4/80, macrophage metalloelastase, and MPS-1 at d 7 of postlactational involution (clusters 13 and 14) (29, 30, 31, 32, 33, 34), as well as by a number of genes whose products have been postulated to mediate the detection and engulfment of apoptotic cells and debris, including C1q (chains A, B, and C), CR3ß, Gas-6, Axl, Cd68, and Abc1, consistent with a role for macrophages in clearing apoptotic epithelial cells after lactation (35). Together, our data indicate that microarray analysis can detect altered cell numbers as well as suggest likely functional roles within the context of a whole organ.
Related to this, endothelial-specific genes, including Pecam, VE-cadherin, and endoglin, were found to exhibit a sharp peak of expression in glands harvested from 2-wk-old animals (Fig. 2
, cluster 10) (36), as were a number of genes involved in the coordination of angiogenesis, such as Flk-1, Tie1, Tie2/Tek, and Vegf (Fig. 2
, clusters 10 and 15) (37). These data suggest active formation of the microvasculature during early postnatal development. Further examination reveals that genes encoding basement membrane proteins, including the
-1 and
-2 chains of type IV collagen, heparan sulfate proteoglycan 2, and the laminin B2 chain, also demonstrate peak expression in the mammary glands of 2-wk-old animals (Fig. 2
, clusters 10 and 11). This observation is intriguing in light of in vitro work demonstrating that endothelial-specific markers can be detected in cultured mammary stromal cells capable of differentiating into adipocytes or forming capillary-like structures in a Vegf-dependent manner when plated on basement membrane components (38). Our results suggest that basement membrane protein expression may provide a similar signal for endothelial differentiation of mammary stromal cells during early postnatal development. Alternately, the coordinate expression of endothelial markers and basement membrane proteins may reflect the simultaneous formation of an extensive capillary network characteristic of white adipose tissue along with deposition of basement membrane associated with each adipocyte. While further work is required to test these hypotheses, the ability to conceptually link disparate pathways on the basis of gene expression patterns during organogenesis demonstrates the utility of this approach for suggesting novel hypotheses regarding vertebrate ontogeny.
The Mammary Gland Contains Brown Adipose Tissue
Having successfully used a statistical approach to identify known developmental processes in the mammary gland, we next examined our data set for groups of genes with previously unsuspected patterns of coordinate developmental regulation. For example, we found that genes encoding "Fatty acid metabolism" enzymes were statistically associated (P < 0.02) with expression clusters manifesting peak expression in the mammary glands of 2-wk-old animals (clusters 5 and 10). This pattern was confirmed by Northern hybridization analysis (Fig. 4A
, Dci and Acadvl). The eight enzymes in the "Fatty acid metabolism" category that map to clusters 5 and 10 catalyze sequential steps in the ß-oxidation of both saturated and unsaturated fatty acids, including the transport of fatty acyl-coenzyme A (CoA) across the inner mitochondrial membrane via the carnitine shuttle, the breakdown of very-long-/long-/medium-/short-chain fatty acyl-CoA, and the isomerization of cis-
3-enoyl CoA (Fig. 4B
) (39, 40). Eight of nine clustered genes in this category exhibit peak expression in the mammary gland at 2 wk of age, emphasizing the striking coordinate regulation of the majority of enzymes in this metabolic pathway.

View larger version (89K):
[in this window]
[in a new window]
|
Figure 4. Coordinate Regulation of Genes Involved in Fatty Acid Metabolism
A, Northern blot blots showing expression levels of Ucp1, PPAR , dodecenoyl-CoA -isomerase (Dci), very-long-chain acyl-CoA dehydrogenase (Acadvl), and Gata-3 across mammary gland development are shown. Gata-3 is expressed in mammary epithelial cells rather than the mammary fat pad (data not shown). B, Coregulation of fatty acid metabolism genes. Genes that encode enzymes necessary for the ß-oxidation of saturated and unsaturated fatty acids and that are up-regulated in early nulliparous development (2-wk) are indicated: 1) Carnitine acetyltransferase (Crat), 2) Carnitine palmitoyltransferase II (Cpt2), 3) Short chain acyl-CoA dehydrogenase (Acads), 4) Medium chain acyl-CoA dehydrogenase (Acadm), 5) Long chain acyl-CoA dehydrogenase (Acadl), 6) Very-long-chain acyl-CoA dehydrogenase (Acadvl), 7) 3-hydroxyacyl CoA dehydrogenase (Hadh), 8) Dodecenoyl-CoA -isomerase (Dci). C, Multilocular adipose tissue is seen in hematoxylin and eosin-stained sections of the nulliparous no. 4 murine mammary gland at 2 wk of age. A region of the mammary fat pad without significant multilocular adipose tissue is shown (white adipose tissue; left) along with a region containing large amounts of multilocular adipose tissue (brown adipose tissue; right). The lymph node contained within the no. 4 mammary gland is also seen. Low power. D, High-power view of mammary adipose tissue shown in C. E, In situ hybridization of serial sections demonstrating the expression of Ucp1 (left panels) and aP2 (right panels) in the no. 4 mammary gland (MG), interscapular fat pad (IFP), and liver (L) of the 2-wk-old female FVB mouse. Regions of the mammary gland with high (MG, upper) and low (MG, lower) Ucp1 expression are shown. Brown (BAT) and white (WAT) adipose tissue are indicated in sections from interscapular fat pad.
|
|
Further examination of developmental clusters 5 and 10 yielded the surprising finding that uncoupling protein 1 (Ucp1), which encodes a protein that dissipates the proton gradient across the mitochondrial inner membrane, is also maximally expressed in the murine mammary gland during early postnatal development (Fig. 4A
and data not shown). This finding was unexpected given that Ucp1 is specifically expressed in brown adipose tissue (41), whereas the adipocyte compartment of the mammary gland has long been described as being composed of white adipose tissue (42, 43, 44, 45, 46). As such, the expression of Ucp1 suggests the possibility that, contrary to previous observations, the mammary gland contains brown adipose tissue. Such a compartment has not previously been recognized within the mammary gland.
Because brown adipose tissue plays a central role in adaptive thermogenesis, we further considered the possibility that the mammary gland may play a role in this process. Moreover, because adaptive thermogenesis requires the ß-oxidation of fatty acids within brown adipose tissue (47), a role for the mammary gland in nonshivering thermogenesis might account for the developmental coregulation of "Fatty acid metabolism" genes and Ucp1. To begin to test this hypothesis, in situ hybridization for Ucp1 expression was performed on serial sections from 2-wk-old mammary gland. This analysis demonstrated that Ucp1 is expressed at high levels in defined, lobular regions of eosinophilic, multilocular adipose tissue located within the mammary gland that are histologically characteristic of brown fat (Fig. 4
, CE). Notably, these Ucp1-expressing regions encompass only a subset of adipocytes within the mammary gland because the adipose-specific marker, aP2, is expressed in additional regions of the gland that do not express Ucp1 (Fig. 4E
). Similarly, aP2 expression can be seen in both interscapular brown fat as well as the surrounding white adipose tissue, whereas Ucp1 expression in the interscapular fat pad is restricted to the regions of brown fat.
Further examination of serial sections from this analysis revealed that genes (Acadvl, Dci) encoding enzymes that catalyze the ß-oxidation of fatty acids are expressed in a spatially restricted pattern that is identical to that of Ucp1 (Fig. 5
). Moreover, expression of the nuclear hormone receptor, PPAR
, which has been shown to regulate the expression of both Ucp1 (48) and fatty acid metabolic genes (49 ; for review, see Ref. 50), was found to exhibit developmental regulation (cluster 15) and spatial expression strikingly similar to that of Ucp1, Acadvl, and Dci (Figs. 4
and 5
). As such, the spatial and temporal coregulation of Ucp1, PPAR
, and genes whose products are involved in the ß-oxidation of fatty acids further supports our hypothesis that the statistical association between "Fatty Acid Metabolism" genes and developmental expression profiles with peak expression in the neonatal mammary gland may be related to the PPAR
-mediated coordinate activation of gene regulatory networks involved in adaptive thermogenesis.

View larger version (90K):
[in this window]
[in a new window]
|
Figure 5. Spatial Colocalization of Fatty Acid Metabolism Genes with Ucp1
In situ hybridization of serial sections was used to demonstrate the colocalization of Ucp1, PPAR , Very-long-chain acyl-CoA dehydrogenase (Acadvl), and Dodecenoyl-CoA -isomerase (Dci). Ob, which is expressed in both brown and white adipose tissue, is shown.
|
|
A Potential Role for the Mammary Gland in Adaptive Thermogenesis
If the above findings in fact represent the identification of a developmentally regulated pathway in the mammary gland similar to that involved in adaptive thermogenesis, we reasoned that the mammary gland would be capable of responding to cold exposure. It has previously been shown that the PPAR
coactivator, PGC-1, is induced in response to cold and activates adaptive thermogenesis by up-regulating Ucp1 (51). Furthermore, PPAR
and PGC-1 can interact to up-regulate the mitochondrial transcription of fatty acid metabolic genes (52). Therefore, to determine if a cold-inducible signaling pathway is present within the mammary gland, 2- or 10-wk-old female FVB mice were exposed to cold (4 C) for 3 h before tissue harvest. Northern hybridization analysis demonstrated a 6-fold up-regulation of PGC-1 in the mammary glands of 2-wk-old mice exposed to cold, compared with animals housed at room temperature (Fig. 6
, A and B). In contrast, while an increase in PGC-1 mRNA levels was observed within the glands of 10-wk-old mice after cold exposure, the absolute level of expression observed in the cold-induced state was comparable to that seen in 2-wk-old animals housed at room temperature. These data indicated that the classically described coordinate transcriptional response to cold characteristic of brown adipose tissue is present in the neonatal mammary gland.

View larger version (49K):
[in this window]
[in a new window]
|
Figure 6. PGC-1 Induction upon Cold Exposure.
Two- or 10-wk-old female FVB mice were housed at room temperature or 4 C for 3 h before tissue harvest. A, Northern hybridization of PGC-1 to RNA from 10-wk liver, 2-wk mammary gland, or 10-wk mammary gland from mice housed at room temperature (W) or 4 C (C). B, PhosphorImager quantitation. Expression levels were normalized to ß-actin, and the mean expression of room temperature 2-wk mammary gland was set to 1. C, Cold-induced changes in mammary gene expression. Two- or 10-wk old FVB female mice were kept at room temperature (W) or 4 C (C) for 3 h before harvest. In situ hybridization was performed with antisense probes to PGC-1, Ucp1, and PPAR . D, Progressive induction of Ucp1 expression by cold exposure. Two- or 10-wk old FVB female mice were housed at 4 C for the indicated time period (2-wk: 0, 3 h; 10-wk: 0, 1, 3, 6, 10 h) before harvest. Low-power views of in situ hybridization slides reveal the spatial pattern of Ucp1 expression across the entire no. 4 mammary gland.
|
|
The increase in PGC-1 expression observed in the mammary glands of cold-induced 10-wk-old mice suggested the presence of brown adipose tissue or other cells capable of regulating transcription in response to cold-induced signals. Although the absolute levels of cold-induced PGC-1 expression within the adult mammary gland were found to be lower than those seen in the mammary glands of 2-wk-old mice, we considered the possibility that this reflected a developmentally regulated decrease in the amount of brown adipose tissue within the mammary gland similar to that described for other white adipose depots (53) rather than an intrinsic inability of brown adipose tissue within the gland to respond to cold-induced signals. To address this hypothesis, we examined PGC-1, PPAR
, and Ucp1 expression by in situ hybridization in mammary glands of 2-wk-old and 10-wk-old mice after cold exposure. Analysis of Ucp1-expressing cells in the mammary glands of 2-wk-old mice confirmed that brown adipose tissue comprises a significant fraction of the neonatal gland (Fig. 6
, C and D). This analysis further demonstrated that PGC-1, PPAR
, and Ucp1 expression colocalize and that are each up-regulated in response to cold exposure in restricted regions of the adult (10-wk-old) gland to levels that are comparable to those seen in 2-wk-old mice (Fig. 6C
). These findings suggest that while basal levels of Ucp1 and PGC-1 are barely detectable in the mammary glands of 10-wk-old mice, the adult gland remains capable of responding to cold by activating transcriptional programs characteristic of adaptive thermogenesis. Interestingly, while PGC-1 expression was induced by exposure to cold in the 2-wk-old gland (Fig. 6
, AC), Ucp1 and PPAR
were not observed to increase beyond their already elevated levels of expression (Fig. 6C
). Thus, while cells within the neonatal mammary gland respond appropriately to cold-induced signals by inducing PGC-1, this signal does not appear to influence expression of the gene encoding the key thermogenic mediator, Ucp1.
Finally, having demonstrated that the adult gland remains capable of inducing pathways required for adaptive thermogenesis, we extended our study to 10-wk-old animals housed at 4 C for increasing periods of time before harvest. This analysis demonstrated that, whereas Ucp1 is expressed at high levels throughout the mammary glands of 2-wk-old mice, only a small fraction of adipocytes within the mammary glands of 10-wk-old animals are capable of expressing Ucp1 in response to acute cold exposure. Furthermore, although increasing periods of cold exposure progressively induce Ucp1 expression in the 10-wk gland, the absolute expression of Ucp1 on a per-cell basis is comparable to that seen within the 2-wk-old gland either at room temperature or following cold exposure (Fig. 6D
). These data demonstrate that both the amount of brown adipose tissue and the transcriptional regulation of the mammary glands adaptive response to cold are developmentally regulated within the murine mammary gland.
 |
DISCUSSION
|
---|
We have systematically analyzed global gene expression patterns in the murine mammary gland during puberty, pregnancy, lactation, and postlactational involution using oligonucleotide microarrays. After clustering developmental expression profiles using SOMs to identify predominant patterns of genes expression, we tested for statistical associations between clustered gene expression profiles and biological annotations derived from SWISS-PROT. This approach permitted the identification of cellular processes and pathways of known importance in mammary development as well as shifts in the relative abundance of different cell types within the gland. Moreover, by identifying genetic pathways with unpredicted patterns of developmental regulation, this approach was shown to facilitate the generation of novel hypotheses. Specifically, we have shown that: the mammary gland contains an abundant, previously unrecognized stromal compartment composed of brown adipose tissue; that expression of PPAR
, Ucp1, and genes involved in the ß-oxidation of fatty acids colocalize within this compartment and are coordinately regulated during development; that the mammary gland plays a functional role in adaptive thermogenesis; and that transcriptional control of the adaptive response to cold within mammary brown adipose tissue is developmentally regulated. These findings suggest a previously unrecognized, developmentally regulated role for the murine mammary gland in adaptive thermogenesis. Notably, despite the fact that current wisdom regarding mammary development considers the neonatal period to be relatively quiescent, studies clearly demonstrate that novel and interesting events occur in the mammary gland during this stage of development.
A central feature of the analytical approach described in this paper is the identification of statistically significant relationships between prospectively annotated functional gene categories and clustered gene expression profiles. This approach significantly decreases the analytical burden imposed by examination of large numbers of genes and avoids several potential sources of selection bias in pathway identification. In the present study, we have demonstrated the validity of this approach for studying vertebrate development by confirming its ability to identify known biological processes during mammary gland development. The automated identification of statistical relationships between functional gene categories (e.g. cell cycle, cell division, milk protein, protein biosynthesis, protein transport, thiol protease, and zymogen genes) and specific stages of mammary gland development proves that this approach is able to reliably identify biologically relevant pathways within a complex mixture of cell types during organogenesis in the absence of prior information about the developmental process analyzed. These data confirm the utility of this technique for studying mammary development and suggest its broader applicability to other developmental systems.
We have further demonstrated the ability of microarray approaches to detect shifts in the relative abundance of adipocytes and epithelial cells within the mammary gland, as well as cell types that constitute only a small percentage of the gland such as macrophages, endothelial cells, and lymphocytes. While our data show that it is clearly possible to detect coordinate pathway regulation in complex mixtures of cell types during organogenesis, interpretation of expression profiles from intact organs will undoubtedly be aided by the identification of genes that mark the contribution of specific cell types. As such, establishing baseline expression profiles for diverse cell types should permit a more detailed description of changes in tissue composition as well as facilitate the identification of developmentally regulated genes by identifying gene expression patterns within a compartment that depart from that compartments canonical baseline expression profile.
Adaptive Thermogenesis in the Mammary Gland
Perhaps the greatest value in pursuing a nondirected analytical approach lies in the ability to detect unexpected patterns of developmental regulation. In this regard, it is noteworthy that the statistical approach employed in this study was able to identify the presence of an abundant cellular compartment within the mammary gland that direct examination over several decades failed to reveal. Our identification of the coordinate up-regulation of genes involved in fatty acid oxidation, coupled with the recognition that this pattern reflects the presence of brown adipose tissue within the mammary gland, demonstrates that this methodology can identify novel and unpredicted features of vertebrate development.
One potential explanation for the observed elevation in fatty acid metabolism gene expression in the mammary glands of 2-wk-old mice is that it reflects the increased use of ß-oxidative pathways to provide the energy requirements of the neonatal gland. This explanation seems unlikely, however, because cells within the neonatal gland are relatively quiescent when compared with the rapid epithelial proliferation that occurs during puberty and pregnancy, or to the large-scale milk production that occurs during lactation. A second hypothesis to explain this pattern of gene expression is that it represents a metabolic response to changes in dietary composition. Inasmuch as 2-wk-old mouse pups are nursing on diet of high-fat rodent milk, it is likely that the mammary gland adapts to utilize the predominant source of available energy. In fact, regulation of ß-oxidative pathways during the neonatal period has been described in other tissues (54). Thus, this developmental expression pattern may represent an adaptation of the organism to its environment rather than a local phenomenon in the mammary gland.
A third hypothesis is suggested by our observation that Ucp1 is expressed in the mammary gland and is developmentally regulated in a manner similar to that observed for genes involved in fatty acid oxidation. Although the mammary fat pad is composed of white adipose tissue, the expression of a gene whose product dissipates the proton gradient across the mitochondrial inner membrane suggested the presence of cells that perform thermogenic functions characteristic of brown fat. Because thermogenesis also requires the ß-oxidation of fatty acids (47), it is possible that this pathway is up-regulated in the mammary gland in response to a need to maintain body temperature during the neonatal period. This interpretation is supported by the spatial coregulation of Ucp1 and genes encoding enzymes involved in fatty acid metabolism. The increased rate of heat loss in neonatal mice due to their comparatively high surface:volume ratio and relative lack of fur suggests the potential advantage of utilizing the mammary fat pad for thermogenesis during early postnatal development. As would be predicted by this model, the mammary glands of male neonatal mice show a similar pattern of Ucp1 expression (data not shown). This model is also consistent with the described developmental regulation of the amount of brown fat in the interscapular fat pad in rodents, with the induction of brown-fat-like regions within white adipose depots, and with the emerging recognition that brown fat and white fat may be interconvertable under some experimental conditions (53, 55, 56). This phenomenon has not been described in the mammary gland, however, validating the ability of this approach to identify unexpected aspects of mammary development.
In addition to demonstrating that the murine mammary gland contains functional depots of brown adipose tissue that respond appropriately to cold, our studies further demonstrate that this response is developmentally regulated. Specifically, the mammary glands of 10-wk-old mice do not express significant levels of Ucp1, PGC-1, or PPAR
when housed at room temperature, but coordinately induce each of these genes in response to cold exposure. In contrast, Ucp1 and PPAR
are each maximally up-regulated in neonatal mice housed at room temperature. These findings indicate that the stimuli that trigger this coordinate transcriptional response differ between neonatal and adult animals. Consistent with this, cold exposure caused an up-regulation of the transcriptional coactivator PGC-1 in the mammary glands of 2-wk-old mice, yet this increase was not accompanied by an increase in the expression of its transcriptional target, Ucp1. This suggests that the up-regulation of PGC-1 is not rate-limiting for Ucp1 expression at this developmental stage. These findings both highlight the developmental regulation of this system and raise the question of what signal is responsible for the up-regulation of Ucp1 expression in neonatal mice housed at room temperature. Potentially, the difficulty that neonatal mice have in maintaining body temperature may itself provide this signal. That is, even when housed at room temperature neonatal mice may exist in what is essentially a cold-induced state. Previous work in rat models supports this hypothesis (57). The elevated expression of Vegf observed in 2-wk-old mice (data not shown) is also consistent with this model because brown adipose tissue has been shown to produce Vegf in response to adrenergic stimulation (58). To our knowledge, this aspect of the developmental regulation of the transcriptional response to cold exposure has not been previously described.
Further consideration of the potential mechanisms underlying this effect suggests alternative developmental signals for mammary brown fat activity. Given that PGC-1 has been shown to act as a transcriptional coactivator for PPAR
(52, 48), and that high levels of PPAR
mRNA are found in the 2-wk-old mammary gland, it is possible that basal levels of PGC-1 and PPAR
expression are sufficient for high levels of expression of Ucp1 in the mammary glands of neonatal mice housed at room temperature. Notably, a high-fat diet has been shown to induce PPAR
expression in the liver (59), raising the possibility that maternal milk provides a developmental signal that induces a protective, thermogenic response in suckling pups. It is also possible, however, that other factors may contribute to the up-regulation of Ucp1 at this developmental stage. Further experiments will be required to dissect the developmental and environmental stimuli responsible for these observations.
In aggregate, our findings demonstrate that statistical methods for detecting associations between functional gene categories and patterns of developmental gene expression allow the rapid, automated identification of key pathways and processes used during mammary development. Additionally, our observation that genes involved in the ß-oxidation of fatty acids are developmentally coregulated with Ucp1 and PPAR
underscores the ability of this technique to suggest novel hypotheses regarding developmental processes. We anticipate that further characterization of clustered gene expression profiles during normal mammary development will provide a framework for identifying and classifying a broad array of pathways that are altered in transgenic and knockout animals exhibiting mammary developmental phenotypes. When combined with studies to determine the spatial localization of gene expression, the parallel analysis of developmental gene expression patterns should facilitate the elucidation of regulatory networks guiding vertebrate development.
 |
MATERIALS AND METHODS
|
---|
RNA Isolation
The nos. 35 mammary glands were harvested from FVB mice at the indicated time points, and the lymph node within the no. 4 gland was visualized and removed. Glands were snap frozen, and total RNA was isolated as described (60). RNA was pooled from approximately 20 mice for each developmental time point to provide an average result for each developmental time point that would be relatively independent of a particular stage of the estrus cycle. For pregnancy, lactation, and early-mid (d 2 and d 7) involution time points, mice were mated at 14 wk. Pups were forcibly weaned after 9 d of lactation in preparation for harvest of early-mid (d 2 and d 7) involution glands. To match the ages of adult nulliparous and parous samples, mice used to generate late (d 28) involution samples were mated at 4 wk of age, and pups were allowed to wean themselves naturally. Thus, late involution samples are age-matched (15 wk) with nulliparous adult animals. For adaptive thermogenesis experiments, three to six mice were housed together at room temperature or 4 C for the indicated times immediately before harvest.
Oligonucleotide Microarray Hybridization
Approximately 30 µg total RNA was used at each developmental time point. Biotinylated cRNA was generated essentially as described (61). Hybridization to a set of four Affymetrix Mu6500 GeneChips was performed overnight per manufacturers instructions. After staining and washing, chips were scanned using a Hewlett-Packard Co. GeneArray scanner. Grid alignment and raw data generation was performed using Affymetrix GeneChip 3.1 Software. Raw expression values, representing the average difference in hybridization intensity between oligonucleotides that perfectly match the transcript sequence and oligonucleotides containing single base pair mismatches, were measured. To scale between chips, these expression values were rank ordered, and the median 96% were averaged. Chips were scaled relative to each other to equalize this average value. A noise value (Q) based on the variance of low-intensity probe cells was used to calculate a minimum expression threshold (2.1*Q) for each chip, and the highest value across development for each of the four Mu6500 subarrays was used to set the minimum value for that subarray at all time points. To estimate the total number of genes being assayed, all probe set accession numbers present in GenBank were matched against Unigene. A total of 5165 probe sets matched 4097 different Unigene clusters, and the remaining 1182 probe sets, which represent 1175 distinct accession numbers, did not match any Unigene entry.
Multidimensional Scaling
Scaled distances between developmental time points were calculated using a weighted modification of the Pearson correlation coefficient.
where wi is a weighting factor for the number of times a given gene is represented (wi = psi-1, where psi is the total number of probe sets that detect the same gene as probe set i), s is the total number of probe sets measured, xi and yi are expression levels of probe set i from respective time points, and
and
are the scaled mean values, e.g.
Multidimensional scaling coordinates were calculated using the "R" software package, and the resulting projection into three-dimensional space was visualized using Blender (Not a Number B.V., Amsterdam, The Netherlands).
Cluster Analysis
Expression profiles with a relative difference greater than 3-fold and an absolute difference greater than 100 between lowest and highest expression during development were selected for clustering. For comparison, the average expression value of all probe sets at a given time point ranged from 93.6134.4. A total of 1,312 genes met both the relative and absolute difference criteria and were selected for cluster analysis. Individual developmental profiles were normalized such that mean expression = 0 and variance = 1. The resulting data were clustered by a two-dimensional SOM using GeneCluster software (2). An 8-by-5 grid was chosen empirically to balance the goals of minimizing the variance of individual clusters while minimizing the risk of splitting closely associated profiles into multiple clusters. A list of genes assigned to each of the 40 clusters is published in the supporting data.
Identification and Statistical Analysis of Gene Categories
SWISS-PROT version 38 and TrEMBL (including updates before 20-Aug-1999) were used for the purposes of this analysis (62). For 2,665 expressed sequence tags, Affymetrix provided the nearest SWISS-PROT accession number, and this was used to extract keyword annotations from SWISS-PROT 38. For remaining sequences, the accession number was searched against SWISS-PROT to find an explicit match. If no such match existed, the nearest SWISS-PROT/TrEMBL sequences were determined by alignment using BLAST 2.0 (63). If the closest match had greater than 95% degree of identity, its keywords were extracted. If the closest match was an unannotated TrEMBL entry and an annotated SWISS-PROT entry also showed greater than 95% homology, the SWISS-PROT keywords were used. If the highest degree of identity was 8595%, the closest match was used regardless of number of keyword annotations. The resulting gene categories were then tested for statistically significant association with SOM clusters (64). The cumulative hypergeometric distribution was used to determine the probability of chance association between a given cluster and a given gene category. To account for the independent testing of 394 categories that appear more than one time in the set of clustered genes, this value was required to be
 |
Visual inspection of the 8 by 5 SOM revealed that several adjacent clusters have similar expression profiles (Fig. 2
). As such, we reasoned that some gene categories might not reach the threshold of statistical significance due to the partitioning of genes with similar expression profiles into adjacent, related clusters. Furthermore, because the precise location and organization of SOM clusters depends on the random seeding of points used in the analysis, subtle shifts in cluster boundaries among independent runs could theoretically impair the ability of this analytical approach to detect biologically relevant associations. Therefore, more sensitive detection of associations was achieved by pooling overlapping sets of significant genes from 30 independent clustering runs. Briefly, for each individual run with a statistically significant (as defined above) association between some cluster and a specific gene category, a list of genes within that cluster and category was generated. Overlapping contigs of the resulting gene lists were generated, and these final groups of genes represent connected regions of expression space associated with a given gene category. A single SOM was chosen for reference purposes, and a cluster within this map that contains 40% or more of the grouped genes was designated the primary cluster for this association. Similarly, any cluster containing 25%40% of the genes was designated as a secondary cluster. To confirm the statistical significance of associations identified by this approach, gene identifications were randomly permuted with expression data to generate 100 data sets for statistical characterization. Each permutation was analyzed as described for the original data set, and categories that exceeded the statistical significance threshold were tabulated. The resulting data were fit to binomial distributions and indicated that a conservative estimate for the probability of an individual keyword associating by chance with one or more clusters of gene expression was less than 0.03. For gene categories containing less than 10 members, the probability of association was individually estimated.
Statistically significant associations were individually examined for redundant gene representation, and significance was recalculated where appropriate. Categories that no longer exceeded the defined statistical threshold after removal of duplicate entries were removed from further consideration. Additionally, significant associations between clusters and muscle-specific gene categories were considered to represent a variable admixture of mammary gland and skeletal muscle across development and were therefore eliminated from further consideration.
Morphological Analysis
For mammary gland whole mounts and in situ hybridization, no. 4 mammary glands were harvested from FVB mice. Whole mounts were generated by mounting the gland and staining with carmine alum as previously described (12). In situ hybridization of sections from paraffin-embedded glands was performed as previously described (60).
Experimental Animals
All animal experimentation described was conducted in accord with accepted standards of humane animal care. Protocols for animal work were approved by the University of Pennsylvania institutional committee on animal care.
 |
ACKNOWLEDGMENTS
|
---|
 |
FOOTNOTES
|
---|
This research was supported in part by NIH Grants CA-92910 and PO1-CA-77596 from the National Cancer Institute, U.S. Army Breast Cancer Research Program Grants DAMD17-98-1-8226, DAMD17-01-1-0364, DAMD17-98-1-8227 (to C.M.D.), DAMD17-00-1-0401 (to S.E.M.) and grants from the Concert for the Cure.
Abbreviations: CoA, Coenzyme A; Ob, leptin; PGC-1, peroxisome proliferator activated receptor
coactivator 1; PPAR
, peroxisome proliferator activated receptor
; SOM, self-organizing map; Ucp1, uncoupling protein 1; Vegf, vascular endothelial growth factor.
Received for publication February 12, 2002.
Accepted for publication March 18, 2002.
 |
REFERENCES
|
---|
-
Eisen MB, Spellman PT, Brown PO, Botstein D 1998 Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci USA 95:1486314868[Abstract/Free Full Text]
-
Tamayo P, Slonim D, Mesirov J, Zhu Q, Kitareewan S, Dmitrovsky E, Lander ES, Golub TR 1999 Interpreting patterns of gene expression with self-organizing maps: methods and application to hematopoietic differentiation. Proc Natl Acad Sci USA 96:29072912[Abstract/Free Full Text]
-
Iyer VR, Eisen MB, Ross DT, Schuler G, Moore T, Lee JC, Trent JM, Staudt LM, Hudson Jr J, Boguski MS, Lashkari D, Shalon D, Botstein D, Brown PO 1999 The transcriptional program in the response of human fibroblasts to serum. Science 283:8387[Abstract/Free Full Text]
-
Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, Bloomfield CD, Lander ES 1999 Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science 286:531537[Abstract/Free Full Text]
-
Alizadeh AA, Eisen MB, Davis RE, Ma C, Lossos IS, Rosenwald A, Boldrick JC, Sabet H, Tran T, Yu X, Powell JI, Yang L, Marti GE, Moore T, Hudson J, Lu L, Lewis DB, Tibshirani R, Sherlock G, Chan WC, Greiner TC, Weisenburger DD, Armitage JO, Warnke R, Levy R, Wilson W, Grever MR, Byrd JC, Botstein D, Brown PO, Staudt LM 2000 Distinct types of diffuse large B-cell lymphoma identified by gene expression profiling. Nature 403:503511[CrossRef][Medline]
-
Perou CM, Jeffrey SS, van de Rijn M, Rees CA, Eisen MB, Ross DT, Pergamenschikov A, Williams CF, Zhu SX, Lee JC, Lashkari D, Shalon D, Brown PO, Botstein D 1999 Distinctive gene expression patterns in human mammary epithelial cells and breast cancers. Proc Natl Acad Sci USA 96:92129217[Abstract/Free Full Text]
-
Ross DT, Scherf U, Eisen MB, Perou CM, Rees C, Spellman P, Iyer V, Jeffrey SS, Van de Rijn M, Waltham M, Pergamenschikov A, Lee JC, Lashkari D, Shalon D, Myers TG, Weinstein JN, Botstein D, Brown PO 2000 Systematic variation in gene expression patterns in human cancer cell lines. Nat Genet 24:227235[CrossRef][Medline]
-
White KP, Rifkin SA, Hurban P, Hogness DS 1999 Microarray analysis of Drosophila development during metamorphosis. Science 286:21792184[Abstract/Free Full Text]
-
Neville MC, Daniel CW, eds. 1987 The mammary gland: development, regulation, and function. New York: Plenum Press
-
Daniel CW, Silberstein GB 1987 Postnatal development of the rodent mammary gland. In: Neville MC, Daniel CW, eds. The mammary gland: development, regulation, and function. New York: Plenum Press; 336
-
Bittner M, Meltzer P, Chen Y, Jiang Y, Seftor E, Hendrix M, Radmacher M, Simon R, Yakhini Z, Ben-Dor A, Sampas N, Dougherty E, Want E, Marincola F, Gooden C, Lueders J, Glatfelter A, Pollock P, Carpten J, Gillanders E, Leja D, Dietrich K, Beaudry C, Berens M, Alberts D, Sondak V, Hayward N, Trent J 2000 Molecular classification of cutaneous malignant melanoma by gene expression profiling. Nature 406:536540[CrossRef][Medline]
-
Gardner HP, Belka GK, Wertheim GBW, Hartman JL, Ha SI, Gimotty PA, Marquis ST, Chodosh LA 2000 Developmental role of the SNF1-related kinase Hunk in pregnancy-induced changes in the mammary gland. Development 127:44934509[Abstract/Free Full Text]
-
Kearsey SE, Labib K 1998 MCM proteins: evolution, properties, and role in DNA replication. Biochim Biophys Acta 1398:113136[Medline]
-
Derenzini M, Sirri V, Trere D, Ochs RL 1995 The quantity of nucleolar proteins nucleolin and protein B23 is related to cell doubling time in human cancer cells. Lab Invest 73:497502[Medline]
-
Rapoport TA, Matlack KE, Plath K, Misselwitz, B, Staeck O 1999 Posttranslational protein translocation across the membrane of the endoplasmic reticulum. Biol Chem 380:11431150[Medline]
-
Greasley PJ, Bonnard C, Amati B 2000 Myc induces the nucleolin and BN51 genes: possible implications in ribosome biogenesis. Nucleic Acids Res 28:446453[Abstract/Free Full Text]
-
Neville MC, Picciano MF 1997 Regulation of milk lipid secretion and composition. Annu Rev Nutr 17:159183[CrossRef][Medline]
-
Bansal MP, Medina D 1993 Expression of fatty acid-binding proteins in the developing mouse mammary gland. Biochem Biophys Res Commun 191:6169[CrossRef][Medline]
-
Todorov PT, McDevitt TM, Meyer DJ, Ueyama H, Ohkubo I, Tisdale MJ 1998 Purification and characterization of a tumor lipid-mobilizing factor. Cancer Res 58:23532358[Abstract]
-
Lund LR, Romer J, Thomasset N, Solberg H, Pyke C, Bissell MJ, Dano K, Werb Z 1996 Two distinct phases of apoptosis in mammary gland involution: proteinase- independent and -dependent pathways. Development 122:181193[Abstract/Free Full Text]
-
Furth PA 1999 Introduction: mammary gland involution and apoptosis of mammary epithelial cells. J Mammary Gland Biol Neoplasia 4:123127[CrossRef][Medline]
-
Gigliotti AP, DeWille JW 1998 Lactation status influences expression of CCAAT/enhancer binding protein isoforms mRNA in the mouse mammary gland. J Cell Physiol 174:232239[CrossRef][Medline]
-
Roberts LR, Adjei PN, Gores GJ 1999 Cathepsins as effector proteases in hepatocyte apoptosis. Cell Biochem Biophys 30:7188[Medline]
-
Benaud C, Dickson RB, Thompson EW 1998 Roles of the matrix metalloproteinases in mammary gland development and cancer. Breast Cancer Res Treat 50:97116[CrossRef][Medline]
-
Gregoire FM, Smas CM, Sul HS 1998 Understanding adipocyte differentiation. Physiol Rev 78:783809[Abstract/Free Full Text]
-
Danesch U, Hoeck W, Ringold GM 1992 Cloning and transcriptional regulation of a novel adipocyte-specific gene, FSP27. CAAT-enhancer-binding protein (C/EBP) and C/EBP-like proteins interact with sequences required for differentiation-dependent expression. J Biol Chem 267:71857193[Abstract/Free Full Text]
-
Taylor-Papadimitriou J, Lane EB 1987 Keratin expression in the mammary gland. In: Neville MC, Daniel CW, eds. The mammary gland: development, regulation, and function. New York: Plenum Press; 181215
-
Smith GH, Mehrel T, Roop, DR 1990 Differential keratin gene expression in developing, differentiating, preneoplastic, and neoplastic mouse mammary epithelium. Cell Growth Differ 1:161170[Abstract]
-
Cherayil BJ, Weiner SJ, Pillai S 1989 The Mac-2 antigen is a galactose-specific lectin that binds IgE. J Exp Med 170:19591972[Abstract]
-
Atkinson PG, Blackwell JM, Barton CH 1997 Nramp1 locus encodes a 65 kDa interferon-
-inducible protein in murine macrophages. Biochem J 325:779786[Medline]
-
McKnight AJ, Gordon S 1998 The EGF-TM7 family: unusual structures at the leukocyte surface. J Leukoc Biol 63:271280[Abstract]
-
Belaaouaj A, Shipley JM, Kobayashi DK, Zimonjic DB, Popescu N, Silverman GA, Shapiro SD 1995 Human macrophage metalloelastase. Genomic organization, chromosomal location, gene linkage, and tissue-specific expression. J Biol Chem 270:1456814575[Abstract/Free Full Text]
-
Spilsbury K, OMara MA, Wu WM, Rowe PB, Symonds G, Takayama Y 1995 Isolation of a novel macrophage- specific gene by differential cDNA analysis. Blood 85:16201629[Abstract/Free Full Text]
-
Strange R, Li F, Saurer S, Burkhardt A, Friis RR 1992 Apoptotic cell death and tissue remodeling during mouse mammary gland involution. Development 115:4958[Abstract]
-
Fadok VA 1999 Clearance: the last and often forgotten stage of apoptosis. J Mammary Gland Biol Neoplasia 4:203211[CrossRef][Medline]
-
Bazzoni G, Dejana E, Lampugnani MG 1999 Endothelial adhesion molecules in the development of the vascular tree: the garden of forking paths. Curr Opin Cell Biol 11:573581[CrossRef][Medline]
-
Risau W 1997 Mechanisms of angiogenesis. Nature 386:671674[CrossRef][Medline]
-
Zangani D, Darcy KM, Masso-Welch PA, Bellamy ES, Desole MS, Ip MM 1999 Multiple differentiation pathways of rat mammary stromal cells in vitro: acquisition of a fibroblast, adipocyte or endothelial phenotype is dependent on hormonal and extracellular matrix stimulation. Differentiation 64:91101[CrossRef][Medline]
-
Bieber LL 1988 Carnitine. Annu Rev Biochem 57:261283[CrossRef][Medline]
-
Salway JG 1994 Metabolism at a glance. Oxford: Blackwell Science
-
Ricquier D, Bouillaud F 2000 The uncoupling protein homologues: UCP1, UCP2, UCP3, StUCP and AtUCP. Biochem J 345:161179[CrossRef][Medline]
-
Hoshino K 1967 Transplantability of mammary gland in brown fat pads of mice. Nature 213:194195[Medline]
-
Elias JJ, Pitelka DR, Armstrong RC 1973 Changes in fat cell morphology during lactation in the mouse. Anat Rec 177:533548[Medline]
-
Bani-Sacchi T, Bianchi S, Bani G, Bigazzi M 1987 Ultrastructural studies of white adipocyte differentiation in the mouse mammary gland following estrogen and relaxin. Acta Anat (Basel) 129:19[Medline]
-
Sakakura T 1991 New aspects of stroma-parenchyma relations in mammary gland differentiation. Int Rev Cytol 125:165203[Medline]
-
Hovey RC, McFadden TB, Akers RM 1999 Regulation of mammary gland growth and morphogenesis by the mammary fat pad: a species comparison. J Mammary Gland Biol Neoplasia 4:5368[CrossRef][Medline]
-
Guerra C, Koza RA, Walsh K, Kurtz DM, Wood PA, Kozak LP 1998 Abnormal nonshivering thermogenesis in mice with inherited defects of fatty acid oxidation. J Clin Invest 102:17241731[Abstract/Free Full Text]
-
Barberá MJ, Schlüter A, Padraza N, Inglesias R, Villarroya F, Giralt M 2001 Peroxisome proliferator-activated receptor
activates transcription of the brown fat uncoupling protein-1 gene. J Biol Chem 276:14861493[Abstract/Free Full Text]
-
Aoyama T, Peters JM, Iritani N, Nakajima T, Furihata K, Hashimoto T, Gonzalez FJ 1998 Altered constitutive expression of fatty acid-metabolizing enzymes in mice lacking the peroxisome proliferator-activated receptor
(PPAR
). J Biol Chem 273:56785684[Abstract/Free Full Text]
-
Desvergne B, Wahli W 1999 Peroxisome proliferator-activated receptors: nuclear control of metabolism. Endocr Rev 20:649688[Abstract/Free Full Text]
-
Puigserver P, Wu Z, Park CW, Graves R, Wright M, Spiegelman BM 1998 A cold-inducible coactivator of nuclear receptors linked to adaptive thermogenesis. Cell 92:829839[Medline]
-
Vega RB, Huss JM, Kelly DP 2000 The coactivator PGC-1 cooperates with peroxisome proliferator-activated receptor
in transcriptional control of nuclear genes encoding mitochondrial fatty acid oxidation enzymes. Mol Cell Biol 20:18681876[Abstract/Free Full Text]
-
Nedergaard J, Connolly E, Canon B 1986 Brown adipose tissue in mammalian neonate. In: Trayhurn P, Nicholls DG, eds. Brown adipose tissue. London: Edward Arnold; 152213
-
Girard J, Ferré P, Pégorier J-P, Duée P-H 1992 Adaptations of glucose and fatty acid metabolism during perinatal period and suckling-weaning transition. Physiol Rev 72:507562[Free Full Text]
-
Cousin B, Cinti S, Morroni M, Raimbault S, Ricquier D, Pénicaud L, Casteilla 1992 Occurrence of brown adipocytes in rat white adipose tissue: molecular and morphological characterization. J Cell Sci 103:931942[Abstract/Free Full Text]
-
Tsukiyama-Kohara K, Poulin F, Kohara M, DeMaria CT, Cheng A, Wu Z, Gingras AC, Katsume A, Elchebly M, Spiegelman BM, Harper ME, Tremblay ML, Sonenberg N 2001 Adipose tissue reduction in mice lacking the translational inhibitor 4E-BP1. Nat Med 7:11281132[CrossRef][Medline]
-
Obregon MJ, Jacobsson A, Kirchgessner T, Schotz MC, Cannon B, Nedergaard J 1989 Postnatal recruitment of brown adipose tissue is induced by the cold stress experienced by the pups. An analysis of mRNA levels for thermogenin and lipoprotein lipase. Biochem J 259:341346[Medline]
-
Asano A, Morimatsu M, Nikami H, Yoshida T, Saito M 1997 Adrenergic activation of vascular endothelial growth factor mRNA expression in rat brown adipose tissue: implication in cold-induced angiogenesis. Biochem J 328:179183[Medline]
-
Bonilla S, Redonnet A, Noel-Suberville C, Pallet V, Garcin H, Higueret P 2000 High-fat diets affect the expression of nuclear retinoic acid receptor in rat liver. Br J Nutr 83:665671[Medline]
-
Marquis ST, Rajan JV, Wynshaw-Boris A, Xu J, Yin GY, Abel KJ, Weber BL, Chodosh LA 1995 The developmental pattern of Brca1 expression implies a role in differentiation of the breast and other tissues. Nat Genet 11:1726[Medline]
-
Lockhart DJ, Dong H, Byrne MC, Follettie MT, Gallo MV, Chee MS, Mittmann M, Wang C, Kobayashi M, Horton H, Brown EL 1996 Expression monitoring by hybridization to high-density oligonucleotide arrays. Nat Biotechnol 14:16751680[Medline]
-
Bairoch A, Apweiler R 1999 The SWISS-PROT protein sequence data bank and its supplement TrEMBL in 1999. Nucleic Acids Res 27:4954[Abstract/Free Full Text]
-
Altschul SF, Madden TL, Schaffer AA, Zhang Z, Miller W, Lipman DJ 1997 Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25:33893402[Abstract/Free Full Text]
-
Tavazoie S, Hughes JD, Campbell MJ, Cho RJ, Church GM 1999 Systematic determination of genetic network architecture. Nat Genet 22:281285[CrossRef][Medline]