Article Text
Abstract
Objectives To identify key predictors of general practitioner (GP) consultations for allergic rhinitis (AR) using meteorological and environmental data.
Design A retrospective, time series analysis of GP consultations for AR.
Setting A large GP surveillance network of GP practices in the London area.
Participants The study population was all persons who presented to general practices in London that report to the Public Health England GP in-hours syndromic surveillance system during the study period (3 April 2012 to 11 August 2014).
Primary measure Consultations for AR (numbers of consultations).
Results During the study period there were 186 401 GP consultations for AR. High grass and nettle pollen counts (combined) were associated with the highest increases in consultations (for the category 216-270 grains/m3, relative risk (RR) 3.33, 95% CI 2.69 to 4.12) followed by high tree (oak, birch and plane combined) pollen counts (for the category 260–325 grains/m3, RR 1.69, 95% CI 1.32 to 2.15) and average daily temperatures between 15°C and 20°C (RR 1.47, 95% CI 1.20 to 1.81). Higher levels of nitrogen dioxide (NO2) appeared to be associated with increased consultations (for the category 70–85 µg/m3, RR 1.33, 95% CI 1.03 to 1.71), but a significant effect was not found with ozone. Higher daily rainfall was associated with fewer consultations (15–20 mm/day; RR 0.812, 95% CI 0.674 to 0.980).
Conclusions Changes in grass, nettle or tree pollen counts, temperatures between 15°C and 20°C, and (to a lesser extent) NO2 concentrations were found to be associated with increased consultations for AR. Rainfall has a negative effect. In the context of climate change and continued exposures to environmental air pollution, intelligent use of these data will aid targeting public health messages and plan healthcare demand.
- epidemiology
- public health
- primary care
- allergy
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
Strengths and limitations of this study
This study quantifies the impact of ambient concentrations of individual types of pollen and meteorological factors on general practitioner consultations for allergic rhinitis.
This study contributes to the small but growing body of evidence quantifying the role of pollen and air pollutants in the burden of allergic rhinitis, controlling for the effects of other environmental and meteorological factors.
This study demonstrates the potential utility of syndromic surveillance systems in supporting the monitoring of the health impact of pollen exposure.
The findings should be interpreted in the context of potential exposure measurement errors that are inherent in these types of studies.
There were insufficient data to include putative predictors (or confounders) such as sulfur dioxide, aggregate PM10, other pollen species and fungal spore measurements in the model.
Background
Allergic rhinitis (AR) is a global health problem affecting an estimated 500 million people worldwide.1 The symptoms can affect individuals’ sleep,2 performance at school3 and work,4 and AR has substantial financial costs to society; in Sweden, AR has been estimated to cost the economy approximately €2.7 billion annually.5 The prevalence of AR is increasing worldwide, with some countries describing over 50% of adolescents self-reporting compatible symptoms.6 In the UK, prevalence among adults has been estimated at 26% (95% CI 20.3% to 30.7%).7 Rhinitis is an inflammation of the lining of the nose characterised by nasal symptoms such as anterior or posterior rhinorrhoea, nasal blockage, sneezing and itching of the nose, and occurring during two or more consecutive days for more than 1 hour on most days. AR is the most common form of non-infectious rhinitis and the symptoms are caused by an IgE-mediated immune response to allergens. AR has a multifactorial aetiology that consists of genetic and environmental factors. Some of these triggers have been well described; and include outdoor allergens like pollen grains and moulds, and indoor allergens like mites, animal dander, insects and moulds.1
Further to the described association with allergens there is increasing evidence for the association between AR and three major air pollutants (namely PM10, sulfur dioxide (SO2) and nitrogen dioxide (NO2)),8–10 although this relationship has not been consistently demonstrated across all relevant studies. In 2001, Hajat et al8 investigated the relationship between air pollution and daily general practice consultations for AR in London using data from 1992 to 1994, and found an association between increasing consultations and ambient concentrations of SO2 and O3, which both contributed independently to the increase.
Despite the growing evidence for the role of air pollutants, the contribution of individual pollen types to consultations for AR has been less well described. Breton et al11 found a significant association between outpatient visits for AR and Ambrosia pollens in Montreal (Quebec, Canada). A Canadian study12 found an association between increasing levels of fungal spores or ragweed pollen and increases in paediatric emergency hospital visits for conjunctivitis and AR. They found that an increase of 551 basidiomycetes spores (fungus)/m3, or of 72 ragweed pollen grains/m3 was associated with a 10% increase in visits for conjunctivitis and rhinitis. In Beijing, Zhang et al13 found an association between high pollen counts and general practitioner (GP) consultations for AR, however, the authors measured a combined pollen count rather than measuring counts of individual pollen species. Less is known about the relative contribution of individual pollen types to the burden of GP consultations for AR.
In the wider context of a changing climate, continued environmental exposure to air pollution and the socioeconomic burden of AR, greater understanding of the main drivers for this illness is necessary to inform appropriate action at all levels. In this study, we aimed to estimate the contribution of putative meteorological and environmental factors on GP consultations for AR.
Methods
Study design
This was a retrospective, time series analysis of daily GP consultations for AR. The study population was all persons who presented to general practices participating in the Public Health England (PHE) national GP in-hours surveillance system within London over the study period (between 3 April 2012 and 11 August 2014). A case of AR was defined as a GP consultation episode where the GP assigned a Read code (the classification system currently used in the UK practice14) consistent with AR symptoms (codes included: AR, allergic rhinosinutis and hay fever (5-Byte Read codes H17; Read clinical terms V.3 XE0Y5)) and included in the daily extract used for PHE syndromic surveillance.
Data sources
GP consultation data
Within the UK, primary care is led by GPs, community-based doctors with generalist training. For patients, general practice is normally the first point of medical contact within the healthcare system; GP-led services provide open and unlimited access for patients, dealing with a wide spectrum of health problems. Services are available mainly through scheduled appointments (face to face or via telephone), with urgent unscheduled consultations also available.
The PHE GPIH (General Practitioner In Hours) system has a network of over 3500 GP practices across the UK covering approximately 35 million people (about 53% of the population) in England1 15 (within London the GPIH system provides coverage in 30 of the 31 administrative districts).16 Daily GPIH consultation data are routinely used to monitor the impact of public health events (eg, seasonal influenza) using clinical indicators (eg, influenza-like illness, diarrhoea, vomiting).
For this study, AR consultation data (daily numbers of consultations and denominator GP practice patient populations) were extracted for London for the period 3 April 2012 to 11 August 2014.
Meteorological data
Daily maximum, minimum and mean temperature (°C), and total precipitation (mm/day) were taken from a 5 km gridded dataset17 provided by the Met Office (the UK’s national weather service). The 5 km grid square containing the London King’s College pollen monitoring station was used in this study.
Pollen data
Daily counts of different pollen species; Betula (birch tree), Quercus (oak tree), Platanus (plane tree), Poaceae (grasses) and Urtica (nettle) measured as grains/m3 (g/m3) provided by the Met Office, for the London site of King’s College (latitude 51.510732, longitude −0.116938). The three tree pollens peak during a broadly similar time period and were aggregated into one variable, and grass/nettle pollens were similarly aggregated into one variable for inclusion in the model. The variables were combined as they broadly peaked during similar times and the combination reduced the covariance. Importantly, ‘out of season’ pollen species levels are considered negligible, and outside of these periods were not collected at the collecting stations and recorded. Periods when data were collected from the London site of King’s College (latitude 51.510732, longitude −0.116938) measuring station were 03/04/2012–19/08/2012, 14/04/2013–02/09/2013, 01/04/2014–11/08/2014. During 2012, earlier data were available from the pollen collecting station, but not included as GP data in the same format were unavailable.
Air pollutant data
Daily measurements of PM2.5 (microns), NO2 (ppb), ozone (O3) were extracted from the Department for Environment Food & Rural Affairs.18 More than 10% of data were missing for SO2 and PM10 so were not included. Data were taken from the London North Kensington monitoring station and used as a proxy for London. Less than 10% of data were missing for PM2.5 values and were filled using scaled data from the Bloomsbury urban background monitoring station (latitude 51.522290, longitude −0.125889), t-tests indicated that there were significant differences between stations and so missing data were imputed using linear regression to model the relationship between pollutant values at different stations. Residual analysis was used to ensure that there were no temporal changes in the relationships between sites.
Data analysis
Rates and proportion calculations
The expected number of AR cases E(Yt) at time t=1,…,T was modelled using a distributed lag nonlinear model.19 20 To account for potential over dispersion in the data, we tested negative binomial and quasi-maximum likelihood Poisson models and selected the model specification with the lowest dispersion parameter. To ensure the robustness of our results, we fitted our models using generalised linear modelling and generalised additive modelling approaches. The general algebraic definition for the negative binomial and Poisson models is given by:
where log(μt) is a logarithmic link function of the expectation E(Yt≡μt), with Yt as the time series of AR counts. The parameter α denotes the intercept; t corresponds to a linear function of time to account for potential long-term trends in the data. Seasonal trends are modelled using Fourier terms with 4 df/year. Long-term and seasonal trends were controlled for because they may be related to factors other than those considered in this study21 and has been used in pollen research previously.22 The terms w, b and h denote Boolean variables for the weekends (w), bank holidays (UK public holidays) (b) and school holidays (h). f(xt), g(zt) and k(pt) indicate the nonlinear and delayed effects of the meteorological, pollen and pollutant predictors, respectively. The relationship between AR and the meteorological, pollen and pollutant predictors is defined in two different ways. The first one corresponds to a 5-day moving average of each predictor time series in the set. The second one consisted of cross-basis functions describing simultaneously the nonlinear form of the relationship between AR and the predictor and its distributed lag effects.19 20 Cross-basis functions are implemented here using a distributed-lag nonlinear model which effectively avoids the collinearity issues resulting from the inclusion of a parameter for each lag of the exposure variable; algebraic definitions of the cross-basis functions are described in detail elsewhere.19 20
Results
Descriptive epidemiology
During the study period there were 186 401 GP consultations for AR. The mean week-day count of AR consultations was 218 with the range of AR consultations between 0 and 2046. Two distinct seasonal AR peaks were observed: the first peaking during May and the second, larger peak during June (figure 1).
Mean daily number of GP consultations per week, GPIH System, London, 3 April 2012 to 3 August 2014. GP, general practitioner.
Pollen species were plotted as counts (grains/m3)/day over the study period. Tree pollen species (plane, oak, birch) were more prevalent during the early summer months (April–May) while grass and nettle pollen were higher later in the year (July–August). The early AR peak coincided with tree pollens, with the late AR peak at broadly similar time to the grass and nettle pollen (figure 2). The rise in GP consultation rates appeared at broadly similar times to rises in pollen types during each year.
Daily mean pollen counts/week (by species, grains/m3), London, 3 April 2012 to 3 August 2014.
Modelling
After testing 10 different model specifications, a negative binomial model incorporating daily mean temperature, daily total rainfall, daily aggregated tree pollen, daily aggregated grass and nettle pollen, daily NO2, ozone and PM2.5 were selected. The model explained ~93% of the deviance in the GP consultations for AR and showed a dispersion parameter close to the unity (1.04), indicating that the negative binomial specification was adequate.
The average estimates of the overall cumulative exposure–response relationships estimated by the model were plotted (figure 3). The relationships shown were interpreted as the accumulated effect of the predictors on AR over the current and previous 4 days under the assumption that the predictor remained constant over the entire 5-day period.
Estimated relationships between allergic rhinitis (AR) GP consultations, temperature, rainfall, aggregated tree and grass/nettle pollen and air pollution from a distributed lag nonlinear model. The y axis (RR or relative risk) represents the excess in AR attendances with respect to the mean values of each covariate. The thick smooth lines indicate the maximum likelihood estimates, and the grey shaded areas represent the 95% CIs. GP, general practitioner.
The results illustrated an overall picture of the delayed effect of temperature on AR compared with a reference value of 10°C. The relative risk (RR) of an AR consultation gradually increased with increasing temperatures up to approximately 17.5°C, after which the RR gradually decreased, becoming non-significant at temperatures above 21°C. The overall mean estimated RR versus the 10°C reference value was 0.55 (95% CI 0.51 to 0.61) for temperatures below 10°C, and 1.36 (95% CI 1.13 to 1.66) for temperatures above 10°C. The estimates of the RR (and their 95% CIs) of AR were all significant at all levels of temperature except for temperatures between 20°C and 25°C (table 1).
Daily mean number of AR consultations and RR by metrological, pollen and pollutant variables
A delayed negative association between rainfall and AR consultation was found, with the RR decreasing as rainfall increased. The mean RR of AR at values of rainfall below the 10 mm day threshold was 1.29 (95% CI 1.16 to 1.45) and for values above this threshold the RR was 0.87 (95% CI 0.78 to 0.97). The effect of rainfall on AR was statistically significant at all levels evaluated in this study (table 1).
There was a positive relationship between pollen counts (combined tree, grass and nettle) and AR consultations with the number of GP consultations increasing with increasing pollen concentrations. Tree pollen concentrations below 30 grains/m3 produced an average RR of 0.95 (95% CI 0.91 to 0.98), whereas tree pollen concentrations above the same threshold generated an average RR of 1.41 (95% CI 1.18 to 1.70). Grass/nettle pollen had an average RR of 0.86 (95% CI 0.83 to 0.90) at concentrations below 30 g/m3 and of 2.29 (95% CI 1.95 to 2.71) at concentrations above that threshold. The effect of concentrations of grass/nettle pollen above 90 grains/m3 produced a considerable larger effect on AR than that generated by tree pollen (table 1).
The magnitude of the associations between AR consultations and air pollution differed depending on the pollutant examined. There was a direct relationship between NO2 and AR consultations with values between 30 and 45 ppb having an RR of 1.08 (95% CI 1.01 to 1.15) rising to an RR of 1.33 (95% CI 1.03 to 1.71), and the RR for NO2 levels below the reference value (30 ppb) having had an average RR of 1.03 (95% CI 00.93 to 1.16).
The effects of NO2 on AR became non-significant at values 85 times above the reference value. The effects of ozone on AR were non-significant. Finally, the effects of PM2.5 on AR were significant only at values lower than the 30 times the reference (table 1) with an average RR of 1.23 (95% CI 1.07 to 1.41).
The model was then re-fitted after disaggregating the two aggregated pollen variables (tree and grass/nettle) into their sub-components (figure 4). The estimated relationships were very similar to those obtained with the aggregated pollen data for most pollens although there were wider CIs around some of the point estimates (eg, Betula); as was expected with a smaller sample size. However, the association between AR consultations and the Oak tree pollen (Quercus) was inverse, a finding that was different to that observed with the other pollen species in the analysis. A closer examination revealed that high concentrations of Quercus pollen always coincided with low concentrations of other tree pollens, which may account for the individual observation.
Estimated relationships between allergic rhinitis (AR) GP consultations, temperature, rainfall, tree, grass and nettle pollen, and air pollution from a distributed lag nonlinear model. The y axis (RR or relative risk) represents the excess in AR attendances with respect to the mean values of each covariate. The thick smooth lines indicate the maximum likelihood estimates, and the grey shaded areas represent the 95% CIs.
Discussion
This study has assessed the relationship between GP consultations for AR and several different variables that are implicated in the overall burden. The observed seasonal pattern of AR consultations is broadly consistent with other studies,14 with a large peak occurring in June, preceded by an earlier, smaller peak, usually during April or May. It is apparent from both existing literature and from the results of this study that the fluctuations in AR consultations coincide broadly with pollen peak levels, and that the observed seasonality is predominantly driven by weather and climate drivers, and differences between the flowering times of different types of trees/grasses/weeds23 (also known as the ‘pollen calendar’).
This study has quantified the association between individual pollen species and AR consultations. The first peak of GP AR consultations was associated with increased tree pollen levels (Platanus, Betula, Quercus), with the second peak of consultations appearing to be associated with increased grass (Poaceae) and nettle (Urtica) levels. These results are broadly consistent with the known seasonal distribution of pollen, where in general tree pollens are higher during the earlier part of the summer season, with grass/nettle pollens appearing later in the season.23 Interestingly, when similar amounts of tree or combined grass and nettle pollens are present, the RR of a GP consultation is higher when exposed to the grass/nettle pollen. The biological reasons for this may include differential effects of some pollens (grasses have been implicated in asthma exacerbations leading to hospitalisations, probably via allergenic asthma hospitalisations).22 and there are variations in allergenicity of pollen between trees and grasses.24 A further finding was in the disaggregated analysis, nettle (Urtica) pollen had a significant effect on GP consultations; these are not usually considered strong allergens, and further work is necessary to better understand this relationship.
Importantly, the peaks in the AR consultations appeared broadly at similar times to the rises in the observed pollen species. Other studies have described differing lag times between exposure to a putative exposure and occurrence of health effects. Hajat et al8 identified significant increases in consultations 3 and 4 days following exposure to SO2 and O3, and Breton et al11 demonstrated strong correlations between pollen counts and higher medical consultation for AR on days 0–3 and 5 days after elevated pollen exposure, as did Osborne et al.22 The spikes in the number of AR consultations appear to precede the rises in the observed pollen species, our findings did not have a lag (day 0). There are several possible reasons for this including the incomplete or delayed capture of pollen at the monitoring sites; and unmeasured variables (which may rise at the same time as the measured variables). Other explanations include potential effectiveness of alerts or media messaging reporting on pollen encouraging patients to seek medical help in advance of symptoms, or weather forecasts/prevailing conditions initiating ‘prophylactic’ visits to a GP to obtain medication in advance of symptoms. It is likely that these would be recorded with similar diagnostic codes to trigger a consultation being coded as AR in the syndromic surveillance system.
Despite the limitations that this may present in interpreting the results of this study, importantly, from a healthcare systems perspective, this raises the possibility that syndromic surveillance may provide an early warning increasing healthcare activity during periods of rising pollen to the wider healthcare system enabling both timely public health messaging and proactive healthcare planning.
The interpretation of the impact of air pollutants on AR consultations is complex. Pollutants such as carbon monoxide and nitrogen oxide have been reported to contribute to the burden of both AR and asthma.9 25 Our study finds this relationship to be direct and significant for NO2 when accounting for pollen levels, although the impact is clearly less than of the pollen. However, there is an unusual relationship observed between consultations and PM2.5; with a significant effect observed at very low levels, but not replicated at higher loads; the lack of both biological plausibility and dose–response might suggest that finding is an artefact. The lack of effect of ozone has been observed in other studies,26 but is complicated by ozone levels often being low in areas where primary emission concentrations are high, and it has been postulated that no (or even a protective) effect from ozone should not be concluded: instead providing further evidence for the impact of other air pollutants.26
The observed relationship with temperatures above 20°C not being independently associated with increased consultation rates warrants further investigation to understand the reasons behind it. Potentially behavioural factors, where there is a higher propensity for people to be outdoors in warmer temperatures, but as they increase over a certain level move indoors to avoid may have a role, or there could biological factors linked to the release of pollens at different temperatures, such as different species of grass pollen being released at hotter temperatures.27
A further interesting finding is that rainfall has an inverse relationship with AR—as rainfall increases AR consultations decreases. The influence of rainfall on pollen concentrations was initially considered to be negative; rainfall washing out airborne pollens from the atmosphere28 and has been replicated recently,29 however, other studies have found the situation to be more complex30 31; and the influence of rainfall in removing pollen from the air is perhaps related to the size of pollen grains for different taxa.32 It has been postulated that the absorption of water by pollen grains can lead to them shattering and releasing starch granules, leading to greater exposure with smaller more numerous particles being breathed deeper into the lungs, as in the hypotheses surrounding thunderstorm asthma. Behavioural factors may also be influential; with increasing rainfall people are less likely to spend significant periods of time outside.
There were a number of limitations to this study. First, the results need to be considered with the possibility of measurement bias arising from errors in measuring the exposures. Second, the geographical area used to gather data on GP consultations (upper tier local authority) was much larger than area covered by the few monitoring sites for the environmental exposures; this was a pragmatic approach as only one monitoring site could provide daily pollen count data, but we recognise that the actual individuals’ exposure from pollen and air pollutants may have been very localised. Third, GP consultations may be biased towards younger children who are able to obtain free medication, and patients presenting with more severe or persistent symptoms. It is likely that most AR cases in the community will ‘self-treat’ or ‘self-medicate’ with antihistamines available without prescription from a pharmacy, and therefore this work may underestimate the overall burden of disease. While this may affect the external validity of our study findings (ie, generalisability to population groups not captured in GP AR data) it is, however, unlikely to affect internal validity. Furthermore, pollen concentration is considered low outside of the usual season and is not measured outside the season, so we were unable to include any out-of-season pollen measurements in our analysis. Finally, there were insufficient data to include either sulfur dioxide, aggregate PM10, other pollen species (in particular Alder and Hazel), and fungal spore measurements in our model; potential independent predictors (and/or confounders) that we were unable to account for in this study, as well as potential synergistic effects of the exposures which we did not explore.
Future work building on the results of this study could incorporate data on other pollen species or fungal spores into the model, and develop forecasting models to provide near real-time predictions of healthcare burden and improve the timeliness of public health messaging during the AR seasons. Recent work by Brennan et al27 shows that cutting edge DNA ‘barcoding’ techniques are now being used to detect individual grass species in the air. Given the greater effect of AR GP consultations found in our study, the possibility to repeat this analysis with counts of individual grass species could lead to new findings about the health impact of specific species of grasses. In the context of climate change and continued exposure to environmental pollutants, intelligent use of these data will be crucial to target public health messaging and plan healthcare demand.
Acknowledgments
We thank and acknowledge the University of Nottingham, ClinRisk and the contribution of EMIS and EMIS practices; TPP, ResearchOne and the SystmOne GP practices contributing to this surveillance system. We also thank the PHE Real-time Syndromic Surveillance Team for technical support, in particular Paul Loveridge. We acknowledge the use of Met Office meteorological data, the London North Kensington and Bloomsbury measuring stations (Information provided by the National Meteorological Library and Archive—Met Office, UK). Observational pollen information provided by the National Meteorological Library and Archive—Met Office, UK. Observed pollen data were collected by London King’s College and we thank all the individual collectors at the site for this. We also thank the Met Office who have operated the UK pollen network since 2011.
References
Footnotes
Twitter @nosborne36, @SotirisVard
Collaborators Paul Loveridge.
Contributors DT co-ordinated the project, provided initial analysis and drafted the manuscript. FdJCG completed the statistical analysis. RMo and AC also provided statistical support. SH, SK, NJO, RMc, SV and KE provided expertise in environmental factors. RMc organised data acquisition for environmental variables. OE, GS and AJE provided expertise in Syndromic Surveillance, and organised data acquisition for syndromic data. SK, GS, DT and AJE were involved in the project initiation. All authors contributed to, and commented on the protocol, methodology and final manuscript.
Funding The research was funded through PHE core activities, and supported in part by the National Institute for Health Research Health Protection Research Unit (NIHR HPRU) in Environmental Change and Health (SH, SK, NJO, RMc) at the London School of Hygiene and Tropical Medicine in partnership with PHE, and in collaboration with the University of Exeter, University College London, and the Met Office. FdJCG, RMo, OE, GS and AJE were supported by the NIHR HPRU in Emergency Preparedness and Response at King’s College London in partnership with PHE. Funding for open access publication was through the NIHR HPRU in Environmental Change and Health. The views expressed are those of the authors not necessarily those of the NHS, the NIHR, the Department of Health and Social Care or Public Health England.
Competing interests None declared.
Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.
Patient consent for publication Not required.
Provenance and peer review Not commissioned; externally peer reviewed.
Data availability statement Data may be obtained from a third party and are not publicly available. Syndromic surveillance data used in this manuscript originate from the University of Nottingham, ClinRisk EMIS and EMIS practices; TPP, ResearchOne and the SystmOne GP practices. Metrological and Pollen information was provided by the National Meteorological Library and Archive—Met Office, UK.