Article Text
Abstract
Background The relationship between atherosclerosis and endotypes of myocardial ischaemia with no obstructive coronary artery disease (INOCA) is unclear. We investigated potential associations between cumulative atherosclerotic plaque burden quantified using the Gensini score, novel invasive indices of coronary microvascular function (microvascular resistance reserve (MRR); resistive reserve ratio (RRR)) and related INOCA endotypes.
Methods Coronary angiography and invasive coronary function tests were simultaneously acquired in the CorMicA cohort. A comprehensive physiological assessment was performed using both a thermodilution-based diagnostic guidewire and intracoronary acetylcholine provocation testing. Angiograms were examined for luminal stenosis in each segment of the SYNTAX coronary model. Cumulative plaque burden was quantified using the Gensini score, which incorporated both the number of diseased coronary segments and stenosis severity. Results were compared with indices of microvascular function and INOCA endotypes. Angiographic analyses were performed blind to coronary physiology findings.
Results In 151 participants (median age 61 years; 73.5% female) without flow-limiting coronary artery disease, medical history included 41.7% smoking, 63.6% hypertension and 19.2% diabetes mellitus. The left anterior descending artery underwent diagnostic guidewire testing in 85.4%, and 55.0% of participants had abnormal coronary flow reserve (CFR) and/or Index of Microcirculatory Resistance (IMR). The median Gensini score was 6.0 (IQR 2.5–11.0). CFR (p=0.012), MRR (p=0.026) and RRR (p=0.026), but not IMR (p=0.445), were univariably associated with raised Gensini scores. These significant effects persisted in multivariable models controlling for potential confounders. Considering INOCA endotypes, Gensini scores differed among participants with microvascular angina (MVA) (7.0 (2.5–11.0)), vasospastic angina (VSA) (4.5 (2.0–10.0)), mixed MVA/VSA (9.0 (5.0–11.5)) and non-cardiac symptoms (3.5 (1.5–8.0)); Kruskal-Wallis p=0.030.
Conclusions Reduced CFR, MRR and RRR, and MVA were associated with increased coronary atherosclerotic plaque burden, as evidenced by higher Gensini scores. These novel findings provide a mechanistic link between INOCA and cardiovascular events, reinforcing the importance of antiatherosclerosis therapy in patients with MVA.
- Atherosclerosis
- Coronary Angiography
- Microvascular Angina
- Cardiovascular Diseases
- Angina Pectoris
Data availability statement
Data are available upon reasonable request.
This is an open access article distributed in accordance with the Creative Commons Attribution 4.0 Unported (CC BY 4.0) license, which permits others to copy, redistribute, remix, transform and build upon this work for any purpose, provided the original work is properly cited, a link to the licence is given, and indication of whether changes were made. See: https://creativecommons.org/licenses/by/4.0/.
Statistics from Altmetric.com
WHAT IS ALREADY KNOWN ON THIS TOPIC
Epicardial coronary atherosclerosis of varying severity is prevalent among patients with ischaemia with no obstructive coronary artery disease (INOCA), although the pathophysiological associations between atherosclerosis and microvascular dysfunction remain uncertain.
WHAT THIS STUDY ADDS
In this cross-sectional study of prospectively recruited individuals, clinical endotypes of INOCA were defined by coronary angiography, fractional flow reserve, thermodilution tests of microvascular function and responses to intracoronary infusion of acetylcholine. This was then correlated with the Gensini score, a quantitative measure of coronary plaque burden.
Novel measures of coronary microvascular dysfunction, including reduced microvascular resistance reserve and resistive reserve ratio, and microvascular angina, were associated with higher Gensini scores.
HOW THIS STUDY MIGHT AFFECT RESEARCH, PRACTICE OR POLICY
Our study identifies mechanistic associations between coronary atherosclerosis burden and microvascular dysfunction in INOCA, thereby reinforcing the importance of antiatherosclerosis therapy in this population.
Introduction
Many patients with suspected or proven myocardial ischaemia are found to have no flow-limiting disease in the epicardial circulation (ischaemia with no obstructive coronary artery disease (INOCA)).1 The Coronary Vasomotion Disorders International Study Group (COVADIS) have established criteria for the diagnosis microvascular and vasospastic angina (VSA).2 3 These conditions are prognostically important1 4 5 and, potentially, a therapeutic target.6 Guidelines recommend coronary function testing to aid clinical diagnosis.1 7 In the CorMicA trial,8 patients underwent invasive coronary function testing, to allow endotype classification according to COVADIS criteria2 3: microvascular angina (MVA), VSA, both, or non-cardiac chest pain.
Epicardial coronary atherosclerosis of varying severity is prevalent among patients with INOCA,9 although the pathophysiological associations between atherosclerosis and microvascular dysfunction remain uncertain. Echavarria-Pinto et al9 studied 78 patients with ‘intermediate’ coronary artery disease. Coronary flow reserve (CFR) and Index of Microcirculatory Resistance (IMR) differed according to the distribution of coronary artery disease (focal vs diffuse). However, this cohort was relatively small, 41% of patients had flow-limiting atherosclerosis (fractional flow reserve (FFR) ≤0.80), acetylcholine vasoreactivity testing was not included, and substantive analyses in the subgroup of patients with non-flow-limiting atherosclerosis (FFR >0.80; n=54 arteries) were uncertain. While CFR reflects the combined vasodilator reserve of the epicardial coronary artery and microcirculation, specific associations between atherosclerosis burden and novel microvascular-specific indices are uncertain.
To address these knowledge gaps, we aimed to investigate the relationships between atherosclerosis and microvascular function in a well-defined, prospectively recruited population. The microvascular resistance reserve (MRR)10 11 and resistive reserve ratio (RRR)12 13 are novel indices. Unlike CFR, they are specific measures of microvascular function and independent of epicardial coronary atherosclerosis. We investigated potential associations between cumulative atherosclerotic plaque burden quantified using the Gensini score, novel invasive indices of coronary microvascular function and related INOCA endotypes. We hypothesised MRR and RRR are inversely correlated with atherosclerosis burden in this population.
Methods
Study design
This was a retrospective analysis of the CorMicA trial population. The trial design, 6- and 12-month results have been published.8 14 15 The study is an investigator-led, parallel-group, randomised, controlled trial with a blinded outcome assessment. Participants with no obstructive coronary artery disease (<50% luminal diameter stenosis on angiography, or FFR >0.80) underwent an interventional diagnostic procedure at the time of invasive coronary angiography.
Participants
Participants were recruited from elective, adult referrals to two regional cardiac centres in the west of Scotland (population of 2.5 million). All outpatients aged ≥18 years who were undergoing clinically indicated invasive coronary angiography were eligible. Exclusion criteria were non-coronary indications for angiography (eg, valvular disease) and inability to provide informed consent. Only participants with non-flow limiting coronary artery disease were included in this study.
Interventional diagnostic procedure and indices of microvascular function
The full details of the IDP protocol are available in online supplemental file 1.14 It comprised two components: assessment with a diagnostic guidewire and vasospasm provocation with a pharmacological agent.
Supplemental material
First, a combined pressure/temperature diagnostic guidewire (PressureWire X, Abbott Cardiovascular) was advanced to the distal (>6 cm from ostium) coronary artery. The left anterior descending (LAD) was selected as default unless technical or clinical reasons (eg, tortuous anatomy) advocated choosing a different coronary artery. Intracoronary short-acting nitrate (GTN 200 µg) was administered to attenuate any resting spasm. Using the principles of thermodilution, rapid boluses of 3 mL room temperature normal saline were performed to obtain transit times at rest. Intravenous adenosine infusion (140 µg/kg/min) was commenced to induce systemic hyperaemia and thermodilution injections repeated. Haemodynamic responses were recorded using linked software (CoroFlow V.3.01, CoroVentis) and used to calculate CFR (abnormal <2.0) and IMR (abnormal ≥25).
The MRR and RRR were retrospectively calculated for all subjects, using previously obtained physiology data. The MRR is specific for the microvasculature, independent of autoregulation and myocardial mass.11 It is calculated by dividing the CFR by FFR and correcting for driving pressures. RRR is a marker of the capacity of the coronary microcirculation to reduce resistance in response to hyperaemia. It is the ratio of the resting to hyperaemic resistance index.13
Finally, pharmacological provocation testing for spasm was undertaken with sequential 2 minute intracoronary acetylcholine infusions of incremental concentrations (10−6 M, 10−5 M, 10−4 M), followed by a bolus (100 µg for left or 50 µg for right coronary artery (RCA)).
Diagnosis of clinical endotypes
Each INOCA endotype is characterised by distinct functional or pathophysiological mechanisms. A reduced CFR (<2.0) reflected impaired vasodilator reserve (functional microvascular dysfunction), whilst an elevated IMR (≥25) reflected microvascular dysfunction due to structural remodelling. The vasoactive response to intracoronary infusion of acetylcholine disclosed epicardial coronary and/or microvascular spasm, or a normal response.
According to COVADIS criteria,2 MVA is diagnosed with an abnormal response during either the diagnostic guidewire test (CFR and/or IMR) and/or pharmacological provocation (slow flow phenomenon, ischaemic ECG changes and chest pain). VSA is supported by epicardial coronary spasm (>90% diameter reduction) with chest pain and ECG changes.3 Mixed MVA and VSA is diagnosed if features from both conditions are present. Non-cardiac symptoms are diagnosed if the diagnostic guidewire and pharmacological tests are normal.
Angiographic analysis
The spatial resolution of invasive coronary angiography is approximately 0.3 mm. Clinician visual assessment of the angiogram is the standard approach for quantifying coronary atherosclerotic plaque. While less sensitive than cross-sectional imaging, detection of luminal stenosis is nonetheless specific for atherosclerotic plaque.16 17
The coronary angiograms were analysed (Carestream Vue PACS, V.12.2.2.1025) by two clinicians trained in angiographic interpretation (DTYA and JC). Stenosis severity was determined by visual assessment, from at least two angiographic views at least 30 degrees projection apart, and scored based on the highest luminal stenosis identified in that segment. The assessment was made for all coronary segments sized >2.0 mm diameter, with each segment defined according to the SYNTAX coronary anatomy map.18
Interobserver disagreements were resolved through a joint review of images and/or adjudication by a third clinician (CB). All angiographic assessments and subsequent Gensini scoring were performed while blind to the results of the coronary physiology findings (ie, CFR, IMR, MRR and RRR).
Assessing atherosclerotic plaque distribution and burden: the Gensini score
First, the atherosclerotic plaque distribution was defined according to the aggregated number of SYNTAX coronary segments with luminal stenosis ≥1% evident on the angiogram. Then, the Gensini score was calculated as follows. The Gensini score was first described in 198319 and later summarised by Rampidis et al,20 as a cumulative quantification of both plaque distribution and severity. The Gensini score considers the full spectrum of atherosclerotic stenosis, with increasing points allocated based on stenosis severity (score 0 for no stenosis, 1 for 1–25% stenosis, 2 for 26–50%, 4 for 51–75%, 8 for 76–90%, 16 for 91–99%, 32 for chronic total occlusion).
As all participants in this study had no obstructive coronary artery disease, the maximum unweighted score anticipated per coronary segment was 4 (moderate luminal stenosis with FFR >0.80). Next, a weighting factor was applied to the initial score according to the coronary segment involved (eg, left main stem x5, proximal circumflex x3.5 if left dominant circulation or x2.5 if right dominant, x0.5 for posterior left ventricular branch).20 Finally, the Gensini score was the sum of the weighted stenosis scores for all coronary segments. These scores were investigated for associations with indices of microvascular function (CFR, IMR, MRR, RRR) and compared between the diagnosed clinical endotypes.
Statistical methods
Continuous variables are presented as mean+/-SD or median (IQR), according to their normal or non-normal distribution. Categorical variables are presented as numbers and percentages. Distributions were assessed using normal Q-Q plots. Continuous variables were compared with the t-test, Mann-Whitney U-test, or Kruskal-Wallis H-test, as appropriate. Categorical variables were compared with Fisher’s exact test. Correlation between indices of microvascular function was assessed using Pearson’s or Spearman’s coefficient, according to their normal or non-normal distribution. Univariable and multivariable linear regression were conducted to identify predictors of the Gensini score. Covariates were chosen on clinical grounds and included in all models regardless of statistical significance. Modelling assumptions were checked by visual inspection of residuals, using Q-Q plots. Models were checked for multicollinearity by examination of variance inflation factors (acceptable if VIF <4). A p value of <0.05 was considered significant. Statistical analyses were performed using SPSS software V.29.0.1.0 (IBM, Armonk, New York, USA).
Patient and public involvement
It was not possible to involve patients or the public in the design, conduct, reporting, or dissemination plans of our research.
Results
Between 25 November 2016 and 12 November 2017, 151 participants with symptoms but no obstructive coronary artery disease underwent the IDP according to the CorMicA protocol.
Baseline demographics
Clinical characteristics are summarised in table 1. The median age was 61 years (IQR 53–68), and 73.5% (n=111) were female. The calculated ASSIGN score for the 10-year likelihood of a cardiovascular event was similar between coronary endotype groups.
Global atherosclerotic plaque burden
The median number of SYNTAX coronary segments with ≥1% stenosis was 4 (IQR 2–7). The median number of proximal coronary segments (defined as: left main stem, proximal LAD, proximal circumflex or proximal RCA) with ≥1% luminal stenosis was 1 (IQR 1–3).
In this population without flow-limiting atherosclerosis, the median Gensini score was 6.0 (IQR 2.5–11.0). When divided according to coronary artery (excluding the left main stem), the median Gensini scores were 3.5 (IQR 1.5–5.0) in the LAD, 2.0 (IQR 0.0–3.5) in the left circumflex and 1.0 (IQR 0.0–2.0) in the RCA.
Atherosclerotic plaque burden according to indices of microvascular function
The LAD underwent diagnostic guidewire testing in 85.4% (n=129), the circumflex in 5.3% (n=8) and RCA in 9.3% (n=14).
Indices of microvascular function
The invasive physiology results are summarised in table 2. In the absence of flow-limiting epicardial coronary disease, the CFR and IMR reflect vasoactive function, flow reserve and resistance of the microvascular bed. The CFR, IMR, MRR and RRR (all p<0.001), but not FFR (p=0.207), differed according to INOCA endotypes, reflecting the diagnostic criteria used. However, the CFR, MRR and RRR were higher and IMR lower in patients with isolated VSA, which was not significantly different from patients with non-cardiac symptoms (all p>0.22).
There were strong, positive monotonic correlations between CFR and MRR (rs=0.945, p<0.001), CFR and RRR (rs=0.940, p<0.001) and MRR and RRR (rs=0.985, p<0.001). Conversely, the relationships between CFR and IMR (p=0.364), MRR and IMR (p=0.134), and RRR and IMR (p=0.274) were not significant.
Plaque distribution and Gensini score
Abnormal CFR <2.0 (n=45; p=0.006), MRR ≤3.0 (n=70; p=0.028) and RRR <3.5 (n=85; p=0.044), but not IMR ≥25 (n=52; p=0.477) were associated with more diffuse atherosclerotic plaque distribution, defined by the aggregated number of SYNTAX coronary segments with luminal stenosis ≥1%.
Similarly, the univariable analysis revealed that abnormal CFR <2.0 (p=0.002), MRR ≤3.0 (p=0.025) and RRR <3.5 (p=0.043), but not IMR ≥25 (p=0.445) were associated with greater cumulative plaque burden, as quantified by the Gensini score.
Predictors of raised Gensini score
The factors associated with a raised Gensini score are summarised in table 3. Age, smoking history, hypertension, previous myocardial infarction and measured FFR (all p<0.04) were significant determinants of the Gensini score on univariable regression analysis. Abnormal CFR (p=0.012), MRR (p=0.026) and RRR (p=0.026), but not IMR, were each independently associated with higher Gensini scores. These effects persisted in multivariable models controlling for the potential confounders above (table 3). There was no evidence of multicollinearity in any analyses (all VIFs <1.2) and residual distributions were approximately normal.
Atherosclerotic plaque burden according to INOCA endotype
In total, 55.0% (n=83) of individuals had abnormal CFR and/or IMR. When considering acetylcholine provocation results, 72.2% (n=109/151) met the criteria for MVA and 37.1% (n=56/151) had features of epicardial VSA. The final diagnosis was MVA in 51.7% (n=78), VSA in 16.6% (n=25), mixed microvascular and VSA in 20.5% (n=31) and reclassification to non-cardiac symptoms in 11.3% (n=17).
Participants with coronary microvascular dysfunction as defined by an abnormal CFR <2.0 and/or IMR ≥25 (n=83) had a higher number of diseased coronary segments versus those with normal CFR and IMR (median 5 (IQR 2–7) vs 41–6 ; p=0.039), table 4. Accordingly, the cumulative atherosclerotic plaque burden as assessed by the Gensini score was also higher in participants with versus without microvascular dysfunction by the same criteria (median 7.5 (IQR 3.5–11.5) vs 5.3 (2.5–10.0); p=0.028).
There was a significant difference in the total number of diseased segments when comparing participants with different INOCA endotypes: MVA (median four segments (IQR 2–7)), VSA (4 (1–6)), mixed microvascular and VSA (6 (3–8)) and non-cardiac symptoms (3 (1–5)); p=0.032 (table 4).
Cumulative plaque burden as reflected by the Gensini score also differed by INOCA endotype: MVA (median Gensini 7.0 (2.5–11.0)), VSA (4.5 (2.0–10.0)), mixed microvascular and VSA (9.0 (5.0–11.5)) and non-cardiac symptoms (3.5 (2.0–8.0)); Kruskal-Wallis p=0.030 (figure 1).
Discussion
We evaluated coronary atherosclerotic plaque burden and microvascular function simultaneously in a prospectively recruited cohort of participants with suspected INOCA. We identified novel associations between indices of microvascular function (CFR, MRR and RRR) and the global burden of coronary atherosclerosis, as quantified by the Gensini score. Our findings address a knowledge gap in INOCA; individuals with coronary microvascular dysfunction, particularly those with dual pathologies (mixed MVA and VSA) have substantially greater plaque burden compared with those with non-cardiac symptoms.
Our study has several novel findings. First, we provide a systematic global plaque assessment of all coronary segments in the SYNTAX classification using the Gensini score, a validated measure of coronary atherosclerosis burden. Second, we explored the associations between plaque burden and novel markers of microvascular function (MRR, RRR). Third, our cohort participants underwent comprehensive assessment with both bolus thermodilution-based indices, as well as intracoronary acetylcholine provocation testing. This allowed for stratification of participants by specific mechanisms of vasomotor dysfunction and differentiates our study from prior research.9 21
AlBadri et al21 previously reported increased coronary plaque burden assessed using single coronary artery intravascular ultrasound in patients with microvascular dysfunction, as identified using a Doppler diagnostic guidewire. While our study agrees with their finding of increased plaque burden in patients with versus without coronary microvascular dysfunction, they found an association with the Doppler-derived hyperaemic resistance ratio (HMR), which is conceptually similar, but distinct from IMR. The differences in our results may be explained by the differing methodology for plaque assessment (global analysis from angiography vs single-vessel intravascular imaging) and measured indices (IMR vs HMR). Other published works demonstrated strong links between reduced CFR and MRR with adverse cardiovascular outcomes.1 4 5 22 However, the underlying mechanisms remain incompletely understood.23
We describe correlations between CFR, MRR, RRR and the Gensini score. The higher Gensini scores in participants with MVA are driven by: a greater number of diseased coronary segments, reflecting more diffuse disease; more severe luminal stenosis, indicating higher plaque burden per segment; and more frequent involvement of proximal coronary arteries. Our findings provide a mechanistic link between these indices and the higher rates of acute coronary syndromes seen with INOCA, via increased coronary plaque burden.23
In the absence of flow-limiting epicardial coronary artery disease, an abnormal CFR reflects abnormalities in the coronary microcirculation,24 supporting a diagnosis of MVA. Our findings provide objective, quantitative evidence of an association between coronary microvascular dysfunction and diffuse atherosclerosis in a population with INOCA.
Importantly, reduced CFR, MRR and RRR remained independently associated with higher Gensini scores after controlling for potential confounders. Our multivariable modelling controlled for FFR, the reference physiological index for assessment of epicardial disease. Furthermore, while CFR is influenced by epicardial stenosis, indices specific to microvascular function like MRR and RRR were also independently associated with a raised Gensini score.
Strong correlations were observed between CFR, MRR and RRR. However, IMR was not associated with cumulative plaque burden. Since IMR may be a surrogate measure for microvascular structural remodelling in chronic coronary syndromes, the lack of association between IMR and atherosclerosis may reflect distinct vascular pathophysiology.
Clinical implications
Our findings have implications for patient care. Treatments that stabilise plaques and potentiate regression, importantly statins, may reduce cardiovascular risk in INOCA, particularly in patients with impaired CFR and/or MRR. MVA associated with cardiovascular risk factors that mediate epicardial atherosclerosis.25 26 Similarly, VSA is linked to smoking history.27 These associations support the hypothesis that microvascular dysfunction and epicardial coronary artery disease have a shared pathophysiological basis. Ultimately, our findings reinforce the importance of antiatherosclerosis therapies in INOCA. The potential prognostic benefits of intensive control of cardiovascular risk factors are currently being assessed in the WARRIOR trial.28
Study limitations
Atherosclerotic changes may only be evident on invasive coronary angiography after plaque burden within the vessel wall exceeds >40% by cross-sectional area.29 This in turn may reduce the sensitivity of the Gensini score. Invasive angiography remains the reference standard for assessing coronary luminal stenoses. The angiographic views (>30 degrees projection apart, ideally orthogonal) reflect standard clinical practice, still allowing for assessment of eccentric plaque.
Our study did not employ intravascular imaging. While standard angiographic analysis is less sensitive for detecting atherosclerotic plaque, it remains highly specific30 and has identified novel associations. Future studies could incorporate multi-vessel assessment with intracoronary imaging (intravascular ultrasound or optical coherence tomography). CT coronary angiography could be an alternative, but its spatial resolution limits the analysis of distal coronary branches. Furthermore, CT does not provide a measure of vasomotor function. In our study, coronary plaque burden and microvascular function assessments were acquired simultaneously.
The participants in the CorMicA study were selected, undergoing clinically indicated invasive investigation.14 The true incidence of microvascular dysfunction and plaque burden in the general population remain ill defined. Our study had a cross-sectional design and future cohort studies should clarify longer-term prognostic implications of atherosclerosis in INOCA.
Conclusions
We provided novel insights into relations between global atherosclerotic plaque burden defined in all SYNTAX coronary segments using the Gensini score, and INOCA endotypes as stratified according to comprehensive invasive coronary function testing.
Abnormal CFR, MRR and RRR, and MVA were associated with increased atherosclerotic plaque burden, evidenced by higher Gensini scores. These novel findings provide a mechanistic link between INOCA and cardiovascular events, reinforcing the importance of anti-atherosclerosis therapy in patients with MVA.
Data availability statement
Data are available upon reasonable request.
Ethics statements
Patient consent for publication
Ethics approval
The West of Scotland Research Ethics Committee approved the CorMicA study (REC 1 reference 16/WS/0192). Participants gave informed consent to participate in the study before taking part.
Acknowledgments
The authors acknowledge the patients, the staff members and the British Heart Foundation, who supported the CorMicA study.
References
Supplementary materials
Supplementary Data
This web only file has been produced by the BMJ Publishing Group from an electronic file supplied by the author(s) and has not been edited for content.
Footnotes
X @DanielTYAng, @_RobSykes
Contributors Drs DTYA, JC, TJF and CB were responsible for study design, data collection and analysis. Dr DTYA was responsible for manuscript preparation and revision. Dr AM provided statistical support for data analysis. All other authors assisted with data collection, interpretation of data, revision of the manuscript for important intellectual content and final approval of the version to be published. All authors agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. Drs DTYA and CB act as guarantors of study data integrity. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.
Funding The CorMicA trial was an investigator-initiated clinical trial that was funded by the British Heart Foundation (PG/17/2532884; RE/13/5/30177; RE/18/6/34217).
Competing interests Dr TF received consulting fees from BioExcel; honoraria from Abbott, Boehringer, Novartis. Dr SW received honoraria from Abbott, AstraZeneca, Biosensors, Boston Scientific, Daiichi Sankyo, GE Healthcare and ShockWave Medical. Dr KR received honoraria from AstraZeneca and Abbott. Dr MM received consulting fees from Abbott, Boston Scientific and ShockWave Medical; honoraria from International Medical Device Solutions and Medtronic. Dr KGO received honoraria from Abbott, Biosensors International and Boston Scientific; full-time employee of Biosensors International since May 2020. Dr CB receives research funding from the British Heart Foundation grant (RE/18/6134217, PG/19/28/34310), Chief Scientist Office, EPSRC (EP/R511705/1, EP/S030875/1) and Medical Research Council (MR/S018905/1). None of the declared interests regards the submitted work. All other authors have nothing to disclose.
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.
Provenance and peer review Not commissioned; externally peer reviewed.
Supplemental material This content has been supplied by the author(s). It has not been vetted by BMJ Publishing Group Limited (BMJ) and may not have been peer-reviewed. Any opinions or recommendations discussed are solely those of the author(s) and are not endorsed by BMJ. BMJ disclaims all liability and responsibility arising from any reliance placed on the content. Where the content includes any translated material, BMJ does not warrant the accuracy and reliability of the translations (including but not limited to local regulations, clinical guidelines, terminology, drug names and drug dosages), and is not responsible for any error and/or omissions arising from translation and adaptation or otherwise.