Joint hypermobility is one visible manifestation of familial connective tissue variants that can impact organ function throughout the body. Joint hypermobility often results in troublesome joint pain and stiffness,Reference Keer and Grahame1 yet remains under-recognised and poorly managed. The diagnosis of joint hypermobility syndrome (HMS), as defined by the Brighton criteria, requires joint hypermobility plus musculoskeletal or connective tissue symptoms (e.g. prolapse, easy bruising, dislocations).Reference Grahame, Bird and Child2 This classification has now been superseded by the term Hypermobility Spectrum Disorder.Reference Castori, Tinkle, Levy, Grahame, Malfait and Hakim3 Rates of anxiety are markedly higher (odds ratio of 4.39) among hypermobile individuals.Reference Smith, Easton, Bacon, Jerman, Armon and Poland4 There is an overrepresentation of individuals with hypermobility in people with anxiety-related conditions and in presentations in which anxiety frequently co-occurs, including neurodevelopmental conditions such as autism and attention-deficit hyperactivity disorder.Reference Cederlöf, Larsson, Lichtenstein, Almqvist, Serlachius and Ludvigsson5,Reference Csecs, Iodice, Rae, Brooke, Simmons and Quadt6 Furthermore, individuals with hypermobility often experience symptoms of dysautonomia, such as postural tachycardia syndrome.Reference Hakim7 One explanation, connecting joint hypermobility, autonomic dysfunction and anxiety, proposes that relative inelasticity of connective tissue within peripheral vasculature compromises vasoconstriction and reduces venous return during standing through venous pooling. Compensatory autonomic responses, including increased sympathetic activity, increases physiological arousal including heart rate.Reference Bohora8–Reference Sharp, Critchley and Eccles11
In addition, hypermobile individuals may also show differences in interoceptive attention and sensitivity (increased sensing of changes from within the bodyReference Domschke, Stevens, Pfleiderer and Gerlach12), reflecting the experience of greater interoceptive surprise through less predictable (more imprecise) afferent visceral signals.Reference Paulus and Stein13 Increased attention may amplify interoceptive prediction error signals that contribute to the feeling of anxiety. Correspondingly, within the brain, hypermobile individuals are reported to show heightened reactivity in response to affective stimulation, both in regions responsible for interoceptive representation (insular cortex) and for emotional processing (amygdala),Reference Mallorquí-Bagué, Garfinkel, Engels, Eccles, Pailhez and Bulbena14 where structural differences are reported even in people with subclinical hypermobile features.Reference Eccles, Beacher, Gray, Jones, Minati and Harrison15
The aim of this study is to use functional neuroimaging to explore the neural basis for the relationship between joint hypermobility and clinical anxiety, building on this earlier work.Reference Mallorquí-Bagué, Garfinkel, Engels, Eccles, Pailhez and Bulbena14,Reference Eccles, Beacher, Gray, Jones, Minati and Harrison15 We hypothesised first that participants with anxiety would exhibit heightened insula and amygdala reactivity when processing social emotional stimuli (facial expressions), replicating prior findings.Reference Paulus and Stein13,Reference Klumpp, Post, Angstadt, Fitzgerald and Phan16–Reference Stein18 No previous functional imaging work has specifically addressed the link between hypermobility and clinical anxiety; we therefore additionally hypothesised that the reactivity of the amygdala and insula in conjunction with engagement of other ‘body-related’ brain regions would vary according to the presence and absence of hypermobility and anxiety, thereby illuminating neural substrates underlying their interaction.Reference Mallorquí-Bagué, Garfinkel, Engels, Eccles, Pailhez and Bulbena14,Reference Eccles, Beacher, Gray, Jones, Minati and Harrison15
Method
Participants
Sixty-three participants were recruited to the study. People volunteered to participate in response to an advertisement from Sussex Partnership NHS Trust either after inclusion in a linked study or via electronic bulletin boards. Members of the non-clinical comparison group were recruited via electronic bulletin boards. Of the 63 participants, 30 (47.6%) participants (age; mean ± s.e.m.) = 42.93 ± 2.24 yrs, 18 female, 12 male) met the threshold for generalised anxiety disorder (GAD), and 33 (52.4%) participants (age; mean ± s.e.m. = 37.42 ± 2.28 yrs, 16 female, 17 male) were healthy controls. There were no statistically significant differences in age or gender between the two groups. Of the people with GAD, 18 (60%) were classified as having joint hypermobility syndrome, and 7 (21.2%) of the non-clinical comparison group met the diagnostic threshold for joint hypermobility syndrome. See Table 1 for participant demographic details and clinical features.
HMS, hypermobility syndrome.
Inclusion criteria for people included the DSM-IV diagnosis of GAD. Participants in the non-clinical comparison group were required to be free from any history of psychiatric disorder. General exclusion criteria included magnetic resonance imaging (MRI) incompatibility, presence of neurological illness and presence of diagnosed psychiatric illness other than anxiety or comorbid depression in people. Written informed consent was obtained from all participants. The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2013. All procedures involving human participants including individuals were approved by the National Research Ethics Service – South East Coast (Brighton and Sussex; REC reference 12/LO/1942; IRAS registration number 115219).
All participants underwent assessment for GAD and hypermobility. DSM-IV diagnosis of GAD was confirmed or refuted using the Mini International Neuropsychiatric Interview (MINI).Reference Sheehan, Lecrubier, Sheehan, Amorim, Janavs and Weiller19 The presence or absence of generalised joint laxity was established through physical examination of joints using the Beighton Scale Reference Beighton, Solomon and Soskolne20 where a cut-off of 4 out of 9 was used in line with the UK literature, e.g. Clinch et al.Reference Clinch, Deere, Sayers, Palmer, Riddoch and Tobias21 All participants were assessed by the same clinician (J.A.E.). The presence of HMS was confirmed or excluded using Brighton Criteria.Reference Grahame, Bird and Child2
Anxiety level was measured using the Beck Anxiety Inventory (BAI).Reference Beck, Epstein, Brown and Steer22There was significant group difference in anxiety levels (BAI; mean ± s.e.m.: anxious group 22.9 ± 1.78 versus non-anxious group 5.48 ± 0.70; t(61) = 9.11, P < 0.001). See Table 1 for anxiety levels measured using the BAI across groups. There was no significant difference in anxiety levels between those with HMS and those without, either in the anxious (mean ± s.e.m.: anxious group: HMS 23.28 ± 2.37 versus anxious group: non-HMS 22.42 ± 2.82; t(28) = −0.23, P = 0.818) or non-anxious group (non-anxious group: HMS 5.38 ± 0.74 versus non-anxious group: non-HMS 5.86 ± 1.92; t(31) = −0.27, P = 0.786). All clinical and demographic data were analysed with IBM SPSS Statistics v. 29.23
Emotional faces task
An emotional faces task was modified from Umeda and colleagues,Reference Umeda, Harrison, Gray, Mathias and Critchley24 wherein five classes of images of emotional faces from the Karolinska Directed Emotional Faces set (KDEF) (classes: angry, afraid, disgusted, happy and neutral)Reference Goeleven, De Raedt, Leyman and Verschuere25 were presented in a randomised order. Null events were pseudo-randomised and presented as fixation cross. These were also included to facilitate the identification of haemodynamic responses to stochastically ordered stimuli. There were 15 trials of each emotion category, and 21 null events, each lasting 4 s. During each face presentation, participants were asked to make an incidental judgement of whether they could see teeth or not (index or middle finger button press with right hand) to ensure attention to the stimuli.
MRI acquisition
Functional MRI data were acquired on a Siemens Avanto 1.5 Tesla with a 32-channel head coil (T2*-weighted echo planar images, repetition time 2520 ms, echo time 43 ms, 34 interleaved slices 3 mm thick with 0.6 mm interslice gap, in-plane resolution 3 × 3 mm). A T1 structural was acquired for registration (repetition time 2730 ms, echo time 3.57 ms, 1 × 1 × 1 mm resolution).
fMRI preprocessing
Functional MRI data were preprocessed and analysed using Statistical Parametric Mapping (SPM12)Reference Penny, Friston, Ashburner, Kiebel and Nichols26 running in MATLAB.27 Preprocessing was performed using default options, including realignment to the mean image, slice-time correction to the 6th slice, co-registration to the T1 structural image and normalisation to Montreal Neurological Institute space, as well as smoothing at 8 mm Gaussian smoothing kernel. To account for head motion, framewise displacement (FD) values were calculated from the six motion parameters (rp_.txt) generated in SPM during realignment using the FD_conn method in the CONN functional connectivity toolbox running in SPM12 (version 22.a).Reference Whitfield-Gabrieli and Nieto-Castanon28 Participants with excessive motion (framewise displacement > 0.5 mm) were flagged for closer inspection. Two participants exceeded this threshold (with framewise displacement values of 0.63 mm and 0.60 mm respectively) but were retained in the analysis to preserve the representativeness and generalisability of the findings within the clinical population.
Statistical analysis
First-level general linear model
Task events were modelled in a general linear model, with five regressors representing the onset and duration of presentation of angry, afraid, disgusted, neutral and happy faces respectively. To account for head motion, six nuisance regressors modelled head movement using the motion parameters calculated during realignment. Single-regressor T-contrasts were generated for viewing (a) angry, (b) afraid, (c) disgusted, (d) neutral and (e) happy faces by assigning a contrast weight of 1 to each of the five experimental conditions, with the intertrial interval fixation cross representing an implicit baseline. These T-contrasts were entered into a full factorial second-level analysis.
Second-level general linear model
A second-level full factorial model contained HMS (non-HMS, HMS) and anxiety (non-anxious, anxious) as an independent (between-subjects) factor, and facial expression (angry, afraid, disgusted, neutral and happy) as a within-subject factor. In addition, two covariates were entered for (a) gender (male, female) and (b) age.
F-contrasts were generated testing for: all effects, main effect of HMS, main effect of anxiety, main effect of task and interactions between the factors. Individual group effects for viewing faces (compared with implicit baseline) were examined using T-contrasts: HMS > non-HMS; HMS (anxious > non-anxious); non-HMS > HMS; and anxious (HMS > non-HMS). A series of two further second-level models included the Beighton score and anxiety level (BAI score) as additional covariates. In the Beighton second-level model, the Beighton score was used as a covariate so that the main effect of hypermobility symptoms could be modelled along with the interaction of hypermobility symptoms with the anxiety factor (i.e. presence of GAD or not). In the anxiety second-level model, the anxiety level (BAI score) was used as a covariate, so that the main effect of the anxiety level could be modelled along with the interaction of the anxiety level (BAI score) with the factor of HMS. All covariates were mean centred around zero.
Statistical images were thresholded at a cluster-forming threshold of P < 0.001 for family-wise error rate correction (FWEc) for multiple comparisons at P < 0.05.Reference Eklund, Nichols and Knutsson29 Significant clusters were localised according to the Anatomy toolbox running in SPM12 (v 3.0).Reference Eickhoff, Paus, Caspers, Grosbras, Evans and Zilles30
Psychophysiological interactions
We performed a series of psychophysiological interaction analyses to investigate how brain activity in response to emotional faces, within regions identified in the above univariate analyses, changed in their functional connectivity to other regions of the brain as a function of hypermobility and anxiety status. On the basis of the univariate fMRI results, we identified three regions from which to seed these functional connectivity analyses: (a) left amygdala (centred on x −32, y 0, z −16); (b) right mid insula (centred on x 34, y −2, z −6); and (c) left inferior frontal gyrus (x −46, y 34, z 2). First, we extracted the first eigenvariate (weighted mean of blood-oxygen-level-dependent (BOLD) time series) for each region by thresholding three contrasts at P < 1 for each participant: (a) the interaction between anxiety and HMS (for left amygdala region of interest (ROI)); (b) the main effect of HMS for (right mid insula ROI); and (c) the main effect of HMS (for left inferior frontal gyrus ROI). Then, an F-contrast was computed for each subject, representing all effects (angry, afraid, disgusted, neutral and happy: ‘eye:5’). In the three contrasts given above, we then extracted a 10 mm sphere of voxels for each ROI, adjusting for the F-contrast of all effects.
Next, the psychophysiological interaction term was calculated according to the main effect of the task (contrast weights: 1 for angry, 1 for afraid, 1 for disgusted, 1 for neutral and 1 for happy) and the BOLD time series for each ROI. These psychophysiological interaction terms were each entered into a first-level model for each participant, alongside a regressor representing the BOLD activity of the ROI (psychophysiological interation, original region of interest eigenvariate (PPI.Y)) and the main effect of the task (psychophysiological interation, attention–no attention task vector (PPI.P)). Single regressor T-contrasts were generated for the psychophysiological interaction term using a single contrast weight to investigate positive changes in the regression slope of voxels elsewhere in the brain relative to the seed ROI during task events relative to baseline.
The first-level T-contrasts were then entered into a series of second-level models that examined the psychophysiological interaction between the seed ROI and voxels across the brain using a full factorial second-level analysis, with HMS (non-HMS, HMS) and anxiety (non-anxious, anxious) as an independent (between-participant) factor, and the first-level T-contrasts representing functional connectivity when viewing faces as a repeated measures (within-participant) factor.
In these second-level models, as with the univariate functional MRI analysis, age and gender were entered as covariates. Contrasts were thresholded at a cluster-forming threshold of P < 0.001 for FWEc at P < 0.05. Significant clusters were localised with reference to the SPM Anatomy toolbox (v 2.2b).Reference Eickhoff, Paus, Caspers, Grosbras, Evans and Zilles30
Transparency declaration
We affirm that the manuscript is an honest, accurate and transparent account of the study being reported, and that no important aspects of the study have been omitted.
Results
Univariate functional MRI
Main effects
As anticipated, there was a significant main effect of anxiety within the left amygdala and right mid insula, with post hoc T-contrasts revealing greater activity in these regions in anxious versus non-anxious participants (Supplementary Material, Tables 4 and 5 available at https://doi.org/10.1192/bjp.2024.259).
We also observed significant main effects of HMS in the left inferior frontal gyrus, precuneus and pre-supplementary motor area (SMA), right mid insula, right posterior and left anterior mid cingulate gyrus and the left supramarginal gyrus (Supplementary Material, Table 1). Post hoc T-contrasts revealed that non-HMS participants showed greater activity in the left inferior frontal gyrus, precuneus, left pre-SMA, right posterior and left anterior mid cingulate gyrus and the left supramarginal gyrus when viewing emotional faces compared with participants with HMS (Fig. 1(a)). Post hoc T-contrasts for the right mid insula were not statistically significant.
We also observed the known main effect of viewing emotional faces, associated with activation of large areas of the occipital lobe, right middle and left superior frontal gyrus (Supplementary Material, Table 3).
Interaction: hypermobility × anxiety
Furthermore, there was a significant interaction between HMS and anxiety in the left amygdala, left hippocampus, right paracingulate gyrus and right mid insula (Supplementary Material, Table 6). Post hoc T-contrasts revealed that: (a) in the group with HMS, there was greater activation in the left amygdala and the right paracingulate gyrus in anxious compared with non-anxious participants (Fig. 1(b)); (b) findings for the left hippocampus were not significant; (c) in participants with anxiety, there was greater activation in the right paracingulate gyrus and right mid-insula in the group with HMS compared with the non-HMS group (Fig. 1(c)).
Interaction: Beighton score × anxiety
The interaction between the number of hypermobile joints (Beighton score) and anxiety status, i.e. the interaction testing for regions in which activation was more positively correlated with the Beighton score for anxious compared with non-anxious participants, showed activation in the left anterior insula (Fig. 2(a)).
Interaction: Beck anxiety score × hypermobility
The interaction between anxiety severity (BAI score) and HMS, i.e. the interaction testing for regions in which activation was more positively correlated with the anxiety score for people with HMS compared with non-HMS participants, showed activation in the left putamen (Supplementary Material, Table 10).
Psychophysiological interactions
Three second-level models examined changes in functional connectivity with the (a) left amygdala, (b) right mid insula and (c) left inferior frontal gyrus, depending on the psychological context of viewing emotional faces. In the left amygdala psychophysiological interaction, there was no significant effect of HMS (F-contrast; HMS, non-HMS). However, the post hoc T-contrast testing for HMS versus non-HMS (T-contrast; HMS > non-HMS [−1 −1 1 1]) revealed that participants with HMS, when viewing faces, showed greater functional connectivity between the left amygdala and the left precuneus (Fig. 2(b)).
In the right mid insula psychophysiological interaction, there was a significant interaction effect between HMS and anxiety severity in the right supramarginal gyrus (Fig. 2(c)) and the right occipital cortex. However, the post hoc T-contrasts to explore the effect did not remain significant after correction for multiple comparisons. The left inferior frontal gyrus psychophysiological interaction did not produce statistically significant results.
Discussion
Here, using functional neuroimaging during the incidental processing of emotional faces, we identify putative neural substrates underpinning the association between hypermobility and anxiety. Hypermobile participants showed reduced activation in discrete areas of the association cortex, notably the prefrontal and parietal regions. Anxious participants, in line with several previous imaging studies, showed amplified reactivity of the amygdala and insula during a socio-emotional challenge. For the first time, we were able to investigate the neural interaction between hypermobility and anxiety. We first confirmed that anxious participants with HMS showed greater activity in the amygdala and insula than non-anxious participants with HMS. Furthermore, the degree of hypermobility (as measured by the Beighton score) was more strongly correlated with insula activity in anxious participants with HMS than non-anxious participants with HMS. Finally, there was a general effect of hypermobility (regardless of anxiety) on functional connectivity between the amygdala and precuneus. Collectively, these findings identify a network of amygdala, insula and association cortices that link anxiety and hypermobility. Greater connectivity and activity within this network may underpin the increased prevalence of anxiety in individuals with HMS.
Interoceptive pathways in the brain
The brain continuously receives sensory information from the visceral organs and peripheral tissues via ascending nerves that enter the brainstem.Reference Craig31 From here, interoceptive signals relaying the state of the body are conveyed to the thalamus, and ultimately the posterior insula cortex. Here, viscerotopic representations of these afferent signals are believed to support the cognitive perception of bodily feelings (e.g. heart rate, respiration, gastric sensations). This viscerotopic information is then re-represented, and integrated, more anteriorly in the insula lobe, underpinning our experience of broader affective states.Reference (Bud) Craig32,Reference Critchley and Garfinkel33
Putatively, hypermobility may render individuals more prone to anxious affective experiences via heightened signalling of interoceptive signals relaying dysautonomic states. Because of changes in the connective tissue of the vasculature, hypermobile people may experience abnormal peripheral vasoconstriction.Reference Csecs, Dowell, Savage, Iodice, Mathias and Critchley34 Specifically, reduced venous return during standing because of venous pooling may be responsible for an increased sympathetic state and autonomic hyperactivity.Reference Bohora8,Reference Mathias, Low, Iodice, Owens, Kirbis and Grahame9 The insular cortex is an important central substrate for receiving this autonomic hyperactivity information.
It is therefore particularly intriguing that the insula was not only identified as overactive in our anxious sample, but more specifically, as more active in anxious versus non-anxious participants with HMS. We also saw that the degree of hypermobility (Beighton score) was more strongly correlated with insula activity in anxious than non-anxious participants with HMS. This identifies the insula as a nexus of affective experience in anxious hypermobile people. This finding extends similar observations reported in other affective conditionsReference Förster, Maliske, Schurz, Henneberg, Dannlowski and Kanske35 (e.g. increased emotional reactivity (hyperactivity of salience-processing regions) in bipolar disorder) and previous findings from our group that have examined the functional activity in non-anxious hypermobile participants.Reference Mallorquí-Bagué, Garfinkel, Engels, Eccles, Pailhez and Bulbena14 Within this context, interoception, i.e. the dynamic signalling, neural and perceptual representation of internal physiological states of the body, is a likely unifying factor. Additionally, participants with bipolar disorder have demonstrated abnormal insular functional connectivity, possibly modulated by inflammatory markers.Reference Chen, Chen, Chen, Zhong, Gong and Zhong36 Similar mechanisms may underpin the heightened insula activity observed in anxious hypermobile participants, potentially mediated by dysautonomic states.
Amygdala interactions with hypermobility
In addition to reactivity differences in the insular cortex, we also observed interesting findings in the amygdala. Affective tasks performed during fMRI reliably engage the amygdala. This activation is typically amplified in anxious participantsReference Fonzo and Etkin37 and is also enhanced following interoceptive stressors such as immune challenges.Reference Davies, Cooper, Voon, Tibble, Cercignani and Harrison38 It is thus noteworthy that the amygdala was identified as a functional neural centre for interaction between anxiety and hypermobility, wherein the amygdala showed even greater activity in anxious participants with HMS than the non-anxious group with HMS. The amygdala is a critical region supporting the detection and perception of threat through associative integration of external and interoceptive information,Reference Bach, Hurlemann and Dolan39 which may underpin its role in affective experience.Reference Garfinkel, Minati, Gray, Seth, Dolan and Critchley40 The greater reactivity of the amygdala in anxious hypermobile individuals may also reflect previously identified differences of amygdalar structure and function in hypermobile people.Reference Mallorquí-Bagué, Garfinkel, Engels, Eccles, Pailhez and Bulbena14,Reference Eccles, Beacher, Gray, Jones, Minati and Harrison15 In parallel, the amygdala response in anxious participants with HMS may reflect the dynamic contributions (e.g. to behavioural/autonomic response and subjective feelings) of the amygdala and the insula within a wider affective network; however, our functional connectivity analyses did not identify a dependent association between these two regions.
Functional connectivity analyses did, however, identify stronger functional coupling between the amygdala and precuneus in hypermobile participants (regardless of anxiety). This may suggest that hypermobile individuals have a tendency towards hyper-connectivity of affective regions, which may have consequences for onward processing of information.
Affective tasks as a probe for functional anatomy
We selected an emotional face processing task as a vehicle for probing affective responses in participants with GAD. A wide range of literature has previously taken such an approach to confirm the involvement of regions such as the amygdala and insula in a variety of anxiety conditions,Reference Engel, Bandelow, Gruber and Wedekind41 as well as experimental interoceptive challenges.Reference Davies, Cooper, Voon, Tibble, Cercignani and Harrison38 Here, we were able to extend this literature to understand the impact of hypermobility on these processes.
In our analyses, we collapsed across the five stimuli types to ask the fundamental question of how viewing social affective stimuli can provoke a neural response. In future work, it would be interesting to further tease apart the nuances of different affective cues on the hypermobile brain. For example, do responses to anger differ from fear?
Limitations and future directions
Despite our implementation of a well-established paradigm to invoke reliable activations in affective brain regions, we acknowledge that fMRI tasks constrain one's investigative potential to the circumscribed set of regions that the task recruits. Other approaches, including resting-state fMRI, provide complementary information that can be leveraged to examine brain-wise network interactions. We did not acquire such data within this study, yet this remains an important avenue for future work to better understand the neural characteristics of hypermobility.
We recruited a community sample that enabled us to screen participants and place them into one of four categories according to anxiety and hypermobility status. However, since hypermobility is a risk factor for anxiety, with an odds ratio of 4.39 for suffering from anxiety if hypermobile,Reference Smith, Easton, Bacon, Jerman, Armon and Poland4 recruiting large numbers of hypermobile individuals who are not anxious (and anxious individuals who are not hypermobile) was a challenge. This means our sample sizes for the non-anxious HMS and anxious non-HMS groups were smaller than the other two groups.
Future work should capitalise on advances in biofeedback interventions (such as aligning dimensions of interoceptive experienceReference Quadt, Garfinkel, Mulcahy, Larsson, Silva and Jones42 and altering dynamics of autonomic processing therapiesReference Davies, Csecs, Ball, Dare, Bremner and Hosking43) to target anxiety in hypermobile individuals, perhaps using neurofeedback to specifically down-regulate insular and amygdala reactivity. As in the broader mental health space, it is becoming increasingly clear that a ‘one-size-fits-all’ approach does not work for many groups experiencing anxiety; the present study will potentially inform personalised treatment approachesReference Williams, Carpenter, Carretta, Papanastasiou and Vaidyanathan44 for a group of individuals who have previously perhaps been dismissed or overlooked.Reference Halverson, Penwell and Francomano45
Supplementary material
Supplementary material is available online at https://doi.org/10.1192/bjp.2024.259
Data availability
The analytic code (SPM batches) that was used for the neuroimaging analysis for this study is available at OSF: osf.io/tcemx/. The MRI acquisition sequence information, demographic and clinical data, and participant mean framewise displacement values are available at OSF: osf.io/tcemx/. The neuroimaging data that support the findings of this study (unthresholded statistic images for every contrast reported) are openly available at https://identifiers.org/neurovault.collection:16863, reference number 16863.Reference Hill, Laird, Marcus, Gorgolewski, Varoquaux and Rivera46
Acknowledgements
We would like to thank the radiography team and administrative staff at the Clinical Imaging Sciences Centre, Brighton and Sussex Medical School at the University of Sussex, UK for all their support throughout this project.
Author contributions
C.N.K.: methodology, formal analysis, visualisation, writing – original draft, writing – review & editing; C.L.R.: methodology, formal analysis, supervision, writing – original draft, writing – review & editing; C.G.V.P.: methodology, writing – review & editing; N.A.H.: writing – review & editing; S.N.G.: writing – review & editing; H.D.C.: conceptualisation, supervision, writing – review & editing; J.A.E.: conceptualisation, project administration, methodology, investigation, data curation, formal analysis, supervision, writing – original draft, writing – review & editing.
Declaration of interest
None.
Funding
Funding for this project came via a fellowship to J.A.E. (MRC MR/K002643/1). J.A.E. was also supported by MQ Transforming Mental Health and Versus Arthritis (MQF 17/19).
eLetters
No eLetters have been published for this article.