Understanding how high stocking densities and concurrent limited oxygen availability drive social cohesion and adaptive features in regulatory growth, antioxidant defense and lipid metabolism in farmed gilthead sea bream (Sparus aurata).

The study combined the use of biometric, behavioral, physiological and external tissue damage scoring systems to better understand how high stocking densities drive schooling behavior and other adaptive features during the finishing growing phase of farmed gilthead sea bream in the Western Mediterranean. Fish were grown at three different final stocking densities (LD, 8.5 kg/m3; MD, 17 kg/m3; HD, 25 kg/m3). Water oxygen concentration varied between 5 and 6 ppm in LD fish to 3-4 ppm in HD fish with the summer rise of water temperature from 19°C to 26°C (May-July). HD fish showed a reduction of feed intake and growth rates, but they also showed a reinforced social cohesion with a well-defined endogenous swimming activity rhythm with feeding time as a main synchronization factor. The monitored decrease of the breathing/swimming activity ratio by means of the AEFishBIT data-logger also indicated a decreased energy partitioning for growth in the HD environment with a limited oxygen availability. Plasma glucose and cortisol levels increased with the rise of stocking density, and the close association of glycaemia with the expression level of antioxidant enzymes (mn-sod, gpx4, prdx5) in liver and molecular chaperones (grp170, grp75) in skeletal muscle highlighted the involvement of glucose in redox processes via rerouting in the pentose-phosphate-pathway. Other adaptive features included the depletion of oxidative metabolism that favored lipid storage rather than fatty acid oxidation to decrease the oxygen demand as last electron acceptor in the mitochondrial respiratory chain. This was coincident with the metabolic readjustment of the Gh/Igf endocrine-growth cascade that promoted the regulation of muscle growth at the local level rather than a systemic action via the liver Gh/Igf axis. Moreover, correlation analyses within HD fish displayed negative correlations of hepatic transcripts of igf1 and igf2 with the data-logger measurements of activity and respiration, whereas the opposite was found for muscle igf2, ghr1 and ghr2. This was indicative of a growth-regulatory transition that supported a proactive instead of a reactive behavior in HD fish, which was considered adaptive to preserve an active and synchronized feeding behavior with a minimized risk of oxidative stress and epidermal skin damage.

SAMBA: Structure-Learning of Aquaculture Microbiomes Using a Bayesian Approach.

Gut microbiomes of fish species consist of thousands of bacterial taxa that interact among each other, their environment, and the host. These complex networks of interactions are regulated by a diverse range of factors, yet little is known about the hierarchy of these interactions. Here, we introduce SAMBA (Structure-Learning of Aquaculture Microbiomes using a Bayesian Approach), a computational tool that uses a unified Bayesian network approach to model the network structure of fish gut microbiomes and their interactions with biotic and abiotic variables associated with typical aquaculture systems. SAMBA accepts input data on microbial abundance from 16S rRNA amplicons as well as continuous and categorical information from distinct farming conditions. From this, SAMBA can create and train a network model scenario that can be used to (i) infer information of how specific farming conditions influence the diversity of the gut microbiome or pan-microbiome, and (ii) predict how the diversity and functional profile of that microbiome would change under other variable conditions. SAMBA also allows the user to visualize, manage, edit, and export the acyclic graph of the modelled network. Our study presents examples and test results of Bayesian network scenarios created by SAMBA using data from a microbial synthetic community, and the pan-microbiome of gilthead sea bream (Sparus aurata) in different feeding trials. It is worth noting that the usage of SAMBA is not limited to aquaculture systems as it can be used for modelling microbiome-host network relationships of any vertebrate organism, including humans, in any system and/or ecosystem.

Client Applications and Server-Side Docker for Management of RNASeq and/or VariantSeq Workflows and Pipelines of the GPRO Suite.

The GPRO suite is an in-progress bioinformatic project for -omics data analysis. As part of the continued growth of this project, we introduce a client- and server-side solution for comparative transcriptomics and analysis of variants. The client-side consists of two Java applications called "RNASeq" and "VariantSeq" to manage pipelines and workflows based on the most common command line interface tools for RNA-seq and Variant-seq analysis, respectively. As such, "RNASeq" and "VariantSeq" are coupled with a Linux server infrastructure (named GPRO Server-Side) that hosts all dependencies of each application (scripts, databases, and command line interface software). Implementation of the Server-Side requires a Linux operating system, PHP, SQL, Python, bash scripting, and third-party software. The GPRO Server-Side can be installed, via a Docker container, in the user's PC under any operating system or on remote servers, as a cloud solution. "RNASeq" and "VariantSeq" are both available as desktop (RCP compilation) and web (RAP compilation) applications. Each application has two execution modes: a step-by-step mode enables each step of the workflow to be executed independently, and a pipeline mode allows all steps to be run sequentially. "RNASeq" and "VariantSeq" also feature an experimental, online support system called GENIE that consists of a virtual (chatbot) assistant and a pipeline jobs panel coupled with an expert system. The chatbot can troubleshoot issues with the usage of each tool, the pipeline jobs panel provides information about the status of each computational job executed in the GPRO Server-Side, while the expert system provides the user with a potential recommendation to identify or fix failed analyses. Our solution is a ready-to-use topic specific platform that combines the user-friendliness, robustness, and security of desktop software, with the efficiency of cloud/web applications to manage pipelines and workflows based on command line interface software.

Genetics and Nutrition Drive the Gut Microbiota Succession and Host-Transcriptome Interactions through the Gilthead Sea Bream (Sparus aurata) Production Cycle.

Fish genetically selected for growth (GS) and reference (REF) fish were fed with CTRL (15% FM, 5-7% FO) or FUTURE (7.5% FM, 10% poultry meal, 2.2% poultry oil + 2.5% DHA-algae oil) diets during a 12-months production cycle. Samples from initial (t0; November 2019), intermediate (t1; July 2020) and final (t2; November 2020) sampling points were used for Illumina 16S rRNA gene amplicon sequencing of the adherent microbiota of anterior intestine (AI). Samples from the same individuals (t1) were also used for the gene expression profiling of AI by RNA-seq, and subsequent correlation analyses with microbiota abundances. Discriminant analyses indicated the gut bacterial succession along the production cycle with the proliferation of some valuable taxa for facing seasonality and different developmental stages. An effect of genetic background was evidenced along time, decreasing through the progression of the trial, namely the gut microbiota of GS fish was less influenced by changes in diet composition. At the same time, these fish showed wider transcriptomic landmarks in the AI to cope with these changes. Our results highlighted an enhanced intestinal sphingolipid and phospholipid metabolism, epithelial turnover and intestinal motility in GS fish, which would favour their improved performance despite the lack of association with changes in gut microbiota composition. Furthermore, in GS fish, correlation analyses supported the involvement of different taxa with the down-regulated expression of pro-inflammatory markers and the boosting of markers of extracellular remodelling and response to bacterium. Altogether, these findings support the combined action of the gut microbiome and host transcriptionally mediated effects to preserve and improve gut health and function in a scenario of different growth performance and potentiality.

Diet and Host Genetics Drive the Bacterial and Fungal Intestinal Metatranscriptome of Gilthead Sea Bream.

The gut microbiota is now recognised as a key target for improving aquaculture profit and sustainability, but we still lack insights into the activity of microbes in fish mucosal surfaces. In the present study, a metatranscriptomic approach was used to reveal the expression of gut microbial genes in the farmed gilthead sea bream. Archaeal and viral transcripts were a minority but, interestingly and contrary to rRNA amplicon-based studies, fungal transcripts were as abundant as bacterial ones, and increased in fish fed a plant-enriched diet. This dietary intervention also drove a differential metatranscriptome in fish selected for fast and slow growth. Such differential response reinforced the results of previously inferred metabolic pathways, enlarging, at the same time, the catalogue of microbial functions in the intestine. Accordingly, vitamin and amino acid metabolism, and rhythmic and symbiotic processes were mostly shaped by bacteria, whereas fungi were more specifically configuring the host immune, digestive, or endocrine processes.

Cross-Talk Between Intestinal Microbiota and Host Gene Expression in Gilthead Sea Bream (Sparus aurata) Juveniles: Insights in Fish Feeds for Increased Circularity and Resource Utilization.

New types of fish feed based on processed animal proteins (PAPs), insect meal, yeast, and microbial biomasses have been used with success in gilthead sea bream. However, some drawback effects on feed conversion and inflammatory systemic markers were reported in different degrees with PAP- and non-PAP-based feed formulations. Here, we focused on the effects of control and two experimental diets on gut mucosal-adherent microbiota, and how it correlated with host transcriptomics at the local (intestine) and systemic (liver and head kidney) levels. The use of tissue-specific PCR-arrays of 93 genes in total rendered 13, 12, and 9 differentially expressed (DE) genes in the intestine, liver, and head kidney, respectively. Illumina sequencing of gut microbiota yielded a mean of 125,350 reads per sample, assigned to 1,281 operational taxonomic unit (OTUs). Bacterial richness and alpha diversity were lower in fish fed with the PAP diet, and discriminant analysis displayed 135 OTUs driving the separation between groups with 43 taxa correlating with 27 DE genes. The highest expression of intestinal pcna and alpi was achieved in PAP fish with intermediate values in non-PAP, being the pro-inflammatory action of alpi associated with the presence of Psychrobacter piscatorii. The intestinal muc13 gene was down-regulated in non-PAP fish, with this gene being negatively correlated with anaerobic (Chloroflexi and Anoxybacillus) and metal-reducing (Pelosinus and Psychrosinus) bacteria. Other inflammatory markers (igm, il8, tnfα) were up-regulated in PAP fish, positively correlating the intestinal igm gene with the inflammasome activator Escherichia/Shigella, whereas the systemic expression of il8 and tnfα was negatively correlated with the Bacilli class in PAP fish and positively correlated with Paracoccus yeei in non-PAP fish. Overall changes in the expression pattern of il10, galectins (lgals1, lgals8), and toll-like receptors (tlr2, tlr5, tlr9) reinforced the anti-inflammatory profile of fish fed with the non-PAP diet, with these gene markers being associated with a wide range of OTUs. A gut microbiota-liver axis was also established, linking the microbial generation of short chain fatty acids with the fueling of scd1- and elovl6-mediated lipogenesis. In summary, by correlating the microbiome with host gene expression, we offer new insights in the evaluation of fish diets promoting gut and metabolism homeostasis, and ultimately, the health of farmed fish.

Targeting the Mild-Hypoxia Driving Force for Metabolic and Muscle Transcriptional Reprogramming of Gilthead Sea Bream (Sparus aurata) Juveniles.

On-growing juveniles of gilthead sea bream were acclimated for 45 days to mild-hypoxia (M-HYP, 40-60% O2 saturation), whereas normoxic fish (85-90% O2 saturation) constituted two different groups, depending on if they were fed to visual satiety (control fish) or pair-fed to M-HYP fish. Following the hypoxia conditioning period, all fish were maintained in normoxia and continued to be fed until visual satiation for 3 weeks. The time course of hypoxia-induced changes was assessed by changes in blood metabolic landmarks and muscle transcriptomics before and after exhaustive exercise in a swim tunnel respirometer. In M-HYP fish, our results highlighted a higher contribution of aerobic metabolism to whole energy supply, shifting towards a higher anaerobic fitness following normoxia restoration. Despite these changes in substrate preference, M-HYP fish shared a persistent improvement in swimming performance with a higher critical speed at exercise exhaustion. The machinery of muscle contraction and protein synthesis and breakdown was also largely altered by mild-hypoxia conditioning, contributing this metabolic re-adjustment to the positive regulation of locomotion and to the catch-up growth response during the normoxia recovery period. Altogether, these results reinforce the presence of large phenotypic plasticity in gilthead sea bream, and highlights mild-hypoxia as a promising prophylactic measure to prepare these fish for predictable stressful events.

Effects of genetics and early-life mild hypoxia on size variation in farmed gilthead sea bream (Sparus aurata).

The present study evaluated, in an 18-month gilthead sea bream trial, the time course effects of genetics on individual size variation and growth compensation processes in families selected by heritable growth in the PROGENSA® breeding program. Families categorized as fast, intermediate, and slow growing had different growth trajectories with a more continuous growth in fast growth families. This feature was coincident with a reduced size variation at the beginning of the trial that clustered together the half-sib families sharing the same father. Regression analysis evidenced that the magnitude of compensatory growth was proportional to the initial size variation with no rescaling of families at this stage. By contrast, the finishing growth depensation process can mask, at least partially, the previous size convergence. This reflects the different contribution across the production cycle of genetics in growth. How early-life experiences affect growth compensation at juvenile stages was also evaluated in a separate cohort, and intriguingly, a first mild-hypoxia pulse at 60-81 days post-hatching (dph) increased survival rates by 10%, preventing growth impairment when fish were exposed to a second hypoxia episode (112-127 dph). The early hypoxia experience did not have a negative impact on growth compensatory processes at juvenile stages. By contrast, a diminished capacity for growth compensation was found with repeated or late hypoxia experiences. All this reinforces the use of size variation as a main criterion for improving intensive fish farming and selective breeding.

Genetic selection for growth drives differences in intestinal microbiota composition and parasite disease resistance in gilthead sea bream.

BACKGROUND: The key effects of intestinal microbiota in animal health have led to an increasing interest in manipulating these bacterial populations to improve animal welfare. The aquaculture sector is no exception and in the last years, many studies have described these populations in different fish species. However, this is not an easy task, as intestinal microbiota is composed of very dynamic populations that are influenced by different factors, such as diet, environment, host age, and genetics. In the current study, we aimed to determine whether the genetic background of gilthead sea bream (Sparus aurata) influences the intestinal microbial composition, how these bacterial populations are modulated by dietary changes, and the effect of selection by growth on intestinal disease resistance. To that aim, three different groups of five families of gilthead sea bream that were selected during two generations for fast, intermediate, or slow growth (F3 generation) were kept together in the same open-flow tanks and fed a control or a well-balanced plant-based diet during 9 months. Six animals per family and dietary treatment were sacrificed and the adherent bacteria from the anterior intestinal portion were sequenced. In parallel, fish of the fast- and slow-growth groups were infected with the intestinal parasite Enteromyxum leei and the disease signs, prevalence, intensity, and parasite abundance were evaluated. RESULTS: No differences were detected in alpha diversity indexes among families, and the core bacterial architecture was the prototypical composition of gilthead sea bream intestinal microbiota, indicating no dysbiosis in any of the groups. The plant-based diet significantly changed the microbiota in the intermediate- and slow-growth families, with a much lower effect on the fast-growth group. Interestingly, the smaller changes detected in the fast-growth families potentially accounted for more changes at the metabolic level when compared with the other families. Upon parasitic infection, the fast-growth group showed significantly lower disease signs and parasite intensity and abundance than the slow-growth animals. CONCLUSIONS: These results show a clear genome-metagenome interaction indicating that the fast-growth families harbor a microbiota that is more flexible upon dietary changes. These animals also showed a better ability to cope with intestinal infections. Video Abstract.

Stearoyl-CoA desaturase (scd1a) is epigenetically regulated by broodstock nutrition in gilthead sea bream (Sparus aurata).

The aim of this study was to generate new knowledge on fish epigenetics, assessing the effects of linolenic acid (ALA) conditioning of broodstock in the offspring of the marine fish Sparus aurata. Attention was focused on gene organization, methylation signatures and gene expression patterns of fatty acid desaturase 2 (fads2) and stearoyl-CoA desaturase 1a (scd1a). Blat searches in the genomic IATS-CSIC database (www.nutrigroup-iats.org/seabreamdb) highlighted a conserved exon-intron organization, a conserved PUFA response region, and CG islands at the promoter regions of each gene. The analysed CpG positions in the fads2 promoter were mostly hypomethylated and refractory to broodstock nutrition. The same response was achieved after conditioning of juvenile fish to low water oxygen concentrations, thus methylation susceptibility at individual CpG sites seems to be stringently regulated in fish of different origin and growth trajectories. Conversely, the scd1a promoter was responsive to broodstock nutrition and the offspring of parents fed the ALA-rich diet shared an increased DNA-methylation, mainly in CpG sites neighbouring SP1 and HNF4α binding sites. Cytosine methylation at these sites correlated inversely with the hepatic scd1a expression of the offspring. Co-expression analyses supported that the HNF4α-dependent regulation of scd1a is affected by DNA methylation. The phenotypic output is a regulated liver fat deposition through changes in scd1 expression, which would also allow the preservation of fatty acid unsaturation levels in fish fed reduced levels of n-3 LC-PUFA. Collectively, these findings reveal a reliable mechanism by which parent's nutrition can shape scd1a gene expression in the fish offspring.

Sex, Age, and Bacteria: How the Intestinal Microbiota Is Modulated in a Protandrous Hermaphrodite Fish.

Intestinal microbiota is key for many host functions, such as digestion, nutrient metabolism, disease resistance, and immune function. With the growth of the aquaculture industry, there has been a growing interest in the manipulation of fish gut microbiota to improve welfare and nutrition. Intestinal microbiota varies with many factors, including host species, genetics, developmental stage, diet, environment, and sex. The aim of this study was to compare the intestinal microbiota of adult gilthead sea bream (Sparus aurata) from three groups of age and sex (1-year-old males and 2- and 4-year-old females) maintained under the same conditions and fed exactly the same diet. Microbiota diversity and richness did not differ among groups. However, bacterial composition did, highlighting the presence of Photobacterium and Vibrio starting at 2 years of age (females) and a higher presence of Staphylococcus and Corynebacterium in 1-year-old males. The core microbiota was defined by 14 Operational Taxonomic Units (OTUs) and the groups that showed more OTUs in common were 2- and 4-year-old females. Discriminant analyses showed a clear separation by sex and age, with bacteria belonging to the phyla Firmicutes, Proteobacteria and Actinobacteria driving the separation. Pathway analysis performed with the inferred metagenome showed significant differences between 1-year-old males and 4-year-old females, with an increase in infection-related pathways, nitrotoluene degradation and sphingolipid metabolism, and a significant decrease in carbohydrate metabolism pathways with age. These results show, for the first time, how intestinal microbiota is modulated in adult gilthead sea bream and highlight the importance of reporting age and sex variables in these type of studies in fish.

Acting locally - affecting globally: RNA sequencing of gilthead sea bream with a mild Sparicotyle chrysophrii infection reveals effects on apoptosis, immune and hypoxia related genes.

BACKGROUND: Monogenean flatworms are the main fish ectoparasites inflicting serious economic losses in aquaculture. The polyopisthocotylean Sparicotyle chrysophrii parasitizes the gills of gilthead sea bream (GSB, Sparus aurata) causing anaemia, lamellae fusion and sloughing of epithelial cells, with the consequent hypoxia, emaciation, lethargy and mortality. Currently no preventive or curative measures against this disease exist and therefore information on the host-parasite interaction is crucial to find mitigation solutions for sparicotylosis. The knowledge about gene regulation in monogenean-host models mostly comes from freshwater monopysthocotyleans and almost nothing is known about polyopisthocotyleans. The current study aims to decipher the host response at local (gills) and systemic (spleen, liver) levels in farmed GSB with a mild natural S. chrysophrii infection by transcriptomic analysis. RESULTS: Using Illumina RNA sequencing and transcriptomic analysis, a total of 2581 differentially expressed transcripts were identified in infected fish when compared to uninfected controls. Gill tissues in contact with the parasite (P gills) displayed regulation of fewer genes (700) than gill portions not in contact with the parasite (NP gills) (1235), most likely due to a local silencing effect of the parasite. The systemic reaction in the spleen was much higher than that at the parasite attachment site (local) (1240), and higher than in liver (334). NP gills displayed a strong enrichment of genes mainly related to immune response and apoptosis. Processes such as apoptosis, inflammation and cell proliferation dominated gills, whereas inhibition of apoptosis, autophagy, platelet activation, signalling and aggregation, and inflammasome were observed in spleen. Proteasome markers were increased in all tissues, whereas hypoxia-related genes were down-regulated in gills and spleen. CONCLUSIONS: Contrasting forces seem to be acting at local and systemic levels. The splenic down-regulation could be part of a hypometabolic response, to counteract the hypoxia induced by the parasite damage to the gills and to concentrate the energy on defence and repair responses. Alternatively, it can be also interpreted as the often observed action of helminths to modify host immunity in its own interest. These results provide the first toolkit for future studies towards understanding and management of this parasitosis.

Somatotropic Axis Regulation Unravels the Differential Effects of Nutritional and Environmental Factors in Growth Performance of Marine Farmed Fishes.

The Gh/Prl/Sl family has evolved differentially through evolution, resulting in varying relationships between the somatotropic axis and growth rates within and across fish species. This is due to a wide range of endogenous and exogenous factors that make this association variable throughout season and life cycle, and the present minireview aims to better define the nutritional and environmental regulation of the endocrine growth cascade over precisely defined groups of fishes, focusing on Mediterranean farmed fishes. As a result, circulating Gh and Igf-i are revitalized as reliable growth markers, with a close association with growth rates of gilthead sea bream juveniles with deficiency signs in both macro- or micro-nutrients. This, together with other regulated responses, promotes the use of Gh and Igf-i as key performance indicators of growth, aerobic scope, and nutritional condition in gilthead sea bream. Moreover, the sirtuin-energy sensors might modulate the growth-promoting action of somatotropic axis. In this scenario, transcripts of igf-i and gh receptors mirror changes in plasma Gh and Igf-i levels, with the ghr-i/ghr-ii expression ratio mostly unaltered over season. However, this ratio is nutritionally regulated, and enriched plant-based diets or diets with specific nutrient deficiencies downregulate hepatic ghr-i, decreasing the ghr-i/ghr-ii ratio. The same trend, due to a ghr-ii increase, is found in skeletal muscle, whereas impaired growth during overwintering is related to increase in the ghr-i/ghr-ii and igf-ii/igf-i ratios in liver and skeletal muscle, respectively. Overall, expression of insulin receptors and igf receptors is less regulated, though the expression quotient is especially high in the liver and muscle of sea bream. Nutritional and environmental regulation of the full Igf binding protein 1-6 repertoire remains to be understood. However, tissue-specific expression profiling highlights an enhanced and nutritionally regulated expression of the igfbp-1/-2/-4 clade in liver, whereas the igfbp-3/-5/-6 clade is overexpressed and regulated in skeletal muscle. The somatotropic axis is, therefore, highly informative of a wide-range of growth-disturbing and stressful stimuli, and multivariate analysis supports its use as a reliable toolset for the assessment of growth potentiality and nutrient deficiencies and requirements, especially in combination with selected panels of other nutritionally regulated metabolic biomarkers.