BACKGROUND
for the majority of intensively-reared finfish species including the european sea bass , diets have traditionally been based on fish meal  and fish oil .
however, the decline in worldwide supplies of marine oils and fish meal  <cit>  has led the industry and several research initiatives to investigate the possibility of using plant proteins and vegetable oils as alternatives to marine fishery-derived proteins and oils.
nevertheless, the use of such plant products is recognised to have several disadvantages, particularly related to their protein contents, amino acid profiles and unsaturated fatty acid imbalances, but also including endogenous anti-nutritional factors.
taking into account these limits and the dietary needs of different fish species, efforts have been made over the last decade to develop diets with a low content in fish resources.
this has been done by using a mixture of vegetable meals and oils  <cit> , resulting in the successful reduction of both fm and fo in the feeds for several species  <cit> .
much progress has indeed been made in the substitution of fm and fo with plant products in feeds for salmonids as well as marine fish, in the recent past  <cit> .
while several studies performed on salmonids indicate that total replacement of fish meal by plant ingredients leads to decreased growth rate  <cit> , kaushik et al.
<cit>  showed that it was possible to almost totally replace fish meal with a mixture of plant protein sources for european sea bass without reducing growth performance.
the same authors did, however, note a significant increase in fat content and a decrease in plasma cholesterol concentrations for sea bass fed with plant protein, suggesting altered regulation of lipid metabolic pathways.
for the replacement of fish oil, it is well established that freshwater or anadromous fish species such as salmonids have higher tolerance to vegetable oil compared with marine fish species.
thus, for atlantic salmon  and rainbow trout , the total replacement of fish oil with a blend of vegetable oils poor in highly unsaturated fatty acids  did not result in diminished growth performance, feed conversion or development of histopathology, despite an increase of polyunsaturated fatty acid  deposition in liver and muscle  <cit> .
in some studies, a high or total substitution of fish oil by linseed and soybean oils for several months induced decreases in growth rate of gilthead sea bream  and european sea bass   <cit> .
some other studies undertaken with gilthead sea bream showed that while there were no differences in growth of fish fed high levels of vegetable oil mixtures, there were possibly other metabolic consequences  <cit> .
this lower adaptation of marine fish species to vegetable oil has been suggested to be linked to their lower efficiency at synthesizing lc-pufa from n- <dig> and n- <dig> precursors present in plants  <cit> .
a recent study performed on european sea bass indicates that the limiting step for lc-pufa synthesis could be linked to a deficiency in the stimulation of delta-6-desaturase  activity in fish fed vegetable oil  <cit> .
the resulting low tissue levels of lc-pufa in marine fish fed vegetable oil could impact fish health, since lc-pufa are not only important as structural components of cell membranes but also as precursors of eicosanoids.
eicosanoids are involved in many physiological processes, including osmoregulation, immune responses, blood coagulation and reproduction  <cit> .
moreover, lowered eicosapentaenoic acid  and docosahexaenoic acid  content in the flesh of marine fish fed a vegetable diet diminishes their nutritional value for consumers.
recent studies on salmonids have suggested there is genetic variability for ability to utilize plant-based diets  <cit> .
interestingly, some genotype-diet interactions for growth have also been recently demonstrated in european sea bass fed on a plant-based diet  <cit> .
the existence of such interactions suggests that it could be possible to select fish, and particularly sea bass, with a better ability to grow on plant-based feeds.
however, the genetic factors and related metabolic or physiological pathways responsible for these advantageous capacities are still unknown.
to our knowledge, studies on the total replacement of both fm and fo have not been undertaken in a marine fish species until now, except for the afore mentioned work by le boucher et al.
<cit> .
moreover, investigations on the impact of fm and fo substitution with plant products for marine fish species have only been performed using molecular and/or biochemical approaches focused on selected target metabolic pathways or physiological functions.
such dedicated approaches do not allow an exhaustive and global overview of the molecular mechanisms underlying tissue and organism response to diet substitution.
in order to gain a fuller picture of the effects of total substitution of both fm and fo, the present study primarily aimed to characterise the regulation of the liver transcriptome in european sea bass fed on a fish-free diet for  <dig> months, using an oligonucleotide microarray recently developed for this species  <cit> .
this investigation was performed on liver because this organ plays a key role in intermediary metabolism, integrates a large part of nutrient uptake and affects a wide range of functions in the body, including plasma protein synthesis, hormone production and detoxification.
the present study was undertaken using two half-sibfamilies that exhibited similar growth rates when they were fed a fm-fo diet, but different growth rates when they were fed a plant protein-vegetable oil based diet.
the second aim of this work was to pinpoint genes and related metabolic and physiological pathways that could explain the different adaptation of these two half-sibfamilies of european sea bass to a plant-based diet.
the hepatic transcriptomes and flesh lc-pufa profiles were, therefore, compared between these half-sibfamilies.
methods
diets and fish
two practical iso-energetic and iso-nitrogenous diets were formulated .
the first, a fish-based diet , was composed of fish meal, wheat gluten and fish oil whereas the second, a vegetable-based diet , was devoid of ingredients of fish origin and composed of plant protein sources and vegetable oil .
the fatty acid composition of the two diets is given in table  <dig>
ingredient composition , amino acid profiles  and chemical composition  of the two diets fed to european sea bass.
1attractant mix contained: glycine , alanine , betaine , taurine  and glucosamine
<dig>  <dig> as per nrc  <cit>
all procedures concerning the animals and their handling were conducted in accordance with the code of ethics of the world medical association .
the study was performed under licence no.
<dig>  of the french department of veterinary services  to conduct experimental protocols and samplings on fish.
the present study focused on fish of two half-sibfamilies , which exhibited a similar daily growth coefficient  when they were fed on a fish-based diet, but had significantly different dgcs when they were fed an all-plant diet.
the two half-sibfamilies of fish were produced from a crossing design between  <dig> females and  <dig> males, which was carried out with european sea bass individuals from a mediterranean stock held at the experimental station of palavas-les-flots   <cit> .
rearing conditions have already been described by le boucher et al.
<cit> : fish were reared in two tanks per diet condition, supplied with recirculated seawater  at a constant temperature of 21°c, and subjected to a photoperiod of  <dig> h light: <dig> h dark.
fish were fed on a commercial diet  until they reached the mean weight of  <dig> g. before the beginning of the nutritional challenge, fish were individually tagged, genotyped for microsatellite markers to infer parentage  <cit> , and distributed randomly into two tanks per dietary condition , with  <dig> fish per tank.
after an acclimation period of  <dig> weeks, fish were hand fed to satiation  for a period of  <dig> months on the experimental diets .
analysis of experimental data was done using the following formulae:
hepatosomatic index%:hsi=100×liver weightg×body weightg-1viscerosomatic index%:vsi=100×carcass weightg×body weightg- <dig>
daily growth coefficient : dgc =  <dig> ×  1/ <dig> - initial individual weight  1/3)/days .
feed efficiency : /feed ration
at the end of the growth trial, liver, muscle and plasma were sampled from  <dig> fishes per half-sibfamily and per dietary treatment.
muscular lc-pufa profiles and real-time pcr investigations were performed on these  <dig> sampled individuals per group.
microarray analysis of hepatic rna and investigation of immune parameters in plasma were investigated on  <dig> to  <dig> sampled individuals per group in order to respect similar sex ratio, sampling time and rna quality  > 9).
lipid extraction and fatty acid analysis
one gram of white muscle was dissected from  <dig> fish per dietary treatment and immediately frozen in liquid nitrogen.
total lipids were extracted according to the method of folch et al.
<cit> , by accelerated solvent extraction  <dig>  with dichloromethane/methanol  containing  <dig> % butylatedhydroxutoluene  as antioxidant.
lipids were extracted at  <dig> bars, 100°c, with a  <dig> min precallingphase,  <dig> min static phase, and 60% flush for  <dig> sec .
the separation of neutral and polar lipids was performed according to the procedure described by juaneda and roquelin  <cit> .
the total lipid  extracts were fractionated on silica cartridges , neutral lipids  were eluted with chloroform and polar lipids  with methanol.
fatty acid methyl esters  were prepared from total lipids by acid-catalysed trans-esterification.
fame were quantified by gas-liquid chromatography  with a bpx <dig> column of  <dig> m length and  <dig>  mm i.d.
hydrogen was used as carrier gas and temperature programming was from 50°c to 180°c at 20°c/min and then to 220°c at 3°c/min.
individual methyl esters were identified by comparison with known standards.
the fatty acid analysis was performed on one sample per fish.
rna extraction and real-time quantitative pcr analysis
total mrna of liver was extracted using trizol reagent  and quantified by measuring absorbance at  <dig> nm in a spectrophotometer .
the reverse transcription was performed using the quantitect® reverse transcription kit , including a genomic dna elimination reaction.
reactions were carried out in a volume of  <dig> μl, containing  <dig> μg of total rna,  <dig> unit of quantiscript reverse transcriptase,  <dig> μl of quantiscript rt buffer  and  <dig> μm primer mix.
seven genes involved in metabolic and/or immune pathways of interest , whose the oligonucleotides were spotted on the chip, were analysed by real-time pcr in order to validate the gene expression patterns obtained through the microarray approach.
the relative mrna levels were automatically normalized with housekeeping elongation factor  <dig>  gene expression and measured by bio-rad iq <dig> software using ΔΔct method :
gene normalized expression in sample  <dig> = / with relative quantity in sample  <dig> for a gene i = e ct  - ct
ef <dig> was chosen to provide an internal control for real-time pcr, since contrary to 18s rrna and actin initially tested, we did not observe any significant difference  between ct values for ef <dig> between the dietary groups .
its stability was also assessed by a low coefficient of variation over all the samples .
specific primers  were designed from european sea bass sequences of fads <dig> , fabp <dig> , hmgcr , angptl <dig> , cxcl <dig> , gck , lpl  and ef <dig>  .
all primers used for real-time quantitative pcr analysis were defined with the primer <dig> software http://frodo.wi.mit.edu/primer3/ in order to respect an annealing temperature of 60°c.
all pcr reactions were performed with an efficiency of 100% .
the pcr reactions were carried out in an i-cycler with an optical module , in a final volume of  <dig> μl containing  <dig>  μl sybr green supermix ,  <dig>  μl of each primer  and  <dig> μl of cdna .
the pcr program consisted of an initial dna denaturation of 94°c for  <dig> s, followed by  <dig> cycles at 95°c for  <dig> s and 60°c for  <dig> s. a triplicate of the amplification reaction was realised for each sample.
plasma lysozyme concentration
plasma lysozyme activity was determined at ambient temperature using a turbidimetric assay  <cit> , adapted to microtitration plates.
briefly, a bacterial suspension of micrococcus lysodeikticus  was prepared at a concentration of  <dig>  g.l- <dig> in a  <dig>  m sodium-phosphate buffer, ph  <dig> .
fifty microlitres of the samples were placed in 96-well microtitration plates.
the reaction was initiated in a labsystems-iems analyser, by addition of  <dig> μl well-1m.
lysodeikticus suspension using an automatic dispenser.
the optic density  reading was taken at a wavelength of  <dig> nm every  <dig> s for  <dig> min, the plate being shaken before each reading.
lysozyme values for samples were converted to mg.ml- <dig>  using a reference curve established with hen egg white lysozyme .
plasma alternative complement pathway activity
determination of the alternative pathway of plasma complement activity was carried out at 4°c through a haemolytic assay with rabbitred blood cells  as described by yano  <cit> , adapted to microtitration plates.
sea bass samples, diluted to 1/ <dig> in egta-mg-gvb buffer to avoid natural haemolytic activity, were added in increasing amounts, from  <dig> to  <dig> μl well- <dig>  wells were filled with egta-mg-gvb buffer to a final volume of  <dig> μl.
finally,  <dig> μl of 2% rrc  suspension was added to each of the wells.
control values of 0% and 100% haemolysis were obtained using  <dig> μl of egta-mg-gvb buffer and  <dig> μl of non-decomplemented trout haemolytic serum at 1/ <dig> in ultrapure water, respectively.
the samples were then incubated for  <dig> h at 20°c.
the microplates were centrifuged  and  <dig> μl of supernatant from each well was then transferred into another 96-well microplate with  <dig> μl of phosphate buffered saline .
the absorbance  was read in a labsystems-iems analyser and the number of ach <dig> units per ml of plasma was determined by reference to the 50% haemolysis.
dicentrarchus labrax oligonucleotide microarray
gene expression profiles were investigated using the agilent- <dig> dicentrarchus labrax oligo microarray .
this platform represents  <dig>   <dig> unique transcripts of the european sea bass.
two non-overlapping probes were designed for each transcript for a total  <dig>   <dig> oligonucleotide probes  synthesized onto the array .
all sequences  are publicly available in a dedicated database  <cit> , together with associated annotations, go entries and putative homologous genes in fish model species.
microarrays were synthesized in situ using agilent non-contact ink-jet technology with a  <dig> ×  <dig> k format, and included default positive and negative controls.
microarray analysis was based on a single color  design.
a mixture of  <dig> different viral polyadenilated rnas  was also added to each rna sample to monitor labelling and hybridization quality as well as microarray analysis work-flow.
sample labelling and hybridization were performed as reported in ferraresso et al.
<cit> .
briefly, for each sample,  <dig> ng of total rna were linearly amplified and labelled with cy3-dctp according to the agilent one-color microarray-based gene expression analysis protocol.
a mixture of  <dig> different viral polyadenilated rnas  was also added to each rna sample to monitor labelling and hybridization quality as well as microarray analysis work-flow.
after fragmentation, a total of  <dig>   <dig> ng of labelled crna were dispensed in the gasket slide and assembled to the microarray slide .
slides were incubated for  <dig> h at 65°c in an agilent hybridization oven and washed following manufacturer's instructions.
data acquisition and normalization
hybridized slides were scanned at  <dig> μm resolution using an agilent g2565ba dna microarray scanner.
default settings were modified to scan the same slide twice at two different sensitivity levels .
the two linked images generated were analyzed together; data were extracted and background subtracted using the standard procedures contained in the agilent feature extraction  software version  <dig> .
<dig>  spike-in probe intensities were used to assess the performance of the normalization procedure for each dataset.
data normalization was performed using r statistical software http://www.r-project.org; microarray data were normalized across all arrays using the cyclic loess approach  <cit> .
fold changes  were calculated for each gene by finding the average value for each group .
raw and normalized fluorescence data of the present microarray experiment have been deposited in the geo database under accession number .
statistical analysis
all the results are presented as mean values with standard deviations .
daily growth coefficient  was studied using a model accounting for diet as a fixed effect and tank, sire, dam, sire*diet and dam*diet as random effects, using sas-glm.
effects of diet and half-sibfamily factors on biometry , fatty acid composition, gene expression , plasma lysozyme concentration and alternative complement pathway activity were tested by two-way anovas  using statistica biosoft  <dig> .
the microarray data were analysed by two-way anova using tmev  statistical software, and gene expression was considered significantly different when p <  <dig> .
significant enrichment of go biological process categories were tested for using ease software  with p <  <dig>   <cit> .
RESULTS
growth and biometry
after  <dig> months of the feeding trial, european sea bass fed vd exhibited significantly lower dgc than those fd .
in addition, the fish of half-sibfamily g fed the vd had a significantly higher dgc than fish of half-sibfamily g fed vd, while there was no difference between these two half-sibfamilies when they were both fed fd .
the hepatosomatic index  was regulated by diet and genetic factors while the viscerosomatic index  was only regulated by the genetic factor.
there were no significant interactions between dietary and genetic factors for these two parameters.
the feed efficiency  in the duplicate tanks was  <dig>  and  <dig>  for the fd and  <dig>  and  <dig>  for the vd, respectively.
mortality was not significantly different between dietary treatment and half-sibfamilies .
initial and final weights and lengths , hsi  and vsi .
effects of diet factor and genetic factor on the biometric parameters are determined by two-way anova.
results are expressed as mean +/- s.d.
and significant differences are indicated by the p value .
lipid and fatty acid compositions of the fillet
the flesh lipid composition was significantly affected by dietary treatment .
feeding vd significantly increased the percentage of saturated lipids in both the neutral lipids and phospholipids.
the α-linolenic acid  and linoleic acid  contents were respectively 10-fold and 3-fold higher, in both lipid classes, when fish were fed vd.
in addition, aa, epa and dha contents were around  <dig> -fold lower in flesh of fish fed vd.
the Σn- <dig> pufa/Σn- <dig> pufa ratio decreased in both the neutral lipid  and phospholipid  fractions in the flesh of european sea bass fed vd.
microarray gene expression profiling
a list of  <dig>   <dig> significant probes was obtained for the effect of diet factor, corresponding to  <dig> unique transcripts with gene ontology  annotation.
among these regulated transcripts,  <dig> exhibited higher levels in fish fed vd while  <dig> were over-expressed in the liver of fish fed fd.
for the family factor effect, total of  <dig> significant probes were revealed corresponding to  <dig> unique transcripts with go annotation.
among these,  <dig> exhibited higher levels in half-sibfamily g while  <dig> were more abundant in half-sibfamily g. in addition,  <dig> probes related to  <dig> genes with functional annotation exhibited significant diet × genetic interactions.
the main biological processes enriched out of those associated with genes that were over-expressed in fish fed vd were physiological process , metabolism , rna splicing , protein catabolism , aerobic respiration , sterol metabolism , carboxylic metabolism , amino acid metabolism , blood coagulation  and hexose catabolism  .
the genes involved in carboxylic acid and sterol metabolism included fatty acid desaturase  <dig>  steroyl-coa desaturase  <dig>  nadh-cytochrome b <dig> reductase, isopentenyl-diphosphate delta-isomerase  <dig>  lanosterol 14-alpha demethylase, farnesyl pyrophosphate synthase, c- <dig> methylsterol oxidase and 3-hydroxy-3-methylglutaryl-coenzyme a reductase.
apolipoprotein a-i, apolipoprotein b- <dig> and lipoprotein lipase, implicated in lipid transport, were also found to be up-expressed in fish fed vd .
similarly, some genes involved in carbohydrate metabolism , such as glucose-6-phosphate 1-dehydrogenase, 6-phosphogluconate dehydrogenase and fructose- <dig>  6-bisphosphatase  <dig> were also expressed at a higher level in the fish fed vd .
expression levels of genes involved in protein metabolism and amino acid metabolism were also increased in fish fed vd .
the fold-changes  are indicated considering fd as the reference group.
.
in contrast, the main biological processes enriched associated with the genes that were lower-expressed in fish fed vd were particularly related to cellular process , cell communication  and cell proliferation  .
in addition, some genes involved in the immune function were less expressed in fish fed vd .
the comparison of hepatic transcriptomes between the two half-sibfamilies indicated that genes involved in immune function, such as complement component c <dig>  c <dig> and c <dig> and mannan-binding lectin serine protease  <dig>  were more greatly expressed in half-sibfamily g than in half-sibfamily g .
inversely, expression levels of nadh dehydrogenase genes , involved in the electron transport to the respiratory chain, were significantly lower in half-sibfamily g compared with half-sibfamily g. half-sibfamily g was also characterised by lower expression level of genes implicated in atp production  and protein synthesis, such as ribosomal subunits  than the half-sibfamily g .
among the  <dig> genes exhibiting an interaction between half-sibfamily and diet factors,  <dig> were involved in metabolism.
however, only the processes related to aromatic amino acid family  and nucleotide metabolism  were found to be over-represented among these genes .
the fold-changes  are indicated considering the half-sibfamily g as the reference group.
.
in order to validate the accuracy of the microarray data, the fads <dig>  hmgcr, fabp <dig>  angptl <dig>  cxcl <dig>  gck and lpl genes, which were spotted on the microarray, were also investigated by means of real-time pcr.
the comparison of the gene expression pattern obtained through the real-time pcr and microarray approaches, revealed a correlation greater than  <dig>  .
immune parameters
lysozyme activity was significantly lower  in fish fed vd than in fish fed fd, while the alternative complement activity was  <dig> -times higher  .
there was no effect of the half-sibfamily factor on these activities.
discussion
the present work is the first investigation into the effect of an exclusively vegetable diet on the hepatic transcriptome in a marine fish species.
it is also the first study to have explored the transcriptome of two half-sibfamilies of european sea bass exhibiting different capacities to grow on such a diet.
the replacement of fm and fo with increasing levels of plant protein and oil sources for marine fish species can modify feed intake and conversion, which should be the major reason for associated growth delay  <cit> .
in the present study, there was a tendency for higher fe in fish fed fd  compared with fish fed vd .
however, this difference could not be statistically tested since fish were reared in only two tanks per diet condition.
a vegetable diet is also known to potentially impact fish metabolism through regulation of gene expression, especially in the liver  <cit> .
analysis of the oligo-dna microarray data by two-way anova indicated that several hundred genes were differentially regulated according to diet or/and half-sibfamily factors.
the accuracy of the present microarray data is validated by the similar gene pattern expression obtained from different oligonucleotides representing the same genes spotted on the array , as well as by the correlation shown between results of microarray and qpcr approaches.
lc-pufa metabolism
metabolism-related biological processes constitute the largest group among the go terms associated with the genes regulated in vd-fed fish.
among these, genes involved in lipid metabolism and, particularly, in lc-pufa biosynthesis, were found to be up-regulated in fish fed vd.
not surprisingly, delta-6-desaturase , steroyl-coa desaturase  <dig>  and nadh-cytochrome b <dig> reductase , involved in long chain fatty acid desaturation and/or elongation, were up-regulated in vd-fed fish.
it is indeed established that in most fish species that the lc-pufa biosynthetic pathway is positively regulated in response to the use of a diet poor in lc-pufa but rich in pufa, although this regulation depends on fish species, degree of fish oil substitution, nature of the vegetable oil, and environmental parameters   <cit> .
the positive regulation of the lc-pufa biosynthetic pathway is in agreement with results obtained by transcriptomic approaches in salmonids fed vegetable oil  <cit> .
altogether, the expression data obtained in marine fish species and salmonids indicate that the different capacity of marine and fresh water species to grow on a lc-pufa-deprived diet does not seem to be due to a different transcriptional regulation of key genes involved in lipid synthesis, such as fads <dig> or scd <dig>  indeed, the level of induction of fads <dig> expression found in this experiment  is of similar amplitude to that observed in the liver and intestine of atlantic salmon  <cit> .
the stimulation of this biosynthetic pathway in fish fed a diet poor in lc-pufas can be explained by the fact that lc-pufas play several key physiological roles in vertebrates, particularly in fish.
for example, fish are poikilothermic organisms and therefore need a high degree of unsaturation of lc-pufa included in membrane phospholipids to maintain phospholipid bilayer fluidity at reduced temperature  <cit> .
lc-pufa, especially arachidonic acid  and eicosapentaenoic acid , are also precursors of eicosanoids, which are involved in pro and anti-inflammatory pathways.
this hypothesis is reinforced by our data, indicating a stimulation of genes involved in phospholipids biosynthesis  when fish were fed the vd.
despite this stimulation of lc-pufa and the phospholipid biosynthesis pathway at the transcriptional level, our investigation of fatty acid profiles indicated that the amounts of lc-pufa, particularly eicosapentaenoic acid  and docosahexaenoic acid , were still considerably lower in the flesh of fish from both half-sibfamilies fed vd in comparison with fish fed fd.
this finding is in agreement with those previously obtained  <cit> , which revealed that the stimulation of fads <dig> expression in fish fed a vegetable diet was not associated with an induction of its enzymatic activity, suggesting a post-transcriptional regulation of fads <dig> expression.
such epa and dha deficiency can notably explain the growth deficiency observed in fish fed vd, as well as effects observed on immune function .
concerning the genetic aspect, the comparable overall expression pattern of genes involved in lc-pufa synthesis in the liver, associated with similar lc-pufa profile in muscle, in both half-sibfamilies suggests that the differing capacities of these european sea bass half-sibfamilies to grow on a vegetable diet are not due to differing capacities to synthesize lc-pufa.
sterol metabolism
the present microarray data also indicate an increase in expression levels of genes involved in sterol metabolism in vd-fed fish.
among these, isopentenyl-diphosphate delta-isomerase  <dig> , lanosterol 14-alpha demethylase , farnesyl pyrophosphate synthase , c- <dig> methylsterol oxidase  and 3-hydroxy-3-methylglutaryl-coenzyme a reductase  genes  <cit>  are known to be implicated in the cholesterol metabolic pathway .
more particularly, hmgcr, a transmembrane glycoprotein involved in the rate-limiting step of sterol biosynthesis, is increased, as shown in mammals  <cit> .
the stimulation of cholesterol biosynthesis in fish fed vd could be related to the difference in sterol composition between diets.
indeed, while the fish diet is rich in cholesterol, the vegetable diet used in this experiment contains exclusively plant sterols, which have been shown to affect membrane properties by decreasing permeability and fluidity, and modifying phospholipid order in mammals  <cit> .
as a consequence, the increase in cholesterol biosynthesis could be a metabolic response to its deficiency in the diet, as well as a way to restore membrane properties by incorporation of endogenous cholesterol.
since we did not measure the cholesterol content in the liver, flesh or blood, it is not possible for us to assess the capacity of european sea bass to compensate for possible dietary deficiencies in cholesterol through a regulation of its biosynthesis.
moreover, similarly to the lc-pufa pathway, no significant difference of cholesterol biosynthetic regulation was observed between the half-sibfamilies.
interestingly, several vd-stimulated genes  involved in the lipogenic pathway  are known to be molecular targets of sterol regulatory element binding proteins , which are key regulators of fatty acid and cholesterol synthesis  <cit> .
recent data indicating an up-regulation of the srebp- <dig> gene expression in european sea bass fed a vegetable diet  <cit>  could thus be due to such stimulations.
lipid and sterol transport
the present microarray data indicate that the stimulation of genes involved in fatty acid and cholesterol synthesis in vd-fed fish was associated with an over-expression of genes involved in their transport, such as apolipoproteins apoa <dig> and apob <dig>  which are the major protein constituents of high and low density lipoprotein , respectively .
the ldl, including apob <dig>  are involved in the transport of cholesterol and lipids from the liver to other tissues.
thus, up-regulation of apob <dig> combined with the induction of the expression of lipoprotein lipase , a key enzyme involved in the hydrolysis of triglyceride, suggests an increase in lipid transport and metabolism from the liver to tissues in fish fed vd.
the decrease in angiopoietin-related protein  <dig>  that we observed in fish fed vd reinforces this idea since angptl <dig> suppresses ldl clearance via the inhibition of lpl activity  <cit> .
in parallel, we observed that the reverse transport of cholesterol not used by tissues via hdl to the liver was also stimulated in fish fed vd.
indeed, apoa <dig>  which participates in the transport of cholesterol to the liver by promoting cholesterol efflux from tissues and by acting as a cofactor for the lecithin cholesterol acyltransferase , exhibited higher transcript levels in fish fed vd.
altogether these results reveal that another major response in the liver of european sea bass fed a vegetable diet is the stimulation of cholesterol synthesis and transport, irrespective of the half-sibfamily considered.
carbohydrate metabolism
lc-pufa and cholesterol biosynthesis require reducing power in the form of nadph.
it is well documented in vertebrates, including fish, that nadph required for malonyl-coa synthesis is mainly supplied by the dehydrogenases of the pentose phosphate shunt  <cit> .
interestingly, our transcriptomic data indicate that the use of the vd induced a significant increase in the level of glucose-6-phosphate dehydrogenase .
g6pd catalyses nadp+-linked oxidation of d-glucose-6-phosphate and has been shown to be a major contributor of nadph production for lipogenesis in atlantic salmon   <cit>  and european sea bass  <cit> .
moreover, our data indicate an increase in the expression of hexose-6-phosphate dehydrogenase  and phosphogluconate dehydrogenase , enzymes of the pentose phosphate pathway that generate nadph, in fish fed vd  <cit>  .
once synthesized, the resulting pentose sugar intermediate generated by the pentose phosphate pathway can be reconverted to intermediates of the glycolysis/gluconeogenesis pathway such as glyceraldehyde 3p or fructose 6p.
in the liver of fish, it is known that glycolysis provides essential precursors for biosynthesis rather than pyruvate for oxidation  <cit> .
thus, the stimulation of fructose- <dig>  6-bisphosphatase  <dig>  and aldolase  expression that we observe in fish fed vd could provide high levels of fructose-6-phosphate from glyceraldehyde 3p, then glucose 6p that serves as substrate for repeated passage in the pentose phosphate shunt .
protein/amino-acid metabolism and atp synthesis
our data revealed over-expression of genes involved in proteolysis and, more particularly, in proteasome activity  and ubiquitin activity  in fish fed vd , which is in total agreement with proteomic data obtained in rainbow trout, indicating a stimulation of proteolysis in fish fed vegetable diets  <cit> .
in our study, the stimulation of proteolysis in the fish fed the vegetable diet was associated with the induction of  <dig> genes involved in amino acid metabolism and, more importantly,  <dig> genes involved in glutamine metabolism.
in addition, gmps, aadat, got <dig> and tat genes, which are implicated in transamination, were also stimulated in fish fed vd.
the processes related to amino acid metabolism and, especially transamination, are important steps in the synthesis of some non-essential amino acids such as α-ketoglutarate.
for example, the synthesis of α-ketoglutarate through transamination reactions could be used in the tca cycle to provide energy.
interestingly, we found over-expression of genes encoding enzymes involved in the tca cycle, such as succinate dehydrogenase , isocitrate dehydrogenase  and malate dehydrogenase , in fish fed vd.
this stimulation of the tca cycle could be related to the higher levels of atp required for lc-pufa and cholesterol biosynthesis in fish fed vd.
since marine fish have a low capacity to digest complex carbohydrates, in contrast to mammals  <cit> , the use of proteins as an essential source of energy can thus explain the stimulation of the amino acid metabolism in fish fed vd.
as shown in rainbow trout fed on a vegetable-based diet  <cit> , the lower growth rate in fish fed vd in the present study could be associated with higher proteolytic activity compared with fish fed fd.
interestingly, while both half-sibfamilies g and g exhibited similar proteolysis regulation, the expression of several genes involved in macromolecule biosynthesis, and particularly in protein biosynthesis , were up-regulated in half-sibfamily g. this result, suggesting a higher protein turnover in half-sibfamily g compared with half-sibfamily g when fish were fed vd, could be related to the higher growth rate observed in half-sibfamily g fed vd.
as protein biosynthesis requires energy from atp hydrolysis, the higher protein biosynthesis in half-sibfamily could be related to a higher activity of mitochondrial atp production.
accordingly, genes involved in atp biosynthesis  and atp synthesis-coupled electron transport  were found up-regulated in half-sibfamily g. the diet × half-sibfamily interaction that we found for the expression of genes involved in aromatic amino acid metabolism reinforces the difference in protein metabolism between the two half-sibfamilies.
immune function
in the present study, the vegetable diet used, based on linseed oil, was characterised by a very low ara content and poor levels of n- <dig> lc-pufa.
eicosanoids derived from ara are known to be involved in the proliferation of hepatocytes and immune cells  <cit> .
as a consequence, the lower hepatosomatic index  measured in fish fed vd could be linked to a lower hepatocyte proliferation due to a deficiency in ara, which was supported by down-expression of  <dig> genes involved in cell proliferation in this dietary group.
moreover, fatty acid imbalances can induce an immune deficiency in all vertebrates including fish  <cit> .
in particular, the Σ n-3/Σ n- <dig> fatty acid ratio is considered as a key element regulating immune cell structure, cell signalling, and eicosanoid production.
the present microarray data revealed genes of the immune system, particularly the innate immune response, exhibiting lower expression in fish fed with vd.
interleukin  <dig>  and c-x-c motif chemokine  <dig> , which are chemotactic factors for granulocytes and monocytes, respectively, were found to be less expressed in fish fed vd, suggesting a down regulation of the innate defence system and the pro-inflammatory pathway in this group.
deficiency of the inflammatory response was also in agreement with the higher levels of transcripts of fatty acid binding protein  <dig> , whose expression in mammals has been shown to be restricted to the kupffer cells  <cit> , and the down expression of the c-reactive protein , an acute phase protein synthesised by hepatocytes, in fish fed vd.
decrease in inflammatory response can also be related to the low level of ara in the fish fed vd, which induces a reduction of prostaglandin synthesis derived from this fatty acid.
our microarray data indeed show that prostaglandin e synthase  <dig> , involved in the synthesis of pro-inflammatory prostaglandin e <dig>  is down-regulated in fish fed vd, while prostaglandin e synthase  <dig> , which has anti-inflammatory properties, exhibited higher messenger levels in fish fed vd.
this depression of innate immune system, particularly pro-inflammatory activity, could also be partially explained by a defect in membrane properties in fish fed vd, as revealed by the down-regulation of a large number of genes  related to cell communication, including factors such as cytokine receptor common subunit gamma , receptor-type tyrosine-protein phosphatase f  or integrin beta  <dig> , which are cell-surface receptor binding proteins and/or cell adhesion receptors involved in immune response.
the depression of the innate immune response in fish fed vd was confirmed by the lower plasmatic lysozyme concentration and lower-expression of lysozyme g  gene.
surprisingly, the alternative complement pathway activity involved in the innate immune response, which we assessed by analysis of plasma parameters, showed a significantly higher level in fish fed vd.
such an opposite regulation of the immune pathway revealed that different components of the immune systems can be regulated in opposite directions.
interestingly, processes related to the humoral immune response were also over-represented among the genes up-regulated in half-sibfamily g. indeed, complement component c <dig>  c <dig> and c <dig> genes showed higher expression levels in half-sibfamily g. however, the up-regulation of genes involved in the alternative complement pathway cannot be associated with an increase of the plasma alternative complement pathway activity, probably due to the complexity of factors and regulation levels  involved in the regulation of this pathway  <cit> .
moreover, the higher expression of masp <dig>  tnrfrf <dig>  c <dig> and c <dig> genes involved in the inflammatory response might reflect higher inflammatory states in half-sibfamily g, which could be associated with a decrease in growth rate, as demonstrated in chicken  <cit> .
blood coagulation
blood coagulation is another process involved in the innate immune system.
lc-pufa and, more specifically, epa, dha and ara are precursors for eicosanoid synthesis involved in the control of the blood coagulation  <cit> .
as mentioned above, the use of a diet composed of vegetable protein and oil induces modifications in the membrane phospholipid composition, with possible consequences for eicosanoid production and the blood coagulation process.
we found that the use of the vegetable diet induced an increase in the expression of genes involved in the blood coagulation pathway .
among these genes, prothrombin , coagulation factor × , fibrinogen beta chain , fibrinogen gamma chain  and coagulation factor vii  were positively involved in the blood coagulation process.
on the basis of these results, the use of a vd seems to cause pro-coagulant action by the stimulation of the blood coagulation pathway, which is in agreement with our visual observation of plasma clotting .
in agreement with these results, tavares-dias  <cit>  showed that, contrary to a vegetable diet, dietary enrichment in long-chain n- <dig> fatty acids has a strong hypocoagulant action.
in addition, the regulation of the coagulation pathway is complex and under the control of several negative factors that maintain a physiological homeostasis.
the pro-coagulation effect in response to a vegetable diet is notably reinforced by some of these genes , which also exhibited higher expression in fish fed vd .
CONCLUSIONS
the nutrigenomic approach used here has revealed several new genes and related biological processes regulated by a vegetable diet.
in particular, genes involved in lipid metabolism, protein/amino acid metabolism, carbohydrate metabolism, immune function, blood coagulation and the rna splicing process were expressed at a higher level in fish fed with vd.
the comparison of transcriptomic response in two half-sibfamilies of fish exhibiting different growth rates when fed the vegetable diet also revealed some biological processes related to protein turnover and immune response, potentially due to better adaptation to this diet.
finally, in the context of developing novel diets for aquaculture and selecting fish families exhibiting higher adaptation to fish oil and meal substitution, this work enabled us to pinpoint potentially useful new molecular markers for identifying the physiological effects of a vegetable diet, as well as a family exhibiting a phenotype of interest.
list of abbreviations
ala: α-linolenic acid; ara: arachidonic acid; dgc: daily growth coefficient; dha: docosahexaenoic acid; dm: dry matter; ef1: elongation factor 1; epa: eicosapentaenoic acid; fd: fish diet; fe: feed efficiency; go: gene ontology; hdl: high density lipoprotein; hsi: hepatosomaticindex; la: linoleic acid; lc-pufa: long chain polyunsaturated fatty acid; ldl: low density lipoprotein; nl: neutral lipid; pl: polar lipid; pufa: polyunsaturated fatty acid; tl: total lipid; vd: vegetable diet; vsi: viscerosomatic index.
authors' contributions
fg and dm co-ordinated the sampling procedures, biochemical studies and microarray data analysis, and drafted the manuscript.
sf and lb carried out the sample preparation, hybridizations and microarray data normalizations.
sk formulated the experimental diet.
mv participated in the design of the study and performed analysis of growth parameters.
cc and jlz conceived and designed the study and led its coordination.
cq co-ordinated dedicated experiments on immune function.
all authors read and approved the final manuscript.
supplementary material
additional file 1
comparison of ct values for ef <dig> gene between the dietary groups.
quantitative pcr did not reveal any significant difference  between ct values for ef <dig> between the dietary groups .
click here for file
additional file 2
fatty acid composition in muscle of two european sea bass half-sibfamilies fed fd or vd.
composition in terms of the main fatty acids  in neutral lipid and phospholipid fractions in muscle of each of the half-sibfamilies  of european sea bass fed fd or vd.
effects of diet factor and half-sibfamily factor on fatty acid composition were determined by two-way anova.
results are expressed as mean +/- s.d.
and significant differences are indicated by the p value .
click here for file
additional file 3
significantly enriched biological processes associated with genes regulated by diet and half-sibfamily factors .
the main biological processes enriched out of those associated with genes that were over-expressed in fish fed vd were related to physiological process, metabolism , rna splicing, protein catabolism, aerobic respiration, blood coagulation and hexose catabolism.
in contrast, the main biological processes associated with the genes lower-expressed in fish fed vd were related to cellular process, cell communication and cell proliferation.
regarding half-sibfamily factors, biological process related to humoral immune response was shown to be over-represented within genes up-expressed in half-sibfamily g while processes related to energy pathways  were enriched within genes up-expressed in half-sibfamily g.
click here for file
acknowledgements
this work was supported by marine genomics europe, europôle mer axis 1: "genomics and 'blue chemistry' ".
the authors thank h. le delliou for fatty acid analysis, f. ruelle for european sea bass breeding, f. terrier for making the experimental diets, f. lamour, c. huelvan, e. desbruyères for the measurement of immune parameters on plasma, dr. a. henaut for discussion on microarrays analysis and h. mccombie for english revisions.

 for the majority of intensively-reared finfish species including the european sea bass, diets have traditionally been based on fish meal and fish oil. however, the decline in worldwide supplies of marine oils and fish meal has led the industry and several research initiatives to investigate the possibility of using plant proteins and vegetable oils as alternatives to marine fishery-derived proteins and oils. nevertheless, the use of such plant products is recognised to have several disadvantages, particularly related to their protein contents, amino acid profiles and unsaturated fatty acid imbalances, but also including endogenous anti-nutritional factors. taking into account these limits and the dietary needs of different fish species, efforts have been made over the last decade to develop diets with a low content in fish resources. this has been done by using a mixture of vegetable meals and oils, resulting in the successful reduction of both fm and fo in the feeds for several species. much progress has indeed been made in the substitution of fm and fo with plant products in feeds for salmonids as well as marine fish, in the recent past. while several studies performed on salmonids indicate that total replacement of fish meal by plant ingredients leads to decreased growth rate, kaushik et al. showed that it was possible to almost totally replace fish meal with a mixture of plant protein sources for european sea bass without reducing growth performance. the same authors did, however, note a significant increase in fat content and a decrease in plasma cholesterol concentrations for sea bass fed with plant protein, suggesting altered regulation of lipid metabolic pathways. for the replacement of fish oil, it is well established that freshwater or anadromous fish species such as salmonids have higher tolerance to vegetable oil compared with marine fish species. thus, for atlantic salmon and rainbow trout, the total replacement of fish oil with a blend of vegetable oils poor in highly unsaturated fatty acids did not result in diminished growth performance, feed conversion or development of histopathology, despite an increase of polyunsaturated fatty acid deposition in liver and muscle. in some studies, a high or total substitution of fish oil by linseed and soybean oils for several months induced decreases in growth rate of gilthead sea bream and european sea bass. some other studies undertaken with gilthead sea bream showed that while there were no differences in growth of fish fed high levels of vegetable oil mixtures, there were possibly other metabolic consequences. this lower adaptation of marine fish species to vegetable oil has been suggested to be linked to their lower efficiency at synthesizing lc-pufa from n- and n- precursors present in plants. a recent study performed on european sea bass indicates that the limiting step for lc-pufa synthesis could be linked to a deficiency in the stimulation of delta-6-desaturase activity in fish fed vegetable oil. the resulting low tissue levels of lc-pufa in marine fish fed vegetable oil could impact fish health, since lc-pufa are not only important as structural components of cell membranes but also as precursors of eicosanoids. eicosanoids are involved in many physiological processes, including osmoregulation, immune responses, blood coagulation and reproduction. moreover, lowered eicosapentaenoic acid and docosahexaenoic acid content in the flesh of marine fish fed a vegetable diet diminishes their nutritional value for consumers. recent studies on salmonids have suggested there is genetic variability for ability to utilize plant-based diets. interestingly, some genotype-diet interactions for growth have also been recently demonstrated in european sea bass fed on a plant-based diet. the existence of such interactions suggests that it could be possible to select fish, and particularly sea bass, with a better ability to grow on plant-based feeds. however, the genetic factors and related metabolic or physiological pathways responsible for these advantageous capacities are still unknown. to our knowledge, studies on the total replacement of both fm and fo have not been undertaken in a marine fish species until now, except for the afore mentioned work by le boucher et al.. moreover, investigations on the impact of fm and fo substitution with plant products for marine fish species have only been performed using molecular and/or biochemical approaches focused on selected target metabolic pathways or physiological functions. such dedicated approaches do not allow an exhaustive and global overview of the molecular mechanisms underlying tissue and organism response to diet substitution. in order to gain a fuller picture of the effects of total substitution of both fm and fo, the present study primarily aimed to characterise the regulation of the liver transcriptome in european sea bass fed on a fish-free diet for months, using an oligonucleotide microarray recently developed for this species. this investigation was performed on liver because this organ plays a key role in intermediary metabolism, integrates a large part of nutrient uptake and affects a wide range of functions in the body, including plasma protein synthesis, hormone production and detoxification. the present study was undertaken using two half-sibfamilies that exhibited similar growth rates when they were fed a fm-fo diet, but different growth rates when they were fed a plant protein-vegetable oil based diet. the second aim of this work was to pinpoint genes and related metabolic and physiological pathways that could explain the different adaptation of these two half-sibfamilies of european sea bass to a plant-based diet. the hepatic transcriptomes and flesh lc-pufa profiles were, therefore, compared between these half-sibfamilies. methods diets and fish two practical iso-energetic and iso-nitrogenous diets were formulated. the first, a fish-based diet, was composed of fish meal, wheat gluten and fish oil whereas the second, a vegetable-based diet, was devoid of ingredients of fish origin and composed of plant protein sources and vegetable oil. the fatty acid composition of the two diets is given in table ingredient composition, amino acid profiles and chemical composition of the two diets fed to european sea bass. 1attractant mix contained: glycine, alanine, betaine, taurine and glucosamine as per nrc all procedures concerning the animals and their handling were conducted in accordance with the code of ethics of the world medical association. the study was performed under licence no. of the french department of veterinary services to conduct experimental protocols and samplings on fish. the present study focused on fish of two half-sibfamilies, which exhibited a similar daily growth coefficient when they were fed on a fish-based diet, but had significantly different dgcs when they were fed an all-plant diet. the two half-sibfamilies of fish were produced from a crossing design between females and males, which was carried out with european sea bass individuals from a mediterranean stock held at the experimental station of palavas-les-flots. rearing conditions have already been described by le boucher et al.: fish were reared in two tanks per diet condition, supplied with recirculated seawater at a constant temperature of 21°c, and subjected to a photoperiod of h light: h dark. fish were fed on a commercial diet until they reached the mean weight of g. before the beginning of the nutritional challenge, fish were individually tagged, genotyped for microsatellite markers to infer parentage, and distributed randomly into two tanks per dietary condition, with fish per tank. after an acclimation period of weeks, fish were hand fed to satiation for a period of months on the experimental diets. analysis of experimental data was done using the following formulae: hepatosomatic index%:hsi=100×liver weightg×body weightg-1viscerosomatic index%:vsi=100×carcass weightg×body weightg- daily growth coefficient: dgc = × 1/ - initial individual weight 1/3)/days. feed efficiency: /feed ration at the end of the growth trial, liver, muscle and plasma were sampled from fishes per half-sibfamily and per dietary treatment. muscular lc-pufa profiles and real-time pcr investigations were performed on these sampled individuals per group. microarray analysis of hepatic rna and investigation of immune parameters in plasma were investigated on to sampled individuals per group in order to respect similar sex ratio, sampling time and rna quality > 9). lipid extraction and fatty acid analysis one gram of white muscle was dissected from fish per dietary treatment and immediately frozen in liquid nitrogen. total lipids were extracted according to the method of folch et al., by accelerated solvent extraction with dichloromethane/methanol containing % butylatedhydroxutoluene as antioxidant. lipids were extracted at bars, 100°c, with a min precallingphase, min static phase, and 60% flush for sec. the separation of neutral and polar lipids was performed according to the procedure described by juaneda and roquelin. the total lipid extracts were fractionated on silica cartridges, neutral lipids were eluted with chloroform and polar lipids with methanol. fatty acid methyl esters were prepared from total lipids by acid-catalysed trans-esterification. fame were quantified by gas-liquid chromatography with a bpx column of m length and mm i.d. hydrogen was used as carrier gas and temperature programming was from 50°c to 180°c at 20°c/min and then to 220°c at 3°c/min. individual methyl esters were identified by comparison with known standards. the fatty acid analysis was performed on one sample per fish. rna extraction and real-time quantitative pcr analysis total mrna of liver was extracted using trizol reagent and quantified by measuring absorbance at nm in a spectrophotometer. the reverse transcription was performed using the quantitect® reverse transcription kit, including a genomic dna elimination reaction. reactions were carried out in a volume of μl, containing μg of total rna, unit of quantiscript reverse transcriptase, μl of quantiscript rt buffer and μm primer mix. seven genes involved in metabolic and/or immune pathways of interest, whose the oligonucleotides were spotted on the chip, were analysed by real-time pcr in order to validate the gene expression patterns obtained through the microarray approach. the relative mrna levels were automatically normalized with housekeeping elongation factor gene expression and measured by bio-rad iq software using ΔΔct method: gene normalized expression in sample = / with relative quantity in sample for a gene i = e ct - ct ef was chosen to provide an internal control for real-time pcr, since contrary to 18s rrna and actin initially tested, we did not observe any significant difference between ct values for ef between the dietary groups. its stability was also assessed by a low coefficient of variation over all the samples. specific primers were designed from european sea bass sequences of fads, fabp, hmgcr, angptl, cxcl, gck, lpl and ef. all primers used for real-time quantitative pcr analysis were defined with the primer software http://frodo.wi.mit.edu/primer3/ in order to respect an annealing temperature of 60°c. all pcr reactions were performed with an efficiency of 100%. the pcr reactions were carried out in an i-cycler with an optical module, in a final volume of μl containing μl sybr green supermix, μl of each primer and μl of cdna. the pcr program consisted of an initial dna denaturation of 94°c for s, followed by cycles at 95°c for s and 60°c for s. a triplicate of the amplification reaction was realised for each sample. plasma lysozyme concentration plasma lysozyme activity was determined at ambient temperature using a turbidimetric assay, adapted to microtitration plates. briefly, a bacterial suspension of micrococcus lysodeikticus was prepared at a concentration of g.l- in a m sodium-phosphate buffer, ph. fifty microlitres of the samples were placed in 96-well microtitration plates. the reaction was initiated in a labsystems-iems analyser, by addition of μl well-1m. lysodeikticus suspension using an automatic dispenser. the optic density reading was taken at a wavelength of nm every s for min, the plate being shaken before each reading. lysozyme values for samples were converted to mg.ml- using a reference curve established with hen egg white lysozyme. plasma alternative complement pathway activity determination of the alternative pathway of plasma complement activity was carried out at 4°c through a haemolytic assay with rabbitred blood cells as described by yano, adapted to microtitration plates. sea bass samples, diluted to 1/ in egta-mg-gvb buffer to avoid natural haemolytic activity, were added in increasing amounts, from to μl well- wells were filled with egta-mg-gvb buffer to a final volume of μl. finally, μl of 2% rrc suspension was added to each of the wells. control values of 0% and 100% haemolysis were obtained using μl of egta-mg-gvb buffer and μl of non-decomplemented trout haemolytic serum at 1/ in ultrapure water, respectively. the samples were then incubated for h at 20°c. the microplates were centrifuged and μl of supernatant from each well was then transferred into another 96-well microplate with μl of phosphate buffered saline. the absorbance was read in a labsystems-iems analyser and the number of ach units per ml of plasma was determined by reference to the 50% haemolysis. dicentrarchus labrax oligonucleotide microarray gene expression profiles were investigated using the agilent- dicentrarchus labrax oligo microarray. this platform represents unique transcripts of the european sea bass. two non-overlapping probes were designed for each transcript for a total oligonucleotide probes synthesized onto the array. all sequences are publicly available in a dedicated database, together with associated annotations, go entries and putative homologous genes in fish model species. microarrays were synthesized in situ using agilent non-contact ink-jet technology with a × k format, and included default positive and negative controls. microarray analysis was based on a single color design. a mixture of different viral polyadenilated rnas was also added to each rna sample to monitor labelling and hybridization quality as well as microarray analysis work-flow. sample labelling and hybridization were performed as reported in ferraresso et al.. briefly, for each sample, ng of total rna were linearly amplified and labelled with cy3-dctp according to the agilent one-color microarray-based gene expression analysis protocol. a mixture of different viral polyadenilated rnas was also added to each rna sample to monitor labelling and hybridization quality as well as microarray analysis work-flow. after fragmentation, a total of ng of labelled crna were dispensed in the gasket slide and assembled to the microarray slide. slides were incubated for h at 65°c in an agilent hybridization oven and washed following manufacturer's instructions. data acquisition and normalization hybridized slides were scanned at μm resolution using an agilent g2565ba dna microarray scanner. default settings were modified to scan the same slide twice at two different sensitivity levels. the two linked images generated were analyzed together; data were extracted and background subtracted using the standard procedures contained in the agilent feature extraction software version. spike-in probe intensities were used to assess the performance of the normalization procedure for each dataset. data normalization was performed using r statistical software http://www.r-project.org; microarray data were normalized across all arrays using the cyclic loess approach. fold changes were calculated for each gene by finding the average value for each group. raw and normalized fluorescence data of the present microarray experiment have been deposited in the geo database under accession number. statistical analysis all the results are presented as mean values with standard deviations. daily growth coefficient was studied using a model accounting for diet as a fixed effect and tank, sire, dam, sire*diet and dam*diet as random effects, using sas-glm. effects of diet and half-sibfamily factors on biometry, fatty acid composition, gene expression, plasma lysozyme concentration and alternative complement pathway activity were tested by two-way anovas using statistica biosoft. the microarray data were analysed by two-way anova using tmev statistical software, and gene expression was considered significantly different when p <. significant enrichment of go biological process categories were tested for using ease software with p <.  growth and biometry after months of the feeding trial, european sea bass fed vd exhibited significantly lower dgc than those fd. in addition, the fish of half-sibfamily g fed the vd had a significantly higher dgc than fish of half-sibfamily g fed vd, while there was no difference between these two half-sibfamilies when they were both fed fd. the hepatosomatic index was regulated by diet and genetic factors while the viscerosomatic index was only regulated by the genetic factor. there were no significant interactions between dietary and genetic factors for these two parameters. the feed efficiency in the duplicate tanks was and for the fd and and for the vd, respectively. mortality was not significantly different between dietary treatment and half-sibfamilies. initial and final weights and lengths, hsi and vsi. effects of diet factor and genetic factor on the biometric parameters are determined by two-way anova. results are expressed as mean +/- s.d. and significant differences are indicated by the p value. lipid and fatty acid compositions of the fillet the flesh lipid composition was significantly affected by dietary treatment. feeding vd significantly increased the percentage of saturated lipids in both the neutral lipids and phospholipids. the α-linolenic acid and linoleic acid contents were respectively 10-fold and 3-fold higher, in both lipid classes, when fish were fed vd. in addition, aa, epa and dha contents were around -fold lower in flesh of fish fed vd. the Σn- pufa/Σn- pufa ratio decreased in both the neutral lipid and phospholipid fractions in the flesh of european sea bass fed vd. microarray gene expression profiling a list of significant probes was obtained for the effect of diet factor, corresponding to unique transcripts with gene ontology annotation. among these regulated transcripts, exhibited higher levels in fish fed vd while were over-expressed in the liver of fish fed fd. for the family factor effect, total of significant probes were revealed corresponding to unique transcripts with go annotation. among these, exhibited higher levels in half-sibfamily g while were more abundant in half-sibfamily g. in addition, probes related to genes with functional annotation exhibited significant diet × genetic interactions. the main biological processes enriched out of those associated with genes that were over-expressed in fish fed vd were physiological process, metabolism, rna splicing, protein catabolism, aerobic respiration, sterol metabolism, carboxylic metabolism, amino acid metabolism, blood coagulation and hexose catabolism. the genes involved in carboxylic acid and sterol metabolism included fatty acid desaturase steroyl-coa desaturase nadh-cytochrome b reductase, isopentenyl-diphosphate delta-isomerase lanosterol 14-alpha demethylase, farnesyl pyrophosphate synthase, c- methylsterol oxidase and 3-hydroxy-3-methylglutaryl-coenzyme a reductase. apolipoprotein a-i, apolipoprotein b- and lipoprotein lipase, implicated in lipid transport, were also found to be up-expressed in fish fed vd. similarly, some genes involved in carbohydrate metabolism, such as glucose-6-phosphate 1-dehydrogenase, 6-phosphogluconate dehydrogenase and fructose- 6-bisphosphatase were also expressed at a higher level in the fish fed vd. expression levels of genes involved in protein metabolism and amino acid metabolism were also increased in fish fed vd. the fold-changes are indicated considering fd as the reference group.. in contrast, the main biological processes enriched associated with the genes that were lower-expressed in fish fed vd were particularly related to cellular process, cell communication and cell proliferation. in addition, some genes involved in the immune function were less expressed in fish fed vd. the comparison of hepatic transcriptomes between the two half-sibfamilies indicated that genes involved in immune function, such as complement component c c and c and mannan-binding lectin serine protease were more greatly expressed in half-sibfamily g than in half-sibfamily g. inversely, expression levels of nadh dehydrogenase genes, involved in the electron transport to the respiratory chain, were significantly lower in half-sibfamily g compared with half-sibfamily g. half-sibfamily g was also characterised by lower expression level of genes implicated in atp production and protein synthesis, such as ribosomal subunits than the half-sibfamily g. among the genes exhibiting an interaction between half-sibfamily and diet factors, were involved in metabolism. however, only the processes related to aromatic amino acid family and nucleotide metabolism were found to be over-represented among these genes. the fold-changes are indicated considering the half-sibfamily g as the reference group.. in order to validate the accuracy of the microarray data, the fads hmgcr, fabp angptl cxcl gck and lpl genes, which were spotted on the microarray, were also investigated by means of real-time pcr. the comparison of the gene expression pattern obtained through the real-time pcr and microarray approaches, revealed a correlation greater than. immune parameters lysozyme activity was significantly lower in fish fed vd than in fish fed fd, while the alternative complement activity was -times higher. there was no effect of the half-sibfamily factor on these activities. discussion the present work is the first investigation into the effect of an exclusively vegetable diet on the hepatic transcriptome in a marine fish species. it is also the first study to have explored the transcriptome of two half-sibfamilies of european sea bass exhibiting different capacities to grow on such a diet. the replacement of fm and fo with increasing levels of plant protein and oil sources for marine fish species can modify feed intake and conversion, which should be the major reason for associated growth delay. in the present study, there was a tendency for higher fe in fish fed fd compared with fish fed vd. however, this difference could not be statistically tested since fish were reared in only two tanks per diet condition. a vegetable diet is also known to potentially impact fish metabolism through regulation of gene expression, especially in the liver. analysis of the oligo-dna microarray data by two-way anova indicated that several hundred genes were differentially regulated according to diet or/and half-sibfamily factors. the accuracy of the present microarray data is validated by the similar gene pattern expression obtained from different oligonucleotides representing the same genes spotted on the array, as well as by the correlation shown between results of microarray and qpcr approaches. lc-pufa metabolism metabolism-related biological processes constitute the largest group among the go terms associated with the genes regulated in vd-fed fish. among these, genes involved in lipid metabolism and, particularly, in lc-pufa biosynthesis, were found to be up-regulated in fish fed vd. not surprisingly, delta-6-desaturase, steroyl-coa desaturase and nadh-cytochrome b reductase, involved in long chain fatty acid desaturation and/or elongation, were up-regulated in vd-fed fish. it is indeed established that in most fish species that the lc-pufa biosynthetic pathway is positively regulated in response to the use of a diet poor in lc-pufa but rich in pufa, although this regulation depends on fish species, degree of fish oil substitution, nature of the vegetable oil, and environmental parameters. the positive regulation of the lc-pufa biosynthetic pathway is in agreement with results obtained by transcriptomic approaches in salmonids fed vegetable oil. altogether, the expression data obtained in marine fish species and salmonids indicate that the different capacity of marine and fresh water species to grow on a lc-pufa-deprived diet does not seem to be due to a different transcriptional regulation of key genes involved in lipid synthesis, such as fads or scd indeed, the level of induction of fads expression found in this experiment is of similar amplitude to that observed in the liver and intestine of atlantic salmon. the stimulation of this biosynthetic pathway in fish fed a diet poor in lc-pufas can be explained by the fact that lc-pufas play several key physiological roles in vertebrates, particularly in fish. for example, fish are poikilothermic organisms and therefore need a high degree of unsaturation of lc-pufa included in membrane phospholipids to maintain phospholipid bilayer fluidity at reduced temperature. lc-pufa, especially arachidonic acid and eicosapentaenoic acid, are also precursors of eicosanoids, which are involved in pro and anti-inflammatory pathways. this hypothesis is reinforced by our data, indicating a stimulation of genes involved in phospholipids biosynthesis when fish were fed the vd. despite this stimulation of lc-pufa and the phospholipid biosynthesis pathway at the transcriptional level, our investigation of fatty acid profiles indicated that the amounts of lc-pufa, particularly eicosapentaenoic acid and docosahexaenoic acid, were still considerably lower in the flesh of fish from both half-sibfamilies fed vd in comparison with fish fed fd. this finding is in agreement with those previously obtained, which revealed that the stimulation of fads expression in fish fed a vegetable diet was not associated with an induction of its enzymatic activity, suggesting a post-transcriptional regulation of fads expression. such epa and dha deficiency can notably explain the growth deficiency observed in fish fed vd, as well as effects observed on immune function. concerning the genetic aspect, the comparable overall expression pattern of genes involved in lc-pufa synthesis in the liver, associated with similar lc-pufa profile in muscle, in both half-sibfamilies suggests that the differing capacities of these european sea bass half-sibfamilies to grow on a vegetable diet are not due to differing capacities to synthesize lc-pufa. sterol metabolism the present microarray data also indicate an increase in expression levels of genes involved in sterol metabolism in vd-fed fish. among these, isopentenyl-diphosphate delta-isomerase, lanosterol 14-alpha demethylase, farnesyl pyrophosphate synthase, c- methylsterol oxidase and 3-hydroxy-3-methylglutaryl-coenzyme a reductase genes are known to be implicated in the cholesterol metabolic pathway. more particularly, hmgcr, a transmembrane glycoprotein involved in the rate-limiting step of sterol biosynthesis, is increased, as shown in mammals. the stimulation of cholesterol biosynthesis in fish fed vd could be related to the difference in sterol composition between diets. indeed, while the fish diet is rich in cholesterol, the vegetable diet used in this experiment contains exclusively plant sterols, which have been shown to affect membrane properties by decreasing permeability and fluidity, and modifying phospholipid order in mammals. as a consequence, the increase in cholesterol biosynthesis could be a metabolic response to its deficiency in the diet, as well as a way to restore membrane properties by incorporation of endogenous cholesterol. since we did not measure the cholesterol content in the liver, flesh or blood, it is not possible for us to assess the capacity of european sea bass to compensate for possible dietary deficiencies in cholesterol through a regulation of its biosynthesis. moreover, similarly to the lc-pufa pathway, no significant difference of cholesterol biosynthetic regulation was observed between the half-sibfamilies. interestingly, several vd-stimulated genes involved in the lipogenic pathway are known to be molecular targets of sterol regulatory element binding proteins, which are key regulators of fatty acid and cholesterol synthesis. recent data indicating an up-regulation of the srebp- gene expression in european sea bass fed a vegetable diet could thus be due to such stimulations. lipid and sterol transport the present microarray data indicate that the stimulation of genes involved in fatty acid and cholesterol synthesis in vd-fed fish was associated with an over-expression of genes involved in their transport, such as apolipoproteins apoa and apob which are the major protein constituents of high and low density lipoprotein, respectively. the ldl, including apob are involved in the transport of cholesterol and lipids from the liver to other tissues. thus, up-regulation of apob combined with the induction of the expression of lipoprotein lipase, a key enzyme involved in the hydrolysis of triglyceride, suggests an increase in lipid transport and metabolism from the liver to tissues in fish fed vd. the decrease in angiopoietin-related protein that we observed in fish fed vd reinforces this idea since angptl suppresses ldl clearance via the inhibition of lpl activity. in parallel, we observed that the reverse transport of cholesterol not used by tissues via hdl to the liver was also stimulated in fish fed vd. indeed, apoa which participates in the transport of cholesterol to the liver by promoting cholesterol efflux from tissues and by acting as a cofactor for the lecithin cholesterol acyltransferase, exhibited higher transcript levels in fish fed vd. altogether these results reveal that another major response in the liver of european sea bass fed a vegetable diet is the stimulation of cholesterol synthesis and transport, irrespective of the half-sibfamily considered. carbohydrate metabolism lc-pufa and cholesterol biosynthesis require reducing power in the form of nadph. it is well documented in vertebrates, including fish, that nadph required for malonyl-coa synthesis is mainly supplied by the dehydrogenases of the pentose phosphate shunt. interestingly, our transcriptomic data indicate that the use of the vd induced a significant increase in the level of glucose-6-phosphate dehydrogenase. g6pd catalyses nadp+-linked oxidation of d-glucose-6-phosphate and has been shown to be a major contributor of nadph production for lipogenesis in atlantic salmon and european sea bass. moreover, our data indicate an increase in the expression of hexose-6-phosphate dehydrogenase and phosphogluconate dehydrogenase, enzymes of the pentose phosphate pathway that generate nadph, in fish fed vd. once synthesized, the resulting pentose sugar intermediate generated by the pentose phosphate pathway can be reconverted to intermediates of the glycolysis/gluconeogenesis pathway such as glyceraldehyde 3p or fructose 6p. in the liver of fish, it is known that glycolysis provides essential precursors for biosynthesis rather than pyruvate for oxidation. thus, the stimulation of fructose- 6-bisphosphatase and aldolase expression that we observe in fish fed vd could provide high levels of fructose-6-phosphate from glyceraldehyde 3p, then glucose 6p that serves as substrate for repeated passage in the pentose phosphate shunt. protein/amino-acid metabolism and atp synthesis our data revealed over-expression of genes involved in proteolysis and, more particularly, in proteasome activity and ubiquitin activity in fish fed vd, which is in total agreement with proteomic data obtained in rainbow trout, indicating a stimulation of proteolysis in fish fed vegetable diets. in our study, the stimulation of proteolysis in the fish fed the vegetable diet was associated with the induction of genes involved in amino acid metabolism and, more importantly, genes involved in glutamine metabolism. in addition, gmps, aadat, got and tat genes, which are implicated in transamination, were also stimulated in fish fed vd. the processes related to amino acid metabolism and, especially transamination, are important steps in the synthesis of some non-essential amino acids such as α-ketoglutarate. for example, the synthesis of α-ketoglutarate through transamination reactions could be used in the tca cycle to provide energy. interestingly, we found over-expression of genes encoding enzymes involved in the tca cycle, such as succinate dehydrogenase, isocitrate dehydrogenase and malate dehydrogenase, in fish fed vd. this stimulation of the tca cycle could be related to the higher levels of atp required for lc-pufa and cholesterol biosynthesis in fish fed vd. since marine fish have a low capacity to digest complex carbohydrates, in contrast to mammals, the use of proteins as an essential source of energy can thus explain the stimulation of the amino acid metabolism in fish fed vd. as shown in rainbow trout fed on a vegetable-based diet, the lower growth rate in fish fed vd in the present study could be associated with higher proteolytic activity compared with fish fed fd. interestingly, while both half-sibfamilies g and g exhibited similar proteolysis regulation, the expression of several genes involved in macromolecule biosynthesis, and particularly in protein biosynthesis, were up-regulated in half-sibfamily g. this result, suggesting a higher protein turnover in half-sibfamily g compared with half-sibfamily g when fish were fed vd, could be related to the higher growth rate observed in half-sibfamily g fed vd. as protein biosynthesis requires energy from atp hydrolysis, the higher protein biosynthesis in half-sibfamily could be related to a higher activity of mitochondrial atp production. accordingly, genes involved in atp biosynthesis and atp synthesis-coupled electron transport were found up-regulated in half-sibfamily g. the diet × half-sibfamily interaction that we found for the expression of genes involved in aromatic amino acid metabolism reinforces the difference in protein metabolism between the two half-sibfamilies. immune function in the present study, the vegetable diet used, based on linseed oil, was characterised by a very low ara content and poor levels of n- lc-pufa. eicosanoids derived from ara are known to be involved in the proliferation of hepatocytes and immune cells. as a consequence, the lower hepatosomatic index measured in fish fed vd could be linked to a lower hepatocyte proliferation due to a deficiency in ara, which was supported by down-expression of genes involved in cell proliferation in this dietary group. moreover, fatty acid imbalances can induce an immune deficiency in all vertebrates including fish. in particular, the Σ n-3/Σ n- fatty acid ratio is considered as a key element regulating immune cell structure, cell signalling, and eicosanoid production. the present microarray data revealed genes of the immune system, particularly the innate immune response, exhibiting lower expression in fish fed with vd. interleukin and c-x-c motif chemokine, which are chemotactic factors for granulocytes and monocytes, respectively, were found to be less expressed in fish fed vd, suggesting a down regulation of the innate defence system and the pro-inflammatory pathway in this group. deficiency of the inflammatory response was also in agreement with the higher levels of transcripts of fatty acid binding protein, whose expression in mammals has been shown to be restricted to the kupffer cells, and the down expression of the c-reactive protein, an acute phase protein synthesised by hepatocytes, in fish fed vd. decrease in inflammatory response can also be related to the low level of ara in the fish fed vd, which induces a reduction of prostaglandin synthesis derived from this fatty acid. our microarray data indeed show that prostaglandin e synthase, involved in the synthesis of pro-inflammatory prostaglandin e is down-regulated in fish fed vd, while prostaglandin e synthase, which has anti-inflammatory properties, exhibited higher messenger levels in fish fed vd. this depression of innate immune system, particularly pro-inflammatory activity, could also be partially explained by a defect in membrane properties in fish fed vd, as revealed by the down-regulation of a large number of genes related to cell communication, including factors such as cytokine receptor common subunit gamma, receptor-type tyrosine-protein phosphatase f or integrin beta, which are cell-surface receptor binding proteins and/or cell adhesion receptors involved in immune response. the depression of the innate immune response in fish fed vd was confirmed by the lower plasmatic lysozyme concentration and lower-expression of lysozyme g gene. surprisingly, the alternative complement pathway activity involved in the innate immune response, which we assessed by analysis of plasma parameters, showed a significantly higher level in fish fed vd. such an opposite regulation of the immune pathway revealed that different components of the immune systems can be regulated in opposite directions. interestingly, processes related to the humoral immune response were also over-represented among the genes up-regulated in half-sibfamily g. indeed, complement component c c and c genes showed higher expression levels in half-sibfamily g. however, the up-regulation of genes involved in the alternative complement pathway cannot be associated with an increase of the plasma alternative complement pathway activity, probably due to the complexity of factors and regulation levels involved in the regulation of this pathway. moreover, the higher expression of masp tnrfrf c and c genes involved in the inflammatory response might reflect higher inflammatory states in half-sibfamily g, which could be associated with a decrease in growth rate, as demonstrated in chicken. blood coagulation blood coagulation is another process involved in the innate immune system. lc-pufa and, more specifically, epa, dha and ara are precursors for eicosanoid synthesis involved in the control of the blood coagulation. as mentioned above, the use of a diet composed of vegetable protein and oil induces modifications in the membrane phospholipid composition, with possible consequences for eicosanoid production and the blood coagulation process. we found that the use of the vegetable diet induced an increase in the expression of genes involved in the blood coagulation pathway. among these genes, prothrombin, coagulation factor ×, fibrinogen beta chain, fibrinogen gamma chain and coagulation factor vii were positively involved in the blood coagulation process. on the basis of these results, the use of a vd seems to cause pro-coagulant action by the stimulation of the blood coagulation pathway, which is in agreement with our visual observation of plasma clotting. in agreement with these results, tavares-dias showed that, contrary to a vegetable diet, dietary enrichment in long-chain n- fatty acids has a strong hypocoagulant action. in addition, the regulation of the coagulation pathway is complex and under the control of several negative factors that maintain a physiological homeostasis. the pro-coagulation effect in response to a vegetable diet is notably reinforced by some of these genes, which also exhibited higher expression in fish fed vd.  the nutrigenomic approach used here has revealed several new genes and related biological processes regulated by a vegetable diet. in particular, genes involved in lipid metabolism, protein/amino acid metabolism, carbohydrate metabolism, immune function, blood coagulation and the rna splicing process were expressed at a higher level in fish fed with vd. the comparison of transcriptomic response in two half-sibfamilies of fish exhibiting different growth rates when fed the vegetable diet also revealed some biological processes related to protein turnover and immune response, potentially due to better adaptation to this diet. finally, in the context of developing novel diets for aquaculture and selecting fish families exhibiting higher adaptation to fish oil and meal substitution, this work enabled us to pinpoint potentially useful new molecular markers for identifying the physiological effects of a vegetable diet, as well as a family exhibiting a phenotype of interest. list of abbreviations ala: α-linolenic acid; ara: arachidonic acid; dgc: daily growth coefficient; dha: docosahexaenoic acid; dm: dry matter; ef1: elongation factor 1; epa: eicosapentaenoic acid; fd: fish diet; fe: feed efficiency; go: gene ontology; hdl: high density lipoprotein; hsi: hepatosomaticindex; la: linoleic acid; lc-pufa: long chain polyunsaturated fatty acid; ldl: low density lipoprotein; nl: neutral lipid; pl: polar lipid; pufa: polyunsaturated fatty acid; tl: total lipid; vd: vegetable diet; vsi: viscerosomatic index. authors' contributions fg and dm co-ordinated the sampling procedures, biochemical studies and microarray data analysis, and drafted the manuscript. sf and lb carried out the sample preparation, hybridizations and microarray data normalizations. sk formulated the experimental diet. mv participated in the design of the study and performed analysis of growth parameters. cc and jlz conceived and designed the study and led its coordination. cq co-ordinated dedicated experiments on immune function. all authors read and approved the final manuscript. supplementary material additional file 1 comparison of ct values for ef gene between the dietary groups. quantitative pcr did not reveal any significant difference between ct values for ef between the dietary groups. click here for file additional file 2 fatty acid composition in muscle of two european sea bass half-sibfamilies fed fd or vd. composition in terms of the main fatty acids in neutral lipid and phospholipid fractions in muscle of each of the half-sibfamilies of european sea bass fed fd or vd. effects of diet factor and half-sibfamily factor on fatty acid composition were determined by two-way anova. results are expressed as mean +/- s.d. and significant differences are indicated by the p value. click here for file additional file 3 significantly enriched biological processes associated with genes regulated by diet and half-sibfamily factors. the main biological processes enriched out of those associated with genes that were over-expressed in fish fed vd were related to physiological process, metabolism, rna splicing, protein catabolism, aerobic respiration, blood coagulation and hexose catabolism. in contrast, the main biological processes associated with the genes lower-expressed in fish fed vd were related to cellular process, cell communication and cell proliferation. regarding half-sibfamily factors, biological process related to humoral immune response was shown to be over-represented within genes up-expressed in half-sibfamily g while processes related to energy pathways were enriched within genes up-expressed in half-sibfamily g. click here for file acknowledgements this work was supported by marine genomics europe, europôle mer axis 1: "genomics and'blue chemistry' ". the authors thank h. le delliou for fatty acid analysis, f. ruelle for european sea bass breeding, f. terrier for making the experimental diets, f. lamour, c. huelvan, e. desbruyères for the measurement of immune parameters on plasma, dr. a. henaut for discussion on microarrays analysis and h. mccombie for english revisions.
biological processes related to protein catabolism, amino acid transaminations, rna splicing and blood coagulation were also found to be regulated by diet, while the expression of genes involved in protein and atp synthesis differed between the half-sibfamilies. efforts towards utilisation of diets without fish meal or fish oil in finfish aquaculture have been being made for more than two decades. we report the first results obtained using a transcriptomic approach on the liver of two half-sibfamilies of the european sea bass that exhibit similar growth rates when fed a fish-based diet, but significantly different growth rates when fed an all-plant diet. no significant diet-family interaction for the regulation of lc-pufa biosynthesis pathways could be detected by microarray analysis. statistical analysis identified unique annotated genes differentially expressed between groups of fish fed the two diets, genes regulated by genetic factors, and genes that exhibited diet-family interactions. metabolic responses to substitution of fishery products have been shown to impact growth performance and immune system of fish as well as their subsequent nutritional value, particularly in marine fish species, which exhibit low capacity for biosynthesis of long-chain poly-unsaturated fatty acids. in addition, the combination of our transcriptomic data with an analysis of plasmatic immune parameters revealed a stimulation of complement activity associated with an immunodeficiency in the fish fed vd, and different inflammatory status between the two half-sibfamilies. overall, the combined gene expression, compositional and biochemical studies demonstrated a large panel of metabolic and physiological effects induced by total substitution of both fm and fo in the diets of european sea bass and revealed physiological characteristics associated with the two half-sibfamilies. overall gene expression was analysed using oligo dna microarrays. the main objective of the present study was to analyse the effects of a plant-based diet on the hepatic transcriptome of european sea bass.
 efforts towards utilisation of diets without fish meal or fish oil in finfish aquaculture have been being made for more than two decades. metabolic responses to substitution of fishery products have been shown to impact growth performance and immune system of fish as well as their subsequent nutritional value, particularly in marine fish species, which exhibit low capacity for biosynthesis of long-chain poly-unsaturated fatty acids. the main objective of the present study was to analyse the effects of a plant-based diet on the hepatic transcriptome of european sea bass.  we report the first results obtained using a transcriptomic approach on the liver of two half-sibfamilies of the european sea bass that exhibit similar growth rates when fed a fish-based diet, but significantly different growth rates when fed an all-plant diet. overall gene expression was analysed using oligo dna microarrays. statistical analysis identified unique annotated genes differentially expressed between groups of fish fed the two diets, genes regulated by genetic factors, and genes that exhibited diet-family interactions. the expression of several genes involved in the lc-pufa and cholesterol biosynthetic pathways was found to be up-regulated in fish fed vd, suggesting a stimulation of the lipogenic pathways. no significant diet-family interaction for the regulation of lc-pufa biosynthesis pathways could be detected by microarray analysis. this result was in agreement with lc-pufa profiles, which were found to be similar in the flesh of the two half-sibfamilies. in addition, the combination of our transcriptomic data with an analysis of plasmatic immune parameters revealed a stimulation of complement activity associated with an immunodeficiency in the fish fed vd, and different inflammatory status between the two half-sibfamilies. biological processes related to protein catabolism, amino acid transaminations, rna splicing and blood coagulation were also found to be regulated by diet, while the expression of genes involved in protein and atp synthesis differed between the half-sibfamilies.  overall, the combined gene expression, compositional and biochemical studies demonstrated a large panel of metabolic and physiological effects induced by total substitution of both fm and fo in the diets of european sea bass and revealed physiological characteristics associated with the two half-sibfamilies. 