BACKGROUND
sheep milk fat is rich in medium-chain triglycerides , which are composed of fatty acids with a carbon chain of 6– <dig> carbon atoms.
mcts are of special therapeutic interest because of their specific metabolism and their consequent application in certain types of metabolic illnesses  <cit> .
by contrast, medium-chain saturated fatty acids, mainly c12: <dig> and c14: <dig>  are considered to have a negative effect on human health when consumed in excess  <cit> .
other fatty acids , such as oleic , linolenic , eicosapentaenoic acid , docosahexaenoic acid  and conjugated linoleic acid , have positive effects on human health  <cit> .
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  <cit> , as well as in sheep  <cit> .
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  <cit> .
studies that have examined the effects of diet on gene expression in sheep and goats under the same dietary treatment are limited  <cit> .
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  <cit>  or the use of high-, medium- and low-quality diets  <cit>  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  <cit> .
total fa was reduced by over 50% during hay making with a greater loss of linolenic acid   <cit> .
to further address these issues, we studied the effect of forage preservation  fed to ewes on the milk fatty acid profile and the expression of  <dig> key genes implicated in fat metabolism within the mammary gland: lipoprotein lipase , acetyl-coa carboxylase α , fatty acid synthase , fatty acid binding protein  <dig>  and  <dig> , diglyceride acyltransferase  <dig> , stearoyl-coa desaturase  and carnitine palmitoyltransferase 1b  <cit> .
lipoprotein lipase  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  and 2-monoacylglycerols.
the resulting ffas can be utilised in different tissues, such as the mammary gland.
the enzyme acetyl-coenzyme a carboxylase α , which is encoded by the acetyl-coa carboxylase α gene , 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  by the enzyme fatty acid synthase .
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  <dig> and  <dig>  which are encoded by the fabp <dig> and fabp <dig> 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 , 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  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  from c18: <dig> trans- <dig> in animal tissues.
the dgat <dig> gene encodes the enzyme acyl-coa:diacylglycerol acyltransferase <dig>  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 peroxisome proliferator – activated receptor gamma  and alpha , and caat/ enhancer- binding protein beta , also studied in this work.
furthermore, the stability of housekeeping genes within the mammary gland for use in gene expression studies was assessed.
methods
experimental site and sheep breed
the experiment was conducted at la garcipollera research station, in the pyrenees , during autumn  <dig>  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  <dig> ,  <dig> ,  <dig>  and  <dig> °c, while the precipitation in these months was  <dig> ,  <dig> ,  <dig>  and  <dig>  mm, respectively.
the breed used in this study was the churra tensina sheep, an endangered local coarse-wooled hardy breed  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  or light lamb .
animal management and experimental design
twenty-four multiparous single-bearing churra tensina ewes were used  of  <dig>  ±  <dig>  kg and a body condition score  of  <dig>  ±  <dig> , on a  <dig> to  <dig> scale  <cit> , at breeding).
the flock was permanently grazing on medium-ranged mountain pastures  during spring .
in summer , the flock grazed on high mountain pastures  .
the mating period was in spring from  <dig> may to  <dig> june, and the lambing was registered, on average, on  <dig> october ± 8 days.
four churra tensina rams were used for mating with the ewes.
to evaluate the effect of the diet  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  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:  <dig> ewes were fed on mountain pastures ad libitum .
during the pre-partum period, snow limited pasture availability for 11 days , and ewes were offered 1 kg/day/ewe of alfalfa  pellets .
· hay treatment:  <dig> ewes were fed pasture hay ad libitum indoors .
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 , 68% grass  and 10% other species   <cit> .
outdoor grazing ewes had free access to a sheltered area.
measurements
the bw of ewes and lambs was recorded at weekly intervals at 8 a.m. with an electronic balance .
the bcs was measured according to russel et al.
<cit>  by two trained technicians at lambing and 5 weeks post-partum .
after lambing, milk production was recorded weekly by the oxytocin technique proposed by doney et al.
<cit>  with machine milking with hand finishing .
ewes were injected with 5 iu oxytocin in the jugular vein prior to machine milking and hand finishing at 08: <dig> 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 .
individual milk samples  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 .
the animal handling and procedures were in accordance with current european legislation   <cit>  and supervised by the animal welfare committee of the institution .
ewes were positioned on a lateral recumbence.
the biopsy site of the mammary gland was prepared by cleaning with 10% povidone-iodine solution , and local anaesthesia was administered by subcutaneous injection of  <dig>  ml lidocaine chlorhydrate .
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  to protect the wound from dirt.
immediately after biopsy, the mammary tissue sample was stored in rna later solution  and then frozen and stored at −80°c.
feed and milk fatty acid analysis
for the fa determination, feed samples were soxhlet-extracted  <cit> , and milk samples were submitted to a lipid separation as described by luna et al.
<cit> .
both fa groups were converted to methyl esters by base-catalysed methanolysis of the glycerides with koh according to the une-en iso 5509: <dig> methods  <cit> .
fatty acid methyl esters were separated on a capillary column  and detected with a flame ionisation detector .
the carrier gas was helium, and the flow rate was 2 ml/min.
the temperatures of the inlet and the detector were maintained at  <dig> 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  <dig> °c/min to 170°c , then programmed to increase at  <dig> °c/min to 180°c , and finally to increase at 10°c/min to 215°c .
the injection volume was  <dig>  μl.
fatty acids were identified based on comparison with the retention times of a standard fa mixture .
the individual fatty acid contents were expressed as weight percentages .
after determining individual fas, the sum of saturated fatty acids , short-chain fa , medium-chain fa , long-chain fa , total unsaturated fa , mono-unsaturated fa  and poly-unsaturated fa  were calculated, according to de la fuente et al.
<cit> .
moreover, the ufa/sfa and n-6/n- <dig> ratios were determined.
rna extraction and cdna synthesis
total rna was prepared from the mammary gland biopsy with tri@reagent  according to the manufacturer's instructions.
the concentration and quality of the rna were determined by nanophotometric analysis .
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  following the manufacturer’s recommendations.
real-time polymerase chain reaction analysis
gene expression levels were determined by rt-pcr on an abi prism  <dig> platform .
primers for rt-pcr were designed with primer <dig> software .
the sequences of the primers and the rt-pcr conditions are described in dervishi et al.
<cit> .
in this study, the fabp <dig> 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  and slope of 60°c,  <cit> ,  <dig>  and − <dig> , 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 prism <dig> , and standard protocols.
homology searches were performed with 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 .
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,  <dig> candidate housekeeping genes were tested.
in this study, the housekeeping genes analysed included those for ovine: tyrosine 3-monooxygenase , ribosomal protein l <dig>  glyceraldehyde-3-phosphate dehydrogenase , glucose-6-phosphate dehydrogenase , beta actin , ubiquitin c , beta-2-microglobulin , and succinate dehydrogenase .
the sequences of the primers and the rt-pcr conditions are described in dervishi et al.
<cit> .
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 .
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  <dig> cycles of 95°c for 15 sec and  <dig> 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.
<cit> .
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  <dig>  pcr-normalised data are presented as the n-fold relative change.
statistical analyses
statistical analyses were performed with the spss statistical software package, version  <dig> .
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 , 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  <cit> .
the following equation was used for the model: yijk = μ + aj + b + )jk + eijk, where yijk = dependent variable; μ = overall mean; aj = the effect of type of forage fed to the ewe; b = body weight at lambing; b = linear regression coefficient; )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  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 <cit>  and normfinder <cit> .
RESULTS
dietary composition
during the study period , pasture biomass decreased steadily from  <dig>  ±  <dig>  to  <dig>  ±  <dig>  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 , whereas the diet of the hay treatment  was pasture hay made in late spring from the same pasture paddocks used for autumn grazing .
the daily voluntary hay intake during the last 5 weeks of pregnancy was, on average,  <dig>  ±  <dig>  kg/ewe , and this increased to  <dig>  ±  <dig>  kg/ewe  during the 5 weeks of lactation.
1sfas  were calculated as the sum of c12:0 + c14:0 + c16:0 + c17:0 + c18:0 + c20:0; mcfas  were calculated as the sum of c12:0 + c14:0; lcfas  were calculated as the sum of c16:0 + c18:0 + c20:0; ufas  were calculated as the sum of ; mufas  were calculated as the sum of ; pufas  were calculated as the sum of .
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  <dig>  the estimated dietary metabolisable energy  content according to cannas et al.
<cit>  was  <dig>  and  <dig>  mj me/kg dm for pasture and hay, respectively.
in addition, the alfalfa pellets supplied  <dig>  mj me/kg dm to grazing ewes when snow prevented pasture grazing.
the higher cp and lower fibre  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  <dig>  and  <dig>  ±  <dig>  in hay and gre ewes, respectively, while at drying-off, it was  <dig>  and  <dig>  ±  <dig> , respectively.
although the forages supplied  met the theoretical energy and protein requirements for milk production  <cit> , 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  presented a  <dig> % greater overall proportion of c18: <dig> n- <dig> than hay feed.
alfalfa pellets exhibited a greater proportion of c18: <dig> cis- <dig> cis- <dig> and a lower proportion of c18: <dig> trans- <dig>  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  <dig> % greater proportion of c18: <dig> n- <dig> than the other forages.
by contrast, fresh pasture had a  <dig>  and  <dig> % greater proportion of c18: <dig> n- <dig> than hay feed and alfalfa pellets.
these results are in agreement with biondi et al.
<cit> , who also observed that c18: <dig> n- <dig> was predominant in fresh pasture whereas c18: <dig> n- <dig> was more important in hay feed.
the sums of lcfa, total pufa and pufa n- <dig> were  <dig> %,  <dig> % and  <dig> % greater in the feedstuffs used in the gre treatment  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  <cit> .
milk fatty acid profile
the fatty acid composition of milk at the fifth week of lactation is shown in table  <dig>  there were significant effects of forage type on c17: <dig> cis- <dig>  c18: <dig>  c18: <dig> cis- <dig> cis- <dig>  cla and c20: <dig> cis- <dig> cis- <dig>  contents.
the remaining milk fa was not significantly affected by forage type.
ewes belonging to the hay group exhibited a higher content of c17: <dig> cis- <dig> and c18: <dig> cis- <dig> cis- <dig> than the gre group.
by contrast, the c18: <dig>  cla and c20: <dig> contents were higher in the gre group .
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.
<cit>  also observed higher cla and c18: <dig> n- <dig> 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: <dig> n- <dig> 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  <cit> .
the lack of effect of diet on most of the milk fa composition contradicts reports in the literature  <cit> , although most of these studies compared the effect of forage vs. concentrate diets.
1within a row, means without a common superscript differ .
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  in milk and dairy products were observed when animals grazed fresh pastures compared to dried forages  <cit> .
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   <cit> , as well as increasing pufa n- <dig>  <cit> .
in the present study, grazing increased milk cla content.
in this sense, fresh pasture and alfalfa are richer in c18:3n- <dig> , 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 .
forage preservation affected only the lcfa and pufa n- <dig> groups .
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 .
total pufa n- <dig> milk fatty acids were greater in the hay group than in the gre group due to the greater amount of linoleic acid  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.
1within a row, means without a common superscript differ .
se is the standard error.
2hay group: ewes received hay from mountain pastures; gre group: ewes grazed low mountain pastures.
3sfas  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  were calculated as the sum of c4:0 + c6:0 + c8:0 + c10:0; mcfas  were calculates as the sum of c12:0 + c14:0; lcfas  were calculated as the sum of c16:0 + c18:0 + c20:0 + c22:0 + c24:0; ufas  were calculated as the sum of mufa + pufa; mufas  were calculated as the sum of c16:1 + c17:1 + c18:1n-9 + c18:11n-7 + c20:1; pufas  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: <dig>
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 , g6pdh  and rpl <dig> .
when analysed with normfinder, the rpl <dig>  and gapdh  genes still exhibited higher levels of gene expression stability.
to normalise the data, we used the geometric mean of the rpl <dig> 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  and the effect of feeding system on gene expression
in the present study, we analysed the expression of  <dig> genes related to lipid metabolism within the mammary gland.
we found significant differences only in the expression of the cpt1b gene .
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  <dig> -fold lower in hay ewes than in gre ewes.
the enzyme carnitine palmitoyltransferase  is located in the outer mitochondrial membrane and represents the main site of control for the entry of lcfa into the mitochondria  <cit> .
cpt <dig> 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 cpt <dig>  cpt <dig> 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 cpt <dig> in the rat neonatal liver  <cit>  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  <cit> .
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- <dig>   <cit> .
the induction of the cpt1b gene occurs via the fatty acid response element  or peroxisome proliferator-activated receptor  response elements , which bind the transcription factor ppar  <cit> .
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-cpt <dig> is subject to regulation at the transcriptional level and to acute control by malonyl-coa  <cit> .
the n-terminal domain of the m-cpt <dig> enzyme is essential for malonyl-coa inhibition  <cit> .
rasmussen and winder  <cit>  suggested that muscle contraction inhibits acetyl-coa carboxylase α  and would be one mechanism for increased fatty acid oxidation.
a decrease in malonyl-coa  promotes fatty acid oxidation because m-cpt <dig> becomes uninhibited, whereas an increase in malonyl-coa levels  will decrease mitochondrial fatty acid uptake and oxidation  <cit> .
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.
<cit> , 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 .
shi et al.
<cit>  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.
srebf <dig> was also highly expressed, and a central role for this gene as a mediator of fa effects has been outlined  <cit> .
CONCLUSIONS
in conclusion, low mountain pastures increased the milk content of c18: <dig> cis- <dig> trans- <dig>  and long chain fatty acids  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 rpl <dig> 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.
authors’ contributions
ed conducted the research and the statistical analysis.
ed, mj and jhc wrote the manuscript.
as and jar contributed to the statistical analysis and the manuscript revision.
fm performed the feed and milk fatty acid analysis.
mj and jhc designed the research.
mj provided animals.
jhc has primary responsibility for the final content.
all of the authors contributed to the manuscript discussion.
all of the authors read and approved the final manuscript.

 sheep milk fat is rich in medium-chain triglycerides, which are composed of fatty acids with a carbon chain of 6– 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: and c14: are considered to have a negative effect on human health when consumed in excess. other fatty acids, such as oleic, linolenic, eicosapentaenoic acid, docosahexaenoic acid and conjugated linoleic acid, have positive effects on human health. 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, as well as in sheep. 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. studies that have examined the effects of diet on gene expression in sheep and goats under the same dietary treatment are limited. 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 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. to further address these issues, we studied the effect of forage preservation fed to ewes on the milk fatty acid profile and the expression of key genes implicated in fat metabolism within the mammary gland: lipoprotein lipase, acetyl-coa carboxylase α, fatty acid synthase, fatty acid binding protein and, diglyceride acyltransferase, stearoyl-coa desaturase and carnitine palmitoyltransferase 1b. lipoprotein lipase 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 and 2-monoacylglycerols. the resulting ffas can be utilised in different tissues, such as the mammary gland. the enzyme acetyl-coenzyme a carboxylase α, which is encoded by the acetyl-coa carboxylase α gene, 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 by the enzyme fatty acid synthase. 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 and which are encoded by the fabp and fabp 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, 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 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 from c18: trans- in animal tissues. the dgat gene encodes the enzyme acyl-coa:diacylglycerol acyltransferase 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 peroxisome proliferator – activated receptor gamma and alpha, and caat/ enhancer- binding protein beta, also studied in this work. furthermore, the stability of housekeeping genes within the mammary gland for use in gene expression studies was assessed. methods experimental site and sheep breed the experiment was conducted at la garcipollera research station, in the pyrenees, during autumn 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,, and °c, while the precipitation in these months was,, and  mm, respectively. the breed used in this study was the churra tensina sheep, an endangered local coarse-wooled hardy breed 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 or light lamb. animal management and experimental design twenty-four multiparous single-bearing churra tensina ewes were used of  ±   kg and a body condition score of  ± , on a to scale, at breeding). the flock was permanently grazing on medium-ranged mountain pastures during spring. in summer, the flock grazed on high mountain pastures. the mating period was in spring from may to june, and the lambing was registered, on average, on october ± 8 days. four churra tensina rams were used for mating with the ewes. to evaluate the effect of the diet 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 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: ewes were fed on mountain pastures ad libitum. during the pre-partum period, snow limited pasture availability for 11 days, and ewes were offered 1 kg/day/ewe of alfalfa pellets. · hay treatment: ewes were fed pasture hay ad libitum indoors. 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, 68% grass and 10% other species. outdoor grazing ewes had free access to a sheltered area. measurements the bw of ewes and lambs was recorded at weekly intervals at 8 a.m. with an electronic balance. the bcs was measured according to russel et al. by two trained technicians at lambing and 5 weeks post-partum. after lambing, milk production was recorded weekly by the oxytocin technique proposed by doney et al. with machine milking with hand finishing. ewes were injected with 5 iu oxytocin in the jugular vein prior to machine milking and hand finishing at 08: 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. individual milk samples 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. the animal handling and procedures were in accordance with current european legislation and supervised by the animal welfare committee of the institution. ewes were positioned on a lateral recumbence. the biopsy site of the mammary gland was prepared by cleaning with 10% povidone-iodine solution, and local anaesthesia was administered by subcutaneous injection of  ml lidocaine chlorhydrate. 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 to protect the wound from dirt. immediately after biopsy, the mammary tissue sample was stored in rna later solution 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: methods. fatty acid methyl esters were separated on a capillary column and detected with a flame ionisation detector. the carrier gas was helium, and the flow rate was 2 ml/min. the temperatures of the inlet and the detector were maintained at 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 °c/min to 170°c, then programmed to increase at °c/min to 180°c, and finally to increase at 10°c/min to 215°c. the injection volume was  μl. fatty acids were identified based on comparison with the retention times of a standard fa mixture. the individual fatty acid contents were expressed as weight percentages. after determining individual fas, the sum of saturated fatty acids, short-chain fa, medium-chain fa, long-chain fa, total unsaturated fa, mono-unsaturated fa and poly-unsaturated fa were calculated, according to de la fuente et al.. moreover, the ufa/sfa and n-6/n- ratios were determined. rna extraction and cdna synthesis total rna was prepared from the mammary gland biopsy with tri@reagent according to the manufacturer's instructions. the concentration and quality of the rna were determined by nanophotometric analysis. 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 following the manufacturer’s recommendations. real-time polymerase chain reaction analysis gene expression levels were determined by rt-pcr on an abi prism platform. primers for rt-pcr were designed with primer software. the sequences of the primers and the rt-pcr conditions are described in dervishi et al.. in this study, the fabp 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 and slope of 60°c,, and −, 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 prism, and standard protocols. homology searches were performed with 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. 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, candidate housekeeping genes were tested. in this study, the housekeeping genes analysed included those for ovine: tyrosine 3-monooxygenase, ribosomal protein l glyceraldehyde-3-phosphate dehydrogenase, glucose-6-phosphate dehydrogenase, beta actin, ubiquitin c, beta-2-microglobulin, and succinate dehydrogenase. 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. 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 cycles of 95°c for 15 sec and 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 pcr-normalised data are presented as the n-fold relative change. statistical analyses statistical analyses were performed with the spss statistical software package, version. 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, 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 + )jk + eijk, where yijk = dependent variable; μ = overall mean; aj = the effect of type of forage fed to the ewe; b = body weight at lambing; b = linear regression coefficient; )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 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.  dietary composition during the study period, pasture biomass decreased steadily from  ±  to  ±   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, whereas the diet of the hay treatment was pasture hay made in late spring from the same pasture paddocks used for autumn grazing. the daily voluntary hay intake during the last 5 weeks of pregnancy was, on average,  ±   kg/ewe, and this increased to  ±   kg/ewe during the 5 weeks of lactation. 1sfas were calculated as the sum of c12:0 + c14:0 + c16:0 + c17:0 + c18:0 + c20:0; mcfas were calculated as the sum of c12:0 + c14:0; lcfas were calculated as the sum of c16:0 + c18:0 + c20:0; ufas were calculated as the sum of; mufas were calculated as the sum of; pufas were calculated as the sum of. 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  the estimated dietary metabolisable energy content according to cannas et al. was and  mj me/kg dm for pasture and hay, respectively. in addition, the alfalfa pellets supplied  mj me/kg dm to grazing ewes when snow prevented pasture grazing. the higher cp and lower fibre 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 and  ±  in hay and gre ewes, respectively, while at drying-off, it was and  ± , respectively. although the forages supplied 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 presented a % greater overall proportion of c18: n- than hay feed. alfalfa pellets exhibited a greater proportion of c18: cis- cis- and a lower proportion of c18: trans- 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 % greater proportion of c18: n- than the other forages. by contrast, fresh pasture had a and % greater proportion of c18: n- than hay feed and alfalfa pellets. these results are in agreement with biondi et al., who also observed that c18: n- was predominant in fresh pasture whereas c18: n- was more important in hay feed. the sums of lcfa, total pufa and pufa n- were %, % and % greater in the feedstuffs used in the gre treatment 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. milk fatty acid profile the fatty acid composition of milk at the fifth week of lactation is shown in table  there were significant effects of forage type on c17: cis- c18: c18: cis- cis- cla and c20: cis- cis- contents. the remaining milk fa was not significantly affected by forage type. ewes belonging to the hay group exhibited a higher content of c17: cis- and c18: cis- cis- than the gre group. by contrast, the c18: cla and c20: contents were higher in the gre group. 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: n- 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: n- 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, although most of these studies compared the effect of forage vs. concentrate diets. 1within a row, means without a common superscript differ. 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 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, as well as increasing pufa n-. in the present study, grazing increased milk cla content. in this sense, fresh pasture and alfalfa are richer in c18:3n-, 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. forage preservation affected only the lcfa and pufa n- groups. 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. total pufa n- milk fatty acids were greater in the hay group than in the gre group due to the greater amount of linoleic acid 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. 1within a row, means without a common superscript differ. se is the standard error. 2hay group: ewes received hay from mountain pastures; gre group: ewes grazed low mountain pastures. 3sfas 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 were calculated as the sum of c4:0 + c6:0 + c8:0 + c10:0; mcfas were calculates as the sum of c12:0 + c14:0; lcfas were calculated as the sum of c16:0 + c18:0 + c20:0 + c22:0 + c24:0; ufas were calculated as the sum of mufa + pufa; mufas were calculated as the sum of c16:1 + c17:1 + c18:1n-9 + c18:11n-7 + c20:1; pufas 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: 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, g6pdh and rpl. when analysed with normfinder, the rpl and gapdh genes still exhibited higher levels of gene expression stability. to normalise the data, we used the geometric mean of the rpl 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 and the effect of feeding system on gene expression in the present study, we analysed the expression of genes related to lipid metabolism within the mammary gland. we found significant differences only in the expression of the cpt1b gene. 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 -fold lower in hay ewes than in gre ewes. the enzyme carnitine palmitoyltransferase is located in the outer mitochondrial membrane and represents the main site of control for the entry of lcfa into the mitochondria. cpt 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 cpt cpt 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 cpt 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. 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-. the induction of the cpt1b gene occurs via the fatty acid response element or peroxisome proliferator-activated receptor response elements, which bind the transcription factor ppar. 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-cpt is subject to regulation at the transcriptional level and to acute control by malonyl-coa. the n-terminal domain of the m-cpt enzyme is essential for malonyl-coa inhibition. rasmussen and winder suggested that muscle contraction inhibits acetyl-coa carboxylase α and would be one mechanism for increased fatty acid oxidation. a decrease in malonyl-coa promotes fatty acid oxidation because m-cpt becomes uninhibited, whereas an increase in malonyl-coa levels 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. 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. srebf 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: cis- trans- and long chain fatty acids 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 rpl 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. authors’ contributions ed conducted the research and the statistical analysis. ed, mj and jhc wrote the manuscript. as and jar contributed to the statistical analysis and the manuscript revision. fm performed the feed and milk fatty acid analysis. mj and jhc designed the research. mj provided animals. jhc has primary responsibility for the final content. all of the authors contributed to the manuscript discussion. all of the authors read and approved the final manuscript. 
alterations in lipid metabolism occur when animals are exposed to different feeding systems. the aim of this study was to evaluate the effect of different forage diets 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. small changes in diet, such as the forage preservation, can affect the milk fatty acid profile and the expression of the cpt1b gene, which is associated with the oxidation of fatty acids. 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. 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 profile in the mammary gland and aided in the elucidation of the mechanisms of the response to diet. more stable housekeeping genes in the ovine mammary gland that would be appropriate for use in gene expression studies were ribosomal protein l and glyceraldehyde- 3- phosphate dehydrogenase. the most highly expressed gene in the mammary gland at the fifth week of lactation is caat/ enhancer- binding protein beta, possibly due to its role in milk fat synthesis in the mammary gland. when compared to hay fed indoors, grazing fresh low mountain pastures stimulates the milk content of cla and lcfa via mammary uptake. this may suggest that these fas act as regulatory factors for the transcriptional control of the carnitine palmitoyltransferase 1b gene, which was more highly expressed in the grazing group. in this study, the forage type affected the c18: cis- trans- and long-chain saturated fatty acid content, with higher percentages during grazing than during hay feeding.
 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 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 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: cis- trans- and long-chain saturated fatty acid 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 gene, which was more highly expressed in the grazing group. the most highly expressed gene in the mammary gland at the fifth week of lactation is caat/ enhancer- binding protein beta, 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 l and glyceraldehyde- 3- phosphate dehydrogenase.  small changes in diet, such as the forage preservation, 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. 