Effect of dietary restriction and subsequent re-alimentation on the transcriptional profile of hepatic tissue in cattle
© Keogh et al. 2016
Received: 12 November 2015
Accepted: 8 March 2016
Published: 17 March 2016
Compensatory growth (CG) is an accelerated growth phenomenon observed in animals upon re-alimentation following a period of dietary restriction. It is typically utilised in livestock systems to reduce feed costs during periods of reduced feed availability. The biochemical mechanisms controlling this phenomenon, however, are yet to be elucidated. This study aimed to uncover the molecular mechanisms regulating the hepatic expression of CG in cattle, utilising RNAseq. RNAseq was performed on hepatic tissue of bulls following 125 days of dietary restriction (RES) and again following 55 days of subsequent re-alimentation during which the animals exhibited significant CG. The data were compared with those of control animals offered the same diet on an ad libitum basis throughout (ADLIB). Elucidation of the molecular control of CG may yield critical information on genes and pathways which could be targeted as putative molecular biomarkers for the selection of animals with improved CG potential.
Following a period of differential feeding, body-weight and liver weight were 161 and 4 kg higher, respectively, for ADLIB compared with RES animals. At this time RNAseq analysis of liver tissue revealed 1352 significantly differentially expressed genes (DEG) between the two treatments. DEGs indicated down-regulation of processes including nutrient transport, cell division and proliferation in RES. In addition, protein synthesis genes were up-regulated in RES following a period of restricted feeding. The subsequent 55 days of ad libitum feeding for both groups resulted in the body-weight difference reduced to 84 kg, with no difference in liver weight between treatment groups. At the end of 55 days of unrestricted feeding, 49 genes were differentially expressed between animals undergoing CG and their continuously fed counterparts. In particular, hepatic expression of cell proliferation and growth genes were greater in animals undergoing CG.
Greater expression of cell cycle and cell proliferation genes during CG was associated with a 100 % recovery of liver weight during re-alimentation. Additionally, an apparent up-regulation in capacity for cellular protein synthesis during restricted feeding may contribute to and sustain CG during re-alimentation. DEGs identified are potential candidate genes for the identification of biomarkers for CG, which may be incorporated into future breeding programmes.
As feed can account for up to 75 % of the variable costs in beef cattle production systems [1, 2], any means by which these costs may be reduced without compromising overall feed efficiency or animal performance would be of benefit to the beef industry worldwide. Compensatory growth (CG) is defined as a physiological process whereby an animal has the potential, following a period of restricted feed intake, to undergo accelerated growth upon re-alimentation . The CG phenomenon is commonly utilised by cattle producers to reduce the overwintering costs of beef cattle . However, despite extensive utilisation by producers, there is little understanding of the biological and molecular mechanisms regulating the exaggerated growth phenotype typically observed.
Although typically attributing to between 1 and 1.3 % of body-weight, the liver is a major metabolic organ, accounting for on average, 24 % of whole body energy use [5, 6]. The energy requirement arises from activities associated with absorption and transportation of nutrients for subsequent use by other tissues and also a large portion of this energy is used for the maintenance of tissue integrity and mass . Alterations in the size of the liver have been shown to be directly proportional to dietary intake . Indeed previous work, including our own data, has shown a reduction in the weight and metabolic activity of this organ during dietary restriction, which facilitates efficient coping with restricted nutrient availability, primarily through a reduction in its basal metabolic rate [9–12]. It is thought that this reduced metabolic rate may continue into the initial stages of re-alimentation and thus facilitate the CG process . During subsequent re-alimentation induced CG, the liver has been shown to be one of the most responsive tissues to re-alimentation, compensating ahead of other organs and tissues in the body [9–12]. Indeed liver tissue of the compensating animals in the current study were found to have achieved 100 % recovery following 55 days of re-alimentation, whereas overall body-weight CG index for animals undergoing CG during the same time was only 48 %. A previous microarray based examination of hepatic gene expression during feed restriction, followed by early phase re-alimentation of cattle, has been reported by Connor et al. . The authors noted alterations in the expression of genes associated with cellular division and mitochondrial function during early CG. Next generation RNAseq technology has distinct advantages over microarray technology, including sensitive unbiased detection of all expressed genes without the requirement to generate an array of probes based on known sequence as well as having a much greater dynamic range . Indeed, previous work from our own group using hepatic tissue comparing both microarray and RNAseq datasets pertaining to the same biological samples identified a greater number of DEGs through utilisation of next generation sequencing technology . Furthermore, in the study of Connor et al.  differential expression of genes was evaluated within the first 2 weeks of re-alimentation, which potentially may have been too early to identify genes associated with more sustained and lasting CG as genes identified as differentially expressed by Connor et al.  may have reflected latent effects of the previous dietary restriction phase. Therefore the objective of the current study was to examine the differential expression of hepatic genes in cattle following an industry typical period of restricted feeding (125 days) and subsequent CG using RNAseq technology. The liver was chosen as a target tissue of interest, as it is a highly metabolic organ and is clearly physically affected by restricted feeding and subsequent re-alimentation induced CG [9, 11]. Our efforts during re-alimentation were focussed within the first 60 days as this is the period where the greatest increment of overall body CG is typically observed .
All procedures involving animals were approved by the University College Dublin, Animal Research Ethics Committee and licensed by the Irish Department of Health and Children in accordance with the European Community Directive 86/609/EC.
This study was conducted as part of a larger research programme designed to physiologically characterise the effect of restricted growth and subsequent re-alimentation in Holstein Friesian bulls [9, 10]. Briefly, sixty purebred Holstein Friesian bulls with a mean (SEM) age of 479 (15) days and body-weight 370 (35) kg were blocked according to weight, age, sire and a pre-trial body-weight gain into one of two groups: (i) restricted feed allowance for 125 days (RES; n = 30) followed by ad libitum access to feed for a further 55 days or (ii) ad libitum access to feed throughout (ADLIB; n = 30). The first 125 days was denoted as Period 1 and the subsequent 55 days, Period 2. Period 1 was designed to reflect an industry typical period of dietary restriction of 125 days, whereas 55 days of re-alimentation in Period 2 was designed to capture the peak of CG expression . All animals were offered a total mixed ration diet consisting of 70 % concentrate and 30 % grass silage on a dry matter basis. All animals received the same diet throughout each period, but with different proportions offered depending on treatment group. Diets were offered individually, with the proportion of feed required based on each animal’s own individual body-weight. Animals were weighed on two days at the start of the study, at the end of Period 1 and again at the end of Period 2. Additionally, throughout the study, animals were weighed every 2 weeks during Period 1 and every week during Period 2. Weighing was at the same time each morning before fresh feed was offered. During Period 1 RES animals were managed to achieve a target mean daily growth rate of 0.6 kg/day, based on dietary energy calculations using NRC . At the end of this period 15 animals from each treatment were slaughtered. All remaining animals were slaughtered at the end of Period 2. At each time point slaughter order was randomised to account for potential confounding effects on treatment outcomes.
Hepatic tissue collection
All animals were slaughtered in an EU licensed abattoir (Euro Farm Foods, Duleek, Co. Meath). Hepatic tissue was sampled from all animals within 30 min of slaughter. All tissue samples were sampled from the same location in each liver. All surgical instruments used for tissue collection were sterilized and treated with RNA Zap prior to use (Ambion, Applera Ireland, Dublin, Ireland). Samples were washed thoroughly with sterile DPBS and immediately snap frozen in liquid nitrogen before subsequent storage at −80 °C.
RNA isolation and purification
Total RNA was isolated from liver tissue samples using the Qiagen RNeasy mini kit (Qiagen), according to the manufacturer’s instructions. Approximately 60 mg of frozen tissue was used for RNA extraction. The quantity of the RNA isolated was determined by measuring the absorbance at 260 nm using a Nanodrop spectrophotometer ND-1000 (Nanodrop Technologies, DE, USA). RNA quality was assessed on the Agilent Bioanalyser 2100 using the RNA 6000 Nano Lab Chip kit (Agilent Technologies Ireland Ltd., Dublin, Ireland). RNA quality was also verified by ensuring all RNA samples had an absorbance (A260/280) of between 1.8 and 2. RNA samples with 28S/18S ratios ranging from 1.8 to 2.0 and an RNA integrity number of between 8 and 10 were deemed to be of sufficiently high quality. High quality RNA samples were selected from 10 representative animals from each treatment within each period.
cDNA library preparation and sequencing
cDNA libraries were prepared from high quality RNA using an Illumina TruSeq RNA sample prep kit following the manufacturer’s instructions (Illumina, San Diego, CA, USA). For each sample, 3 μg of total RNA was used for cDNA preparation. Briefly, mRNA was purified from total RNA and then fragmented. First strand cDNA synthesis was performed using SuperScript II Reverse Transcriptase (Applied Biosystems Ltd.) subsequently synthesising the second strand using components of the Illumina TruSeq RNA samples prep kit. Adaptors were ligated to the cDNA which was then enriched by PCR. Final individual cDNA libraries were validated on the Agilent Bioanalyser 2100 using the DNA 1000 Nano Lab Chip kit, ensuring that library fragment size was ~260 bp and library concentration was >30 ng/μl. After quality control procedures, individual RNAseq libraries were pooled based on their respective sample-specific-6 bp adaptors and sequenced at 100 bp/sequence single-end reads using an Illumina HiSeq 2000 sequencer. Approximately 16 million sequences per sample (Mean ± SD = 15,964,874 ± 1,903,207) were generated.
RNAseq data analyses
Raw sequence reads were first checked for quality using FASTQC software (version 0.10.0). Input reads were then aligned to the bovine reference genome (UMD3.1) using TopHat (v2.0.9). The software package HTSeq (v0.5.4p5) (http://pypi.python.org/pypi/HTSeq) was employed to calculate the number of sequence reads overlapping all protein coding genes from the ENSEMBL v74 annotation of the bovine genome. The number of read counts mapping to each annotated gene from HTSeq was then collated into a single file and used for subsequent differential gene expression. Only uniquely mapped reads were used for subsequent differential gene expression analysis. The R (v2.14.1) Bioconductor package, EdgeR (v3.4.1), which uses a negative binomial distribution model to account for both biological and technical variation, was applied to identify statistically significant differentially expressed genes (DEGs). Reads were first filtered before subsequent differential gene expression analysis, a gene was deemed to be expressed if the number of reads per gene per animal was ≥4. The analysis was undertaken using moderated tagwise dispersions. DEGs are defined as having a Benjamini and Hochberg false discovery rate of < 0.05 % and a fold change cut-off of 1.25 was used for each time-point.
In RNAseq experiments the differences in transcript length can yield different levels of total reads, even if transcripts are expressed at the same level. GOseq is an application for performing gene ontology analysis on RNAseq data while appropriately incorporating the effect of this transcript length selection bias . Biological pathways that were over-represented (P < 0.05) among DEGs were identified using the GOseq software package (v.1.14.0 and Kyoto Encyclopaedia of Genes and Genomes (KEGG)). Pathways were deemed over-represented when there were more DEGs in the pathway than would be expected given the size and gene length distribution . Due to the incomplete functional annotation of the bovine genome, to facilitate GOseq analysis, the online tool BioMart (www.ensembl.org/biomart/martview) was used to convert bovine gene IDs to their human orthologs. The resultant set of DEGs was then applied to test KEGG pathways (http://www.genome.jp/kegg/pathway.html) for over- or under-representation. The significant KEGG pathway maps were examined for significant DEGs. To examine the molecular functions and genetic networks, the RNAseq data were further analysed using Ingenuity Pathway Analysis (v. 8.8, Ingenuity Systems, Mountain View, CA; http://www.ingenuity.com), a web-based software application that enables identification of over-represented biological mechanisms, pathways and functions most relevant to experimental datasets or genes of interest [18–21].
qRT-PCR validation of RNAseq data
The RNAseq results were validated against gene expression values obtained from the same animals used in the current study on component genes of the somatotropic axis which has been described previously by Keogh et al. . These genes represented genes that were identified as both significantly differentially expressed as well as those not affected by either dietary restriction and subsequent re-alimentation induced CG. Briefly, using the same RNA samples that were analysed in the current RNAseq study, cDNA was synthesised and the expression of genes of the somatotropic axis examined using qRT-PCR following both dietary restriction and subsequent re-alimentation. Expression levels of candidate genes (SOCS3, JAK2, STAT5B, IGF1, IGFBP1-6, ALS, GHR1A) were normalised against expression values of selected hepatic reference genes (ACTB, CAP1). Gene expression data were checked for normality using the UNIVARIATE procedure of SAS (SAS Inst. Inc., Cary, NC). Where necessary, data were transformed using the Transreg procedure by raising values to the power of λ. Data were analysed using mixed models methodology (PROC MIXED, SAS). The Tukey critical difference test was performed to determine the existence of statistical differences between treatment mean values. The CORR procedure of SAS was used to determine correlations between RNAseq and qRT-PCR data. Pearson correlation coefficients were estimated for each individual gene across all animals. A P < 0.05 was considered to be statistically significant.
mRNA read alignment and differential gene expression
Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathways that were significantly over-represented in hepatic tissue in restricted fed animals compared with ad libitum control animals following a period of dietary restriction at the end of Period 1
Enriched KEGG pathways
Over represented P value
Steroid hormone biosynthesis
Arginine and proline metabolism
Glycine, serine and threonine metabolism
Starch and sucrose metabolism
PPAR1 signalling pathway
Protein processing in endoplasmic reticulum
MAPK signalling pathway
Valine, leucine and isoleucine degradation
Insulin signalling pathway
Differentially expressed hepatic genes involved in nutrient transport in RES compared with ADLIB animals at the end of Period 1
Solute carrier family 36 (proton/amino acid symporter), member 1
Solute carrier family 38, member 2
Solute carrier family 38, member 4
Solute carrier family 6 (amino acid transporter), member 14
Solute carrier family 7 (cationic amino acid transporter, y + system), member 2
P < 0.001
Solute carrier family 7 (amino acid transporter light chain, bo,+ system), member 9
Solute carrier family 27 (fatty acid transporter), member 4
P < 0.001
Solute carrier family 2 (facilitated glucose/fructose transporter), member 5
P < 0.001
Solute carrier family 5 (sodium/glucose cotransporter), member 1
Solute carrier family 37 (glucose-6-phosphate transporter), member 4
Solute carrier family 30, member 10
Solute carrier family 30 (zinc transporter), member 6
Solute carrier family 41 (magnesium transporter), member 2
Hepatic genes involved in the cell cycle differentially expressed following a period of dietary restriction (Period 1) and a subsequent period of re-alimentation and compensatory growth (Period 2)
Cdk5 and Abl enzyme substrate 1
Cell division cycle 7
Cyclin-dependent kinase 11A
Cyclin-dependent kinase 12
Cyclin-dependent kinase 2 associated protein 1
Cyclin-dependent kinase 2 associated protein 2
Cyclin-dependent kinase inhibitor 1B (p27, Kip1)
CDC28 protein kinase regulatory subunit 2
DNA-damage-inducible transcript 3
Growth arrest and DNA-damage-inducible, gamma interacting protein 1
Leucine zipper, putative tumor suppressor 2
NIMA-related kinase 9
NUF2, NDC80 kinetochore complex component
DNA-damage-inducible transcript 3
Hepatic genes involved in cell growth and proliferation differentially expressed following a period of dietary restriction (Period 1) and a subsequent period of re-alimentation and compensatory growth (Period 2)
Dual-specificity tyrosine-(Y)-phosphorylation regulated kinase 1A
Dual-specificity tyrosine-(Y)-phosphorylation regulated kinase 1B
Epidermal growth factor receptor
Fibroblast growth factor receptor 4
Inhibitor of DNA binding 2, dominant negative helix-loop-helix protein
Inhibin, beta C
Mesencephalic astrocyte-derived neurotrophic factor
Phosphatidylinositol-4-phosphate 3-kinase, catalytic subunit type 2 gamma
Zinc finger protein 516
Mesencephalic astrocyte-derived neurotrophic factor
Secreted protein, acidic, cysteine-rich (osteonectin)
qRT-PCR validation of RNAseq data
Genes validated between qRT-PCR and RNAseq methodologies, including fold change (FC), P-values and correlation coefficients (R)
The phenomenon of CG following a period of restricted feeding has been associated with an improved efficiency and utilisation of feed, most notably through greater overall gain in body-weight for a given level of feed intake. Indeed this was observed in the animals used in the current study, whereby animals undergoing CG were more feed efficient . Several studies have investigated the physiological control of restricted feeding and CG [9–13, 24–28] however, understanding of the molecular control of CG still remains to be elucidated. This is important as genes contributing to the phenomenon may represent putative biological markers for improved growth and feed efficiency potential in cattle. The liver was chosen as a target tissue of interest, as it is a highly metabolic organ and is clearly physically affected by restricted feeding and subsequent re-alimentation [9, 11]. Indeed liver tissue of the compensating animals in the current study were found to have achieved 100 % recovery following 55 days of re-alimentation, whereas overall body-weight and carcass CG indexs  for RES animals were 48 and 32 % respectively, during the same time.
Metabolism and nutrient transporters
As the liver is the central organ responsible for intermediary metabolism within the body , it is not surprising that a reduction in dietary intake would incur a reduction in its size and metabolic activity . Following digestion of food, the chemical constituents of the feed are transferred to the liver from the small intestine where they are further metabolised. The liver can subsequently orchestrate flux and inter-conversion of nutrients and metabolites to support the changes in demand and supply of nutrients during periods of restrictive feeding. This may be controlled through alterations in gene expression, enzyme activities and the resultant nutrient fluxes which are essential for optimal liver function and nutrient inter-conversion [31, 32]. Furthermore, as a highly metabolic organ, the liver has a substantial basal energy demand, responsible for between 18 and 25 % of the total oxygen consumption in cattle . In order to utilise energy more effectively and efficiently during restricted feeding, the liver has the capacity to regulate its size and metabolic activity so as to reduce energy requirements during times of limited nutrient availability . Indeed we and others have shown that liver size is reduced during periods of restricted feeding [9, 11, 13]. Furthermore, in the current study we observed genes associated with amino acid, lipid and carbohydrate metabolism to be differentially expressed, with a large number of these genes down-regulated following a period of dietary restriction at the end of Period 1 (Fig. 4). Following 55 days of re-alimentation, lipid and carbohydrate metabolic processes were significantly affected by re-alimentation, with the majority of DEGs involved in both lipid and carbohydrate metabolism up-regulated in RES compared with ADLIB. These results indicate an acquired greater capacity for hepatic metabolic processes during CG, which is not surprising given the documented greater feed intake per unit of bodyweight during re-alimentation which was evident in the animals used in the current study as outlined by Keogh et al. . Indeed, Burrin et al.  observed that metabolic processes associated with dietary energy intake were enhanced during CG in sheep. This was particularly evident in that study through a comparison of liver energy consumption in rams undergoing either maintenance (22 % energy consumption of whole body energy use) or CG (41 % energy consumption of whole body energy use).
This lower nutrient requirement of the liver to process nutrients was evidenced in the current study through down-regulation in the expression of a number of nutrient transporter genes in animals following feed restriction compared with their ad libitum fed contemporaries. Nutrient flux through the plasma membrane is facilitated by nutrient transporters. These trans-membrane proteins are substrate specific and are differentially expressed between different tissues to aid in the partitioning of nutrients. Following a period of feed restriction, genes associated with the transport of nutrients such as lipids, amino acids and carbohydrates as well as mineral transporters were differentially expressed. All DEGs related to nutrient transport displayed reduced expression at the end of Period 1 in RES animals compared with ADLIB animals, with the exception of SLC38A2 and SLC6A14. Expression of SLC6A14 has previously been shown to be up-regulated in the duodenal epithelia of cows which displayed greater feed and production efficiencies . A reduction in hepatic expression of nutrient transporters following restricted feed intake has also been observed in chickens . Additionally, alterations in solute carrier transporters was observed in dairy cows during negative energy balance, a period of time after calving when energy consumption is typically less than requirements . Amino acid, sugar and mineral transporters were down regulated due to alterations in energy partitioning during the early post-partum period in that study, reinforcing the role of transporters in energy partitioning in highly metabolic organs such as in the liver.
Cellular proliferation and growth
Effects of plane of nutrition on liver size and growth may be due to alterations in cellular proliferation in addition to an overall metabolic activity or workload [13, 14, 25]. We observed a number of growth and proliferative genes to be differentially expressed between feed restricted and non-restricted animals. Of note, all of these genes had lower transcript abundance in the feed restricted animals. These included genes encoding cell receptors involved in cellular growth including epidermal growth factor receptor (EGFR)  and fibroblast growth factor receptor (FGFR4) . The transcription factor ZNF516 which is involved in cell proliferation was also down regulated in RES animals at the end of Period 1. Signalling processes involved in cellular proliferation and growth such as the dual specificity tyrosine phosphorylation regulated kinases DYRK1A and DYRK1B  and PI3-kinase signalling including PIK3C2G  all had reduced expression at the end of a period of restricted feeding. The TGF-beta superfamily is involved in many cellular processes including cell growth and differentiation . A gene which codes for the beta C chain of inhibin (INHBC), which is a member of this superfamily, displayed reduced expression in the feed restricted animals. Additionally, there was up-regulation of genes associated with the inhibition of cellular proliferation, including: MANF and ID2. A previous examination of liver gene expression using microarray technology also showed alterations in genes associated with cellular growth and proliferation in hepatic tissue of steers following restricted feeding . However, there was no consistency in the specific genes identified between the two studies. This may be due to differences in experimental design between the two studies in addition to utilising varying technologies, or alternatively, it may be due to sampling tissue at different time points, or stringency in data analysis, ultimately confounding comparison of the outcomes. However, in that study, the authors subsequently noted greater expression of hepatic genes associated with cellular proliferation and growth in compensating animals early into re-alimentation (days 1 and 14 of re-alimentation). A similar finding was observed in our study with cell cycle and growth processes up-regulated at the end of Period 1 (Fig. 4). One gene in particular SPARC, which codes for a cysteine-rich acidic matrix-associated protein, was up-regulated in both hepatic tissue studies, during the initial stages of accelerated growth  and also following 55 days of re-alimentation in the current study. This gene appears to regulate cell growth through interactions with the extracellular matrix and cytokines . SPARC may be a potential genomic target for enhanced CG or improved feed efficiency potential in cattle particularly as it was differentially expressed during both the initiation of re-alimentation  and also during more sustained CG in the current study.
Continued increased expression of cellular proliferation genes by 55 days into the re-alimentation period was somewhat unexpected as by then, at least on a weight recovery basis, the liver appeared to have compensated fully (displayed a 100 % weight recovery index). However, as only 2 genes associated with proliferation and growth were up-regulated by day 55 of re-alimentation compared with a larger number observed earlier into re-alimentation in other studies , this suggests that the accelerated growth of this organ had declined by day 55 of re-alimentation which is consistent with overall recovery of this tissue. However, it must be noted that although no difference was apparent in the weight of the liver between treatment groups at the end of Period 2, a return to equal mass of the liver may not reflect a return to equal function. Thus, further evaluations on the functional control of hepatic tissue during CG is warranted.
Genes involved in the cell cycle were also differentially expressed in liver of animals undergoing feed restriction and subsequent compensation, relative to their ad libitum fed contemporaries. Genes coding for proteins important to the G1/S transition of the cell cycle were also down regulated in RES animals. For example, NEK9, which is a regulator of mitotic progression, participating in the control of spindle dynamics and chromosome separation  was down regulated in RES animals. Additionally, CDC7 which encodes a cell division cycle protein with kinase activity that phosphorylates critical substrates regulating the G1/S phase transition . Down regulation of genes involved in the G1/S transition occurred in parallel with up regulation of GADD45GIP1- a nuclear-localised protein that may be induced by p53 and regulates the cell cycle by inhibiting G1 to S phase progression . The encoded protein acts as a negative regulator of G1 to S phase progression by inhibiting cyclin-dependent kinases.
Genes coding for structural components of the cell cycle were also affected by the feed restriction regimen employed here. For example, NUF2 which codes for a component of the essential kinetochore-associated NDC80 complex, which is required for chromosome segregation and spindle checkpoint activity was down regulated in the restricted animals. The protein encoded by this gene is required for kinetochore integrity and the organisation of stable microtubule binding sites in the outer plate of the kinetochore . The effect of restricted feeding on hepatic cell cycle progression was further established through up-regulation of genes associated with cell cycle inhibition including, DDIT3 and LZTS2. DDIT3, a transcription factor that induces cell cycle arrest and apoptosis , was up regulated in feed restricted animals at the end of Period 1. Differences in cell cycle genes followed a similar pattern to those involved in cellular proliferation. In particular, with respect to the up regulation of proliferative genes 55 days into re-alimentation, where DDIT3, was subsequently down regulated in animals undergoing CG.
Cyclins are a family of proteins that control the progression of the cell cycle by activating cyclin-dependent kinase enzymes . Cyclin G2 (CCNG2), cyclin-dependent kinase 11a (CDK11A), cyclin-dependent kinase 12 (CKS2) and a cyclin regulatory subunit (CDK12) were all down regulated in restricted animals following 125 days of feed restriction. Additionally, CABLES1, which encodes a protein involved in regulating the cell cycle through interactions with cyclin-dependent kinases was also down-regulated in RES animals at the end of Period 1. Down regulation of these cell cycle progression genes further implies less cell cycle division and replication taking place in hepatic tissue of the restricted animals compared with the ad libitum fed control animals. However, CCDN3 which forms a complex with, and functions as a regulatory subunit of cyclin dependent kinase 4 and 6, whose activity is required for the G1/S cell cycle transition, was in fact up-regulated in animals following a period of restricted feeding. This result suggests that the GADD45GIP1 gene was in fact expressed in order to cause a disruption to the G1/S phase transition in the cell cycle. Indeed inhibitors of cyclin activity were identified as differentially expressed following a period of restricted feeding. Proteins encoded by CDK2AP1 and CDK2AP2 are both thought to function as negative regulators of cyclin dependent kinase 2, during S phase of the cell cycle . Similarly, CDKN1B which encodes a cyclin-dependent kinase inhibitor was also up regulated in animals fed a restricted diet. The encoded protein binds to and prevents the activation of cyclin E-CDK2 or cyclin D-CDK4 complexes and thus controls the cell cycle progression at G1 and is ultimately involved in G1 arrest. Overall, it is apparent that the documented difference in liver weight and volume at the end of a restricted feeding regimen may have been due to a reduction in the occurrence of hepatic cell division. As only one gene associated with cell division was differentially expressed 55 days into re-alimentation, together with the observed full recovery of liver weight, suggests that most, if not all, of hepatic tissue compensation had occurred at this time.
Perhaps the most striking result from this dataset was the large number of DEGs with denoted ribosomal functions (Additional file 1: Table S1). This was also evident through KEGG pathway analysis, where at the end of Period 1, the ribosome was identified as the second most significantly over-represented pathway (Table 1). The ribosome is a large and complex molecular machine that serves as the primary site of biological protein synthesis . This organelle works to link amino acids together in the order specified by messenger RNA molecules and consists of two components, the small ribosomal subunit (40S) which reads the RNA, and the large subunit (60S) which joins amino acids to form a polypeptide chain . During differential feeding in Period 1, 28 genes coding for components of the ribosomal 40S subunit and 35 genes of the 60S ribosome were all up-regulated (Additional file 1: Table S1). Greater expression of these genes in Period 1 coincided with greater expression of a number of genes involved in amino acid synthesis and protein processing. Additionally, a number of tRNA (transfer RNA) genes were also up regulated following a period of restricted feeding. Amino acids are selected, collected and carried to the ribosome by tRNA, which enter the ribosome and bind to the mRNA chain. Increased protein synthesis has been documented previously in rodents when examining the effect of feed restriction on ageing properties [51–53]. Greater expression of genes related to protein synthesis in RES animals suggests a greater efficiency and utilisation of diet derived nutrients in hepatic tissue in these animals during dietary restriction.
A greater degree of protein deposition occurs within the body during the initial stages of CG [54–56]. This is thought to act through a necessity to increase the metabolic capacity of organs such as the liver and gastrointestinal tract in order to be able to process the greater quantities of nutrients available during re-alimentation. There is potential, given the up-regulation of hepatic genes associated with protein synthesis during restricted feeding, that this may continue on into the re-alimentation period and ultimately contribute to the occurrence of CG through repletion of metabolically important tissues. In their examination of hepatic DEGs, Connor et al.  observed up-regulation of ribosomal genes on the first day of re-alimentation. Of the ribosomal genes identified in that study, 19 genes coding for subunits of the large ribosome, whilst 15 encoded subunits of the small 40S subunit were up-regulated in animals undergoing CG on the first day of re-alimentation. In our own study we identified 14 60S ribosomal genes which were in agreement with the findings of Connor et al. , whilst 14 40S ribosomal genes were also common between the two studies. By day 55 of re-alimentation, only one ribosomal gene, RPS27, was differentially expressed between RES and ADLIB animals and this was down-regulated in RES. Taken together with evidence presented earlier, it is appropriate to suggest that hepatic compensation was complete at this stage of re-alimentation. When metabolically important tissues have been restored fully, there is an apparent increase in adipose deposition as opposed to protein tissue deposition [55, 57–59]. An indication towards an increase in adipose deposition occurring during CG was obtained through up-regulation of genes involved in adipose deposition including FADS1 and SREBF1on day 55 of re-alimentation. These results potentially indicate an increase in adipose deposition together with a decrease in protein deposition in hepatic tissue coinciding with a decrease in overall body growth rate at the end of Period 2.
Potential molecular biomarkers
The CG phenomenon is utilised in beef production systems worldwide [60, 61]. However knowledge of the underlying molecular control regulating the expression of CG is lacking. A greater understanding of the genetic basis for CG is critical to the future effective exploitation of the trait and may lead to the discovery of DNA-based biomarkers which could be incorporated into genomic selection breeding programmes to select animals with a greater propensity to display CG following prior dietary restriction. Furthermore, as CG is associated with an improvement in feed efficiency, differentially expressed genes identified in this study may contribute to breeding protocols for the selection of animals with improved feed efficiency. An examination of the findings of both the current study and those of Connor et al. , show that a number of DEGs were observed to be in agreement during restricted feeding in the current dataset and during early CG  including: AKR1C3; INSIG1; SELK and UBL5. These genes may be potential targets for the accelerated growth observed during the early stages of CG. Additionally, hepatically expressed genes which had greater transcript abundance in feed efficient (low residual feed intake) animals including GOLTA1, IDH2, INHBA, PSPH, PYCR1, RPS4X and STEAP4  were also up-regulated in the present study during restricted feeding and may represent markers for feed efficiency. During CG, SPARC was up-regulated in both the current study and the data of Connor et al. . Only one common gene between these two studies during the CG phase may be due to different sampling time points, with samples in the current study taken on day 55 of re-alimentation and samples taken by Connor et al.  on the first day of re-alimentation. Additionally, the gene SREBF1 was identified as up-regulated in the current hepatic data-set on day 55 of re-alimentation, as well as on day 15 of re-alimentation in skeletal muscle in the same animals undergoing CG as in the current study . SREBF1 codes for a sterol regulator element-binding transcription factor and has a crucial role in energy homeostasis through promotion of glycolysis, lipogenesis and adipogenesis [64–66]. Further investigation is warranted to determine if these results can be utilised as potential molecular markers for CG and feed efficiency in cattle.
During dietary restriction evidence of reduced metabolic activity of the liver was apparent through less mRNA abundance of nutrient transporters. Additionally, through an examination of differential transcript abundance we observed evidence for a reduction in cell cycle processes as well as a reduction in cellular proliferation and growth. For both processes this was manifested as a down-regulation in the expression of promoter genes coincident with an increase in transcripts coding for inhibitor proteins. Our data also suggests that an increased capacity for protein synthesis following feed restriction may sustain into CG in re-alimentation. This may then allow for the preferential deposition of protein ahead of adipose deposition in order to aid in the recovery of metabolically active organs such as the liver and gastrointestinal tract, which are required to facilitate the increased feed intake and nutrient availability associated with re-alimentation. Our results further suggest that hepatic protein synthesis may have subsided by 55 days into re-alimentation, in favour of adipose deposition and may signal full hepatic tissue recovery at that stage. In the context of the current study, 55 days into re-alimentation may have been too late to identify large differences in pathways and genes regulating CG, as although an overall body compensatory index of only 48 % was achieved in this time, the liver had fully recovered by day 55. Potentially sampling hepatic tissue earlier into re-alimentation for example following 1 month of re-alimentation may yield further information on the underlying biological control regulating the expression of CG. However, genes differentially expressed by day 55 in the current study may represent a more sustained or prolonged CG. It must be noted, however, that a return to equal mass of the liver following 55 days of re-alimentation and CG may not reflect a return to equal function, thus, further evaluations on the functional control of hepatic tissue during CG is warranted. The new knowledge generated in this study offer further insights into some of the molecular processes underlying restricted and CG in cattle. Furthermore, differential gene expression patterns provide data which may be further integrated and used for the selection of robust biomarkers to identify animals with superior genetic potential for CG and feed efficiency.
The authors wish to acknowledge funding from Science Foundation Ireland (09/RFP/GEN2447).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- NRC. Nutrient requirements of beef cattle. 7th ed. Washington, DC: National Academy Press; 2000.Google Scholar
- Finneran E, Crosson P, O’Kiely P, Shalloo L, Forristal D, Wallace M. Simulation modeling of the cost of producing and utilizing feeds for ruminants on Irish farms. J Farm Management. 2010;14:95–116.Google Scholar
- Hornick JL, Van Eenaeme C, Gerard O, Dufrasne I, Istasse L. Mechanisms of reduced and compensatory growth. Domest Anim Endocrinol. 2000;19(2):121–32.View ArticlePubMedGoogle Scholar
- Keane MG, Drennan MJ. Effects of winter supplementary concentrate level on the performance of steers slaughtered immediately or following a period at pasture. Ir J Agric Res. 1994;33(2):111–9.Google Scholar
- McBride BW, Kelly JM. Energy cost of absorption and metabolism in the ruminant gastrointestinal tract and liver: a review. J Anim Sci. 1990;68:2997–3010.PubMedGoogle Scholar
- Caton JS, Bauer ML, Hidari H. Metabolic components of energy expenditure in growing beef cattle –review. Asian-Australas J Anim Sci. 2000;13:702–10.View ArticleGoogle Scholar
- Baldwin VIRL, McLeod KR, Klotz JL, Heitmann RN. Rumen development, intestinal growth and hepatic metabolism in the pre- and postweaning ruminant. J Dairy Sci. 2004;87:E55–65.View ArticleGoogle Scholar
- Johnson DE, Johnson KA, Baldwin RL. Changes in liver and gastrointestinal tract energy demands in response to physiological workload in ruminants. J Nutr. 1990;120:649–55.PubMedGoogle Scholar
- Keogh K, Waters SM, Kelly AK, Kenny DA. Feed restriction and compensatory growth in Holstein Friesian bulls: I. Effect on animal performance, muscle, fat and linear body measurements and slaughter characteristics. J Anim Sci. 2015;93(7):3578–89.View ArticlePubMedGoogle Scholar
- Keogh K, Waters SM, Kelly AK, Wylie ARG, Sauerwein H, Sweeney T, Kenny DA. Feed restriction and compensatory growth in Holstein Friesian bulls: I. Effect on blood pressure and circulating levels of metabolites and metabolic hormones. J Anim Sci. 2015;93(7):3590–601.View ArticlePubMedGoogle Scholar
- Yambayamba ESK, Price MA, Jones SDM. Compensatory growth of carcass tissues and visceral organs in beef heifers. Livest Prod Sci. 1996;46(1):19–32.View ArticleGoogle Scholar
- Yambayamba ES, Price MA, Foxcroft GR. Hormonal status, metabolic changes, and resting metabolic rate in beef heifers undergoing compensatory growth. J Anim Sci. 1996;74(1):57–69.PubMedGoogle Scholar
- Ryan WJ, Williams IH, Moir RJ. Compensatory growth in sheep and cattle II. Changes in body composition and tissue weights. Aust J Agric Res. 1993;44(7):1623–33.View ArticleGoogle Scholar
- Connor EE, Kahl S, Elsasser TH, Parker JS, Li RW, Van Tassell CP, Baldwin RL VI, Barao SM. Enhanced mitochondrial complex gene function and reduced liver size may mediate improved feed efficiency of beef cattle during compensatory growth. Funct Integr Genomics. 2010;10(1):39–51.View ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.View ArticlePubMedPubMed CentralGoogle Scholar
- McCabe M, Waters S, Morris D, Kenny D, Lynn D, Creevey C. RNA-seq analysis of differential gene expression in liver from lactating dairy cows divergent in negative energy balance. BMC Genomics. 2012;13:193–203.View ArticlePubMedPubMed CentralGoogle Scholar
- Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-SEQ: accounting for selection bias. Genome Biol. 2010;11(2):R14.View ArticlePubMedPubMed CentralGoogle Scholar
- Barilli A, Rotoli BM, Visigalli R, Bussolati O, Gazzola GC, Kadija Z, Rodi G, Mariani F, Ruzza ML, Luisetti M, Dall’Asta V. In lysinuric protein intolerance system y + L activity is defective in monocytes and in GM-CSF-differentiated macrophages. Orphanet J Rare Dis. 2010;5:32.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu X, Milo M, Lawrence ND, Rattray M. Probe-level measurement error improves accuracy in detecting differential gene expression. Bioinformatics. 2006;22(17):2107–13.View ArticlePubMedGoogle Scholar
- Loor JJ, Everts RE, Bionaz M, Dann HM, Morin DE, Oliveira R, Rodriquez-Zas SL, Drackly JK, Lewin HA. Nutrition-induced ketosis alters metabolic and signalling gene networks in liver of periparturient dairy cows. Physiol Genomics. 2007;32(1):105–16.View ArticlePubMedGoogle Scholar
- Selvey S, Thompson EW, Matthaei K, Lea RA, Irving MG, Griffiths LR. Beta-actin –an unsuitable internal control for RT-PCR. Mol Cell Probes. 2001;15(5):307–11.View ArticlePubMedGoogle Scholar
- Keogh K, Waters SM, Kelly AK, Wylie ARG, Kenny DA. Effect of feed restriction and subsequent re-alimentation on hormones and genes of the somatotropic axis in cattle. Physiol Genomics. 2015;47(7):264–73.View ArticlePubMedGoogle Scholar
- Goeman JJ, Buhlmann P. Analyzing gene expression data in terms of gene sets: methodological issues. Bioinformatics. 2007;23(8):980–7.View ArticlePubMedGoogle Scholar
- Hayden JM, Williams JE, Collier RJ. Plasma growth hormone, insulin-like growth factor, insulin, and thyroid hormone association with body protein and fat accretion in steers undergoing compensatory gain after dietary energy restriction. J Anim Sci. 1993;71(12):3327–38.PubMedGoogle Scholar
- Hornick JL, Van Eenaeme C, Clinquart A, Diez M, Istasse L. Different periods of feed restriction before compensatory growth in Belgian Blue bulls: I. animal performance, nitrogen balance, meat characteristics, and fat composition. J Anim Sci. 1998;76(1):249–59.PubMedGoogle Scholar
- Hornick JL, Van Eenaeme C, Diez M, Minet V, Istasse L. Different periods of feed restriction before compensatory growth in Belgian Blue bulls: II. Plasma metabolites and hormones. J Anim Sci. 1998;76(1):260–71.PubMedGoogle Scholar
- Ryan WJ, Williams IH, Moir RJ. Compensatory growth in sheep and cattle I. Growth pattern and feed intake. Aust J Agric Res. 1993;44(7):1609–21.View ArticleGoogle Scholar
- Sainz RD, De la Torre F, Oltjen JW. Compensatory growth and carcass quality in growth-restricted and refed beef steers. J Anim Sci. 1995;73(10):2971–9.PubMedGoogle Scholar
- Stangassinger M, Giesecke D. Splanchnic metabolism of glucose and related energy substrates. In: Milligan LP, Grovum WL, Dobson A, editors. Control of digestion and metabolism in ruminants. Englewood Cliffs, NJ: Pretice-Hall; 1986. p. 347–66.Google Scholar
- Burrin DG, Ferrell CL, Britton RA, Bauer M. Level of nutrition and visceral organ size and metabolic activity in sheep. Br J Nutr. 1990;64:439–48.View ArticlePubMedGoogle Scholar
- Loor JJ, Dann HM, Guretzky NA, Everts RE, Oliveira R, Green CA, Litherland NB, Rodriquez-Zas SL, Lewin HA, Drackly JK. Plane of nutrition prepartum alters hepatic gene expression and function in dairy cows as assessed by longitudinal transcript and metabolic profiling. Physiol Genomics. 2006;27(1):29–41.View ArticlePubMedGoogle Scholar
- Reynolds CK, Aikman PC, Lupoli B, Humphreys DJ, Beever DE. Splanchnic metabolism of dairy cows during the transition from late gestation through early lactation. J Dairy Sci. 2003;86(4):1201–17.View ArticlePubMedGoogle Scholar
- Burrin DG, Ferrell CL, Eisemann JH, Britton RA, Nienaber JA. Effect of level of nutrition on splanchnic blood flow and oxygen consumption in sheep. BR J Nutr. 1989;62:23–34.View ArticlePubMedGoogle Scholar
- Waters SM, Keogh K, Buckley F, Kenny DA. Effect of genotype expression of nutrient transporter genes in dairy cows. J Anim Sci Biotechnol. 2013;4:49.View ArticlePubMedPubMed CentralGoogle Scholar
- Humphrey BD, Stephensen CB, Calvert CC, Klasing KC. Lysine deficiency and feed restriction independently alter cationic amino acid transporter expression in chickens (Gallus gallus comesticus). Comp Biochem Physiol A Mol Integr Physiol. 2006;143(2):218–27.View ArticlePubMedGoogle Scholar
- McCarthy SD, Waters SM, Kenny DA, Diskin MG, Fitzpatrick R, Patton J, Wathes DC, Morris DG. Negative energy balance and hepatic gene expression patterns in high-yielding dairy cows during the early postpartum period: a global approach. Physiol Genomics. 2010;42(3):188–99.View ArticleGoogle Scholar
- Herbst RS. Review of epidermal growth factor receptor biology. Int J Radiat Oncol Biol Phys. 2004;59(2):21–6.View ArticlePubMedGoogle Scholar
- Eswarakumar VP, Lax I, Schlessinger J. Cellular signaling by fibroblast growth factor receptors. Cytokine Growth Factor Rev. 2005;16(2):139–49.View ArticlePubMedGoogle Scholar
- Becker W. Emerging role of DYRK family protein kinases as regulators of protein stability in cell cycle control. Cell Cycle. 2012;11(18):3389–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Foster FM, Traer CJ, Abraham SM, Fry MJ. The phosphoinositide (PI) 3-kinase family. J Cell Sci. 2003;116:3037–40.View ArticlePubMedGoogle Scholar
- Huang SS, Huang JS. TGF-beta control of cell proliferation. J Cell Biochem. 2005;96(3):447–62.View ArticlePubMedGoogle Scholar
- Brekken RA, Sage EH. SPARC, a matricellular protein: at the crossroads of cell-matrix communication. Matrix Biol. 2001;19(8):816–27.View ArticlePubMedGoogle Scholar
- Yang SW, Gao C, Chen L, Song YL, Zhu JL, Qi ST, Jiang ZZ, Wang ZW, Lin F, Huang H, Xing FQ, Sun QY. Nek9 regulates spindle organization and cell cycle progression during mouse oocyte meiosis and its location in early embryo mitosis. Cell Cycle. 2012;11(23):4366–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Kumagai H, Sato N, Yamada M, Mahony D, Seghezzi W, Lees E, Arai K, Masai H. A novel growth- and cell cycle-regulated protein, ASK, activates human Cdc7-related kinase and is essential for G1/S transition in mammalian cells. Mol Cell Biol. 1999;19(7):5083–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Nakayama K, Nakayama N, Wang TL, Skik leM. NAC-1 controls cell growth and survival by repressing transcription of Gadd45GIP1, a candidate tumor suppressor. Cancer Res. 2007;67(17):8058–64.View ArticlePubMedGoogle Scholar
- Sundin LJ, Guimaraes GJ, Deluca JG. The NDC80 complex proteins Nuf2 and Hec1 make distinct contributions to kinetochore-microtubule attachment in mitosis. Mol Biol Cell. 2011;22(6):759–68.View ArticlePubMedPubMed CentralGoogle Scholar
- Bento C, Andersson MK, Aman P. DDIT3/CHOP and the sarcoma fusion oncoprotein FUS-DDIT3/TLS-CHOP bind cyclin-dependent kinase 2. BMC Cell Biol. 2009;10:89.View ArticlePubMedPubMed CentralGoogle Scholar
- Nigg EA. Cyclin-dependent protein kinases: key regulators of the eukaryotic cell cycle. Bioessays. 1995;17(6):471–80.View ArticlePubMedGoogle Scholar
- Sun M, Jiang R, Wang G, Zhang C, Li J, Jin C, Zhang X. Cyclin-dependent kinase 2-assoicated protein 1 suppresses growth and tumorigenesis of lung cancer. Int J Oncol. 2013;42(4):1376–82.PubMedGoogle Scholar
- Alberts BA, Johnson A, Lewis J, Martin R, Roberts K, Walter P: Molecular biology of the cell. 4th edition. Garland Publishing, Inc., New York, NY, USA.Google Scholar
- Lakowski B, Hekimi S. The genetics of caloric restriction in Caenorhabditis elegans. Proc Natl Acad Sci USA. 1998;95(22):13091–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Lewis SE, Goldspink DF, Phillips JG, Merry BJ, Holehan AM. The effects of aging and chronic dietary restriction on whole body growth and protein turnover in the rat. Exp Gerontol. 1985;20(5):253–63.View ArticlePubMedGoogle Scholar
- Ward W, Richardson A. Effect of age on liver protein synthesis and degradation. Hepatology. 1991;14(5):935–48.View ArticlePubMedGoogle Scholar
- Carstens GE, Johnson DE, Ellenberger MA, Tatum JD. Physical and chemical components of the empty body during compensatory growth in beef steers. J Anim Sci. 1991;69(8):3251–64.PubMedGoogle Scholar
- Rompala RE, Jones SDM, Buchanan-Smith JG, Bayley HS. Feedlot performance and composition of gain in late-maturing steers exhibiting normal and compensatory growth. J Anim Sci. 1985;61:637–46.Google Scholar
- Ryan WJ. Compensatory growth in cattle and sheep. Nutr Abstr Rev. 1990;60(9):653–64.Google Scholar
- Drew KR, Barry TN, Duncan SJ, Kleim C. Compensatory growth and reproductive performance in young sheep given differential nutrition from 9 to 13 months of age. New Zealand J Exper Agric. 1973;1(2):109–14.View ArticleGoogle Scholar
- Hogg BW: Effects of growth patterns on body composition and compensatory growth in sheep. PhD thesis. University of Melbourne, Australia; 1977.Google Scholar
- Wright IA, Russel AFJ. Changes in the body composition of beef cattle during compensatory growth. Anim Prod. 1991;52(1):105–13.View ArticleGoogle Scholar
- Keane MG, Drennan MJ. Effects of winter supplementary concentrate level on the performance of steers slaughtered immediately or following a period at pasture. Ir J Agric Res. 1994;33(2):111–9.Google Scholar
- Ashfield A, Wallace M, McGee M, Crosson P. Bioeconomic modelling of compensatory growth for grass-based dairy calf-to-beef production systems. J Agric Res. 2014;152:805–16.Google Scholar
- Chen Y, Gondro C, Quinn K, Herd RM, Parnell PF, Vanselow B. Global gene expression profiling reveals genes expressed differentially in cattle with high and low residual feed intake. Anim Genet. 2011;42(5):475–90.View ArticlePubMedGoogle Scholar
- Keogh K, Kenny DA, Cormican P, Kelly AK, Waters SM. Effect of dietary restriction and subsequent re-alimentation on the transcriptional profile of bovine skeletal muscle. Plos one. 2015;11(2), e0149373.View ArticleGoogle Scholar
- Hoashi S, Ashida N, Ohsaki H, Utsugi T, Sasazaki S, Taniguchi M, Oyama K, Mukai F, Mannen H. Genotype of bovine sterol regulatory element binding protein-1 (SREBP-1) is associated with fatty acid composition in Japanese Black cattle. Mamm Genome. 2007;18(12):880–6.View ArticlePubMedGoogle Scholar
- Bhuiyan MSA, Yu SL, Jeon JT, Yoon D, Cho YM, Park EW, Kim NK, Kim KS, Lee JH. DNA polymorphisms in SREBF1 and FASN genes affect fatty acid composition in Korean cattle (Hanwoo). Asian-Aust J Anim Sci. 2009;22(6):765–73.View ArticleGoogle Scholar
- Ntambi JM. Regulation of stearoyl-CoA desaturase by polyunsaturated fatty acids and cholesterol. J Lipid Res. 1999;40(9):1549–58.PubMedGoogle Scholar