Objectives To investigate associations between daily concentrations of air pollution and myocardial infarction (MI), ST-elevation MI (STEMI) and non-ST-elevation MI (NSTEMI).
Methods Modelled daily ground-level gaseous, total and speciated particulate pollutant concentrations and ground-level daily mean temperature, all at 5 km×5 km horizontal resolution, were linked to 202 550 STEMI and 322 198 NSTEMI events recorded on the England and Wales Myocardial Ischaemia National Audit Project (MINAP) database. The study period was 2003–2010. A case-crossover design was used, stratified by year, month and day of the week. Data were analysed using conditional logistic regression, with pollutants modelled as unconstrained distributed lags 0–2 days. Results are presented as percentage change in risk per 10 µg/m3 increase in the pollutant relevant metric, having adjusted for daily mean temperature, public holidays, weekly influenza consultation rates and a sine-cosine annual cycle.
Results There was no evidence of an association between MI or STEMI and any of O3, NO2, PM2.5, PM10 or selected PM2.5 components (sulfate and elemental carbon). For NSTEMI, there was a positive association with daily maximum 1-hour NO2 (0.27% (95% CI 0.01% to 0.54%)), which persisted following adjustment for O3 and adjustment for PM2.5. The association appeared to be confined to the midland and southern regions of England and Wales.
Conclusions The study found no evidence of an association between the modelled pollutants (including components) investigated and STEMI but did find some evidence of a positive association between NO2 and NSTEMI. Confirmation of this association in other studies is required.
- Air pollution
This is an Open Access article distributed in accordance with the terms of the Creative Commons Attribution (CC BY 4.0) license, which permits others to distribute, remix, adapt and build upon this work, for commercial use, provided the original work is properly cited. See: http://creativecommons.org/licenses/by/4.0/
Statistics from Altmetric.com
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.
What is already known about this subject?
Evidence from epidemiological studies examining associations between short-term variations in concentrations of ambient air pollution and acute myocardial infarction (MI) is mixed.
Most studies have focused on urban populations—few have explored the wider geographical coverage (especially for exposure to PM2.5 and its components) offered by the use of atmospheric chemistry transport models (ACTM).
What does this study add?
Using ACTM data, this study found no evidence of an association between the pollutants investigated (O3, NO2, PM2.5, PM10 and the PM2.5 components sulfate and elemental carbon) and MI or ST-elevation MI.
However, there was evidence of a positive association between NO2 (a traffic related pollutant) and non-ST-elevation MI which appeared to be independent of PM2.5 and O3.
How might this impact on clinical practice?
Our findings add to the growing epidemiological literature investigating the acute effects of air pollution on cardiovascular health.
The identification of susceptible population subgroups most at risk on high pollution days would enable appropriate advice to be formulated and communicated in a timely manner to reduce the risk to health.
Air pollution has been associated with adverse cardiovascular heath events in long-term exposure (cohort) studies and in short-term (time-series or case-crossover) studies, and there is mounting evidence from clinical investigations as to potential mechanisms.1–4
In terms of the more specific outcome of acute myocardial infarction (MI), epidemiological studies examining the short-term effects of outdoor ambient air pollution have varied in their findings.5–10 While a minority of studies suggest an increase in the risk of an MI following exposure to higher concentrations of ozone (O3), evidence from the literature of an increase in risk following a rise in exposure to the mass of particulates with an aerodynamic diameter <10 µm (PM10), or <2.5 µm (PM2.5), nitrogen dioxide (NO2), carbon monoxide (CO) and sulfur dioxide (SO2) is more consistent.5–10 Part of the reason for inconsistencies between studies may be due to the different outcome measures used (from emergency department visits and hospital admissions to mortality), the accuracy of diagnosis or different assumptions about the nature of any association (linear or non-linear in the log relative risk). Results may also be influenced by the relative proportion of MI subtypes in the study population, as well as the source of gasses and the type and source of particle species. The few studies that have presented results for ST-elevation MI (STEMI) and non-ST-elevation MI (NSTEMI) separately have produced conflicting results.11–13 An analysis of data (based on 452 343 MI events between 2003 and 2009) from the England and Wales Myocardial Ischaemia National Audit Project (MINAP)14 reported small positive associations of NO2 (unconstrained distributed lag model (UDLM) lags 0–4) with all MI and NSTEMI,11 while a much smaller study in the USA reported a positive association of PM2.5 (lag 1 hour) with STEMI.12
Previous studies have also varied in the accuracy, completeness and representativeness of the pollution data, the choice of pollution metric and exposure period/lag, as well as to the choice and modelling of potential covariates. Most short-term studies to date, including those based on MINAP data, have used pollution monitoring networks to provide their daily exposure assessments. Few have so far explored the potential advantages of using modelled data that have complete coverage over urban and rural areas, diversity (in terms of particle species components) and enable air quality-related policy scenarios to be investigated. The EMEP4UK atmospheric chemistry transport model (ACTM)15 ,16 is a high-resolution regional application of the well-established European Monitoring and Evaluation Programme (EMEP) MSC-W model.17 It simulates the evolution of ambient pollution concentrations using official pollutant emission inventories, relevant natural emissions and driving meteorology, through a detailed treatment of atmospheric chemistry and physics. In this paper, we use ACTM data at 5 km by 5 km spatial resolution to explore associations between modelled ground-level gaseous, total and speciated particulate pollutant concentrations linked to MI, STEMI and NSTEMI events recorded on the England and Wales MINAP database between 2003 and 2010.
The outcome data for the study period of 2003–2010 come from MINAP. This is a register of hospital admissions for acute coronary syndromes (ACS) covering all acute National Health Service hospitals in England and Wales. In addition to discharge diagnosis, the register contains over 100 separate fields of patient-level data, including demographic information (eg, sex, age, smoking status, ethnicity), medical history (eg, history of MI, cerebrovascular disease, peripheral vascular disease, COPD), drug treatment prior to and during the admission and clinical findings (eg, ECG results, symptoms).14
For the purpose of this study, we excluded patients with missing geocodes, insufficient information on date of event, missing information on discharge diagnosis or linked to areas outside England and Wales. This left 630 116 events occurring during 2003–2010 and recorded on the MINAP database, of which 203 804 had a discharge diagnosis of STEMI and 323 999 had a discharge diagnosis of NSTEMI (including troponin positive ACS). New left bundle branch block was treated as synonymous with ST-elevation in the diagnosis of STEMI.18
Model pollution and weather data
The pollution data are surface daily outputs (derived from hourly outputs) at 5 km×5 km spatial resolution from an ACTM: EMEP4UK version rv4.3.15 ,16 High resolution is achieved through a nested approach whereby a 5 km×5 km (inner) domain over the British Isles is nested within, and takes boundary and initial conditions from, a larger 50 km×50 km (outer) domain over Europe. The meteorology driving the ACTM comes from the Weather Research and Forecasting (WRF) model version 3.1.119 which is constrained to contemporary meteorological observations, ensuring that the applied meteorology is representative of the real weather conditions prevailing throughout the simulated period. Pollutant and pollutant–precursor emissions over the UK were taken from the UK National Atmospheric Emissions Inventory20 and for the outer domain from EMEP estimates provided by the Centre for Emission Inventories and Projections.21
The pollutant metrics investigated in this study were daily means of PM2.5, PM10, sulfate (SO42−) and elemental carbon (EC) and daily maximum 8-hour-running mean for O3 and daily maximum 1-hour mean for NO2. The weather metric was daily mean temperature (see below).
Monitor pollution data
Monitor data were used only in the assessment of model performance (see below and online supplementary table S1). For this purpose, daily maximum 8-hour running mean O3, daily maximum 1-hour NO2, daily mean PM2.5 and daily mean PM10 were calculated for 2001–2010 using a 75% data capture threshold on hourly data from urban background and rural monitoring sites of the Automatic Urban and Rural Network (AURN) of the UK Department for Environment, Food and Rural Affairs.22 PM2.5 data were not available for the full period of interest, as the monitoring of this pollutant began only in the latter part of the decade. Between 2001 and 2010 there were changes in the instrumentation used to monitor PM10.
Validation of the pollution and weather models
The EMEP4UK model has undergone extensive validation,15 ,16 ,23 and their potential for use in epidemiological analyses is beginning to be explored.24 For 2001–2010 the average Pearson correlation over time between urban background monitored pollution concentrations at AURN monitoring sites,22 and their equivalents for the EMEP4UK model grid incorporating the monitor, was relatively high for daily maximum 8-hour O3 (no. of sites (n)=63; =0.76; SD (r)=0.04) and daily mean PM2.5 (n=39; =0.69; SD (r)=0.09) and lower for daily mean PM10 (n=57; =0.50; SD (r)=0.07) and maximum 1-hour NO2 (n=75; =0.54; SD (r)=0.10).
Imprecision in the estimation of daily pollution concentrations whether modelled or measured may on average lead to some attenuation (ie, bias towards the null) in estimates of the log relative risk obtained from epidemiological analyses. Simple predictions as to the level of attenuation (ie, percentage bias towards the null) that might be expected due to the use of ACTM data are explored in online supplementary table S1.
Since monitored temperature data are available for only a small subset of the model grid boxes, we used the WRF model 2-metre temperature for covariate adjustment. As discussed above, the WRF model was nudged with reanalysis data every 6 hours to closely represent observations such as the surface temperature.
For daily mean temperature, the average model-monitor correlation over time was very high (n=93; =0.98; SD (r)=0.01) and plots of the relationship between MI and modelled temperature (see online supplementary figure S1) are similar to those previously published by Bhaskaran et al,25 using the MINAP database and monitored temperature.
Each MI event was linked to the modelled weather and pollutant exposure data in the 5 km grid closest to the output area (OA) of the patient's residence (using OA centroids rounded to 1000 m to avoid personal identification).
The analysis was conducted at the level of the individual using a time-stratified case-crossover analysis.26 For each case we defined the index day as the day of the MI and the referent or control days as those days within the same month and on the same day of the week as the event day.26 Within each individual we then compared the modelled pollutant exposures between the index and referent days as in a 1:M matched case–control study. The analysis was conducted in STATA V.12 (StataCorp: Stata Statistical Software: Release 12. College Station, TX: StataCorp LP; 2011) using conditional logistic regression. In this way each subject acted as their own control, automatically adjusting for potential non-time varying/time-insensitive confounders such as sex, age, smoking status, socioeconomic status etc. The matched sets or strata were defined in time (ie, year, month, day of the week) to remove trend and seasonal pattern and additional covariates were added to the regression models to adjust for temperature, public holidays, influenza epidemics and residual seasonality. The primary regression model included: two natural cubic splines (each with 5 degrees of freedom (df)) representing mean daily temperature averaged across the day and the day before (mean lag 0–1) and mean daily temperature averaged across days 2–6 before (mean lag 2–6); a binary indicator variable for public holidays; the Royal College of General Practitioners (RCGP) England and Wales weekly consultation rate for influenza-like illness for the week of the event;27 and sine and cosine terms representing a simple annual cycle. A priori, pollutants were included in analyses as unconstrained distributed lags 0–2 days (UDLM 0–2). Under the rare disease assumption, ORs from conditional logistic regression were interpreted as relative risks and are presented as such in tables and plots with their 95% CIs.
Possible effect modification by season (autumn=September–November, winter=December–February, spring=March–May; and summer=June–August) and, where applicable, by sex and age group (≤64, 65–74, 75–84, ≥85) was investigated by including appropriate interaction terms in the regression models and testing for an improvement in fit using likelihood ratio tests. However, when investigating effect modification by Government Office Region (10 in England and Wales), a two-stage analysis was employed whereby ORs from region-specific conditional logistic regressions were meta-analysed using METAN in STATA V.12, (StataCorp LP; 2011) to obtain overall relative risks (across all regions), subtotal relative risks for each of three broader areas or ‘super regions’ (ie, the North, Midlands and the South) and tests of heterogeneity between regions and between ‘super regions’. This approach facilitated region-specific covariate adjustment.
Underlying our analyses is the assumption that any association between MI and pollution is approximately log-linear. This assumption was investigated in sensitivity analysis by fitting natural cubic splines (each with 2 df) to simple pollutant averages (averaged over lags 0–2) and testing for non-linearity using the Wald χ2 tests (df=1).
Having excluded records with missing data on exposures or covariates, our main analysis was based on 626 239 events occurring during the study period (ie, 2003–2010) and recorded on the MINAP database. The number of control days ranged from 1 to 4 per case but with <1% of cases matched to only 1 or 2 controls. Of these 626 239 events, 202 550 were diagnosed on discharge as STEMI, 322 198 as NSTEMI and 37 579 as ACS (troponin negative). The remaining 63 912 had discharge diagnoses, including threatened MI, chest pain of uncertain cause, unconfirmed MI, Takotsubo cardiomyopathy, ACS (troponin unspecified), PCI-related MI and other (ie, admitted with clinical suspicion of cardiac pain but diagnosis other than cardiac ischaemia confirmed). Our definition of ALL MI includes all 626 239 events and is therefore broad and inclusive. The median (IQR) age, in years, of cases was 70.6 (59.1–80.0) for ALL MI, 66.2 (56.0–76.5) for STEMI and 74.2 (63.2–82.4) for NSTEMI and the percentage (n/N) male was 65% (404 374/624 926), 70% (141 971/202 024) and 62% (199 216/321 633), respectively.
Association with pollution
In our primary logistic regression models adjusted for temperature lags 0–1 and 2–6 days, influenza, public holidays and residual seasonality (table 1), there was no evidence of an association of ALL MI or STEMI with any of O3, NO2, PM10 or PM2.5. For NSTEMI, there was a statistically significant (p=0.043) positive association with NO2 suggestive of a 0.27% (95% CI 0.01% to 0.54%) increase in risk per 10 µg/m3 increase in the daily 1-hour maximum concentration. This association did not differ significantly by season (p=0.571), age-group (p=0.986) or sex (p=0.859). It also persisted following adjustment for O3, and increased in magnitude following adjustment for PM2.5 (p=0.008) (table 1: two pollutant models).
Effect modification by season
Returning to single pollutant models, evidence of modification by season (table 2) was observed for PM2.5 in relation to ALL MI and for O3 in relation to STEMI. For O3 there was evidence of a negative association with STEMI in the autumn and for PM2.5 there was evidence of a positive association with ALL MI in the autumn. However, in both cases, there was no consistency in the direction of relative risk estimates across the four seasons.
Effect modification by region
We next investigated differences in association between Government Office Regions, having adjusted for covariates at the regional level (table 3). With respect to NO2, PM10 and PM2.5, regression models included additional regional-level adjustment for ozone (lags 0, 1 and 2).
For ALL MI and for STEMI there was no evidence of heterogeneity in association between regions or between ‘super regions’ (North, Midlands, South) whether in relation to O3, NO2, PM2.5 or PM10. In contrast for NO2 and NSTEMI, subtotal relative risks for the Midlands and the South were significant and >1, suggestive of increases of 0.79% (0.13% to 1.45%) and 0.51% (0.09% to 0.94%) per 10 µg/m3, respectively, whereas in the North the subtotal relative risk was <1 and non-significant (test for heterogeneity between ‘super regions’, p=0.022). For O3 and NSTEMI, associations differed significantly between regions (p=0.006) but not between ‘super regions’.
When we compared the overall estimates (ie, over all regions) in table 3 with their corresponding values in table 1, it appeared that adjustment for region had little effect.
Table 4 displays the results from investigating the associations between our three outcomes and two selected PM2.5 components. The two components tabulated are EC (emitted from combustion) and sulfate (SO42−, formed within the atmosphere from SO2). The conditional logistic regression models adjust for temperature at lags 0–1 and 2–6 days, public holidays, influenza and residual seasonality. Neither EC nor SO42− was associated with ALL MI, STEMI or NSTEMI either before or after additional adjustment for O3 at lags 0, 1 and 2.
As part of our sensitivity analyses, we returned to our single pollutant regressions (table 1) and plotted the individual components of each UDLM as the percentage change in risk per 10 µg/m3. From figure 1, it appeared that the association between NO2 and NSTEMI was primarily due to concentrations on the day (lag 0) rather than on the 1 day (lag 1) or 2 days prior (lag 2).
Next we looked for evidence of non-linearity in the associations between simple pollution averages (across lags 0, 1 and 2) and the log relative risks of ALL MI, STEMI and NSTEMI. Having adjusted for covariates, we found no evidence of non-linearity based on Wald tests with 1 df (all p>0.19).
Finally, to investigate the sensitivity of our findings to the level of control for seasonality, we re-ran our single pollutant regression models without the sine and cosine terms and observed little change in risk estimates and CIs (see online supplementary table S2).
We investigated associations between gaseous and particulate air pollution and three outcomes, and found no evidence of associations except for one positive association between modelled NO2 and NSTEMI, suggestive of a 0.27% increase in NSTEMI per 10 µg/m3 increase in daily maximum 1-hour NO2 (UDLM 0–2). This association increased in magnitude and statistical significance following adjustment for PM2.5. While not significant in any one season, the relative risk for NO2 and NSTEMI in our single pollutant model was consistently >1 in the autumn, winter, spring and summer months (table 2). However, this consistency was not replicated in our regional analysis, which suggested that any positive association between NO2 and NSTEMI was confined to the Midlands and the South (table 3). The lack of association with EC and SO42− is noteworthy given their novelty and what they represent.
Comparison with other studies
A recent systematic review of air pollution and MI10 identified 34 time-series/case-crossover studies published between 1997 and 2011 for inclusion in meta-analyses. The studies were mainly from the USA and Europe (including two from the UK) and varied in their source of outcome data (eg, hospital admissions, death registries, MI registries). In common with our results, the review found no evidence of an association between MI and O3 (lag 0). However, it reported statistically significant increases in the risk of MI of 1.1%, 0.6% and 2.5% per 10 µg/m3 increases in NO2 (lag 1), PM10 (lag 0) and PM2.5 (lag 1), respectively. These represent very different relative risks to those observed (UDLM 0–2) in our study (0.09%, −0.20%, −0.04%, respectively). In particular, the increase in risk for PM2.5 is not apparent in our study. One possible explanation for the difference between our finding and those from the review (besides obvious differences of time-period, geography and outcome measure) is our use of modelled rather than monitored pollutant data. We, therefore, compared our findings with those of a previously published analysis of the MINAP database (2003–2009) by Milojevic et al,11 which investigated associations between MI and nearest monitor pollution concentrations (up to 50 km distant from residential address). In common with our own results, they found no evidence of positive associations with O3, PM2.5 or PM10 but a significant positive association of daily mean NO2 (lags 0–4) with NSTEMI, although their estimated increase in risk was more marked (0.68% per 10 µg/m3 increase in NO2 (95% CI 0.07% to 1.3%)). One possible explanation for this difference is the different uncertainties in the monitor-based pollution data used by Milojevic et al (measurements made up to 50 km from residential address) and the model-based exposure assessment used here (see online supplementary table S1).
Strengths and limitations
One of the strengths of using MINAP data is that it provides us with a large number of events and hence good statistical power for detecting associations with ALL MI but also with the MI subtypes of STEMI and NSTEMI. However, given the consideration of three outcomes and six pollutants in our analyses, the possibility of obtaining a spurious significant association cannot be discounted.
We acknowledge the potential for modelled pollution concentrations based on any ACTM to introduce more classical measurement error into a time series analysis than monitor data and hence lead to greater attenuation in relative risk estimates. However, based on our comparison of model and monitor data in online supplementary table S1 and our previous statistical simulation study,24 we expect this to be less of an issue for O3 and for urban PM2.5. Measurement error (classical and Berkson) may also lead to a reduction in statistical power, although, for ACTMs, the ability to provide more complete pollutant time series may partly compensate by limiting any reduction in sample size due to missing data.
One advantage of using modelled pollution data is that we can obtain geographically complete data, as measurements of PM2.5 and other species are limited (in space and time) over the UK. In addition, models can provide particulate matter component data, which, in England and Wales, are sparsely monitored, or restricted to intensive field campaigns, and hence for which little historical data are available. While, in this study, we found no evidence of associations with the EC or sulfate components of PM2.5, there may be associations with other particle component matter not yet investigated.
Critical exposure period
If there is a detrimental effect of some pollutants, it is not clear whether these effects are immediate (in the few hours prior to the event), delayed or cumulative. A small study of MI and particulate air pollution in New York reported a positive association of STEMI with PM2.5 exposure in the 1 hour prior to event, with ORs for the 3, 12 and 24 hours prior to event, >1 but non-significant.12 Similarly, a large, though exploratory, study of air pollution and acute MI in Alberta, Canada, reported positive associations of NSTEMI in the 65 and over age group with average NO2 exposures in the 6, 12 and 24 hours prior to event (suggesting elevations in risk of 3.7%, 3.8% and 3.7% per 10 µg/m3, respectively).13
In terms of MINAP, a previous study of the database (2003–2006) by Bhaskaran et al,28 found evidence of positive associations between MI and NO2 and PM10 exposure in the 1–6 hours prior to event (suggesting elevations in risk of 1.1% and 1.2% per 10 µg/m3, respectively) but no evidence of associations with exposures in the 7–12, 13–18, 19–24 or 25–72 hours prior. While, in contrast, the Milojevic et al,11 study found stronger associations of all MI and NSTEMI with NO2 at lags 0–4 days rather than NO2 at lags 0–1 days, our sensitivity analyses with modelled pollutant data (figure 1) suggested stronger associations with the more recent exposure (ie, lag 0 days) than lags 1 or 2 days.
Much has been written and hypothesised as to the various possible mechanisms by which an increase in the level of exposure to certain pollutants may increase the risk of an acute coronary event, including by increasing blood pressure, blood viscosity, promoting an inflammatory response, interfering with heart rhythm or promoting vasoconstriction.3 ,4 Of interest here is how these mechanisms may differentially impact on STEMI and NSTEMI. The total occlusion of a coronary artery characteristic of STEMI is more likely to result from plaque rupture and thrombus formation than the partial occlusion characteristic of NSTEMI.12 ,29 Gardner et al,12 therefore, suggested that whether exposure to a given pollutant (in their case PM2.5) was more important in the aetiology of STEMI than NSTEMI may depend on whether or how that pollutant influences the processes leading to thrombus formation and/or thrombus dissolution.
In common with a previous study of the MINAP database,11 but using a different data set with notably wider geographical coverage, modelled temperature and pollution data modelled by an atmospheric chemistry transport model, we found some evidence of a positive association between exposure to NO2 and the risk of NSTEMI. The appropriate exposure period or lag to consider is still unclear and more immediate effects of other pollutants in the hours rather than the day of and or days prior to event may have been missed. Given our study detected only one positive association, it is possible that this finding is simply an artefact. Confirmation of any such association in databases other than MINAP is therefore required.
This study includes data collected on behalf of the British Cardiovascular Society under the auspices of the National Institute for Cardiovascular Outcomes Research (NICOR). The use of monitoring data from the UK Department for Environment, Food and Rural Affairs Automated Urban and Rural Network (AURN) for air pollutant monitoring (© Crown 2015 copyright Defra via uk-air.defra.gov.uk, licenced under the Open Government Licence (OGL)) is acknowledged. The use of emissions data from the National Atmospheric Emissions Inventory (© Crown 2015 copyright Defra and DECC via naei.defra.gov.uk, licenced under the Open Government Licence (OGL)) is acknowledged. Also licenced under the OGL, the authors acknowledge the use of data on RCGP influenza-like illness (ILI) consultation rates per 100 000 population, England and Wales (Source: the Office for National Statistics and the Royal College of General Practitioners Research and Surveillance Centre).
Contributors BKB conducted the main statistical analysis and took the lead in drafting the paper. RWA, MRH, RMD, BGA and PW contributed to the design and concept of the study. AM pre-analysed the MINAP data and constructed the linked database. CL and MRH assembled the monitor-model comparison data sets and along with BGA and BKB were involved in producing the model performance statistics. IAM, MV and RMD generated the ACTM data and processed it for use in this study. MV is the main developer of the regional EMEP4UK atmospheric chemistry-transport model based on the widely used EMEP European regional model. All authors contributed to the drafting of the paper and the interpretation of results.
Funding This work was performed under the cross UK-research councils Environmental Exposure and Human Health Initiative (EEHI) grants: NE/I007865/1, NE/I007938/1 and NE/I008063/1. Other than to review the study design as part of the grant giving process, the funders did not have any involvement in the collection, analysis and interpretation of data, in the writing of the report or in the decision to submit the article for publication. The EMEP4UK ACTM is also supported by the UK Department for the Environment, Food and Rural Affairs (Defra) and, the NERC Centre for Ecology and Hydrology (CEH).
Competing interests BKB owns shares in Royal Dutch Shell and Scottish and Southern Energy and her spouse has a deferred Shell pension. Her work on this project was funded by a grant from the UK-research councils Environmental Exposure and Human Health Initiative (EEHI). RWA and RMD report grants from the UK-research councils EEHI, during the conduct of the study. MRH, IAM, CL and PW report grants from the Natural Environment Research Council, during the conduct of the study.
Ethics approval Ethics approval was obtained from the National Research Ethics Service (15/SC/0494).
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement No additional data are available.