|
|
||||||||
Mammalian NutriPhysioGenomics, Department of Animal Sciences and Division of Nutritional Sciences, University of Illinois, Urbana 61801
1 Corresponding author: jloor{at}illinois.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: transcriptomics nuclear receptor transcription regulator
| INTRODUCTION |
|---|
|
|
|---|
Studies in nonlactating sheep and lactating cows showed that short-term (hours to a few days) manipulation of nutrients [e.g., propionate or conjugated linoleic acid (CLA)] available to the animal can have substantial effects on subcutaneous adipose tissue (SUBQ) gene expression (Lee and Hossner, 2002; Harvatine et al., 2009). Based on those data, it is likely that long-term nutritional regulation of milk fat synthesis during lactation involves not only transcriptomic adaptations in mammary tissue but also in adipose tissue. Milk production, fat composition and yield, and SUBQ rates of lipogenesis and esterification are affected by stage of lactation and supplemental lipid (McNamara et al., 1995).
Metabolic regulation in complex organisms relies partly on transcriptional control as a long-term mechanism affecting the level of expression of several key enzymes (Desvergne et al., 2006). A cellular gene network can be defined as a collection of DNA segments that interact with a regulator such as a transcription factor or nuclear receptor, but also with each other through their RNA and protein products and with other molecules in the cell (Wittkopp, 2007). These "global" interactions, thus, can govern the rates at which genes in the network are transcribed into mRNA. In rodents, there is high correlation between mRNA expression of target genes and recruitment of lipogenic transcription factors or nuclear receptors and their co-regulatory proteins to promoter regions, suggesting that gene expression analysis is useful for inferring transcriptional activity (Bennett et al., 2008). Transcriptional regulation of gene networks involved in hepatic lipogenesis and adipogenesis in adipose tissue of rodents is under control of sterol regulatory element binding factor 1 (SREBF1) (Horton et al., 2002) and the ligand-activated nuclear receptor PPAR
(PPARG), respectively (Fernyhough et al., 2007).
As a means to evaluate long-term metabolic regulation, our group has recently studied transcriptional networks involved in mammary (Bionaz and Loor, 2008a,b) or skeletal muscle tissue (Graugnard et al., 2009) during lactation or rapid growth. In the course of similar studies of long-term nutritional regulation of milk fat synthesis (Thering, 2008), we sought to examine transcriptional changes of adipogenic and lipogenic networks in SUBQ. Our hypothesis was that long-term diet-induced milk fat depression in response to a blend of fish and soybean oil would be associated with greater mRNA abundance of tail-head SUBQ lipogenic or adipogenic gene networks (Figure 1). Specific objectives were to measure mRNA abundance of 20 genes encoding proteins required for fatty acid uptake and activation, intracellular fatty acid transport, de novo fatty acid synthesis, esterification, desaturation, transcriptional regulation of adipogenesis and differentiation, and lipid droplet formation.
|
| MATERIALS AND METHODS |
|---|
|
|
|---|
Diets were based on corn silage and alfalfa silage (Thering, 2008). Treatments included a control diet with no added lipid (n = 5 cows), a milk fat-depressing diet containing a blend of fish oil (10 g/kg; Omega Protein Inc., Houston, TX) and soybean oil (25 g/kg; Archer Daniels Midland, Decatur, IL; n = 5 cows, FSO), and a diet containing a saturated lipid supplement (Energy Booster 100, Milk Specialties, Dundee, IL; 35 g/kg of DM; n = 6 cows, EB100). Diets had a forage:concentrate ratio of 70:30 (DM basis), with lipid supplements in FSO and EB100 replacing soybean hulls. Ingredient and chemical composition of experimental diets were reported previously (Thering, 2008). Intake of NEL was calculated as DMI (kg/d) x NEL (Mcal/kg of DM) content of each diet (1.55, 1.66, and 1.67, respectively, for control, FSO, and EB100). Milk NEL output (Mcal/d) was calculated as milk yield (kg) x [0.0929 x (fat %) + 0.0563 x (true protein %) + 0.0395 x (lactose %)] (NRC, 2001). Diets were fed ad libitum as a TMR.
Individual BW were recorded weekly; individual feed intake was recorded daily. Samples of fresh TMR, feed ingredients, and feed refusals were collected weekly and stored at –20°C. At the end of the experiment, feed samples were oven-dried (110°C) for 24 h and DM percentage determined.
Milk yields were electronically recorded at each milking. Samples of milk for the determination of fat, protein, and lactose were collected from each cow at 0500 and 1600 h each day (total of 56 individual samples). Composite samples of milk (28 total) from a.m. and p.m. milkings were obtained by combining a volume of a.m. or p.m. sample based on the proportion of milk produced at each milking. Samples were stored with Bronopol preservative (D&F Control Systems Inc., Dublin, CA) at 4°C before milk composition analysis by mid-infrared methods (AOAC, 1995) at a DHIA laboratory (Dairy One, Ithaca, NY).
Blood samples from the coccygeal artery/vein were collected before the morning feeding on d 21 of feeding the diets. Samples were analyzed to determine concentrations of NEFA, BHBA, TAG, and glucose using commercial kits in an autoanalyzer at the University of Illinois Veterinary Diagnostic Laboratory (Urbana, IL).
Biopsies
On d 21 relative to start of treatments, and following an approved Institutional Animal Care and Use Committee protocol, the tail-head area of all cows was evaluated visually and by palpation to determine if sufficient amounts of SUBQ could be harvested without the need for extensive dissection around the incision, which may have resulted in excessive trauma to the chosen area or excessive pain/distress to the animal. An additional concern was the fact that, in our experience, the yield of total RNA per gram of SUBQ is markedly lower compared with that from tissues such as liver and mammary. It was not expected that the cow selection approach would confound our analysis because the chosen cows responded to dietary treatments in the same fashion (Table 1) as the entire group (Thering, 2008). Based on the above criteria, it was deemed that only 3 cows per treatment were adequate for biopsies. Before biopsies, the hair on the tail-head and to one side of the tail-head was clipped closely and thoroughly scrubbed with surgical soap. Local anesthesia was applied over the area between the point of the ischium and coccygeal vertebra. A 6- to 8-cm incision was made and the skin pulled back using sterile forceps, exposing the SUBQ. Samples were taken using a sterile scalpel blade and forceps. After the sample was taken, pressure was applied with sterile gauze to stop any external bleeding. The incision was closed with 8 to 12 surgical staples (#89063337, Appose ULC Skin Stapler, 35 wide, Henry Schein Inc., Melville, NY) and iodine ointment was applied. Biopsied tissue (3–5 g) was weighed and stored in liquid N2 before RNA extraction.
|
Complementary DNA was synthesized using 100 ng of RNA, 1 µg of dT18 (Operon Biotechnologies, Huntsville, AL), 1 µL of 10 mmol/L dNTP mix (Invitrogen Corp., Carlsbad, CA), 1 µL of random primers (Invitrogen Corp.), and 10 µL of DNase/RNase–free water. The mixture was incubated at 65°C for 5 min and kept on ice for 3 min. A total of 6 µL of master mix composed of 4.5 µL of 5x First-Strand Buffer, 1 µL of 0.1 M dithithreitol, 0.25 µL (50 U) of SuperScript III RT (Invitrogen Corp.) and 0.25 µL of RNase Inhibitor (10 U, Promega, Madison, WI) was added. The reaction was performed in an Eppendorf Mastercycler Gradient using the following temperature program: 25°C for 5 min, 50°C for 60 min, and 70°C for 15 min. Sample cDNA was then diluted 1:4 with DNase/RNase–free water, and cDNA for the standard curve was diluted 1:3 to ensure that it encompassed all levels of expression.
Quantitative PCR was performed using 4 µL of diluted cDNA combined with 6 µL of a mixture composed of 5 µL of 1 x SYBR Green master mix (Applied Biosystems, Foster City, CA), 0.4 µL each of 10 µM forward and reverse primers, and 0.2 µL of DNase/RNase–free water in a MicroAmp Optical 384-Well Reaction Plate (Applied Biosystems). Quantitative PCR performance for all genes including internal controls is shown in Supplemental Table 1 (http://jds.fass.org/content/vol92/issue9/). Each sample was run in triplicate and a 5-point relative standard curve plus the nontemplate control were used (User Bulletin #2, Applied Biosystems). The reactions were performed in an ABI Prism 7900 HT SDS instrument (Applied Biosystems) using the following conditions: 2 min at 50 °C, 10 min at 95 °C, 40 cycles of 15 s at 95 °C (denaturation) and 1 min at 60 °C (annealing + extension). The presence of a single PCR product was verified by the dissociation protocol using incremental temperatures to 95°C for 15 s plus 65°C for 15 s. Data were calculated with the 7900 HT Sequence Detection Systems Software (version 2.2.1, Applied Biosystems). The final data were filtered to remove unreliable data and subsequently were normalized (5 to 9 replicates per diet) using the geometric mean of ubiquitously expressed transcript (UXT), eukaryotic translation initiation factor 3, subunit K (EIF3K), and mitochondrial ribosomal protein L39 (MRPL39), which were identified as suitable internal controls among several tested (Kadegowda et al., 2009a).
Details on primer design, testing, sequencing, and actual primer features were reported previously (Bionaz and Loor, 2008a,b; Graugnard et al., 2009). Briefly, primers were designed using Primer Express 2.0 with minimum amplicon size of 80 bp (when possible amplicons of 100 to 150 bp were chosen) and limited 3' G+C (Applied Biosystems). When possible, primer sets were designed to fall across exon–exon junctions. Primers were aligned against publicly available databases using BLASTN at National Center of Biotechnology Information (National Center for Biotechnology Information, 2008) and UCSCs Cow (Bos taurus) Genome Browser Gateway (Genome Browser Gateway, 2008).
Gene Network Analysis
Summary networks among the chosen genes (Figure 1) were developed using the web-based software package Ingenuity Pathway Analysis (www.ingenuity.com; Redwood City, CA). The networks were generated using the respective gene identifiers and not the actual fold-changes in expression. Connections among genes were based on known relationships available in the Ingenuity Pathway Analysis knowledge base. This is a proprietary, manually curated database based on the published literature primarily in rodents and humans.
Statistical Analysis
All data were analyzed using the MIXED procedure in SAS (SAS Institute Inc., Cary, NC). Normalized mRNA abundance data were log-transformed before statistical analysis. The model to examine treatment differences on mRNA abundance and blood serum metabolites consisted of treatment as fixed effect and cow within treatment as random effect. A similar model including treatment, time, and treatment x time as fixed effects and cow within treatment as random effect was used to determine differences in daily DMI, NEL intake, milk NEL yield, surplus NEL intake, weekly BW, milk production, milk composition, and milk component yield. Least squares means (LSM) were separated using the PDIFF statement in SAS.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
|
|
. Dietary fish oil increases the total fatty acid concentration in TAG of very low density lipoproteins (VLDL; Offer et al., 2001). However, because the saturated lipid supplement (EB100) in our study did not affect milk fat percentage or yield from controls (Table 2) it is likely that greater adipose LPL mRNA with FSO was driven by reduced use of circulating TAG for milk fat production (approximately 26% lower fat yield vs. control or EB100; Table 2, Figure 2); that is, more VLDL was available for use by SUBQ. Under that scenario, there is no reason to believe that the change in LPL mRNA was not gradual and opposite to the pattern of milk fat and energy yield (Table 2, Figure 2); that is, over the 21-d period there was more TAG available for SUBQ utilization that might, in and of itself, have served as feed-forward mechanism for upregulation of LPL. Lipoprotein lipase activity and volume of adipocytes are substantially lower in SUBQ from lactating compared with dry/fattening cows (Chilliard and Robelin, 1985). The volume of subcutaneous adipocytes was not affected greatly by dietary fat after peak lactation but it increased over the course of lactation (McNamara et al., 1995). Further, the capacity of SUBQ for lipogenesis from acetate and esterification of long-chain fatty acids also should be greater past peak lactation (McNamara et al., 1995). Thus, it is intuitive that a reduction in milk fat production that leads to a surplus of dietary energy past peak lactation would divert use of substrates for lipogenesis (acetate, BHBA) or direct esterification (VLDL-derived long-chain fatty acids) to adipose tissues.
Despite often reducing DMI, FSO increased ruminal organic matter and fiber digestibility (Wonsil et al., 1994; Doreau and Chilliard, 1997) as well as ruminal propionate concentration (Doreau and Chilliard, 1997). In addition, reduced milk fat synthesis with FSO previously led to lower mammary uptake of BHBA (approximately 60% from controls) without affecting its arterial plasma concentration (Loor et al., 2005), which also implies that this lipogenic precursor would be available to SUBQ for lipogenesis. In the present study, we did not observe changes in arterial concentrations of BHBA, NEFA, glucose, or TAG (Table 2). Although not statistically significant, the numerically greater ACSS2 and FASN (Figure 4) as well as the greater capacity of SUBQ for lipogenesis at >120 DIM (McNamara et al., 1995) suggests that de novo synthesis might have been increased in cows fed FSO. The ACSS2 gene encodes a cytosolic protein with high affinity for acetate and it channels acetate toward fatty acid synthesis in nonruminants (Luong et al., 2000). We previously observed a marked increase in mammary ACSS2 mRNA abundance at the onset and during lactation (Bionaz and Loor, 2008b), and its transcript pattern corresponded with bovine mammary acetyl-CoA production throughout lactation (Mellenberger et al., 1973).
The expression of classical lipogenic genes in SUBQ increased in response to milk fat depression induced by short-term exogenous trans-10,cis-12–18:2 (Harvatine et al., 2009). Infusion of CLA led to greater mRNA of LPL, FASN, and SCD as well as the transcription regulators SREBF1, THRSP, and PPARG (Harvatine et al., 2009). We found that some of these genes were also greater statistically (LPL, SCD) or numerically (FASN, THRSP) by feeding FSO (Figure 4). This response, along with greater (P
0.07) ADFP and LPIN1 (Figure 4), both of which would allow for TAG (Reue and Zhang, 2008) and lipid droplet formation (McManaman et al., 2007), are suggestive of greater lipid deposition. Volume of SUBQ adipocytes increased beyond 120 DIM in cows fed control or lipid-supplemented diets (McNamara et al., 1995). Further, the capacity for esterification at >120 DIM is also expected to be greater in cows fed supplemental lipid. Thus, it would appear that the total pool size of long-chain fatty acids available to SUBQ is a likely driving force for greater lipid filling. However, it remains to be determined if very long chain polyunsaturated fatty acids (PUFA) from fish oil were directly responsible for the changes in gene expression observed. Previous work from our laboratory showed that 20:5 versus 16:0 reduced markedly the mRNA of LPL, SCD, and SREBF1 in MacT cells (Kadegowda et al., 2009b), all of which would be required for TAG synthesis. Palmitate induced intracellular lipid droplet formation but 20:5 did not. Further studies with isolated adipocytes would have to be conducted to clarify the tissues response to specific fatty acids during lactation.
Contrary to a previous short-term (4 d) study of CLA-induced milk fat depression (Harvatine et al., 2009), we did not observe greater SREBF1 in adipose tissue from cows fed FSO despite lower milk fat production (Figure 4). This was not unexpected because SREBF1 is acutely regulated by insulin and activates the lipogenic transcriptional cascade primarily in rodent liver (Horton et al., 2002). Furthermore, at least in rodents, SREBF1 is a well-established PPARG target gene (Kast-Woelbern et al., 2004) and we have found evidence in vitro (Kadegowda et al., 2009b) and in vivo (Graugnard et al., 2009) of the same in cattle. Our temporal data from longissimus muscle of rapidly growing cattle clearly showed that PPARG upregulation preceded any increases in SREBF1, and PPARG responses were paralleled by a battery of genes required for TAG synthesis (Graugnard et al., 2009).
There is no evidence, to our knowledge, that long-chain fatty acids directly bind to or promote maturation of SREBF1 to elicit a change in its mRNA expression; that is, most of the negative effects of fatty acids (primarily by PUFA) on SREBF1 mRNA abundance or transcription activation are via increases in decay or decreases in stability of the mRNA (Postic et al., 2007). Induction of SREBF1 in liver activates expression of ATP citrate lyase, glucose-6-phosphate dehydrogenase, ACACA, FASN, SCD, and glycerol-3-phophate acyltransferase, mitochondrial (GPAM) leading to synthesis of palmitic acid, oleic acid, and formation of TAG. Although SREBF1 activation is associated with greater lipogenic gene expression in rodent liver and adipose tissue (Horton et al., 2003; Li et al., 2003), recent data (Sekiya et al., 2007) found that lipogenic gene expression in adipose was independent of SREBF1. Thus, it is doubtful that SREBF1 is essential for regulation of adipocyte lipogenic/adipogenic gene expression in SUBQ of lactating cows as we have proposed recently in rapidly growing cattle (Graugnard et al., 2009).
It is noteworthy that cows fed FSO also had numerically greater insulin-induced gene 1 (INSIG1) mRNA expression (Figure 4). In murine adipose tissue INSIG1 mRNA is induced by activation of PPAR
but it occurs relatively late in the adipogenic program, preceded by peak of PPARG and SREBF1 mRNA expression (Kast-Woelbern et al., 2004). Interestingly, the primary function of INSIG1 in adipose tissue or liver is to block processing of SREBF1 (Li et al., 2003), and upregulation of INSIG1 transcription in adipocytes effectively downregulated expression of PPARG and SREBF1. The net result of INSIG1 action is the control of preadipocyte differentiation (Li et al., 2003), thus providing a feedback signaling mechanism to restrict both lipogenesis and adipogenesis. Although limited by number of animals, our data do not seem to support such a mechanism in lactating cow SUBQ because INSIG1 mRNA expression was numerically greater in cows fed FSO, which also had greater mRNA of LPL, SCD, ADFP, and LPIN1. Greater INSIG1 expression together with its apparent regulation via PPAR
suggest that the protein encoded by this gene might be pro-lipogenic rather than anti-lipogenic in cattle (Graugnard et al., 2009) as inferred in the mouse (Kast-Woelbern et al., 2004).
Expression of FABP3 mRNA was lower in cows fed the saturated lipid diet (EB100) compared with control or FSO (Figure 4). Fatty acid binding protein and acyl-CoA binding protein (ACBP or DBI) are the main intracellular fatty acid transporters in nonruminant cells (McArthur et al., 1999). The former has high affinity for free long-chain fatty acids but can also bind acyl-CoA (Whetstone et al., 1986; Frolov et al., 1997). We previously observed the presence of mRNA of all FABP isoforms except FABP2 mRNA in bovine mammary tissue, with greater abundance and upregulation of FABP3 mRNA during lactation (Bionaz and Loor, 2008a). In a subsequent study, we showed that FABP3 was the second most abundant transcript (~16%) in mammary tissue among 45 measured, in accordance with the large cytosolic content of its protein in mammary epithelium (Whetstone et al., 1986). However, nonruminant data (Rosen and MacDougald, 2006) have revealed that FABP4 (another key PPARG target gene) is the most abundant and biologically relevant FABP in adipose tissue and our previous work with cattle longissimus muscle tissue provides evidence of the same in cattle. Thus, the biological role of lower FABP3 due to EB100 in the present study is not readily apparent. Cows with milk fat depression during a short-term trans-10,cis-12–18:2 abomasal infusion had greater FABP4 in SUBQ (Harvatine et al., 2009), suggesting that along with lipogenic genes (FASN, THRSP, LPL) it helps coordinate adipogenesis. The lower FABP3 mRNA due to EB100 could be related to the fact that these cows were diverting dietary fatty acids toward the mammary gland for milk fat synthesis (Table 2, Figure 2).
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
|
|
|---|
Received for publication December 27, 2008. Accepted for publication May 26, 2009.
| REFERENCES |
|---|
|
|
|---|
Related articles in JDS:
activation and long-chain fatty acids alter lipogenic gene networks in bovine mammary epithelial cells to various extents
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |