|
|
||||||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
1 Department of Physiology and Biophysics, University of Colorado Health Sciences Center at Fitzsimmons, Aurora, Colorado
2 Department of Obstetrics and Gynecology, University of Colorado Health Sciences Center at Fitzsimmons, Aurora, Colorado
3 Department of Anesthesiology, University of Colorado Health Sciences Center at Fitzsimmons, Aurora, Colorado
4 Department of Pharmacology, University of Colorado Health Sciences Center at Fitzsimmons, Aurora, Colorado
5 Department of Pathology, University of Colorado Health Sciences Center at Fitzsimmons, Aurora, Colorado
6 Division of Cancer Sciences and Molecular Pathology, Western Infirmary, Glasgow, United Kingdom
| ABSTRACT |
|---|
|
|
|---|
lipid synthesis; microarray; metabolomics; dietary lipid; magnetic resonance spectroscopy
| INTRODUCTION |
|---|
|
|
|---|
19 days including an intensive proliferative phase followed by a differentiation phase marked by an increase in milk protein gene expression, lipid droplet formation, and stromal adaptations (42, 43). A fall in progesterone
day 18 initiates secretory activation, a programmed series of changes in the epithelium that leads to the copious secretion of very rich milk consisting of
12% protein, 30% lipid, and 5% lactose. Lipid synthesis is particularly remarkable: The mammary gland of the FVB mouse must synthesize an amount of triacylglycerol (TAG) equivalent to the entire weight of the mouse during a 20-day lactation, generally while the mouse is eating a diet containing <8% of the calories as fat (37). Coordinate transcriptional regulation of many of the lipid synthesis enzymes occurs in other organs such as the liver (12, 17) and adipose tissue (39) as well as pancreatic ß-cells (2). We hypothesized that coordinate regulation of the same pathways might be integral to the remarkable increase in lipid and lactose synthesis during the activation of secretion (sometimes called lactogenesis II) in the mammary gland. We wished to identify both the key enzymes that change during the initiation of this program and the transcriptional regulators involved.
Our strategy was first to utilize microarray technology to analyze gene expression in the mammary glands of FVB mice during secretory activation. Having organized the metabolic genes by functional category, we identified genes whose expression changed significantly between the critical time points of pregnancy day 17 (P17) and lactation day 2 (L2). One advantage of these time points is the adipocyte compartment is quite small (41), so that we have some confidence that changes seen occur in the epithelial compartment. In addition, we have previously shown that adipocyte genes have a specific signature (42) that allows them to be identified as such.
As a reference we compared expression levels in the mammary gland to those in the liver at mature lactation (day 9). The synthetic picture that emerged from this analysis was then amplified by quantitating relevant metabolic intermediates using nuclear magnetic resonance spectroscopy (MRS). Finally, we identified potential regulatory genes by further mining both the original data and data from an experiment in which dietary lipid was switched from 8 to 40% (in kcal) in lactating mice. This approach provided a global picture of the multiple points at which the transcriptome of the enzymes of lipid synthesis is regulated and produced several candidate regulatory molecules. We found that mammary gene expression is altered in a unique manner to provide a framework for a balanced upregulation of lactose and lipid synthesis at secretory activation.
| EXPERIMENTAL PROCEDURES |
|---|
|
|
|---|
Isolation and Validation of Total RNA
Total RNA was isolated and purified following the Qiagen RNA extraction/clean-up protocol (including DNase treatment). Purity, concentration, and integrity of total RNA from each sample were verified using a NanoDrop spectrophotometer (NanoDrop Technologies, Wilmington, DE) and the RNA 6000 Nano Assay (Agilent Technologies, Palo Alto, CA).
Microarray Analysis
Total RNA was acquired, verified as above, amplified, biotin labeled (Enzo, Farmingdale, NY), and fragmented following the 2002 protocol for eukaryotic target preparation (Affymetrix, Santa Clara, CA). Acceptable samples were hybridized to Affymetrix Mu74Av2 microarray chips. Raw data were gathered from scanned array chips using Affymetrix Microarray Suite (MAS) version 5.0. Quality control chip statistics were evaluated from each array, including background intensity, noise, scaling factor, 3'5' ratios of the "housekeeping" genes, % distribution of calls, and presence of "spiked" control genes. Compiled data in the form of 28 individual CEL files, the primary output of scanned Mu74Av2 microarray chips using the contemporary GeneChip Operating Software v1.1 (Affymetrix, Santa Clara, CA), were imported to GeneSpring (Silicon Genetics, Redwood City, CA) for analysis using the native probe level GC-Robust Multi-array Average (GC-RMA) algorithm. The complete results of this analysis are available at http://www.ncbi.nlm.nih.gov/geo under the series record GSE4222.
Data Analysis
Genes were removed that scored below 15.0 RAW intensity data units in four or more of the conditions analyzed. Genes that were expressed within 15% of their mean across all time points were also removed giving a filtered list of 6,560 genes. To define genes that vary significantly over either the time course or the diet experiment, the Kruskal-Wallis nonparametric statistical test was run on the 6,560 gene list using a false discovery rate of 0.05.
Quantitative Real-time RT-PCR
We utilized real-time RT-PCR (QRT) to amplify the RNA for 15 lipid synthesis genes as well as forms of the sterol regulatory element binding proteins not available on the chip and four genes commonly used as references. In some samples a multiplexing strategy was used with MHC gene ß2-microglobulin. Copy number per initial 1 µg of RNA was calculated from the concentration of the amplicons used to generate the standard curve. Because the differences between results obtained at P17 and P18 were, in general, not significant (see Table S3A), the data were combined for comparison with P17 data from array analysis.
Primer and Probe Design
Commercial software and web-based algorithms were used in the design of all primer and probe sets, TaqMan Primer Express Applied Biosystems software (http://www.appliedbiosystems.com/) and Integrated DNA Technologies (http://biotools.idtdna.com/primerquest/), respectively. Generally, we constructed a primer/probe construct to target a splice boundary within a message sequence so that genomic DNA was not amplified. All primers, labeled probes, and synthetic amplicon sequences were purchased from either Applied Biosystems or Integrated DNA Technologies. Sequences are given in Supplementary Table S3B. The primers and probes for independent targets were combined at proper concentrations prior to addition of the PCR master mix for both the gene investigated and ß2-microglobulin. For each investigated gene in the 25-µl total volume reaction, a 20-µl master mix was prepared, including 12.5 µl PCR mix, 2.5 µl GeneA primer/probe solution (900 nM forward, 200 nM probe, 900 nM reverse), 0.25 µl of 100 units/µl AmpErase-UNG (a uracil-N-glycosylase that prevents carry-over amplification from templates containing deoxyuridine, i.e., RNA), 0.125 µl of 50 units/µl AmpliTaqGold polymerase. This mix with 2.125 µl of nuclease-free water was loaded to the 96-well plate. To the 20 µl of master mix, 5 µl of 1:10 diluted cDNA template was added. cDNA was synthesized in a 20-µl reaction volume using 1.0 µg of total RNA in 8.5 µl of nuclease free water, 2.5 µl of random hexamers (50 µM), 4 µl of 5x 1st-strand buffer, 2 µl DTT (0.1 M), 1 µl dNTP blend (10 mM), 1 µl RNase inhibitor (40 U/µl), and 1 µl superscript III (200 u/µl); reagents were incubated at 42 degrees for 60 min and enzymes were inactivated at 75°C degrees for 15 min. All real-time PCR data were captured using an Opticon Monitor II real-time PCR detection system plus the accompanying software version 2.02.24 (MJ Research, Bio-Rad Laboratories, Hercules, CA). Copy numbers of the investigated genes were calculated from the respective standard curves using the formula 10^[(Ct minus Y intercept)/(slope)], where Ct is cycle time. Finally, for each sample the ratio of the copy number of the investigated gene to the copy number of the 18S RNA times 106 produced normalized data.
Mammary Gland Extraction for Metabolomic Analysis
For MRS analysis, snap-frozen mammary glands (
0.2 g) were extracted using the dual-phase methanol-chloroform procedure. Briefly, frozen tissues were weighed, powdered in the presence of liquid nitrogen, and homogenized in 4 ml of ice-cold methanol. Ice-cold chloroform (4 ml) was added, samples vortexed, and an additional 4 ml of ice-cold water were added. After storage at 4°C overnight, they were centrifuged for 30 min at 1,400 g at 4°C. The upper methanol phase (water-soluble extract) was transferred into a lyophilizing glass and freeze-dried overnight then dissolved in 0.5 ml of deuterium oxide (D2O) prior to [1H]MRS analysis.
Quantitative MRS Metabolic Analysis
All one- and two-dimensional [1H]MR spectra were obtained on a Bruker 500 MHz DRX spectrometer (Bruker Bisopin, Fremont, CA) using an inverse 5-mm TXI probe. For metabolite identification in water-soluble mammary gland extracts, a two-dimensional (2D)-H, C-heteronuclear single quantum correlation technique was used. The experiments were acquired with 512 increments and 160 scans per increment, using 90° pulse and a recovery delay of 0.5 s. The spectral width was 10 ppm in the proton dimension and 140 ppm in the 13C dimension. Lactate (Lac3, CH3) was used as a chemical shift reference for both carbon (21 ppm) and proton (1.32 ppm) axes. Tissue metabolites were identified based on the results from our chemical shift data base. For metabolite quantification, a standard water presaturation pulse program "zgpr" was used to suppress water residue signal (44). The total number of acquisitions was 40, the pulse delay of 12.8 s was applied between acquisitions for fully relaxed [1H]NMR spectra (calculated as 5*T1). An external standard, trimethylsilylpropionic acid (TMSP-d4, 0.6 mmol/l) was used for metabolite quantification and as a chemical shift reference (0 ppm). After performing Fourier transformation and making phase and baseline corrections, we integrated each identified 1H peak using 1D WINNMR program (Bruker Biospin). The absolute concentrations of single metabolites were then referred to the TMSP integral and calculated according to the Eq. 1:
![]() | (1) |
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
62% of the values obtained with microarray probably due to differences in normalization. We obtained expression levels of three genes that are often used for normalizing RNA data in the mammary gland, GAPDH, cytokeratin 18, ß2-macroglobulin, as well as the 18S ribosomal RNA (see Table S3A). Because the 18S ribosomal RNA varied the least we used it for normalization, slightly reducing the variance between individual samples. The dilution of the mRNA by milk protein mRNA is likely to further dilute the values obtained from the lactating gland, probably explaining at least part of the observed decrease in GAPDH, keratin 18, and ß2-macroglobulin at lactation when normalized to S18. In the following discussion, we make a notation where marked differences between either the Stein data or the QRT data were found.
|
160 transporters for these substances are available on the Affymetrix chip used, we catalogued only transporters for those substances with reasonably high levels of mammary expression (raw expression values >40) in Fig. 1. Glucose transporters that are present from the expression data are GLUT-1 (SLC2a1), GLUT-4 (SLC2a4), GLUT-8 (SLC2a8), and the sodium-dependent glucose transporter SLC5a1. GLUT-1 and SLC5a1 mRNAs were present at much higher levels in the mammary gland than the liver, and GLUT-1 gene expression increased at least 1.7-fold after parturition by both array and real-time RT-PCR (Fig. 3, Table S3C). GLUT-4 (SLC2a4) mRNA is also expressed in the mammary gland, but it is localized in the adipose compartment and shows the decrease during pregnancy and lactation expected of adipose genes (data not shown and Ref. 41). GLUT-12, known to be expressed in the mammary gland (30), is not present on the Mu74Av2 chip. Transcripts for four genes associated with amino acid transport were identified from the array data. SLC1a4, a system alanine-serine cysteine (ASC) transporter, is a sodium-coupled transporter with specificity for alanine, serine, and cysteine capable of producing high gradients of these amino acids between the cell and the interstitial fluid (24). The L-system neutral amino acid transporter consists of a heteromer of SLC3a2 (4F2hc) and SLC7a5 (LAT1) with a broad specificity for neutral amino acids. It is thought to operate in the exchange mode only (36), possibly using amino acid gradients produced by the ASC and other transporters to increase cellular concentrations of such amino acids as valine, serine, and threonine. SLC1a4 and SLC7a5 mRNAs were upregulated 1.5-fold or more at the initiation of lactation and are expressed at levels many times higher than the liver, suggesting mammary specificity.
Genes for three members of the SLC27 family of fatty acid transporters are expressed in the mammary gland: SLC27a1, SLC27a3, and SLC27a4. SLC27a1 shows a temporal expression pattern similar to adipose tissue genes (data not shown)(42). SLC27a3 is expressed at a level about eight times that in the liver and, as we shall see, is subject to dietary regulation. SLC27a4 is expressed at a level similar to that of the liver, but its mRNA is upregulated at least 1.5-fold at the onset of lactation. CD36 is a multifunctional scavenger receptor that is also thought to transport fatty acids (8); although not upregulated at the onset of lactation, it is expressed at a level much higher than the liver.
Together these data suggest that expression of the RNA for three classes of transporters is significantly increased at the onset of lactation, glucose, amino acids, and fatty acids. Since substrate transport into the cell likely represents a potent regulation point for synthesis of milk components, this increase in gene expression is appropriate.
Lactose synthesis.
Lactose synthesis from glucose and UDP-galactose takes place in the trans-Golgi catalyzed by lactose synthase, an enzyme complex made of up the 1,4 ß-galactosyl transferase (ß1,4 GalT) and the mammary-specific protein
-lactalbumin. The presence of the cofactor,
-lactalbumin, is rate limiting for lactose synthesis (47), and its gene expression increases 5- to 15-fold between P17 and L2 (Figs. 1, 3; Table S. 3C).
Glycolysis.
Surprisingly, mRNA levels for the enzymes conducting the initial two steps of glycolysis, hexokinase and glucose phosphate isomerase (GPI), were significantly downregulated at the initiation of lactation and GPI was expressed at levels significantly lower than the liver (Figs. 1 and 2). Because the Km of lactose synthase for glucose is in the millimolar range (26), a decrease in hexokinase activity would effectively partition glucose toward lactose production. A decrease in GPI would help shunt glucose-6-PO4 to the pentose phosphate shunt (see below; Fig. 2). Phosphofructokinase, aldolase C, and pyruvate kinase were expressed at levels much higher than their liver counterparts, probably because the isoforms present in the mammary gland differ from those in the liver. Aldolase C is expressed mainly in the brain where it likely contributes to myelin synthesis (1); its substantial upregulation by both array and PCR analysis at the onset of lactation is consistent with previous observations that the enzyme localizes to mammary alveoli and is responsive to prolactin (31). It is likely important in generation of the glycerol-phosphate backbone for TAG synthesis.
Pentose phosphate shunt.
In lipid synthesizing tissues, the pentose phosphate shunt provides a substantial fraction of the reducing equivalents needed for NADPH synthesis. From our data this increase in activity would be accomplished both by downregulating GPI (see above) and possibly increasing expression of the gene for glucose-6-phosphate dehydrogenase (Figs. 1 and 2), although a large increase was not confirmed in either the Stein data set or by QRT; mRNA expression for other enzymes in this pathway, phosphogluconase dehydrogenase was increased significantly in both arrays, and this and other enzymes and this as well as transketolase and transaldolase were expressed at levels three times those in the liver. At secretory activation in the rat mammary gland the activity of this pathway was increased 22 times (9). The fact that none of the increases in expression for the enzymes in this pathway was close to 22-fold, suggests that the increased flux is likely also due increased substrate availability or other modes of regulation.
Gluconeogenesis.
The data in Fig. 1 confirm the expectation that gluconeogenesis is not a function of the murine mammary epithelium. Mammary expression of two of the key enzymes of gluconeogenesis, PEPCK and fructose-bis-phosphatase, is <2% of that of the liver, suggesting that carbon backbones transported into the mammary gland in the form of amino acids are not significantly converted into glucose in this organ, a point of importance when we consider citrate synthesis.
Citric acid cycle.
Gene expression for many enzymes of the citric acid cycle is significantly downregulated at the onset of lactation; these genes are also are expressed at levels significantly lower than the liver (Fig. 1). The exceptions are pyruvate carboxylase, some components of pyruvate dehydrogenase, and citrate synthase. All are expressed at levels at least threefold higher than the liver, and most are significantly upregulated at secretory activation. In addition, there is an impressive increase in the tricarboxylic acid (citrate) transporter, three- to fourfold, by both microarray and real-time RT-PCR (Figs. 1 and 3). If these adaptations are translated to the protein level they would serve to focus mitochondrial activity on the generation and export of citrate to the cytoplasm where it would be converted to substrates for fatty acid and cholesterol synthesis as well as for the generation of NADPH via the malate shuttle.
Management of plasma fatty acids.
Plasma TAG are hydrolyzed in the capillary by lipoprotein lipase, and their products are transferred across the cell membrane by fatty acid and glycerol transporters, thought to be specific aquaporins. Unfortunately these aquaporins were not represented on the chip. Lipoprotein lipase has long been known to be expressed at very high levels in the lactating mammary gland (15), and the array results are consistent with this finding. Once hydrolysis is complete, lipid transporters, described above, transfer the fatty acids to the cytoplasm where they must be bound to fatty acid binding proteins (FABP) or acylated and bound to fatty acyl binding proteins. As shown in Fig. 1, the mRNA for FABP3, the heart-type binding protein, is expressed at very high levels in the mammary gland (4).
Acylated fatty acids can be metabolized by the enzymes of the ß-oxidation pathway or incorporated into TAG for secretion into milk. In the mammary gland the mRNA for most enzymes of the ß-oxidation pathway falls two- to threefold at secretory activation (Fig. 1, Table S3C) directing fatty acids to TAG synthesis. One unexpected finding was that cytoplasmic thioesterase 1 (Cte-1), an epithelial enzyme that specifically hydrolyzes long-chain saturated acyl-CoAs (20), is expressed at very high levels in pregnancy (see Table S1), falling nearly 26-fold at the onset of lactation. In the liver acyl-CoA thioesterases are thought to be critical for maintaining the balance between energy storage (TAG synthesis) or energy utilization (ß-oxidation) (21). The fall in the mRNA for Cte-1 suggests at least a gate-keeping function of this enzyme for entry of fatty acids into the ß-oxidation pathway; in addition the fall would prevent hydrolysis of acyl-CoAs, which are then available for TAG synthesis.
Fatty acid synthesis.
Citrate, generated in the mitochondria and transported to the cytoplasmic compartment, is hydrolyzed by ATP citrate lyase to provide substrate for both NADPH and malonyl CoA synthesis (3). As shown in Figs. 1 and 2 and Table S3C, the mRNA of most of the enzymes involved in fatty acid synthesis is upregulated at secretory activation, and the final expression level of most of them is two- to threefold above the liver. ATP citrate lyase in particular was upregulated at least 2.7-fold by microarray and QRT, suggesting that its activity is a major metabolic control point. Acetyl-CoA carboxylase has been identified as a major regulator of fatty acid synthesis in other species, but it was not present on the MU74AV2A chip used in these experiments; preliminary assays using real-time QRT (not shown) showed a 1.5-fold change in gene expression in the murine mammary gland. However, the existence of several uncharacterized isoforms in the mouse make this conclusion unreliable. Fatty acid synthase gene expression was increased at least 1.3-fold by microarray and real-time RT-PCR at the initiation of lactation (Figs. 1 and 3).
The mRNAs for SCD2 and Evolv 1 are expressed at much higher levels in the mammary gland than the liver as is the
5 fatty acid desaturase 1. These enzymes are responsible in part for formation of very-long-chain polyunsaturated fatty acids; the high expression of their mRNAs, along with the higher expression of mRNA for other enzymes for fatty acid synthesis, suggests most of the long-chain polyunsaturated fatty acids present in mouse milk are formed in the mammary gland and not transported from the liver.
TAG synthesis.
To synthesize TAGs, both fatty acyl-CoAs and glycerol 3-phosphate must be available. Glycerol can enter the mammary alveolar cell from the plasma to be phosphorylated by glycerol kinase, expressed at a level three times that of the liver and significantly upregulated at secretory activation in the Stein data set (Table S3C). Prior to TAG synthesis, fatty acids must be acetylated; the most highly expressed of the enzymes that do this is long-chain acyl-CoA synthetase 1, moderately upregulated at the onset of lactation. Finally the acyl-CoA groups must be added to glycerol phosphate. In the mammary gland the most highly expressed enzyme in this process is 1-acylglycerol-3-phosphate O-acyltransferase (AgPAT).
Cholesterol synthesis.
Many of the enzymes of cholesterol synthesis are present on the microarray chip, and a good number were significantly upregulated between P17 and L2 (Figs. 1 and 2), suggesting a high degree of regulation in the mammary gland. However, many fall by L9 (see Table S1), when they are expressed at high levels in the liver. This finding is consistent with the data from goats that only
20% of the cholesterol in milk is synthesized in the mammary gland, the remainder being imported from the plasma (29).
Metabolomic Analysis of Substrate Concentrations
Nuclear MRS was used to assess concentrations of water soluble metabolites in mammary tissue at P17 and L2. Representative MRS spectra, with major identified peaks indicated are shown in Fig. 4A. The spectra, which are clearly different in pregnancy and lactation, were quantified as described in methods, and the results are summarized in the bar graph in Fig. 4B, which shows only metabolites that changed significantly (P < 0.05) at secretory activation. Levels of all the metabolites that could be quantified from these spectra are shown in Table S4; metabolites with concentrations below 10 µM cannot be quantified by this method. Both the
- and ß-isoforms of glucose, galactose, and lactose can be discerned in MRS spectra; however, only the ß-isoforms of the sugars have peaks with differing chemical shifts allowing the concentrations of the individual sugars to be measured. As expected, the concentrations of glucose, lactose, and galactose were all increased in glands from lactating animals. ß-Glucose increased more than twofold between pregnancy and lactation as expected from the increase in transport in the face of a fall in hexokinase. The change is consistent with previous measures of glucose concentration in the lactating gland (27). Both galactose and lactose were nearly undetectable in pregnancy but rose to millimolar levels in the lactating gland. Sialic acid (neuraminic acid or NANA), utilized in the formation of the complex carbohydrates present in milk, was detectable only in the lactating gland. Inositol, which plays a role in the maintenance of osmotic pressure as well as being an important component of phospholipids, also increased significantly.
|
The third group of metabolites in Fig. 4B includes citrate, reflecting the increase in citrate synthesis and transport from the mitochondrion; betaine, a choline derivative important in homocysteine metabolism; as well as creatine and hydroxybutyrate. The final group includes choline, precursor to phosphocholine (PC), and glycerolphosphocholine (GPC), the major choline metabolites in rat milk (40). Both are precursors of the major membrane phospholipid, phosphatidyl choline found in the milk fat globule membrane and other membrane fractions of milk. No changes in lactate were observed between pregnancy and lactation, explained by the rapid utilization of pyruvate for citrate formation.
Analysis of Potential Regulatory Genes
Having determined a number of functionally important genes and metabolites that are regulated in a tissue-specific manner at the onset of lactation, we asked the question: How is the expression of these genes coordinately regulated? Here we must examine two types of regulation: 1) positive upregulators of lipid and lactose synthetic processes and 2) negative regulators whose expression is high in pregnancy to inhibit milk synthesis and which decrease at the onset of lactation to release this repression. For the first category we examined the temporal profiles and expression levels of categories of genes known to regulate lipid metabolism in other tissues including the Akts, C/EBPs, IGFs, PPARs, LXRs, SREBPs, and STATs (17, 25, 46) (Fig. 5). Many of these genes were expressed at levels much lower in the mammary gland than in the liver (shown in blue in MG/liver column), suggesting that they may not have a role in the regulation of metabolism in the lactating mammary gland. The exceptions were AKT1, C/EBP
, LXRß, PPAR
, STAT5a, and STAT5b, all of which were expressed at levels at least threefold above those in the liver. In addition, although the prolactin receptor (PRLR) and sterol regulatory element binding protein-1 (SREBP-1) were expressed at levels approximately equal to those of the liver, these two genes, along with Akt1 and STAT5b, were significantly upregulated at the onset of lactation.
|
For the SREBPs, however, the microarray analysis was inadequate because it didnt distinguish between the splice forms SREBP-1a and SREBP-1c, which have different roles in the fatty acid synthesis; furthermore, SREBP-2 was not represented on the chip. We, therefore, used real-time RT-PCR to quantitate the expression of SREBP family members during the initiation of lactation (Fig. 3), finding that SREBP-1c was the most highly expressed and nearly doubled at the onset of lactation. SREBP-1c has been implicated as a metabolic regulator in many organs including the liver (11), adipocytes (28), muscle (14), and pancreas (2), as well as the mammary gland (5, 6). In pregnancy SREBP-1a was expressed at a level about equal to SREBP-1c and the regulator of cholesterol synthesis, SREBP-2, was expressed at a threefold lower level. However, only SREBP1-c increased at the onset of lactation. During secretory activation, expression of 38 of the 102 metabolic genes examined was upregulated >1.4-fold. Of these genes, 17 corresponded to genes whose activity has been suggested to be upregulated by SREBPs on the basis of gene expression studies in the liver (17, 19). This observation plus the increase in SREBP-1c at the initiation of lactation (Fig. 3 and Table S3A) suggest that SREBP-1c may be a prominent regulator of lipid synthesis in the mammary gland.
A number of the highly expressed metabolic regulators were significantly downregulated at the onset of lactation including AKT2, C/EBP
, PPAR
, STAT3. Of these C/EBP
and STAT3 have been implicated in the apoptotic processes that follow weaning in the mammary gland (10, 49); a fall in these factors in the context of the initiation of lactation suggests protection from apoptosis. The role of AKT2 and PPAR
in the mammary gland is unknown but they are highly expressed in adipose tissue and show expression patterns characteristic of that compartment (data not shown and Ref. 41).
We have mentioned negative regulators above and have previously implicated transforming growth factor (TGF)-ß2, TGF-ß3, Wnt5b, and IGFBP5 (41). As shown in Fig. 5, all were downregulated in this analysis; however, none has been linked extensively to regulation of metabolic processes and they will not be discussed further here. On the other hand a role for negative metabolic regulators should also be considered. In the liver, expression of the genes for fatty acid metabolism enzymes appears to be regulated by two reciprocal systems involving PPAR and LXR nuclear hormone receptors. Competitive dimerization of these proteins with RXRs (influenced by binding of their respective ligands) balances the cellular choice between lipogenesis and ß-oxidation (22). Possibly the isomers of PPAR (PPAR
) and LXR (LXRß), present at appreciable levels in the mammary gland (see Table S1) mediate these reactions (23).
Effect of Changing Dietary Fat Content on Gene Expression
Because fatty acid synthesis in the mammary gland is known to be regulated by dietary lipid (35), we sought to confirm a role for some of the regulatory molecules described above. We fed FVB mice diets containing 8 or 40% of their caloric content as lipid, substituting lipid for carbohydrate, from L5 to L10. Gene expression was examined in the mammary gland using the same Affymetrix gene chips. We found that expression of 85 genes was significantly altered by the high fat diet in the mammary gland, compared with 760 in the liver (Table S6), reflecting the much more complex role of the liver in metabolic regulation. By comparison the lactating mammary gland is devoted to one process, the synthesis and secretion of milk. Although the mammary gland is known to decrease fatty acid synthesis in response to a high fat diet (35), there are few other diet-sensitive changes. Because of space limitation we confine this discussion to the metabolic and regulatory genes previously analyzed in this paper.
In the liver four genes for fatty acid degradation were upregulated, but expression in the mammary gland remained low (Fig. 6). Expression of many hepatic cholesterol synthesis genes in the liver was also affected by the increased dietary lipid although the direction of the changes was inconsistent. In the mammary gland the major effects were to increase expression of the mRNA for the major fatty acid transporter FAT 3, and decrease the mRNA for enzymes of the pentose phosphate shunt, the mitochondrial citrate transporter, and fatty acid synthesis. Of the regulatory genes, expression of C/EBPß and C/EBP
as well as, most significantly, SREBP-1 was decreased. Of the 12 metabolic genes whose mammary expression was significantly altered, eight have been shown in numerous experiments to be regulated by SREBP-1c (18), firmly implicating this transcription factor in regulation of fatty acid synthesis in the mammary gland. Surprisingly, expression of TGF-ß3 and IGFBP5 was increased, possibly suggesting that high-fat, low-carbohydrate diets lead to a degree of apoptosis of the mammary epithelium. Linkages between energy metabolism and cell survival have been made in many systems (7, 38) and possibly warrant additional investigation in the mammary gland. C/EBPß has been found to be necessary for ductal morphogenesis, lobuloalveolar differentiation, and functional differentiation of the mammary gland (13). We found that it is significantly decreased by diets with fat content >40%, suggesting that it may respond to changes in dietary lipid.
|
-lactalbumin (
-lact) leading to highly increased lactose synthesis in the lactating gland. The high levels of lactose observed by MRS are also consistent. The data also lead to multiple hypotheses that will need to be tested in subsequent experiments examining both substrate fluxes and enzyme activity. These are:
, LXRß, and fatty acids generated by CTE-1 activity, playing inhibitory roles during pregnancy.
These studies indicate that significant regulation of lactose and lipid synthesis in the mammary gland occurs at the level of mRNA expression. From studies in other tissues we surmise that much of this regulation is at the level of gene transcription and, indeed, a number of activating transcription factors including SREBP-1c, STAT5a and STAT5b, and C/EBPß, as well as inhibitory factors, PPAR
and LXRß, are likely involved. This study of mRNA levels and prominent metabolites has allowed us to obtain a global picture of the adjustments in gene expression necessary to turn a quiescent organ into an incredibly efficient machine for the massive synthesis and secretion of a complex mixture that includes highly regulated concentrations of lactose and lipid. These data provide a comprehensive framework for testing of the hypotheses generated, testing that will have to be carried out at the protein and functional levels.
| GRANTS |
|---|
|
|
|---|
| FOOTNOTES |
|---|
Article published online before print. See web site for date of publication (http://physiolgenomics.physiology.org).
| REFERENCES |
|---|
|
|
|---|
-lactalbumin transgene in
-lactalbumin-deficient mice can rescue lactation. In vivo relationship between bovine
-lactalbumin expression content and milk composition. J Dairy Res 64: 145148, 1997.[CrossRef][Web of Science][Medline]This article has been cited by other articles:
![]() |
Q. K. Kleerekoper and J. A. Putkey PEP-19, an Intrinsically Disordered Regulator of Calmodulin Signaling J. Biol. Chem., March 20, 2009; 284(12): 7455 - 7464. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. D. Maningat, P. Sen, M. Rijnkels, A. L. Sunehag, D. L. Hadsell, M. Bray, and M. W. Haymond Gene expression in the human mammary epithelium during lactation: the milk fat globule transcriptome Physiol Genomics, March 3, 2009; 37(1): 12 - 22. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. J. Loor and W. S. Cohick ASAS Centennial Paper: Lactation biology for the twenty-first century J Anim Sci, February 1, 2009; 87(2): 813 - 824. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Bobrovnikova-Marjon, G. Hatzivassiliou, C. Grigoriadou, M. Romero, D. R. Cavener, C. B. Thompson, and J. A. Diehl PERK-dependent regulation of lipogenesis during mouse mammary gland development and adipocyte differentiation PNAS, October 21, 2008; 105(42): 16314 - 16319. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. R. Andrechek, S. Mori, R. E. Rempel, J. T. Chang, and J. R. Nevins Patterns of cell signaling pathway activation that characterize mammary development Development, July 15, 2008; 135(14): 2403 - 2413. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Bionaz and J. J. Loor ACSL1, AGPAT6, FABP3, LPIN1, and SLC27A6 Are the Most Abundant Isoforms in Bovine Mammary Tissue and Their Expression Is Affected by Stage of Lactation J. Nutr., June 1, 2008; 138(6): 1019 - 1024. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. E. Bauman, J. W. Perfield II, K. J. Harvatine, and L. H. Baumgard Regulation of Fat Synthesis by Conjugated Linoleic Acid: Lactation and the Ruminant Model J. Nutr., February 1, 2008; 138(2): 403 - 409. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. D. Russell, C. A. Palmer, D. J. Orlicky, E. S. Bales, B. H.-J. Chang, L. Chan, and J. L. McManaman Mammary glands of adipophilin-null mice produce an amino-terminally truncated form of adipophilin that mediates milk lipid droplet formation and secretion J. Lipid Res., January 1, 2008; 49(1): 206 - 216. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. N. VanHouten, M. C. Neville, and J. J. Wysolmerski The Calcium-Sensing Receptor Regulates Plasma Membrane Calcium Adenosine Triphosphatase Isoform 2 Activity in Mammary Epithelial Cells: A Mechanism for Calcium-Regulated Calcium Transport into Milk Endocrinology, December 1, 2007; 148(12): 5943 - 5954. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. D. Russell, C. A. Palmer, D. J. Orlicky, A. Fischer, M. C. Rudolph, M. C. Neville, and J. L. McManaman Cytoplasmic lipid droplet accumulation in developing mammary epithelial cells: roles of adipophilin and lipid metabolism J. Lipid Res., July 1, 2007; 48(7): 1463 - 1475. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Visit Other APS Journals Online |