Injecting epidemiology into population viability analysis: Avian cholera transmission dynamics at an arctic seabird colony

Infectious diseases have the potential to spread rapidly and cause high mortality within populations of immunologically naïve hosts. The recent appearance of avian cholera, a highly virulent disease of birds caused by the bacterium Pasteurella multocida, at remote Arctic seabird colonies is an emerging conservation concern. Determining disease risk to population viability requires a quantitative understanding of transmission potential and the factors that regulate epidemic persistence. Estimates of the basic (R0) and real-time (Rt) reproductive number are critical in this regard - enumerating the number of secondary infections caused by each primary infection in a newly invaded host population and the decline in transmission rate as susceptible individuals are removed via mortality or immunized recovery. Here, we use data collected at a closely monitored common eider (Somateria mollissima) breeding colony located in the Canadian Arctic to examine transmission and host population dynamics. Specifically, we infer epidemic curves from daily mortality observations and use a likelihood-based procedure to estimate changes in the reproductive number over a series of annual outbreaks. These data are interpreted in relation to concurrent changes in host numbers to assess local extinction risk. Consistent with expectations for a novel pathogen invasion, case incidence increased exponentially during the initial wave of exposure (R0 = 2·5; generation time = 6·5 days ± 1·1 SD). Disease conditions gradually abated, but only after several years of smouldering infection (Rt ≈ 1). In total, 6194 eider deaths were recorded during outbreaks spanning eight consecutive breeding seasons. Breeding pair abundance declined by 56% from the pre-outbreak peak; however, a robust population of >4000 pairs remained intact upon epidemic fade-out. Overall, outbreak patterns were consistent with herd immunity acting as a mitigating factor governing in the extent and duration of mortality. Disease mortality is frequently modelled as a form of stochastic catastrophe in wildlife population assessments, whereas our approach gives shape to the functional response between transmission and host population dynamics. We conclude that increased emphasis on