Alterations in lipid metabolism occur when animals are exposed to different feeding systems. In the last few decades, the characterisation of genes involved in fat metabolism and technological advances have enabled the study of the effect of diet on the milk fatty acid (FA) profile in the mammary gland and aided in the elucidation of the mechanisms of the response to diet. The aim of this study was to evaluate the effect of different forage diets (grazing vs. hay) near the time of ewe parturition on the relationship between the fatty acid profile and gene expression in the mammary gland of the Churra Tensina sheep breed.
In this study, the forage type affected the C18:2 cis-9 trans-11 (CLA) and long-chain saturated fatty acid (LCFA) content, with higher percentages during grazing than during hay feeding. This may suggest that these FAs act as regulatory factors for the transcriptional control of the carnitine palmitoyltransferase 1B (CPT1B) gene, which was more highly expressed in the grazing group (GRE). The most highly expressed gene in the mammary gland at the fifth week of lactation is CAAT/ enhancer- binding protein beta (CEBPB), possibly due to its role in milk fat synthesis in the mammary gland. More stable housekeeping genes in the ovine mammary gland that would be appropriate for use in gene expression studies were ribosomal protein L19 (RPL19) and glyceraldehyde- 3- phosphate dehydrogenase (GAPDH).
Small changes in diet, such as the forage preservation (grazing vs. hay), can affect the milk fatty acid profile and the expression of the CPT1B gene, which is associated with the oxidation of fatty acids. When compared to hay fed indoors, grazing fresh low mountain pastures stimulates the milk content of CLA and LCFA via mammary uptake. In this sense, LCFA in milk may be acting as a regulatory factor for transcriptional control of the CPT1B gene, which was more highly expressed in the grazing group.
Sheep milk fat is rich in medium-chain triglycerides (MCT), which are composed of fatty acids with a carbon chain of 6–10 carbon atoms. MCTs are of special therapeutic interest because of their specific metabolism and their consequent application in certain types of metabolic illnesses . By contrast, medium-chain saturated fatty acids, mainly C12:0 and C14:0, are considered to have a negative effect on human health when consumed in excess . Other fatty acids (FAs), such as oleic (C18:1 n-9), linolenic (C18:3 n-3), eicosapentaenoic acid (EPA), docosahexaenoic acid (DHA) and conjugated linoleic acid (CLA, C18:2 cis-9 trans-11), have positive effects on human health [3, 4]. Milk fat composition is affected by factors such as: nutrition, lactation stage and genetics. With diet having a high and immediate impact on milk FA composition in cows and goats [5, 6], as well as in sheep [7–10].
In the last few decades, the characterisation of genes involved in fat metabolism and technological advances have enabled the study of the effect of diet on the milk FA profile within the mammary gland and aided in the elucidation of the mechanisms of the response to diet [11, 12]. Studies that have examined the effects of diet on gene expression in sheep and goats under the same dietary treatment are limited [13, 14]. In most studies, the feedstuffs used in animal diets have differed to relate the effects of diet to gene expression or enzymatic activity. Accordingly, grazing versus concentrate-rich diets [15, 16] or the use of high-, medium- and low-quality diets  have been compared. However, few studies have explored the role of forage preservation on milk or meat FA profiles and the expression of genes related to fat metabolism. The hay making process results in the loss of fatty acids which act as precursors for CLA synthesis in the rumen and mammary gland . Total FA was reduced by over 50% during hay making with a greater loss of linolenic acid (C18:3 n-3) . To further address these issues, we studied the effect of forage preservation (grazing vs. hay) fed to ewes on the milk fatty acid profile and the expression of 8 key genes implicated in fat metabolism within the mammary gland: lipoprotein lipase (LPL), acetyl-CoA carboxylase α (ACACA), fatty acid synthase (FASN), fatty acid binding protein 3 (FABP3) and 4 (FABP4), diglyceride acyltransferase 1 (DGAT1), stearoyl-CoA desaturase (SCD) and carnitine palmitoyltransferase 1B (CPT1B). Lipoprotein lipase (LPL) is a rate-limiting enzyme in the hydrolysis of triglycerides circulating in the form of chylomicrons and very low-density lipoproteins into free fatty acids (FFAs) and 2-monoacylglycerols. The resulting FFAs can be utilised in different tissues, such as the mammary gland. The enzyme acetyl-coenzyme A carboxylase α (ACC), which is encoded by the acetyl-CoA carboxylase α gene (ACACA), catalyses the ATP-dependent carboxylation of acetyl-CoA to form malonyl-CoA, which is the substrate for the synthesis of palmitic acid and long-chain fatty acids (acyl-CoA > C22:0) by the enzyme fatty acid synthase (FAS). In the presence of NADPH, FAS, which is encoded by the FASN gene, catalyses the synthesis of long-chain saturated fatty acids from palmitate derived from acetyl-CoA and malonyl-CoA. Fatty acid binding proteins 3 and 4, which are encoded by the FABP3 and FABP4 genes, respectively, supply long-chain fatty acids as an important energy source and have been assumed to be mammary gland- and adipose tissue-specific, respectively. Carnitine palmitoyltransferase I (M- CPT 1), which is encoded by the CPT1B gene, is part of the mitochondrial transport system and is a key enzyme in the control of long-chain fatty acid oxidation. Stearoyl-CoA desaturase (SCD) is the rate-limiting enzyme that converts palmitoyl- and stearoyl-CoA to palmitoleoyl- and oleoyl-CoA, respectively. In ruminants, SCD is encoded by the SCD gene and participates in the formation of conjugated fatty acid (CLA) from C18:1 trans-11 in animal tissues. The DGAT1 gene encodes the enzyme acyl-CoA:diacylglycerol acyltransferase1, which plays a central role in the synthesis of triglycerides. It catalyses the reaction of diacylglycerol and fatty acyl-CoA to form triglycerides. Most of these key enzymes are subject to both acute and chronic control, and their gene expression is regulated by important transcription factors such as sterol regulatory element binding protein (SREBP1)peroxisome proliferator – activated receptor gamma (PPARG) and alpha (PPARA), and CAAT/ enhancer- binding protein beta (CEBPB), also studied in this work. Furthermore, the stability of housekeeping genes within the mammary gland for use in gene expression studies was assessed.
Experimental site and sheep breed
The experiment was conducted at La Garcipollera Research Station, in the Pyrenees (north-eastern Spain, 4237 N, 030 W, 945 m above sea level, a.s.l.), during autumn 2008. The average annual rainfall is bimodally distributed with peaks in spring and autumn, with dry summers and some precipitation in the form of snow in autumn and winter. The mean temperature during September, October, November and December was 16.2, 11.6, 5.0 and 3.0°C, while the precipitation in these months was 31.8, 120.6, 94.0 and 59.2 mm, respectively.
The breed used in this study was the Churra Tensina sheep, an endangered local coarse-wooled hardy breed (approximately 7,000 head) belonging to the Churra group that is raised for lamb production in the mountain area of the southern Pyrenees. The main Churra product is meat lamb, which can be of two commercial categories, suckling (10–12 kg LW) or light lamb (22–24 kg LW).
Animal management and experimental design
Twenty-four multiparous single-bearing Churra Tensina ewes were used (body weight (BW) of 46.5 ± 0.66 kg and a body condition score (BCS) of 3.0 ± 0.03, on a 1 to 5 scale , at breeding). The flock was permanently grazing on medium-ranged mountain pastures (920–1500 m a.s.l.) during spring (68-72% neutral detergent fibre, NDF; 10-14% crude protein, CP). In summer (early to mid-pregnancy), the flock grazed on high mountain pastures (1500–2200 m a.s.l.) (68-75% NDF, 7-10% CP). The mating period was in spring from 5 May to 15 June, and the lambing was registered, on average, on 19 October ± 8 days. Four Churra Tensina rams were used for mating with the ewes.
To evaluate the effect of the diet (grazing vs. hay) on the relationship between the fatty acid profile and gene expression within the mammary gland of the Churra Tensina sheep breed, the experimental study began when pregnant ewes were in the last third of pregnancy (last 5 weeks) and finished at 5 weeks of lactation, when lambs reached the target BW of 10–12 kg for inclusion in the suckling lamb commercial category. Therefore, 5 weeks prior to parturition, ewes were divided into two groups:
GRE treatment: 12 ewes were fed on mountain pastures ad libitum (43–46% NDF and 18–20% CP). During the pre-partum period, snow limited pasture availability for 11 days (15th-19th and 22 th -27th October), and ewes were offered 1 kg/day/ewe of alfalfa (Medicago sativa) pellets (35–37% NDF, 14–17% CP).
HAY treatment: 12 ewes were fed pasture hay ad libitum indoors (49–52% NDF, 11–13% CP). This hay was made in late spring from the same pasture paddocks used by those grazed in autumn.
The pasture was composed of 22% legumes (mainly Trifolium repens), 68% grass (the main species were Festuca arundinaceaFestuca pratensis and Dactylis glomerata) and 10% other species (mainly Rumex acetosa and Ranunculus bulbosus) . Outdoor grazing ewes had free access to a sheltered area.
The BW of ewes and lambs was recorded at weekly intervals at 8 a.m. with an electronic balance (0.1-kg precision). The BCS was measured according to Russel et al.  by two trained technicians at lambing and 5 weeks post-partum (drying off).
After lambing, milk production was recorded weekly by the oxytocin technique proposed by Doney et al.  with machine milking with hand finishing (4 h interval). Ewes were injected with 5 IU oxytocin in the jugular vein prior to machine milking and hand finishing at 08:00 and 12:00 h. Ewes were returned to their paddock between the two milkings, while lambs were confined. The milk obtained in the second milking was weighed, and the yield was extrapolated to the daily period (daily production = milk obtained x 6). Individual milk samples (50 ml) were preserved by adding potassium dichromate and stored at −20°C until future milk fatty acid analysis during the subsequent 3 months. The sample that was taken on the last day of lactation was used to study the relationship between milk FA composition and gene expression in the mammary gland.
Mammary gland biopsy
At the fifth week of lactation, mammary gland tissue was collected with a biopsy needle (Tru-Cut®, CareFusion, France). The animal handling and procedures were in accordance with current European legislation (directive 86/609/EEC)  and supervised by the Animal Welfare Committee of the institution (protocol number 2009-01_MJT). Ewes were positioned on a lateral recumbence. The biopsy site of the mammary gland was prepared by cleaning with 10% Povidone-iodine solution (Lainco, S.A. Barcelona, Spain), and local anaesthesia was administered by subcutaneous injection of 2.5 ml lidocaine chlorhydrate (Xilocaina ovejero 2%, Ovejero Laboratorios, León, Spain). A 2-mm incision was made approximately 2 cm above the nipple to facilitate the insertion of the biopsy needle. The disposable needle had a specimen notch and a sheath over the needle to cut the tissue, which was necessary to withdraw the needle and sheath with the specimen from the mammary gland. After the biopsy, the mammary gland surface was sprayed with Aludemin (Chemical Ibérica, Salamanca, Spain) to protect the wound from dirt. Immediately after biopsy, the mammary tissue sample was stored in RNA later solution (QIAGEN, Madrid, Spain) and then frozen and stored at −80°C.
Feed and milk fatty acid analysis
For the FA determination, feed samples were Soxhlet-extracted , and milk samples were submitted to a lipid separation as described by Luna et al. . Both FA groups were converted to methyl esters by base-catalysed methanolysis of the glycerides with KOH according to the UNE-EN ISO 5509:2000 methods . Fatty acid methyl esters were separated on a capillary column (HP-88 100 m x 0.25 mm ID and 0.20 μm film thickness, Agilent Technologies, Waldbronn, Germany) and detected with a flame ionisation detector (FID, gas chromatograph HP-6890). The carrier gas was helium, and the flow rate was 2 ml/min. The temperatures of the inlet and the detector were maintained at 250 and 300°C, respectively. The temperature program was as follows: the initial temperature was held at 100°C after injection, then programmed to increase at 1.5°C/min to 170°C (held there for 15 min), then programmed to increase at 0.5°C/min to 180°C (held there for 2 min), and finally to increase at 10°C/min to 215°C (held there for 20 min). The injection volume was 1.0 μl. Fatty acids were identified based on comparison with the retention times of a standard FA mixture (Sigma-Aldrich, Madrid, Spain). The individual fatty acid contents were expressed as weight percentages (g/100 g of total FA).
After determining individual FAs, the sum of saturated fatty acids (SFA), short-chain FA (SCFA; C4-C10), medium-chain FA (MCFA; C12-C14), long-chain FA (LCFA; C16-C24), total unsaturated FA (UFA), mono-unsaturated FA (MUFA) and poly-unsaturated FA (PUFA) were calculated, according to De la Fuente et al. . Moreover, the UFA/SFA and n-6/n-3 ratios were determined.
RNA extraction and cDNA synthesis
Total RNA was prepared from the mammary gland biopsy with TRI@REAGENT (Sigma-Aldrich, Madrid, Spain) according to the manufacturer's instructions. The concentration and quality of the RNA were determined by nanophotometric analysis (Implen, Madrid, Spain). To exclude possible amplification of contaminating genomic DNA, samples were treated with DNAse. Single-stranded cDNA was synthesised from 1 μg RNA with the SuperScript®III Reverse Transcriptase kit (Invitrogen, Carlsbad, CA, USA) following the manufacturer’s recommendations.
Gene expression levels were determined by RT-PCR on an ABI Prism 7500 platform (Applied Biosystems, Madrid, Spain). Primers for RT-PCR were designed with Primer3 software (http://frodo.wi.mit.edu/primer3/). The sequences of the primers and the RT-PCR conditions are described in Dervishi et al.. In this study, the FABP3 gene was also studied. The primer sequences were as follow: 5’- TTCAAGCTGGGAGTCGAGTT-3’ and 5’- TGTCCATTCCACTTCTGCAC-3’
RT-PCR was performed under the same conditions as the other genes with an annealing temperature, primer concentrations, correlation (R2) and slope of 60°C, [900, 900], 0.99 and −3.35, respectively. Before performing the RT-PCR reactions, a conventional PCR was performed for all genes to test the primers and verify the amplified products. The PCR products were sequenced to confirm gene identity with an ABI Prism3700 (Applied Biosystem, Madrid, Spain), and standard protocols. Homology searches were performed with BLAST (National Center for Biotechnology Information: http://www.ncbi.nlm.nih.gov/BLAST/) to confirm the identity of the amplified fragments. The PCR reaction was performed in a 10-μl PCR total reaction mixture containing SYBR Green PCR Master Mix (Applied Biosystem, Madrid, Spain). Each reaction was performed in triplicate, and the average was used to calculate the relative amount of the target gene. To normalise the results of the target genes, 8 candidate housekeeping genes were tested. In this study, the housekeeping genes analysed included those for ovine: tyrosine 3-monooxygenase (YWHAZ), ribosomal protein L19(RPL19), glyceraldehyde-3-phosphate dehydrogenase (GAPDH), glucose-6-phosphate dehydrogenase (G6PDH), beta actin (ACTB), ubiquitin C (UBC), beta-2-microglobulin (B2M), and succinate dehydrogenase (SDHA). The sequences of the primers and the RT-PCR conditions are described in Dervishi et al.. Standard curves for genes were generated to calculate the amplification efficiency. The efficiency of PCR amplification for each gene was calculated with the standard curve method (E = 10-1/slope -1). The standard curves for each gene were generated by 5-fold serial dilution of pooled cDNA. The amplification conditions were an initial step of 10 min at 95°C, followed by 40 cycles of 95°C for 15 sec and 59 or 60°C for 30 sec. The specificity of the amplification products and the lack of primer dimers were confirmed by melting curve analysis in all cases. To quantify the relative gene expression, the standard curve method was used according to the recommendation of Larionov et al. . Normalised RT-PCR data were transformed to the fold-change relative to the control group, GRE; thus, data in the GRE group were transformed to obtain a perfect mean of 1. PCR-normalised data are presented as the n-fold relative change.
Statistical analyses were performed with the SPSS statistical software package, version 19.0. To evaluate the effects of the forage type fed to the ewe on milk fatty acid content and gene expression in the mammary gland, the analyses were performed with a General Linear Model (GLM), in which the diet of the ewe was included as a fixed effect and the ewe body weight at lambing was included as a covariable, due to its possible relevance on the studied parameters during lactation . The following equation was used for the model: Yijk = μ + Aj + b(Ek) + (A*b(E))jk + eijk, where Yijk = dependent variable; μ = overall mean; Aj = the effect of type of forage fed to the ewe; b(Ek) = body weight at lambing; b = linear regression coefficient; (A*b(E))jk = interaction between the forage type and body weight at lambing; eijk = residual error. Fatty acid results were expressed as least square means ± the standard error (SE) values, and the relative differences in gene expression between the groups were calculated and defined as the relative increase, setting the control means at 100%. To test the stability of housekeeping genes, we used geNorm and NormFinder.
Results and discussion
During the study period (70 days), pasture biomass decreased steadily from 14.4 ± 4.2 to 5.1 ± 0.6 cm. In the final pre-partum period, there was snow for 11 days during which 1 kg of alfalfa pellets was fed to each grazing ewe per day. Permanent pasture or alfalfa pellets were the dietary components of the grazing treatment (GRE), whereas the diet of the hay treatment (HAY) was pasture hay made in late spring from the same pasture paddocks used for autumn grazing (Table 1). The daily voluntary hay intake during the last 5 weeks of pregnancy was, on average, 1.38 ± 0.31 kg/ewe (as-fed basis), and this increased to 1.55 ± 0.31 kg/ewe (as-fed basis) during the 5 weeks of lactation.
Chemical composition (% of dry matter, DM) and fatty acid content (in % of total fatty acid methyl esters, FAME) of pasture hay, fresh pasture grazed and alfalfa pellets fed to Churra Tensina ewes during the peripartum period1
fresh pasture grazed
Dry matter M (DM, % of fresh matter)
Crude protein, CP
Neutral detergent fibre, NDF
Acid detergent fibre, ADF
C18:1 trans-11 (VA)
C18:2 cis-9 cis-12
1SFAs (saturated fatty acids) were calculated as the sum of C12:0 + C14:0 + C16:0 + C17:0 + C18:0 + C20:0; MCFAs (medium-chain fatty acids) were calculated as the sum of C12:0 + C14:0; LCFAs (long-chain fatty acids) were calculated as the sum of C16:0 + C18:0 + C20:0; UFAs (unsaturated fatty acids) were calculated as the sum of (MUFA + PUFA); MUFAs (monounsaturated fatty acids) were calculated as the sum of (C16:1 + C17:1 + C18:1n-9 + C18:11n-7 + C20:1); PUFAs (polyunsaturated fatty acids) were calculated as the sum of (C18:2 t-t + C18:2c-c + C18:3n-6 + C18:3n-3 + CLA + C20:2 + C20:3 + C20:4 + C20:5 + C22:5 + C22:6).
2Supplement during 11 days in the pre-partum grazing ewes, when snow limited pasture availability.
The dietary chemical composition and fatty acid profile are shown in Table 1. The estimated dietary metabolisable energy (ME) content according to Cannas et al.  was 9.1 and 10.2 MJ ME/kg DM for pasture and hay, respectively. In addition, the alfalfa pellets supplied 10.2 MJ ME/kg DM to grazing ewes when snow prevented pasture grazing. The higher CP and lower fibre (NDF and acid detergent fibre, ADF) contents in fresh pasture compared to hay indicate that grazed herbage had an early stage of growth and maintained vegetative characteristics during the experiment. Both feeding systems offered rather similar estimated energy contents, and crude protein was higher in GRE than in HAY.
The ewes BCS at lambing was 2.69 and 2.61 ± 0.07 in HAY and GRE ewes, respectively, while at drying-off, it was 2.32 and 2.45 ± 0.08, respectively. Although the forages supplied (fresh or hay) met the theoretical energy and protein requirements for milk production , there was BCS loss between lambing and drying-off in both treatments, suggesting that the voluntary intake of ewes did not meet the requirements for avoiding a negative energy balance.
The feedstuffs used in the GRE treatment (fresh pasture and alfalfa pellets) presented a 31.4% greater overall proportion of C18:3 n-3 than HAY feed. Alfalfa pellets exhibited a greater proportion of C18:2 cis-9 cis-12 and a lower proportion of C18:1 trans-11 (VA) than the other forages, but it is expected to exert only a minor influence due to the limited amount of intake. Hay feed had a 78.4% greater proportion of C18:3 n-6 than the other forages. By contrast, fresh pasture had a 36.05 and 16.66% greater proportion of C18:3 n-3 than hay feed and alfalfa pellets. These results are in agreement with Biondi et al. , who also observed that C18:3 n-3 was predominant in fresh pasture whereas C18:3 n-6 was more important in hay feed.
The sums of LCFA, total PUFA and PUFA n-3 were 11.2%, 35.6% and 37.6% greater in the feedstuffs used in the GRE treatment (fresh pasture and alfalfa pellets) than in the HAY treatment. Forage FA composition depends on the species and variety of the plant, the phenological phase of the plant and the preservation of the plant [34, 35].
Milk fatty acid profile
The fatty acid composition of milk at the fifth week of lactation is shown in Table 2. There were significant effects of forage type on C17:1 cis-10, C18:0, C18:2 cis-9 cis-12, CLA and C20:2 cis-11 cis-14 (P < 0.01) contents. The remaining milk FA was not significantly affected by forage type. Ewes belonging to the HAY group exhibited a higher content of C17:1 cis-10 and C18:2 cis-9 cis-12 than the GRE group. By contrast, the C18:0, CLA and C20:2 contents were higher in the GRE group (P < 0.01). The animal diets in both treatments were based exclusively on forage and the only difference was forage preservation, linked to plant phenological phase. Cabiddu et al.  also observed higher CLA and C18:3 n-3 contents in milk when ewes grazed pasture in the vegetative phase, whereas the CLA content decreased considerably when the same pasture was in the reproductive phase. In the present study, the most remarkable difference among forages was the greater C18:3 n-3 content in the fresh pasture compared to hay feed, in agreement with the loss of FA precursors of CLA caused by the hay-making process . The lack of effect of diet on most of the milk FA composition contradicts reports in the literature [7, 8, 19, 37], although most of these studies compared the effect of forage vs. concentrate diets.
Fatty acid content (in %) in milk of Churra Tensina ewes at the fifth week of lactation, according to forage type1–2
C18:2 trans-9 trans-12
C18:2 cis-9 cis-12
CLA C18:2 cis-9 trans-11
C20:2 cis-11 cis-14
C20:3 cis-8 cis-11 cis-14
1Within a row, means without a common superscript differ (P < 0.05). SE is the standard error.
2HAY group: ewes received hay from mountain pastures; GRE group: ewes grazed low mountain pastures.
Milk from grazing ewes contained higher levels of CLA than those fed hay indoors. Similarly, higher contents of CLA and vaccenic acid (VA, C18:1 trans-11) in milk and dairy products were observed when animals grazed fresh pastures compared to dried forages . Grazing pasture may enhance the growth of specific bacteria in the rumen, stimulating the production of CLA and/or blocking the final reduction of VA to stearic acid (C18:0) , as well as increasing PUFA n-3 [5, 34]. In the present study, grazing increased milk CLA content. In this sense, fresh pasture and alfalfa are richer in C18:3n-3 (21.89 and 18.90 vs. 14.00 for GRE and alfalfa vs. HAY, respectively), which can stimulate the CLA production.
Milk FAs were grouped according to the length and degree of saturation from a health perspective or according to ratios related to the human health (Table 3). Forage preservation affected only the LCFA and PUFA n-6 groups (P < 0.05). GRE ewes presented a significantly greater proportion of LCFA in milk than the HAY group, which agrees with their higher proportion in grazed pasture as well in dehydrated alfalfa compared to hay feed (Table 1). Total PUFA n-6 milk fatty acids were greater in the HAY group than in the GRE group due to the greater amount of linoleic acid (C18:2 cis-9 cis-12) in milk from HAY ewes than in milk from GRE ones. Collectively, the present study demonstrated that even small changes in diet, such as the forage preservation form, can affect the milk fatty acid profile.
Main FA categories content (%) in milk of Churra Tensina ewes at the fifth week of lactation, according to forage type1–3
1Within a row, means without a common superscript differ (P < 0.05). SE is the standard error.
2HAY group: ewes received hay from mountain pastures; GRE group: ewes grazed low mountain pastures.
3SFAs (saturated fatty acids) were calculated as the sum of C4:0 + C6:0 + C8:0 + C10:0 + C12:0 + C14:0 + C16:0 + C17:0 + C18:0 + C20:0 + C22:0 + C24:0; SCFAs (short-chain fatty acids) were calculated as the sum of C4:0 + C6:0 + C8:0 + C10:0; MCFAs (medium-chain fatty acids) were calculates as the sum of C12:0 + C14:0; LCFAs (long-chain fatty acids) were calculated as the sum of C16:0 + C18:0 + C20:0 + C22:0 + C24:0; UFAs (unsaturated fatty acids) were calculated as the sum of MUFA + PUFA; MUFAs (monounsaturated fatty acids) were calculated as the sum of C16:1 + C17:1 + C18:1n-9 + C18:11n-7 + C20:1; PUFAs (polyunsaturated fatty acids) were calculated as the sum of C18:2 t-t + C18:2c-c + C18:3n-6 + C18:3n-3 + CLA + C20:2 + C20:3 + C20:4 + C20:5 + C22:5 + C22:6.
Selection of optimal housekeeping genes
To normalise the mammary gland gene expression results, we examined the stability of eight housekeeping genes with two different programs, geNorm and NormFinder. More stable genes analysed with the geNorm program were GAPDH (M = 0.60), G6PDH (M = 0.651) and RPL19 (M = 0.825). When analysed with NormFinder, the RPL19 (stability value = 0.11) and GAPDH (stability value = 0.11) genes still exhibited higher levels of gene expression stability. To normalise the data, we used the geometric mean of the RPL19 and GAPDH genes. Both programs determined the most stable housekeeping genes, and we decided to use the most stable genes determined with the NormFinder program because this program considers the type of treatment, which was one of the most important factors.
Real-time polymerase chain reaction analysis (RT-PCR) and the effect of feeding system on gene expression
In the present study, we analysed the expression of 12 genes related to lipid metabolism within the mammary gland. We found significant differences only in the expression of the CPT1B gene (Figure 1). We did not observe differences in the expression of the rest of genes, which may be due to the similarity of the diets. The CPT1B gene, which is associated with the oxidation of fatty acids, is expressed 0.34-fold lower in HAY ewes than in GRE ewes. The enzyme carnitine palmitoyltransferase (CPT1) is located in the outer mitochondrial membrane and represents the main site of control for the entry of LCFA into the mitochondria . CPT1 catalyses the synthesis of long-chain acylcarnitines from long-chain CoA esters, and its activity has a high degree of control over the rate of fatty acid oxidation. Furthermore, in the liver, the utilisation of fatty acids for ketogenesis is also controlled by CPT1. CPT1 regulation of oxidation would suggest that the treatments altered energy balance more on the GRE ewes, which may be driven by exercise or lower intake. Furthermore, the promoter of CPT1B could be regulated by LCFA levels, and it has been previously reported that LCFAs are important regulatory factors for transcriptional control of CPT1 in the rat neonatal liver  and in rat primary cardiac myocytes. Furthermore, a response element regulated by LCFA levels has been reported for both the rat and human CPT1B promoters [41–43]. In the present work, the forage type only affected the LCFA content of milk, demonstrating that LCFA mammary gland uptake was stimulated in grazing compared to hay-fed ewes. In this sense, LCFA may be acting as a regulatory factor for transcriptional control of the CPT1B gene. As this gene lacks TATA boxes, the baseline transcription regulation relies on transcription factor specificity protein-1 (Sp1) . The induction of the CPT1B gene occurs via the fatty acid response element (FARE) or peroxisome proliferator-activated receptor (PPAR) response elements (PPREs), which bind the transcription factor PPAR [41, 45]. However, in our study, we did not observe significant differences in the expression of either the PPARA or PPARG gene in the mammary gland at the fifth week of lactation. Therefore, LCFA may regulate CPT1B gene expression via another mechanism. M-CPT1 is subject to regulation at the transcriptional level and to acute control by malonyl-CoA . The N-terminal domain of the M-CPT1 enzyme is essential for malonyl-CoA inhibition . Rasmussen and Winder  suggested that muscle contraction inhibits acetyl-CoA carboxylase α (ACC) and would be one mechanism for increased fatty acid oxidation. A decrease in malonyl-CoA (during starvation) promotes fatty acid oxidation because M-CPT1 becomes uninhibited, whereas an increase in malonyl-CoA levels (during a fed state) will decrease mitochondrial fatty acid uptake and oxidation . In the HAY group, ewes were fed pasture hay indoors ad libitum, while GRE group ewes were grazing fresh pasture throughout the experiment. In this sense, the treatments altered energy balance more on the GRE ewes, which may be driven by exercise. However, intake differences and other compounds such as protein and the plant secondary metabolites could also influence the results. Future studies are being planned to verify these factors not studied in the present work.
These results are in complete agreement with those obtained previously by Dervishi et al. , in which grazing systems promoted higher levels of CPT1B gene expression in the semitendinous muscle and fatty acid oxidation for energy production. These results suggest that grazing promotes higher levels of CPT1B gene expression in the mammary gland at the fifth week of lactation.
The results also revealed that CEBPB is the most highly expressed gene among all of the genes studied in the mammary gland at the fifth week of lactation (data not shown). Shi et al.  have also suggested a role for this factor in milk-fat synthesis. The activity of this factor may be a key regulator of the nutritional control of ACACA expression in lipogenic tissues. SREBF1 was also highly expressed, and a central role for this gene as a mediator of FA effects has been outlined .
In conclusion, low mountain pastures increased the milk content of C18:2 cis-9 trans-11 (CLA) and long chain fatty acids (LCFA) via mammary uptake, relative to hay fed indoors. CPT1B gene was more highly expressed in the grazing group possibly due to an altered energy balance by reducing energy intake and increasing exercise. Furthermore, LCFA in milk may be also acting as a regulatory factor for the transcriptional control of the CPT1B gene. More stable housekeeping genes in the ovine mammary gland for gene expression studies were RPL19 and GAPDH. It was demonstrated that small changes in diet, such as forage preservation, can cause differences in the milk fatty acid profile and in the expression of CPT1B.
The authors wish to thank the staff of ‘La Garcipollera’ and of ‘CITA de Aragón’ for their assistance in sample collection and analysis. Special thanks to I. Blasco, R. Ripoll-Bosch and G. Estopañan for their assistance and suggestions. E. Dervishi was supported by a doctoral grant from AECI (Agencia Española de Cooperación Internacional), and J. Alvarez-Rodriguez from INIA-DGA. This study was partially financed by the research projects INIA-RTA-2008-0098, INIA-RTA-2009-0091-C02 and INIA-RZP-2009-005 and the Research Group Funds of the Aragón Government (A49 and A11).
Unidad de Tecnología en Producción Animal, CITA
Departament de Producció Animal, Universitat de Lleida
Chilliard YA, Ferlay R, Mansbridge M, Doreau M: Ruminant milk fat plasticity: nutritional control of saturated, polyunsaturated, trans and conjugated fatty acids.Ann Zootech 2000, 49:181–205.View Article
Williams CM: Dietary fatty acids and human health.Ann Zootech 2000, 49:165–180.View Article
Massaro M, Carluccio MA, de Caterina R: Direct vascular antiatherogenic effects of oleic acid: a clue to the cardioprotective effects of the Mediterranean diet.Cardiologia 1999, 44:507–513.PubMed
Pariza MW, Park Y, Cook ME: The biologically active isomers of conjugated linoleic acid.Prog Lipid Res 2001, 40:283–298.PubMedView Article
Chilliard Y, Ferlay A, Doreau M: Effect of different types of forages, animal fat or marine oils in cow's diet on milk fat secretion and composition, especially conjugated linoleic acid (CLA) and polyunsaturated fatty acids.Livest Prod Sci 2001, 70:31–48.View Article
Chilliard Y, Ferlay A: Dietary lipids and forage interactions on cow and goat milk fatty acid composition and sensory properties.Reprod Nut Devel 2004, 44:467–492.View Article
Scerra M, Caparra P, Foti F, Galofaro V, Sinatra MC, Scerra V: Influence of ewe feeding systems on fatty acid composition of suckling lambs.Meat Sci 2007, 76:390–394.PubMedView Article
Joy M, Alvarez-Rodriguez J, Revilla R, Delfa R, Ripoll G: Ewe metabolic performance and lamb carcass traits in pasture and concentrate-based production systems in Churra Tensina breed.Small Rumin Res 2008, 75:24–35.View Article
Joy M, Ripoll R, Molino F, Dervishi E, Alvarez-Rodriguez J: Influence of the type of forage supplied to ewes in pre- and post-partum periods on the meat fatty acids of suckling lambs.Meat Sci 2012, 90:775–782.PubMedView Article
Tsiplakou E, Mountzouris KC, Zervas G: Concentration of conjugated linoleic acid in grazing sheep and goat milk fat.Livest Sci 2006, 103:74–84.View Article
Bernard L, Leroux C, Faulconnier Y, Durand D, Shingfield KJ, Chilliard Y: Effect of sunflower-seed oil or linseed oil on milk fatty acid secretion and lipogenic gene expression in goats fed hay-based diets.J Dairy Res 2009, 76:241–248.PubMedView Article
Bionaz M, Loor JJ: Gene networks driving bovine milk fat synthesis during the lactation cycle.BMC Genomics 2008, 9:366.PubMedView Article
Tsiplakou E, Flemetakis E, Kalloniati C, Papadomichelakis G, Katinakis P, Zervas G: Sheep and goats differences in CLA and fatty acids milk fat content in relation with mRNA stearoyl-CoA desaturase and lipogenic genes expression in their mammary gland.J Dairy Res 2009,76(4):392–401.PubMedView Article
Tsiplakou E, Flemetakis E, Kalloniati C, Zervas G: Differences in mRNA lipogenic gene expression in the subcutaneous adipose tissue of sheep and goats under the same dietary treatments.Small Rumin Res 2011,99(2–3):110–115.View Article
Tsiplakou E, Kominakis A, Zervas G: The interaction between breed and diet on CLA and fatty acids content of milk fat of four sheep breeds kept indoors or at grass.Small Rumin Res 2008,74(1–3):179–187.View Article
Dervishi E, Serrano C, Joy M, Serrano M, Rodellar C, Calvo JH: The effect of feeding system in the expression of genes related with fat metabolism in semitendinous muscle in sheep.Meat Sci 2011, 89:91–97.PubMedView Article
Reverter A, Byrne KA, Bruce HL, Wang YH, Dalrymple BP, Lehnert SA: A mixture model- based cluster analysis of DNA microarray gene expression data on Brahman and Brahman composite steers fed high-, medium-, and lowquality diets.J Anim Sci 2003, 81:1900–1910.PubMed
Doreau M, Poncet C: Ruminal biohydrogenation of fatty acids originating from fresh or preserved grass.Reprod Nutr Dev 2000, 40:201.
Valvo MA, Lanza M, Bella M, Fasone V, Scerra M, Biondi L: Effect of ewe feeding system (grass v. concentrate) on intramuscular fatty acids of lambs raised exclusively on maternal milk.Anim Sci 2005, 81:431–436.View Article
Russel AJF, Doney JM, Gunn RG: Subjective assessment of body fat in live sheep.J Agr Sci 1969, 72:451–454.View Article
Alvarez-Rodriguez J, Sanz A, Delfa R, Revilla R, Joy M: Performance and grazing behaviour of Churra Tensina sheep stocked under different management systems during lactation on Spanish mountain pastures.Livest Sci 2007, 107:152–161.View Article
Doney JM, Peart JN, Smith WF, Louda F: Consideration of the techniques for estimation of milk-yield by suckled sheep and a comparison of estimates obtained by 2 methods in relation to the effect of breed, level of production and stage of lactation.J Agr Sci 1979, 92:123–132.View Article
European Union (EU) Directive No. 86/609/CEE: Council Directive of 24 November 1986 on the approximation of laws, regulations and administrative provisions of the Member States regarding the protection of animals used for experimental and other scientific purposes.Offic J Europ Commun 1986,1986(358):1–32.
Sukhija PS, Palmquist DL: Rapid method for determination of total fatty acid content and composition of feedstuffs and faeces.J Agric Food Chem 1988, 36:1202–1206.View Article
Luna P, Juárez M, de la Fuente MA: Validation of rapid milk fat separation method to determine the fatty acid profile by gas cromatography.J Anim Sci 2005, 88:3377–3381.
BOE: Aceites y grasas de origen animal y vegetal. Preparación de ésteres metílicos de ácidos grasos.Resolución de 7 de mayo de 2001, de la Dirección General de Política Tecnológica, por la que se publica la relación de normas UNE aprobadas por AENOR durante el mes de febrero de 2001 2001, 136:20223–20227.
De La Fuente LF, Barbosa E, Carriedo JA, Gonzalo C, Arenas R, Fresno JM, San Primitivo F: Factors influencing variation of fatty acid content in ovine milk.J Dairy Sci 2009, 92:3791–3799.PubMedView Article
Larinov A, Krause A, Miller W: A standard curve based method for relative real time PCR data processing.BMC Bioinformatics 2005, 6:62.View Article
Andersen CL, Jensen JL, Orntoft TF: Normalization of Real- Time quantitative reverse transcription- PCR data: A model- based variance estimation approach to identify genes suited for normalization, applied to bladder and colon cancer data set.Cancer Res 2004, 64:5245–5250.PubMedView Article
Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes.Genome Biol 2002,3(7):0034.View Article
Cannas A, Tedeschi LO, Fox DG, Pell AN, Van Soest PJ: A mechanistic model for predicting the nutrient requirements and feed biological values for sheep.J Anim Sci 2004, 82:149–169.PubMed
INRA: Alimentation des bovins, ovins et caprins. Besoin des animaux, valeurs des aliments. Quae Editions edition. Versailles, France; 2007.
Biondi L, Valvo MA, Di Gloria M, Tenghi ES, Galofaro V, Priolo A: Changes in ewe milk fatty acids following turning out to pasture.Small Rumin Res 2008, 75:17–23.View Article
Dewhurst RJ, Shingfield KJ, Lee MRF, Scollan ND: Increasing the concentrations of beneficial polyunsaturated fatty acids in milk produced by dairy cows in high-forage systems.Anim Feed Sci Tech 2006,131(3–4):168–206.View Article
Woods VB, Fearon AM: Dietary sources of unsaturated fatty acids for animals and their transfer into meat, milk and eggs: A review.Livest Sci 2009,126(1–3):1–20.View Article
Cabiddu A, Decandia M, Addis M, Piredda G, Pirisi A, Molle G: Managing Mediterranean pastures in order to enhance the level of beneficial fatty acids in sheep milk.Small Rumin Res 2005,59(2–3):169–180.View Article
Atti N, Rouissi H, Othmane MH: Milk production, milk fatty acid composition and conjugated linoleic acid (CLA) content in dairy ewes raised on feedlot or grazing pasture.Livest Sci 2006, 104:121–127.View Article
Nudda A, McGuire MA, Battacone G, Pulina G: Seasonal variation inconjugated linoleic acid and vaccenic acid in milk fat of sheep and its transfer to cheese and ricotta.J Dairy Sci 2005, 88:1311–1319.PubMedView Article
McGarry JD, Woeltje KF, Kuwajima M, Foster DW: Regulation of ketogenesis and the renaissance of carnitine palmitoyltransferase.Diabetes Metab Rev 1989, 5:271–284.PubMedView Article
Yu GS, Lu YC, Gulick T: Co-regulation of tissue-specific alternative human carnitine palmitoyltransferase Iβ gene promoters by fatty acid enzyme substrate.J Biol Chem 1998, 273:32901–32909.PubMedView Article
Brandt JM, Djouadi F, Kelly DP: Fatty acids activate transcription of muscle carnitine palmitoyltransferase 1 gene in cardiac myocytes via peroxisome proliferator activated receptor alpha.J Biol Chem 1998, 273:23786–23792.PubMedView Article
Moore ML, Wang GL, Belaguli NS, Schwartz RJ, McMillin JB: GATA-4 and serum response factor regulate transcription of the muscle-specific carnitine palmitoyltransferase I beta in rat heart.J Biol Chem 2001, 276:1026–1033.PubMedView Article
Steffen ML, Harrison WR, Elder FF, Cook GA, Park EA: Expression of the rat liver carnitine palmitoyltransferase I (CPT-Ialpha) gene is regulated by Sp1 and nuclear factor Y: chromosomal localization and promoter characterization.Biochem J 1999, 340:425–432.PubMedView Article
Djouadi F, Brandt JM, Weinheimer CJ, Leone TC, Gonzalez FJ, Kelly DP: The role of the peroxisome proliferator-activated receptor alpha (PPAR alpha) in the control of cardiac lipid metabolism.Prostaglandins Leukotrienes Essent Fatty Acids 1999, 60:339–343–443–460.View Article
McGarry JD, Leatherman GF, Foster DW: Carnitine palmitoyltransferase I: the site of inhibition of hepatic fatty acid oxidation by malonyl-CoA.J Biol Chem 1978, 253:4128–4136.PubMed
Rasmussen BB, Winder WW: Effect of exercise intensity on skeletal muscle malonyl- CoA and acetyl- CoA carboxylase.Am Physiol Soc 1997, 83:1104–1109.
Bonnefont JP, Djouadi F, Prip-Buus C, Gobin S, Munnich A, Bastin J: Carnitine palmitoyltransferases 1 and 2: biochemical, molecular and medical aspects.Mol Aspects Med 2004, 25:495–520.PubMedView Article
Shi X, Liu S, Metges CC, Seyfert HM: C/EBP-beta drives expression of the nutritionally regulated promoter IA of the acetyl-CoA carboxylase-alpha gene in cattle.Biochim Biophys Acta 2010, 1799:561–567.PubMedView Article
Harvatine KJ, Bauman DE: SREBP1 and Thyroid Hormone Responsive Spot 14 (S14) are involved in the regulation of bovine mammary lipid synthesis during diet-induced milk fat depression and treatment with CLA.J Nutr Biochem 2006,136(10):2468–2473.