Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention16485486336763705-059610.3201/eid1112.050596DispatchDispatchPhocine Distemper Outbreak, the Netherlands, 2002Phocine Distemper Outbreak, the Netherlands, 2002RijksJolianne M.*Van de BildtMarco W.G.JensenTrine§PhilippaJoost D.W.OsterhausAlbert D.M.E.*KuikenThijs*Dutch Wildlife Health Center, Rotterdam, the NetherlandsFaculty of Veterinary Medicine, Utrecht, the NetherlandsErasmus Medical Center, Rotterdam, the NetherlandsRoyal Veterinary and Agricultural University, Frederiksberg, DenmarkAddress for correspondence: Thijs Kuiken, Institute of Virology, Erasmus MC Rotterdam, PO Box 1738, 3000 DR Rotterdam, the Netherlands; fax: +31-10-408-9485; email: t.kuiken@erasmusmc.nl122005111219451948

During the 2002 phocine distemper epidemic, 2,284 seals, primarily harbor seals (Phoca vitulina), were found stranded along the Dutch coast. Stranding pattern varied with age, sex, state of decomposition, wind, and location. Cumulative proportion of deaths (54%) was comparable to that in the first reported epidemic in 1988.

Keywords: EpidemiologyPhocine distemper virusPhoca vitulinathe Netherlandsdispatch

Marine mammal morbilliviruses are among the most pathogenic infectious agents to emerge in wildlife. Phocine distemper virus (PDV) infection (13) was considered responsible for the deaths of ≈18,000 seals in Europe in the first recorded outbreak in 1988 (4), and of ≈22,000 seals in the second outbreak in 2002 (5,6). We examined the effect of different variables on the dynamics of the 2002 PDV epidemic in the Netherlands. This epidemic started 6 weeks after the first cases were noted on Anholt Island, Denmark (5). Subsequently, the disease spread east to Germany and Denmark, and west to Belgium, France, the United Kingdom, and Ireland (6). We also compared the epidemiologic characteristics of the 1988 and 2002 PDV epidemics in the Netherlands.

Seal strandings were reported to a central telephone service. Live stranded seals were rehabilitated or euthanized. Dead stranded seals were collected for necropsy during which species, sex, standard body length, and state of decomposition were determined. Seals were divided into age categories, based on sex and standard body length (7): male juveniles (age <1 yr; length <95 cm), subadults (1 yr < age <4 yr; 95 cm < length <140 cm), or adults (age >4 yr; length >140 cm); female juveniles (age <1 yr; length <90 cm), subadults (1 yr < age <3 yr; 90 cm < length <130 cm), or adults (age >3 yr; length >130 cm) (Table 1).Of 1,315 seals that underwent necropsy, complete data were obtained for 1,096 harbor seals (Phoca vitulina) (Table 1). These seals originated from the entire Dutch coast, except from the islands Rottumeroog and Rottumerplaat, where they were buried; from the island of Texel, where they were collected for a different study; and from the mainland coasts of North Holland and South Holland, where only a few seals were submitted for necropsy (Table 1, Figure 1A). Because seals on which a necropsy was performed represented 56%–73% of the stranded seals in the remaining locations (Table 1, Figure 1A) and had a similar-shaped epidemic curve to that of stranded seals (Figure A1), they were considered representative of stranded seals. The daily wind factor was calculated by multiplying average daily wind force at Den Helder, North Holland (obtained from the Royal Netherlands Meteorological Institute [KNMI]), with its coefficient. Coefficients were positive for winds north of the line west-southwest–east-northeast, negative for winds south of this line, and ranged from 0, when the wind direction was parallel to this line, to 4, when at right angles to it. To analyze the effect of spring tide, the number of strandings on the day of spring tide and the 2 subsequent days was compared to the number of strandings on other days. We used the χ2 test for categorical comparisons and linear trends, and Mann-Whitney U and Kruskal-Wallis tests for temporal scales, with pair-wise comparison for the variables that showed significant overall effect (SPSS for Windows, SPSS Inc., Chicago, IL, USA). For stranded seals with missing observations, age category, sex, and state of decomposition were imputed by using data on seals that underwent necropsy and had been stranded in the same location and on the same or closest weekly date.

Age and sex distribution of harbor seals stranded during 2002 phocine distemper virus epidemic in the Netherlands
Stranding locationCoastline (km)*Harbor seals
No. strandedNo. that underwent necropsy†Juvenile‡
Subadult‡
Adult‡
MFMFMF
Texel57281
Vlieland413031982535861231123
Terschelling673381904236921341717
Ameland49279156222287963418
Schiermonnikoog293312193331821015232
Rottumeroog and Rottumerplaat14172
Friesland mainland9321715982157662738
Groningen mainland80194138183428502440
North Holland12497
South Holland3618
Zealand1165136
Total7062,2841,096148179432570165168

*Source: ref (8).

†These are harbor seals for which complete data were obtained. Seals from Texel, Rottumeroog and Rottumerplaat, North Holland, and South Holland were excluded because of the small proportions that underwent necropsy.

‡Age category and sex for stranded seals with missing observations were imputed by using those of seals that were stranded in the same location and on the same or closest weekly date.

Spatial and temporal distribution of seal strandings in the Netherlands during the 2002 phocine distemper virus epidemic. A) Spatial distribution of seal strandings and proportion of seals necropsied at each location. The diameter of each pie chart corresponds to the number of seals stranded at a particular location. The names of the Wadden Sea islands have been abbreviated (Tx, Texel; V, Vlieland; Ts, Terschelling; A, Ameland; S, Schiermonnikoog; R, Rottumeroog and Rottumerplaat). B) Weekly stranding rate of all stranded harbor seals and effect of age category on weekly stranding rate.

Between June 16, 2002, when the Dutch index case was found on Vlieland, and the end of November 2002, when the stranding rate returned to preepidemic levels, 2,284 seals (2,154 dead, 130 live) were stranded along the Dutch coast (Figure 1B, Table 1). Almost all (2,279 of 2,284) were identified as harbor seals, and the remaining 5 as gray seals (Halichoerus grypus), despite recently increased gray seal numbers in the Netherlands and their likely exposure to PDV. This finding is consistent with experimental findings that PDV infection is more pathogenic for harbor seals than for gray seals (9). At gross necropsy, ≈80% of harbor seals had pulmonary consolidation consistent with PDV infection, while about 50% had either immunoglobulin M to morbillivirus by serology, morbillivirus-specific nucleic acid by reverse-transcription polymerase chain reaction (PCR), or both (unpub. data). This PCR fragment corresponded to that of PDV by phylogenetic analysis (5). Together these results confirm PDV infection as the primary cause of the epidemic. The rapid course of the epidemic, high cumulative proportion of deaths, and involvement of all age categories (Table 1) fit with a virgin soil epidemic and correspond with lack of preexisting specific immunity to PDV in most of the seals (5,6).

Age and sex affected temporal distribution of strandings. The median stranding date varied significantly among age categories (p<0.001). The date was significantly earlier for subadults than for juveniles and adults (p<0.05; Figure 1B). Subadults display considerably more social play than seals in other age categories, especially in early summer (10). In contrast, juveniles and their mothers are relatively more separated from other seals during the lactation period, are more sedentary, and have fewer new contacts (11,12). The median stranding date for males was significantly earlier than that for females in juveniles (p<0.001), subadults (p<0.001), and adults (p<0.001) (Figure A1 and A2). The average change in the number of individually identified seals hauled out between consecutive days is significantly higher for males than for females (12), and both subadult and adult males have the longest and most aggressive interactions with each other (13). These behavioral differences suggest that contact rates and intensity of contact with other seals, including seals with PDV infection, were higher for subadults at the start of the epidemic than for juveniles and adults, and higher for males than females, thus increasing the risk and severity of infection. Alternatively, the above patterns may be linked to age-related and sex-related differences in the effects of contaminants. The contaminant levels in the tissues of seals that died in the 1988 PDV epidemic were considered sufficiently high to cause immunosuppression and thus to increase the severity of the PDV outbreak (14). Pollutant levels in tissues of seals that died during the 2002 PDV epidemic have yet to be reported.

Age affected geographic distribution of strandings. The proportion of stranded seals of each age category varied significantly among Wadden Sea locations (p<0.001). The highest proportions of juveniles and adults stranded at mainland Groningen (Table 1), which includes Eemsmond, a core breeding area. The highest proportion of subadults stranded on Vlieland (Table 1) in the western part of the Dutch Wadden Sea, an area assumed to have an influx of migrating young seals (15). The number of seals stranded per kilometer of coastline varied significantly among locations for juveniles (p<0.001), subadults (p<0.001), and adults (p<0.001), with 2.2 to 3.1 more seals stranded per kilometer of coastline on Schiermonnikooog, an island in the eastern part of the Dutch Wadden Sea, than would be expected had the seals been evenly distributed per km coast (Table 1). This coincides with the summer distribution of harbor seals in the Dutch Wadden Sea, which is highly skewed toward the east (15,16). Within each age category, the proportion of males to females varied significantly among locations only for adults (p<0.001). Ameland had the highest proportion of adult males, and Vlieland the lowest (Table 1).

Location affected the temporal distribution of strandings. The median stranding date varied significantly among locations (p<0.001); that for Zealand (week 39) was significantly later than that for all Wadden Sea locations (weeks 35–37) (Figure A2). This is likely because seals in Zealand are fewer and more widely dispersed than in the Wadden Sea, so the chance of the virus spreading is lower.

Wind appeared to have a confounding effect on stranding rate: periods of southerly wind corresponded with decreased overall stranding rates, e.g., in weeks 33 and 36, and the opposite for northerly winds (Figure 2A). This is probably because dead seals floated in the top water layer, which shows parallel drift to surface winds. A similar effect of wind on strandings has been shown for seabirds (17). Spring tide did not affect stranding rate (p>0.05).

Effects of environmental variables on seal strandings in the Netherlands during the 2002 phocine distemper virus epidemic. A) Effect of wind direction and force on temporal distribution of stranded seals. Stranding rate of seals is expressed as number of seals reported per day. The wind factor is a function of wind force and wind direction. Negative wind factors correspond to southerly winds. B) Effect of state of decomposition on temporal distribution of stranded harbor seals, overall and per location. Percentages of decomposed seals are expressed per month.

State of decomposition (as a measure of length of time between death of a seal and its detection) also had a confounding effect on stranding rate. From July to October, the overall proportion of decomposed seals differed significantly among months (p<0.001) and increased significantly with time (p<0.001) (Figure 2B). This finding is probably because recovery of seal carcasses was not 100% so that, as the epidemic progressed, a higher proportion of stranded carcasses consisted of seals that had died before the previous shore survey. These findings show that over time stranding rate became a less accurate estimate of mortality rate, as observed in 1988 (18). The proportion of decomposed carcasses varied significantly by location (p<0.001) (Figure A3), with high proportions of decomposed carcasses on the mainland coasts of Friesland and Groningen and on Schiermonnikoog (Figure A4).

The timing of deaths and cumulative proportion of deaths of the 2002 PDV epidemic were similar to those characteristics in 1988 (Table 2). A difference, however, was that the index case was detected ≈1 month later in 2002 than in 1988. The similarity between estimated cumulative proportion of deaths in 1988 (53% of the population) and 2002 (54%) suggests that the pathogenicity of PDV for the harbor seal population has not changed noticeably. However, more detailed examination of the genetic composition of both the virus and the harbor seal is needed to exclude changes in the host-pathogen relationship.

Comparison of overall characteristics of the 1988 and 2002 phocine distemper virus epidemics in the Netherlands
Variable19882002
Date index caseMay 22*Jun 16
Date median caseSep 4*Sep 2
Central epoch (d)115*93
No. found stranded417*2,284
No. counted in preepidemic year966†3,595‡
No. counted in postepidemic year535†2,365‡
Average annual population growth in preepidemic years (%)19¶
Estimated cumulative proportion of deaths (%)#5354

*Source: ref (4).

†Source: Report "Compilatie van gegevens over zeehonden en zeehondenopvang in de Nederlandse Waddenzee" (available at http://www.waddenzee.nl/fileadmin/content/Dossiers/Natuur_en_Landschap/pdf/zeehondenplatform_compilatierapport.pdf).

‡Source: Wadden Sea Newsletters 2001-3 and 2003-2 (available at http://www.waddensea-secretariat.org).

§Source: ref (19).

¶Source: ref (15).

#Calculated as follows: [(number of live seals counted in preepidemic year + average annual growth in preepidemic years) – (number of live seals counted in postepidemic year – average annual growth in preepidemic years)]/(number of live seals counted in preepidemic year + average annual population growth in preepidemic years).

Suggested citation for this article: Rijks JM, Van de Bildt MWG, Jensen T, Philippa JDW, Osterhaus ADME, Kuiken T. Phocine distemper outbreak, the Netherlands, 2002. Emerg Infect Dis [serial on the Internet]. 2005 Dec [date cited]. http://dx.doi.org/10.3201/eid1112.050596

Acknowledgments

We thank the volunteers who searched the Dutch coast daily for stranded seals; the staff at the Seal Research and Rehabilitation Center in Pieterburen for logistical support; V.O.P. Containers for providing the location to perform necropsies of the seals; the staff of the Dutch Ministry of Agriculture, Nature and Food Quality (LNV-Noord) for providing access to and help with the centralized seal registration data; the Common Wadden Sea Secretariat and the Trilateral Seal Expert Group for international coordination of the outbreak; Fiona Read for her help with geographic coordinates; Xavier Harduin for help in entering and analyzing data; Robin Huisman for help with the Geographic Information Systems program; and Hans Heesterbeek, Rik de Swart, and Cock van Duijn for reviewing the manuscript at different stages of development.

This study received financial support from the Dutch Ministry of Agriculture, Nature and Food Quality.

Weekly stranding rate of all recorded stranded seals compared to that of seals that underwent necropsy.

Weekly stranding rate of stranded harbor seals in each age category, by sex.

Effect of location on temporal distribution of stranded harbor seals. Strandings per location are expressed as a relative cumulative frequency curve. The 50% value of each curve corresponds with the median stranding date for a particular location. Note that strandings at Zealand occur ≈1 month later than at Wadden Sea locations.

Effect of state of decomposition of stranded harbor seals. The diameter of each pie chart corresponds to the number of seals stranded per km of coastline at a particular location.

Ms Rijks is a veterinarian at the Erasmus Medical Center in Rotterdam, the Netherlands. Her research interests include the effect of infectious diseases on wild animal populations and the relationship between wildlife health and development aid.

ReferencesOsterhaus ADME, Vedder EJ Identification of virus causing recent seal deaths. Nature. 1988;335:20 10.1038/335020a03412456Mahy BWJ, Barrett T, Evans S, Anderson EC, Bostock CJ Characterization of a seal morbillivirus. Nature. 1988;336:115 10.1038/336115a03185731Cosby SL, McQuaid S, Duffy N, Lyons C, Rima BK, Allan GM, Characterization of a seal morbillivirus. Nature. 1988;336:1156 10.1038/336115b03185731Dietz R, Heide-Jørgensen M-P, Härkönen T Mass deaths of harbor seals (Phoca vitulina) in Europe. Ambio. 1989;18:25864Jensen T, van de Bildt M, Dietz HH, Andersen TH, Hammer AS, Kuiken T, Another phocine distemper outbreak in Europe. Science. 2002;297:209 10.1126/science.107534312114617Müller G, Wohlsein P, Beineke A, Haas L, Greiser-Wilke I, Siebert U, Phocine distemper in German seals, 2002. Emerg Infect Dis. 2004;10:723515200869McLaren IA Growth in pinnipeds. Biol Rev Camb Philos Soc. 1993;68:179 10.1111/j.1469-185X.1993.tb00731.x8457634Camphuysen CJ Oil pollution and oiled seabirds in the Netherlands, 1969–97: signals of a cleaner sea. Sula. 1997;11:41156Harder T, Willhaus T, Frey H-R, Liess B Morbillivirus infections of seals during the 1988 epidemic in the Bay of Heligoland: III. transmission studies of cell culture-propagated phocine distemper virus in harbour seals (Phoca vitulina) and a grey seal (Halichoerus grypus): clinical, virological and serological results. J Vet Med B. 1990;37:64150 10.1111/j.1439-0450.1990.tb01109.x2267888Renouf D, Lawson JW Quantitative aspects of harbour seal (Phoca vitulina) play. J Zool (Lond). 1987;212:26773 10.1111/j.1469-7998.1987.tb05989.xTraut IM, Ries EH, Donat B, Vareschi E Spacing among harbour seals (Phoca vitulina vitulina) on haul-out sites in the Wadden Sea of Niedersachsen. Z Saugetierkd. 1999;64:513Godsell J Herd formation and haul-out behaviour in harbour seals (Phoca vitulina). J Zool (Lond). 1988;215:8398 10.1111/j.1469-7998.1988.tb04886.xSullivan RM Agonistic behavior and dominance relationships in the harbor seal, Phoca vitulina. J Mammal. 1982;63:55469 10.2307/1380260Ross PS, De Swart RL, Reijnders PJH, Van Loveren H, Vos JG, Osterhaus ADME Contaminant-related suppression of delayed-type hypersensitivity and antibody responses in harbor seals fed herring from the Baltic Sea. Environ Health Perspect. 1995;103:1627 10.1289/ehp.951031627737064Ries EH, Traut IM, Brinkman AG, Reijnders PJH Net dispersal of harbour seals within the Wadden Sea before and after the 1988 epizootic. J Sea Res. 1999;41:23344 10.1016/S1385-1101(98)00052-5Leopold MF, van der Werf B, Ries EH, Reijnders PJH The importance of the North Sea for winter dispersal of harbour seals Phoca vitulina from the Wadden Sea. Biol Conserv. 1997;81:97102 10.1016/S0006-3207(96)00092-4Camphuysen CJ, Heubeck M Marine oil pollution and beached bird surveys: the development of a sensitive monitoring instrument. Environ Pollut. 2001;112:44361 10.1016/S0269-7491(00)00138-X11291451Hall AJ, Pomeroy PP, Harwood J The descriptive epizootiology of phocine distemper in the UK during 1988/89. Sci Total Environ. 1992;115:3144 10.1016/0048-9697(92)90030-V1594933Reijnders PJH, Ries EH, Tougaard S, Nørgaard N, Heidemann G, Schwarz J, Population development of harbour seals Phoca vitulina in the Wadden Sea after the 1988 virus epizootic. J Sea Res. 1997;38:1618 10.1016/S1385-1101(97)00031-2