Email updates

Keep up to date with the latest news and content from BMC Veterinary Research and BioMed Central.

Open Access Research article

Risk map and spatial determinants of pancreas disease in the marine phase of Norwegian Atlantic salmon farming sites

Saraya Tavornpanich1*, Mathilde Paul2, Hildegunn Viljugrein13, David Abrial4, Daniel Jimenez1 and Edgar Brun1

Author Affiliations

1 Section for Epidemiology, Department of Health Surveillance, Norwegian Veterinary Institute, Ullevålsveien 68, Pb 750 Sentrum, N-0106, Oslo, Norway

2 UR AGIRs, Centre de coopération internationale en recherche agronomique pour le développement (CIRAD), TA C22/E, Campus international de Baillarguet, 34398, Montpellier cedex 5, France

3 Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biology, University of Oslo, P.O. Box 1066, Blindern, N-0316, Oslo, Norway

4 UR 346, Institut National de la Recherche Agronomique (INRA), 63122, Saint-Genès-Champanelle, France

For all author emails, please log on.

BMC Veterinary Research 2012, 8:172  doi:10.1186/1746-6148-8-172


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1746-6148/8/172


Received:29 June 2012
Accepted:19 September 2012
Published:24 September 2012

© 2012 Tavornpanich et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Outbreaks of pancreas disease (PD) greatly contribute to economic losses due to high mortality, control measures, interrupted production cycles, reduced feed conversion and flesh quality in the aquaculture industries in European salmon-producing countries. The overall objective of this study was to evaluate an effect of potential factors contributing to PD occurrence accounting for spatial congruity of neighboring infected sites, and then create quantitative risk maps for predicting PD occurrence. The study population included active Atlantic salmon farming sites located in the coastal area of 6 southern counties of Norway (where most of PD outbreaks have been reported so far) from 1 January 2009 to 31 December 2010.

Results

Using a Bayesian modeling approach, with and without spatial component, the final model included site latitude, site density, PD history, and local biomass density. Clearly, the PD infected sites were spatially clustered; however, the cluster was well explained by the covariates of the final model. Based on the final model, we produced a map presenting the predicted probability of the PD occurrence in the southern part of Norway. Subsequently, the predictive capacity of the final model was validated by comparing the predicted probabilities with the observed PD outbreaks in 2011.

Conclusions

The framework of the study could be applied for spatial studies of other infectious aquatic animal diseases.

Keywords:
Pancreas disease; Aquatic epidemiology; Spatial analysis; Disease mapping; Bayesian modeling

Background

Pancreas disease (PD) is a viral disease affecting Atlantic salmon (Salmon salar L.,) and rainbow trout (Oncorhynchus mykiss). A necessary cause of the disease is salmon pancreas disease virus, an alphavirus often referred as salmonid alphavirus (SAV). The onset of a PD outbreak is favored by unspecified influential environmental factors. Diagnostic criteria for PD outbreak were described by Taksdal et al., 2007 [1]. So far, 6 subtypes of SAV have been reported worldwide, whereas, only subtypes 2 and 3 have been detected in Norwegian salmonid farming sites [2, Hilde Sindre per communications]. The disease has been a significant problem in salmonid farming industries in Scotland, Ireland, and Norway [2]. An economic modeling study estimated that the direct cost of PD accounting for biological losses, disease controls and treatments, and insurance payment in a Norwegian salmonid site with an approximate of 500,000 smolts to be at least 14.4 million NOK in 2007 [3].

The characteristics of SAV and PD outbreaks have been thoroughly reviewed [4-11]. In Norway, PD outbreaks have been reported only during the marine phase. The disease is regarded as endemic in the area along the western coast of Norway from Rogaland county to a firewall at Hustadvika (Figure 1) in Møre and Romsal with few sporadic outbreaks occurring north of this border. This firewall is a 10 nautical mile exposed costal area without any aquaculture activity. Veterinary authorities together with salmon industries have launched different prevention and control strategies to control the situation in the south and to keep the northern area free of PD [12]. All findings of epidemiological studies of PD outbreaks support the hypothesis that the virus can spread through the water media, causing a PD outbreak both with and without human intervention. This aspect of the disease spreading via water media is a crucial element of aquatic infectious diseases [9,13]. Therefore, it is worthwhile to evaluate the spatial clustering of infected sites in epidemiological studies of PD outbreaks. In this study, we initially evaluated potential risk factors of PD in marine phase of Atlantic salmon farming sites, accounting for spatial congruity of infected sites. Then, the predicted probabilities of PD occurrence based on the model results were displayed in an interpolated map covering the southern part of Norwegian marine aquaculture.

thumbnailFigure 1. Geographical distribution of 359 active Atlantic salmon farming sites located near 6 counties of Norway (Vest-Agder, Rogaland, Hordaland, Sogn & Fjordane, Møre & Romsdal, SørTrøndelag) from 1 January 2009 to 31 December 2010.

Methods

Data

In this study, we included Atlantic salmon farming sites located in 6 Norwegian counties (Vest-Agder, Rogaland, Hordaland, Sogn & Fjordane, Møre & Romsdal, SørTrøndelag), where most of the PD outbreaks were reported so far. To be considered as a PD at-risk site, the included sites were required to rear fish on the site for at least 6 consecutive months during the study period (1 January 2009 to 31 December 2010). All operators of Norwegian salmonid farming industries are required to register production statistical data to responsible authorities on a monthly basis. These statistics are linked to a farm site identity, which is geo-referenced in the Directorate of Fisheries (www. fiskeridir.no). According to the inclusion criteria, 359 Atlantic salmon farming sites were qualified for model estimation of this study. The geographical distribution of the study sites is presented in Figure 1. In addition, a second dataset of active farming sites during a period of 2011 was created for validating the predictive capacity of the model. The second dataset contained 418 farming sites located throughout the study area.

PD infection status

All farming sites were monitored by the farmer for unusual increases in fish mortality and activities (e.g. runt, moribund). If PD was suspected, samples were routinely sent for laboratory confirmation using histopathological examination and PCR testing. The study sites were classified as PD positive if at least 1 PD outbreak were confirmed during the study period, and PD negative, if otherwise.

Site latitude

In the present study, we evaluated if the site location, based on the latitude, had an influence on the PD occurrence. In this study, the latitudes of the study sites cover the area of 63°39'5,481"N (the farthest north) to 58°12'33,449"N (the farthest south). The value of site latitude was converted to meters prior to analyses.

Site density

Site density for a given site was the number of farming sites (excluding the given site) that located within 10-km seaway distance of the given site. In this study, the site density ranged from 1 to14, had mean of 4, and a median of 4.5.

Smolt cohort

A smolt cohort is the group of smolts that were put to sea and reared in sea cages until slaughtering. We identified the time when the smolt cohorts were put to sea. If the cohort was put to sea in the period of August-February it was classified as ‘a cohort of autumn smolt’; otherwise, if the cohort was put to sea in the period of March-July, it was classified as ‘a cohort of spring smolt’. Of the 359 study sites, 168 (46.8%) sites had only autumn cohorts; 170 (47.4%) sites had only spring cohorts; and 21 (5.8%) sites had mixed cohorts (both cohorts of autumn and spring smolts).

PD history in 2008

Previous recent history of PD in a nearby area is likely to have an impact on a future probability of PD occurrence. In this study, we evaluated this effect by classifying the study sites based on whether or not they were located within 10-km seaway distance of sites where PD was suspected or confirmed in 2008. The 10-km seaway distance was used as the cutoff because it was estimated that sites locating within 10-km seaway distance of PD infected sites had a large impact on PD risk [14]. Of 359 study sites, 225 (63%) sites were located within 10-km radius of PD suspected/confirmed sites in 2008. Furthermore, among these 225 sites, 110 (49%) sites were classified as PD confirmed cases during the study period.

Local biomass density

The local biomass density (LBD) indicates the potential of infection endangerment at a given site. This factor depends on average biomass of the surrounding sites to a given site (excluding LBD of the given site). The LBD data were estimated for all registered salmon fish farming sites in Norway and the estimation method was described by Jansen et al (2012) [15]. In the present study, we evaluated the effect of LBD on the probability of PD occurrence at each farming site. The LBD of the study sites ranged from 4.87e + 02 to 4.62e + 08; had a mean of 2.19e + 08; and median of 2.11e + 08. For the analysis purpose (facilitating convergence of the model), the LBD was transformed using log10, and then back transformed for interpretation.

Statistical analysis

Univariate and multivariate analysis

All regression analyses were performed using WinBUGS (windows version of Bayesian inference using Gibbs Sampling) open source software [16]. First, we performed univariate logistic regression to individually evaluate the effect of each potential risk factor on PD occurrence. Records of PD confirmed cases from 1 January 2009 to 31 December 2010 were used as a dependent variable of the model. The evaluated risk factors were site latitude, site density within 10-km seaway distance, smolt cohort, PD history in 2008, and LBD. We explored the assumption of linear relationships with generalized additive models, using the mgcv package in R software. Variables with a 80% posterior probability interval including zero were excluded from further analyses. Then, multivariate logistic regression models were performed with and without a spatial component. A stepwise selection based on the Deviance Information Criterion (DIC) was used for selecting the combination of the covariates; the model with lower DIC is preferred to the model with larger DIC [17]. DIC was also used for evaluating whether adding spatial random effect would improve model fit in the final model. Residuals of all models were in turn checked for spatial dependency, as any pattern present in the residuals of a statistical model may indicate that results are biased [18]. To do so, we applied a Monte Carlo method [19,20] using the geoRglm package of the R software. This consists in comparing the observed variogram with variogram ‘envelopes’ that were computed by simulating 999 permutations of the data values across locations.

Spatial dependency modeling

We grouped fish farming sites into hexagonal subunits which were included in the model as a spatial area-level random effect [21]. Nineteen hexagons (width of 80 km, corresponding to 5546.6 km2) were created; they contained between 1 and 52 farming sites (mean = 35). The area-level effect was modeled using a Conditional Autoregressive (CAR) prior structure, in which an adjacency matrix was specified with a weight of 1 given to adjacent hexagons and a weight of 0 given to nonadjacent hexagons.

Briefly, each site in this study was labeled with i = 1…19 (the hexagon in which the site is localised), and j = 1…359 (the unique site identifier). The observed case of PD outbreaks at each site (y[ij]) during the study period followed a Bernoulli distribution. The model can be written as the following equations:

<a onClick="popup('http://www.biomedcentral.com/1746-6148/8/172/mathml/M1','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1746-6148/8/172/mathml/M1">View MathML</a>

(1)

<a onClick="popup('http://www.biomedcentral.com/1746-6148/8/172/mathml/M2','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1746-6148/8/172/mathml/M2">View MathML</a>

(2)

<a onClick="popup('http://www.biomedcentral.com/1746-6148/8/172/mathml/M3','MathML',630,470);return false;" target="_blank" href="http://www.biomedcentral.com/1746-6148/8/172/mathml/M3">View MathML</a>

(3)

β is the vector of coefficients, x[ij] is the covariates matrix observed on site j, v[i] is the random area-level effect which is spatially structured according to the model of Besag et al (1991) [22]. The mean of this model is driven by the average risk level of the adjacent hexagon set δ[i] of the hexagon i. The precision (σ2) of v[i] follow a gamma distribution. Priors for model coefficients were based on a normal distribution (with a mean of 0.5 and a precision of 5e-04). We implemented the models in WinBUGS, and estimated the model parameters using Gibbs sampling, which is the method for simulating a new value for the model parameters from its full posterior conditional distribution, given the current values for the remaining parameters [23]. In this study, the model parameters were estimated based on 45000 effective iterations. An additional first 5000 iterations were discarded as a burn-in period. Geweke and Heidelberger-Welch tests were used to assess the convergence of the models [16].

Interpolated map of PD predicted probability

The median value (50th) of the posterior distribution of the probability estimate of PD occurrence, obtained from the final Bayesian logistic model, was used for creating a PD risk map. First, the empirical semi-variogram of the point estimates was plotted. Several statistical models (exponential, Gaussian, spherical) were investigated to identify the best-fit variogram model. Those parameters were then used to interpolate values between all the 359 farming sites considering a fixed radius of 40 km and a cell size of 500 m. An ordinary spatial kriging was used and allowed to produce a raster map; the variance of output raster was also mapped. Kriging is a group of geostatistical techniques for interpolating the value of a random field at an unobserved location based on observed values at nearby locations; it has been widely used in veterinary epidemiology [20]. Kriging and mapping were performed with the gstat package of the R software and the spatial analyst extension of ArcGIS software v.9.3 (ESRI Inc. Redlands, CA USA).

Sensitivity analysis

The threshold of 80 km, which was used to define the width of the spatial grid, was selected in order to have an acceptable number of farming sites per hexagon. To see whether the size of the grid would influence the model outputs, we carried out a sensitivity analysis by running the final model with a grid size of 40 km. As this size produced hexagons without any farming sites, the empty adjacent hexagons were merged. This resulted in a final grid of 39 polygons, which were used in the model. DIC and effect estimates of this model were compared to those of the model of originally selected grid size.

Model validation

Relative operation characteristic (ROC) statistical methodology was used to validate the accuracy of the model. The approach has been widely used to access performance of diagnostic tests of animal diseases, and has been applied for validation of models predicting spatial distributions [24-27]. For the present study, the ROC curve presents the plot of the “sensitivity” (true-positive value) versus “specificity” (true-negative value) at each plotted point. The sensitivity is the probability that the model predicts that a site is infected when that site is truly infected. The specificity is the probability that the model predicts that a site is not infected when that site is truly not infected. Each plotted point reflects the average probability that infection will be present on a site. We plotted ROC and calculated an area under the curve (AUC) to check the model fit by comparing the estimated probability of PD occurrence and the observed PD status during the study period. Also, we compared the predicted probabilities based on the final model with the observed PD status in 2011 in order to evaluate if this model could be used for predicting future outbreaks. If the PD occurrence is 100% predicted by the model, AUC of 1 is obtained [24]. All those performance measures were calculated using the pROC package of the R software [28].

Results

Univariate and multivariate analyses

A total of 359 salmon fish farming sites were included in the present study based on the inclusion criteria. Of this study population, 127 sites (35%) were classified as PD cases during the study period. Geographical distribution of the study sites classified by PD status is presented in Figure 1. DIC and coefficient estimates from the univariate logistic regression analysis are presented in Table 1. According to our significance criteria, all covariates except the smolt cohort were further evaluated in multivariate analysis. Table 2 presents the model covariates and DIC of various models depending on combinations of covariates.

Table 1. Univariate analysis using a Bayesian approach of potential risk factors for PD occurrence of 359 Atlantic salmon farming sites located in the southern part of Norway from 2009 to 2010

Table 2. Multivariate analysis using a Bayesian approach of potential risk factors for PD occurrence of 359 Atlantic salmon farming sites located in the southern part of Norway from 2009 to 2010

There was no evidence of multicollinearity among the covariates. Linear correlation between LBD and site density was moderate (Pearson coefficient equals to 0.63 (p < 0.05). A generalized additive model of LBD as a nonlinear function of site density showed a linear positive association for small site densities, the association leveling off for the higher site densities. The final model consisted of site latitude, site density, PD history, and LBD. Although there was an evidence of a spatial cluster of PD infected sites, which was indicated by a decrease of DIC after adding a spatial component to a model with an intercept alone (Table 1), the PD predicted probability was well explained by the covariates of the final model and the spatial component was finally removed. A plot of empirical semivariogram of model residuals and simulation envelopes (data not shown here) confirmed that the spatial pattern of PD cases was well accounted for by the final model. According to the final model, the site latitude alone had a negative effect on the PD probability (median = -0.42 with 95% probability interval between -0.73 and -0.15) meaning that the risk is increasing as one goes south. Having a history of surrounding sites with PD suspected/confirmed the previous year would increase the risk of getting PD (median = 1.22, 95% probability interval between 0.58 and 1.89). Increasing LBD 10 times of an average LBD would increase the risk of getting PD by 1.24 (95% a probability interval between 0.04 and -2.59). At last, increasing site density by adding one more site in a 10-km seaway distance would slightly increase the risk of getting PD (median = 0.08 with a 95% probability interval between -0.03 and 0.18).

Probability estimates of PD occurrence

Based on the results of the final model, probability estimates of PD occurrence of the study sites ranged from 0 to 0.75 with a median of 0.4, and about 34% of the study sites had the probability estimates of ≥ 0.5. An interpolated map of the PD predicted probability based on the final model covering the southern part of Norway is presented in Figure 2. To improve the legibility of the map, we classified the predicted probability into 6 categories ranging from 0 to 1.0. The variance of output raster, which helps assessing the uncertainty of the PD predicted probability estimates produced by kriging, ranged from 1-5%. The predictions are more precise in the center of the circles than at their edges (Appendix 1).

thumbnailFigure 2. An interpolated map presenting PD predicted probability based on the final model of the present study.

Model validation

Figure 3 shows the ROC curves for validating the capacity of the final model for: a) model estimation using observed PD data (359 sites) from 2009 to 2010, b) model prediction using observed PD data of 2011 (438 sites). A reasonable predictive capacity of the final model was observed. The AUC were 0.76 (95%CI: 0.71-0.81), and 0.72 (95%CI: 0.66-0.78) for the model estimation, and the model prediction, respectively.

thumbnailFigure 3. ROC curves for a) validation of the model estimates with observed PD occurrence in time period 2009-2010 (AUC = 0.76, 95%CI: 0.71-0.81), b) validation of the model prediction with observed PD occurrence in 2011 (AUC = 0.72, 95%CI: 0.66-0.78). Each plotted point reflects the average probability that infection will be present on a site.

Sensitivity analysis

For the model with only a spatial component, changing the grid size yielded different estimate of model DICs. The models corresponding to a width of 40 km and 80 km yielded DICs of 424, and 440, respectively. However, changing grid size presented similar DICs and model estimates in the full model (Table 3). This indicates that the effect of spatial congruity on PD occurrence, if occurred, was not influenced by changing of grid size when other significant explanatory variables were included in the model.

Table 3. Sensitivity analysis of the effect of changing grid size on DIC and coefficient estimates of the final model

Discussion

In the present study, we created quantitative PD risk maps of Norwegian Atlantic salmon farming sites accounting for influential factors of PD occurrence at site level. The study was based on the observed data from 1 January 2009 to 31 December 2010. The influential factors included site latitude, site density, PD history, and local biomass density. Subsequently, the probabilities of PD occurrence at site level were displayed in a map covering the southern part of Norwegian marine aquaculture, and validated with observed data of 2011.

The present study included sites located in the southern part of Norway where PD is endemic as well as sites in the southernmost part of the non-endemic area bordering the endemic zone. By including this non-endemic area, the strength of the model to predict the development of PD in this area could be observed. PD strategies in the PD endemic area mainly focus on control and prevention of the disease spreading outside of the affected sites and the endemic area. In the endemic area, with movement restriction, the PD-affected sites have been allowed to keep their fish until slaughtering time. Vaccination has been practiced in the majority of the sites in this PD endemic area; however, the efficacy and effectiveness of the currently used vaccine is debatable [29]. Conversely, vaccination is not a common practice in the PD non-endemic area (the northern part of Norway), but a stamping out policy is regulated and followed by a fallowing period [12]. So far, relatively few cases of PD have been detected in the non-endemic area; therefore, controlling PD outbreak in this area has been manageable financially. Nonetheless, the confining of PD outbreaks within the endemic area may not be feasible in the near future. A great number of financial supports and a strong commitment of fish farming industries would be needed if more stringent control strategies would be applied to prevent the disease spreading northward and to minimize the number of PD outbreaks in the endemic area. An economic model accessing a merit of new disease control strategies would be an assistant tool for decision makers.

The present study showed that the recent history of PD was a useful predictive factor of upcoming PD outbreaks. We found a positive association between the occurrence of PD outbreaks in 2009-2010 and the sites located within 10-km seaway distance of PD suspected/confirmed sites in 2008. This result supports the findings that increasing risk of getting PD is associated with a close distance to infected farming sites [7,14]. We explored the influence of LBD and site density on PD occurrence, and found that these variables also yielded a positive association with increasing risk of getting PD after adjusting for the effects of the other variables. Both variables (LBD, site density) present a concentration of fish; however, they could yield different interpretations. A large site density could have a low LBD if the biomass of neighboring sites is low. The large site density means that there are many sites located in a certain area (in this case 10 km seaway distance). Therefore, frequent visits or multiple activities, such as movement of live fish by a well boat, leading to a higher chance of introducing and spreading the disease agent. On the other hand, a small number of sites each with a large biomass could result in a large overall value of LBD being likely to affect the spreading of virus more than introducing it.

The present study found a non-significant association between the time when smolt were put to sea and the PD occurrence. It has been hypothesized that the cohorts of autumn smolts were more susceptible to PD occurrence than cohorts of spring smolts due to a smaller starting weight of the fish, a lower seawater temperature and a shorter day length, etc [7]. The previous study used a fish cohort as a unit of analysis, and included the cohorts since 2002, which could introduce confounding due to the changes in site management and environment over a period of time. Further works are needed to explore in more details the effect of time when smolt cohorts are put to sea on PD occurrence.

Our results showed a significant effect of latitude, with the risk of getting PD increasing as one goes south. Interpretation of this result is sensitive as latitude can be correlated with several unmeasured variables related to PD risk, including seawater temperature. Exploring the effect of seasonality and seawater temperature on the PD occurrence would be of great interest. However, a cyclical pattern and large missing records of this variable in our study population made its integration into our spatial model challenging. An evaluation of this variable using a subset of the study population was attempted by assessing the proportion of months when the seawater temperature is suitable for survival of SAV [30]. Partly due to a similar pattern of seawater temperature of the subset of the study population we did not find this variable significant (data not shown here). A further evaluation of the effect of this variable on the PD occurrence would require a selection of farming sites with a difference in the pattern of seawater temperature, good-quality and comprehensive data, and would be facilitated by a more dynamic time-scale incorporated into the spatial analysis.

Other research groups have conducted epidemiological studies to determine risk factors of PD transmission [5,31,32]. A close proximity to infected sites, sharing the same water current network, sharing the same equipment with infected sites, sharing of farm ownerships, movement of fish and well boats, inadequate fallowing period, and abundance of sea lice were suggested as potential PD risk factor. This knowledge is invaluable; however, the accessibility and data quality for an extensive area is limited. Another challenging aspect is to clearly distinguish between the risk factors influencing SAV transmission and those inducing PD outbreak.

Peeler and Taylor (2011) published a thorough review emphasizing importance of applying epidemiology tools in a study of aquatic animal diseases, and mentioned a little use of tools describing a spread of infectious disease [33]. In the case of PD, a close proximity of PD-infected sites is clearly a risk factor for the disease transmission [7,14], and sea currents can greatly contribute to spread of the virus between sites [9,13]. Spatial autocorrelation often occurs in outbreaks of infectious diseases in which a farming site located closer to other infected sites is likely to experience an increased chance of acquiring the disease. If this spatial correlation exists and is not accounted for in the model, it could lead to biased parameter estimates and overly optimistic standard errors [20,34]. In our case, despite the spatial clustering of PD infected sites, we did not find a significant improve by adding a spatial component to the final model. This indicates that the predictor variables appropriately accounted for the local spatial correlation.

The use of a Bayesian modeling approach in this work instead of a frequentist approach for estimating the probability of PD occurrence offers a few advantages. It is a flexible tool for accounting for hierarchical levels, such as spatial dependencies [35]. The method allows incorporating previous information providing an appropriate setting for complex models and missing data problem, and the approach yields results in a form of probability distribution that is interpretable intuitively. In our study, non-informative priors were used throughout due to uncertainty of previous information of the parameter estimates. However, the model estimates obtained from this present study could be later used as the prior knowledge to update the future prediction of PD occurrence. The outputs of the Bayesian method gave the posterior distribution of monitoring parameters. We extracted the median of the posterior distribution of the probability of PD occurrence, and displayed in a so-called “PD risk map”. The raster map of PD risk that we produced made it possible to identify local “cold” (green color) and “hot” spots (red color) of PD risk in the area where PD outbreaks area commonly occurring. Based on the PD risk map, the PD predicted probability was high (> 0.5) in the area of Hordaland and some area of Rogaland [Figure 2. In contrast, the probability in Møre and Romsdal was moderate. Without considering the site density and LBD, one could mistakenly believe that there is a higher PD risk in Møre and Romsdal based solely on observed data.

Area under the curve indicated that the raster map produced from data collected in 2009-2010 had a reasonable capacity in discriminating infected farming sites from non-infected ones for the year 2011. This provided a considerable confidence that the risk map could be used as a tool to select areas with an acceptable probability of PD occurrence for salmon farming sites, and to focus surveillance and control measures on high-risk areas. The spatial congruity of infectious aquatic animal diseases is undeveloped and could be explored more in the future by including the use of time-dynamic spatial risk models, and the investigation of the spread of pathogenic agents (with no clinical outbreak) using intensive screening programs. The framework of the study could be applied for spatial studies of other infectious aquatic animal diseases.

Conclusions

This study presents quantitative risk maps of pancreas disease of Atlantic salmon in Norway. It also provides an insight into the relationship between the risk of PD outbreaks in farming sites and several spatial determinants, including latitude, site density, recent history of PD in a nearby area and local biomass density. The framework of the study could be applied for spatial studies of other infectious aquatic animal diseases.

Appendix 1

Variance of PD predicted probability.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

ST coordinated the work, performed all statistical analyses of the present study including ROC/AUC, and drafted the manuscript. MP contributed to development of the Bayesian code, participated in the analyses of the spatial aspect of the study, and assisted in manuscript writing. HV and EB were involved in the study design, analyses of potential risk factors, and assisted in manuscript writing. DA contributed to development of the Bayesian approach. DJ carried out the handling of the raw data, and assisted in performing descriptive statistics of each variable evaluated in the study. All authors read and approved the final manuscript.

Acknowledgements

The study was financially supported by the Research Council of Norway (RCN): TRISAV and SIP RCN project number 190484 and 186907. We thank Anja B. Kristoffersen and Brit Bang Jensen for a critical review of the manuscript. We are also grateful to Matthew Farnsworth and Michael Ward for an intensive training in spatial epidemiology.

References

  1. Taksdal T, Olsen AB, Bjerkås I, Hjortaas MJ, Dannevig BH, Graham DA, McLoughlin MF: Pancreas disease in farmed Atlantic salmon, Salmo salar L., and rainbow trout, Oncorhynchus mykiss (Walbaum) in Norway.

    J Fish Dis 2007, 30:545-558. PubMed Abstract | Publisher Full Text OpenURL

  2. Graham DA, Frost P, McLoughlin K, Rowley HM, Gabestad I, Gordon A, McLoughlin MF: A comparative study of marine salmonid alphavirus subtype 1-6 using an experimental cohabitation challenge model.

    J Fish Dis 2011, 33:123-135. OpenURL

  3. Aunsmo A, Valle PS, Sandberg M, Midtlyng PJ, Bruheim T: Stochastic modeling of direct costs of pancreas disease (PD) in Norwegian farmed Atlantic salmon (Salmo salar L.).

    Prev Vet Med 2010, 93:233-241. PubMed Abstract | Publisher Full Text OpenURL

  4. Poppe T, Rimstad E, Hyllseth B: Pancreas disease of Atlantic salmon Salmo salar L. post-smolts infected with infectious pancreatic necrosis virus (IPNV).

    Bull Euro Asso Fish Pathol 1989, 9:83-85. OpenURL

  5. McLoughlin MF, Peeler EJ, Foyle KL, Rodger H, O’Ceallachain D, Geoghegan F: An epidemiological investigation of the re-emergence of pancreas disease in Irish farmed atlantic salmon (Salmo salar L.) in 2002.

    Marine Environment and Health Series 2003, 14:1-14. OpenURL

  6. McLoughlin MF, Graham DA: Alphavirus infections in salmonids – a review.

    J Fish Dis 2007, 33:511-531. OpenURL

  7. Kristoffersen AB, Viljugrein H, Kongtorp RT, Brun E, Jansen PA: Risk factors for pancreas disease (PD) outbreaks in farmed Atlantic salmon and rainbow trout in Norway during 2003–2007.

    Prev Vet Med 2009, 90:127-136. PubMed Abstract | Publisher Full Text OpenURL

  8. Fringuelli E, Rowley HM, Wilson JC, Hunter R, Rodger H, Graham DA: Phylogenetic analyses and molecular epidemiology of European salmonid alphaviruses (SAV) based on partial E2 and nsP3 gene nucleotide sequences.

    J Fish Dis 2008, 31:811-823. PubMed Abstract | Publisher Full Text OpenURL

  9. Viljugrein H, Staalstrøm A, Molvær J, Urke HA, Jansen PA: Integration of hydrodynamics into a statistical model on the spread of pancreas disease (PD) in salmon farming.

    Dis Aquat Organ 2009, 88:35-44. PubMed Abstract | Publisher Full Text OpenURL

  10. Kongtorp RT, Stene A, Andreassen PA, Aspehaug V, Graham DA, Lyngstad TM, Olsen AB, Olsen RS, Sandberg M, Santi N, Wallace C, Breck O: Lack of evidence for vertical transmission of SAV 3 using gametes of Atlantic salmon, Salmo salar L., exposed by natural and experimental routes.

    J Fish Dis 2010, 33:879-888. PubMed Abstract | Publisher Full Text OpenURL

  11. Jansen MD, Taksdal T, Wasmuth MA, Gjerset B, Brun E, Olsen AB, Breck O, Sandberg M: Salmonid alphavirus (SAV) and pancreas disease (PD) in Atlantic salmon, Salmo salar L., in freshwater and seawater sites in Norway from 2006 to 2008.

    Prev Vet Med 2010, 33:391-402. OpenURL

  12. Anonymous. Norwegian PD regulations 2007-11-20 number 1315: regulations relating to zone to prevent infection and fight pancreas disease in aquaculture. Norwegian: PD regulations 2008-06-17 number 819: regulations relating to the sale of aquaculture animals and products of aquaculture animals; 2008. OpenURL

  13. Stene A, Viljugren H, Jansen P-A: Transmission of SAV3 by water transport? Current patterns and infection pressure between farms. Belfast: PD-trination meeting; 2011. OpenURL

  14. Aldrin M, Storvik B, Frigessi A, Viljugrein H, Jansen PA: A stochastic model for the assessment of the transmission pathways of heart and skeleton muscle inflammation, pancreas disease and infectious salmon anemia in marine fish farms in Norway.

    Prev Vet Med 2010, 93:51-61. PubMed Abstract | Publisher Full Text OpenURL

  15. Jansen PA, Kristoffersen AB, Viljugrein H, Jimenez D, Aldrin M, Stien A: Sea lice as a density dependent constraint to salmonid farming.

    Proc R S B 2012,  .

    published online. 1-9

    Publisher Full Text OpenURL

  16. Lunn DJ, Thomas A, Best N, Spiegelhalter DJ: WinBUGS – a Bayesian modelling framework: concepts, structure, and extensibility.

    Stat Comput 2000, 10:325-337. Publisher Full Text OpenURL

  17. Spiegelhalter DJ, Best N, Carlin BP, Van Der Linde A: Bayesian measures of model complexity and fit.

    J R Stat Soc Series B Stat Methodol 2002, 64:583-639. Publisher Full Text OpenURL

  18. Dormann CF, McPherson JM, Araujo MB, Bivand R, Bolliger J, Carl G, Davies RG, Hirzel A, Jetz W, Kissling WD, Kuhn I, Ohlemuller R, Peres-Neto PR, Reineking B, Schroder B, Schurr FM, Wilson R: Methods to account for spatial autocorrelation in the analysis of species distribution data: a review.

    Ecography 2007, 30:609-628. Publisher Full Text OpenURL

  19. Diggle PJ, Ribeiro PJ: Model-based Geostatistics. Springer, New York, USA: Springer Series in Statistics; 2007. OpenURL

  20. Pfeiffer DU, Robinson TP, Stevenson M, Stevens KB, Rogers DJ, Clements ACA: Spatial Analysis in Epidemiology. New York: Oxford University Press; 2008. OpenURL

  21. Pullan RL, Kabatereine NB, Bukirwa H, Staedke SG, Brooker S: Heterogeneities and Consequences of Plasmodium Species and Hookworm Coinfection: a Population Based Study in Uganda.

    J Infect Dis 2011, 203:406-417. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Besag J, York J, Mollié A: Bayesian image restoration, with two applications in spatial statistics.

    Ann Inst Stat Math 1991, 43:1-59. Publisher Full Text OpenURL

  23. Mollié A: From Bayesian and Empirical Bayes approaches to disease mapping. In Disease mapping and risk assessment for public health. Edited by Lawson A, Biggeri A, Bohning D, Lesaffre E, Viel JF, Bertollini R. Chichester, UK: John Wiley & Sons Ltd; 1999:15-29. OpenURL

  24. Greiner M, Pfeiffer D, Smith RD: Principles and practical application of the receiver-operating characteristic analysis for diagnostic tests.

    Prev Vet Med 2000, 45:23-41. PubMed Abstract | Publisher Full Text OpenURL

  25. Fosgate GT, Adesiyun AA, Herd DW, Johnson WO, Hietala SK, Schurig GG, Ryan J: Estimation of receiver-operating characteristic curves to determine accuracy of a competitive enzyme-linked immunosorbent assay for the serodiagnosis of Brucella infection in domestic water buffalo (Bubalus bubalis) and cattle.

    Am J Vet Res 2003, 64:57-64. PubMed Abstract | Publisher Full Text OpenURL

  26. Pontius RG, Schneider LC: Land-cover change model validation by an ROC method for the Ipswich watershed, Massachusetts, USA.

    Agric Ecosyst Environ 2001, 85:239-248. Publisher Full Text OpenURL

  27. Arekhi S: Modelling spatial pattern of deforestation using GIS and logistic regression: a case study of northern llam forests, llam province, Iran.

    Afr J Biotechnol 2011, 10:16236-16249. OpenURL

  28. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, Müller M: pROC: an open-source package for R and S + to analyze and compare ROC curves.

    BMC Bioinforma 2011, 12:77. BioMed Central Full Text OpenURL

  29. Karlsen M, Tingbø T, Solbakk I, Evensen Ø, Furevik A, Aas-Eng A: Efficacy and safety of an inactivated vaccine against Salmonid alphavirus (family Togaviridae).

    Journal of Vaccines 2012.

    In press

    OpenURL

  30. Graham DA, Staples C, Wilson CJ, Jewhurst H, Cherry K, Gordon A, Rowley HM: Biophysical properties of salmonid alphaviruses: influence of temperature and pH on virus survival.

    J Fish Dis 2007, 30:533-543. PubMed Abstract | Publisher Full Text OpenURL

  31. Wheatly SB, McLaughlin MF, Menzies FD, Goodall EA: Site management factors influencing mortality rates in Atlantic salmon (Salmo salar L.) during marine production.

    Aquaculture 1995, 136:195-207. Publisher Full Text OpenURL

  32. Rodge H, Mitchell S: Epidemiological observations of pancreas disease of farmed Atlantic salmon, Salmo salar L., in Ireland.

    J Fish Dis 2007, 30:157-167. PubMed Abstract | Publisher Full Text OpenURL

  33. Muller EJ, Taylor NG: The application of epidemiology in aquatic animal health-opportunities and challenges.

    Vet Res 2011, 42:94. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  34. Firestone SM, Christley RM, Ward MP, Dhand NK: Adding the spatial dimension to the social network analysis of an epidemic: investigation of the 2007 outbreak of equine influenza in Australia.

    Prev Vet Med 2012.

    In press

    OpenURL

  35. Torabi M, Rosychuk RJ: Hierarchical Bayesian Spatiotemporal Analysis of Childhood Cancer Trends.

    Geogr Anal 2012, 44:109-120. Publisher Full Text OpenURL