Unraveling the genotype by environment interaction in a thermosensitive fish with a polygenic sex determination system.

In most animals, sex determination occurs at conception, when sex chromosomes are segregated following Mendelian laws. However, in multiple reptiles and fishes, this genetic sex can be overridden by external factors after fertilization or birth. In some species, the genetic sex may also be governed by multiple genes, further limiting our understanding of sex determination in such species. We used the European sea bass (Dicentrarchus labrax) as a model and combined genomic (using a single nucleotide polymorphism chip) and transcriptomic (RNA-Sequencing) approaches to thoroughly depict this polygenic sex determination system and its interaction with temperature. We estimated genetic sex tendency (eGST), defined as the estimated genetic liability to become a given sex under a liability threshold model for sex determination, which accurately predicts the future phenotypic sex. We found evidence that energetic pathways, concerning the regulation of lipids and glucose, are involved in sex determination and could explain why females tend to exhibit higher energy levels and improved growth compared to males. Besides, early exposure to high-temperature up-regulated sox3, followed by sox9a in individuals with intermediate eGST, but not in individuals showing highly female-biased eGST, providing the most parsimonious explanation for temperature-induced masculinization. This gonadal state was maintained likely by DNA methylation and the up-regulation of several genes involved in histone modifications, including jmjd1c Overall, we describe a sex determination system resulting from continuous genetic and environmental influences in an animal. Our results provide significant progress in our understanding of the mechanisms underlying temperature-induced masculinization in fish.

In a marine teleost, the significance of oxygen supply for acute thermal tolerance depends upon the context and the endpoint used.

Eight juvenile European seabass were exposed to two thermal ramping protocols with different levels of aerobic activity and tolerance endpoint: the critical thermal maximum for swimming (CTSmax) while exercising aerobically until fatigue and the critical thermal maximum (CTmax) under static conditions until loss of equilibrium (LOE). In the CTSmax protocol, warming caused a profound increase in the rate of oxygen uptake (MO2), culminating in a gait transition from steady aerobic towards unsteady anaerobic swimming, then fatigue at 30.3±0.4°C (mean±s.e.m.). Gait transition and fatigue presumably indicate an oxygen limitation, which reflects the inability to meet the combined demands of swimming plus warming. The CTmax protocol also elicited an increase in MO2, culminating in LOE at 34.0±0.4°C, which is significantly warmer than fatigue at CTSmax. The maximum MO2 achieved in the CTmax protocol was, however, less than 30% of that achieved in the CTSmax protocol. Therefore, the static CTmax did not exploit full cardiorespiratory capacity for oxygen supply, indicating that LOE was not caused by systemic oxygen limitation. Consequently, systemic oxygen supply can be significant for tolerance of acute warming in seabass but this depends upon the physiological context and the endpoint used.

Whole-genome sequencing identifies interferon-induced protein IFI6/IFI27-like as a strong candidate gene for VNN resistance in European sea bass.

BACKGROUND: Viral nervous necrosis (VNN) is a major disease that affects European sea bass, and understanding the biological mechanisms that underlie VNN resistance is important for the welfare of farmed fish and sustainability of production systems. The aim of this study was to identify genomic regions and genes that are associated with VNN resistance in sea bass. RESULTS: We generated a dataset of 838,451 single nucleotide polymorphisms (SNPs) identified from whole-genome sequencing (WGS) in the parental generation of two commercial populations (A: 2371 individuals and B: 3428 individuals) of European sea bass with phenotypic records for binary survival in a VNN challenge. For each population, three cohorts were submitted to a red-spotted grouper nervous necrosis virus (RGNNV) challenge by immersion and genotyped on a 57K SNP chip. After imputation of WGS SNPs from their parents, quantitative trait loci (QTL) were mapped using a Bayesian sparse linear mixed model (BSLMM). We found several QTL regions that were specific to one of the populations on different linkage groups (LG), and one 127-kb QTL region on LG12 that was shared by both populations and included the genes ZDHHC14, which encodes a palmitoyltransferase, and IFI6/IFI27-like, which encodes an interferon-alpha induced protein. The most significant SNP in this QTL region was only 1.9 kb downstream of the coding sequence of the IFI6/IFI27-like gene. An unrelated population of four large families was used to validate the effect of the QTL. Survival rates of susceptible genotypes were 40.6% and 45.4% in populations A and B, respectively, while that of the resistant genotype was 66.2% in population B and 78% in population A. CONCLUSIONS: We have identified a genomic region that carries a major QTL for resistance to VNN and includes the ZDHHC14 and IFI6/IFI27-like genes. The potential involvement of the interferon pathway, a well-known anti-viral defense mechanism in several organisms (chicken, human, or fish), in survival to VNN infection is of particular interest. Our results can lead to major improvements for sea bass breeding programs through marker-assisted genomic selection to obtain more resistant fish.

Natural cortisol production is not linked to the sexual fate of European sea bass.

In this study, we aimed to investigate the relationship between cortisol and the determination of sexual fate in the commercially important European sea bass (Dicentrarchus labrax). To test our hypothesis, we designed two temperature-based experiments (19 ℃, 21 ℃ and 23 ℃, experiment 1; 16 ℃ and 21 ℃, experiment 2) to assess the effects of these thermal treatments on European sea bass sex determination and differentiation. In the fish from the first experiment, we evaluated whether blood cortisol levels and expression of stress key regulatory genes were different between differentiating (149 to 183 dph) males and females. In the second experiment, we assessed whether cortisol accumulated in scales over time during the labile period for sex determination as well as the neuroanatomical localisation of brain cells expressing brain aromatase (cyp19a1b) and corticotropin-releasing factor (crf) differed between males and females undergoing molecular sex differentiation (117 to 124 dph). None of the gathered results allowed to detect differences between males and females regarding cortisol production and regulatory mechanisms. Altogether, our data provide strong physiological, molecular and histochemical evidence, indicating that in vivo cortisol regulation has no major effects on the sex of European sea bass.

The genetic correlation between feed conversion ratio and growth rate affects the design of a breeding program for more sustainable fish production.

BACKGROUND: Most fish breeding programs aim at improving growth rate and include feed conversion ratio (FCR) neither in the breeding goal nor in the selection index, although decreasing FCR is known to increase farm profit and decrease environmental impacts. This is because FCR is difficult to measure in fish that live in groups and FCR is assumed to have a favourable (negative) genetic correlation with growth, although the magnitude of this correlation is unknown. We investigated the effect of the genetic correlation between growth and FCR on the economic and environmental responses of a two-trait breeding goal (growth and FCR), compared to a single-trait breeding goal (growth only). Next, we evaluated the weights to assign to growth and FCR in a two-trait breeding goal to maximize sustainability of fish production. METHODS: We used pseudo-best linear unbiased prediction (BLUP) index calculations to simulate a breeding program for sea bass. For the single-trait breeding goal, the trait in the breeding goal and in the index was thermal growth coefficient (TGC) and for the two-trait breeding goal, the traits in the breeding goal were TGC and FCR and the traits in the index were TGC and percentage of fat in the dorsal muscle (an indirect measure of FCR). We simulated responses to selection for genetic and phenotypic correlations between TGC and FCR ranging from 0 to - 0.8. Then, in the two-trait breeding goal, we calculated the economic return and the change in eutrophication when using economic values (EV) or environmental values (ENV). RESULTS: When the genetic correlation between TGC and FCR was lower than - 0.45, we found major differences in economic returns and in eutrophication between single and two-trait breeding programs. At a correlation of - 0.25, the two-trait breeding goal based on EV increased economic return by 25% compared to the single-trait breeding goal, while using ENV decreased eutrophication by 1.34% per ton of fish produced after one generation of selection. CONCLUSIONS: The genetic correlation between TGC and FCR affects the magnitude of economic losses due to omitting FCR in the breeding program. In addition, the genetic correlation affects the importance of choosing EV or ENV to reduce eutrophication and increase profit.

Population-specific variations of the genetic architecture of sex determination in wild European sea bass Dicentrarchus labrax L.

Polygenic sex determination (PSD) may show variations in terms of genetic and environmental components between populations of fish species exposed/adapted to different environments. The European sea bass (Dicentrarchus labrax) is an interesting model, combining both a PSD system and a genetic subdivision into an Atlantic and a Mediterranean lineage, with genetic substructures within the Mediterranean Sea. Here, we produced experimental progeny crosses (N = 927) from broodstock sampled in four wild populations (North Atlantic, NAT; Western Mediterranean, WEM; North-Eastern Mediterranean, NEM; South-Eastern Mediterranean, SEM). We found less females than males in the progeny, both in the global dataset (32.5%) and within each paternal group (from 25.1% for NEM to 39.0% for WEM), with significant variation among populations, dams, and sires. Sex, body weight (BW), and body length (BL) showed moderate heritability (0.52 ± 0.17, 0.46 ± 0.17, 0.34 ± 0.15, respectively) and sex was genetically correlated with BW and BL, with rAsex/BW = 0.69 ± 0.12 and rA sex/BL = 0.66 ± 0.13. A weighted GWAS performed both on the global dataset and within each paternal group revealed a different genetic architecture of sex determination between Atlantic and Mediterranean populations (9 QTLs found in NAT, 7 in WEM, 5 in NEM, and 4 in SEM, with a cumulated variance explained of 27.04%, 21.87%, 15.89%, and 12.10%, respectively) and a more similar genetic architecture among geographically close populations compared to geographically distant populations, consistent with the hypothesis of a population-specific evolution of polygenic sex determination systems in different environments.

Quantifying the genetic parameters of feed efficiency in juvenile Nile tilapia Oreochromis niloticus.

BACKGROUND: Improving feed efficiency in fish is crucial at the economic, social and environmental levels with respect to developing a more sustainable aquaculture. The important contribution of genetic improvement to achieve this goal has been hampered by the lack of accurate basic information on the genetic parameters of feed efficiency in fish. We used video assessment of feed intake on individual fish reared in groups to estimate the genetic parameters of six growth traits, feed intake, feed conversion ratio (FCR) and residual feed intake in 40 pedigreed families of the GIFT strain of Nile tilapia, Oreochromis niloticus. Feed intake and growth were measured on juvenile fish (22.4 g mean body weight) during 13 consecutive meals, representing 7 days of measurements. We used these data to estimate the FCR response to different selection criteria to assess the potential of genetics as a means of increasing FCR in tilapia. RESULTS: Our results demonstrate genetic control for FCR in tilapia, with a heritability estimate of 0.32 ± 0.11. Response to selection estimates showed FCR could be efficiently improved by selective breeding. Due to low genetic correlations, selection for growth traits would not improve FCR. However, weight loss at fasting has a high genetic correlation with FCR (0.80 ± 0.25) and a moderate heritability (0.23), and could be an easy to measure and efficient criterion to improve FCR by selective breeding in tilapia. CONCLUSION: At this age, FCR is genetically determined in Nile tilapia. A selective breeding program could be possible and could help enabling the development of a more sustainable aquaculture production.

Genome-wide association and genomic prediction of resistance to viral nervous necrosis in European sea bass (Dicentrarchus labrax) using RAD sequencing.

BACKGROUND: European sea bass (Dicentrarchus labrax) is one of the most important species for European aquaculture. Viral nervous necrosis (VNN), commonly caused by the redspotted grouper nervous necrosis virus (RGNNV), can result in high levels of morbidity and mortality, mainly during the larval and juvenile stages of cultured sea bass. In the absence of efficient therapeutic treatments, selective breeding for host resistance offers a promising strategy to control this disease. Our study aimed at investigating genetic resistance to VNN and genomic-based approaches to improve disease resistance by selective breeding. A population of 1538 sea bass juveniles from a factorial cross between 48 sires and 17 dams was challenged with RGNNV with mortalities and survivors being recorded and sampled for genotyping by the RAD sequencing approach. RESULTS: We used genome-wide genotype data from 9195 single nucleotide polymorphisms (SNPs) for downstream analysis. Estimates of heritability of survival on the underlying scale for the pedigree and genomic relationship matrices were 0.27 (HPD interval 95%: 0.14-0.40) and 0.43 (0.29-0.57), respectively. Classical genome-wide association analysis detected genome-wide significant quantitative trait loci (QTL) for resistance to VNN on chromosomes (unassigned scaffolds in the case of 'chromosome' 25) 3, 20 and 25 (P < 1e06). Weighted genomic best linear unbiased predictor provided additional support for the QTL on chromosome 3 and suggested that it explained 4% of the additive genetic variation. Genomic prediction approaches were tested to investigate the potential of using genome-wide SNP data to estimate breeding values for resistance to VNN and showed that genomic prediction resulted in a 13% increase in successful classification of resistant and susceptible animals compared to pedigree-based methods, with Bayes A and Bayes B giving the highest predictive ability. CONCLUSIONS: Genome-wide significant QTL were identified but each with relatively small effects on the trait. Tests of genomic prediction suggested that incorporating genome-wide SNP data is likely to result in higher accuracy of estimated breeding values for resistance to VNN. RAD sequencing is an effective method for generating such genome-wide SNPs, and our findings highlight the potential of genomic selection to breed farmed European sea bass with improved resistance to VNN.

Gene-centromere mapping in meiotic gynogenetic European seabass.

BACKGROUND: Fully isogenic lines in fish can be developed using "mitotic" gynogenesis (suppression of first zygotic mitosis following inactivation of the sperm genome). However, genome-wide verification of the steps in this process has seldom been applied. We used ddRADseq to generate SNP markers in a meiotic gynogenetic family of European seabass (Dicentrarchus labrax): (i) to verify the lack of paternal contribution in a meiotic gynogenetic family; (ii) to generate a gene-centromere map from this family; (iii) to identify telomeric markers that could distinguish mitotic gynogenetics from meiotic gynogenetics, which sometimes arise spontaneously in mitotic gynogenetic families. RESULTS: From a single meiotic gynogenetic family consisting of 79 progeny, 42 million sequencing reads (Illumina, trimmed to 148 bases) resolved 6866 unique RAD-tags. The 340 male-informative SNP markers that were identified confirmed the lack of paternal contribution. A gene-centromere map was constructed based on 804 female-informative SNPs in 24 linkage groups (2n = 48) with a total length of 1251.02 cM (initial LG assignment was based on the seabass genome assembly, dicLab v1). Chromosome arm structure could be clearly discerned from the pattern of heterozygosity in each linkage group in 18 out of 24 LGs: the other six showed anomalies that appeared to be related to issues in the genome assembly. CONCLUSION: Genome-wide screening enabled substantive verification of the production of the gynogenetic family used in this study. The large number of telomeric and subtelomeric markers with high heterozygosity values in the meiotic gynogenetic family indicate that such markers could be used to clearly distinguish between meiotic and mitotic gynogenetics.

How could fully scaled carps appear in natural waters in Madagascar?

The capacity of organisms to rapidly evolve in response to environmental changes is a key feature of evolution, and studying mutation compensation is a way to evaluate whether alternative routes of evolution are possible or not. Common carps (Cyprinus carpio) carrying a homozygous loss-of-function mutation for the scale cover gene fgfr1a1, causing the 'mirror' reduced scale cover, were introduced in Madagascar a century ago. Here we show that carps in Malagasy natural waters are now predominantly covered with scales, though they still all carry the homozygous mutation. We also reveal that the number of scales in mutated carps is under strong polygenic genetic control, with a heritability of 0.49. As a whole, our results suggest that carps submitted to natural selection could evolve a wild-type-like scale cover in less than 40 generations from standing polygenic genetic variation, confirming similar findings mainly retrieved from model organisms.

A new SNP-based vision of the genetics of sex determination in European sea bass (Dicentrarchus labrax).

BACKGROUND: European sea bass (Dicentrarchus labrax) is one of the most important farmed species in Mediterranean aquaculture. The observed sexual growth and maturity dimorphism in favour of females adds value towards deciphering the sex determination system of this species. Current knowledge indicates the existence of a polygenic sex determining determination system that interacts with temperature. This was explored by restriction-site associated DNA (RAD) marker analysis in a test panel of 175 offspring that originated from a factorial cross between two dams and four sires from a single full-sib family. RESULTS: The first high-density single nucleotide polymorphism (SNP) based linkage map for sea bass was constructed, consisting of 6706 SNPs on 24 linkage groups. Indications for putative sex-determining QTL (quantitative trait loci) that were significant at the genome-wide threshold were detected on linkage groups 6, 11 and 18 to 21, although a genome-wide association study (GWAS) did not identify individual significant SNPs at a genome-wide threshold. A preliminary genomic prediction approach that tested the efficiency of SNP-based selection for female sea bass showed a slight advantage compared to traditional pedigree-based selection. However, when the same models were tested on the same animals for selection for greater length, a clear advantage of the SNP-based selection was observed. CONCLUSIONS: Overall, the results of this study provide additional support to the polygenic sex determination hypothesis in sea bass. In addition, identification of sex-ratio QTL may provide new opportunities for sex-ratio control in sea bass.

Physiological mechanisms underlying individual variation in tolerance of food deprivation in juvenile European sea bass, Dicentrarchus labrax.

Although food deprivation is a major ecological pressure in fishes, there is wide individual variation in tolerance of fasting, whose mechanistic bases are poorly understood. Two thousand individually tagged juvenile European sea bass were submitted to two 'fasting/feeding' cycles each comprising 3 weeks of food deprivation followed by 3 weeks of ad libitum feeding at 25°C. Rates of mass loss during the two fasting periods were averaged for each individual to calculate a population mean. Extreme fasting tolerant (FT) and sensitive (FS) phenotypes were identified that were at least one and a half standard deviations, on opposing sides, from this mean. Respirometry was used to investigate two main hypotheses: (1) tolerance of food deprivation reflects lower mass-corrected routine metabolic rate (RMR) in FT phenotypes when fasting, and (2) tolerance reflects differences in substrate utilisation; FT phenotypes use relatively less proteins as metabolic fuels during fasting, measured as their ammonia quotient (AQ), the simultaneous ratio of ammonia excretion to RMR. There was no difference in mean RMR between FT and FS over 7 days fasting, being 6.70±0.24 mmol h(-1) fish(-1) (mean ± s.e.m., N=18) versus 6.76±0.22 mmol h(-1) fish(-1) (N=17), respectively, when corrected to a body mass of 130 g. For any given RMR, however, the FT lost mass at a significantly lower rate than FS, overall 7-day average being 0.72±0.05 versus 0.90±0.05 g day(-1) fish(-1), respectively (P<0.01, t-test). At 20 h after receiving a ration equivalent to 2% body mass as food pellets, ammonia excretion and simultaneous RMR were elevated and similar in FT and FS, with AQs of 0.105±0.009 and 0.089±0.007, respectively. At the end of the period of fasting, ammonia excretion and RMR had fallen in both phenotypes, but AQ was significantly lower in FT than FS, being 0.038±0.004 versus 0.061±0.005, respectively (P<0.001, t-test). There was a direct linear relationship between individual fasted AQ and rate of mass loss, with FT and FS individuals distributed at opposing lower and upper extremities, respectively. Thus the difference between the phenotypes in their tolerance of food deprivation did not depend upon their routine energy use when fasting. Rather, it depended upon their relative use of tissue proteins as metabolic fuels when fasting, which was significantly lower in FT phenotypes.

Thermal modulation of energy allocation during sex determination in the European sea bass (Dicentrarchus labrax).

Water temperature governs physiological functions such as growth, energy allocation, and sex determination in ectothermic species. The European sea bass (Dicentrarchus labrax) is a major species in European aquaculture, exhibiting early dimorphic growth favoring females. The species has a polygenic sex determination system that interacts with water temperature to determine an individual's sex, with two periods during development that are sensitive to temperature. The current study investigated the influence of water temperature on energy allocation and sex-biased genes during sex determination and differentiation periods. RNA-Sequencing and qPCR analyses were conducted in two separate experiments, of either constant water temperatures typical of aquaculture conditions or natural seasonal thermal regimes, respectively. We focused on eight key genes associated with energy allocation, growth regulation, and sex determination and differentiation. In Experiment 1, cold and warm temperature treatments favored female and male proportions, respectively. The RNA-seq analysis highlighted sex-dependent energy allocation transcripts, with higher levels of nucb1 and pomc1 in future females, and increased levels of egfra and spry1 in future males. In Experiment 2, a warm thermal regime favored females, while a cold regime favored males. qPCR analysis in Experiment 2 revealed that ghrelin and nucb1 were down-regulated by warm temperatures. A significant sex-temperature interaction was observed for pank1a with higher and lower expression for males in the cold and warm regimes respectively, compared to females. Notably, spry1 displayed increased expression in future males at the all-fins stage and in males undergoing molecular sex differentiation in both experimental conditions, indicating that it provides a novel, robust, and consistent marker for masculinization. Overall, our findings emphasize the complex interplay of genes involved in feeding, energy allocation, growth, and sex determination in response to temperature variations in the European sea bass.

APIS: an updated parentage assignment software managing triploids induced from diploid parents.

In aquaculture, sterile triploids are commonly used for production as sterility gives them potential gains in growth, yields and quality. However, they cannot be reproduced, and DNA parentage assignment to their diploid or tetraploid parents is required to estimate breeding values for triploid phenotypes. No publicly available software has the ability to assign triploids to their parents. Here, we updated the R package APIS to support triploids induced from diploid parents. First, we created new exclusion and likelihood tables that account for the double allelic contribution of the dam and the recombination that can occur during female meiosis. As the effective recombination rate of each marker with the centromere is usually unknown, we set it at 0.5 and found that this value maximises the assignment rate even for markers with high or low recombination rates. The number of markers needed for a high true assignment rate did not strongly depend on the proportion of missing parental genotypes. The assignment power was however affected by the quality of the markers (minor allele frequency, call rate). Altogether, 96 to 192 SNPs were required to have a high parentage assignment rate in a real rainbow trout dataset of 1232 triploid progenies from 288 parents. The likelihood approach was more efficient than exclusion when the power of the marker set was limiting. When more markers were used, exclusion was more advantageous, with sensitivity reaching unity, very low False Discovery Rate (<0.01) and excellent specificity (0.96-0.99). Thus, APIS provides an efficient solution to assign triploids to their diploid parents.

An accurate formula to calculate exclusion power of marker sets in parentage assignment.

In studies on parentage assignment with both parents unknown, the exclusion power of a marker set is generally computed under the hypothesis that the potential families tested are independent and unrelated samples. This tends to produce overly optimistic exclusion power estimates. In this work, we have developed a new formula that gives almost unbiased results at the population level.

Effects of the total replacement of fish-based diet with plant-based diet on the hepatic transcriptome of two European sea bass (Dicentrarchus labrax) half-sibfamilies showing different growth rates with the plant-based diet.

BACKGROUND: Efforts towards utilisation of diets without fish meal (FM) or fish oil (FO) in finfish aquaculture have been being made for more than two decades. Metabolic responses to substitution of fishery products have been shown to impact growth performance and immune system of fish as well as their subsequent nutritional value, particularly in marine fish species, which exhibit low capacity for biosynthesis of long-chain poly-unsaturated fatty acids (LC-PUFA). The main objective of the present study was to analyse the effects of a plant-based diet on the hepatic transcriptome of European sea bass (Dicentrarchus labrax). RESULTS: We report the first results obtained using a transcriptomic approach on the liver of two half-sibfamilies of the European sea bass that exhibit similar growth rates when fed a fish-based diet (FD), but significantly different growth rates when fed an all-plant diet (VD). Overall gene expression was analysed using oligo DNA microarrays (GPL9663). Statistical analysis identified 582 unique annotated genes differentially expressed between groups of fish fed the two diets, 199 genes regulated by genetic factors, and 72 genes that exhibited diet-family interactions. The expression of several genes involved in the LC-PUFA and cholesterol biosynthetic pathways was found to be up-regulated in fish fed VD, suggesting a stimulation of the lipogenic pathways. No significant diet-family interaction for the regulation of LC-PUFA biosynthesis pathways could be detected by microarray analysis. This result was in agreement with LC-PUFA profiles, which were found to be similar in the flesh of the two half-sibfamilies. In addition, the combination of our transcriptomic data with an analysis of plasmatic immune parameters revealed a stimulation of complement activity associated with an immunodeficiency in the fish fed VD, and different inflammatory status between the two half-sibfamilies. Biological processes related to protein catabolism, amino acid transaminations, RNA splicing and blood coagulation were also found to be regulated by diet, while the expression of genes involved in protein and ATP synthesis differed between the half-sibfamilies. CONCLUSIONS: Overall, the combined gene expression, compositional and biochemical studies demonstrated a large panel of metabolic and physiological effects induced by total substitution of both FM and FO in the diets of European sea bass and revealed physiological characteristics associated with the two half-sibfamilies.

Parental selection for growth and early-life low stocking density increase the female-to-male ratio in European sea bass.

In European sea bass (Dicentrarchus labrax), as in many other fish species, temperature is known to influence the sex of individuals, with more males produced at relatively high temperatures. It is however unclear to what extent growth or stress are involved in such a process, since temperature is known to influence both growth rate and cortisol production. Here, we designed an experiment aiming at reducing stress and affecting early growth rate. We exposed larvae and juveniles originating from both captive and wild parents to three different treatments: low stocking density, food supplemented with tryptophan and a control. Low stocking density and tryptophan treatment respectively increased and decreased early growth rate. Each treatment influenced the stress response depending on the developmental stage, although no clear pattern regarding the whole-body cortisol concentration was found. During sex differentiation, fish in the low-density treatment exhibited lower expression of gr1, gr2, mr, and crf in the hypothalamus when compared to the control group. Fish fed tryptophan displayed lower crf in the hypothalamus and higher level of serotonin in the telencephalon compared to controls. Overall, fish kept at low density produced significantly more females than both control and fish fed tryptophan. Parents that have been selected for growth for three generations also produced significantly more females than parents of wild origin. Our findings did not allow to detect a clear effect of stress at the group level and rather point out a key role of early sexually dimorphic growth rate in sex determination.

Sex dimorphism in European sea bass (Dicentrarchus labrax L.): New insights into sex-related growth patterns during very early life stages.

The European sea bass (Dicentrarchus labrax) exhibits female-biased sexual size dimorphism (SSD) early in development. New tagging techniques provide the opportunity to monitor individual sex-related growth during the post-larval and juvenile stages. We produced an experimental population through artificial fertilization and followed a rearing-temperature protocol (~16°C from hatching to 112 days post-hatching, dph; ~20°C from 117 to 358 dph) targeting a roughly balanced sex ratio. The fish were tagged with microchips between 61 and 96 dph in five tagging trials of 50 fish each; individual standard length (SL) was recorded through repeated biometric measurements performed between 83 to 110 dph via image analyses. Body weight (BW) was modelled using the traits measured on the digital pictures (i.e. area, perimeter and volume). At 117 dph, the fish were tagged with microtags and regularly measured for SL and BW until 335 dph. The experiment ended at 358 dph with the sexing of the fish. The sex-ratio at the end of the experiment was significantly in favor of the females (65.6% vs. 34.4%). The females were significantly longer and heavier than the males from 103 dph (~30 mm SL, ~0.44 g BW) to 165 dph, but the modeling of the growth curves suggests that differences in size already existed at 83 dph. A significant difference in the daily growth coefficient (DGC) was observed only between 96 and 103 dph, suggesting a physiological or biological change occurring during this period. The female-biased SSD pattern in European sea bass is thus strongly influenced by very early growth differences between sexes, as already shown in previous studies, and in any case long before gonadal sex differentiation has been started, and thus probably before sex has been determined. This leads to the hypothesis that early growth might be a cause rather than a consequence of sex differentiation in sea bass.

Optimization of Genomic Selection to Improve Disease Resistance in Two Marine Fishes, the European Sea Bass (Dicentrarchus labrax) and the Gilthead Sea Bream (Sparus aurata).

Disease outbreaks are a major threat to the aquaculture industry, and can be controlled by selective breeding. With the development of high-throughput genotyping technologies, genomic selection may become accessible even in minor species. Training population size and marker density are among the main drivers of the prediction accuracy, which both have a high impact on the cost of genomic selection. In this study, we assessed the impact of training population size as well as marker density on the prediction accuracy of disease resistance traits in European sea bass (Dicentrarchus labrax) and gilthead sea bream (Sparus aurata). We performed a challenge to nervous necrosis virus (NNV) in two sea bass cohorts, a challenge to Vibrio harveyi in one sea bass cohort and a challenge to Photobacterium damselae subsp. piscicida in one sea bream cohort. Challenged individuals were genotyped on 57K-60K SNP chips. Markers were sampled to design virtual SNP chips of 1K, 3K, 6K, and 10K markers. Similarly, challenged individuals were randomly sampled to vary training population size from 50 to 800 individuals. The accuracy of genomic-based (GBLUP model) and pedigree-based estimated breeding values (EBV) (PBLUP model) was computed for each training population size using Monte-Carlo cross-validation. Genomic-based breeding values were also computed using the virtual chips to study the effect of marker density. For resistance to Viral Nervous Necrosis (VNN), as one major QTL was detected, the opportunity of marker-assisted selection was investigated by adding a QTL effect in both genomic and pedigree prediction models. As training population size increased, accuracy increased to reach values in range of 0.51-0.65 for full density chips. The accuracy could still increase with more individuals in the training population as the accuracy plateau was not reached. When using only the 6K density chip, accuracy reached at least 90% of that obtained with the full density chip. Adding the QTL effect increased the accuracy of the PBLUP model to values higher than the GBLUP model without the QTL effect. This work sets a framework for the practical implementation of genomic selection to improve the resistance to major diseases in European sea bass and gilthead sea bream.

Corrigendum: Optimization of Genomic Selection to Improve Disease Resistance in Two Marine Fishes, The European Sea Bass (Dicentrarchus labrax) and the Gilthead Sea Bream (Sparus aurata).

[This corrects the article DOI: 10.3389/fgene.2021.665920.].