Article Text

Download PDFPDF

Aetiological differences in neuroanatomy of the vegetative state: insights from diffusion tensor imaging and functional implications
  1. Virginia F J Newcombe1,2,
  2. Guy B Williams2,
  3. Daniel Scoffings3,
  4. Justin Cross3,
  5. T Adrian Carpenter2,
  6. John D Pickard2,4,
  7. David K Menon1,2
  1. 1University Division of Anaesthesia, University of Cambridge, Cambridge, UK
  2. 2Wolfson Brain Imaging Centre, Department of Clinical Neurosciences, University of Cambridge, Cambridge, UK
  3. 3Department of Radiology, University of Cambridge, Addenbrooke's Hospital, Cambridge, UK
  4. 4Academic Neurosurgery Unit, Department of Clinical Neurosciences, University of Cambridge, Cambridge, UK
  1. Correspondence to Dr Virginia Newcombe, University Division of Anaesthesia, University of Cambridge, Box 93, Addenbrooke's Hospital, Hills Road, Cambridge CB2 2QQ, UK; vfjn2{at}


Background An improved in vivo understanding of variations in neuropathology in the vegetative state (VS) may aid diagnosis, improve prognostication and help refine the selection of patients for particular treatment regimes. The authors have used diffusion tensor imaging (DTI) to characterise the extent and location of white matter loss in VS secondary to traumatic brain injury (TBI) and ischaemic–hypoxic injury.

Methods Twelve patients with VS (seven TBI, five ischaemic/hypoxic injuries) underwent MRI including DTI at a minimum of 3 months postinjury. Mean apparent diffusion coefficient, fractional anisotropy and eigenvalues were obtained for whole-brain grey and white matter, the pons, thalamus, ventral midbrain, dorsal midbrain and the corpus callosum. DTI measures of supratentorial damage were compared with a summed measure from the JFK modified Coma Recovery Scale (CRS-R) and with a three-point scale of functional magnetic resonance imaging (fMRI) response to an auditory paradigm to assess whether residual integrity of supratentorial white matter connectivity correlated with cortical processing.

Results Conventional radiological approaches did not detect lesions in regions where quantitative DTI demonstrated abnormalities. There was evidence of marked, broadly similar, abnormalities in the supratentorial grey- and white-matter compartments from both aetiologies. In contrast, discordant findings were found in the infratentorial compartment, with DTI abnormalities in the brainstem confined to the TBI group. Supratentorial DTI abnormalities correlated with the CRS-R as well as responses to an fMRI paradigm that detected convert cognitive processing.

Conclusions DTI may help to characterise differences in patients in VS. These findings may have implications for response to therapies, and should be taken into account in trials of interventions aimed at arousal in VS.

  • Vegetative state
  • traumatic brain injury
  • hypoxic brain injury
  • diffusion tensor imaging
  • MRI
  • neuroradiology
  • persistent vegetative state

Statistics from

Request Permissions

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.


The vegetative state (VS) encompasses a spectrum of patients who have emerged from coma to be awake, but show no indication of awareness. The aetiology of these clinical outcomes is varied and includes processes as disparate as traumatic brain injury (TBI), ischaemic/hypoxic injuries (IHI), intracerebral haemorrhage and infection. Modern imaging tools continue to provide more data on the neuroanatomical substrate of these conditions in vivo, and such information may be important to gain further insights into the disease process.

The commonest aetiologies of VS are TBI and IHI. While it is widely accepted that both these conditions arise from either damage to the brainstem/midbrain structures and/or diffuse axonal/cortical damage, most data come from histopathological studies. Such studies show that thalamic damage is common in the vegetative state secondary to both IHI (neuronal loss secondary to hypoxia) and TBI (retrograde thalamic degeneration and/or neuronal loss secondary to hypoxia).1 2 While brainstem lesions have been recognised in IHI, these appear to affect predominantly brainstem grey matter nuclei, whereas more extensive white-matter damage is found in TBI (see Kinney and Samuels3 for an early review of histopathological studies). In TBI patients, focal lesions in the corpus callosum (diffuse axonal injuries (DAI) grade 2), or combined lesions in the corpus callosum and the dorsolateral brainstem (DAI grade 3), occurred in 71% of a series of 25 patients.1 In the same series, the majority of patients with IHI exhibited diffuse damage in the neocortex and hippocampus.1 Some of these patients also exhibited minor abnormalities in the brainstem, largely confined to the substantia nigra or other motor nuclei. These reports, using classical staining techniques, did not find extensive brainstem injuries for either aetiology. However, β-amyloid-precurser protein (β-APP), a stain more sensitive to axonal injury,4 has been used in another series.5 This study, which focused on neuropathology in the pons, found axonal injury present in 80–100% of patients who survived more than 3 h after TBI. Axonal injury was also found in patients post-IHI, although it was noted that the intensity of β-APP accumulation was less than in cases of traumatic injury.5 These studies therefore suggest that classical neuropathological approaches may underestimate the extent and distribution of axonal injury, and that there may be subtle (but potentially important) differences in VS arising from different aetiologies. The ability to characterise ante mortem correlates of these pathological findings using modern imaging techniques has distinct advantages. A better understanding of the variations in neuropathology between patients and aetiologies could result in greater diagnostic precision and improved prognostication, and refine the selection and stratification of patients entered into clinical trial, in both the acute and chronic/rehabilitative stages.

Diffusion-weighted imaging is ideally placed to provide such insights into both the macroscopic and the microstructural damage found in VS. This technique characterises the natural displacements of water molecules. The diffusion tensor can be used to represent the magnitude of water diffusion (apparent diffusion coefficient, ADC), how directional the diffusion is (anisotropy) and the orientation of that direction (eigenvectors/eigenvalues). The diffusion process may be affected by various tissue structures and environments, including pathological processes such as demyelination, axonal truncation or microglia activation. In this study, we used diffusion tensor imaging to characterise and differentiate the extent and location of white-matter loss in patients with disorders of consciousness secondary to traumatic brain injury and ischaemic–hypoxic injury.


These studies were conducted on a bank of data from patients with VS, which have been used for a number of behavioural and imaging analyses. The data in this manuscript are based on 12 (seven TBI, five IHI), who underwent MR imaging at a minimum of 3 months postinjury using a 3 T Siemens Magnetom Total Imaging Matrix (TIM) Trio system. Patients were characterised as VS based on a full assessment using the JFK modified Coma Recovery Scale (CRS-R).6 7 Patients included in this analysis were selected primarily on the basis of DTI data quality, and no dataset was excluded once initial DTI processing revealed data of adequate quality. Images were usually excluded on the basis of excessive movements. In order to avoid bias, all such exclusions preceded placement of regions of interest and measurement of DTI variables. Informed assent for acquisition of imaging data for research was obtained from the next of kin in all cases. Ethical approval was obtained from the Local Research Ethics Committee. Thirty-two age-matched controls (healthy volunteers) underwent an identical imaging protocol which included a 3D T1 weighted structural sequence, magnetisation prepared rapid gradient echo (MPRAGE), functional magnetic resonance imaging (fMRI) and diffusion tensor imaging. The large size compared with the patient population was due to two main reasons: first, the patient group is a difficult group to image successfully, and as DTI parameters in controls are well known to have great variability, it was decided to make the group large to account for this. The DTI parameters were as follows; 12 non-collinear directions, 5 b values ranging from 338 to 1588 s/mm2, 5 b=0 images, field of view 192 mm×192 mm, acquisition matrix 96×96, 63 axial slices, 2 mm slice thickness, TR=8300 ms, TE=98 ms. All images were inspected to ensure that movement during the scan was within acceptable limits (<4 mm). Patients were also scanned with conventional MR sequences (proton density and T2-weighted spin echo (PD-T2), gradient echo (GE), fluid attenuated inversion recovery (FLAIR)) and these were inspected by two neuroradiologists (JC and DS). The presence of lesions that may confound the study was noted (table 1). None of the patients had ventricular shunts, cranioplasty plates or any other metallic implants that would cause artefacts.

Table 1

Categorisation of functional magnetic resonance imaging (fMRI) responses to auditory stimuli8 9

Regions of interest (ROIs) were manually drawn using Analyze 7.0 (Bioimaging Imaging Resource, Rochester, Minnesota, USA)10 in MNI125 space using Colin2711 as a high-resolution, high-signal-to-noise template and included the pons, ventral midbrain, dorsal midbrain, thalamus and corpus callosum (genu and splenium) (figure 1). A supratentorial whole-brain central white matter (WBCWM) ROI was created by subtracting the cortical grey matter, the transitional zone and infratentorial WM from a whole-brain white-matter mask segmented using FMRIB's Automated Segmentation Tool (FAST).12 The segmentation also produced a whole-brain grey matter mask (WBGM). These ROIs were transformed into the normalised space (see below for methods). All coregistered images were visually inspected to ensure that ROIs corresponded to the regions specified. Manual adjustment was performed using Analyze 7.010 if they did not. Care was taken to ensure that the ROIs did not contain any lesions, as blood breakdown products such as haemosiderin cause signal drop-out with DTI sequences. Thus, all ROIs contained only tissue that appeared ‘normal’ on conventional MR. Any misclassified tissue in the WBCWM and WMGM masks was manually removed using Analyze 7.0.10 The mean ADC, fractional anisotropy (FA) and eigenvalues for the different ROIs were calculated. Such an ROI-based approached was chosen over voxel-based methods including tract-based spatial statistics,13 as the often marked atrophy of the brains which can cause coregistration difficulties meant that it was important to be able to manually adjust the regions to ensure like was compared with like.

Figure 1

Examples of regions of interest used. From left to right, whole-brain grey matter, supratentorial white matter, the pons (left and right); midbrain (anterior and posterior); and thalamus, anterior corpus callosum, and posterior corpus callosum respectively.

The DTI data underwent eddy current correction, and FA, ADC and eigenvalue maps were created using the FMRIB's Diffusion Toolbox (FDT)14 To aid coregistration, the skull and extracranial soft tissue were stripped from the MPRAGE images using the Brain Extraction Tool.15 The diffusion-weighted data were normalised using a two-step approach. First, all patient and control MPRAGE images were coregistered to the MNI152 template using the vtkCISG normalised mutual information algorithm.16 The b=0 image was subsequently coregistered to the subject's own MPRAGE. The transformation matrix normalising the MPRAGE image was then applied to the b=0 image.

DTIquery17 was used in the subject's native space to create whole-brain tractography using the FACT algorithm,18 a variable-step streamline tracking method. Seed points for tracts were placed evenly throughout the brain as a grid with spacing of 4 mm. The tractography parameters were: a step-size of 2 mm, deflection angle of 45° and FA termination threshold of 0.15. Using the supratentorial WBCWM mask, to reduce errors due to noise, the number of ‘tracts’ greater than 3 cm were counted to create a surrogate marker of the burden of white matter injury.

In order to assess the clinical relevance of the DTI parameters, we compared them with a summary measure of functional status—the summed highest score derived from the JFK CRS-R.6 7 We also wished to explore whether the residual integrity of supratentorial white matter connectivity correlated with cortical processing as demonstrated by fMRI, independent of the ability to demonstrate behavioural outputs. We therefore compared DTI measures of supratentorial damage with a four-point scale of fMRI responses with an auditory fMRI paradigm based on a hierarchical pattern of speech processing described in Rodd et al,19 and previously described in MCS and VS by Coleman et al8 9(see table 1) with level 0 being no cortical fMRI response to auditory stimuli; level 1, activations characteristic of response to sound only; level 2, including activations characteristic of response to speech; and level 3 retrieval of semantic information. Statistical analyses were conducted using SPSS (SPSS 14.0; SPSS, Chicago, Illinois), and graphs were produced using StatView. Due to the small sample size in the HIH group all analyses, used non-parametric statistics. Mann–Whitney U (MHU) was used for unpaired tests, Spearman ρ was used for correlations, and the Jonckheere–Terpstra test was used to test for trend. Statistical significance was accepted if the p values were <0.05, corrected for multiple comparisons between the three groups (ie, p<0.017 was considered significant).


Five patients (four female, median age 41 (range 33–49) years) with ischaemic brain injuries were compared with seven TBI (two female, median age 39 (range 21–67) years) (table 2). All were classified as VS according to the CRS-R. These patients were compared with 32 controls (seven female, median age 25.2 (range 25–70) years). No significant differences were found between groups for age (p=0.536), or in the time to scan from ictus in the patient populations (IHI median 325 (range 209–1518) days, TBI median 198 (range 105–681) days, p=0.172). Table 3 shows the radiological findings on structural MRI sequences (MPRAGE, PDT2, FLAIR and gradient echo sequences) in the patient group. It is important to note that DTI was more sensitive than conventional radiology—where quantitative DTI measures were calculated, these consistently demonstrated abnormalities in regions where no lesions were detectable by experienced neuroradiologists (tables 3, 4).

Table 2

Demographic and clinical characteristics of the study subjects

Table 3

Conventional imaging characteristics of the patients with percentage of patients showing lesions in the sites shown

Table 4

Conventional imaging characteristics of the patients and abnormalities detected using diffusion tensor imaging

Figure 2 shows whole-brain tractography from a representative control subject, TBI patient and IHI patient. These qualitative appearances were underpinned by quantitative measures of the DTI parameters FA and ADC as well as ‘track counts’ in both supratentorial and infratentorial structures derived from the source DTI images (table 5). It is important to note that the lack of brainstem visualisation was not due to a lack of MR coverage of this area, but the consequence of a decrease in FA, which meant that the tracts could not be traced using standard tracking parameters. ADC and FA were compared within each ROI for males and females in the control group ADC, and FA values for left and right sides were compared for all ROIs in both the control and patient groups separately. No significant differences were found for any of these comparisons, and for subsequent analyses, data from right and left sides were combined.

Figure 2

Fluid-attenuated inversion recovery (FLAIR) and ‘global’ tractography with representative examples from the control, traumatic brain injury (TBI) and ischaemic–hypoxic injury (IHI) groups providing a qualitative image of white matter integrity. This approach is illustrative, and was not used to quantify diffusion tensor imaging (DTI) parameters. It can be clearly seen that the amount of ‘paths’ able to be traced is far less in the patients than the controls. Despite relatively unremarkable appearances on conventional structural imaging (see TBI patient above) the DTI images showed consistent and striking loss of white matter. Also note the preservation of brainstem white matter tracts in IHI when compared with TBI (see arrows; the brainstem in the control subject is visualised despite the presence of a more substantial overlapping temporal lobe). For ease of visualisation, tracts with lengths less than 3 cm are not shown, and the coronal images are overlaid on fractional anisotropy maps.

Table 5

Comparison of patients with vegetative state secondary to hypoxic brain injury to those caused by traumatic brain injury and controls

In both patient groups, FA was significantly decreased and ADC increased compared with controls in the whole-brain white matter, anterior corpus callosum and posterior corpus callosum (figures 3, 4). Both patient groups had higher ADCs than controls in the WBGM and thalamic ROIs. DTI measures were broadly similar for the two aetiologies in these supratentorial ROIs. However, the TBI patients showed significantly lower FA values than both the control group and the IHI in the brainstem regions (dorsal and ventral midbrain and pons) which corresponded to an increase in ADC. Most of the changes in FA and ADC were driven by increases in radial diffusivity (table 4), although less consistent increases in axial diffusivity were noted in some brain regions, particularly in the TBI group.

Figure 3

Fractional anisotropy values for the regions of interest studied. The central lines in the boxes denote the median values, the upper and lower edges the 75th and 25th percentiles, the error bars the 90th and 10th percentiles, and the closed circles the data outside these percentiles. *p<0.017; ** p<0.01; ***p<0.001. IHI, ischaemic/hypoxic injuries; NS, non-significant; TBI, traumatic brain injuries; WBCWM, whole-brain central white matter; WBGM, whole-brain grey matter.

Figure 4

Apparent diffusion coefficient (ADC) values for the regions of interest studied. The central lines in the boxes denote the median values, the upper and lower edges the 75th and 25th percentiles, the error bars the 90th and 10th percentiles, and the closed circles the data outside these percentiles. *p<0.017; **p<0.01; ***p<0.001. IHI, ischaemic/hypoxic injuries; NS, non-significant; TBI, traumatic brain injuries; WBCWM, whole-brain central white matter; WBGM, whole-brain grey matter.

In order to determine whether our measures of supratentorial white-matter damage related to preserved cortical processing, we correlated DTI parameters in the WBCWM and WBGM with clinical parameters and to fMRI responses to an auditory fMRI paradigm, graded on a three-point scale. We found significant trends of decreasing ADC in the WBGM and increasing FA in WBCWM for all patients with a higher fMRI language response (WBGM vs. ADC: J-T statistic=−2.6, p=0.009; WBCWM vs FA: J-T statistic=2.1, p=0.033). A higher CRS-R score had a positive correlation with central white matter FA (0.756, p=0.004) and inverse correlations with central ADC (−0.871, p<0.001) (figure 5). These relationships were not different for the two aetiological categories. No significant correlations were found with age at scan or time from injury to scan.

Figure 5

Relationship of JFK modified Coma Recovery Scale (CRS-R) to the burden of supratentorial central whole-brain white matter damage (fractional anisotropy; r=0.756 (p=0.004), apparent diffusion coefficient; r=−0.871 (p<0.001)). The green line represents the IQR for the controls. IHI, ischaemic/hypoxic injuries; TBI, traumatic brain injuries.


We provide the first DTI comparison of the two most common etiologies of VS, namely ischaemic hypoxic brain injury and traumatic brain injury. Despite the sensitivity of conventional MRI in this setting,20 we confirm that quantitative DTI demonstrates abnormalities missed by conventional radiological assessment. Both TBI and IHI showed evidence of marked, and concordant, abnormalities in the supratentorial GM and WM compartments. These correlated with fMRI responses to a paradigm that detected covert cortical cognitive processing,19 and to functional status as defined by the summed CRS-R. In contrast, the two aetiologies exhibited markedly discordant findings in the infratentorial compartment, where DTI abnormalities in the brainstem were confined to the TBI group. An eigenvalue analysis suggested that FA changes were consistently associated with increased radial diffusivity, with less consistent changes in axial diffusivity.

Few studies have used DTI to assess patients in chronic VS. A cross-sectional study of 52 patients with impaired consciousness (19 VS and 33 severe disability) secondary to TBI correlated motor-evoked potentials with corticospinal tract damage as found with DWI.21 One case study of a TBI patient who recovered expressive language after 19 years in the minimally conscious state (MCS) raised the intriguing possibility of DTI being able to detect neuronal plasticity,22 since this apparent functional improvement was associated with increases in anisotropy. However, the complexity of microstructural changes and the confounders caused by crossing fibres in data-processing algorithms makes it difficult to interpret these findings with confidence without larger numbers and imaging at more time points. In addition, most previous studies have classified patients based on functional status (ie, VS vs MCS) and have tended to focus on a single aetiological group. We believe that an analysis that takes account of the concordant and discordant findings in different aetiologies may provide further insights into the anatomical substrates and pathophysiology of these conditions.

The basis of VS is complex in both aetiologies. Structural and functional imaging studies, including FDG PET, have illustrated widespread damage to areas that are part of cortical networks implicated in consciousness and its constituent processes, including memory, attention and language.23 Marked structural loss has been found using both conventional MRI sequences24 and voxel-based morphometry.25 In a small sample size of five patients, with IHI VS, decreased volumes were found in the inferior parietal and superior/medial frontal cortices, insula, temporal lobes, cingulum, midbrain, dorsal pons, caudate, cerebellum and thalamus.25 FDG-PET in the same set of patients demonstrated widespread hypometabolism in the parietal and parietoccipital cortices, the cingulum, bilateral thalami, and frontal, medial and precentral gyri. ‘Islands’ of preserved metabolism have also been seen in some subjects,26 a finding which may help account for the small numbers of patients who, despite no other behavioural evidence of complex cortical processing, may retain some residual fragmentary verbal output. In other patients, relatively complex cognitive processing may continue, but without any behavioural evidence, and may only be detected using fMRI.8 27 28 It appears, therefore, that VS may be the consequence not just of focal injury to discrete sites but also due to abnormal connectivity in distributed neuronal networks. Given this background, our demonstration of white-matter abnormalities in VS may provide an additional neuroanatomical substrate for the deficits observed. In addition, the comparison of differences in the two aetiologies may provide additional insights.

The extent of damage seen after IHI is determined by many factors which include the duration and degree of ischaemia, blood glucose levels and brain temperature during the period of insult.29 Certain brain regions exhibit a selective vulnerability to the effects of ischaemia including the cerebral cortex, striatum, hippocampus and cerebellar purkinje cells.1 2 Thalamic damage may be the predominant injury found at autopsy postcardiorespiratory arrest.30 The exact reasons for such increased regional susceptibility are unknown but may relate to higher metabolic rates in these regions as well as an increased proportion of excitotoxic glutaminergic cells.29 31 Our study demonstrates in vivo evidence for increased damage to particular areas including the thalamus and cerebral cortex, with the brainstem relatively preserved, which supports previous histopathological evidence32–35 as well as the widespread early (within the first 2 weeks of injury) patterns of restricted diffusion cortically, in the cerebellum and basal ganglia.32 35–37

The shearing forces experienced in TBI, and especially those with a component of traumatic axonal injury, mean that both the mechanism and the distribution of injury are different to IHI. Of course, TBI patients make also have a component of ischaemic injury, but in general the areas damaged are distinct histologically for TBI and IHI. Areas of damage in TBI are classically located in the parasagittal WM, corpus callosum (particularly the splenium/posterior aspect), and the corona radiata of the frontal and temporal lobes.38 In our study, thalamic damage was present in both patient populations, as evidenced by higher ADCs. The higher thalamic FA in hypoxic patients compared with TBI may represent the histopathological differences of predominantly neuronal loss versus retrograde degeneration, but true correlation with post-mortem studies would be required to see if this is indeed the case. In addition to abnormalities in these supratentorial regions, additional areas of damage are classically located in the dorsolateral brainstem.

These differences in the location of lesions in the two aetiologies are precisely reflected by our DTI findings. In contrast to the clear commonality of pathology in the supratentorial compartment in the two aetiologies that we studied, we found that patients with TBI exclusively demonstrated evidence of white-matter disconnection in the brainstem (figure 2). These differences underline the mechanistic heterogeneity that may result in VS. In addition, these findings may have important implications for selecting patients for interventions, or for stratification in clinical trials of novel interventions. There is increasing interest in using a range of therapies to enhance arousal, awareness and behavioural recovery in patients who are clinically categorised as VS. It seems self evident that not all patients are likely to respond to either pharmacological therapies or devices meant to enhance endogenous ascending neurochemical arousal systems. While functional imaging may (uniquely) provide the only evidence of preserved cortical processing in many such patients, detection of such processing may not predicate a higher likelihood of responding to therapy. Indeed, the patient who shows complex volitional control of cortical activation27 but is unable to manifest such activity through behavioural outputs may be least likely to respond to ‘arousal’ therapies, however intrinsically effective. Such a patient has adequate arousal; the problem is with a block in generating behavioural outputs. On the other hand, patients who show no evidence of cortical activity but have relatively preserved cortical structural integrity coupled with damage to key (brainstem) arousal systems and pathways may arguably benefit most from pharmacological39–42 or device-based enhancement43 of ascending arousal systems. This is a complex issue, and the inferences we outline above clearly remain unproven. However, these data do set novel hypotheses that can be tested. Although we have detected differences based on aetiology, findings in individual subjects, rather than aetiological groups, are the more important issue. However, the discussion above makes a case for detailed imaging of white-matter anatomy in studies that assess the impact of such therapies.

Our study has clear limitations. Despite a relatively large patient cohort in the context of previous studies, our sample size remains small, and our findings require confirmation. Further, longitudinal studies could assess the relationship between behavioural deficits and cognitive recovery on one hand, and structural and functional plasticity on the other. MRI in these patients can be challenging, and distorted anatomy may hinder traditional analysis methods—the degree of atrophy and distortion mean that careful ROI placement was required so that the chosen ROIs were easily identifiable even in abnormal brains.

Despite these limitations, the combination of conventional MRI, DTI and functional imaging may be important in this patient population. The demonstration of covert perceptual and cognitive processing using fMRI is now well established in selected patients,8 27 28 but it remains unclear why only some patients exhibit such responses. The presence or absence of such functional imaging responses may be able to be explained by the distribution and/or extent of the structural damage detected using DTI techniques. In addition, as we discuss above, DTI may potentially be used to help stratify patients into those likely to respond and may allow clinicians to select subjects for more careful assessment using fMRI.


We are very grateful to M Coleman for access to the imaging and clinical data used in this study.



  • Linked articles 205815.

  • Funding This research was supported by: UK National Institute of Health Research Biomedical Research Centre at Cambridge Technology Platform funding provided by the UK Department of Health. VFJN was supported by the Gates Cambridge Trust and an Overseas Research Studentship. JC and DS are supported by the NIHR Cambridge Biomedical Research Centre.

  • Competing interests None.

  • Patient consent Obtained.

  • Ethics approval Ethics approval was provided by Cambridgeshire 2 Research Ethics Committee.

  • Provenance and peer review Not commissioned; externally peer reviewed.

Linked Articles