- C Reactive Protein Coronary Heart Disease Genetics Collaboration (CCGC)
- Correspondence to:
- Accepted 9 November 2010
Objective To use genetic variants as unconfounded proxies of C reactive protein concentration to study its causal role in coronary heart disease.
Design Mendelian randomisation meta-analysis of individual participant data from 47 epidemiological studies in 15 countries.
Participants 194 418 participants, including 46 557 patients with prevalent or incident coronary heart disease. Information was available on four CRP gene tagging single nucleotide polymorphisms (rs3093077, rs1205, rs1130864, rs1800947), concentration of C reactive protein, and levels of other risk factors.
Main outcome measures Risk ratios for coronary heart disease associated with genetically raised C reactive protein versus risk ratios with equivalent differences in C reactive protein concentration itself, adjusted for conventional risk factors and variability in risk factor levels within individuals.
Results CRP variants were each associated with up to 30% per allele difference in concentration of C reactive protein (P<10−34) and were unrelated to other risk factors. Risk ratios for coronary heart disease per additional copy of an allele associated with raised C reactive protein were 0.93 (95% confidence interval 0.87 to 1.00) for rs3093077; 1.00 (0.98 to 1.02) for rs1205; 0.98 (0.96 to 1.00) for rs1130864; and 0.99 (0.94 to 1.03) for rs1800947. In a combined analysis, the risk ratio for coronary heart disease was 1.00 (0.90 to 1.13) per 1 SD higher genetically raised natural log (ln) concentration of C reactive protein. The genetic findings were discordant with the risk ratio observed for coronary heart disease of 1.33 (1.23 to 1.43) per 1 SD higher circulating ln concentration of C reactive protein in prospective studies (P=0.001 for difference).
Conclusion Human genetic data indicate that C reactive protein concentration itself is unlikely to be even a modest causal factor in coronary heart disease.
Persistent inflammation has been implicated in the pathogenesis of coronary heart disease, but causality has not been established for any specific inflammatory mediator.1 C reactive protein, an acute phase protein produced by the liver, is the most extensively studied systemic marker of inflammation.2 Observational epidemiological studies have shown that C reactive protein concentration is log linearly related to risk of subsequent coronary heart disease, though this association depends considerably on levels of conventional risk factors.3 C reactive protein binds to low density lipoproteins4 and is present in atherosclerotic plaques.5 There is, therefore, interest in whether long term average (“usual”) concentration of C reactive protein is itself causally relevant to coronary heart disease.6 7 Randomised trials of interventions specific to C reactive protein, however, have not yet been done in relation to vascular disease outcomes.
In the absence of such trials, focused genetic studies can be used to help judge causality. This approach is known as “mendelian randomisation” because it is based on Mendel’s second law, which states that alleles of different genes assort independently of one another during gamete formation. Consequently, mendelian randomisation analyses are based on Mendel’s observation that inheritance of one trait should be independent of inheritance of other traits.8 For the causal assessment of C reactive protein, a mendelian randomisation analysis should reduce confounding, provided the genetic variants used as proxies for concentration of C reactive protein are unrelated to conventional vascular risk factors and other disease markers. Such studies should also avoid distortions caused by factors occurring later in life (such as the onset of disease) because genetic variants are fixed at conception.8 Hence, mendelian randomisation analyses should confer certain design advantages akin to those in randomised trials. When applied to other risk factors in coronary heart disease, this approach has previously confirmed the causal relevance of low density lipoprotein cholesterol,9 increased the likelihood of causality for Lp(a) lipoprotein,10 and reduced the likelihood of causality for fibrinogen.11
Findings from previous human genetic studies have reduced the likelihood of a major causal role for C reactive protein concentration in coronary heart disease.12 13 14 15 16 17 18 As most known genetic variants related to C reactive protein account for a relatively small portion of the variability in its concentration, however, even larger and more detailed analyses are needed to reliably assess the possibility of any moderate causal role. We report such an analysis based on individual data from 194 418 participants in 47 epidemiological studies. We studied these genetic variants in relation to C reactive protein concentration, other risk factors, and risk of coronary heart disease.
Design and rationale
The study had five inter-related components. Firstly, we selected, a priori, a panel of four single nucleotide polymorphisms that explain 98% of the variation in the CRP gene in populations of European descent. These variants have been shown to influence concentration of C reactive protein without affecting its protein sequence.19 Secondly, we studied whether these polymorphisms are likely to be exclusively associated with C reactive protein concentration by evaluating them in relation to a range of other risk factors. Thirdly, we calculated risk ratios for coronary heart disease with genetically raised concentration using information on these CRP variants. Fourthly, we calculated risk ratios for coronary heart disease with circulating C reactive protein concentration after adjustment for conventional risk factors and variability in risk factors within individuals. Fifthly, we compared risk ratios for coronary heart disease with genetically raised concentration of C reactive protein versus risk ratios seen with equivalent differences in circulating concentrations.
We used detailed information about the composition of the CRP gene to select a parsimonious set of “tagging” single nucleotide polymorphisms (rs3093077, rs1205, rs1130864, and rs1800947) that fully covers the common variations of this gene in populations of European descent (that is, minor allele frequency ≥0.05 and an r2 threshold of ≥0.8).19 20 Data available on 36 further single nucleotide polymorphisms enabled use of proxy variants when principal polymorphisms were not measured. To enhance power, we also studied combinations of alleles inherited together, or “haplotypes” (see table A in appendix 1 on bmj.com).
Details of the C Reactive Protein Coronary Heart Disease Genetics Collaboration (CCGC) have been described previously.20 Tables B-E in appendix 1 on bmj.com provides details of contributing studies, and appendix 2 lists study acronyms. Studies were identified through registry approaches and systematic searches of the literature (see fig A in appendix 3, and appendix 4 on bmj.com). Individual data were supplied on 194 418 participants, including 46 557 with incident or prevalent coronary heart disease, in 47 studies. Studies used different genotyping platforms: 23 used TaqMan assays, three used KASPAR technology (KBioscience), three used restriction fragment length polymorphism, 10 used the ITMAT-Broad-CARe 50K SNP array, and eight used other multiplex methods. Thirty five studies measured C reactive protein with high sensitivity assays, directly or indirectly standardised on the International Reference Standard for C reactive protein immunoassay (WHO 85/506). The outcome was defined as fatal coronary heart disease (based on International Classification of Diseases codes), non-fatal myocardial infarction (using WHO criteria), or coronary stenosis (>50% narrowing of one of more coronary arteries assessed by angiography). All study participants provided written informed consent for use of their DNA for genetic testing.
Appendix 5 on bmj.com provides details of the statistical methods. Levels of C reactive protein and other positively skewed variables were natural log transformed. Principal analyses assumed additive effects (per allele associations), with subsidiary analyses of other genetic models. All analyses of circulating C reactive protein and other risk factors susceptible to reverse association biases were limited to participants without known coronary heart disease at time of blood sampling. We calculated study specific associations of baseline (and repeat) values of C reactive protein concentration and other characteristics with a linear regression model, adjusted for sex and ethnicity and combined across studies using random effects meta-analysis to obtain summary estimates of single nucleotide polymorphism and C reactive protein correlates.21 We calculated risk ratios for coronary heart disease per addition of a prespecified risk allele of each single nucleotide polymorphism using a two stage approach, stratified, where appropriate, by sex, ethnicity, and adjustment for age. We calculated natural log (ln) risk ratios separately within each study before pooling them across studies by random effects meta-analysis.3 For prospective cohort studies, we used stratified Cox proportional hazard regression models to calculate hazard ratios for incident coronary heart disease risk. For retrospective studies (and for prevalent coronary heart disease cases in prospective studies), we used either conditional or unconditional logistic regression as appropriate to calculate odds ratios. Hazard ratios and odds ratios were assumed to approximate the same risk ratios.
To correct risk ratios for coronary heart disease for potential bias caused by measurement error and variability in risk factors within an individual (“regression dilution bias”), we made allowance for regression dilution bias in both C reactive protein concentration and potential confounding factors. Regression dilution ratios were obtained by regressing 93 998 serial measurements (mean interval 2.9 (SD 1.9) years) from 35 023 participants on baseline levels of the relevant characteristic, adjusted for conventional risk factors plus baseline ln C reactive protein concentration and duration of follow-up.22 Risk ratios for coronary heart disease, adjusted for conventional risk factors, were calculated per 1 SD higher usual ln concentration in participants. This difference corresponds to an approximately threefold difference on the original scale of C reactive protein measurement, as the pooled standard deviation for baseline ln C reactive protein is 1.05 (that is, e1.05). To estimate the effect of genetically predicted C reactive protein on coronary heart disease, we used a hierarchical Bayesian meta-analysis framework.23 Using information on baseline ln C reactive protein concentration within each genotype or haplotype group, we obtained a summary estimate from meta-analysis of study specific risk ratios for coronary heart disease per 1 SD higher genetically predicted C reactive protein concentration. We used I2 statistic to assess heterogeneity.24 Diversity at the study level was investigated by grouping studies by relevant characteristics and by meta-regression. Effect modification by several prespecified variables was investigated by formal tests of interaction. Conventional analyses were conducted in Stata v 11.0 and Bayesian analyses in WinBUGS.
Mean age at entry was 59 years (SD 10), 89% of participants were of European descent, and 44% were women. Minor allele frequencies ranged from 0.06 to 0.33 for the principal single nucleotide polymorphisms. For any given polymorphism, minor allele frequencies were similar across studies. Mean baseline concentrations of C reactive protein varied across studies, though standard deviations were broadly similar (see fig B in appendix 3 on bmj.com), with an overall median of 1.7 (5th, 95th centile 0.3, 12.7) mg/L. The regression dilution ratio of ln C reactive protein, adjusted for age and sex, was 0.57 (95% confidence interval 0.51 to 0.64), similar to those observed here for systolic blood pressure (0.51, 0.48 to 0.54), and total cholesterol (0.55, 0.52 to 0.60).
CRP variants, C reactive protein concentration, and levels of other risk factors
Each of these CRP variants was associated with baseline C reactive protein, with per allele differences in C reactive protein concentration of 23% (95% confidence interval 19% to 27%) for rs3093077, 19% (17% to 21%) for rs1205, 14% (12% to 16%) for rs1130864, and 30% (26% to 34%) for rs1800947 (P<0.001; fig 1⇓) (see also fig C in appendix 3 and table F in appendix 1 on bmj.com). These associations were consistent over multiple repeat measures of circulating C reactive protein taken several years apart (see fig D in appendix 3 on bmj.com). With the exception of ethnicity (with which the principal variants were strongly correlated; P<0.001), single genetic variants associated with C reactive protein were not associated with conventional risk factors or other inflammatory markers (fig 1⇓). With haplotype 1 as a reference, CRP haplotypes were associated with differences in ln C reactive protein concentration ranging from 0.1 to 0.34 mg/L and were not associated with conventional risk factors and other characteristics (see figs E-G in appendix 3 and table F in appendix 1 on bmj.com). By contrast, baseline C reactive protein concentration itself was associated with conventional risk factors, inflammatory markers, and other characteristics (see table G in appendix 1 on bmj.com).
Risk ratios for coronary heart disease with CRP genetic variants
Risk ratios for coronary heart disease per addition of a risk (“C reactive protein raising”) allele were 0.93 (0.87 to 1.00) for rs3093077, 1.00 (0.98 to 1.02) for rs1205, 0.98 (0.96 to 1.00) for rs1130864, and 0.99 (0.94 to 1.03) for rs1800947 (fig 2⇓). There was modest heterogeneity in these estimates (I2 values for those risk ratios were 0%, 26%, 0%, and 0%, respectively, see fig C in appendix 3 on bmj.com) with similar findings under a range of circumstances (see table H in appendix 1 on bmj.com). For haplotype analyses, risk ratios for coronary heart disease per copy of haplotype (relative to two copies of haplotype 1) were 1.01 (0.97 to 1.04) with haplotype 2, 0.98 (0.92 to 1.03) with haplotype 3, 0.99 (0.96 to 1.03) with haplotype 4, and 0.98 (0.91 to 1.05) with haplotype 5 (fig 3⇓) (see also fig F in appendix 3 on bmj.com). There was little heterogeneity in these risk ratios (see fig G in appendix 3 on bmj.com). Data were insufficient to investigate effects in different ethnic groups.
Risk ratios for coronary heart disease with usual concentrations of C reactive protein
In analyses restricted to 27 long term prospective studies comprising 124 931 participants and 10 981 incident cases of coronary heart disease, there was an approximately log linear association between C reactive protein concentration and risk of coronary heart disease (see fig H in appendix 3 on bmj.com). In these studies, the risk ratio for coronary heart disease, adjusted for age, sex, and ethnicity only, was 1.49 (1.40 to 1.59) per 1 SD higher “usual” ln C reactive protein concentration (that is, a risk ratio that has made allowance for regression dilution) (fig 4⇓). The risk ratio for coronary heart disease was 1.33 (1.23 to 1.43) after further adjustment for smoking status, history of diabetes mellitus, and usual levels of systolic blood pressure, body mass index (BMI), non-high density lipoprotein cholesterol, high density lipoprotein cholesterol, and triglyceride concentrations. There was little evidence of heterogeneity (see fig I in appendix 3 on bmj.com). Risk ratios with higher C reactive protein concentration were broadly similar under a range of circumstances (see fig J in appendix 3 on bmj.com). Multivariable adjusted risk ratios for coronary heart disease with C reactive protein concentration weakened after further adjustment for fibrinogen or interleukin 6 (see table I in appendix 1 on bmj.com).
Integration of data on CRP variants and C reactive protein concentration
Risk ratios for coronary heart disease were 1.00 (0.90 to 1.13) per 1 SD higher genetically raised C reactive protein as determined by all four principal single nucleotide polymorphisms (163 174 participants, 37 736 cases, 44 studies). This corresponds to a risk ratio of 1.00 (0.97 to 1.02) per 20% lower C reactive protein, which is equivalent to about 0.38 mg/L lower C reactive protein when the population mean is 1.88 mg/L. From information on all common CRP haplotypes in populations of European descent, the risk ratio was 1.00 (0.89 to 1.12; 152 678 participants, 33 589 cases, 39 studies; fig 4⇑) (see also fig K in appendix 3 on bmj.com).
We observed qualitatively similar results to those reported above in analyses that used different genetic models (see fig L in appendix 3 on bmj.com); excluded variants or studies that deviated from Hardy-Weinberg equilibrium or were judged to be insufficiently strong genetic instruments25; used fixed effect meta-analysis models (see figs C, G, I, and K in appendix 3 on bmj.com); omitted the 11 734 participants who reported using cardiovascular drugs (including statins) at baseline; omitted the 18 198 participants in clinical trials; included only people of European descent (see fig M in appendix 3 on bmj.com); compared prospective and retrospective studies (fig K in appendix 3 on bmj.com); and compared larger studies versus smaller studies (available on request).
Implications for disease aetiology
Using the principle of mendelian randomisation, we used CRP genetic variants as proxies to help judge whether usual C reactive protein concentration is causally relevant to coronary heart disease. Our results indicate that genetically raised concentrations of C reactive protein are unrelated to conventional risk factors and risk of coronary heart disease. Given the power of our study, these results suggest that C reactive protein concentration is unlikely to have even a modest causal role in coronary heart disease. Furthermore, we found no material associations between CRP genetic variants and fibrinogen, interleukin 6, and leucocyte count. Because these results suggest that C reactive protein concentration is unlikely to be causally relevant to these inflammatory markers, the findings imply that analyses of C reactive protein concentration and disease outcomes should adjust for these inflammatory markers (that is, treat them as genuine confounders). Indeed, risk ratios for coronary heart disease with circulating C reactive protein were further weakened when we adjusted for such inflammatory markers. This finding further decreases the likelihood that C reactive protein concentration is causally relevant.3 Collectively, therefore, our results suggest that studies seeking to test the inflammation hypothesis in coronary heart disease should examine inflammatory mediators other than C reactive protein.1 26 27
Our analysis builds on previous studies, combining several major advantages.12 13 14 15 16 17 18 Firstly, we carefully selected and assessed specific polymorphisms to capture the full range of common variability in the CRP gene. Variants in the CRP gene related to C reactive protein concentration are more likely to be exclusively associated with C reactive protein concentration than variants in other genes.18 Secondly, our study involved over 90 000 more participants than the previous largest study, substantially enhancing statistical power. Thirdly, the availability of data from individual participants enabled a comprehensive and detailed range of analyses not possible with aggregated data (such as haplotype analyses and instrumental variables analyses23). For example, compared with analysis of single nucleotide polymorphisms, haplotype analysis allows more complete consideration of genetic variation at a locus. Fourthly, we had information on 94 000 serial measurements, enabling correction of risk ratios for regression dilution bias22 and yielding results supporting the idea that CRP variants are related to lifetime variation in C reactive protein concentration. Finally, the generalisability of our results is supported by their consistency in 47 studies in 15 countries.24
For studies that did not directly measure the principal CRP single nucleotide polymorphisms we investigated, we used alternative variants as proxies. There is the possibility of residual confounding by unrecognised effects of genotypes on other risk factors and by adaptation during early life to compensate for genetically raised concentrations of C reactive protein,8 though there is no evidence of their impact in the current context. Further study is needed in non-white people, particularly as the genetic architecture of CRP differs by ethnicity.28
Our results generally agree with those from previous studies that used other methods, such as the null findings of in vivo studies of atherosclerosis that have involved either injection of C reactive protein in different species or transgenic expression of C reactive protein in mice and rabbits.29 30 Our results also agree with well designed in vitro studies, in which uncontaminated preparations of C reactive protein were used, which have generally failed to produce persuasive evidence for pro-atherothrombotic effects of C reactive protein on various cell types.31 32 Irrespective of the causal relevance of C reactive protein itself to coronary heart disease, however, there is considerable evidence that persistent inflammation might contribute to coronary heart disease.1 33 Hence, there is a need to identify specific genetic, biochemical, and environmental determinants responsible for these associations. Anti-inflammatory interventions have not yet been adequately studied to test the inflammation hypothesis in coronary heart disease.1 26 27
Furthermore, as distinct from possible associations of underlying usual C reactive protein concentration with later development of coronary heart disease, there is interest in the possibility that massively raised concentrations of C reactive protein at the time of acute ischaemic events could contribute to severity and outcome of ischaemic lesions.34 Our findings also do not address the separate issue of the value of measurement of circulating C reactive protein in prediction of long term vascular risk.35
What is already known on this topic
Blood concentrations of C reactive protein are strongly and continuously associated with future risk of coronary heart disease, though it is unknown whether this correlation reflects cause and effect
Genetic variants related to C reactive protein can be used as proxies for C reactive protein concentration to help judge causality (“mendelian randomisation analyses”)
Previous studies have been insufficiently powerful and detailed to evaluate the possibility of any moderate causal role for C reactive protein in coronary heart disease
What this study adds
With individual data from almost 200 000 people (including almost 47 000 with coronary heart disease), the current study has shown that genetically raised concentration of C reactive protein is unrelated to conventional risk factors and risk of coronary heart disease
Human genetic data indicate that C reactive protein concentration itself is unlikely to be even a modest causal factor in coronary heart disease
Cite this as: BMJ 2011;342:d548
Frances Wensley, Pei Gao (University of Cambridge); Stephen Burgess (MRC Biostatistics Unit, Cambridge,); Stephen Kaptoge, Emanuele Di Angelantonio (University of Cambridge); Tina Shah (University College London); James C Engert (McGill University, Canada); Robert Clarke (University of Oxford); George Davey-Smith (University of Bristol); Børge G Nordestgaard (University of Copenhagen, Denmark); Danish Saleheen (University of Cambridge); Nilesh J Samani (University of Leicester); Manjinder Sandhu (University of Cambridge); Sonia Anand (McMaster University, Canada); Mark B Pepys (University College London); Liam Smeeth, John Whittaker, Juan Pablo Casas (London School of Hygiene and Tropical Medicine, UK); Simon G Thompson (MRC Biostatistics Unit, Cambridge); Aroon D Hingorani (University College London); John Danesh (University of Cambridge).
C reactive Protein Coronary Heart Disease Genetics Collaboration (CCGC) Investigators
AGES: G Eiriksdottir, T B Harris, L J Launer, V Gudnason; ARIC: A R Folsom, G Andrews, C M Ballantyne; BHF-FHS: N J Samani, A S Hall, P S Braund, A J Balmforth; BRHS: P H Whincup, R Morris; BWHHS: D A Lawlor, G D O Lowe, N Timpson, S Ebrahim; CAPS: Y Ben-Shlomo, G Davey-Smith, N Timpson; CCHS, CGPS, and CIHDS: B G Nordestgaard, A Tybjærg-Hansen, J Zacho; CHAOS: M Brown, M Sandhu, S L Ricketts, S Ashford; CHS: L Lange, A Reiner, M Cushman, R Tracy (see www.chs-nhlbi.org for acknowledgments); CRPHANS: C Wu, J Ge, Y Zou, A Sun; CUDAS and CUPID: J Hung, B McQuillan, P Thompson, J Beilby, N Warrington, L J Palmer; 4D: C Wanner, C Drechsler, M M Hoffmann; EAS: F G R Fowkes, G D O Lowe, I Tzoulaki; ELSA: M Kumari, M Miller, M Marmot; EPICNL: C Onland-Moret, Y T van der Schouw, J M Boer, C Wijmenga; EPICNOR: S L Ricketts, S Ashford, M Sandhu, K-T Khaw; FRAMOFF: R S Vasan, R B Schnabel, J F Yamamoto, E J Benjamin; GERMIFS: H Schunkert, J Erdmann, I R König, C Hengstenberg; GISSI-P: B Chiodini, M G Franzosi, S Pietri, F Gori; HEALTHABC: M Rudock; Y Liu, K Lohman, T B Harris; HIFMECH: S E Humphries, A Hamsten; HIMS: P E Norman, G J Hankey, K Jamrozik, L J Palmer; HPFS: E B Rimm, J K Pai; HVHS: B M Psaty, S R Heckbert, J C Bis; INTERHEART: S Yusuf, S Anand, J C Engert, C Xie; ISIS: R Collins, R Clarke, D Bennett; LOLIPOP: J Kooner, J Chambers, P Elliott; LURIC: W März, M E Kleber, B O Böhm, B R Winkelmann; MALMO: O Melander, G Berglund; MONICAKORA: W Koenig, B Thorand, J Baumert, A Peters; NHS: E B Rimm, J Manson, J K Pai; NPHSII: S E Humphries, J A Cooper, P J Talmud; NSC: P Ladenvall, L Johansson, J-H Jansson, G Hallmans; PENNCATH: M P Reilly, L Qu, M Li, D J Rader; PROCARDIS: H Watkins, R Clarke, J Hopewell; PROMIS: D Saleheen, J Danesh, P Frossard; PROSPER: N Sattar, M Robertson, J Shepherd, E Schaefer; ROTTERDAM: A Hofman, J C M Witteman, I Kardys, A Dehghan; SHEEP: U de Faire, A Bennet, B Gigante, K Leander; SPEED: Y Ben-Shlomo, G Davey-Smith, N Timpson; UCP: B Peters, A H Maitland-van der Zee, A de Boer, O Klungel; WHIOS: A Reiner, J Manson, P Greenland, J Dai, S Liu; WHITEII: M Kumari, E Brunner, M Kivimaki, M Marmot; WOSCOPS: N Sattar, D O’Reilly, I Ford, C J Packard.
We thank C Hennekens for helpful comments.
Contributors: Statistical team: S G Thompson (co-chair), J Whittaker (co-chair), S Burgess, P Gao, S Kaptoge, C Verzilli, and F Wensley (coordinator). Data management: M Walker and S Watson. Coordinating centre: F Wensley (coordinator), S Anand, S Burgess, J P Casas, E Di Angelantonio, J Engert, P Gao, S Kaptoge, MB Pepys, T Shah, L Smeeth, SG Thompson, C Verzilli, M Walker, S Watson, J Whittaker, S Yusuf, A D Hingorani (co-principal investigator), and J Danesh (co-principal investigator). P Gao, S Burgess, S Kaptoge, and E Di Angelantonio contributed equally to this work, as did AD Hingorani and J Danesh. J Danesh is guarantor.
Funding: The coordinating centre is supported by grant SP/08/007 from the British Heart Foundation. Various sources have supported recruitment, follow-up, and laboratory measurements in the contributing studies. Investigators of several of these studies have contributed to a list naming some of these funding sources, which can be found at http://ceu.phpc.cam.ac.uk/research/ccgc/studies/. AD Hingorani acknowledges support from the Rosetrees Foundation and the British Heart Foundation (FS05/125).
Competing interests: All authors have completed the Unified Competing Interest form at www.icmje.org/coi_disclosure.pdf (available on request from the corresponding author) and declare: no support from any organisation for the submitted work; no financial relationships with any organisations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
Ethical approval: This study was approved by the Cambridgeshire ethics review committee, and was conducted independently from its funders.
Data sharing: No additional data available.
This is an open-access article distributed under the terms of the Creative Commons Attribution Non-commercial License, which permits use, distribution, and reproduction in any medium, provided the original work is properly cited, the use is non commercial and is otherwise in compliance with the license. See: http://creativecommons.org/licenses/by-nc/2.0/ and http://creativecommons.org/licenses/by-nc/2.0/legalcode.