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. So far, 6 subtypes of SAV have been
reported worldwide, whereas, only subtypes 2 and 3 have been detected in Norwegian
salmonid farming sites, Hilde Sindre per communications]. The disease has been a
significant problem in salmonid farming industries in Scotland, Ireland, and Norway. 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.
The characteristics of SAV and PD outbreaks have been thoroughly reviewed. 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 Mre 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. 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. 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.
Methods
Data
In this study, we included Atlantic salmon farming sites located in 6 Norwegian counties (Vest-Agder, Rogaland, Hordaland, Sogn & Fjordane, Mre & Romsdal, SrTrndelag), 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 georeferenced in the Directorate of Fisheries (DEF; 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 6339'5,481"N (the farthest north) to 5812'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. 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). 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. 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. 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. To do so, we applied a Monte Carlo method 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. 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 = 119 (the hexagon in which the site is
localised), and j = 1359 (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:
B 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). 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. 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.
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. 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. 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. All those performance measures were calculated using the pROC package of the R software.
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.
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
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).
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.
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)
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 mode.
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. 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. 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. 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. 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. 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. 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. In the case of PD, a close proximity of PDinfected
sites is clearly a risk factor for the disease transmission, and sea currents can
greatly contribute to spread of the virus between sites. 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. 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. 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 Mre
and Romsdal was moderate. Without considering the site density and LBD, one could
mistakenly believe that there is a higher PD risk in Mre 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.
Further ReadingYou can view the full report by clicking here. |
December 2012