# Transcriptomic response to aquaculture intensification in Nile tilapia

Deiene Rodriguez-Barreto, Olivier REY, Tamsyn Uren-Webster, Giovanni Castaldo, Sonia Consuegra, Carlos Garcia de Leaniz

► To cite this version:

Deiene Rodriguez-Barreto, Olivier REY, Tamsyn Uren-Webster, Giovanni Castaldo, Sonia Consuegra, et al.. Transcriptomic response to aquaculture intensification in Nile tilapia. Evolutionary Applications, 2019, 12 (9), pp.1757-1771. 10.1111/eva.12830. hal-04080755

HAL Id: hal-04080755

https://hal.science/hal-04080755

Submitted on 10 Oct 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

ORIGINAL ARTICLE

<figure><img src="image_3.png" /></figure>

WILEY

# Transcriptomic response to aquaculture intensification in Nile tilapia

Deiene Rodriguez-Barreto<fnref n="1" /> |

Olivier Rey<fnref n="1" />,<fnref n="2" /> |

Tamsyn M. Uren-Webster<fnref n="1" /> |

Giovanni Castaldo<fnref n="1" />,<fnref n="3" /> |

Sonia Consuegra<fnref n="1" /> ID |

Carlos Garcia de Leaniz<fnref n="1" /> ID

**Correspondence**

Carlos Garcia de Leaniz, Centre for Sustainable Aquatic Research (CSAR), College of Science, Swansea University, Swansea SA2 8PP, UK.

Email: c.garciadeleaniz@swansea.ac.uk

**Funding information**

NRN-LCEE AquaWales; ERDF WEFO
SMARTAQUA

**Abstract**

To meet future global demand for fish protein, more fish will need to be farmed using fewer resources, and this will require the selection of nonaggressive individuals that perform well at high densities. Yet, the genetic changes underlying loss of aggression and adaptation to crowding during aquaculture intensification are largely unknown. We examined the transcriptomic response to aggression and crowding in Nile tilapia, one of the oldest and most widespread farmed fish, whose social structure shifts from social hierarchies to shoaling with increasing density. A mirror test was used to quantify aggression and skin darkening (a proxy for stress) of fish reared at low and high densities, and gene expression in the hypothalamus was analysed among the most and least aggressive fish at each density. Fish reared at high density were darker, had larger brains, were less active and less aggressive than those reared at low density and had differentially expressed genes consistent with a reactive stress-coping style and activation of the hypothalamus-pituitary-interrenal (HPI) axis. Differences in gene expression among aggressive fish were accounted for by density and the interaction between density and aggression levels, whereas for nonaggressive fish differences in gene expression were associated with individual variation in skin brightness and social stress. Thus, the response to crowding in Nile tilapia is context dependent and involves different neuroendocrine pathways, depending on social status. Knowledge of genes associated with the response to crowding may pave the way for more efficient fish domestication, based on the selection of nonaggressive individuals with increasing tolerance to chronic stress necessary for aquaculture intensification.

**KEYWORDS**

aggression, aquaculture intensification, crowding, fish domestication, gene expression, HPI axis, stress response

<footnote n="1">¹Centre for Sustainable Aquatic Research (CSAR), College of Science, Swansea University, Swansea, UK</footnote>

<footnote n="2">²Université de Perpignan Via Domitia, Perpignan, France</footnote>

<footnote n="3">³Systemic Physiological and Ecotoxicological Research, Department of Biology, University of Antwerp, Antwerp, Belgium</footnote>
# 1 | INTRODUCTION

To meet global fish demand, aquaculture will need to produce more fish with less food, less water and less space in the future (Godfray et al., 2010; Goldburg & Naylor, 2005). This will require the selection of fish that can thrive under crowded conditions and adapt well to life in captivity (Huntingford, 2004; Huntingford et al., 2006) while maintaining high welfare standards (Ashley, 2007; FAWC, 2014). But to select fish that perform well at high densities requires knowledge on the genetic basis of social behaviours, which for most farmed fish is lacking. Fish domestication involves profound changes in social behaviour (Huntingford et al., 2006), but knowledge on the expression of genes underlying social behaviours has mostly focused on model or "simple" organisms (Sokolowski, 2010) or in relation to the production of terrestrial livestock (Mormède, 2005).

The genetic basis of social behaviours has historically been difficult to study (Blumstein et al., 2010), as these encompass complex phenotypic traits that often depend on genotype by environment interactions (Komers, 1997; Robinson, Fernald, & Clayton, 2008). Aggression is one of the most ubiquitous social behaviours (Maxson & Canastar, 2005), but also one of the most labile ones, as it is easily influenced by the environment experienced during early development (Fernald, 2012; Marks, West, Bagatto, Moore, & Taylor, 2005; Maruska, 2015; Trainor, Lin, Finy, Rowland, & Nelson, 2007). Individuals resort to aggressive behaviour to protect themselves and their progeny, to compete for resources and mating partners or to establish a social rank that might accrue future benefits (Nelson & Trainor, 2007). A loss of aggression is possibly the most pronounced behavioural change that accompanies animal domestication (Price, 1999), but there is controversy about the underlying molecular mechanisms (Wilkins, Wrangham, & Fitch, 2014). In cultured fish, rearing density has a marked influence on aggression, but its effects are complex and species-specific (Huntingford, 2004; Huntingford et al., 2006). For example, a breakdown of agonistic behaviours (i.e., social behaviour associated with aggression including threats, displays, retreats, placation and/or conciliation; Barrows, 2000) is commonly observed among nonsocial fish reared at high densities, and this can be explained by a trade-off between maintaining a social rank and effectively competing for resources (Ellis et al., 2002). Yet, captive bred fish can either be more or less aggressive than wild counterparts, depending on species (Huntingford, 2004). While crowding may increase aggression in some species, it can suppress it in others (Martins et al., 2012). Thus, understanding the molecular basis of such contrasting responses of fish to crowding is key for selecting individuals that adapt well to captivity during aquaculture intensification.

The advent of genomic tools offers new opportunities to decipher the molecular mechanisms underlying aggression and other complex behaviours affected by domestication (Robinson et al., 2008). These studies have shown that some behaviours are orchestrated by neurohormonal and gene expression regulatory networks which are largely conserved across vertebrates (Freudenberg,

Carreño Gutierrez, Post, Reif, & Norton, 2016; Goodson, 2005; O'Connell & Hofmann, 2012). Yet, one outstanding challenge is to document how the environment interacts with molecular pathways to shape variation in individual behaviours (Maruska & Fernald, 2014). For example, studies in model organisms indicate that the display of aggression depends on genotype × social environment interactions, and that these drive individual variation in agonistic behaviours (Gallardo-Pujol, Andrés-Pueyo, & Maydeu-Olivares, 2013; Rohde, Gaertner, Wards, Sørensen, & Mackay, 2017). In some fish, changes in social status can alter the expression of specific genes in the brain, and these can modify skin darkening (a proxy for social stress; Höglund, Balm, & Winberg, 2000), dopamine pathways (Weitekamp, Nguyen, & Hofmann, 2017) and complex behaviours (Fernald, 2012). This suggests that changes in social stress, brought about by crowding, might be reflected in changes in gene expression, which might in turn alter behaviour.

We therefore investigated how aggression, social stress and gene expression changed during aquaculture intensification in Nile tilapia (Oreochromis niloticus, Linnaeus, 1758), one of the oldest and most extensively farmed fish worldwide (FAO, 2016). As most cichlids, this species is structured into social hierarchies maintained by agonistic interactions in the wild (El-Sayed, 2006), but under aquaculture conditions rearing density has a profound effect on the morphology, physiology and behaviour of tilapia (Barcellos, Nicolaiewsky, De Souza, & Lulhier, 1999; El-Sayed, 2002; Fessehaye, Kabir, Bovenhuis, & Komen, 2006). In particular, high stocking density causes a shift from antagonistic (aggressive) to shoaling behaviour (Gall & Bakar, 1999), and more generally, from a proactive to a reactive stress-coping style (Champneys, Castaldo, Consuegra, & Garcia de Leaniz, 2018) that has implications for welfare and disease resistance (Ellison et al., 2018). Yet, the underlying molecular mechanisms of such a dramatic behavioural shift remain largely unknown. We compared patterns of gene expression at two contrasting densities in order to disentangle the molecular pathways responsible for the behavioural changes that accompany crowding during aquaculture intensification. Our hypothesis was that rearing density would modulate the frequency of agonistic interactions, and that individuals with different aggression levels and stress-coping styles would differ in the expression in the hypothalamus of key genes involved in the stress response.

# 2 | MATERIAL AND METHODS

## 2.1 | Origin and rearing of fish

A total of 360 mixed-sex three-week-old Nile tilapia (O. niloticus, silver strain) were sourced from a commercial supplier (Fishgen Ltd) which employs communal tank spawning, typically involving four sires and 12-15 dams per spawning tank. The fish were acclimatized at a recirculation fish facility for 3 days and were then randomly distributed into nine identical 20-L tanks (40L × 30W × 22H cm). Six tanks were stocked with 20 fry/tank (initial biomass = 0.21 g/L, final biomass = 25.5 g/L), and three tanks were stocked with 80 fry/tank (initial biomass = 0.86 g/L, final biomass = 95 g/L). These densities are
commercially relevant for tilapia farming and representative of low and high densities in recirculation aquaculture systems (Conte, 2004). The higher number of replicates used for the low density treatment (n = 6) compared to high density (n = 3) was motivated by the need to mark the same number of fish for testing (see below), but tank effects were explicitly taken into consideration in the statistical analysis.

After 10 weeks of rearing, 120 fish (10 from each of the six low density tanks and 20 fish from each of the three high density tanks) were individually marked with intraperitoneal PIT tags (7 × 1.35 mm, 30 mg, Loligo Systems), returned to their original tanks and allowed to recover for 10 days before the start of the behavioural screening. Fish were fed twice daily (Skretting), progressively reducing the ration from 20% to 5% body weight per day and increasing pellet size following commercial guidelines as per feed manufacturers' recommendations. Rearing conditions and water quality were maintained within the optimal range for the species (El-Sayed, 2006); temperature: 25-27.5°C; dissolved oxygen > 75%; photoperiod 12D:12L).

## 2.2 | Behavioural screening

We quantified the behaviour of tilapia towards their own mirror image using the mirror image stimulation (MIS) test (Balzarini, Taborsky, Wanner, Koch, & Frommen, 2014) to assess individual variation in aggression of fish reared at low and high densities. The lack of visual self-recognition in fish supports the use of MIS to quantify aggression under repeatable conditions, unaffected by chemical cues from conspecifics (Balzarini et al., 2014; Barreto et al., 2009). Fish were assessed individually in four experimental test tanks (60L × 30W × 30H cm) consisting of an acclimation area (15L × 30W cm), fitted with a remotely operated gate door, and a mirror at the opposite end of the tank. Three lines drawn on the bottom of the tanks delimited four equal zones at varying distances from the

mirror to help us assess activity levels (Figure 1). Tanks were separated by opaque walls to prevent test fish from seeing each other.

Fish were introduced singly in the acclimation zone, and after 10 min, the gate door was opened remotely and their behaviour was recorded for 20 min using two video cameras (A-Tech Sony EFFIO 580TVL CCD Outdoor Camera) mounted above and to the side. We measured two behaviours in each fish: (a) "activity," defined as the number of crosses between zones, and (b) "aggression," defined by the number of agonistic actions (nips and charges) directed towards their mirror image (Barreto et al., 2009). At the end of each test, the fish were returned to their original rearing tanks, and the test tanks were washed with 90% ethanol and rinsed with distilled water to prevent the build-up of stress hormones that might affect the behaviour of subsequent fish (Roberts, Taylor, & Garcia de Leaniz, 2011). After 10 days—to give fish time to recover from any stress associated with testing—all individually tagged fish used in the mirror test were humanely euthanized by an overdose of anaesthesia, sexed by visual inspection of the gonads and weighed. Their brains were rapidly dissected, weighed and stored in RNAlater at -20°C for subsequent analyses of gene expression.

## 2.3 | Assessment of skin brightness

To assess the extent of skin brightness, photographs of each fish were taken underwater during the acclimation period using a Canon EOS 400D Digital camera and a white background fitted with a colour standard (Classic Target-X-rite-Color Checker). Colour standardization and analysis were performed using GIMP 2.8.16 (Solomon, 2009) as per Clarke and Schlüter (2011). Greyscale filtered values (0-255) of R, G and B were measured along the fish flank, between the beginning and end of the dorsal fin, and compared to the background colour standard. From this, background-corrected average greyscale values

<figure><img src="image_4.png" /><figcaption>**FIGURE 1** Experimental set up for mirror image stimulation used to assess aggression in Nile tilapia</figcaption></figure>
were converted to luminance (brightness) in HSV space using the rg-btohsv function in the grDevices R base package (R Core Team, 2017).

## 2.4 | Transcriptomic analysis

For transcriptomic analysis, we chose for each density the six most aggressive and six least aggressive individuals on the basis of their MIS scores (n = 24) and selected only males to reduce unwanted variability resulting from sex differences in gene expression (Trainor & Hofmann, 2007; Zabegalov et al., 2019). The hypothalamus of each fish was detached under a dissecting microscope; we chose this brain region because previous studies had shown it is involved in the control of social behaviours and aggression in fish (Filby, Paull, Hickmore, & Tyler, 2010; Goodson, 2005). Total RNA was extracted using the AllPrep DNA/RNA Mini Kit (Qiagen) following manufacturer's instructions. The final product was eluted in 40 µl RNase-free water. RNA quality (quantity, purity and integrity) was checked using a NanoDrop NS-100 Spectrophotometer (NanoDrop Technologies) and an Agilent 2100 Bioanalyzer (Agilent Technologies). All RNA used for library construction was of high quality, having 260/230 and 260/230 ratios > 1.8 and RIN scores > 8. Individual cDNA libraries were prepared using Illumina TruSeq RNA sample preparation kit (1ug of total RNA; 8p.m. final concentration) and quantified using a Qubit Fluorometer (Invitrogen). The resulting 24 libraries were indexed, pooled and sequenced on an Illumina HiSeq2500 platform (2 lanes-2 × 126 bp).

Low-quality reads and Illumina Truseq adaptors were filtered out using Trimmomatic v 0.33 (Bolger, Lohse, & Usadel, 2014), excluding reads that were <36 bases long. After a quality check using FastQC v0.11.2 (Andrews, 2010), we mapped mRNA-seq reads using the Tophat 2.1.1 and Cufflinks 2.2.1 (Linux ×86-64) pipeline (Trapnell et al., 2012) to the available tilapia (O. niloticus) genome assembly (Orenil 1.0, accession number PRJNA5957, Brawand et al., 2014. Sequenced and assembled at the Broad Institute from a female Nile tilapia originating from a clonal line provided by D. Penman, Institute of Aquaculture, Stirling, UK). Between 15 and 20 million reads per sample were mapped to the genome, representing 87%-98% of all generated reads, and 70% of these were mapped to Ensembl annotated coding regions. Using uniquely mapped reads from Tophat, read counts per exon were obtained using the summarize overlaps function from the Genomic Alignments package (Lawrence et al., 2013), predefining gene models, grouping exons by gene for counting reads using Ensembl Tilapia's GTF file.

## 2.5 | Statistical analyses

All statistical analyses were conducted with R 3.3.3 (R Core Team, 2017). To model the effects of rearing density on fish mass, we used a linear mixed model approach (LMM) with rearing density and sex as fixed factors and tank identity as a random factor. To achieve model simplification, we started with a full model with all main effects and interactions and selected the model with the lowest AIC value via backward selection using the step function in the lmerTest package (Kuznetsova, Brockhoff, & Christensen, 2017) which was then

refitted via restricted maximum likelihood or as a linear model when the random component was not significant when compared to the fixed effects only model by the log likelihood ratio (LRT). The same approach was used to assess the effects of rearing density on brain weight and skin brightness, using rearing density, body mass and sex as predictors, and tank identity as a random effect.

To model the effect of rearing density on activity and agonistic interactions (both count data), we employed a generalized linear mixed model (GLMM) with a Poisson log-link while statistically controlling for the effects of body mass and sex. We corrected for overdispersion using fish identity nested within tank to generate a random effect with one level per observation using the lme4 R package (Bates, Maechler, Bolker, & Walker, 2014), as suggested by Harrison (2014). To achieve model simplification, we started with a full model with all main effects and selected the model with the lowest AICc value using the dredge function in the MuMln R package (Barton, 2015). We used the ANOVA command to compare the null model with only a random effect grouping structure and the most plausible mixed model with both random and fixed effects.

Read counts per gene were analysed with the DESeq2 R package, based on the negative binomial distribution (Love, Huber, & Anders, 2014). Counts were prefiltered with a threshold of >20 reads per sample in more than one sample via the independent filtering component of DeSeq2, using default parameters to optimize model sensitivity and improve computational speed (Love et al., 2014). As the number of attacks and extent of body brightness differed significantly between densities, these were also included as predictors in the analysis. Initial inspection of the data by PCA and hierarchical clustering indicated the existence of two extreme values (one from each density; Figure S1a). Although not a single standard exists for dealing with outliers in RNA-seq analysis (Conesa et al., 2016), as these two outliers were outside two standard deviations of the median of the first two principal components they were excluded from analysis (Ellis et al., 2013). Inspection of variance inflation factors (VIF) indicated low multicollinearity for both aggressive (VIF on PC1 = 2.7, 3.2 and 5.9) and nonaggressive fish (VIF on PC1 = 1.4, 1.6 and 1.1).

After removing the two outliers and filtering low counts, 18,963 Ensembl genes were used for exploratory and differential expression analysis. PCA inspection showed no obvious clustering by density, skin brightness or aggression intensity after outlier removal (Figure S1b), but revealed a much higher variation for nonaggressive fish than for aggressive ones (Figure S2). Different intragroup variability is a well-known problem in RNA-seq data analysis, as it can inflate per gene dispersion estimates and reduce statistical power for detecting differentially expressed genes, which can go undetected (Landau & Liu, 2013). Therefore, as recommended for such cases within DESeq2 analysis (Love et al., 2014), we modelled gene expression separately for the two heterogenous groups (most aggressive and least aggressive fish) using density, number of attacks and skin brightness (and their interactions) as predictors. For comparisons, we also include the pooled results using aggression as a categorical variable in a 2-way ANOVA (i.e., gene expression ~ density + aggression group + density: aggression group) as supplementary
material (Tables S1-S3). We used a false discovery rate (FDR) of *p* corrected values < 0.05 for differential gene expression (Ellison et al., 2018). Functional annotation and enrichment analysis, based on *Gene Ontology* terms, was performed using DAVID 6.8 (Huang et al., 2007).

## 3 | RESULTS

### 3.1 | Body mass

Body mass was unaffected by rearing density (*p* = 0.204) or the interaction between density and sex (*p* = 0.827), but males were 27% larger than females (LM estimate for males = 8.30, *SE* = 2.16, *t*<sub>102</sub> = 3.84, *p* < 0.001). The random component (tank identity) was not significant (LRT = 2.189, *df* = 1, *p* = 0.139) and was dropped from the model.

### 3.2 | Brain weight

Brain weight increased with body mass (*p* < 0.001) and was proportionally larger at high density than at low density (LM estimate for low density = -2.49 × e-03, *SE* = 9.82 × e-04, *t*<sub>94</sub> = -2.54, *p* = 0.013, Figure 2). There were no sex differences in relative brain size (*p* = 0.722) and none of the interactions were significant (*p* > 0.5); variation among tanks (random component) was also not significant (LRT = 1.802, *df* = 1, *p* = 0.179).

### 3.3 | Skin brightness

Tilapia reared at high density were significantly darker (i.e., had lower values of brightness) than fish at low density (LM estimate for low density = 0.072, *SE* = 0.017, *t*<sub>105</sub> = 4.287, *p* < 0.001; Figure 3a). Skin brightness was unaffected by body mass (*p* = 0.609), sex (*p* = 0.454) or any of the interactions (*p* > 0.3). As for brain size, the random component (tank identity) was not significant (LRT = 2.274e-13, *df* = 1, *p* = 1.00) and was dropped from the model.

<figure><img src="image_5.png" /><figcaption>**FIGURE 2** Relationship between body weight and brain weight in Nile tilapia reared at low and high density (density effect *t*<sub>94</sub> = -2.54, *p* = 0.013)</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**FIGURE 3** Effect of rearing density on variation of (a) skin brightness (0 = black, 1 = white), (b) activity (no. of zone crosses in the test arena) and (c) aggression (no. of agonistic interactions against a mirror image)</figcaption></figure>

### 3.4 | Activity

Activity levels, measured as number of crosses between zones, was unrelated to body mass (*p* = 0.952) or sex (*p* = 0.967), but was strongly influenced by rearing density (Figure 3b). Fish reared at low density were significantly more active (mean number of crosses = 19.55, *SE* = 3.33) than those reared at high density (mean number of crosses = 7.43, *SE* = 1.40), while statistically controlling for tanks effects (GLMM estimate for low density = 0.850, *SE* = 0.299, *z* = 2.841, *p* = 0.004).

### 3.5 | Aggression

Aggression, measured as number of agonistic interactions against the mirror, was unrelated to body mass (*p* = 0.750) or sex (*p* = 0.514), but fish reared at low density were significantly more aggressive (mean number of interactions = 28.63, *SE* = 5.93) than those reared at high density (mean number of interactions = 7.11, *SE* = 2.09), while statistically controlling for tanks effects (GLMM estimate for low density = 3.103, *SE* = 0.887, *z* = 3.498, *p* < 0.001; Figure 3c).
## 3.6 | Transcriptional analysis (RNA-seq)

We examined the effects of rearing density, aggression intensity (number of attacks) and skin darkening (and their interactions) as predictors of gene expression in the brain of the most and least aggressive individuals at each density. After correcting for multiple testing, we found 25 responsive genes among the most aggressive fish and 41 genes among the least aggressive fish (Table 1), with only one gene (fosab) being responsive among both groups. Most of the variation in gene expression among aggressive fish was accounted for by density (eight genes) and the interaction between density and aggression intensity (19 genes), whereas for nonaggressive fish most of the variation in gene expression was associated with variation in skin brightness (38 genes).

The eight genes that were differentially expressed between high and low density among aggressive fish included somatostatin (sst1), a gene that encodes for a hormone involved in numerous cellular process, jarid2b—involved in transcription and histone methylation regulation, krt5—involved in the immune response and bag6l—in-volved in apoptosis and cellular stress. The expression of 19 genes showed significant aggression × density interactions, including genes involved in organismal response to stress such as v-fos FBJ murine osteosarcoma viral oncogene homolog Ab (fosab; *q* = 0.010; Figure 4a), proopiomelanocortin (pomca)—a gene that encodes a peptide hormone (ACTH) involved in the release of cortisol (*q* = 0.048; Figure 4b), krt5 (skin stress; *q* < 0.01; Figure 4c), nnt (mitochondrial antioxidant defence; *q* < 0.01; Figure 4d), irg1 (inflammation process and response to bacterial infections; *q* < 0.01; Figure 4e) and bag6l (apoptosis and cellular stress; *q* < 0.01; Figure 4f). Other significant interactions in gene expression included density × skin brightness (three genes) and number of attacks × skin brightness (five genes, Table 1).

Among nonaggressive fish, 38 genes were significantly associated with skin brightness, being in most cases upregulated among dark fish and downregulated among pale ones (Table 1). These included genes involved in the stress response and ion homeostasis (oxt, avp), neuropeptides involved in feeding behaviour (grp, hcrt), and several genes that coded for hormone receptors, including the thyrotropin-releasing hormone (trh) and neurotensin (nts) involved in dopamine signalling, as well as homebox transcription factors (Figure 5). Three genes were differentially expressed between rearing densities, including prolactin (prl) and fosab which were also affected by the interaction between density and brightness (*q* = 0.078), as well as one uncharacterized gene (Table 1).

## 4 | DISCUSSION

Fish domestication is thought to be constrained by the capacity of individuals to adapt to high densities in captivity (Ashley, 2007; Huntingford & Adams, 2005), and our study shows that one consequence of crowding in Nile tilapia is the inhibition of aggression, a process that we found was associated with the expression in the

hypothalamus of stress-related genes. Stress has a profound effect on the hypothalamic–pituitary–interrenal (HPI) axis and the neuroendocrine response of zebrafish (Pavlidis, Sundvik, Chen, & Panula, 2011; Pavlidis, Theodoridi, & Tsalafouta, 2015), but the transcriptional response of farmed fish to the increase in density that accompanies aquaculture intensification is not well understood. We employed a mirror test to quantify aggression and skin darkening (a proxy for stress) at low and high density and screened the most and least aggressive individuals at each density for differential gene expression in the hypothalamus in order to examine the transcriptional response to crowding. We found that tilapia reared at high density were darker, less active and less aggressive than fish reared at low density, and that these differences were associated with transcriptional differences in the brain.

Rearing density affects the nature and strength of social interactions in many fish (Ashley, 2007) and in species such as tilapia, where social hierarchies are maintained by agonistic interactions, aggressive behaviour can be used as a predictor of stress (Barreto et al., 2009). Our results suggest the existence of two types of stress in tilapia: social stress caused by the formation of social hierarchies maintained by agonistic interactions at low density and chronic stress caused by crowding at high density.

Although only a relatively small number of genes were differentially expressed in our study, this is common in studies of gene expression in the brain and can be explained by the wide range of specialized neuronal cell types present in the brain, and in the hypothalamus in particular (Machluf, Gutnick, & Levkowitz, 2011), and the tight homeostatic balance of the nervous system (Aubin-Horth, Landry, Letcher, & Hofmann, 2005; Filby et al., 2010). Our results are also consistent with those reported for zebrafish, where 70 genes were differentially expressed in the brain of fighting and nonfighting (isolated) individuals (Malki et al., 2016), seven of which were homologous to those differentially expressed in aggressive and nonaggressive mice.

Of the various genes that were differentially expressed between high and low density among aggressive individuals, somatostatin (sst1) has previously been linked to aggression in fish (Filby et al., 2010; Trainor & Hofmann, 2006), but the link is complex and appears to be species-specific. Thus, while upregulation of sst was found to inhibit aggression in the African cichlid *Astatotilapia burtoni* (Trainor & Hofmann, 2006), a species that switches from territorial to nonterritorial depending on the social environment, in zebrafish, a typically shoaling fish, the opposite appears to be true (Filby et al., 2010). In our study, sst1 was upregulated among fish reared at high density, which were significantly less aggressive than those reared at low density, suggesting an inhibitory role for somatostatin on aggression in Nile tilapia, as seen in other cichlids (Hofmann & Fernald, 2000; Trainor & Hofmann, 2006, 2007). Thus, sst1 appears to be one of the key genes regulating the different response of species to crowding. The two other density-dependent genes differently expressed among aggressive fish, *pomca* and *fosab*, are both involved in the stress response (Eissa & Wang, 2016), and their expression in our study depended on the
<table><caption>**TABLE 1** ENSMBL genes significantly influenced by density (low vs. high), aggression intensity (number of agonistic interaction), skin brightness or their interactions among the six most aggressive and six least aggressive fish at each density</caption><thead><tr><th>Gene ID</th><th>log2 fold change</th><th>SE</th><th>stat</th><th>padj</th><th>Gene name</th><th>Gene description</th><th>GOTerm (BP)</th></tr></thead><tbody><tr><td colspan="8">Aggressive fish</td></tr><tr><td colspan="8">Influenced by density (n = 8)</td></tr><tr><td>ENSONIG00000009137</td><td>-32.47</td><td>6.21</td><td>-5.23</td><td>0.00</td><td>krt5</td><td>keratin 5</td><td></td></tr><tr><td>ENSONIG00000001057</td><td>-11.64</td><td>2.70</td><td>-4.31</td><td>0.03</td><td>sst1.1</td><td>somatostatin 1, tandem duplicate 1</td><td>0010469</td></tr><tr><td>ENSONIG00000002652</td><td>-9.19</td><td>2.17</td><td>-4.23</td><td>0.04</td><td>jarid2b</td><td>junonji, AT-rich interactive domain 2b</td><td>0006355;0007275;0016569;0031061<br>0051574</td></tr><tr><td>ENSONIG00000010340</td><td>8.47</td><td>1.92</td><td>4.40</td><td>0.03</td><td>arhgef1a</td><td>Rho guanine nucleotide exchange factor (GEF) 1a</td><td>0035556;0035023</td></tr><tr><td>ENSONIG00000003670</td><td>14.19</td><td>2.58</td><td>5.50</td><td>0.00</td><td>nnt</td><td>nicotinamide nucleotide transhydrogenase</td><td>0015992;0055114</td></tr><tr><td>ENSONIG00000019114</td><td>19.53</td><td>4.51</td><td>4.33</td><td>0.03</td><td>MFAP4 (one of many)</td><td>Microfibril-associated protein 4</td><td></td></tr><tr><td>ENSONIG00000020827</td><td>25.72</td><td>5.67</td><td>4.54</td><td>0.02</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG00000015953</td><td>35.47</td><td>4.91</td><td>7.22</td><td>0.00</td><td>bag6l</td><td></td><td></td></tr><tr><td colspan="8">Influenced by number of attacks (n = 2)</td></tr><tr><td>ENSONIG0000000379</td><td>0.13</td><td>0.03</td><td>4.58</td><td>0.04</td><td>minpp1a</td><td>multiple inositol-polyphosphate phosphatase 1a</td><td></td></tr><tr><td>ENSONIG00000015953</td><td>0.41</td><td>0.06</td><td>6.67</td><td>0.00</td><td>bag6l</td><td>BCL2-associated athanogene 6 like</td><td></td></tr><tr><td colspan="8">Influenced by brightness (n = 1)</td></tr><tr><td>ENSONIG00000015953</td><td>38.16</td><td>5.95</td><td>6.41</td><td>0.00</td><td>bag6l</td><td>BCL2-associated athanogene 6, like</td><td></td></tr><tr><td colspan="8">Influenced by density × no. of attacks (n = 19)</td></tr><tr><td>ENSONIG00000003670</td><td>0.08</td><td>0.01</td><td>7.17</td><td>0.00</td><td>Nnt</td><td>nicotinamide nucleotide transhydrogenase</td><td>0015992;0055114</td></tr><tr><td>ENSONIG00000004830</td><td>-0.09</td><td>0.02</td><td>-4.11</td><td>0.04</td><td>si:ch73-335121.4</td><td>si:ch73-335121.4</td><td></td></tr><tr><td>ENSONIG00000006273</td><td>-0.05</td><td>0.01</td><td>-5.01</td><td>0.00</td><td>tuft1a</td><td>tuftelin 1a</td><td></td></tr><tr><td>ENSONIG00000006744</td><td>-0.15</td><td>0.03</td><td>-5.61</td><td>0.00</td><td>zgc:162944</td><td>zgc:162944</td><td></td></tr><tr><td>ENSONIG00000006747</td><td>-0.15</td><td>0.03</td><td>-5.68</td><td>0.00</td><td>irg1</td><td>immunoresponsive 1 homolog (mouse)</td><td>0006954;0002376;0045087</td></tr><tr><td>ENSONIG00000007532</td><td>-0.06</td><td>0.01</td><td>-5.45</td><td>0.00</td><td>zgc:122979</td><td>zgc:122979</td><td>0006457</td></tr><tr><td>ENSONIG00000009137</td><td>-0.15</td><td>0.03</td><td>-4.86</td><td>0.00</td><td>krt5</td><td>keratin 5</td><td></td></tr><tr><td>ENSONIG00000009460</td><td>-0.16</td><td>0.04</td><td>-4.09</td><td>0.04</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG00000009551</td><td>-0.06</td><td>0.01</td><td>-4.10</td><td>0.04</td><td>SPIDR</td><td>scaffolding protein involved in DNA repair</td><td>0000724;0006974</td></tr><tr><td>ENSONIG00000010366</td><td>-0.05</td><td>0.01</td><td>-4.89</td><td>0.00</td><td>Txnipa</td><td>thioredoxin-interacting protein a</td><td></td></tr><tr><td>ENSONIG00000012279</td><td>-0.11</td><td>0.03</td><td>-4.18</td><td>0.03</td><td>zmp:0000000801</td><td>zmp:0000000801</td><td>0007186</td></tr><tr><td>ENSONIG00000015140</td><td>-0.13</td><td>0.03</td><td>-4.85</td><td>0.00</td><td>-</td><td>-</td><td></td></tr></tbody></table>

(Continues)
<table><caption>TABLE 1 (Continued)</caption><thead><tr><th>Gene ID</th><th>log2 fold change</th><th>SE</th><th>stat</th><th>padj</th><th>Gene name</th><th>Gene description</th><th>GOTerm (BP)</th></tr></thead><tbody><tr><td>ENSONIG00000015264</td><td>-0.06</td><td>0.01</td><td>-4.49</td><td>0.01</td><td>Fosab</td><td>v-fos FBJ murine osteosarcoma viral oncogene homolog Ab</td><td>0033555;0006355</td></tr><tr><td>ENSONIG00000015953</td><td>0.13</td><td>0.02</td><td>6.50</td><td>0.00</td><td>bag6l</td><td>BCL2-associated athanogene 6, like</td><td></td></tr><tr><td>ENSONIG00000016987</td><td>-0.06</td><td>0.01</td><td>-5.62</td><td>0.00</td><td>ddit4</td><td>DNA-damage-inducible transcript 4</td><td>0009968</td></tr><tr><td>ENSONIG00000016884</td><td>-0.05</td><td>0.01</td><td>-4.98</td><td>0.00</td><td>nt5dc2</td><td>5'-nucleotidase domain containing 2</td><td></td></tr><tr><td>ENSONIG00000018737</td><td>0.05</td><td>0.01</td><td>4.05</td><td>0.05</td><td>Pomca</td><td>proopiomelanocortin a</td><td>0007218;0032400;0033555</td></tr><tr><td>ENSONIG00000019777</td><td>-0.09</td><td>0.02</td><td>-4.21</td><td>0.03</td><td>hbae5</td><td>haemoglobin, alpha embryonic 5</td><td>0015671</td></tr><tr><td>ENSONIG00000021429</td><td>-0.08</td><td>0.02</td><td>-4.62</td><td>0.01</td><td>Cebpd</td><td>CCAAT/enhancer-binding protein (C/EBP), delta</td><td>0006355</td></tr><tr><td colspan="8">Influenced by density × brightness (n = 3)</td></tr><tr><td>ENSONIG0000003670</td><td>19.59</td><td>3.48</td><td>5.63</td><td>0.00</td><td>Nnt</td><td>nicotinamide nucleotide transhydrogenase</td><td>0015992;0055114</td></tr><tr><td>ENSONIG00000010340</td><td>12.50</td><td>2.59</td><td>4.83</td><td>0.01</td><td>arhgef1a</td><td>Rho guanine nucleotide exchange factor (GEF) 1a</td><td>0035556;0035023</td></tr><tr><td>ENSONIG00000015953</td><td>35.63</td><td>6.16</td><td>5.78</td><td>0.00</td><td>bag6l</td><td>BCL2-associated athanogene 6, like</td><td></td></tr><tr><td colspan="8">Influenced by no. of attacks × brightness (n = 5)</td></tr><tr><td>ENSONIG0000006744</td><td>0.69</td><td>0.14</td><td>4.78</td><td>0.01</td><td>zgc:162944</td><td>zgc:162944</td><td></td></tr><tr><td>ENSONIG0000006747</td><td>0.70</td><td>0.14</td><td>4.85</td><td>0.01</td><td>irg1</td><td>immunorespective 1 homolog (mouse)</td><td>0006954;0002376;0045087</td></tr><tr><td>ENSONIG0000009137</td><td>0.69</td><td>0.15</td><td>4.45</td><td>0.03</td><td>krt5</td><td>keratin</td><td></td></tr><tr><td>ENSONIG00000015140</td><td>0.71</td><td>0.14</td><td>4.93</td><td>0.01</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG00000015953</td><td>-0.76</td><td>0.11</td><td>-6.84</td><td>0.00</td><td>bag6l</td><td>BCL2-associated athanogene 6, like</td><td></td></tr><tr><td colspan="8">Nonaggressive fish</td></tr><tr><td colspan="8">Influenced by density (n = 3)</td></tr><tr><td>ENSONIG0000006467</td><td>-24.58</td><td>3.94</td><td>-6.25</td><td>0.00</td><td></td><td></td><td></td></tr><tr><td>ENSONIG0000006469</td><td>-27.49</td><td>5.91</td><td>-4.65</td><td>0.02</td><td>prl</td><td>prolactin</td><td>0010469</td></tr><tr><td>ENSONIG00000015264</td><td>7.48</td><td>1.56</td><td>4.80</td><td>0.01</td><td>fosab</td><td>v-fos FBJ murine osteosarcoma viral oncogene homolog Ab</td><td>0033555;0006355</td></tr><tr><td colspan="8">Influenced by brightness (n = 38)</td></tr><tr><td>ENSONIG0000002382</td><td>-37.14</td><td>9.17</td><td>-4.05</td><td>0.04</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG00000015235</td><td>-24.84</td><td>3.79</td><td>-6.55</td><td>0.00</td><td>oct</td><td>oxytocin</td><td>0050801</td></tr><tr><td>ENSONIG00000015218</td><td>-18.10</td><td>4.39</td><td>-4.12</td><td>0.03</td><td>avp</td><td>arginine vasopressin</td><td>0007165;0010469</td></tr><tr><td>ENSONIG0000003962</td><td>-17.98</td><td>3.15</td><td>-5.70</td><td>0.00</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG0000002871</td><td>-16.93</td><td>3.24</td><td>-5.23</td><td>0.00</td><td>prdm12b</td><td>PR domain containing 12b</td><td>0006351;0006355; 0021521</td></tr><tr><td>ENSONIG0000007689</td><td>-13.67</td><td>2.75</td><td>-4.97</td><td>0.00</td><td>grp</td><td>gastrin-releasing peptide(grp)</td><td>0007218</td></tr><tr><td>ENSONIG00000020662</td><td>-12.86</td><td>2.94</td><td>-4.37</td><td>0.01</td><td>-</td><td>-</td><td></td></tr></tbody></table>

(Continues)
<table><caption>TABLE 1 (Continued)</caption><thead><tr><th>Gene ID</th><th>log2 fold change</th><th>SE</th><th>stat</th><th>padj</th><th>Gene name</th><th>Gene description</th><th>GOTerm (BP)</th></tr></thead><tbody><tr><td>ENSONIG00000005071</td><td>-12.73</td><td>2.92</td><td>-4.37</td><td>0.01</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG00000002896</td><td>-12.49</td><td>3.13</td><td>-3.99</td><td>0.04</td><td>sp8b</td><td>sp8 transcription factor b</td><td>0006351;0006355;0030326</td></tr><tr><td>ENSONIG00000005998</td><td>-12.16</td><td>3.08</td><td>-3.95</td><td>0.04</td><td>sytl3</td><td>synaptotagmin-like 3</td><td>0006886;0006887;0006906</td></tr><tr><td>ENSONIG00000004333</td><td>-10.17</td><td>1.93</td><td>-5.27</td><td>0.00</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG00000001128</td><td>-9.62</td><td>1.96</td><td>-4.91</td><td>0.00</td><td>hcrt</td><td>hypocretin (orexin) neuropeptide precursor</td><td>0007218;0007631;0033555</td></tr><tr><td>ENSONIG00000009224</td><td>-9.23</td><td>1.88</td><td>-4.91</td><td>0.00</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG000000010719</td><td>-8.92</td><td>2.23</td><td>-4.00</td><td>0.04</td><td>si:ch211-236k19.2</td><td>si:ch211-236k19.2</td><td>0006355</td></tr><tr><td>ENSONIG00000004631</td><td>-8.80</td><td>2.11</td><td>-4.18</td><td>0.03</td><td>col8a2</td><td>collagen, type VIII, alpha 2</td><td></td></tr><tr><td>ENSONIG000000014386</td><td>-7.94</td><td>1.46</td><td>-5.44</td><td>0.00</td><td>arxa</td><td>aristaless-related homeobox a</td><td>0003322;0006351;0006355</td></tr><tr><td>ENSONIG000000019324</td><td>-7.58</td><td>1.91</td><td>-3.96</td><td>0.04</td><td>tprg1</td><td>tumour protein p63 regulated 1</td><td></td></tr><tr><td>ENSONIG00000003653</td><td>-7.40</td><td>1.53</td><td>-4.84</td><td>0.00</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG00000001403</td><td>-7.31</td><td>1.86</td><td>-3.93</td><td>0.04</td><td>meis2a</td><td>Meis homeobox 2a</td><td>0006355;0048703;</td></tr><tr><td>ENSONIG00000003843</td><td>-5.75</td><td>1.08</td><td>-5.31</td><td>0.00</td><td>nts</td><td>neurotensin</td><td>0010469</td></tr><tr><td>ENSONIG000000017229</td><td>-4.67</td><td>1.18</td><td>-3.95</td><td>0.04</td><td>thr</td><td>thyrotropin-releasing hormone</td><td>0001692;0007165;0009755</td></tr><tr><td>ENSONIG000000014561</td><td>-4.50</td><td>1.01</td><td>-4.48</td><td>0.01</td><td>tacr3a</td><td>neuromedin-K receptor-like</td><td>0007165;0007186</td></tr><tr><td>ENSONIG000000019723</td><td>-4.20</td><td>1.04</td><td>-4.05</td><td>0.04</td><td>s100a10a</td><td>S100 calcium-binding protein A10a</td><td>0042127</td></tr><tr><td>ENSONIG000000009148</td><td>-3.64</td><td>0.86</td><td>-4.22</td><td>0.02</td><td>vax1</td><td>ventral anterior homeobox 1</td><td>0006351;0006355;0007275</td></tr><tr><td>ENSONIG000000016412</td><td>-2.88</td><td>0.65</td><td>-4.45</td><td>0.01</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG000000007099</td><td>-2.73</td><td>0.66</td><td>-4.13</td><td>0.03</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG000000014947</td><td>-2.68</td><td>0.67</td><td>-3.99</td><td>0.04</td><td>CACNA2D1</td><td>calcium voltage-gated channel auxiliary subunit alpha2delta 1</td><td>1904646;0070588</td></tr><tr><td>ENSONIG000000016069</td><td>-2.59</td><td>0.52</td><td>-5.02</td><td>0.00</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG000000008322</td><td>-2.58</td><td>0.64</td><td>-4.05</td><td>0.04</td><td>dhrs13a.3</td><td>dehydrogenase/reductase (SDR family) member 13a, duplicate 3</td><td></td></tr><tr><td>ENSONIG00000000606</td><td>2.10</td><td>0.54</td><td>3.86</td><td>0.05</td><td>EIF4A1A</td><td>eukaryotic translation initiation factor 4A1A</td><td>0006412;0006446;0016926</td></tr><tr><td>ENSONIG000000014333</td><td>2.27</td><td>0.56</td><td>4.03</td><td>0.04</td><td>syn2b</td><td>synapsin IIb</td><td>0007269</td></tr><tr><td>ENSONIG000000006168</td><td>2.69</td><td>0.67</td><td>3.99</td><td>0.04</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG000000010283</td><td>3.36</td><td>0.77</td><td>4.39</td><td>0.01</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG000000016676</td><td>4.54</td><td>1.14</td><td>3.98</td><td>0.04</td><td>-</td><td>-</td><td></td></tr><tr><td>ENSONIG000000009325</td><td>4.69</td><td>1.20</td><td>3.91</td><td>0.05</td><td>CYR61</td><td>cysteine-rich angiogenic inducer 61</td><td>0001558;0007155;0007165;0007267;0060548</td></tr><tr><td>Influenced by density × brightness (n = 1)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ENSONIG00000006467</td><td>-38.04</td><td>6.79</td><td>-5.60</td><td>0.00</td><td>-</td><td>-</td><td></td></tr></tbody></table>

Note: Positive log2 fold change values denote downregulation at high density (HD) compared to low density (LD), while for continuous predictors they represent rate of unit change.
<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

**FIGURE 4** Relationship between aggression intensity (no. of agonistic interactions) and gene expression (DESeq2 normalized counts) in the six most aggressive fish at each density. Shown are examples of six responsive genes (out of 19) where there were significant interactions between aggression intensity and rearing density ($q < 0.05$; Table 1), suggesting that agonistic behaviour is regulated in a density-dependent manner: (a) v-fos FBJ murine osteosarcoma viral oncogene homolog Ab (fosab), (b) proopiomelanocortin a (pomca), (c) krt5, (d) nnt, (e) irg1 and (f) bag6l

Density Low High

interaction between rearing density and aggression. We found a negative association between *pomca* expression and number of attacks at high density, but a positive association at low density. As *pomca* encodes a preproprotein whose proteolytic products include adrenocorticotropin (ACTH), which stimulates cortisol secretion (Bornstein & Chrousos, 1999), this suggests that elevated levels of aggression may result in increased cortisol production. Indeed, a recent study in a group living cichlid has shown that aggressive behaviour is associated with higher production of cortisol (Culbert, Gilmour, & Balshine, 2018). Another proteolytic product of *pomca* includes the alpha-melanocyte-stimulating hormone, α-MSH, which disperses pigment-containing melanosomes in the pigment cells, thereby making fish look darker (Kobayashi, Mizusawa, Chiba, Tagawa, & Takahashi, 2012; Kobayashi, Mizusawa, Saito, & Takahashi, 2012). We found that body darkening (i.e., lower brightness) was significantly associated with high density, adding support to the idea that body darkening is a good proxy for chronic stress in Nile tilapia (Champneys et al., 2018), as seen in other fish (Höglund et al., 2000). In zebrafish, aggressive individuals are also darker, while those that exhibit fear become paler (Gerlai, Lahav, Guo, & Rosenthal, 2000), which serves to highlight the relationship that exists between social stress, skin darkening and aggression.

An interaction between rearing density and aggression was also found with respect to the expression of fosab among aggressive fish. We found that fosab expression was positively associated with aggression intensity at high density, but negatively associated at low density. Fosab encodes c-Fos, whose upregulation has previously been associated to both stress (Kovács, 1998; Pavlidis et al., 2015) and aggression (Davis & Marler, 2004; Haller, Tóth, Halasz, & De Boer, 2006; Malki et al., 2016), and which also modulates neural plasticity underlining behavioural flexibility (Oliveira, 2012). This suggests that expression of fosab and aggression in tilapia are context dependent and are driven by social stress caused by agonistic interactions at low density and by crowding and chronic stress at high density. Aggression is influenced by genotype × environment interactions, notably social stress, in mammals and fish alike (Zabegalov et al., 2019), and many genes that have been associated with aggression and stress in mammalian studies respond in a similar manner in zebrafish (Freudenberg et al., 2016). These include oxt, avp, *hcrt* and *sst1* which are upregulated in dominant zebrafish, and which we found in our study to be influenced by skin darkening (and hence social stress) and density (*sst1*) in Nile tilapia.

Crowding has previously been found to have a darkening effect on fish, usually associated with increased chronic stress (Brown &
<figure><img src="image_13.png" /><figcaption>**FIGURE 5** Relationship between skin brightness (a proxy for stress) and gene expression (DESeq2 normalized counts) in the six least aggressive fish at each density. Shown are examples of five responsive genes (out of 38) where there was a significant negative relationship between skin brightness and gene expression ($q < 0.05$; Table 1) independently of rearing density, suggesting a common stress response: (a) *grp*, (b) *hcrt*, (c) oxytocin (*oxt*), (d) arginine vasopressin (*avp*) and (e) thyrotropin-releasing hormone (*trh*)</figcaption></figure>

Shahidi, 1997; Van der Salm, Martinez, Flik, & Wendelaar Bonga, 2004; Zeng et al., 2010). Among nonaggressive fish, we found a negative association between brightness and the expression in the hypothalamus of several important neuropeptides, hormones and hormone receptors, including the gastrin-releasing peptide (*grp*), hypocretin/orexin neuropeptide precursor (*hcrt*), oxytocin (*oxt*), vasopressin (*avp*), thyrotropin-releasing hormone (*trh*), parathyroid hormone 1 receptor (*pth1ra*) and the growth hormone-releasing hormone receptor (*ghrhrl*). All these genes have previously been reported to be upregulated in stressed or subordinate fish (Balment, Lu, Weybourne, & Warne, 2006; Bernier, 2006; Jezova, Skultetyova, Tokarev, Bakos, & Vigas, 1995; Pavlidis et al., 2011, 2015). The expression of prolactin (*prl*) and *fosab* also depended on the interaction between density and body darkening, though the effect was not strong. Prolactin is involved in the regulation of body pigmentation (Leclercq, Taylor, & Migaud, 2010; Oshima, Makino, Iwamuro, & Bern, 1996), and a common response to stress is an increase in prolactin secretion and mRNA levels (Auperin, Baroiller, Ricordel, Fostier, & Prunet, 1997; Avella, Schreck, & Prunet, 1991; Pavlidis et al., 2015), and hence darkening. Upregulation of *fosab* has also been associated with an increase in stress and skin darkening in other species (Kovács,

1998; Pavlidis et al., 2015), suggesting that both prolactin and *fosab* are conserved indicators of chronic stress caused by crowding.

The response of fish to stress depends not only on the type of stressor but also on the coping style that characterizes an individual (Koolhaas, De Boer, Coppens, & Buwalda, 2010; Koolhaas et al., 1999). Thus, the response along the HPI axis can be viewed as a dynamic multivariate system, where the effect of each element depends on the intensity of the others (Pavlidis et al., 2011, 2015). In our study, and with one exception (*fosab*), aggressive and nonaggressive fish expressed different genes in relation to rearing density, and among aggressive fish there were significant density-dependent interactions between aggression intensity and gene expression, suggesting that different neuroendocrine pathways may be involved. The existence of significant interactions between rearing density and aggression, and between density and body darkening, serves to highlight the fact that the regulation of aggression not only varies among individuals, but is also context dependent and shaped by social stress. On the other hand, aggression is one of the most repeatable and heritable fish behaviours (Bell, Hankison, & Laskowski, 2009; Way, Ruhl, Snekser, Kiesel, & McRobert, 2015), which might facilitate the selective breeding of nonaggressive individuals with
high tolerance to stress under crowded conditions. For example, several serotonergic genes and genes related to dopamine signalling and to c-Fos found in our study are also differentially expressed between aggressive and nonaggressive in zebrafish (Malki et al., 2016; Zabegalov et al., 2019), suggesting that these could be used as bio-markers for aggression. In addition, there is mounting evidence for an epigenetic regulation of aggression in humans and rodents (Waltes, Chiocchetti, & Freitag, 2016), which given the conserved nature of many cellular networks implicated in aggression across taxa (Malki et al., 2016), might also make it possible to use epigenetic conditioning to reduce aggression and stress under aquaculture intensification, as shown recently for disease resistance (Uren Webster et al., 2018).

Beyond the observed transcriptional and behavioural changes associated with crowding, tilapia reared at low density had smaller brains (relative to their body size) than fish reared at high density. In mammalian systems, chronic stress can impact on brain structures through the production of glucocorticoids which can in turn affect cognition and social interactions (Lupien, McEwen, Gunnar, & Heim, 2009), and it is possible that stress has a similar effect on the brain size of fishes. Fish brain size can evolve very rapidly in response to environmental conditions, reflecting a trade-off between enhanced cognition and reproductive performance (Kotrschal et al., 2013). For example, fish living in the wild develop larger brains compared to hatchery-reared conspecifics living in more impoverished environments (Kihslinger, Lema, & Nevitt, 2006), as do fish living in high predation areas compared to habitats with fewer predators (Kotrschal, Deacon, Magurran, & Kolm, 2017). This suggests that the smaller brain of fish reared at low density may have been the consequence of social stress caused by elevated aggression, and perhaps more generally, a trait associated with a proactive (rather than a reactive) stress-coping style (Champneys et al., 2018).

## 5 | CONCLUSIONS

Our study indicates the crowding inhibits aggressive behaviour in Nile tilapia and results in changes in the expression of stress-related genes that accompany the shift from social hierarchies maintained by agonistic interactions at low density, to shoaling at high density. This suggests that aquaculture can substantially alter the aggression level and stress response of Nile tilapia. Given that loss of aggression and stress tolerance are two of the defining features of animal domestication (Belyaev, 1969; Jensen, 2014), and that our study shows that these were associated with differential gene expression in Nile tilapia, it might be possible to selectively breed fish that perform well under crowded conditions under aquaculture intensification. Some of the genes identified here, such as *sstl* and *fosab*, will be key candidates for this, as they seem to play an important role in the responses of different species to crowding stress.

## ACKNOWLEDGEMENTS

This work was partially funded by NRN-LCEE AquaWales and ERDF WEFO SMARTAQUA Projects.

## CONFLICT OF INTEREST

The authors declare no conflict of interest.

## AUTHOR CONTRIBUTIONS

SC and CGL wrote the grant and secured the funding. OR, CGL and SC designed the study. OR, DRB, GC and TUW collected the data. CGL, DRB and OR wrote the MS with input from all authors. DRB and CGL carried out the statistical analysis. All authors approved the final submission.

## DATA AVAILABILITY STATEMENT

All sequences have been submitted to NCBI SRA (Submission ID: SUB4925804; BioProject: PRJNA510732; BioSample SAMN10614327-SAMN10614330). Phenotypic and behavioural data are available from Figshare (https://doi.org/10.6084/m9.fishare.8198408.v1) Rodriguez-Barreto et al., 2019

## ANIMAL ETHICS

This work was screened and approved by Swansea University's College of Science Ethics Committee (permit BS 09062016).

## ORCID

Sonia Consuegra ID https://orcid.org/0000-0003-4403-2509

Carlos Garcia de Leaniz D https://orcid.org/0000-0003-1650-2729

## REFERENCES

Andrews, S. (2010). *FastQC: A quality control tool for high throughput sequence data*. Retrieved from http://www.bioinformatics.babraham.ac.uk/projects/fastqc1662

Ashley, P. J. (2007). Fish welfare: Current issues in aquaculture. *Applied Animal Behaviour Science*, 104(3), 199-235. https://doi.org/10.1016/j.applanim.2006.09.001

Aubin-Horth, N., Landry, C. R., Letcher, B. H., & Hofmann, H. A. (2005). Alternative life histories shape brain gene expression profiles in males of the same population. *Proceedings of the Royal Society B: Biological Sciences*, 272(1573), 1655-1662. https://doi.org/10.1098/rspb.2005.3125

Auperin, B., Baroiller, J.-F., Ricordel, M.-J., Fostier, A., & Prunet, P. (1997). Effect of confinement stress on circulating levels of growth hormone and two prolactins in freshwater-adapted tilapia (*Oreochromis niloticus*). *General and Comparative Endocrinology*, 108(1), 35-44. https://doi.org/10.1006/gcen.1997.6938

Avella, M., Schreck, C. B., & Prunet, P. (1991). Plasma prolactin and cortisol concentrations of stressed coho salmon, *Oncorhynchus kisutch*, in fresh water or salt water. *General and Comparative Endocrinology*, 81(1), 21-27. https://doi.org/10.1016/0016-6480(91)90121-L

Balment, R., Lu, W., Weybourne, E., & Warne, J. (2006). Arginine vasotocin a key hormone in fish physiology and behaviour: A review with insights from mammalian models. *General and Comparative Endocrinology*, 147(1), 9-16. https://doi.org/10.1016/j.ygcen.2005.12.022

Balzarini, V., Taborsky, M., Wanner, S., Koch, F., & Frommen, J. G. (2014). Mirror, mirror on the wall: The predictive value of mirror tests for
measuring aggression in fish. *Behavioral Ecology and Sociobiology*, 68(5), 871-878. https://doi.org/10.1007/s00265-014-1698-7

Barcellos, L., Nicolaiewsky, S., De Souza, S., & Lulhier, F. (1999). The effects of stocking density and social interaction on acute stress response in Nile tilapia *Oreochromis niloticus* (L.) fingerlings. *Aquaculture Research*, 30(11-12), 887-892. https://doi.org/10.1046/j.1365-2109.1999.00419.x

Barreto, R. E., Volpato, G. L., de Brito Faturi, C., Giaquinto, P. C., Gonçalves de Freitas, E., & Fernandes de Castilho, M. (2009). Aggressive behaviour traits predict physiological stress responses in Nile tilapia (*Oreochromis niloticus*). *Marine and Freshwater Behaviour and Physiology*, 42(2), 109-118.

Barrows, E. M. (2000). *Animal behavior desk reference: A dictionary of animal behavior, ecology, and evolution*. Boca Raton, FL: CRC Press.

Barton, K. (2015). *Package 'MuMln'. Multi-Model inference, version 1*. Retrieved from https://cran.r-project.org/web/packages/MuMln/index.html

Bates, D., Maechler, M., Bolker, B., & Walker, S. (2014). *Ime4: Linear mixed-effects models using Eigen and S4*. R Package v. 1.7.

Bell, A. M., Hankison, S. J., & Laskowski, K. L. (2009). The repeatability of behaviour: A meta-analysis. *Animal Behaviour*, 77, 771-783. https://doi.org/10.1016/j.anbehav.2008.12.022

Belyaev, D. K. (1969). Domestication of animals. *Science*, 5(1), 47-52.

Bernier, N. J. (2006). The corticotropin-releasing factor system as a mediator of the appetite-suppressing effects of stress in fish. *General and Comparative Endocrinology*, 146(1), 45-55. https://doi.org/10.1016/j.ygcen.2005.11.016

Blumstein, D. T., Ebensperger, L., Hayes, L., Vásquez, R. A., Ahern, T. H., Burger, J. R., ... Young, L. J. (2010). Towards an integrative understanding of social behavior: new models and new opportunities. *Frontiers in Behavioral Neuroscience*, 4, 34. https://doi.org/10.3389/fbehav.2010.00034

Bolger, A. M., Lohse, M., & Usadel, B. (2014). Trimmomatic: A flexible trimmer for Illumina sequence data. *Bioinformatics*, 30(15), 2114-2120. https://doi.org/10.1093/bioinformatics/btu170

Bornstein, S., & Chrousos, G. (1999). Adrenocorticotropin (ACTH)-and non-ACTH-mediated regulation of the adrenal cortex: Neural and immune inputs. *The Journal of Clinical Endocrinology & Metabolism*, 84(5), 1729-1736. https://doi.org/10.1210/jcem.84.5.5631

Brawand, D., Wagner, C. E., Li, Y. I., Malinsky, M., Keller, I., Fan, S., ... Di Palma, F. (2014). The genomic substrate for adaptive radiation in African cichlid fish. *Nature*, 513, 375-381. https://doi.org/10.1038/nature13726

Brown, M. J. A. & Shahidi, F. (1997). Effects of stocking density on colour characteristics and deposition of carotenoids in cultured Arctic charr (*Salvelinus alpinus*). *Food Chemistry*, 59(1), 107-114. https://doi.org/10.1016/S0308-8146(96)00205-1

Champneys, T., Castaldo, G., Consuegra, S., & Garcia de Leaniz, C. (2018). Density-dependent changes in neophobia and stress-coping styles in the world's oldest farmed fish. *Royal Society Open Science*, 5(12), 181473. https://doi.org/10.1101/394288

Clarke, J. M., & Schlüter, D. (2011). Colour plasticity and background matching in a threespine stickleback species pair. *Biological Journal of the Linnean Society*, 102, 902-914. https://doi.org/10.1111/j.1095-8312.2011.01623.x

Conesa, A., Madrigal, P., Tarazona, S., Gomez-Cabrero, D., Cervera, A., McPherson, A., ... Mortazavi, A. (2016). A survey of best practices for RNA-seq data analysis. *Genome Biology*, 17(1), 13. https://doi.org/10.1186/s13059-016-0881-8

Conte, F. S. (2004). Stress and the welfare of cultured fish. *Applied Animal Behaviour Science*, 86(3), 205-223. https://doi.org/10.1016/j.appla nim.2004.02.003

Culbert, B. M., Gilmour, K. M., & Balshine, S. (2018). Stress axis regulation during social ascension in a group-living cichlid fish.

*Hormones and Behaviour*, 103, 121-128. https://doi.org/10.1016/j.ybehav.2018.06.007

Davis, E. S., & Marler, C. A. (2004). C-fos changes following an aggressive encounter in female California mice: A synthesis of behavior, hormone changes and neural activity. *Neuroscience*, 127(3), 611-624. https://doi.org/10.1016/j.neuroscience.2004.05.034

Eissa, N., & Wang, H. P. (2016). Transcriptional stress responses to environmental and husbandry stressors in aquaculture species. *Reviews in Aquaculture*, 8(1), 61-88. https://doi.org/10.1111/raq.12081

Ellis, S. E., Gupta, S., Ashar, F. N., Bader, J. S., West, A. B., & Arking, D. E. (2013). RNA-Seq optimization with eQTL gold standards. BMC Genomics, 14(1), 892. https://doi.org/10.1186/1471-2164-14-892

Ellis, T., North, B., Scott, A., Bromage, N., Porter, M., & Gadd, D. (2002). The relationships between stocking density and welfare in farmed rainbow trout. *Journal of Fish Biology*, 61(3), 493-531. https://doi.org/10.1006/jfbi.2002.2057

Ellison, A. R., Uren Webster, T. M., Rey, O., Garcia de Leaniz, C., Consuegra, S., Orozco-terWengel, P., & Cable, J. (2018). Transcriptomic response to parasite infection in Nile tilapia (*Oreochromis niloticus*) depends on rearing density. BMC Genomics, 19, 723. https://doi.org/10.1186/s12864-018-5098-7

El-Sayed, A. F. M. (2002). Effects of stocking density and feeding levels on growth and feed efficiency of Nile tilapia (*Oreochromis niloticus* L.) fry. *Aquaculture Research*, 33(8), 621-626. https://doi.org/10.1046/j.1365-2109.2002.00700.x

El-Sayed, A.-F.-M. (2006). *Tilapia culture*. Oxford, UK: CABI Publishing.

FAO (2016). *The State of World Fisheries and Aquaculture 2016. Contributing to Food Security and Nutrition for All*. Retrieved from Rome.

FAWC (2014). *Opinion on the welfare of farmed fish* (p. 40). London, UK: Farm Animal Welfare Committee. Retrieved from http://www.fawc.org.uk/freedoms.htm

Fernald, R. D. (2012). Social control of the brain. *Annual Review of Neuroscience*, 35, 133-151. https://doi.org/10.1146/annur ev-neuro-062111-150520

Fessehaye, Y., Kabir, A., Bovenhuis, H., & Komen, H. (2006). Prediction of cannibalism in juvenile *Oreochromis niloticus* based on predator to prey weight ratio, and effects of age and stocking density. *Aquaculture*, 255(1), 314-322. https://doi.org/10.1016/j.aquac ulture.2005.11.033

Filby, A. L., Paull, G. C., Hickmore, T. F., & Tyler, C. R. (2010). Unravelling the neurophysiological basis of aggression in a fish model. BMC Genomics, 11(1), 498. https://doi.org/10.1186/1471-2164-11-498

Freudenberg, F., Carreño Gutierrez, H., Post, A. M., Reif, A., & Norton, W. H. (2016). Aggression in non-human vertebrates: Genetic mechanisms and molecular pathways. *American Journal of Medical Genetics Part B: Neuropsychiatric Genetics*, 171, 603-640.

Gall, G. A., & Bakar, Y. (1999). Stocking density and tank size in the design of breed improvement programs for body size of tilapia. *Aquaculture*, 173(1), 197-205. https://doi.org/10.1016/S0044-8486(98)00487-6

Gallardo-Pujol, D., Andrés-Pueyo, A., & Maydeu-Olivares, A. (2013). MAOA genotype, social exclusion and aggression: An experimental test of a gene-environment interaction. *Genes, Brain and Behavior*, 12(1), 140-145.

Gerlai, R., Lahav, M., Guo, S., & Rosenthal, A. (2000). Drinks like a fish: Zebra fish (*Danio rerio*) as a behavior genetic model to study alcohol effects. *Pharmacology Biochemistry and Behavior*, 67, 773-782. https://doi.org/10.1016/S0091-3057(00)00422-6

Godfray, H. C. J., Beddington, J. R., Crute, I. R., Haddad, L., Lawrence, D., Muir, J. F., ... Toulmin, C. (2010). Food security: The challenge of feeding 9 billion people. *Science*, 327(5967), 812-818. https://doi.org/10.1126/science.1185383

Goldburg, R., & Naylor, R. (2005). Future seascapes, fishing, and fish farming. *Frontiers in Ecology and the Environment*, 3(1), 21-28. https://doi.org/10.1890/1540-9295(2005)003[0021:FSFAFF]2.0.CO;2
Goodson, J. L. (2005). The vertebrate social behavior network: Evolutionary themes and variations. *Hormones and Behavior, 48*(1), 11-22. https://doi.org/10.1016/j.yhbeh.2005.02.003

Haller, J., Tóth, M., Halasz, J., & De Boer, S. F. (2006). Patterns of violent aggression-induced brain c-fos expression in male mice selected for aggressiveness. *Physiology & Behavior, 88*(1-2), 173-182. https://doi.org/10.1016/j.physbeh.2006.03.030

Harrison, X. A. (2014). Using observation-level random effects to model overdispersion in count data in ecology and evolution. *PeerJ*, 2, e616. https://doi.org/10.7717/peerj.616

Hofmann, H. A., & Fernald, R. D. (2000). Social status controls somatostatin neuron size and growth. *Journal of Neuroscience*, 20(12), 4740-4744. https://doi.org/10.1523/JNEUROSCI.20-12-04740.2000

Höglund, E., Balm, P. H. M., & Winberg, S. (2000). Skin darkening, a potential social signal in subordinate Arctic charr (*Salvelinus alpinus*): The regulatory role of brain monoamines and pro-opiomelanocortin-derived peptides. *Journal of Experimental Biology*, 203, 1711-1721.

Huang, D. A., Sherman, B. T., Tan, Q., Collins, J. R., Alvord, W. G., Roayaei, J., ... Lempicki, R. A. (2007). The DAVID gene functional classification tool: A novel biological module-centric algorithm to functionally analyze large gene lists. *Genome Biology*, 8(9), R183-R183. https://doi.org/10.1186/gb-2007-8-9-r183

Huntingford, F. A. (2004). Implications of domestication and rearing conditions for the behaviour of cultivated fishes. *Journal of Fish Biology*, 65, 122-142. https://doi.org/10.1111/j.0022-1112.2004.00562.x

Huntingford, F., & Adams, C. (2005). Behavioural syndromes in farmed fish: Implications for production and welfare. *Behaviour*, 142(9-10), 1207-1221.

Huntingford, F. A., Adams, C., Braithwaite, V. A., Kadri, S., Pottinger, T. G., Sandøe, P., & Turnbull, J. F. (2006). Current issues in fish welfare. *Journal of Fish Biology*, 68(2), 332-372. https://doi.org/10.1111/j.0022-1112.2006.00104.x

Jensen, P. (2014). Behavior genetics and the domestication of animals. *Annual Review of Animal Biosciences*, 2(1), 85-104. https://doi.org/10.1146/annurev-animal-022513-114135

Jezova, D., Skultetyova, I., Tokarev, D. I., Bakos, P., & Vigas, M. (1995). Vasopressin and oxytocin in stress. *Annals of the New York Academy of Sciences*, 771(1), 192-203. https://doi.org/10.1111/j.1749-6632.1995.tb44681.x

Kihslinger, R., Lema, S. C., & Nevitt, G. (2006). Environmental rearing conditions produce forebrain differences in wild Chinook salmon *Oncorhynchus tshawytscha*. *Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology*, 145(2), 145-151. https://doi.org/10.1016/j.cbpa.2006.06.041

Kobayashi, Y., Mizusawa, K., Chiba, H., Tagawa, M., & Takahashi, A. (2012). Further evidence on acetylation-induced inhibition of the pigment-dispersing activity of α-melanocyte-stimulating hormone. *General and Comparative Endocrinology*, 176(1), 9-17. https://doi.org/10.1016/j.ygcen.2011.12.001

Kobayashi, Y., Mizusawa, K., Saito, Y., & Takahashi, A. (2012). Melanocortin systems on pigment dispersion in fish chromatophores. *Frontiers in Endocrinology*, 3, 9. https://doi.org/10.3389/fendo.2012.00009.

Komers, P. E. (1997). Behavioural plasticity in variable environments. *Canadian Journal of Zoology*, 75(2), 161-169. https://doi.org/10.1139/z97-023

Koolhaas, J., De Boer, S., Coppens, C., & Buwalda, B. (2010). Neuroendocrinology of coping styles: Towards understanding the biology of individual variation. *Frontiers in Neuroendocrinology*, 31(3), 307-321. https://doi.org/10.1016/j.yfrne.2010.04.001

Koolhaas, J. M., Korte, S. M., De Boer, S. F., Van Der Vegt, B. J., Van Reenen, C. G., Hopster, H., ... Blokhuis, H. J. (1999). Coping styles in animals: Current status in behavior and stress-physiology. *Neuroscience & Biobehavioral Reviews*, 23, 925-935. https://doi.org/10.1016/S0149-7634(99)00026-3

Kotrschal, A., Deacon, A. E., Magurran, A. E., & Kolm, N. (2017). Predation pressure shapes brain anatomy in the wild. *Evolutionary Ecology*, 31(5), 619-633. https://doi.org/10.1007/s10682-017-9901-8

Kotrschal, A., Rogell, B., Bundsen, A., Svensson, B., Zajitschek, S., Brännström, I., ... Kolm, N. (2013). Artificial selection on relative brain size in the guppy reveals costs and benefits of evolving a larger brain. *Current Biology*, 23(2), 168-171. https://doi.org/10.1016/j.cub.2012.11.058

Kovács, K. J. (1998). c-Fos as a transcription factor: A stressful (re) view from a functional map. *Neurochemistry International*, 33(4), 287-297.

Kuznetsova, A., Brockhoff, P. B., & Christensen, R. H. B. (2017). ImerTest package: Tests in linear mixed effects models. *Journal of Statistical Software*, 82(13), 1-26. https://doi.org/10.18637/jss.v082.i13

Landau, W. M., & Liu, P. (2013). Dispersion estimation and its effect on test performance in RNA-seq data analysis: A simulation-based comparison of methods. *PLoS ONE*, 8(12), e81415. https://doi.org/10.1371/journal.pone.0081415

Lawrence, M., Huber, W., Pagès, H., Aboyoun, P., Carlson, M., Gentleman, R., ... Carey, V. (2013). Software for computing and annotating genomic ranges. *PLoS Computational Biology*, 9(8), e1003118. https://doi.org/10.1371/journal.pcbi.1003118

Leclercq, E., Taylor, J. F., & Migaud, H. (2010). Morphological skin colour changes in teleosts. *Fish and Fisheries*, 11(2), 159-193. https://doi.org/10.1111/j.1467-2979.2009.00346.x

Love, M. I., Huber, W., & Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. *Genome Biology*, 15, 550(12).

Lupien, S. J., McEwen, B. S., Gunnar, M. R., & Heim, C. (2009). Effects of stress throughout the lifespan on the brain, behaviour and cognition. *Nature Reviews Neuroscience*, 10(6), 434. https://doi.org/10.1038/nrn2639

Machluf, Y., Gutnick, A., & Levkowitz, G. (2011). Development of the zebrafish hypothalamus. *Annals of the New York Academy of Sciences*, 1220(1), 93-105. https://doi.org/10.1111/j.1749-6632.2010.05945.x

Malki, K., Du Rietz, E., Crusio, W. E., Pain, O., Paya-Cano, J., Karadaghi, R. L., ... Tosto, M. G. (2016). Transcriptome analysis of genes and gene networks involved in aggressive behavior in mouse and zebrafish. *American Journal of Medical Genetics Part B: Neuropsychiatric Genetics*, 171(6), 827-838. https://doi.org/10.1002/ajmg.b.32451

Marks, C., West, T. N., Bagatto, B., Moore, F. B., & Taylor, C. (2005). Developmental environment alters conditional aggression in zebrafish. *Copeia*, 2005(4), 901-908. https://doi.org/10.1643/0045-8511(2005)005[0901:DEACAI]2.0.CO;2

Martins, C. I. M., Galhardo, L., Noble, C., Damsgård, B., Spedicato, M. T., Zupa, W., ... Kristiansen, T. (2012). Behavioural indicators of welfare in farmed fish. *Fish Physiology and Biochemistry*, 38(1), 17-41. https://doi.org/10.1007/s10695-011-9518-8

Maruska, K. P. (2015). Social transitions cause rapid behavioral and neuroendocrine changes. *Integrative and Comparative Biology*, 55(2), 294-306. https://doi.org/10.1093/icb/icv057

Maruska, K. P., & Fernald, R. D. (2014). Social regulation of gene expression in the African cichlid fish *Astatotilapia burtoni*. In T. Canli (Ed.), *The Oxford handbook of molecular psychology* (pp. 52-78). Oxford, UK: Oxford University Press.

Maxson, C. S., & Canastar, A. (2005). Genetic aspects of aggression in nonhuman animals. In R. J. Nelson (Ed.), *Biology of aggression* (pp. 3-19). New York, NY: Oxford University Press.

Mormède, P. (2005). Molecular genetics of behaviour: Research strategies and perspectives for animal production. *Livestock Production Science*, 93(1), 15-21. https://doi.org/10.1016/j.livprodsci.2004.11.002

Nelson, R. J., & Trainor, B. C. (2007). Neural mechanisms of aggression. *Nature Reviews Neuroscience*, 8(7), 536. https://doi.org/10.1038/nrn2174

O'Connell, L. A., & Hofmann, H. A. (2012). Evolution of a vertebrate social decision-making network. *Science*, 336(6085), 1154-1157.
Oliveira, R. F. (2012). Social plasticity in fish: Integrating mechanisms and function. *Journal of Fish Biology*, **81**(7), 2127-2150. https://doi.org/10.1111/j.1095-8649.2012.03477.x

Oshima, N., Makino, M., Iwamura, S., & Bern, H. A. (1996). Pigment dispersion by prolactin in cultured xanthophores and erythrophores of some fish species. *Journal of Experimental Zoology*, **275**(1), 45-52. https://doi.org/10.1002/(SICI)1097-010X(19960501)275:1<45;AID-JEZ7>3.0.CO;2-A

Pavlidis, M., Sundvik, M., Chen, Y.-C., & Panula, P. (2011). Adaptive changes in zebrafish brain in dominant-subordinate behavioral context. *Behavioural Brain Research*, **225**(2), 529-537. https://doi.org/10.1016/j.bbr.2011.08.022

Pavlidis, M., Theodoridi, A., & Tsalafouta, A. (2015). Neuroendocrine regulation of the stress response in adult zebrafish, *Danio rerio*. *Progress in Neuro-Psychopharmacology and Biological Psychiatry*, **60**, 121-131. https://doi.org/10.1016/j.pnpbp.2015.02.014

Price, E. O. (1999). Behavioral development in animals undergoing domestication. *Applied Animal Behaviour Science*, **65**(3), 245-271. https://doi.org/10.1016/S0168-1591(99)00087-8

R Core Team (2017). *R: A language and environment for statistical computing*. Vienna, Austria: R Foundation for Statistical Computing. Retrieved from https://www.R-project.org/

Roberts, L. J., Taylor, J., & Garcia de Leaniz, C. (2011). Environmental enrichment reduces maladaptive risk-taking behavior in salmon reared for conservation. *Biological Conservation*, **144**(7), 1972-1979. https://doi.org/10.1016/j.biocon.2011.04.017

Robinson, G. E., Fernald, R. D., & Clayton, D. F. (2008). Genes and social behavior. *Science*, **322**, 896-900. https://doi.org/10.1126/science.1159277

Rodriguez-Barreto, D., Rey, O., Uren Webster, T. M., Castaldo, G., Consuegra, S., & Garcia de Leaniz, C. (2019). Transcriptomic response to aquaculture intensification in Nile tilapia. *Figshare*. https://doi.org/10.6084/m9 figshare.8198408.v1

Rohde, P. D., Gaertner, B., Wards, K., Sørensen, P., & Mackay, T. F. C. (2017). Genomic analysis of genotype by social environment interaction for *Drosophila aggressive* behavior. *Genetics*, **206**(4), 1969-1984.

Sokolowski, M. B. (2010). Social interactions in “simple” model systems. *Neuron*, **65**(6), 780-794. https://doi.org/10.1016/j.neuron.2010.03.007

Solomon, R. W. (2009). Free and open source software for the manipulation of digital images. *American Journal of Roentgenology*, **192**(6), W330-W334. https://doi.org/10.2214/AJR.OB.2190

Trainor, B. C., & Hofmann, H. A. (2006). Somatostatin regulates aggressive behavior in an African cichlid fish. *Endocrinology*, **147**(11), 5119-5125. https://doi.org/10.1210/en.2006-0511

Trainor, B. C., & Hofmann, H. A. (2007). Somatostatin and somatostatin receptor gene expression in dominant and subordinate males of an African cichlid fish. *Behavioural Brain Research*, **179**(2), 314-320. https://doi.org/10.1016/j.bbr.2007.02.014

Trainor, B. C., Lin, S., Finy, M. S., Rowland, M. R., & Nelson, R. J. (2007). Photoperiod reverses the effects of estrogens on male aggression via genomic and nongenomic pathways. *Proceedings of the National Academy of Sciences of the United States of America*, **104**(23), 9840-9845. https://doi.org/10.1073/pnas.0701819104

Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., ... Pachter, L. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with *TopHat* and *Cufflinks*. *Nature Protocols*, **7**(3), 562. https://doi.org/10.1038/nprot.2012.016

Uren Webster, T. M., Rodriguez-Barreto, D., Martin, S. A. M., van Oosterhout, C., Orozco-terWengel, P., Cable, J., ... Consuegra, S. (2018). Contrasting effects of acute and chronic stress on the transcriptome, epigenome, and immune response of Atlantic salmon. *Epigenetics*, **13**(12), 1191-1207. https://doi.org/10.1080/1559294.2018.1554520

Van der Salm, A. L., Martinez, M., Flik, G., & Wendelaar Bonga, S. E. (2004). Effects of husbandry conditions on the skin colour and stress response of red porgy, *Pagrus pagrus*. *Aquaculture*, **241**(1-4), 371-386. https://doi.org/10.1016/j.aquaculture.2004.08.038

Waltes, R., Chiocchetti, A. G., & Freitag, C. M. (2016). The neurobiological basis of human aggression: A review on genetic and epigenetic mechanisms. *American Journal of Medical Genetics Part B: Neuropsychiatric Genetics*, **171**, 650-675.

Way, G. P., Ruhl, N., Snekser, J. L., Kiesel, A. L., & McRobert, S. P. (2015). A comparison of methodologies to test aggression in zebrafish. *Zebrafish*, **12**, 144-151. https://doi.org/10.1089/zeb.2014.1025

Weitekamp, C. A., Nguyen, J., & Hofmann, H. A. (2017). Social context affects behavior, preoptic area gene expression, and response to D2 receptor manipulation during territorial defense in a cichlid fish. *Genes, Brain and Behavior*, **16**(6), 601-611. https://doi.org/10.1111/gbb.12389

Wilkins, A. S., Wrangham, R. W., & Fitch, W. T. (2014). The “Domestication Syndrome” in mammals: A unified explanation based on neural crest cell behavior and genetics. *Genetics*, **197**(3), 795-808. https://doi.org/10.1534/genetics.114.165423

Zabegalov, K. N., Kolesnikova, T. O., Khatsko, S. L., Volgin, A. D., Yakovlev, O. A., Amstislavskaya, T. G., ... Kalueff, A. V. (2019). Understanding zebrafish aggressive behavior. *Behavioural Processes*, **158**, 200-210. https://doi.org/10.1016/j.beproc.2018.11.010

Zeng, W., Li, Z., Ye, S., Xie, S., Liu, J., Zhang, T., & Duan, M. (2010). Effects of stocking density on growth and skin color of juvenile darkbarbel catfish *Pelteobagrus vachelli* (Richardson). *Journal of Applied Ichthyology*, **26**(6), 925-929. https://doi.org/10.1111/j.1439-0426.2010.01523.x

## SUPPORTING INFORMATION

Additional supporting information may be found online in the Supporting Information section at the end of the article.

**How to cite this article:** Rodriguez-Barreto D, Rey O, Uren-Webster TM, Castaldo G, Consuegra S, Garcia de Leaniz C. Transcriptomic response to aquaculture intensification in Nile tilapia. *Evol Appl.* 2019;12:1757-1771. https://doi.org/10.1111/eva.12830