Global gene expression in muscle from fasted/refed trout reveals up-regulation of genes promoting myofibre hypertrophy but not myofibre production

Compensatory growth is a phase of rapid growth, greater than the growth rate of control animals, that occurs after a period of growth-stunting conditions. Fish show a capacity for compensatory growth after alleviation of dietary restriction, but the underlying cellular mechanisms are unknown. To learn more about the contribution of genes regulating hypertrophy (an increase in muscle fibre size) and hyperplasia (the generation of new muscle fibres) in the compensatory muscle growth response in fish, we used high-density microarray analysis to investigate the global gene expression in muscle of trout during a fasting-refeeding schedule and in muscle of control-fed trout displaying normal growth. The compensatory muscle growth signature, as defined by genes up-regulated in muscles of refed trout compared with control-fed trout, showed enrichment in functional categories related to protein biosynthesis and maturation, such as RNA processing, ribonucleoprotein complex biogenesis, ribosome biogenesis, translation and protein folding. This signature was also enriched in chromatin-remodelling factors of the protein arginine N-methyl transferase family. Unexpectedly, functional categories related to cell division and DNA replication were not inferred from the molecular signature of compensatory muscle growth, and this signature contained virtually none of the genes previously reported to be up-regulated in hyperplastic growth zones of the late trout embryo myotome and to potentially be involved in production of new myofibres, notably genes encoding myogenic regulatory factors, transmembrane receptors essential for myoblast fusion or myofibrillar proteins predominant in nascent myofibres. Genes promoting myofibre growth, but not myofibre formation, were up-regulated in muscles of refed trout compared with continually fed trout. This suggests that a compensatory muscle growth response, resulting from the stimulation of hypertrophy but not the stimulation of hyperplasia, occurs in trout after refeeding. The generation of a large set of genes up-regulated in muscle of refed trout may yield insights into the molecular and cellular mechanisms controlling skeletal muscle mass in teleost and serve as a useful list of potential molecular markers of muscle growth in fish.


Background
Skeletal muscle has a number of important functions: it maintains posture, produces locomotion and is a dominant organ in energy metabolism. Skeletal muscle is a dynamic structure that forms during embryogenesis, grows afterwards and is subjected to remodeling and changes in mass under a variety of physiological conditions such as nutritional perturbations, modifications of activity or ageing. In mammals, embryonic and fetal skeletal muscle forms and grows via the proliferation, differentiation and fusion of myogenic cells, whereas postnatal muscle grows largely through remodelling of pre-existing myofibres [1]. In contrast, fish lastingly combine muscle hyperplasia (generation of new fibres) and hypertrophy (increase in fibre size) to generate indeterminate muscle growth [2,3]. Muscle hyperplasia in fish initially occurs in a discrete layer at the surface of the myotome in late embryos or early larvae. This regionalized phase of myogenesis, termed stratified hyperplasia, results from the differentiation of myogenic progenitor cells originating from the dermomyotome-like epithelium present at the surface of the embryonic myotome [4][5][6]. At post-larval stages, new myofibres are formed between existing muscle fibres throughout the myotome, thus producing a typical mosaic appearance in muscle cross sections [2]. A Myog:GFP transgenic trout model has revealed that mosaic hyperplasia is prevalent in the juvenile stage, but progressively decreases as trout age, and eventually ceases at approximately 18 months postfertilization [7]. Nevertheless, potentially recruitable muscle stem cells are present in muscles of aged trout, as shown by their ability to form myofibres de novo after muscle injury [7,8].
Although the genetic programs determining muscle development are well characterized in vertebrates, including fish [9], little is known regarding the programs involved in muscle remodelling and muscle mass changes occurring throughout life. Among adaptive responses that remodel muscle, there is compensatory growth. Compensatory growth is a phase of rapid growth, greater than the growth rate of control animals, that follows growth depression [10]. Compensatory growth is triggered by refeeding of animals after a period of weight loss induced by a prolonged fasting period. In fish, compensatory growth has been related to an increase in feed intake (hyperphagia) [11] and efficiency of food utilization [10]. Hence, there is considerable interest in exploiting compensatory growth, through the design of feeding schedule, to optimize muscle growth [12,13]. Using nylon macroarrays, we have previously reported a preliminary exploration of the temporal gene expression profiling of trout muscle during a fasting-refeeding schedule [14]. In the present study, using a high-density Agilent-based microarray platform for trout, and considering the gene expression of control-fed trout that displayed usual growth, we exhaustively defined genes specifically associated with the compensatory muscle growth response, and characterized the corresponding functional categories. Additionally, to further characterize the cellular mechanisms involved in the compensatory growth response, we compared the genes composing the compensatory muscle growth signature and the genes up-regulated in hyperplastic growth zones of the post-embryonic trout myotome, as they were previously identified using laser-capture microdissection combined with the same Agilent-based microarray platform [15].

Results
Growth characteristics in fasted, fasted-refed and control (normally) fed trout Changes in body weight and condition factor K, an indicator of body shape, are shown in Fig. 1. This figure shows a slight decrease in total body weight after fasting that was followed by an increase in total body weight after refeeding. Nonetheless, total body weight increase in refed animals was found to be comparable to that observed in control group from day 0 to day 36. The condition factor K decreased after nutrient deprivation and increased after refeeding to reach control level. By contrast, condition factor only moderately decreased at day 36. On the whole, these data show a resumption of growth after refeeding, but comparison of growth curves in refed and control trout did not objectify a burst of muscle growth after refeeding. To further provide experimental evidence of a compensatory growth response in muscle of refed trout, a global gene expression analysis was conducted in muscle of refed and control-fed trout.

Gene expression profiling overview
Anova test was used to define genes whose expression levels in trout muscle were significantly different across the various feeding conditions considered (fasted for 1 month (F0), fasted for 2 month and then refed 4, 11 or 36 days (RF4, RF11 and RF36, respectively) and normally fed (control: C0). This analysis led to the identification of approximately 2300 unique differentially expressed genes. The hierarchical clustering of differentially expressed genes, shown in Fig. 2 and available through the heat map file (Additional file 1) and Java TreeView tool (https://sourceforge.net/projects/jtreeview/files/) resulted in the formation of three major clusters with distinct expression profile. Cluster I contained genes up-regulated in muscles of fasted trout and then down-regulated after refeeding. This cluster contained markers of muscle atrophy such as FBXO32/atrogin1 and trim63 (murf1), and markers of autophagy such as SQSTM1, thus indicating the activation of proteolytic systems in muscles of fasted trout. Clusters IIa and IIb both contained genes up-regulated after refeeding. Cluster IIa included genes whose expression levels after refeeding were restored to values found in control (normally fed) trout, and cluster IIb included genes whose expression levels after refeeding exceeded those in muscles of control-fed trout. Because compensatory growth refers to a faster than usual (control-fed) growth rate that is induced after refeeding of fasted animals [16], we reasoned that the compensatory muscle growth signature, consisted only of genes contained in cluster IIb. The study presented here focused on genes up-regulated after refeeding and more especially on those forming the compensatory growth signature (cluster IIb).
Cluster IIa: Genes whose expression levels after refeeding was restored to values found in control trout Cluster IIa was composed of 422 unique genes whose expression levels during muscle recovery progressively increased to reach those of normally fed trout. Gene ontology analysis was performed to understand the functional significance of the genes contained in cluster IIa. The DAVID analysis of 357 eligible genes specific to this cluster showed enrichment in GO terms linked to cell division (mitosis (P < 1.4.10 −7 ), organelle fission (P < 2.6.10 −7 ) and chromosome condensation (P < 3.9.10 −4 )). In addition, the genes in cluster IIa showed enrichment in GO terms related to cytoskeleton organization (P < 3.7.10 −6 ), components of the sarcomere (P < 2.5.10 −5 , 12 genes including myosins, skeletal muscle actins, troponins and tropomyosins), and proteins of the extracellular matrix (P < 4.10 −11 , 31 genes including collagen chains, laminin subunits, periostin, CILP1, tenascin and many proteoglycans such as syndecan 2, mimecan and keratocan). Overall, cluster IIa showed enrichment in genes involved in cell division and genes encoding structural components of muscle fibres (for details, see Table 1 and Additional file 2 for list of genes that formed cluster IIa).
Cluster IIb: Genes whose expression levels after refeeding exceeded those in control trout Cluster IIb was composed of 1161 unique genes upregulated in muscles of refed trout compared with both fasted and control-fed trout; as such, we identified cluster IIb as the compensatory muscle growth signature. Notably, genes composing cluster IIb clearly displayed an earlier up-regulation after refeeding than genes forming cluster IIa. A DAVID analysis of 960 eligible genes showed enrichment in GO terms linked to transcription, such as RNA processing (P < 5.6.10 −53 ), rRNA metabolic process (P < 7.4.10 −30 ), mRNA metabolic process (p < 2.8.10 −19 ), tRNA metabolic process (P < 1.2.10 −17 ), and ribonucleoprotein complex biogenesis (P < 7.10 −49 ). Other GO terms associated with cluster IIb included translation (P < 4.7.10 −27 ), ribosome biogenesis (P < 9.1.10 −38 ), which determines translation capacity, cellular macromolecular complex assembly (P < 8.4.10 −16 ) and protein folding (P < 7.1.10 −20 ). Enrichment in genes involved in mitochondria biogenesis and activity (P < 2.10 ) was also found in cluster IIb. Finally, cluster IIb included many genes encoding epigenetic transcriptional regulators. Among them were the histone-lysine Nmethyltransferases EZH2, jarid2/Jumonji and RBBP4 which are all components of the PCR2 complex, the histone-lysine N-methyltransferases EHMT2, SETDB1 and WHSC1/NSD2, and many histone modifying enzymes of the protein arginine methyltransferase (PRMT) family such as Prmt1-A, Prmt1-B Prmt3; Prmt5, Prmt6 and Prmt7 (Fig. 3). Overall, the molecular signature of compensatory muscle growth showed enrichment in genes encoding histone methyltransferases, genes regulating protein biosynthesis for cell growth and genes involved in mitochondrion biogenesis for energy supply (for details, see Tables 1, 2, 3 and 4, and Additional file 3 for list of genes that that formed the major functional categories of cluster IIb).

The compensatory muscle growth signature does not include genes involved in the formation of new myofibres
To characterize the genetic mechanisms regulating the formation of new myofibres in fish we previously examined the transcriptome of the superficial hyperplastic growth zones of the late trout embryo myotome by using laser capture microdissection and microarray analysis [15]. In this study, to further characterize the cellular mechanisms involved in the compensatory growth response, we compared its molecular signature with that of the superficial hyperplastic growth zones. A Venn diagram showed that most of the genes (849 out of 1161) of the compensatory muscle growth signature were included in the list of genes up-regulated in superficial hyperplastic growth zones (Fig. 4). As a result, GO categories inferred from genes common to both situations were mostly those inferred from the compensatory growth signature, and were related to protein biosynthesis and maturation such as RNA processing, ribonucleoprotein complex biogenesis, ribosome biogenesis, translation and protein folding. Additionally, several epigenetic factors of the protein arginine N-methyl transferase (PRMT) family (i.e., Prmt1-A, Prmt-B Prmt3; Prmt5, Prmt6 and Prmt7), as well as several histone lysine Nmethyltransferases (EZH2, EHMT2, WHSC1/NSD2 and SETDB1), were found in both situations. In contrast, none of the genes up-regulated in hyperplastic growth zones and that encode SWI/SNF chromatin-remodelling enzymes (i.e., Smarcd1, Smarce1, Smarcb1A, smarca5, smarcad1, Smarcab1, Smarcc1 and Smarca4/BRG1) appeared in the compensatory growth signature. This showed that only a subset of the chromatin-remodelling factors up-regulated in hyperplastic growth zones, was present in the compensatory muscle growth signature. Unexpectedly, GO terms related to DNA replication and cell cycle and that were inferred from genes upregulated in hyperplastic growth zones were not found in the compensatory muscle growth signature. Additionally, transcriptional regulators, immunoglobulin domaincontaining transmembrane proteins and secreted/signalling molecules known or predicted to be important in the differentiation of myogenic cells or their fusion into new myofibres were virtually absent from the compensatory growth signature. For example, neither canonical myogenic transcriptional regulators (Pax3, Pax7 and MRFs such as MyoD1a, MyoD1b and MyoD1c, myf5, myogenin and mrf4), nor genes encoding transcriptional regulators known to modulate MRF activity, such as Tsh3, ARX, meis1 and pbx1, were found to be associated with compensatory growth, whereas they were upregulated in hyperplastic growth zones. Additionally, among the >30 genes found to be specific to the superficial growth areas and to encode homeobox-containing transcriptional regulators or members of the Hairy and enhancer of split family, none were included within the compensatory growth signature. Of note, Hairy and enhancer of split 6 and myogenin were found to exhibit increased expression after refeeding, but their transcription levels were never above those found in control animals. With the notable exception of M-cadherin, hyperplasiacorrelated genes encoding promyogenic membrane receptors such as N-cadherin, Brother of CDO, NCAM and Fig. 2 Hierarchical clustering of differentially expressed genes in muscle during a fasting-refeeding schedule and in control-fed trout displaying usual growth. Hierarchical clustering of differentially expressed genes led to the formation of three distinct clusters: I, IIa and IIb. Cluster IIb, which includes genes up-regulated in muscles of refed trout compared with control-fed trout, defines the specific molecular signature of compensatory muscle growth following refeeding. Each row represents the expression pattern of a single gene, and each column corresponds to a single sample: columns 1 to 5, muscles from fasted trout; columns 6 to 10, 11 to 15 and 16 to 19, muscles from 4-, 11-and 36-days refed trout respectively; columns 20 to 23, muscles of control-fed trout. The expression levels are represented by colored tags, with red representing the highest levels of expression and green representing the lowest levels of expression protogenin [17], or the essential fusion effectors Kin of Irre like-3, and Jamb [18,19] were also absent from the compensatory muscle growth signature. Similarly, with the exception of stromal cell-derived factor 2-like protein 1 and Hepatoma-derived growth factor-related protein, virtually all the hyperplasia-correlated genes encoding secreted factors were absent from the compensatory signature. Finally, because the same Agilent-based microarray was used for gene expression profiling in both experiments, we reliably determined that among the >30 distinct genes encoding myofibrillar proteins (e.g., troponins, myosin chains, myosin binding proteins, tropomyosins, alpha actins, myozenins) that were up-regulated in lasercaptured hyperplastic growth zones, only one (tropomyosin α4-chain) was found in the compensatory muscle growth signature (cluster IIb), and none were found among the numerous myofibrillar protein genes found in cluster IIa. Together, these data showed that genes common to superficial hyperplastic growth zone and to the compensatory muscle growth signature were related to protein biosynthesis for myofibre growth. In contrast, hyperplasia-correlated genes involved in myofibre production or encoding myofibrillar proteins specific to nascent myofibres were virtually excluded from the compensatory growth signature.

Discussion
Until now, very little was known regarding the muscle growth pattern during a fasting-refeeding sequence, essentially because cellularity analysis is difficult to interpret in   Fig. 2 this case. Indeed, small muscle fibres result either from the production of new fibres (muscle hyperplasia) or from the shrinking of pre-existing fibres (muscle fibre atrophy).
To gain insight into the relative contribution of hypertrophy and hyperplasia in the compensatory muscle growth response, we performed transcriptomic analysis in muscles of fasted, refed and control-fed trout and specifically investigate the expression of genes potentially involved in these two processes. Although our experimental protocol was designed to exacerbate compensatory growth (i.e., long fasting period and large excess of food during refeeding), growth curve in refed trout did not show the idealized pattern of growth compensation as presented in reviews in the field [10,20]. Methodological problems in studies on compensatory growth have however been raised when considering growth trajectories [10]. And it is noteworthy that an absence of clear compensatory growth trajectory has been previously reported in fasted/refed trout [21,22] and tilapia [23]. Changes in body shape (condition factor) was more suggestive of a compensatory response in our study, but food present de novo in the digestive tract of refed animals likely induced biased weights and condition factors. On the other hand, one cannot formally exclude that the experimental period chosen in our study was too limited to allow the full expression of the compensatory growth potential. Anyway, that a specific compensatory muscle growth response occurred in trout after refeeding was suggested by the identification in refed animals of a muscle molecular signature enriched in genes promoting muscle growth. Gene expression profiling overview showed that genes up-regulated after refeeding fell into two distinct categories: those whose expression level was restored to the values found in control-fed trout displaying normal growth and those whose expression level exceeded that found in control trout. Given that compensatory growth refers to a faster than usual growth rate [10], we reasoned that only genes in the second category formed the compensatory muscle growth signature, whereas genes of the first category signaled the restoration of cellular mechanisms involved in normal growth. The GO terms associated with genes whose expression was restored only after refeeding included cell division, organelle fission and condensed chromosome, thus suggesting that cell proliferation decreased in muscle during fasting but resumed after refeeding. Whether the resumption of cell proliferation after refeeding targets myogenic cells remains to be determined. However, a stimulation of myogenic cell proliferation has been reported in the Antartic fish Notothenia coriiceps after feeding [24]. In line with results from previous studies on muscle transcriptome dynamics during fasting-induced recovery growth [14,25], we observed, in refed trout, the up-regulation of genes encoding structural components such as sarcomeric proteins and matricial compounds. However, we found that this up-regulation was limited to a restoration of the expression level found during normal growth. Therefore, as in the case of cell cycle regulators, matricial compounds and myofibrillar proteins were excluded from the compensatory growth signature. In sharp contrast, a very large number of genes stimulating ribosome biogenesis or enhancing translational efficiency were up-regulated in compensatory muscle growth compared to normal growth. This finding strongly suggested that the compensatory growth response was associated with an accretion of the protein mass necessary for muscle fibre hypertrophy. In agreement with this finding, a correlation has been recently established between ribosome biogenesis and the magnitude of fibre hypertrophy in overloaded mouse skeletal muscle [26], and emerging evidence supports the view that ribosome biogenesis is a crucial mechanism used by skeletal muscle to regulate protein synthesis and control muscle mass [27]. Our study further showed that the capacity to convert nascent polypeptides into functional three-dimensional structures also increased, as shown by the up-regulation of a large number of genes involved in this process, notably HSP90 and HSP70. HSP90 is required for myofibril assembly in developing zebrafish embryos [28], and its expression has been shown to increase during muscle hypertrophy resulting from functional overload in rats and mice [29]. Recently, Hsp70-null mice have been reported to display a deficit in muscle fibre size [30]. Additionally, in agreement with protein synthesis and cellular growth, which require adjustments in mitochondrial ATP production, many genes involved in mitochondrial biogenesis were found in the compensatory muscle growth signature. This finding is in line with previous morphometric analyses showing an increase in mitochondrial volume density during compensatory muscle hypertrophy produced by tenotomy of the tibialis anterior muscles of rats [31]. Palstra et al. have reported in zebrafish that muscle fibre hypertrophy promoted by swimming-induced exercise is associated with an activation of the myogenic program [32]. In contrast, we found here that the up-regulation of genes promoting myofibre hypertrophy was not associated with an activation of genes involved in myofibre production. In particular, myogenic regulatory factors of the MyoD family (MyoD, myogenin, myf5 and mrf4), which constitute a cross-regulatory transcriptional network at the core of myogenesis [33], were not found in the compensatory growth response. In line with this observation, Johansen and Overturf have reported only very few changes in MyoD1b and myf5 expression in trout muscle during refeeding following starvation [34]. In addition, the two Igdomain transmembrane proteins Kin of Irre like3 and jamb which are essential for myocyte fusion in zebrafish embryo [18,19], and the genes encoding contractile proteins prevalent in hyperplastic growth zones were absent from the compensatory muscle growth signature. Finally, the compensatory muscle growth signature was found to share only a subset of the chromatin-remodelling factors evidenced in hyperplastic growth zones. In particular, the SWI/SNF chromatin-remodelling enzymes up-regulated in hyperplastic growth zones were totally absent from this signature. This finding and the demonstration that Brg1/ smarca4-SWI/SNF complexes are key epigenetic determinants of skeletal muscle differentiation [35] strengthen the notion that the steeper growth rate associated with compensatory growth does not involve an increase in myogenic cell differentiation. Contrasting to SWI/SNF components, many histone N-methyl transferases were both upregulated in hyperplastic growth zones and present in the Fig. 4 The compensatory muscle growth response involves only a subpart of the molecular signature of the hyperplastic growth zone. Venn diagram representing the distribution of genes of the compensatory muscle growth signature and genes up-regulated in the superficial hyperplastic growth zone of the late trout embryo. Functional categories inferred from genes common to the compensatory muscle growth and the hyperplastic growth zone signatures are detailed and major functional categories specific to hyperplastic growth zones are mentioned. The 312 genes specific of the compensatory muscle growth response were mostly related to translation, protein folding, RNA processing and ribosome biogenesis compensatory growth signature. Among them were members of the protein arginine M-methyltransferases family. Nine distinct protein arginine M-methyltransferases have been identified in mammals [36] and only Prmt4/carm1 and Prmt5 have been reported to exert myogenic activity [37,38]. Both Prmt4/carm1 and Prmt5 were up-regulated in hyperplastic growth zones but only Prmt5 was found in the compensatory muscle growth signature. How chromatin-remodelling factors regulate the expression of genes involved in the compensatory muscle growth response will be an interesting topic for futures studies. Together, all these data strongly suggest that the steeper growth rate associated with compensatory growth results from the stimulation of hypertrophy rather than the stimulation of hyperplasia. Of note, the compensatory muscle growth signature strongly differs from the transcriptome associated with regenerating muscle [8]; the latter shows redeployment of the transcriptional program involved in myogenesis or myofibrillogenesis of nascent myotubes [8]. This result indicates that distinct genetic pathways may be activated to drive distinct aspects of the growth and remodeling of adult fish muscle.

Conclusion
Using an Agilent-based microarray platform, we identified genes composing the molecular signature of muscle growth in fasted/refed trout and found that they were relevant to functional categories mostly involved in protein synthesis and accretion. In contrast, genes upregulated in hyperplastic growth zones of trout muscle and that are known or predicted to be important in the formation of new myofibres were virtually absent from this signature. These findings suggest that a burst of muscle growth occurs in refed trout which is mediated by stimulation of myofibre hypertrophy. The generation of a large set of genes composing the compensatory muscle growth signature expands our understanding of the molecular and cellular mechanisms controlling skeletal muscle mass, and provide a useful list of potential molecular markers of muscle growth in fish.

Methods
Animals and experimental design A spring strain of l-year-old rainbow trout (Oncorhynchus mykiss) was used.Trout were continually fed to satiation until the beginning of the experiment, at which time two groups of trout weighing approximately 130 g were constituted and reared in two separate tanks containing 2000 l of freshwater each. The experimental group was deprived of food for 30 days and refed for 36 days at a rate three times higher than satiation ration to exacerbate growth resumption. Trout of the control group were continually fed to satiation until the end of the experiment. Fasted fish and refed fish were taken from the experimental group at 0 day, 4 days, 11 days and 36 days following refeeding (termed F0, RF4, RF11 and RF36, respectively) with 10 fish sampled at each time point. Control-fed trout were taken 30 days (termed C0) and 66 days after the beginning of the experiment with 10 fish sampled at each time point. Fish were exposed to natural photoperiod and fed with a commercial diet (BioMar, Nersac, France). The fish were rapidly anaesthetized with phenoxy-ethanol (Aquaveto, 4 ml per 10 l of fresh water) before total body weighing, length measurement and muscle sampling (10 fish at each time point). The condition factor (an indicator of the body shape) was calculated as follows: K = body weightx100/ body length 3 (the body length did not include the caudal fin length). Growth trajectories and condition factors were calculated from the 10 samples taken at each time point.
RNA extraction, labelled cRNA preparation, hybridisation RNA samples of five distinct trout per feeding conditions (F0, RF4, RF11, RF36 and C0) were used for microarray experiments. A transverse slice of fast muscle situated just beneath the dorsal fin was taken for RNA extraction using TRIzol (Invitrogen, Carlsbad, CA, USA) reagent following the manufacturer's instructions. Cy3labelled cRNA generation and hybridisation were performed as previously described [15]. A control sample that did not give signal on microarray was discarded from the gene expression analysis.

Microarray slides
An Agilent-based microarray platform with 8 × 60 K probes per slide was used. This platform (GEO platform record: GPL15840) that was based on a rainbow trout resource designed by Yao and colleagues [39] was enriched with oligonucleotides designed from recent trout NGS data (http://ngspipelines-sigenae.toulouse.inra.fr:9064/). Microarray data sets have been submitted to the GEO-NCBI with the accession number: GSE91048.

Data acquisition and analysis
After hybridisation, the slides were rinsed and scanned at a 3 μm resolution with the Agilent DNA Microarray Scanner. Fluorescence intensity was calculated using the standard procedures found in the Agilent Feature Extraction (FE) software version 10.7.3.1. The arrays were normalised (scale normalisation) and log-transformed using Gene-Spring software (version 12.6.1). An ANOVA analysis (Benjamini-Hochberg (BH) corrected pval < 0.05) was used to determine the genes that were differentially expressed between feeding conditions. Gene expression