Skip to main content
  • Research article
  • Open access
  • Published:

Seasonality and the effects of weather on Campylobacter infections

Abstract

Background

Campylobacteriosis is a major public health concern. The weather factors that influence spatial and seasonal distributions are not fully understood.

Methods

To investigate the impacts of temperature and rainfall on Campylobacter infections in England and Wales, cases of Campylobacter were linked to local temperature and rainfall at laboratory postcodes in the 30 days before the specimen date. Methods for investigation included a comparative conditional incidence, wavelet, clustering, and time series analyses.

Results

The increase of Campylobacter infections in the late spring was significantly linked to temperature two weeks before, with an increase in conditional incidence of 0.175 cases per 100,000 per week for weeks 17 to 24; the relationship to temperature was not linear. Generalized structural time series model revealed that changes in temperature accounted for 33.3% of the expected cases of Campylobacteriosis, with an indication of the direction and relevant temperature range. Wavelet analysis showed a strong annual cycle with additional harmonics at four and six months. Cluster analysis showed three clusters of seasonality with geographic similarities representing metropolitan, rural, and other areas.

Conclusions

The association of Campylobacteriosis with temperature is likely to be indirect. High-resolution spatial temporal linkage of weather parameters and cases is important in improving weather associations with infectious diseases. The primary driver of Campylobacter incidence remains to be determined; other avenues, such as insect contamination of chicken flocks through poor biosecurity should be explored.

Peer Review reports

Background

The distinctive seasonal pattern of Campylobacter incidence has suggested that seasonal changes in the environment contribute to this pattern. The examination of the numbers of Campylobacter in the small intestine and caeca of chicken over a 12-month period showed significant seasonal variation in the carriage rate which was associated with sunshine and temperature [1]. The analysis of surveillance data from Europe, Canada, Australia and New Zealand found a lack of clear and consistent association with rainfall, humidity, sunshine, and temperature [2]. The analysis of Campylobacter data from Denmark, Sweden, Austria, Finland, New Zealand, Scotland, and Wales showed consistent seasonal patterns of infection with more cases in the summer, but with cases peaking at different times of the year [3]. Campylobacter rates in England and Wales were associated with temperature, sunshine, and rainfall, with temperature being the most significant variable and a 5% increase in Campylobacter cases was associated with a 1 °C rise in temperature to a maximum temperature of 14 °C [4,5,6].

An examination of the relationship between Campylobacter cases (along with Salmonella and pathogenic Escherichia coli) in Canada found an increase in cases of 2.2% for every degree of increase in temperature, with a low temperature threshold of 10 °C; another Canadian study modelled Campylobacter and temperature and found an association with temperature above a threshold of 10 °C [7, 8]. Campylobacter incidence in Massachusetts closely followed the peak in annual temperature with a lag of 2–14 days, and in Georgia (US) a strong association was found between temperature and Campylobacter infections [9, 10].

The impact of weather parameters on Campylobacter differed between temperate and sub-tropical areas of Australia [11]. A spatial examination of Campylobacter in New Zealand found no association with climate variables. Temperature did not appear to drive the seasonality, and there was little association with rainfall [12,13,14]. Jore et al. [15] showed a relationship between temperature with both human infections and chicken colonisation in six European countries. There has also been work in Switzerland suggesting the seasonality of Campylobacter is partly explained by contamination of chicken broilers [16].

While an association with weather in general, and temperature in particular, appears to be a relatively common finding in these studies, nevertheless the link is not as clear, consistent and direct as with Salmonella; nor does it have a well-defined mechanism. Nichols [17] and Ekdahl et al. [18] suggested that flies might explain the epidemiology of Campylobacter, through the direct transfer of contaminated material from raw meats or faeces to ready-to-eat foods, and/or by the seasonal contamination of chicken flocks [19,20,21]. Nichols [17] proposed a fly transmission model that explained the seasonal distribution by the impact of temperature on non-biting fly numbers, as measured by the development time of Musca domestica larvae. Evidence suggested that the use of fly screens on chicken production houses could reduce the seasonal contamination of chickens and subsequent human disease significantly [21]. Similar responses to temperature in both chicken broiler carcase contaminations with Campylobacter and human infections with Campylobacter cases in six different European Countries, implying that the seasonal relationship between human infection and temperature is more likely a result of exposure to chicken [15].

One of the additional complicating factors in the Campylobacter transmission in most developed countries is the significant percentage of cases related to foreign travel, and the seasonality of this additional factor [22, 23]. Hartnack et al. [24] highlighted the necessity of quantifying the transmission of Campylobacter from broiler to humans and to include climatic factors in order to gain further insight in its epidemiology. Sterk et al. [25] used a runoff simulator to analyse the effect of climate change on Campylobacter (along with Cryptosporidium) runoff and human infection risks in the Netherlands; and found that climate change has little overall impact on runoff of Campylobacter (along with Cryptosporidium) from land to the surface waters.

Within climate change research, there are three reasons for undertaking this work. Firstly to demonstrate evidence of the potential short-term consequences on human health from climate change; secondly to understand what factors influence the occurrence of infectious diseases; and thirdly to provide models which can forecast future disease risks. The incidence of Campylobacter in the UK demonstrated both a strong and regular seasonality, and a variation in the number of cases per year [4, 6, 26]. This raises a question of whether the seasonal patterns of Campylobacter have changed over time as a result of climate change and other environmental factors.

Methods

Campylobacter surveillance data were extracted from the Second Generation Surveillance System (SGSS) database of Public Health England from 1989 to 2014. A case was defined as a human faecal sample submitted from a patient that was positive for Campylobacter. All Campylobacter surveillance data were laboratory-confirmed cases and were predominantly based on the isolation of Campylobacter from faeces using selective culture media in a micro aerobic environment; membrane filter, PCR or other approaches made up a very small percentage of cases. A few sentinel laboratory sites conducted screening using a range of PCR targets during the 2012 Olympic Games. Standard antimicrobial testing methods, used in primary diagnostic laboratories, were predominantly disc diffusion methods. Where speciation was reported, it was conducted in diagnostic laboratories using conventional phenotypic methods [6]. Most cases were thought to be symptomatic and included patients with extra intestinal infections [6]. The symptoms were those recorded on the physician request form and generally are in the simplistic form of “D & V” or something equivalent.

The meteorological data were supplied by the Met Office and held on the Medical and Environmental Data a Mash-up Infrastructure (MEDMI) platform (see https://www.data-mashup.org.uk).

Reported cases between 2005 and 2009 were linked to local weather parameters at the laboratory postcode locations; namely: maximum, minimum and average daily temperature and daily rainfall in 30 days before the specimen date using a previously published approach that prevented patient identification [6]. The linkage was conducted on laboratory postcode because this was more complete than patient residential postcode and preserved patient confidentiality. An additional study [27] validated that the weather parameters at the laboratory postcode correlate well with the weather at the domestic postcode of cases.

Comparative conditional incidence

GIS analysis was performed on cases with Lower Layer Super Output Areas (LSOA) using ArcGIS 10.2 for Desktop by ESRI. A map of Campylobacter incidence between 2005 and 2009 from patient’s residential postcode data was produced. Some areas had poor residential postcode completeness resulting in areas with an apparent lower incidence. Comparative conditional incidence (CCI) is a new approach which uses cases and LSOA populations to examine differences between disease incidence at different temperature and rainfall in the 30 days before the specimen date within equal sized datasets. CCI is a sum of cases per week within all LSOAs with same temperature divided by sum of population in the same LSOAs. It is conditional as the weather variables were available for the weeks and LSOAs where there were positive cases and excluded patients without a residential postcode as it was not possible to establish the LSOA.

Wavelet analysis

To detect potential time-varying seasonality in Campylobacter that might result from Salmonella control measures in chicken production, such as improved biosecurity practices, we used wavelet analysis. Campylobacter daily counts between 1989 and 2009 were adjusted from the day of the year using a seven-day moving average and systematic adjustments for the reduced reporting over bank holidays and for long-term trends as previously reported [6]. The analyses were conducted in R using ‘waveletComp’ [28].

Generalized structural time series model

Campylobacter cases across ten geographical regions in England and Wales were examined using the generalized structural time series model (GEST) [29] with a negative binomial distribution, where the natural logarithm of the expected cases of Campylobacter was decomposed into baseline, trend and seasonal components. To investigate the effect of weather variables on the development of Campylobacter cases, explanatory variables of the average temperature (°C) and the total rainfall (mm), of one, two, three and four weeks before diagnosis were included in the model. The temperature variable was examined using fixed and varying-coefficient analyses in order to account for variations in the relationship through the year.

Cluster analysis

We examined the geographic similarity of Campylobacter seasonality in England and Wales in thirty-three sub-regional areas between 1989 and 2009. The fitted seasonality of Campylobacter infections across thirty-three Strategic Health Authorities were plotted in a dendrogram using hierarchical clustering analysis with Ward’s minimum variance. Clusters were drawn up that exhibited a distance of one or more from other clusters; and this was used to examine the seasonality of groups of SHAs. Additional file 1: Section S4 illustrates an example of the GEST analysis of Campylobacter cases in Avon, Gloucestershire and Wiltshire SHA. The clustering was conducted using the ‘hclust’ routine of the stats package in R statistical software. Table 1 illustrates methods used, time period, geography, and reason for choice in the analysis.

Table 1 Methods, time, geography, and reason for choice in the analyses

Results

Incidence was determined from patient residence postcode data (2005–2014) and plotted by LSOA across England and Wales. The geographic representation of cases per 100,000 per week showed areas of the country with an apparently low incidence (Fig. 1a; shown in blue); this is presumed to be where residential postcode data were missing and/or there were under-reporting. Data on individual years of cases showed a gradual improvement in the completeness of reporting from 2005 to 2014 (Fig. 1b).

Fig. 1
figure 1

Incidence. a Average Campylobacter cases per 100,000 population for the 10 years 2005 to 2014 by Lower Layer Super Output Area (LSOA) in England and Wales; b Campylobacter cases per 100,000 per year from 2005 to 2014; c Comparative conditional incidence (CCI) of Campylobacter distributed evenly over ten data points and sorted by maximum, average and minimum temperatures two weeks before. d CCI of Campylobacter on four separate periods plotted against average temperature two weeks before; e CCI of Campylobacter on four-week periods splitting each dataset into two equal parts of temperature two weeks before. f CCI Campylobacter on four-week periods splitting each dataset into two equal parts of average rainfall four weeks before; g R2 values using different numbers of data points based on the average temperature two weeks before

The CCI was estimated using datasets that split cases into ten equal parts by temperature and rainfall. The highest CCI values were found for the minimum temperature (9 °C) two weeks before the specimen date, average temperature (16 °C) two weeks before the specimen date and maximum temperature (21 °C) two weeks before the specimen date (Fig. 1c). The cases at the height of summer when temperatures were highest had lower cases than a few weeks earlier. Between average temperatures of 5 °C and 15 °C, there was a 0.09 case increase in the incidence per 100,000 per week (Fig. 1d). It is unclear how this relates to the actual incidence as the LSOA per weeks that did not have a Campylobacter case which could not be included. Using four week periods, average temperatures two weeks before and splitting the data for each time period into two equal halves by temperature, demonstrated that the relationship between temperature and incidence in each time period differed, forming a cyclic loop (Fig. 1e). For weeks 17 to 20 and 21 to 24 there was an increase in the CCI of 0.175 cases per 100,000 per week for each degree increase in the temperature two weeks before (Fig. 1e). A similar approach using average rainfall over the previous four weeks showed no consistent links with the CCI (Fig. 1f). Adjusting the number of data points to give larger or smaller percentages of the data in each point affected the R2 value of the CCI on all data, with optimum results for 200 data points or less (Fig. 1g).

The number of Campylobacter cases per week from 2005 to 2009 in England and Wales, average rainfall (mm) in the previous four weeks, maximum, minimum and average temperature (°C) in the two weeks before the specimen date are shown in Fig. 2a. Numbers of cases per week and average temperature two weeks previously associated with cases on colour are shown in Fig. 2b.

Fig. 2
figure 2

Campylobacter cases and weather. a Campylobacter cases per week from 2005 to 2009 in England and Wales; average rainfall per day in the previous four weeks; maximum, minimum and average temperature two weeks before the specimen date; b Cases per week and average temperature two weeks before the specimen date associated with cases

Wavelet analysis of Campylobacter cases showed a strong annual periodicity with smaller semi-annual and 4 monthly periodicities, i.e. the ‘harmonics’ corresponding to the red horizontal regions in the power spectrum (Fig. 3a) and the three main peaks in the global power spectrum (Fig. 3d) located at 365, 180 and 121 days. Furthermore, wavelet analysis for temperature showed a strong annual seasonality but was also dominated by semi-annual. However, there were years (e.g. from April 2004 to July 2008) where this component was absent in the temperature, but present in the Campylobacter cases. Rainfall data over the same time period exhibited different patterns, except for the annual seasonality.

Fig. 3
figure 3

Wavelet analysis for Campylobacter cases reported between 1989 and 2009 and for temperature and rainfall in England and Wales; a Wavelet power spectrum of the root transformed time-series of daily Campylobacter cases adjusted using a seven-day rolling mean, removal of bank holiday artefacts and adjusted for long term trend; b Wavelet power spectrum of the root transformed time-series of daily temperature; c of daily Rainfall. Low values of the power spectrum are shown in dark blue and high values in dark red. The black lines show the maxima of the undulations of the wavelet power spectrum. The light white shaded areas identify the region influenced by edge; d Global average wavelet power spectrum of the root transformed time-series of Campylobacter cases, the black dots show the 5% significant levels computed based on 500 bootstrapped series; e as d but for averaged temperature

GEST with single coefficients for temperature and rainfall found a statistically significant association between Campylobacter and temperature of two weeks before the case diagnosis (see Additional file 1: Table S1). The relationship between Campylobacter and temperature was positive, as for every degree increase in average weekly temperature (°C) of two weeks before the specimen date there is an increase of 0.73% in the mean number of Campylobacter cases. For rainfall, there was a slightly significant association between Campylobacter and total rainfall of one week before the case diagnosis. The relationship was negative, as for every millimetre increase in total rainfall of one week before the specimen date there is a decrease of 0.24% in the mean number of Campylobacter cases.

Additional file 1: Table S2 provides clear evidence for an association between Campylobacter and average, maximum and minimum temperatures of two weeks before the case diagnosis in most geographical regions of England. Wales did not show significant association. All the regions showed statistically significant increasing trend in cases except the East of England. Additional file 1: Figure S1 shows the weekly cases of Campylobacter in England and Wales with the fitted mean in red (panel a), and the decomposition of the fitted mean into fitted trend (panel e), fitted seasonality (panel f), fitted linear effect of average temperature (panel c) and fitted linear effect of total rainfall (panel d). Expected cases of Campylobacter increased with the average temperature two weeks before diagnosis above a threshold of 11 °C. The relative change in the expected cases of Campylobacter due to temperature and rainfall in England and Wales were 7.7% and − 1.9% respectively and the effect of seasonality was − 3.4% (see Additional file 1: Table S3).

In GEST with varying coefficients of temperature, using 13 four week periods as factors through the year and a single coefficient for rainfall, the relative change in the expected cases of Campylobacter due to temperature in England and Wales increased to 33.3%, the effect of rainfall was − 1.7%, and the effect of seasonality was − 0.5% (see Additional file 1: Table S4). Figure 4 shows the weekly cases of Campylobacter in England and Wales and the fitted mean with varying coefficients of temperature, rainfall, seasonality and trend, and the decomposition of the fitted mean. The fitted effect of 13 four week periods on temperature is shown in Additional file 1: Figure S2c.

Fig. 4
figure 4

GEST analysis for Campylobacter cases in England and Wales from 2005 to 2009 using temperature and rainfall as explanatory variables, and seasonal hierarchical clustering. a Cases with the fitted mean decomposed by trend, seasonality, temperature and rainfall; b Fitted linear trend in time; c Fitted seasonality; d Fitted varying coefficients for temperature; e Fitted fixed coefficient for rainfall; f Campylobacter cases; g Sub-regional clustering by seasonality using Ward’s Minimum Variance based on the GEST model. h Seasonal distribution of Campylobacter cases in the three main clusters

The Ward’s minimum variance clustering of the fitted seasonality of the GEST showed three main clusters: metropolitan (London and Birmingham), rural (Wales and the North East), and other geographical areas (Fig. 4g). When cases were plotted by week of year for the three clusters and the same time period, the seasonal distribution differed significantly (Fig. 4h).

Discussion

The epidemiology of Campylobacter in humans remains a subject of importance because it is the commonest bacterial cause of diarrhoea in most developed countries, and because intervention requires an understanding of disease transmission pathways and possible association with climate and other environmental change. Within the context of climate change research, the approaches to linking weather parameters with infectious disease data commonly rely on an assumption that the relationships are continuous throughout the year, sometimes assuming a threshold value above which the continuous relationship changes. Here, we used different methods to demonstrate that Campylobacter incidence changes dynamically through the season, and we were able to decouple the relationship of temperature from the seasonality. Most notable was the very steep increase in the number of cases diagnosed between weeks 19 and 24, although there was still a relationship to temperature before this (weeks 1-18) as well as later in the year (weeks 29-53). The distribution was also consistent with the direct and/or indirect fly transmission hypotheses [17,18,19], with the fly populations increasing during the late spring period resulting in transmission passing from faeces or other contaminated material to either food and/or chicken flocks (as a result of poor biosecurity). The expression of weekly incidence by the temperature during that week in an LSOA and averaged for all weeks with one or more Campylobacter cases, gives a means of examining the relationship between case numbers and temperature that is partially decoupled from seasonality and geography.

The GEST provides a useful approach for modelling over-dispersed counts with a negative binomial distribution. It allows for the extraction of hidden trend and seasonality and estimating the effect of weather parameters. The use of multilevel factors for temperature separated into 13 four week periods allowed the contribution of temperature two weeks prior to the diagnosis to be estimated, explaining that 33.3% of the contribution to the annual burden of Campylobacter (calculated as relative change in the predictive mean) was attributable to temperature, compared to 7.7% using a single coefficient value. The remaining regular seasonality across the five years was small, and was limited to holiday periods, perhaps reflecting surveillance bias. The contribution of rainfall to the incidence of Campylobacter was small. While being able to predict cases from the temperature in a much more complete way than the GEST model with a fixed coefficient, the GEST with the varying coefficient model did not provide an explanation of the causes of these changes, and suggested that other drivers linked to temperature were important.

When the GEST model was applied to weekly number of cases of Campylobacter from 1989 to 2009 in thirty three Strategic Health Authorities, it allowed the fitted seasonality to be used to group sub-regional SHAs into three geographical clusters.

The distribution of cases across geographically distinct sub-regions demonstrated strong differences in the strength of the seasonal pattern, but not the week when cases increase. The seasonality was strongest in Wales and the North East, and weakest and a bit more delayed in the large urban areas of London and Birmingham, with other geographic areas distributed in between. This implies that there are local factors that also influence Campylobacter disease occurrence. These could be demographic (e.g. more young people in cities), different risk profiles (e.g. more exposure to manures in the countryside; greater foreign travel in cities), and/or different impacts of weather locally (e.g. fly populations causing direct transmission to food or indirect transmission through the contamination of chicken flocks especially in rural areas) [17, 18, 21, 30].

Wavelet analysis showed that Campylobacter cases could be interpreted as the super-imposition of three main harmonics of annual, semi-annual and 4-months periods. The origin of these harmonics is not clear; however, annual and higher order (1/2 year, 1/3 year, etc.) harmonics, whose relative importance depends on the geographic region, are typical features of a range of climatic variables. Potential other drivers could be other environmental variables, reporting bias, human behaviour, and/or intrinsic seasonal patterns in animal/vector reservoirs (i.e. the contamination of chickens needs to be investigated). In this type of study the quality of data is important. We used laboratory postcode as a surrogate for patient residence postcode to anonymise the data and improve data completeness in linking cases to weather parameters in 30 days before the specimen date. Djennad et al. [27] demonstrated that it is valid to use laboratory postcode as a proxy for patient postcode. The algorithm to locally link case and weather was limited to a five-year period, for three temperature parameters and one rainfall parameter and had a lag of 30 days. Another limitation is that the data linkage was in cases only, with the result that weeks without cases were not recorded.

Being able to link cases to an LSOA allows the CCI to be applied to compare incidence using a range of weather parameters. This method utilised the geographic variation in weather parameters by combining data from small areas and their populations that had the same weather conditions. The CCI is limited by the weather data only being linked to cases only and not to the underlying population. While the approach provides a useful way of examining weather variables that are de-coupled from seasonality by grouping cases per week and their underlying populations at a particular temperature or rainfall level together as a combined incidence, it is clear that some weeks and areas will be left out because there are no Campylobacter cases. It would be expected that a measurement of the true incidence would give a result that differed somewhat from the CCI because the true incidence would include all the weeks with negative results not just those with positive ones. An extended version of this that uses weather parameters for both cases and populations has the potential to provide a more robust representation of the incidence of Campylobacter at different temperatures through the year.

The CCI and GEST approaches gave good results in examining environmental drivers for Campylobacter. Adding other weather parameters to these models should further improve the predictive power of these models or methods. As an example, the development of standard weather database and algorithms through the MEDMI project (https://www.data-mashup.org.uk) facilitates data linkage of environment and health data such that different analyses may be examined in a more systematic way for a wider range of infectious diseases [31]. The PHE Campylobacter dataset is sufficiently large, nevertheless the epidemiology of Campylobacter remains enigmatic and interventions have had limited success to date. It is therefore an important organism to study to elucidate methodologies that have power and utility in explaining infectious disease and climate change epidemiology for other infectious diseases in the future.

Conclusions

The study provides strong association between Campylobacter and temperature. Using a range of statistical methods, the study suggests that temperature and/or rainfall alone cannot explain the entire seasonal variation of Campylobacteriosis risk in England and Wales. Further research should investigate if the temporal dependency of the relationship between Campylobacter incidence and temperature on the week might be driven by other environmental variables, or perhaps by an intrinsic seasonality in the dynamics of the bacterial population in the environment or in the zoonotic reservoir or potential vectors such as flies.

Abbreviations

CCI:

Comparative Conditional Incidence

GEST:

Generalized structural time series model

GIS:

Geographic Information System

LSOA:

Lower Layer Super Output Areas

PHE:

Public Health England

SHA:

Strategic Health Authorities

WMVC:

Ward’s Minimum Variance Clustering

References

  1. Wallace JS, Stanley KN, Currie JE, Diggle PJ, Jones K. Seasonality of thermophilic campylobacter populations in chickens. J Appl Microbiol. 1997;82(2):219–24.

    Article  CAS  Google Scholar 

  2. Kovats RS, Edwards SJ, Charron D, Cowden J, et al. Climate variability and campylobacter infection: an international study. Int J Biometeorol. 2005;49(4):207–14.

    Article  Google Scholar 

  3. Nylen G, Dunstan F, Palmer SR, Andersson Y, et al. The seasonal distribution of campylobacter infection in nine European countries and New Zealand. Epidemiol Infect. 2002;128(3):383–90.

    Article  CAS  Google Scholar 

  4. Louis VR, Gillespie IA, O'Brien SJ, Russek-Cohen E, Pearson AD, Colwell RR. Temperature-driven campylobacter seasonality in England and Wales. Appl Environ Microbiol. 2005;71(1):85–92.

    Article  CAS  Google Scholar 

  5. Tam CC, Rodrigues LC, O'Brien SJ, Hajat S. Temperature dependence of reported campylobacter infection in England, 1989-1999. Epidemiol Infect. 2006;134(1):119–25.

    Article  CAS  Google Scholar 

  6. Nichols GL, Richardson JF, Sheppard SK, Lane C, Sarran C. Campylobacter epidemiology: a descriptive study reviewing 1 million cases in England and Wales between 1989 and 2011. BMJ Open. 2012;2:e001179. https://doi.org/10.1136/bmjopen-2012-001179.

  7. Fleury M, Charron DF, Holt JD, Allen OB, Maarouf AR. A time series analysis of the relationship of ambient temperature and common bacterial enteric infections in two Canadian provinces. Int J Biometeorol. 2006;50(6):385–91.

    Article  Google Scholar 

  8. Allard R, Plante C, Garnier C, Kosatsky T. The reported incidence of campylobacteriosis modelled as a function of earlier temperatures and numbers of cases, Montreal, Canada, 1990-2006. Int J Biometeorol. 2011;55(3):353–60.

    Article  Google Scholar 

  9. Naumova EN, Jagai JS, Matyas B, DeMaria A Jr, MacNeill IB, Griffiths JK. Seasonality in six enterically transmitted diseases and ambient temperature. Epidemiol Infect. 2007;135(2):281–92.

    Article  CAS  Google Scholar 

  10. Weisent J, Seaver W, Odoi A, Rohrbach B. The importance of climatic factors and outliers in predicting regional monthly campylobacteriosis risk in Georgia, USA. Int J Biometeorol. 2014;58(9):1865–78.

    Article  CAS  Google Scholar 

  11. Bi P, Cameron AS, Zhang Y, Parton KA. Weather and notified campylobacter infections in temperate and sub-tropical regions of Australia: an ecological study. J Infect. 2008;57(4):317–23.

    Article  Google Scholar 

  12. Rind E, Pearce J. The spatial distribution of campylobacteriosis in New Zealand, 1997-2005. Epidemiol Infect. 2010;138(10):1359–71.

    Article  CAS  Google Scholar 

  13. Spencer SE, Marshall J, Pirie R, Campbell D, Baker MG, French NP. The spatial and temporal determinants of campylobacteriosis notifications in New Zealand, 2001-2007. Epidemiol Infect. 2012;140(9):1663–77.

    Article  CAS  Google Scholar 

  14. Holland RC, Jones G, Benschop J. Spatio-temporal modelling of disease incidence with missing covariate values. Epidemiol Infect. 2015;143(8):1777–88.

    Article  CAS  Google Scholar 

  15. Jore S, Viljugrein H, Brun E, Heier BT, Borck B, et al. Trends in campylobacter incidence in broilers and humans in six European countries, 1997-2007. Prev Vet Med. 2010;93(1):33–41.

    Article  CAS  Google Scholar 

  16. Wei W, Schüpbach G, Held L. Time-series analysis of campylobacter incidence in Switzerland. Epidemiol Infect. 2015;143(9):1982–9.

    Article  CAS  Google Scholar 

  17. Nichols GL. Fly transmission of campylobacter. Emerg Infect Dis. 2005;11(3):361–4.

    Article  Google Scholar 

  18. Ekdahl K, Normann B, Andersson Y. Could flies explain the elusive epidemiology of campylobacteriosis? BMC Infect Dis. 2005;5:11.

    Article  Google Scholar 

  19. Hald B, Skovgård H, Bang DD, Pedersen K, Dybdahl J, Jespersen JB, Madsen M. Flies and campylobacter infection of broiler flocks. Emerg Infect Dis. 2004;10(8):1490–2.

    Article  Google Scholar 

  20. Hald B, Skovgård H, Pedersen K, Bunkenborg H. Influxed insects as vectors for campylobacter jejuni and campylobacter coli in Danish broiler houses. Poult Sci. 2008;87(7):1428–34.

    Article  CAS  Google Scholar 

  21. Hald B, Sommer HM, Skovgard H. Use of fly screens to reduce campylobacter spp. introduction in broiler houses. Emerg Infect Dis. 2007;13(12):1951–3.

    Article  Google Scholar 

  22. Ekdahl K, Andersson Y. Regional risks and seasonality in travel-associated campylobacteriosis. BMC Infect Dis. 2004;4(1):54.

    Article  Google Scholar 

  23. Ekdahl K, Giesecke J. Travellers returning to Sweden as sentinels for comparative disease incidence in other European countries, campylobacter and giardia infection as examples. Euro Surveill. 2004;9(9):6–9.

    Article  Google Scholar 

  24. Hartnack S, Doherr MG, Alter T, Toutounian-Mashad K, Greiner M. Campylobacter monitoring in German broiler flocks: an explorative time series analysis. Zoonoses Public Health. 2009;56(3):117–28.

    Article  CAS  Google Scholar 

  25. Sterk A, Schijven J, de Roda Husman AM, de Nijs T. Effect of climate change on runoff of campylobacter and Cryptosporidium from land to surface water. Water Res. 2016;95:90–102.

    Article  CAS  Google Scholar 

  26. Gillespie IA, O'Brien SJ, Bolton FJ. Age patterns of persons with campylobacteriosis, England and Wales, 1990-2007. Emerg Infect Dis. 2009;15(12):2046–8.

    Article  Google Scholar 

  27. Djennad A, Lo Iacono G, Sarran C, et al. A comparison of weather variables linked to infectious disease patterns using laboratory addresses and patient residence addresses. BMC Infect Dis. 2018;18(189):2–9.

    Google Scholar 

  28. Rösch A, Schmidbauer H. WaveletComp: Computational Wavelet Analysis. R package version 1.0, https://CRAN.R-project.org/package=WaveletComp. 2014.

  29. Djennad A, Rigby R, Stasinopoulos D, Voudouris V, Eilers P. Beyond location and dispersion models: the generalized structural time series model with applications. Munich Pers REPIC Repos. 2015:1–33.

  30. Guerin MT, Martin SW, Reiersen J, Berke O, McEwen SA, et al. Temperature-related risk factors associated with the colonization of broiler-chicken flocks with campylobacter spp. in Iceland, 2001-2004. Prev Vet Med. 2008;86(1–2):14–29.

    Article  CAS  Google Scholar 

  31. Fleming LE, Haines A, Golding B, Kessel A, et al. Data mashups: potential contribution to decision support on climate change and health. Int J Environ Res Public Health. 2014;11(2):1725–46.

    Article  Google Scholar 

Download references

Acknowledgements

We thank Dr. Robert Rigby, Professor Mikis Stasinopoulos, Professor Ben Armstrong and Dr. Shakoor Hajat for their insightful comments and feedback on statistical analysis and methodology.

Funding

The research was funded in part by the UK Medical Research Council (MRC) and UK Natural Environment Research Council (NERC) for the MEDMI Project; the National Institute for Health Research Health Protection Research Unit (NIHR HPRU) in Environmental Change and Health at the London School of Hygiene and Tropical Medicine in partnership with Public Health England, and in collaboration with the University of Exeter, University College London and the Met Office. It also included collaboration with the NIHR HPRU in Gastrointestinal Infection Programme. The funders had no role in the design of the study; collection, analysis, and interpretation of data; and in writing the manuscript. The views expressed are those of the author(s) and not necessarily those of the NERC, MRC, NHS, the NIHR, the Department of Health or Public Health England.

Availability of data and materials

The datasets on MEDMI are described at https://www.data-mashup.org.uk/data/data-library/ and include the SGSS infectious disease dataset. Permissions are required to access these datasets and users require an account to be set up as described at https://www.data-mashup.org.uk/data/accessing-data/.

Author information

Authors and Affiliations

Authors

Contributions

All authors were involved in the conception and design of the study and interpretation of the results. AD and TB were responsible for the mathematical modelling including time series analysis. GL was responsible for wavelet analysis and interpretation of the results. CS was responsible for temperature and rainfall linkages with Campylobacter. CL was responsible for LSOA linkage from postcode. AD, IL and FCG were responsible for cluster analysis. RE was responsible for producing the incidence map. CH and GN were responsible for cumulative incidence at different temperatures and CCI method. IL, JS, SK were responsible for epidemiological input. LEF and AK were the project consultants and lead. GN was responsible for the conception, design of the article and statistical input. All authors revised the manuscript critically and contributed to the final draft. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Abdelmajid Djennad.

Ethics declarations

Ethics approval and consent to participate

PHE has legal permission, provided by Regulation 3 of The Health Service (Control of Patient Information) Regulations 2002, to process patient confidential information for national surveillance of communicable diseases and as such, individual patient consent was not required.

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S5 shows: a) Average weekly reported Campylobacter cases averaged over 20 years (from 1989 to 2009). All time series were square root transformed and then normalised to sum to unity. c) wavelet power spectrum of the transformed time-series of Campylobacter. Low values of the power spectrum are shown in dark blue, and high values in dark red. The black lines show the maxima of the undulations of the wavelet power spectrum. The light white shaded areas identify the region subjected to errors arising from dealing with a finite-length time series (edge effect). e) global average wavelet power spectrum, the black dots show the 5% significant levels computed based on 100 bootstrapped series g) original and reconstructed time-series according to all harmonics and the selected first 3 harmonics only. b), d), f) h) As in figures a), c) e) and g) but after the time-series of Campylobacter cases were adjusted using a seven day rolling mean, removal of bank holiday artefacts and adjusted for long term trend. (PDF 1032 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Djennad, A., Lo Iacono, G., Sarran, C. et al. Seasonality and the effects of weather on Campylobacter infections. BMC Infect Dis 19, 255 (2019). https://doi.org/10.1186/s12879-019-3840-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-019-3840-7

Keywords