Skip to main content

QTL detection for Aeromonas salmonicida resistance related traits in turbot (Scophthalmus maximus)



Interactions between fish and pathogens, that may be harmless under natural conditions, often result in serious diseases in aquaculture systems. This is especially important due to the fact that the strains used in aquaculture are derived from wild strains that may not have had enough time to adapt to new disease pressures. The turbot is one of the most promising European aquaculture species. Furunculosis, caused by the bacterium Aeromonas salmonicida, produces important losses to turbot industry. An appealing solution is to achieve more robust broodstock, which can prevent or diminish the devastating effects of epizooties. Genomics strategies have been developed in turbot to look for candidate genes for resistance to furunculosis and a genetic map with appropriate density to screen for genomic associations has been also constructed. In the present study, a genome scan for QTL affecting resistance and survival to A. salmonicida in four turbot families was carried out. The objectives were to identify consistent QTL using different statistical approaches (linear regression and maximum likelihood) and to locate the tightest associated markers for their application in genetic breeding strategies.


Significant QTL for resistance were identified by the linear regression method in three linkage groups (LGs 4, 6 and 9) and for survival in two LGs (6 and 9). The maximum likelihood methodology identified QTL in three LGs (5, 6 and 9) for both traits. Significant association between disease traits and genotypes was detected for several markers, some of them explaining up to 17% of the phenotypic variance. We also identified candidate genes located in the detected QTL using data from previously mapped markers.


Several regions controlling resistance to A. salmonicida in turbot have been detected. The observed concordance between different statistical methods at particular linkage groups gives consistency to our results. The detected associated markers could be useful for genetic breeding strategies. A finer mapping will be necessary at the detected QTL intervals to narrow associations and around the closely associated markers to look for candidate genes through comparative genomics or positional cloning strategies. The identification of associated variants at specific genes will be essential, together with the QTL associations detected in this study, for future marker assisted selection programs.


In aquaculture, disease resistance related traits are of particular importance. The reason is that interactions between fish and pathogens, that may be harmless under natural conditions, often result in disease problems in aquaculture systems because of the added stress from biological, physical and chemical factors [1]. This is especially important due to the fact that, in contrast to farm animals, the strains used in aquaculture usually have been very recently derived from wild strains [2] and, therefore, have had little time to adapt to the new disease pressures within the aquaculture environment.

Improvements in the performance of any productive or physiological characteristic of the cultured species can be achieved (if the trait is genetically determined) through artificial selection. Several breeding programmes have been developed for different traits in aquaculture species usually involving growth rate, cold tolerance and disease resistance. However, implementing classical breeding programs focused on disease resistance traits could be highly problematic, since the phenotypic measurement of these traits is often complex and expensive. It is not possible to evaluate specimens to be selected and thus, evaluation has to be performed on relatives. Moreover, the required challenges may cause animal suffering and increase risks of infection at farm facilities [3]. Despite those disadvantages, some populations of Atlantic salmon have been already selected for resistance to bacterial and viral diseases. Selection for resistance to Infectious Pancreatic Necrosis (IPN) virus, based on bath challenge tests of several hundred families of first feeding fry, showed 66.6% and 29.3% mortality for low and high resistant strains, respectively [4]. Also, high resistance to IPN in rainbow trout was achieved by selective breeding [5]. In that study, the commercial strain RT-201 artificially challenged with IPN virus showed a mortality of 4.3%, whereas the highly sensitive controls reached 96.1%. In carps, Schäperclaus [6] found resistance to the dropsy disease, selected lines suffering low mortality (11.5%) well below the unselected ones (57%).

Positive response to selection pressures is possible because resistance against particular diseases affecting aquaculture species often shows moderate to high heritabilities and, thus, there is a large potential for genetic improvement. For example recent heritability estimates for resistance to Aeromonas salmonicida ranged from 0.43 to 0.62 in Atlantic salmon [reviewed by [7]], and it was 0.51 ± 0.03 in Salvelinus fontinalis[8].

Turbot (Scophthalmus maximus) is a flatfish that has been intensively cultured during the last decade due to its great commercial value. Its production in Europe has increased from 3000 Tm in 1996 to 9246 Tm in 2009 [9]. Increasing growth rate, controlling sex ratio (females largely outgrow males) and enhancing disease resistance currently constitute the main goals of genetic breeding programmes in this species. Pathologies constitute one of the main problems of turbot culture. Among these, furunculosis, caused by Aeromonas salmonicida, has produced important losses to turbot industry [10, 11]. Genomic resources of turbot have increased in the last years [12, 13] and an immune-enriched oligo-microarray was designed and applied to identify candidate genes of resistance to A. salmonicida[14]. Also, a microsatellite consensus genetic map including centromere positions was reported in this species [15, 16], and recently, 31 EST-linked microsatellites, particularly useful for comparative genomics, were added to the consensus map [17]. Combining functional genomics strategies with the detection of genomic regions associated to productive characters (using genetic maps) increases the power to identify genes involved in the phenotypic differences occurring within and between families.

Besides an infinitesimal component (due to small effects of a huge number of loci) the variation in quantitative traits may be also controlled by a few genes with larger effects. Genomic regions closely linked to those genes show association with the trait phenotype and are known as quantitative trait loci (QTL) [18]. The effects of allele segregation at molecular markers throughout the genome can be used to determine the number and position of trait-related QTL, as well as the magnitude of their effects [19]. If, eventually, the responsible gene of a large effect on a trait is detected and the causal mutation determined, selection could be exerted directly on the genotype for that locus (Gene Assisted Selection, GAS). Alternatively, genetic maps provide DNA markers tightly linked to genes affecting different traits. Such markers can be used in Marker-Assisted Selection (MAS), selection based partly or fully on DNA marker genotypes. Consequently, disease resistance traits are candidates for the implementation of MAS and, especially, GAS programs, which would allow the evaluation of individuals without exposing them to the pathogen or relying on relatives' information alone.

There have been a certain number of studies on disease resistance QTL in the main aquaculture species, especially focused in rainbow trout and Atlantic salmon within fish. In rainbow trout, several QTL were detected for resistance to IPN and IHN (Infectious Hematopoietic Necrosis) viruses [2023] and to different parasites, including Ceratomyxa shasta and Myxobolus cerebralis[24, 25]. The identified QTL in response to IHN virus were detected in more than one family supporting their consistency [22], and in the case of resistance to Myxobolus cerebralis, a very strong association was detected which explained between 50 and 86% of the phenotypic variance across families [25]. In Atlantic salmon, QTL for disease resistance have been reported for the ISA (Infectious Salmon Anaemia) virus [26], for the bacterium Aeromonas salmonicida responsible of furunculosis [27] and for the ectoparasite Gyrodactylus salaries[28]. Of particular relevance in this species was a recent study on ten full-sib families for mapping QTL for resistance against the IPN virus in post-smolts of Scottish origin, based on data from a field trial [29]. In this study, a major QTL explained up to 21% of the phenotypic variation in the data set and was found to segregate in 7 out of 20 parents investigated. Additionally, this QTL mapped to the same location of a recently detected QTL for IPN-resistance that explained 29% of the phenotypic variance using ten large full-sib families of challenge-tested Norwegian Atlantic salmon. This particular QTL was found to be segregating in 10 out of 20 parents, and a subsequent fine-mapping with additional markers narrowed the QTL peak to a 4 cM region on linkage group 21 [3]. This QTL, detected in two different populations, is now being implemented in within-family selection in both Scotland [30] and Norway [3]. QTL for disease resistance have also been reported in a few cases in non-salmonid teleosts, including those for stress/immune response in tilapia [31], for resistance to pasteurellosis in gilthead sea bream [32] and for resistance to lymphocystis viral disease in Japanese flounder [33]. Other QTL for disease resistance traits in aquaculture species were identified in Crassostrea virginica, Ostrea edulis and Paralichthys olivaceus [see 7 for a review].

To date, only two studies identified QTL in turbot. One detected a QTL for body length highly associated with the marker YSKr51, explaining 12.4% of the phenotypic variance [34], and the other one identified a significant sex-determining QTL highly associated with the SmaUSC-E30 marker, which allowed for correct sex determination in up to 98.4% of the studied individuals [35]. The construction of a genetic map with appropriate marker density is necessary to detect QTL controlling quantitative traits of economic interest in aquaculture [20]. However, genetic linkage maps of most aquaculture species have only recently been available. Bouza et al. [17] reported an updated consensus including a total of 273 microsatellites clustered in 26 linkage groups, comprising 1343.2 cM length, with an average distance between markers of 6.5 ± 0.5 cM.

In the present work, a genome scan for QTL affecting resistance and survival to A. salmonicida in four turbot families was carried out using the reported microsatellite panel. The objectives were: (1) to locate QTL on the available linkage map, (2) to compare QTL obtained through the use of two different methodologies: linear regression and maximum likelihood, and (3) to determine the association between markers and traits.



Four full-sib turbot families were used to identify QTL for resistance and survival to A. salmonicida and to evaluate the association of those traits with the genotyped markers. Three of these families were obtained from the Stolt Sea Farm S.A. breeding program (FamAS-1, FamAS-2 and FamAS-3) and one family was obtained from Insuiña S.A. (FamAP). Both companies have their facilities located in NW Spain.

Since no selected strains for high and low disease resistance exist in turbot, families were selected trying to emphasize genetic divergence between parents to detect as much as possible allelic variants associated with resistance. Thus, families founded with unrelated grandparents and, when possible, from different Atlantic origins, were chosen. Also, we selected families where a three-generation pedigree was available. This enabled us to know the linkage-phase between markers for a more consistent statistical analysis. This process finally led to four independent full-sib pair analyses.

Trait measurement

All offspring of each family (approximately 150 individuals per family) were intracelomically injected at the age of four months (mean weight around 20 g) with a highly virulent A. salmonicida strain. Although A. salmonicida infection was previously reported to last around 21 days [11], the experiment was prolonged to 39 days to increase phenotypic variance for survival. After injection, two disease resistance related traits were evaluated: resistance (Re) and survival (Su). Around one hundred individuals for each family were selected to evaluate both traits (the 50 most resistant and the 50 most sensitive individuals). The dichotomous trait Re was defined as the survival vs non survival status of individuals at the end of the experiment (i. e. day 39 after the injection). The trait Su was defined as the number of days one individual survived (i.e. elapsed days until the individual died or the experiment finished). Consequently, this is a censored trait since all individuals still alive at the end of the study were scored with the same value for survival.

Genetic map

The panel of markers used for QTL identification was reported by Martínez et al. [35] and it is based on the consensus map by Bouza et al. [15], the new EST-linked microsatellites by Bouza et al. [17] and the centromere mapping [16].

Table 1 shows the number of analyzed microsatellites, the number of linkage groups, the map length, the average distance between microsatellites and the average number of microsatellites per linkage group. The average distance between markers in the map ranged between 15.53 and 16.06 cM, being below the minimum distance proposed for QTL detection (< 20 cM) [36].

Table 1 Screening figures in the four families analyzed for QTL identification

QTL analyses

Two programs were used to detect QTLs: GridQTL [37] and QTLMap [38]. Within both methodologies, two approaches were followed. First, a single QTL was assumed at each LG. Afterwards, a two-QTL model was also tested within each LG.


This software implements a linear regression (LR) methodology, considering the linkage phase between markers according to pedigree information. The default regression method [39] was applied, and the genome and chromosome-wide significance thresholds were estimated by implementing bootstrapping at p = 0.05 and 0.01 [40], with a permutation test of 10, 000 iterations [41].


This software detects QTL through interval mapping using a maximum likelihood (ML) test. To determine the significance level 10, 000 simulations were performed for each trait and LG, with a heritability set to 0.10 [42].

For both methodologies, an outbred full-sib model was used and a QTL was considered suggestive when significance was between 5% and 1% at chromosome-wide level, and significant when significance was below 1% at chromosome-wide level or when significance was below 5% at genome-wide level [43]. These thresholds also allow establishing a confidence interval to allocate the detected QTL.

Sex was not included into the model because it was not available. Although there is some debate about the influence of growth-related traits on the resistance to diseases (see for example [44] for positive results, and [45] and [46] for absence of correlations between both types of traits), weight and length were included as covariates within the model to reduce the stochasticity.

Association analyses

A one way analysis of variance (ANOVA) was performed on the phenotypic values (resistance and survival) of the progeny for each family using individual genotypes from markers within the LGs where a significant or suggestive QTL was found. The objective was to detect associations between markers and traits by estimating the between-genotype component of the observed phenotypic variance (i.e. differences attributable to the different marker genotypes). To avoid false positives due to multiple testing, a simple Bonferroni correction was performed for all tests involving those markers within the same LG. Each ANOVA also provided a corrected R2 value useful to estimate the reduction of the overall phenotypic variance of the trait due to the model fitting, thus providing the proportion of the trait variance predictable from the given marker genotypes.

Results and Discussion

Trait values

Table 2 shows the number of observations and the mean for the traits recorded (resistance and survival) in the four families of the experiment. The standard deviation ranged between 13.12 - 15.31 days for survival and 0.44 - 0.50 for resistance. Notice that a greater power for detection is expected in challenges where there is approximately a 50% of survival (LD50), as it occurred in AS families. Although a lower survival (26%) was observed for family AP, this did not preclude the detection of QTL in that family (see Table 3). Mean weight and length ranged between 31.62 - 46.45 g and 10.70 - 12.83 cm, respectively.

Table 2 Statistics of the measured traits in the four families analyzed (± standard deviation)
Table 3 Location of the QTL detected for resistance and survival to A. salmonicida with two different statistical methodologies

QTL analyses

Table 3 shows the location of the detected QTL (LG, estimated position and interval) with both methodologies (LR and ML). There were no QTL at genome-wide level with the LR approach, but two suggestive (LG12 and LG13) and three significant (LG4, LG6 and LG9) QTL were detected for resistance, and four suggestive (LG4, LG11, LG12 and LG13) and two significant (LG6 and LG9) QTL were detected for survival at chromosome-wide level. Using the ML methodology, one QTL for resistance (LG6) and three for survival (LG5, LG6 and LG9) were identified at genome-wide level with significance below 1%. Additionally, three suggestive (LG4, LG12 and LG16) and two significant (LG5 and LG9) QTL for resistance, and five suggestive (LG1, LG4, LG12, LG16 and LG18) QTL for survival were detected at chromosome level.

In our study, four full-sib families from segregating populations were evaluated separately using two different methodologies. Kao [47] investigated the performances of LR and ML methods for QTL detection, suggesting the application of LR as an initial procedure to obtain preliminary results and then use the ML method as a final procedure in order to obtain the most conclusive results. Our data support in part Kao's suggestion, since the highest significance levels were yielded by ML for the same QTL detected with both methodologies. The concordance in QTL location between both methods or between different families (LG16 in FamAS-2 and FamAP; see Table 3) increases the confidence of the results obtained.

QTL were less significant when length and weight were excluded from the model (data not shown) indicating that body size may contribute to or may be correlated with another trait affecting disease resistance. The regression coefficients for the covariates themselves were not significant except for body weight in one of the families (FamAP). Overturf et al. [44] indicated a positive correlation between body size and disease resistance related traits, but other authors detected no significant correlation between both types of traits [45, 46]. Although the inclusion of body size improved QTL detection in our study, the exact relationship with resistance to A. salmonicida is not known.

Classical QTL mapping methods assume that traits follow a normal distribution [48]. However, the categorical (dead or alive) and survival (length of life) data used in this study are not-normally distributed. For such traits, classical QTL detection methods could have a lower power and a bias when estimating the effects and position of QTL. Nevertheless, it has been proposed that classical interval mapping methods using a Gaussian model on censored data, analyzed as if they were uncensored, have not a reduced accuracy on QTL location and estimation of QTL effects [49].

As can be seen in Table 3, several QTL were detected for resistance and survival at close positions within the same LG. One reason for this agreement could be the high relationship expected between both traits, since the resistance concept also includes survival (resistant individuals survive for a long period, in fact the whole experiment). Other reason could be that the same genes were involved in the mechanisms underlying both disease resistance traits.

The power to detect QTL depends on the heritability of the trait, the recombination distance between QTL and markers, the proportion of phenotypic variance explained by the QTL, the QTL allele frequency and the sample size [50]. In addition, the power of QTL analysis is limited, since only QTL segregating in one or both parents can be detected. In fact, there were no QTL detected within FamAS-1. Considering that families come from unrelated and genetically divergent grandparents from natural populations of the Atlantic area, identified QTL in our study could be representative of the genetic architecture of disease resistance related traits in turbot. Anyway, the detected QTL should be verified in other turbot families before their use in breeding programs.

The probability of success and the accuracy of results of QTL studies depend on the precision and density of the available genetic map. Following this idea, the map used in this study is being further refined and, thus, the QTL detection could also be improved. Currently, Martínez et al. (unpublished data) are developing a consensus map comprising 463 markers (microsatellites and EST) with a moderately dense coverage.

Association analyses

All markers within the same linkage group where a suggestive or significant QTL was detected were analyzed. Table 4 shows the results only for markers significantly associated with any of the traits. At least one significantly associated marker was found in all linkage groups where a QTL was detected, except for LGs 1, 12, 13 and 18. However, within these LGs only suggestive QTL were detected (see Table 3). Remarkably, not always the most significantly associated markers were the closest to the estimated position of the QTL (data not shown). The reasons could be: (1) low information content of the closest marker due to the parents genotype; (2) a large extension of the area in linkage disequilibrium with the detected QTL, which could result in positive associations between the trait and the genotypes at several markers; (3) the existence of a secondary segregating QTL, although this situation was not detected when LGs were analyzed for the possibility of carrying two QTL; and finally (4) marker positions in the map are not definitive, which could modify the order and the distance between markers. For example, it should be noticed that the present genetic map of turbot has more linkage groups (26) than chromosomes observed in this species (2n = 44 chromosomes) [51].

Table 4 Proportion of the phenotypic variance explained by the markers significantly associated with the evaluated traits

Markers showing significant association with studied traits explained from 7 to 17% of the phenotypic variance. The marker SmaUSC-E30, significantly associated with resistance and survival traits at LG5, was previously associated to the major QTL for sex determination in turbot [35]. This may suggest an influence of sex in the measured traits that could be related either to the reported growth rate differences between sexes in turbot [52] or to sex-susceptibility/immune differences [53]. A relationship between growth and immune function was previously reported in fish [54, 55]. In the case of LG9, up to three markers showed significant association with both traits (Table 4: SmaUSC-E23, SmaUSC-E41 and Sma-USC21), the two first ones being functionally annotated as immune-related genes [14].

A high explained phenotypic variance was also observed for an IPN-resistance QTL in salmonids 21% [29] and 29% [3]. A QTL for ISA-resistance, also in Atlantic salmon, explained 6 - 9% of the phenotypic variance [26]. In rainbow trout two QTL for IPN-resistance were identified explaining 27 and 34% of the phenotypic variance, respectively [20].

Several associations have been detected in teleosts between disease resistance and candidate genes, mostly involving polymorphism of the Major Histocompatibility Complex (MHC), an essential gene family for adaptive immunity [7]. In rainbow and cutthroat trout, an association was detected with the resistance to IHN virus [56]; in rainbow trout, a suggestive association with the resistance to Bacterial Cold Water Disease (BCWD) [57]; and in Atlantic salmon, associations were detected for resistance to A. salmonicida[27, 58, 59], IHN [60], and ISA [27]. Beyond salmonids, associations between MHC markers and disease resistance have also been reported in other important farmed teleosts like Japanese flounder (resistance to Vibrio anguillarum) [61]; common carp (resistance to Cyprinid Herpesvirus-3) [62], and in the turbot (resistance to Edwardsiella tarda) [63]. Several of these associations involved the MHC class IIB gene, an observation which strongly recommends the mapping of this gene in turbot to cross this information against the QTL locations from this study and to evaluate its possible role in general resistance to bacterial pathogens. Nevertheless, direct selection on MHC variants could be not advisable since balancing selection [57] or overdominance [64] have been suggested to be acting on this polymorphism. The risk of pathogen-specific allelic selection could determine undesired loss of MHC variation, which may increase susceptibility to other pathogens. Furthermore, it has been highlighted that other genes may also play an important role to explain genetic variance for disease resistance [7]. Recently, significant associations with other relevant immune-genes have been reported in common carp (Interleukin 10 (IL-10) with resistance to cyprinid herpesvirus-3[65]) and in grass carp (Toll-like Receptor-3 (TLR3) with resistance to reovirus [66]).

A revision and a refinement of the mapping regions where QTL were located should be performed to further narrow the intervals and facilitate the advance on candidate gene strategies. Also, mapping of differentially expressed (DE) genes identified in response to A. salmonicida infections [14] should be evaluated or acquired to cross this information with QTL position. None of the immune-related genes significantly associated with resistance and survival traits in this study (SmaUSC-E23, SmaUSC-E41) were found to be regulated in response to A. salmonicida infection [14]. Nevertheless, two DE genes detected in response to this pathogen [14] were located in the vicinity of two QTL at LG9 and LG11, respectively. SmaUSC-E41, linked to the immune-related DE ecdysone receptor A (BLAST annotation; E-value: 3E-8), was located at LG9 close to the marker SmaUSC-E5 [4 cM; see [17]]. This gene pertains to a nuclear receptor superfamily present in all Metazoa, containing ligand dependent transcription factors related to immune system and regulation of inflammatory processes [67]. On the other hand, Sma-USC147 was 7 cM apart from the significantly associated SmaUSC-E24 marker at LG11. This marker is linked to a DE gene (annotated as PRA1 family protein 3 using BLAST; E-value: 3E-9), which modulates antiapoptotic activity and immune response [68]. Since these DE markers were not included in the panel of markers for QTL screening in this work, further association analysis across the same QTL families should be carried out. Finally, markers linked to relevant DE genes in the turbot and previously associated with disease resistance in other fish (e.g., TLR3, MHC class IA and MHC IIA) [14, 66, 7], should be mapped. This would enable to test their association within families to provide new candidates for a functional explanation of resistance to A. salmonicida in the turbot.

If the responsible gene and the causal mutation were detected, then, GAS could be performed. Increasing map density should facilitate the detection of more closely linked markers useful to implement MAS. For both strategies, another promising approach is the use of comparative genomics [69]. A study of comparative mapping by Bouza et al. [15] identified syntenic relationships between sequences of turbot and the model fish Tetraodon nigroviridis (Tni): between turbot LG16 and Tni chromosome 19; between turbot LG5 and Tni1; and between turbot LG6 and Tni13. Since in the present study QTL at LG5, LG6 and LG16 were detected, the analysis of genes located at these syntenic blocks in T. nigroviridis through comparative genomics tools with the turbot DE genes in response to A. salmonicida[14] could provide relevant information about candidates underlying the detected QTL.

Finfish aquaculture is still relatively in its first stages in the large-scale animal production sector. However, it can take advantage on the vast amounts of genomic data currently being generated for livestock, as well as for model fish organisms. Many livestock genomics programs have focused on identifying genes influencing similar economic traits as those important in aquaculture, such as growth rate, disease resistance and meat quality traits. Through comparative genomic approaches, sequence information from other vertebrates can be used to quickly isolate homologous genes in finfish. It could be also interesting to extend this study to other pathogens for detection of disease resistance related QTL affecting turbot aquaculture systems in order to investigate more general disease resistance related traits.


Regions controlling resistance-related traits to Aeromonas salmonicida in turbot were detected. Concordance on the detection of QTL between statistical methods gives more consistency to our results. A finer mapping will be necessary on those linkage groups with effect on disease resistance by increasing the density of markers on the position where the QTL were detected. In addition, functional and comparative genomics strategies could be used to confirm the detected QTL and to look for candidate genes. Alternatively, the associated markers could be useful in programs for marker-assisted selection in turbot.


  1. 1.

    Wedemeyer G: Physiology of fish in intensive culture systems. 1996, London, Chapman & Hall

    Google Scholar 

  2. 2.

    Duarte CM, Marbá N, Holmer M: Rapid domestication of marine species. Science. 2007, 316: 383-383. 10.1126/science.1142201.

    Article  Google Scholar 

  3. 3.

    Moen T, Baranski M, Sonesson AK, Kjøglum S: Confirmation and fine-mapping of a major QTL for resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar): population-level associations between markers and trait. BMC Genomics. 2009, 10: 368-10.1186/1471-2164-10-368.

    Article  PubMed  PubMed Central  Google Scholar 

  4. 4.

    Storset A, Strand C, Wetten M, Kjøglum S, Ramstad A: Response to selection for resistance against infectious pancreatic necrosis in Atlantic salmon (Salmo salar L.). Aquaculture. 2007, 272: 62-68.

    Article  Google Scholar 

  5. 5.

    Okamoto N, Tayama T, Kawanobe M, Fujiki N, Yasuda Y, Sano T: Resistance of a rainbow trout strain to infectious pancreatic necrosis. Aquaculture. 1993, 117: 71-76. 10.1016/0044-8486(93)90124-H.

    Article  Google Scholar 

  6. 6.

    Schäperclaus W: Traité de pisciculture en étang. 1962, Paris, Vigot Frères

    Google Scholar 

  7. 7.

    Ødegård J, Baranski M, Gjerde B, Gjedrem T: Methodology for genetic evaluation of disease resistance in aquaculture species: challenges and future prospects. Aquaculture Research. 2011, 42: 103-114.

    Article  Google Scholar 

  8. 8.

    Perry GML, Tarte P, Croisetière S, Belhumeur P, Bernatchez L: Genetic variance and covariance for 0 + brook charr (Salvelinus fontinalis) weight and survival time of furunculosis (Aeromonas salmonicida) exposure. Aquaculture. 2004, 235: 263-271. 10.1016/j.aquaculture.2004.03.002.

    Article  Google Scholar 

  9. 9.

    APROMAR: Informe anual sobre el desarrollo de la acuicultura marina en España y el mundo. 2010, España

    Google Scholar 

  10. 10.

    Ellis AE: Immunization with bacterial antigens: furunculosis. Fish Vaccinology. Edited by: Gudding R, Lillehaug A, Midtlyng PJ, Brown PR. 1997, Karger, Basel, 107-116.

    Google Scholar 

  11. 11.

    Austin B, Austin DA: Bacterial fish pathogens. Diseases of farmed and wild fish. 1999, United Kingdom, Springer-Praxis Publishing Ltd

    Google Scholar 

  12. 12.

    Pardo BG, Fernández C, Millán A, Bouza C, Vázquez-López A, Vera M, Álvarez-Dios JA, Calaza M, Gómez-Tato A, Vázquez M, Cabaleiro S, Magariños B, Lemos ML, Leiro JM, Martínez P: Expressed sequence tags (ESTs) from immune tissues of turbot (Scophthalmus maximus) challenged with pathogens. BMC Vet Res. 2008, 4: 37-10.1186/1746-6148-4-37.

    Article  PubMed  PubMed Central  Google Scholar 

  13. 13.

    Vera M, Álvarez-Dios JA, Millán A, Pardo BG, Bouza C, Hermida M, Fernández C, de la Herrán R, Molina-Luzón J, Martínez P: Validation of single nucleotide polymorphism (SNP) markers from an immune expressed sequence tag (EST) turbot, Scophthalmus maximus, database. Aquaculture. 2011, 313: 31-41. 10.1016/j.aquaculture.2011.01.038.

    CAS  Article  Google Scholar 

  14. 14.

    Millán A, Gómez-Tato A, Fernández C, Pardo BG, Álvarez-Dios JA, Calaza M, Bouza C, Vázquez M, Cabaleiro S, Martínez P: Design and performance of a turbot (Scophthalmus maximus) oligo-microarray based on ESTs from immune tissues. Mar Biotechnol. 2010, 12: 452-465. 10.1007/s10126-009-9231-0.

    Article  PubMed  Google Scholar 

  15. 15.

    Bouza C, Hermida M, Pardo BG, Fernández C, Castro J, Fortes GG, Sánchez L, Presa P, Pérez M, Sanjuán A, Comesaña S, Álvarez-Dios JA, Calaza M, Cal R, Piferrer F, Martínez P: A microsatellite genetic map in the turbot (Scophthalmus maximus). Genetics. 2007, 177: 2457-2467. 10.1534/genetics.107.075416.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Martínez P, Hermida M, Pardo BG, Fernández C, Castro J: Centromere-linkage in the turbot (Scophthalmus maximus) through half-tetrad analysis in diploid meiogynogenetics. Aquaculture. 2008, 280: 81-88. 10.1016/j.aquaculture.2008.05.011.

    Article  Google Scholar 

  17. 17.

    Bouza C, Hermida M, Millán A, Vilas R, Vera M: Characterization of EST- derived microsatellites for gene mapping and evolutionary genomics in turbot. Anim Genet. 2008, 39: 666-670. 10.1111/j.1365-2052.2008.01784.x.

    CAS  Article  PubMed  Google Scholar 

  18. 18.

    Mackay TFC: The genetic architecture of quantitative traits. Ann Rev Genet. 2001, 35: 303-339. 10.1146/annurev.genet.35.102401.090633.

    CAS  Article  PubMed  Google Scholar 

  19. 19.

    Lynch M, Walsh B: Genetics and analysis of quantitative traits. 1998, Sunderland, Sinauer Associates

    Google Scholar 

  20. 20.

    Ozaki A, Sakamoto T, Khoo S, Nakamura K, Coimbra MR, Akutsu T, Okamoto N: Quantitative trait loci (QTLs) associated with resistance/susceptibility to infectious pancreatic necrosis virus (IPNV) in rainbow trout (Oncorhynchus mykiss). Mol Genet Genomics. 2001, 265: 23-31. 10.1007/s004380000392.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Palti Y, Parsons JE, Thorgaard GH: Identification of candidate DNA markers associated with IHN virus resistance in backcrosses of rainbow (Oncorhynchus mykiss) and cutthroat trout (O clarki). Aquaculture. 1999, 173: 81-94. 10.1016/S0044-8486(98)00471-2.

    CAS  Article  Google Scholar 

  22. 22.

    Rodríguez MF, Lapatra S, Williams S, Famula T, May B: Genetic markers associated with resistance to infectious hematopoietic necrosis in rainbow and steelhead trout. Aquaculture. 2004, 241: 93-115. 10.1016/j.aquaculture.2004.08.003.

    Article  Google Scholar 

  23. 23.

    Barroso RM, Wheeler PA, Lapatra SE, Drew RE, Thorgaard GH: QTL for IHNV resistance and growth identified in a rainbow (Oncorhynchus mykiss) × Yellowstone cutthroat (Oncorhynchus clarki bouvieri) trout cross. Aquaculture. 2008, 277: 156-163. 10.1016/j.aquaculture.2008.03.001.

    CAS  Article  Google Scholar 

  24. 24.

    Nichols KM, Bartholomew J, Thorgaard GH: Mapping multiple genetic loci associated with Ceratomyxa shasta resistance in Oncorhynchus mykiss. Dis Aquat Organ. 2003, 56: 145-54.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Baerwald MR, Petersen JL, Hedrick RP, Schisler GJ, May B: A major effect quantitative trait locus for whirling disease resistance identified in rainbow trout (Oncorhynchus mykiss). Heredity. 2010, 106: 920-926.

    Article  PubMed  PubMed Central  Google Scholar 

  26. 26.

    Moen T, Sonesson AK, Hayes B, Lien S, Munck H, Meuwissen THE: Mapping of a quantitative trait locus for resistance against infectious salmon anaemia in Atlantic salmon (Salmo Salar): comparing survival analysis with analysis on affected/resistant data. BMC Genetics. 2007, 8: 53-

    Article  PubMed  PubMed Central  Google Scholar 

  27. 27.

    Grimholt U, Larsen S, Nordmo R, Midtlyng P, Kjoeglum S, Storset A, Saebo S, Stet RJM: MHC polymorphism and disease resistance in Atlantic salmon (Salmo salar); facing pathogens with single expressed major histocompatibility class I and class II loci. Immunogenetics. 2003, 55: 210-219. 10.1007/s00251-003-0567-8.

    CAS  Article  PubMed  Google Scholar 

  28. 28.

    Gilbey J, Verspoor E, Mo TA, Sterud E, Olstad K, Hytterød S, Noble LR, Jones CS: Identification of genetic markers associated with Gyrodactylus salaris resistance in Atlantic salmon (Salmo salar L). Diseases of Aquatic Organisms. 2006, 71: 119-129.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Houston RD, Haley CS, Hamilton A, Guy DR, Tinch AE, Taggart JB, Mcandrew BJ, Bishop SC: Major quantitative traits loci affect resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar). Genetics. 2008, 178: 1109-1115. 10.1534/genetics.107.082974.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Houston RD, Haley CS, Hamilton A, Guy DR, Mota-Velasco JC, Gheyas AA, Tinch AE, Taggart JB, Bron JE, McAndrew BJ, Verner-Jeffreys DW, Paley RK, Rimmer GSE, Tew IJ, Bishop SC: The susceptibility of Atlantic salmon fry to freshwater infectious pancreatic necrosis is largely explained by a major QTL. Heredity. 2009, 105: 318-327.

    Article  PubMed  Google Scholar 

  31. 31.

    Cnaani A, Zilberman N, Tinman S, Hulata G, Ron M: Genome-scan analysis for quantitative trait loci in an F2 tilapia hybrid. Mol Genet Genomics. 2004, 272: 162-172.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    Massault C, Franch R, Haley C, de Koning DJ, Bovenhyis H, Pellizzan C, Patarnello T, Bargelloni L: Quantitative trait loci for resistance to fish pasteurellosis in gilthead sea bream (Sparus aurata). Animal Genetics. 2011, 42: 191-203. 10.1111/j.1365-2052.2010.02110.x.

    CAS  Article  PubMed  Google Scholar 

  33. 33.

    Fuji K, Kobayashi K, Hasegawa O, Coimbra MR, Sakamoto T, Okamoto N: Identification of a single major genetic locus controlling the resistance to lymphocystis disease in Japanese flounder (Paralichthys olivaceus). Aquaculture. 2006, 254: 203-210. 10.1016/j.aquaculture.2005.11.024.

    CAS  Article  Google Scholar 

  34. 34.

    Ruan XH, Wang WJ, Kong J, Yu F, Huang XQ: Genetic linkage mapping of turbot (Scophthalmus maximus L.) using microsatellite markers and its application in QTL analysis. Aquaculture. 2010, 308: 89-100. 10.1016/j.aquaculture.2010.08.010.

    CAS  Article  Google Scholar 

  35. 35.

    Martínez P, Bouza C, Hermida M, Fernández J, Toro MA, Vera M, Pardo BG, Millán A, Fernández C, Vilas R, Viñas A, Sánchez L, Felip A, Piferrer F, Ferreiro I, Cabaleiro S: Identification of the major sex-determining region of turbot (Scophthalmus maximus). Genetics. 2009, 183: 1443-52. 10.1534/genetics.109.107979.

    Article  PubMed  PubMed Central  Google Scholar 

  36. 36.

    Dekkers JCM, Hospital F: The use of molecular genetics in the improvement of agricultural populations. Nat Rev Genet. 2002, 3: 22-32. 10.1038/nrg701.

    CAS  Article  PubMed  Google Scholar 

  37. 37.

    Seaton G, Hernández J, Grunchec JA, White I, Allen J, De Koning DJ, Wei W, Berry D, Haley C, Knott S: GridQTL: A grid portal for QTL mapping of compute intensive datasets. Proceedings of the 8th World Congress on Genetics Applied to Livestock Production. 2006, Belo Horizonte, Brazil

    Google Scholar 

  38. 38.

    Gilbert H, Le Roy P, Moreno C, Robelin D, Elsen JM: QTLMAP, a software for QTL detection in outbred population. Annals of Human Genetics. 2008, 72 (5): 694-

    Google Scholar 

  39. 39.

    Haseman JK, Elston RC: The investigation of linkage between a quantitative trait and a marker locus. Behav Genet. 1972, 2: 3-19. 10.1007/BF01066731.

    CAS  Article  PubMed  Google Scholar 

  40. 40.

    Visscher P, Thompson R, Haley C: Confidence intervals in QTL mapping by bootstrapping. Genetics. 1996, 143: 1013-1020.

    CAS  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. Genetics. 1994, 138: 963-971.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. 42.

    Antonello J, Massault C, Francha R, Haley CS, Pellizzari C, Bovo G, Atarnello Th, de Koning DJ, Bargelloni L: Estimates of heritability and genetic correlation for body length and resistance to fish pasteurellosis in the gilthead sea bream (Sparus aurata L.). Aquaculture. 2009, 298: 29-35. 10.1016/j.aquaculture.2009.10.022.

    Article  Google Scholar 

  43. 43.

    Lander E, Kruglyak L: Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. Nature Genet. 1995, 11: 241-247. 10.1038/ng1195-241.

    CAS  Article  PubMed  Google Scholar 

  44. 44.

    Overturf K, Lapatra S, Towner R, Campbell N, Narum S: Relationships between growth and disease resistance in rainbow trout, Oncorhynchus mykiss (Walbaum). J Fish Dis. 2010, 33: 321-329. 10.1111/j.1365-2761.2009.01124.x.

    CAS  Article  PubMed  Google Scholar 

  45. 45.

    Henryon M, Jokumsenb A, Berga P, Lundb I, Pedersenb PB, Olesenc NJ, Slierendrechtc WJ: Genetic variation for growth rate, feed conversion efficiency, and disease resistance exists within a farmed population of rainbow trout. Aquaculture. 2002, 209: 59-76. 10.1016/S0044-8486(01)00729-3.

    Article  Google Scholar 

  46. 46.

    Silverstein JT, Vallejo RL, Palti Y, Leeds TD, Rexroad CD, Welch TJ, Wiens GD, Ducrocq V: Rainbow trout resistance to bacterial cold-water disease is moderately heritable and is not adversely correlated with growth. J Anim Sci. 2009, 87: 860-867.

    CAS  Article  PubMed  Google Scholar 

  47. 47.

    Kao CH: On the differences between maximum likelihood and regression interval mapping in the analysis of quantitative trait loci. Genetics. 2000, 156: 855-865.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. 48.

    Knott SA, Elsen JM, Haley CS: Methods for multiple marker mapping of quantitative trait loci in half-sib populations. Theor Appl Genet. 1996, 93: 71-80. 10.1007/BF00225729.

    CAS  Article  PubMed  Google Scholar 

  49. 49.

    Moreno CR, Elsen JM, Le Roy P, Ducrocq V: Interval mapping methods for detecting QTL affecting survival and time-to-event phenotypes. Genet Res Camb. 2005, 85: 139-149. 10.1017/S0016672305007366.

    CAS  Article  Google Scholar 

  50. 50.

    Massault C, Bovenhuis H, de Koning DJ: QTL mapping designs for aquaculture. Aquaculture. 2008, 285: 23-29. 10.1016/j.aquaculture.2008.06.040.

    CAS  Article  Google Scholar 

  51. 51.

    Bouza C, Sánchez L, Martínez P: Karyotypic characterisation of turbot (Scophthalmus maximus) with conventional fluorochrome and restriction endonuclease banding techniques. Mar Biol. 1994, 120: 609-613. 10.1007/BF00350082.

    Article  Google Scholar 

  52. 52.

    Imsland AK, Folkvord A, Grung GL, Stefansson SO, Taranger GL: Sexual dimorphism in growth and maturation of turbot, Scophthalmus maximus (Rafinesque, 1810). Aquaculture Research. 1997, 28: 101-114. 10.1111/j.1365-2109.1997.tb01022.x.

    Article  Google Scholar 

  53. 53.

    Cuesta A, Vargas-Chacoff L, García-López A, Arjona FJ, Martínez-Rodríguez G, Meseguer J, Mancera JM, Esteban MA: Effect of sex-steroid hormones, testosterone and estradiol, on humoral immune parameters of gilthead seabream. Fish & Shellfish Immunology. 2007, 23: 693-700. 10.1016/j.fsi.2007.01.015.

    CAS  Article  Google Scholar 

  54. 54.

    Peterson BC, Bosworth BG, Bilodeau AL: Differential expression of IGF-I, IGF-II, and toll-like receptors 3 and 5 mRNA during embryogenesis in hybrid (channel X blue) and channel catfish. Comparative Biochemistry and Physiology Part A. 2005, 141: 42-47. 10.1016/j.cbpb.2005.03.009.

    Article  Google Scholar 

  55. 55.

    Johansen K, Sealey W, Overturf KE: The effects of chronic immune stimulation on muscle growth in rainbow trout. Comparative Biochemistry and Physiology Part B. 2006, 144: 520-531. 10.1016/j.cbpb.2006.05.006.

    Article  Google Scholar 

  56. 56.

    Palti Y, Nichols KM, Waller KI, Parsons JE, Thorgaard GH: Association between DNA polymorphisms tightly linked to MHC class II genes and IHN virus resistance in backcrosses of rainbow and cutthroat trout. Aquaculture. 2001, 194: 283-289. 10.1016/S0044-8486(00)00526-3.

    CAS  Article  Google Scholar 

  57. 57.

    Johnson NA, Vallejo RL, Silverstein JT, Welch TJ, Wiens GD, Hallerman EM, Palti Y: Suggestive association of major histocompatibility IB genetic markers with resistance to bacterial cold water disease in rainbow trout (Oncorhynchus mykiss). Mar Biotechnol. 2008, 10: 429-437. 10.1007/s10126-007-9080-7.

    CAS  Article  PubMed  Google Scholar 

  58. 58.

    Langefors A, Lohm J, Grahn M, Andersen O, von Schantz T: Association between major histocompatibility complex class IIB alleles and resistance to Aeromonas salmonicida in Atlantic salmon. Proc R Soc Lond B Biol Sci. 2001, 268: 479-485. 10.1098/rspb.2000.1378.

    CAS  Article  Google Scholar 

  59. 59.

    Lohm J, Grahn M, Langefors Å, Andersen Ø, Storset A, von Schantz T: Experimental evidence for major histocompatibility complex allele-specific resistance to a bacterial infection. Proc R Soc Lond B Biol Sci. 2002, 269: 2029-2033. 10.1098/rspb.2002.2114.

    CAS  Article  Google Scholar 

  60. 60.

    Miller KM, Winton JR, Schulze AD, Purcell MK, Ming TJ: Major histocompatibility complex loci are associated with susceptibility of Atlantic salmon to infectious hematopoietic necrosis virus. Environ Biol Fishes. 2004, 69: 307-316.

    Article  Google Scholar 

  61. 61.

    Xu T, Chen S, Zhang Y: MHC class IIα gene polymorphism and its association with resistance/susceptibility to Vibrio anguillarum in Japanese flounder (Paralichthys olivaceus). Developmental & Comparative Immunology. 2010, 34: 1042-1050. 10.1016/j.dci.2010.05.008.

    CAS  Article  Google Scholar 

  62. 62.

    Rakus KL, Wiegertjes GF, Adamek M, Siwicki AK, Lepa A, Irnazarow I: Resistance of common carp (Cyprinus carpio L.) to Cyprinid herpesvirus-3 is influenced by major histocompatibility (MH) class II B gene polymorphism. Fish & Shellfish Immunology. 2009, 26: 737-743. 10.1016/j.fsi.2009.03.001.

    CAS  Article  Google Scholar 

  63. 63.

    Xu JY, Chen SL, Ding H: Specific MHC class II B alleles associated with resistance to Edwardsiella tarda in turbot, Psetta maxima (L.). J Fish Dis. 2009, 32: 637-640. 10.1111/j.1365-2761.2009.01028.x.

    CAS  Article  PubMed  Google Scholar 

  64. 64.

    Turner SM, Faisal M, DeWoody JA: Zygosity at the major histocompatibility class IIB locus predicts susceptibility to Renibacterium salmoninarum in Atlantic salmon (Salmo salar L.). Animal Genetics. 2007, 38: 517-519. 10.1111/j.1365-2052.2007.01631.x.

    CAS  Article  PubMed  Google Scholar 

  65. 65.

    Kongchum P, Sandel E, Lutzky S, Hallerman EM, Hulata G, David L, Palti Y: Association between IL-10a SNPs and resistance to Cyprinid herpesvirus-3 infection in common carp (Cyprinus carpio). Aquaculture. 2011, 315: 417-421. 10.1016/j.aquaculture.2011.02.035.

    CAS  Article  Google Scholar 

  66. 66.

    Heng J, Su J, Huang T, Dong J, Chen L: The polymorphism and haplotype of TLR3 gene in grass carp (Ctenopharyngodon idella) and their associations with susceptibility/resistance to grass carp reovirus. Fish & Shellfish Immunology. 2011, 30: 45-50. 10.1016/j.fsi.2010.09.004.

    CAS  Article  Google Scholar 

  67. 67.

    Howarth DL, Law SHW, Barnes B, Hall JM, Hinton DE, Moore L, Maglich JM, Moore JT, Kullman SW: Paralogous vitamin D receptors in teleosts: transition of nuclear receptor function. Endocrinology. 2008, 149: 2411-2422. 10.1210/en.2007-1256.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  68. 68.

    Li LY, Shih HM, Liu MY, Chen JY: The cellular protein PRA1 modulates the anti-apoptotic activity of Epstein-Barr virus BHRF1, a homologue of Bcl-2, through direct interaction. The Journal of Biological Chemistry. 2001, 276: 27354-27362. 10.1074/jbc.M103821200.

    CAS  Article  PubMed  Google Scholar 

  69. 69.

    Li J, Boroevich KA, Koop BF, Davidson WS: Comparative genomics identifies candidate genes for Infectious Salmon Anemia (ISA) resistance in Atlantic salmon (Salmo salar). Mar Biotechnol. 2011, 13: 232-241. 10.1007/s10126-010-9284-0.

    CAS  Article  PubMed  Google Scholar 

Download references


We thank two anonymous referees for helpful comments on the manuscript. This study was supported by the Spanish government projects (Consolider Ingenio Aquagenomics: CSD200700002; MCI: CGL2009-13278-C02 and AGL2009-13273), and Xunta de Galicia and Fondos Feder (Grupos de Referencia Competitiva, 2010/80). B. G. P. was supported by an Isidro Parga Pondal research fellowship from Xunta de Galicia (Spain).

Author information



Corresponding author

Correspondence to Jesús Fernández.

Additional information

Authors' contributions

STRR was responsible of the final version of QTL and marker analysis and the writing of the article. MAT and JF were responsible of the supervision of the QTL and marker analysis and also of the supervision of the article. CB, MH and BGP performed genotyping at all families and contributed to discussion on association analysis. SC was responsible of challenges of fish and measuring of phenotypic traits. PM was responsible of the work and supervised all tasks. JF was responsible for a first version of the QTL analysis. All authors read and approved the final manuscript.

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Rodríguez-Ramilo, S.T., Toro, M.A., Bouza, C. et al. QTL detection for Aeromonas salmonicida resistance related traits in turbot (Scophthalmus maximus). BMC Genomics 12, 541 (2011).

Download citation


  • Quantitative Trait Locus
  • Atlantic Salmon
  • Differentially Express
  • Differentially Express Gene
  • Quantitative Trait Locus Detection