BACKGROUND
the control of infectious diseases is a prime concern in salmon farming due to severe economic losses, reduced animal welfare and challenges to the sustainability of the industry  <cit> .
piscirickettsia salmonis is an intracellular bacterium, causing salmon rickettsial syndrome   <cit>  in salmon species and is considered one of the major pathogens of the salmon farming industry.
during p. salmonis outbreaks fish exhibit pale gills, abdominal swelling and haemorrhages on the base of fins  <cit> .
infected fish often have skin lesions  <cit> .
internally, serosanguinous ascites and swollen kidneys, livers and spleens are common.
p. salmonis is horizontally transmitted through gills or skin  <cit> .
this pathogen has evolved over time, being insidious with each outbreak and becoming refractory to treatments  <cit> .
antibiotic use may inhibit the growth of the pathogen, but treatments have been unsuccessful in stopping disease outbreaks  <cit> .
commercial vaccines have also not proven to be as efficient as needed  <cit> .
in chile alone, p. salmonis is responsible for the 74 % of the infectious-related deaths in atlantic salmon  <cit>  resulting in economic losses estimated at us $ <dig> million  <cit> .
disease resistance has become an increasingly important trait of interest included in the genetic programs of aquatic species  <cit> .
selective breeding for disease resistance in aquaculture species commonly includes survival information from experimental infection of full-sibs of the breeding candidates  <cit> .
the identification of the most resistant families is performed using the information obtained from challenge results, which can be used to identify the genomic regions  involved in the genetic variation of these traits  <cit> .
selection for disease resistance traits is based on estimated breeding values  calculated using phenotypic records and pedigree.
previous studies have estimated low to moderate heritabilities , ranging between  <dig>  and  <dig>  for resistance to p. salmonis  <cit>  in the same commercial population of atlantic salmon used in the present study; indicating a potential for selective breeding for p. salmonis resistance in this population.
the availability of molecular tools like dense panels of snps  in domesticated species has led to the search for genomic regions and causal mutations underlying variation in complex traits through the use of linkage and association mapping .
high-density snp panels have recently been developed for atlantic salmon  <cit> .
in genome-wide association studies  a group of individuals are phenotyped and then genotyped using a large number of snps, in order to detect statistical association between a marker and a trait of interest  <cit> .
data from gwas experiments are commonly analyzed testing one snp at a time using linear models  <cit> .
two limitations of this approach are that it does not account for epistasis effects  and adjustments are needed to control the false positive rate due to the multiple testing problem
the associated markers can be used to accelerate the genetic progress by means of marker assisted selection , if the total variance explained by significant markers is high  <cit> .
mas can use a causative mutation that has been identified in a gene or regulatory region that has a major effect over the trait; or can use snps that are in linkage disequilibrium  with a qtl to increase the response to selection  <cit> .
information from snp panels can be used for the evaluation of breeding candidates through genomic selection .
in gs, sibs of the selection candidates are phenotyped and genotyped, while the selection candidates are only genotyped.
the sibs of the selection candidates are utilized to predict the effect of each marker, and genomic ebvs  are calculated for the candidates using the information of all genotyped snps  <cit> .
knowledge concerning the architecture of the resistance trait is needed before deciding on the most appropriate approach to enhance artificial selection for p. salmonis resistance in atlantic salmon, .
here we present the first genome wide association analysis performed on atlantic salmon to dissect the genetic architecture of resistance to p. salmonis using a 50 k snp array in  <dig> fish.
methods
fish material and records
the origin of the population used in the present study is a chilean salmon strain derived from the irish strain fanad-mowi   <cit> .
fish from the fanad-mowi strain were introduced for farming purposes to chile during the 1990’s through commercial agreements.
in  <dig>  a breeding program for atlantic salmon was started by the company aquachile , aimed at improving economically relevant traits including growth rate, carcass quality traits and late maturation.
the base population of the breeding nucleus was established mainly with fish from the fanad-mowi strain.
more specifically, the population used in the present study corresponds to the year-class  <dig>  which already has four generations of selective breeding in chilean farming conditions.
two thousand six hundred one atlantic salmon smolts weighing on average  <dig>  g  were experimentally challenged against p. salmonis as described previously  <cit> .
briefly, infection was induced through intra-peritoneal  injection with  <dig>  ml of a ld <dig> inoculum of p. salmonis.
post-injection, infected fish were distributed in three different tanks with salt water .
these smolts were the same individuals used in a previous study  <cit>  and were the progeny of  <dig> sires and  <dig> dams from the breeding population of salmones chaicas, x region, chile.
each individual was identified using a pit  tag, and a similar number of fish from each of the  <dig> full sib families were placed in each tank.
the average number of fish per family was  <dig>  ranging from  <dig> to  <dig>  in order to determine the specific response to p. salmonis infection, the absence of other pathogens is required.
fish were found to be negative for infectious salmon anemia virus, infectious pancreatic necrosis virus and renibacterium salmoninarum by rt-pcr and negative for flavobacterium spp.
by culture.
mortalities were recorded daily and all survivors were anesthetized and euthanized at day  <dig>  final body weight was recorded at the day of death for every dead fish or at the end of the challenge for survivors.
fin samples were taken from all fish and preserved in ethanol at −80 ° c until dna extraction.
the challenge test and sampling procedures were approved by the comité de bioética animal, facultad de ciencias veterinarias y pecuarias, universidad de chile .
genotyping
genomic dna was extracted from the fin clips from all challenged fish using a commercial kit , following the manufacturer’s instructions.
fish were genotyped using a 50 k affymetrix® axiom® mydesigntm snp genotyping array designed by aquainnovo and the university of chile.
50 k genotyping array
the markers included in this array were chosen from a 200 k snp array previously developed and validated in six commercial populations  and two wild populations   <cit> .
affymetrix software snpolisher  <cit>  was used to retain  <dig>  poly-high-resolution   <cit>  and no- minor-homozygote   <cit>  snps.
snps were further filtered based on mendelian error in genotypes of  <dig> trios  and one sire-offspring pair from one commercial population.
<dig>  snps with less than  <dig> mendelian error were selected.
snps were further selected on the basis of minor allele frequency  in five commercial populations of european origin.
we gave preference to snps with higher maf in the commercial populations of european origin because the target populations which will be used for further implementation of genomic selection and genome wide association studies have european origin.
for this purpose four related european populations were grouped together as group a and a fifth population was called group b. snps were selected if they met one of the following two criteria: i) maf > = <dig>  in group a and maf > = <dig>  in group b or ii) maf > = <dig>  in group b and maf > = <dig>  in group a. from  <dig>  snps that passed the above criteria,  <dig>  snps were selected so that they are as evenly distributed along the genome as possible.
this was done by selecting snps from windows of equal size across various linkage groups of the genome.
when selecting snps from windows, higher preference was given to snps with lower mendelian error.
pearson correlation was used to estimate ld among the snp pairs belonging to same linkage group.
only one member of each snp pair with r-square > = <dig>  was retained.
for each pair higher preference was given to snps near the end of the linkage group or having lower mendelian error.
from the remaining  <dig>  candidates,  <dig>  snps were selected by removing a snp that had a distance of less than 10kb from its neighbor, so markers will be distributed equally across the genome.
genotypes were obtained for our samples following the best practice analysis workflow from affymetrix  <cit>  and selecting the poly-high-resolution and the no-minor-homozygote snps from snpolisher  <cit> .
quality control was performed on the genotypes to filter based on hardy weinberg equilibrium , minor allele frequency  and call rate for snps and samples .
for the quality control the r statistical software and the genabel library  <cit>  were used.
trait definition and genome wide association analysis
the resistance phenotypes were defined as the time to death , measured in days with values ranging from  <dig> to  <dig>  depending on the day the fish died; and as a binary survival , scored as  <dig> if the fish died during the 40-day challenge and  <dig> if the fish survived until the end of the trial.
fish weight at the end of the experiment and tank designation were also recorded to be included as covariate and factor within the model, respectively.
linear regression and logistic regression were used to identify association between snps and resistance to p. salmonis, for the td and bs traits, respectively, using the mlreg function implemented in genabel  <cit> .
logistic regression was used to account for the binary nature of the bs trait.
a genomic kinship matrix was calculated from snp data using the gkin function.
the inbreeding coefficient  was estimated using plink  <cit>  to determine the extent of inbreeding in our samples.
heritability values were also estimated from genotype data, using a linear mixed model implemented on polygenic function .
the general formula used for the linear regression model is as follows: yi=β0+β1∗snp+β2∗tk+β3∗w+ei
where yi is the phenotypic record ; β <dig> the intercept; β <dig> the effect of each snp; β <dig> is the effect of each tank ; β <dig> is the effect of the weight  and ei the random residual.
the general formula for the logistic regression model for bs is: πx=eβ0+β1∗snp+β2∗tk+β3∗w1+eβ0+β1xβ0+β1∗snp+β2∗tk+β3∗w
where πi is the probability of the random variable being one, β <dig> the intercept; β <dig> the effect of each snp; β <dig> is the effect of each tank  and β <dig> is the effect of the weight .
the -log <dig>  for each snp across the genome was plotted to summarize the gwas results.
the significance threshold was determined using false discovery rate correction  <cit> .
this correction was performed as follows:  we ordered p values p1
≤ p2
≤ … ≤ pk where k is the number of markers tested.
starting with the largest p value, we found the first individual p value  that satisfied: pi ≤ i/k * <dig> , where i = ith observation.
the pi that satisfied the condition above became the critical value for the experiment.
the proportion of the heritability explained by each significant marker was obtained by comparing the heritability estimated with the polygenic function and the heritability estimated with the inclusion of the significant snp genotype as factor  <cit> .
the proportion of the phenotypic variance explained by each marker for both traits was estimated multiplying the heritability by the proportion of heritability explained by each marker.
the approximate effect of each marker was obtained from the linear regression analysis performed in genabel  <cit> .
the level of linkage disequilibrium between significant markers was determined by calculating r <dig> statistics using the r2fast function implemented in genabel  <cit> .
for r <dig> calculations we used all genotyped samples which passed the qc.
RESULTS
the average mortality per family at the end of the test was  <dig>  %  .
we observed considerable variation in the phenotypic survival across families: some families showed more resistance with a low mortality , while other families had high mortality .fig.
<dig> cumulative mortality by family after piscirickettsia salmonis experimental infection.
the average mortality after 40 days of experimental infection against p. salmonis is plotted for each of the  <dig> full-sib families of atlantic salmon included in the present study
final body weight was recorded as a covariate in the statistical models as previously described  <cit> .
average final body weight was  <dig>  g .
fish grew during the trial, and the increased sd and cv reflects a higher level of dispersion of final body weight compared with body weight measured at the beginning of the trial .
this is likely due to varying days of death among animals, as some fish died early in the trial  and some survived the entire trial period.
forty-eight thousand eight hundred eighty-six snps and  <dig> samples passed all quality control criteria.
average f was  <dig>  .
thus, the inbreeding level in this population is considered low and will have a small impact on the power of the association analysis.
furthermore, we accounted for population structure within the sample aimed at decreasing the chances of occurrence of type i and type ii errors in the association analysis.
the heritability estimated using the genomic information was  <dig>  for td and  <dig>  for bs   <cit> , which is consistent with the results obtained in previously reported estimations using pedigree information  <cit> .
in these previous studies when the binary survival  was analysed using a threshold model, a heritability of  <dig>  ±  <dig>  was calculated  <cit> .
when a linear model was used to analyse the day of death, a heritability of  <dig>  ±  <dig>  was estimated  <cit> .
the p values obtained from the gwas indicate evidence of significant associations on atlantic salmon chromosomes ssa <dig> and ssa <dig> .fig.
<dig> gwas for piscirickettsia salmonis resistance measured as time to death: red line indicates false discovery rate significance threshold
fig.
<dig> gwas for piscirickettsia salmonis resistance measured as binary survival: red line indicates false discovery rate significance threshold
results were consistent between models after multiple testing corrections .
the linear model performed slightly better in terms of detecting significant markers associated with td than the logistic model for bs.
a linear model using td appears to be a better fit than a logistic model using bs; which leads to more efficient identification of important snps.
similarly, a slightly higher accuracy of ebvs was obtained when fitting a linear model to analyse day of death when compared to a threshold model to analyse overall test period survival  for p. salmonis resistance in atlantic salmon  <cit> .table  <dig> significant snps found for piscirickettsia salmonis resistance in s. salar detected for time to death and binary survival traits
asignificant after fdr corrections
the logistic regression model revealed two significant snp associated with bs located in one genomic region of chromosome ssa <dig> .
in the linear regression model, five significant snps associated with td were identified, the same two identified by the logistic regression and three additional snps.
one marker is located in the same genomic region of the afore-mentioned snps  and two are located in different regions of chromosomes ssa <dig> and ssa <dig> .
two of these markers are significant for both trait definitions and models used.
markers that were significant only for td were suggestive for bs.
the resistance alleles are the minor alleles, with frequencies ranging between  <dig>  and  <dig>  , indicating that there is still potential to increase their frequencies in the breeding population.table  <dig> minor allele frequency of each marker and effect of each allele on the resistance against piscirickettsia salmonis
a <dig> and a <dig> are the minor and major alleles, respectively.
mean value of the time to death by genotype for each marker.
a1: high resistance allele.
effect of the a <dig> allele and standard error in brackets
it is interesting to note that three of the identified markers are located in a ~660kb region  in chromosome ssa <dig>  this region contains several predicted genes  <cit>  and some of them could play an important role in the resistance against p. salmonis .
aqi_uch- <dig> is located in an intronic region of alpha--fucosyltransferase, which is involved in neutrophil and t cell recruitment and lymphocyte trafficking .
the snp marker located on chromosome ssa <dig> is close to a metalloprotease, described as up-regulated in atlantic salmon infected with p. salmonis  <cit> .
snp aqi_uch- <dig> is located in an intronic region of a predicted gene corresponding to tenascin, which has been described to generate an effective immune response to bacterial lipopolysaccharide  <cit> .
other genes like interleukin receptors are located in the genomic regions identified in the gwas and could be involved in differential susceptibility to p. salmonis infection.
some of these genes have been previously identified in gene expression assays of diseased fish  <cit> , and thus these results support the evidence that they may be involved in atlantic salmon’s response to p. salmonis.
further assays like differential expression or sequencing of selected genes and regions in susceptible and resistance fish could lead to a better understanding of the biological response to this pathogen.table  <dig> summary of the locations and functions of candidate genes that could play an important role in the immune response to p. salmonis infection and that have been identified in regions surrounding the most significant markers
the proportion of the heritability explained by each identified marker ranged between  <dig>  and  <dig>  for td and  <dig>  and  <dig>  for bs .
the proportion of phenotypic variance explained by each marker was small, ranging from  <dig>  to  <dig>  for bs and from  <dig>  to  <dig>  for td.table  <dig> proportion of heritability and phenotypic variance explained by each marker for time to death and binary survival traits
markers located in the same contig  exhibit high level of ld, ranging between  <dig>  and  <dig>  .
the ld between markers in this contig and the marker located in the other contig  in chromosome ssa <dig> is also relatively high , which is consistent with the physical proximity of both contigs.
the marker located on chromosome ssa <dig> showed low ld with other markers.
the closer the markers are located, the higher the value of r <dig> is to be expected  <cit> .
due to the high amount of ld present between the markers identified, their effect is likely to be shared.table  <dig> linkage disequilibrium estimates  between pairs of markers are shown above the diagonal
number of samples used to estimate r <dig> are under the diagonal
many genes usually affect quantitative traits and consequently the benefit of conducting mas depends on the proportion of the variance explained by each qtl  <cit> .
we found five markers significantly associated with resistance to p. salmonis, but the proportion of phenotypic variance and heritability explained are relatively small.
in a typical gwas analysis, markers need to surpass a specific threshold to be considered significant.
in the case of traits that are affected by many genes of small effect, their effect can be underestimated  <cit> .
on the other hand, polygenic approaches like genomic selection take into account the nature of such traits and include the information of all marker genotypes.
although there is no major effect qtl, the marker information can still be used for breeding programs  <cit>  to predict the gebvs of selection candidates  <cit>  through gs.
gs takes into account all markers to estimate the gebvs, without the need to surpass a significance threshold for association with a particular trait  <cit> .
the sib-test design can be used for mas  <cit>  or gs, where the association between markers and phenotypes is estimated in the sibs of the candidates, and the candidates are selected on breeding values that result from summing the estimates of the effects of their marker alleles  <cit> .
different approaches can be used to determine the marker effects to estimate the breeding values of the candidates.
the gs methods facilitates computation of individual breeding values for all genotyped animals and do not require any prior knowledge of the underlying qtl  <cit> .
taking the relatively small proportion of the genetic variance explained by the significant markers into account, the results presented in this study suggest that a gs approach will be the most appropriate way of incorporating molecular information to assist artificial selection for p. salmonis resistance in atlantic salmon.
CONCLUSIONS
this is the first study designed to dissect the genetic basis for the resistance for p. salmonis using a dense snp array in atlantic salmon.
resistance to p. salmonis can be described as a moderate polygenic trait, as there are likely several loci involved but each with a small effect.
five snps could be identified as significantly associated with p. salmonis resistance.
these markers explained a relatively small proportion of the variance for the trait.
availability of supporting data
snp data from the 50 k array used in this study have been deposited on salmondb database  <cit>  .
sequences flanking significant markers can be found as additional file 1: table s <dig>
additional file
additional file 1: table s <dig>  sequences flanking significant markers identified in gwas.
abbreviations
gwasgenome wide association study
ldlinkage disequilibrium
masmarker assisted selection
qtlquantitative trait loci
snpsingle nucleotide polymorphism
srssalmon rickettsial syndrome
gsgenomic selection
competing interests
the authors declare that they have no competing interests.
authors’ contributions
kc performed gwas analysis, discussion and wrote the initial version of the manuscript.
jpl contributed with study design.
mel contributed in the data analysis.
lb managed samples and performed dna extractions.
sn and nd contributed with the filtering of the 200 k array.
adg and am contributed with prediction and functional annotation of s. salar genes, markers anchoring and variant effects.
wd contributed with discussion.
jmy conceived and designed the study, supervised the work of kc and contributed to the analysis, discussion and writing.
all authors have reviewed and approved the manuscript.

 the control of infectious diseases is a prime concern in salmon farming due to severe economic losses, reduced animal welfare and challenges to the sustainability of the industry. piscirickettsia salmonis is an intracellular bacterium, causing salmon rickettsial syndrome in salmon species and is considered one of the major pathogens of the salmon farming industry. during p. salmonis outbreaks fish exhibit pale gills, abdominal swelling and haemorrhages on the base of fins. infected fish often have skin lesions. internally, serosanguinous ascites and swollen kidneys, livers and spleens are common. p. salmonis is horizontally transmitted through gills or skin. this pathogen has evolved over time, being insidious with each outbreak and becoming refractory to treatments. antibiotic use may inhibit the growth of the pathogen, but treatments have been unsuccessful in stopping disease outbreaks. commercial vaccines have also not proven to be as efficient as needed. in chile alone, p. salmonis is responsible for the 74 % of the infectious-related deaths in atlantic salmon resulting in economic losses estimated at us $ million. disease resistance has become an increasingly important trait of interest included in the genetic programs of aquatic species. selective breeding for disease resistance in aquaculture species commonly includes survival information from experimental infection of full-sibs of the breeding candidates. the identification of the most resistant families is performed using the information obtained from challenge results, which can be used to identify the genomic regions involved in the genetic variation of these traits. selection for disease resistance traits is based on estimated breeding values calculated using phenotypic records and pedigree. previous studies have estimated low to moderate heritabilities, ranging between and for resistance to p. salmonis in the same commercial population of atlantic salmon used in the present study; indicating a potential for selective breeding for p. salmonis resistance in this population. the availability of molecular tools like dense panels of snps in domesticated species has led to the search for genomic regions and causal mutations underlying variation in complex traits through the use of linkage and association mapping. high-density snp panels have recently been developed for atlantic salmon. in genome-wide association studies a group of individuals are phenotyped and then genotyped using a large number of snps, in order to detect statistical association between a marker and a trait of interest. data from gwas experiments are commonly analyzed testing one snp at a time using linear models. two limitations of this approach are that it does not account for epistasis effects and adjustments are needed to control the false positive rate due to the multiple testing problem the associated markers can be used to accelerate the genetic progress by means of marker assisted selection, if the total variance explained by significant markers is high. mas can use a causative mutation that has been identified in a gene or regulatory region that has a major effect over the trait; or can use snps that are in linkage disequilibrium with a qtl to increase the response to selection. information from snp panels can be used for the evaluation of breeding candidates through genomic selection. in gs, sibs of the selection candidates are phenotyped and genotyped, while the selection candidates are only genotyped. the sibs of the selection candidates are utilized to predict the effect of each marker, and genomic ebvs are calculated for the candidates using the information of all genotyped snps. knowledge concerning the architecture of the resistance trait is needed before deciding on the most appropriate approach to enhance artificial selection for p. salmonis resistance in atlantic salmon,. here we present the first genome wide association analysis performed on atlantic salmon to dissect the genetic architecture of resistance to p. salmonis using a 50 k snp array in fish. methods fish material and records the origin of the population used in the present study is a chilean salmon strain derived from the irish strain fanad-mowi. fish from the fanad-mowi strain were introduced for farming purposes to chile during the 1990’s through commercial agreements. in a breeding program for atlantic salmon was started by the company aquachile, aimed at improving economically relevant traits including growth rate, carcass quality traits and late maturation. the base population of the breeding nucleus was established mainly with fish from the fanad-mowi strain. more specifically, the population used in the present study corresponds to the year-class which already has four generations of selective breeding in chilean farming conditions. two thousand six hundred one atlantic salmon smolts weighing on average  g were experimentally challenged against p. salmonis as described previously. briefly, infection was induced through intra-peritoneal injection with  ml of a ld inoculum of p. salmonis. post-injection, infected fish were distributed in three different tanks with salt water. these smolts were the same individuals used in a previous study and were the progeny of sires and dams from the breeding population of salmones chaicas, x region, chile. each individual was identified using a pit tag, and a similar number of fish from each of the full sib families were placed in each tank. the average number of fish per family was ranging from to in order to determine the specific response to p. salmonis infection, the absence of other pathogens is required. fish were found to be negative for infectious salmon anemia virus, infectious pancreatic necrosis virus and renibacterium salmoninarum by rt-pcr and negative for flavobacterium spp. by culture. mortalities were recorded daily and all survivors were anesthetized and euthanized at day final body weight was recorded at the day of death for every dead fish or at the end of the challenge for survivors. fin samples were taken from all fish and preserved in ethanol at −80 ° c until dna extraction. the challenge test and sampling procedures were approved by the comité de bioética animal, facultad de ciencias veterinarias y pecuarias, universidad de chile. genotyping genomic dna was extracted from the fin clips from all challenged fish using a commercial kit, following the manufacturer’s instructions. fish were genotyped using a 50 k affymetrix® axiom® mydesigntm snp genotyping array designed by aquainnovo and the university of chile. 50 k genotyping array the markers included in this array were chosen from a 200 k snp array previously developed and validated in six commercial populations and two wild populations. affymetrix software snpolisher was used to retain poly-high-resolution and no- minor-homozygote snps. snps were further filtered based on mendelian error in genotypes of trios and one sire-offspring pair from one commercial population. snps with less than mendelian error were selected. snps were further selected on the basis of minor allele frequency in five commercial populations of european origin. we gave preference to snps with higher maf in the commercial populations of european origin because the target populations which will be used for further implementation of genomic selection and genome wide association studies have european origin. for this purpose four related european populations were grouped together as group a and a fifth population was called group b. snps were selected if they met one of the following two criteria: i) maf > = in group a and maf > = in group b or ii) maf > = in group b and maf > = in group a. from snps that passed the above criteria, snps were selected so that they are as evenly distributed along the genome as possible. this was done by selecting snps from windows of equal size across various linkage groups of the genome. when selecting snps from windows, higher preference was given to snps with lower mendelian error. pearson correlation was used to estimate ld among the snp pairs belonging to same linkage group. only one member of each snp pair with r-square > = was retained. for each pair higher preference was given to snps near the end of the linkage group or having lower mendelian error. from the remaining candidates, snps were selected by removing a snp that had a distance of less than 10kb from its neighbor, so markers will be distributed equally across the genome. genotypes were obtained for our samples following the best practice analysis workflow from affymetrix and selecting the poly-high-resolution and the no-minor-homozygote snps from snpolisher. quality control was performed on the genotypes to filter based on hardy weinberg equilibrium, minor allele frequency and call rate for snps and samples. for the quality control the r statistical software and the genabel library were used. trait definition and genome wide association analysis the resistance phenotypes were defined as the time to death, measured in days with values ranging from to depending on the day the fish died; and as a binary survival, scored as if the fish died during the 40-day challenge and if the fish survived until the end of the trial. fish weight at the end of the experiment and tank designation were also recorded to be included as covariate and factor within the model, respectively. linear regression and logistic regression were used to identify association between snps and resistance to p. salmonis, for the td and bs traits, respectively, using the mlreg function implemented in genabel. logistic regression was used to account for the binary nature of the bs trait. a genomic kinship matrix was calculated from snp data using the gkin function. the inbreeding coefficient was estimated using plink to determine the extent of inbreeding in our samples. heritability values were also estimated from genotype data, using a linear mixed model implemented on polygenic function. the general formula used for the linear regression model is as follows: yi=β0+β1∗snp+β2∗tk+β3∗w+ei where yi is the phenotypic record; β the intercept; β the effect of each snp; β is the effect of each tank; β is the effect of the weight and ei the random residual. the general formula for the logistic regression model for bs is: πx=eβ0+β1∗snp+β2∗tk+β3∗w1+eβ0+β1xβ0+β1∗snp+β2∗tk+β3∗w where πi is the probability of the random variable being one, β the intercept; β the effect of each snp; β is the effect of each tank and β is the effect of the weight. the -log for each snp across the genome was plotted to summarize the gwas results. the significance threshold was determined using false discovery rate correction. this correction was performed as follows: we ordered p values p1  ≤ p2  ≤ … ≤ pk where k is the number of markers tested. starting with the largest p value, we found the first individual p value that satisfied: pi ≤ i/k *, where i = ith observation. the pi that satisfied the condition above became the critical value for the experiment. the proportion of the heritability explained by each significant marker was obtained by comparing the heritability estimated with the polygenic function and the heritability estimated with the inclusion of the significant snp genotype as factor. the proportion of the phenotypic variance explained by each marker for both traits was estimated multiplying the heritability by the proportion of heritability explained by each marker. the approximate effect of each marker was obtained from the linear regression analysis performed in genabel. the level of linkage disequilibrium between significant markers was determined by calculating r statistics using the r2fast function implemented in genabel. for r calculations we used all genotyped samples which passed the qc.  the average mortality per family at the end of the test was  %. we observed considerable variation in the phenotypic survival across families: some families showed more resistance with a low mortality, while other families had high mortality.fig. cumulative mortality by family after piscirickettsia salmonis experimental infection. the average mortality after 40 days of experimental infection against p. salmonis is plotted for each of the full-sib families of atlantic salmon included in the present study final body weight was recorded as a covariate in the statistical models as previously described. average final body weight was  g. fish grew during the trial, and the increased sd and cv reflects a higher level of dispersion of final body weight compared with body weight measured at the beginning of the trial. this is likely due to varying days of death among animals, as some fish died early in the trial and some survived the entire trial period. forty-eight thousand eight hundred eighty-six snps and samples passed all quality control criteria. average f was. thus, the inbreeding level in this population is considered low and will have a small impact on the power of the association analysis. furthermore, we accounted for population structure within the sample aimed at decreasing the chances of occurrence of type i and type ii errors in the association analysis. the heritability estimated using the genomic information was for td and for bs, which is consistent with the results obtained in previously reported estimations using pedigree information. in these previous studies when the binary survival was analysed using a threshold model, a heritability of  ±  was calculated. when a linear model was used to analyse the day of death, a heritability of  ±  was estimated. the p values obtained from the gwas indicate evidence of significant associations on atlantic salmon chromosomes ssa and ssa.fig. gwas for piscirickettsia salmonis resistance measured as time to death: red line indicates false discovery rate significance threshold fig. gwas for piscirickettsia salmonis resistance measured as binary survival: red line indicates false discovery rate significance threshold results were consistent between models after multiple testing corrections. the linear model performed slightly better in terms of detecting significant markers associated with td than the logistic model for bs. a linear model using td appears to be a better fit than a logistic model using bs; which leads to more efficient identification of important snps. similarly, a slightly higher accuracy of ebvs was obtained when fitting a linear model to analyse day of death when compared to a threshold model to analyse overall test period survival for p. salmonis resistance in atlantic salmon.table significant snps found for piscirickettsia salmonis resistance in s. salar detected for time to death and binary survival traits asignificant after fdr corrections the logistic regression model revealed two significant snp associated with bs located in one genomic region of chromosome ssa. in the linear regression model, five significant snps associated with td were identified, the same two identified by the logistic regression and three additional snps. one marker is located in the same genomic region of the afore-mentioned snps and two are located in different regions of chromosomes ssa and ssa. two of these markers are significant for both trait definitions and models used. markers that were significant only for td were suggestive for bs. the resistance alleles are the minor alleles, with frequencies ranging between and, indicating that there is still potential to increase their frequencies in the breeding population.table minor allele frequency of each marker and effect of each allele on the resistance against piscirickettsia salmonis a and a are the minor and major alleles, respectively. mean value of the time to death by genotype for each marker. a1: high resistance allele. effect of the a allele and standard error in brackets it is interesting to note that three of the identified markers are located in a ~660kb region in chromosome ssa this region contains several predicted genes and some of them could play an important role in the resistance against p. salmonis. aqi_uch- is located in an intronic region of alpha--fucosyltransferase, which is involved in neutrophil and t cell recruitment and lymphocyte trafficking. the snp marker located on chromosome ssa is close to a metalloprotease, described as up-regulated in atlantic salmon infected with p. salmonis. snp aqi_uch- is located in an intronic region of a predicted gene corresponding to tenascin, which has been described to generate an effective immune response to bacterial lipopolysaccharide. other genes like interleukin receptors are located in the genomic regions identified in the gwas and could be involved in differential susceptibility to p. salmonis infection. some of these genes have been previously identified in gene expression assays of diseased fish, and thus these results support the evidence that they may be involved in atlantic salmon’s response to p. salmonis. further assays like differential expression or sequencing of selected genes and regions in susceptible and resistance fish could lead to a better understanding of the biological response to this pathogen.table summary of the locations and functions of candidate genes that could play an important role in the immune response to p. salmonis infection and that have been identified in regions surrounding the most significant markers the proportion of the heritability explained by each identified marker ranged between and for td and and for bs. the proportion of phenotypic variance explained by each marker was small, ranging from to for bs and from to for td.table proportion of heritability and phenotypic variance explained by each marker for time to death and binary survival traits markers located in the same contig exhibit high level of ld, ranging between and. the ld between markers in this contig and the marker located in the other contig in chromosome ssa is also relatively high, which is consistent with the physical proximity of both contigs. the marker located on chromosome ssa showed low ld with other markers. the closer the markers are located, the higher the value of r is to be expected. due to the high amount of ld present between the markers identified, their effect is likely to be shared.table linkage disequilibrium estimates between pairs of markers are shown above the diagonal number of samples used to estimate r are under the diagonal many genes usually affect quantitative traits and consequently the benefit of conducting mas depends on the proportion of the variance explained by each qtl. we found five markers significantly associated with resistance to p. salmonis, but the proportion of phenotypic variance and heritability explained are relatively small. in a typical gwas analysis, markers need to surpass a specific threshold to be considered significant. in the case of traits that are affected by many genes of small effect, their effect can be underestimated. on the other hand, polygenic approaches like genomic selection take into account the nature of such traits and include the information of all marker genotypes. although there is no major effect qtl, the marker information can still be used for breeding programs to predict the gebvs of selection candidates through gs. gs takes into account all markers to estimate the gebvs, without the need to surpass a significance threshold for association with a particular trait. the sib-test design can be used for mas or gs, where the association between markers and phenotypes is estimated in the sibs of the candidates, and the candidates are selected on breeding values that result from summing the estimates of the effects of their marker alleles. different approaches can be used to determine the marker effects to estimate the breeding values of the candidates. the gs methods facilitates computation of individual breeding values for all genotyped animals and do not require any prior knowledge of the underlying qtl. taking the relatively small proportion of the genetic variance explained by the significant markers into account, the results presented in this study suggest that a gs approach will be the most appropriate way of incorporating molecular information to assist artificial selection for p. salmonis resistance in atlantic salmon.  this is the first study designed to dissect the genetic basis for the resistance for p. salmonis using a dense snp array in atlantic salmon. resistance to p. salmonis can be described as a moderate polygenic trait, as there are likely several loci involved but each with a small effect. five snps could be identified as significantly associated with p. salmonis resistance. these markers explained a relatively small proportion of the variance for the trait. availability of supporting data snp data from the 50 k array used in this study have been deposited on salmondb database. sequences flanking significant markers can be found as additional file 1: table s additional file additional file 1: table s sequences flanking significant markers identified in gwas. abbreviations gwasgenome wide association study ldlinkage disequilibrium masmarker assisted selection qtlquantitative trait loci snpsingle nucleotide polymorphism srssalmon rickettsial syndrome gsgenomic selection competing interests the authors declare that they have no competing interests. authors’ contributions kc performed gwas analysis, discussion and wrote the initial version of the manuscript. jpl contributed with study design. mel contributed in the data analysis. lb managed samples and performed dna extractions. sn and nd contributed with the filtering of the 200 k array. adg and am contributed with prediction and functional annotation of s. salar genes, markers anchoring and variant effects. wd contributed with discussion. jmy conceived and designed the study, supervised the work of kc and contributed to the analysis, discussion and writing. all authors have reviewed and approved the manuscript. 
candidate genes including interleukin receptors and fucosyltransferase have been found to be physically linked with these genetic markers and may play an important role in the differential immune response against this pathogen. knowledge concerning the architecture of the resistance trait is needed before deciding on the most appropriate approach to enhance artificial selection for p. salmonis resistance in atlantic salmon. linear regression and logistic regression models were tested. atlantic salmon smolts were experimentally challenged against p. salmonis by means of intra-peritoneal injection. pisciricketssia salmonis is the causal agent of salmon rickettsial syndrome, which affects salmon species and causes severe economic losses. the purpose of the study was to dissect the genetic variation in the resistance to this pathogen in atlantic salmon. due to the small amount of variance explained by each significant marker we conclude that genetic resistance to this pathogen can be more efficiently improved with the implementation of genetic evaluations incorporating genotype information from a dense snp array. a genome wide association analysis was performed on data from the challenged fish. genome wide association analysis indicated that resistance to p. salmonis is a moderately polygenic trait. salmon rickettsial syndromegenome wide association analysispathogen resistanceatlantic salmonsingle nucleotide polymorphismissue-copyright-statement© the author 2015 the proportion of the phenotypic variance explained by each marker is small, ranging from to. selective breeding for disease resistance represents one approach for controlling srs in farmed atlantic salmon.
 pisciricketssia salmonis is the causal agent of salmon rickettsial syndrome, which affects salmon species and causes severe economic losses. selective breeding for disease resistance represents one approach for controlling srs in farmed atlantic salmon. knowledge concerning the architecture of the resistance trait is needed before deciding on the most appropriate approach to enhance artificial selection for p. salmonis resistance in atlantic salmon. the purpose of the study was to dissect the genetic variation in the resistance to this pathogen in atlantic salmon. methods atlantic salmon smolts were experimentally challenged against p. salmonis by means of intra-peritoneal injection. these smolts were the progeny of sires and dams from a chilean breeding population. mortalities were recorded daily and the experiment ended at day post-inoculation. fish were genotyped using a 50k affymetrix® axiom® mydesigntm single nucleotide polymorphism genotyping array. a genome wide association analysis was performed on data from the challenged fish. linear regression and logistic regression models were tested.  genome wide association analysis indicated that resistance to p. salmonis is a moderately polygenic trait. there were five snps in chromosomes ssa and ssa significantly associated with the traits analysed. the proportion of the phenotypic variance explained by each marker is small, ranging from to. candidate genes including interleukin receptors and fucosyltransferase have been found to be physically linked with these genetic markers and may play an important role in the differential immune response against this pathogen.  due to the small amount of variance explained by each significant marker we conclude that genetic resistance to this pathogen can be more efficiently improved with the implementation of genetic evaluations incorporating genotype information from a dense snp array. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords salmon rickettsial syndromegenome wide association analysispathogen resistanceatlantic salmonsingle nucleotide polymorphismissue-copyright-statement© the author 2015 