
==== Front
BMJ Open
BMJ Open
bmjopen
bmjopen
BMJ Open
2044-6055
BMJ Publishing Group BMA House, Tavistock Square, London, WC1H 9JR

39231552
10.1136/bmjopen-2024-085814
bmjopen-2024-085814
Original Research
Rheumatology
1732
1506
Quantitative sensory testing and chronic pain syndromes: a cross-sectional study from TwinsUK
https://twitter.com/AmberRhee_EPI
http://orcid.org/0000-0003-4894-6578
Rhee Amber 12arhee6@jhmi.edu

http://orcid.org/0000-0003-4294-3317
Granville Smith Isabelle 1isabelle.smith@kcl.ac.uk

http://orcid.org/0000-0002-7034-1092
Compte Roger 1roger.compte@kcl.ac.uk

http://orcid.org/0000-0003-2804-7399
Vehof Jelle 134jelle.vehof@kcl.ac.uk

http://orcid.org/0000-0002-1993-1345
Nessa Ayrun 1ayrun.nessa@kcl.ac.uk

http://orcid.org/0009-0000-7289-8436
Wadge Samuel 1sam.wadge@kcl.ac.uk

http://orcid.org/0000-0002-1439-6259
Freidin Maxim B 15m.freydin@qmul.ac.uk

http://orcid.org/0000-0002-7996-2696
Bennett David L 6david.bennett@ndcn.ox.ac.uk

http://orcid.org/0000-0002-2998-2744
Williams Frances M K 1frances.williams@kcl.ac.uk

1 Department of Twin Research and Genetic Epidemiology, King's College London, London, UK
2 Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, USA
3 Departments of Ophthalmology and Epidemiology, University of Groningen, Groningen, The Netherlands
4 Department of Ophthalmology, Vestfold Hospital Trust, Tønsberg, Norway
5 Department of Biology, Queen Mary University of London, London, UK
6 Nuffield Department of Clinical Neurosciences, Oxford University, Oxford, UK
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.

None declared.

IsabelleGranville Smith; isabelle.smith@kcl.ac.uk
2024
3 9 2024
14 9 e08581401 3 2024
30 7 2024
Copyright © Author(s) (or their employer(s)) 2024. Re-use permitted under CC BY. Published by BMJ.
2024
https://creativecommons.org/licenses/by/4.0/ 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/.

Abstract

Objective

The chronic pain syndromes (CPS) include syndromes such as chronic widespread pain (CWP), dry eye disease (DED) and irritable bowel syndrome (IBS). Highly prevalent and lacking pathognomonic biomarkers, the CPS are known to cluster in individuals in part due to their genetic overlap, but patient diagnosis can be difficult. The success of quantitative sensory testing (QST) and inflammatory biomarkers as phenotyping tools in conditions such as painful neuropathies warrant their investigation in CPS. We aimed to examine whether individual QST modalities and candidate inflammatory markers were associated with CWP, DED or IBS in a large, highly phenotyped population sample.

Design

Cross-sectional study.

Setting

Community-dwelling cohort.

Participants

Twins from the TwinsUK cohort

Primary and secondary outcome measures

We compared 10 QST modalities, measured in participants with and without a CWP diagnosis between 2007 and 2012. We investigated whether inflammatory markers measured by Olink were associated with CWP, including interleukin-6 (IL-6), IL-8, IL-10, monocyte chemoattractant protein-1 and tumour necrosis factor. All analyses were repeated in DED and IBS with correction for multiple testing.

Results

In N=3022 twins (95.8% women), no association was identified between individual QST modalities and CPS diagnoses (CWP, DED and IBS). Analyses of candidate inflammatory marker levels and CPS diagnoses in n=1368 twins also failed to meet statistical significance.

Conclusion

Our findings in a large population cohort suggest a lack of true association between singular QST modalities or candidate inflammatory markers and CPS.

Chronic Pain
Cross-Sectional Studies
EPIDEMIOLOGIC STUDIES
Irritable Bowel Syndrome
Dry Eye Syndromes
http://dx.doi.org/10.13039/100010269 Wellcome Trust 212904/Z/18/Z
==== Body
pmcSTRENGTHS AND LIMITATIONS OF THIS STUDY

Large sample (N=3022) from a well-characterised, cohort.

Appropriate correction for multiple testing across all analyses.

Adjusted for major covariates (family relatedness, age, body mass index) in all analyses.

Inclusion of participants with common painful conditions other than chronic pain syndromes in analysis.

Variation in sample sizes for each quantitative sensory testing modality.

Introduction

Chronic pain is a major public health burden and leading cause of disability globally.1 Because of its clinical heterogeneity, the necessity of pain phenotyping and measurement tools are widely acknowledged.2 Chronic pain syndromes (CPS) are a recognised cluster of highly prevalent syndromes, including chronic widespread musculoskeletal pain (CWP)/fibromyalgia, dry eye disease (DED) and irritable bowel syndrome (IBS).36 CPS demonstrate genetic and symptomatology overlaps in the TwinsUK cohort, a large population-based cohort of community-dwelling twins, a finding supported in other population-based cohorts and clinical studies.379 Notably, CPS are characterised by lack of pathognomonic tissue injury or clinical biomarkers.10 Phenotyping and stratification, crucial steps to improving therapy success, thus prove difficult in these syndromes.

Quantitative sensory testing (QST) is a collection of psychophysical tests that assess percepts evoked by defined sensory stimuli.11 The tests evaluate a range of sensory modalities (ie, thermal and mechanical) in both the non-noxious and noxious range (ie, warm detection and heat pain thresholds (HPT)), providing insights into somatosensory nervous system functions.12 13 QST protocols have been standardised and are increasingly applied to clinical cohorts for grading neuropathic pain (pain caused by damage to the somatosensory nervous system) and subclassifying chronic pain.1418 For instance, in neuropathic pain, QST can help build a sensory profile via an unbiased clustering algorithm, generating three principal groups: sensory loss, thermal hyperalgesia and mechanical hyperalgesia.19 Broadening the utilisation of QST beyond neuropathic pain is a subject of debate, however, as the extent of insight QST provides into the pathological mechanisms driving chronic pain is unclear.2023

Inflammatory markers have been studied as potential phenotyping tools in chronic pain conditions. For example, osteoarthritis studies have identified biomarker profiles characteristic to the inflammatory osteoarthritis phenotype in synovial fluid.24 While CPS such as CWP were previously thought not inflammatory, some evidence suggests low-grade inflammation in CPS.2528 Unlike QST, inflammatory markers offer direct insight into circulating mediators of disease state.24

This cross-sectional study aimed to examine the potential for QST and inflammatory markers to serve as phenotyping tools in CPS by investigating whether individual QST modalities and candidate inflammatory markers were associated with CPS diagnoses in a population-based cohort.

Methods

Participants

TwinsUK is an adult twin registry with over 15 000 volunteers, recruited from across the UK for research. Established in 1992, the ongoing longitudinal cohort, initially restricted to female recruitment, is predominantly white race/ethnicity and female (82%), with age mean 59.6 years (range 18–90+ years), and comprises monozygotic (MZ) and dizygotic (DZ) twins.29 TwinsUK has one of the largest QST datasets, administered to TwinsUK participants between 2007 and 2012 as part of collaborative studies with Pfizer.30 The cohort has been reported similar to age-matched British women from a singleton population-based cohort with regards to a range of health traits and diseases.31 Phenotypic data and biological specimens are collected from TwinsUK participants through annual questionnaires and approximately quadrennial clinic visits.29

Participants from TwinsUK were included in this cross-sectional study if they previously completed at least one QST measure and at least one CPS questionnaire. Participants were excluded from QST at the time of visit if they had severe skin disease, previous stroke or chemotherapy, likely impaired upper limb neurology, allergy to electrodes, history of melanoma, were pregnant or used painkillers on the day of the test.30

Quantitative sensory testing protocols

In this study, we examined 10 QST modalities: cold intolerable threshold, cold pain threshold, heat pain supra threshold (HPST), HPT, mechanical detection threshold, mechanical pain threshold, skin flare extent, pain during burn induction, punctate hyperalgesia and thermal hyperalgesia. The last four modalities are part of a milder thermal burn protocol than the more commonly utilised protocol.32 QST was administered at a single site during standard TwinsUK visits. Protocols were established in TwinsUK in collaboration with the Stephen McMahon lab, King’s College London under the auspices of Prof D Bennett (coauthor, now at Oxford).32 Detailed descriptions of QST protocols are found in online supplemental file S1. A high degree of standardisation is necessary to perform QST accurately; to achieve this, nurses and research assistants underwent considerable training. Heritability and reliability of QST measures in this particular population have been formally assessed and reported previously, with heritability and inter-rater reliability estimates for each modality ranging from 0.29 to 0.55 and 0.34 to 0.91 respectively.30 32

Candidate inflammatory marker measures

Five ‘candidate’ inflammatory markers were compiled a priori as exposure variables for secondary analysis: interleukin-6 (IL-6), IL-8, IL-10, monocyte chemoattractant protein-1 (MCP-1) and tumour necrosis factor (TNF). Markers were selected from the Olink Target 96 Inflammation panel to reflect both current literature and assay availability.262833 Serum inflammatory marker proteomics were collected and assayed as part of a large proteomics study. A subset of TwinsUK volunteers participated in both the QST and proteomics studies.

Olink uses Proximity Extension Assay to assess multiple proteins and report levels through a preprocessed relative (Normalised Protein eXpression) quantification on a log2 scale.39 40 Proteomics were assayed in two batches in 2019 and 2020, and data from the two batches were combined by author MBF into a single dataset through bridge sample normalisation according to Olink recommendations.41 42 Only samples collected within 2 years of participant QST visit were included in this study. For participants with multiple longitudinal samples, the sample collected on QST visit or closest to the participant’s QST visit was used in analysis.

Chronic pain syndromes

CWP status was ascertained using a modified version of the London Fibromyalgia Epidemiology Study Screening Questionnaire.43 DED classification was determined according to the validated Women’s Health Study questionnaire, while IBS status was determined based on Rome 3 Criteria and, if unavailable, self-report of clinician diagnosis or treatment.44 Questionnaires were administered between 2002 and 2020. Participants were counted as cases if they were ever diagnosed with a CPS during this time.

Statistical analysis

For each of the 10 QST modalities, we conducted a Mann-Whitney U test among participants who completed a CWP questionnaire; we compared scores between participants with CWP and participants without CWP (ie, comparing HPT between participants with CWP and participants without CWP). This was repeated for DED and IBS. With 10 QST modalities, a Bonferroni-correction cut-off was set at p=0.005.

To account for potential confounding due to CPS comorbidities, we conducted Mann-Whitney U tests in a sensitivity analysis comparing QST scores between participants with CWP and a control group of participants without any of the CPS diagnoses (‘true controls’). To address temporality, a secondary sensitivity analysis compared QST scores of participants with prevalent CWP diagnosis at QST date and participants with incident CWP diagnosis after QST date. In consideration of family relatedness and potential confounding by age and body mass index (BMI), we repeated the main analysis using mixed effects logistic regressions of each QST modality (scaled) on CWP diagnosis (ie, regression of HPST (scaled) on CWP diagnosis) with family ID as a random effect and age (scaled) and BMI category (nominal) as fixed effects. We used a BOBYQA (Bound Optimization BY Quadratic Approximation) optimisation technique, using the lme4 package in R.45 BMI categories were defined according to the Centers for Disease Control and Prevention BMI cut-off standards.46 All sensitivity analyses were repeated for DED and IBS.

For each of the five candidate inflammatory markers, we conducted a mixed effects logistic regression of inflammatory marker level on CWP diagnosis (ie, regression of IL-6 on CWP diagnosis) in a subset of participants who had Olink proteomics data collected within 2 years of their QST data collection visit. Model specifications were identical to those of the regression analyses in QST. Fixed effects included age (scaled) and BMI category (nominal). Family ID was included as a random effect to control for twin relatedness. Men were excluded from inflammatory panel analyses due to small sample size. This was repeated for DED and IBS. In total, we considered 15 models, examining five markers in each CPS, and imposed a Bonferroni-corrected p value cut-off of 0.003.

In a sensitivity analysis, we conducted a discordant twin analysis (MZ and DZ) of inflammatory marker levels on CWP diagnosis in twin pairs who were discordant for CWP. Using a conditional logistic regression analysis in the R survival package, associations were adjusted for BMI category (nominal).47 This analysis was restricted to samples collected on the same day as QST visit. We repeated these analyses in DED and IBS and imposed a Bonferroni-corrected p-value cut-off of 0.003.

We analysed all data using R V.4.2.1 (R Foundation for Statistical Computing, Vienna, Austria).

Patient and public involvement

Patients or the public were not involved in the design, or conduct, or reporting or dissemination plans of our research.

Results

Participants with CWP questionnaire data (N=2996) completed at least one QST modality (table 1). Prevalence of CWP was 22.4%; n=564 reported CWP at the time of QST visit and n=106 developed incident CWP after QST visit. More participants with CWP were classified as obese (26.9%) compared with those without CWP (16.9%).

Table 1 Characteristics of participants who completed a CWP questionnaire and QST

	Total* (N=2996)	Participants without CWP (n=2326)	Participants with CWP (n=670)	
Zygosity	
 Dizygotic (%)	1435 (47.9)	1078 (46.3)	357 (53.3)	
 Monozygotic (%)	1550 (51.7)	1239 (53.3)	311 (46.4)	
 Missing (%)	11 (0.4)	9 (0.4)	2 (0.3)	
Sex	
 Female (%)	2872 (95.9)	2220 (95.4)	652 (97.3)	
 Male (%)	124 (4.1)	106 (4.6)	18 (2.7)	
Age (years)	
 Mean (SD)	57.4 (12.2)	56.2 (12.7)	61.3 (9.3)	
BMI category	
 Underweight (%)	35 (1.2)	29 (1.2)	6 (0.9)	
 Healthy (%)	1326 (44.3)	1084 (46.6)	242 (36.1)	
 Overweight (%)	1063 (35.5)	821 (35.3)	242 (36.1)	
 Obese (%)	572 (19.1)	392 (16.9)	180 (26.9)	
* Total participants in TwinsUK who have completed a CWP questionnaire and at least one QST modality.

BMI, body mass indexCWP, chronic widespread pain; QST, quantitative sensory testing

Among participants with DED questionnaire data, N=2583 completed at least one QST modality (table 2). With prevalence of DED at 28.8%, approximately half of the DED cases (n=358) were prevalent at QST visit; n=387 were incident and developed after QST visit.

Table 2 Characteristics of participants who completed a DED questionnaire and QST

	Total* (N=2583)	Participants without DED (n=1838)	Participants with DED (n=745)	
Zygosity	
 Dizygotic (%)	1234 (47.8)	890 (48.4)	344 (46.2)	
 Monozygotic (%)	1340 (51.9)	939 (51.1)	401 (53.8)	
 Missing (%)	9 (0.3)	9 (0.5)	0 (0.0)	
Sex	
 Female (%)	2485 (96.2)	1750 (95.2)	735 (98.7)	
 Male (%)	98 (3.8)	88 (4.8)	10 (1.3)	
Age (years)	
 Mean (SD)	58.3 (11.2)	57.6 (11.8)	60.2 (9.7)	
BMI category	
 Underweight (%)	30 (1.2)	21 (1.1)	9 (1.2)	
 Healthy (%)	1159 (44.9)	793 (43.1)	366 (49.1)	
 Overweight (%)	919 (35.6)	674 (36.7)	245 (32.9)	
 Obese (%)	475 (18.4)	350 (19.0)	125 (16.8)	
* Total participants in TwinsUK who have completed a DED questionnaire and at least one QST modality.

BMI, body mass indexDED, dry eye disease; QST, quantitative sensory testing

Participants with IBS questionnaire data (N=2677) completed at least one QST modality (table 3). Prevalence of IBS was 26.2%, n=368 of which was prevalent at QST visit and n=334 of which was incident after QST visit.

Table 3 Characteristics of participants who completed an IBS questionnaire and QST

	Total* (N=2677)	Participants without IBS (n=1975)	Participants with IBS (n=702)	
Zygosity	
 Dizygotic (%)	1272 (47.5)	945 (47.8)	327 (46.6)	
 Monozygotic (%)	1397 (52.2)	1022 (51.7)	375 (53.4)	
 Missing (%)	8 (0.3)	8 (0.4)	0 (0.0)	
Sex	
 Female (%)	2580 (96.4)	1895 (95.9)	685 (97.6)	
 Male (%)	97 (3.6)	80 (4.1)	17 (2.4)	
Age (years)	
 Mean (SD)	58.0 (11.5)	58.2 (11.4)	57.4 (11.7)	
BMI category	
 Underweight (%)	30 (1.1)	23 (1.2)	7 (1.0)	
 Healthy (%)	1202 (44.9)	873 (44.2)	329 (46.9)	
 Overweight (%)	949 (35.5)	718 (36.4)	231 (32.9)	
 Obese (%)	496 (18.5)	361 (18.3)	135 (19.2)	
* Total participants in TwinsUK who have completed an IBS questionnaire and at least one QST modality.

BMI, body mass indexIBS, irritable bowel syndrome; QST, quantitative sensory testing

In total, N=3022 unique participants were included across analyses. Within QST participants who completed all three CPS questionnaires (n=2502; 82.8%), n=1156 were true controls and never diagnosed with any CPS. Overlap of each analytical group and their CPS diagnoses are demonstrated in online supplemental figure S1.

Most participants completed only heat QST modalities for the Pfizer study (n=2633) with n=365 participants completing both heat and mechanical QST modalities. Sample sizes for each QST modality are available in table 4.

Table 4 QST sample sizes by modality and CPS analytical group

QST	CWP	DED	IBS	Total participants	
CIT	100	95	94	100	
CPT	133	127	126	133	
Flare extent	101	99	98	101	
HPST	2276	2066	2125	2298	
HPT	2973	2561	2656	2999	
MDT	380	360	357	380	
MPT	386	364	362	386	
Pain during burn induction	101	99	98	101	
Punctate hyperalgesia	100	98	97	100	
Thermal hyperalgesia	101	99	98	101	
Each cell represents the number of participants who completed the corresponding QST and CPS questionnaire (ie, participants who completed the HPT test and CWP questionnaire=2973). The total participants’ column indicates the total unique participants across all analytical groups who have completed the corresponding QST test (ie, total participants who completed HPT test and any of the three CPS questionnaires=2999).

CIT, cold intolerable threshold; CPS, chronic pain syndrome; CPT, cold painful threshold; CWP, chronic widespread pain; DED, dry eye disease; HPST, heat pain supra threshold; HPT, heat pain threshold; IBS, irritable bowel syndrome; MDT, mechanical detection threshold; MPT, mechanical painful thresholdQST, quantitative sensory testing

Of QST participants with CWP questionnaire data, N=1342 had data for inflammatory markers collected within 2 years of QST visit after excluding men (n=18) and analysed in mixed effects logistic regressions (online supplemental table S1). Prevalence of CWP was 27.5%; n=117 twin pairs discordant for CWP diagnosis had inflammatory markers collected on QST visit and were examined in the sensitivity analysis (online supplemental table S2).

Of QST participants with DED questionnaire data, N=1211 had data for inflammatory markers collected within 2 years of QST visit after excluding men (n=16) and analysed in mixed effects logistic regressions (online supplemental table S3). Consistent with the main QST sample, prevalence of DED was 30.4% at sample collection. There were n=129 twin pairs discordant for DED diagnosis who had inflammatory markers collected on QST visit and included in the sensitivity analysis (online supplemental table S4).

Of QST participants with IBS questionnaire data, N=1248 had data for inflammatory markers collected within 2 years of QST visit after excluding men (n=15) and analysed in the mixed effects logistic regressions (online supplemental table S5). Prevalence of IBS was 27.0%, similar to the main QST sample. There were n=125 twin pairs discordant for IBS diagnosis with inflammatory markers collected on QST visit who were examined in the sensitivity analysis (online supplemental table S6).

In total, N=1368 unique participants had inflammatory marker data across analyses. Overlap of each analytical group and their CPS diagnoses are displayed in online supplemental figure S2. All inflammatory marker samples were above limit of detection (LOD) for IL-10, MCP-1 and TNF with n=106 samples below LOD for IL-6 and n=9 samples below LOD for IL-8. Many participants (n=1230; 89.9%) had data for inflammatory marker levels collected within a year of their QST visit, with most (n=1147; 83.8%) samples obtained on the day of QST visit.

A flowchart of all study populations is documented in online supplemental figure S3.

QST measures in CPS

We found no differences between the central tendencies of QST scores in participants with and without CWP for all 10 QST modalities (figure 1). Mann-Whitney U test p values ranged from 0.076 to 0.874 with a Bonferroni threshold of p=0.005. This finding was repeated in analyses comparing QST scores in participants with and without DED and in participants with and without IBS. Mann-Whitney U test p values in these CPS ranged from 0.135 to 0.994 and 0.077 to 0.773, respectively. Minimal detectable effect sizes (Cohen’s d) with 80.0% power for each test are found in online supplemental table S7.48 Of the total 30 Mann-Whitney U tests, comparison groups did not meet the unequal variances assumption in nine tests (online supplemental table S7). Thus, further inference about differences between medians in these tests cannot be made.

Figure 1 Heatmap of p values from Mann-Whitney U tests comparing QST scores in participants with and without CWP, DED and IBS. Each cell represents the p value of an individual Mann-Whitney U test for the corresponding QST in the relevant CPS questionnaire population (ie, p value for Mann-Whitney U test comparing CIT scores in participants with CWP and participants without CWP=0.076). Bonferroni-corrected p value threshold=0.005. CIT, cold intolerable threshold; CPT, cold pain threshold; CWP, chronic widespread pain; DED, dry eye disease; HPST, heat pain supra threshold; HPT, heat pain threshold; IBS, irritable bowel syndrome; MDT, mechanical detection threshold; MPT, mechanical pain threshold; QST, quantitative sensory testing.

Sensitivity analyses comparing QST scores of participants with CWP and true controls were consistent with main analyses and not statistically significant (online supplemental figure S4). Comparisons of QST scores in participants with prevalent CWP and participants with incident CWP were also not statistically significant (online supplemental figure S5). Mixed effects regression analyses of QST on CWP, adjusted for twin relatedness, age (scaled) and BMI category (nominal), were also consistent with the main Mann-Whitney U findings and failed to reach statistical significance (online supplemental table S8). These findings were repeated in DED and IBS analyses (online supplemental figures S4 and S5, online supplemental table S8).

Inflammation markers in CPS

In the CWP case–control mixed effects logistic regressions of inflammatory marker levels, no association reached statistical significance after Bonferroni-correction at p=0.003 (table 5). The association between IL-6 and CWP was nominally significant in a univariate model with an OR of 1.31 (95% CI 1.03 to 1.66) and p value of 0.030. All associations, however, were null after adjusting for age (scaled), BMI category (nominal) and twin relatedness.

Table 5 Mixed effects logistic regressions of inflammatory marker level on CWP diagnosis

Inflammatory marker level (NPX)	Univariate model*	Multivariate model†	
OR (95% CI)	P value	OR (95% CI)	P value	
IL-6	1.31 (1.03 to 1.66)	0.030	1.06 (0.82 to 1.38)	0.652	
IL-8	1.28 (0.99 to 1.65)	0.055	1.17 (0.90 to 1.51)	0.237	
IL-10	0.92 (0.69 to 1.22)	0.563	0.86 (0.64 to 1.15)	0.305	
MCP-1	1.05 (0.76 to 1.47)	0.759	0.89 (0.63 to 1.25)	0.488	
TNF	1.11 (0.82 to 1.50)	0.509	0.94 (0.68 to 1.29)	0.692	
Bonferroni-corrected p- value threshold=0.003.

* Random effect: family ID; fixed effects: none.

† Random effect: family ID; fixed effects: age (scaled), BMI category (nominal).

BMIbody mass indexCWP, chronic widespread pain; IL-6, interleukin-6; IL-8, interleukin-8; IL-10, interleukin-10; MCP-1, monocyte chemoattractant protein-1; NPX, Normalised Protein eXpression; TNF, tumour necrosis factor

We found no association between any inflammatory marker and DED diagnosis in mixed effects logistic regressions (table 6). ORs for all inflammatory markers approximated 1.00 with p values ranging from 0.411 to 0.775.

Table 6 Mixed effects logistic regressions of inflammatory marker level on DED diagnosis

Inflammatory marker level (NPX)	Univariate model*	Multivariate model†	
OR (95% CI)	P value	OR (95% CI)	P value	
IL-6	1.11 (0.90 to 1.37)	0.336	1.07 (0.85 to 1.34)	0.584	
IL-8	1.03 (0.82 to 1.28)	0.815	0.97 (0.77 to 1.22)	0.775	
IL-10	1.08 (0.85 to 1.37)	0.528	1.06 (0.83 to 1.35)	0.654	
MCP-1	1.15 (0.86 to 1.54)	0.342	1.06 (0.78 to 1.44)	0.701	
TNF	1.21 (0.93 to 1.57)	0.146	1.12 (0.85 to 1.47)	0.411	
Bonferroni-corrected p- value threshold=0.003.

* Random effect: family ID; fixed effects: none.

† Random effect: family ID; fixed effects: age (scaled), BMI category (nominal).

BMIbody mass indexDED, dry eye disease; IL-6, interleukin-6; IL-8, interleukin-8; IL-10, interleukin-10; MCP-1, monocyte chemoattractant protein-1; NPX, Normalised Protein eXpression; TNF, tumour necrosis factor

In the mixed effects logistic regressions of inflammatory marker levels on IBS diagnosis, the association between IL-8 and IBS diagnosis was nominally significant with an OR of 1.29 (95% CI 1.02 to 1.64) and p value of 0.036 (table 7). With a p value threshold of p=0.003, no association reached statistical significance with all other p values ranging from 0.249 to 0.893.

Table 7 Mixed effects logistic regressions of inflammatory marker level on IBS diagnosis

Inflammatory marker level (NPX)	Univariate model*	Multivariate model†	
OR (95% CI)	P value	OR (95% CI)	P value	
IL-6	0.94 (0.74 to 1.18)	0.580	0.98 (0.77 to 1.26)	0.893	
IL-8	1.20 (0.95 to 1.52)	0.124	1.29 (1.02 to 1.64)	0.036	
IL-10	1.13 (0.88 to 1.45)	0.349	1.15 (0.89 to 1.48)	0.283	
MCP-1	1.00 (0.74 to 1.35)	0.988	1.11 (0.81 to 1.51)	0.522	
TNF	1.07 (0.81 to 1.42)	0.632	1.18 (0.89 to 1.58)	0.249	
Bonferroni-corrected p- value threshold=0.003.

* Random effect: family ID; fixed effects: none.

† Random effect: family ID; fixed effects: age (scaled), BMI category (nominal).

BMIbody mass indexIBS, irritable bowel syndrome; IL-6, interleukin-6; IL-8, interleukin-8; IL-10, interleukin-10; MCP-1, monocyte chemoattractant protein-1; NPX, Normalised Protein eXpression; TNF, tumour necrosis factor

In the discordant twin sensitivity analyses, no associations were detected between intrapair differences in inflammatory marker level and CWP, in agreement with the main analyses (Bonferroni correction at p=0.003; online supplemental table S9). These findings were also repeated in DED and IBS.

Full results of inflammatory marker mixed effects regression analyses are found in online supplemental table S10. Discordant twin analyses can be found in online supplemental table S11.

Discussion

Considerations for QST interpretation in CPS

This study is the first large-scale investigation of the association of individual QST modalities with CPS in a population-based cohort. With high heterogeneity in both presentation of pain and response to common treatments, identifying a patient’s pain phenotype, and optimal treatment, is a necessary next step to improving clinical care.2 QST is an appealing tool to assist in this endeavour because it is a semi-objective, quantitative method to potentially characterise pain.49

Despite its appeal, there is an ongoing debate on how QST should be used in the clinic. While larger university hospitals have the resources to perform multiple QST modalities on patients, many clinical settings are limited to use of one or two QST modalities to measure somatosensory function, due to the expensive equipment and the highly specialised training required for QST implementation.23 50 Notably, in our cohort, no single QST modality was able to distinguish between participants with and without CWP diagnosis, DED diagnosis or IBS diagnosis. This was true with both Mann-Whitney U tests and mixed effects logistic regressions, adjusted for twin relatedness, age and BMI category. We also found no difference between QST scores in participants with prevalent CPS diagnoses at the time of QST measures and participants who were diagnosed with incident CPS later, suggesting that temporality between diagnosis and QST does not impact this outcome. This is in line with previous literature determining a lack of association between QST and migraine diagnosis; migraine, while not part of the genetic CPS cluster, is considered a common overlapping condition.51 52 In a small subset of our sample, we reported associations between presence of DED pain symptoms and heat QST modalities (HPT, HPST); however, this study also did not find significant differences in HPT or HPST between participants with and without a DED diagnosis.53 Thus, while the presence of certain subsets of pain symptoms may be associated with specific QST modalities, the null associations in the present study suggest that single QST modalities are unable to capture the heterogeneity of CPS phenotypes. This highlights the need for careful interpretation of existing QST data in CPS patients and clarification of the utility and limits of QST prior to clinical implementation that requires further exploration in future studies.

One of the strengths of our study is the large participant sample size who undertook QST measures. Pain thresholds for heat stimuli were determined in approximately 3000 participants, while pain thresholds for mechanical stimuli were determined in approximately 380 participants. QST studies are typically conducted in patient cohorts with less than 100 controls. In addition, our participants were sampled from a well-characterised cohort demonstrated to resemble an age-matched, population-based British cohort.31 No association with CPS status was detected with minimal detectable effect sizes of 0.163–0.186 at 80.0% power (1-β) in the heat modalities and 0.421–0.456 in the mechanical modalities (online supplemental table S7). If these associations do exist, they are likely to be small.

Individual inflammatory markers in CPS

This study is, to our knowledge, the largest analysis of IL-6, IL-8, IL-10, MCP-1 and TNF levels in participants with CPS. Selected a priori according to current literature, no inflammatory markers were significantly associated with CPS diagnosis in the case–control mixed effects analysis following adjustment for age, BMI category and twin relatedness; similar results were obtained in the sensitivity analysis in discordant twin pairs. This consistency of results across analyses is significant, considering the advantages of the discordant twin design—primarily the inherent matching for age, genotype (totally for MZ twins, partially for DZ twins) and most socioeconomic and environmental factors across comparison groups without additional adjustment.54

A recent systematic review and meta-analysis of 29 studies (N=2458) reported significant increases of TNF, IL-6, IL-8 and IL-10 in CWP/fibromyalgia patients compared with healthy controls.55 The component studies of this review paint a more complex picture—some reported significant increases in inflammatory marker levels, but others reported significant decreases or no significant differences. Many of the studies did not apply appropriate multiple testing corrections or adjust for age and BMI in their analyses.5559 Our results add to previous research by addressing this limitation and suggest significant increases in levels of candidate markers in CWP patients may be attributable to the inflammation associated with age and BMI than to CWP (online supplemental table S10). A similar review in IBS has pointed to the large overlap of IL-6, IL-8, IL-10 and TNF levels between patients and healthy controls in numerous studies, despite meta-analytic reports of cytokine imbalance.60 Studies included in the IBS meta-analysis also did not adjust for age or BMI. Associations between candidate cytokines and DED have primarily been derived from tear samples.27 36 The failure to replicate these associations may be due to our analysis being conducted in blood serum samples when DED is a tear and ocular surface disease. The systemic inflammation postulated to play a role in CPS may be so low, it is not reflected in levels of individual markers. Future studies examining metabolomic pathway analyses may demonstrate inflammatory pathways are over-represented in CPS patients compared with controls.

The absence of differences in inflammatory marker levels between participants with CPS and control participants does not necessarily indicate their absolute inability to be used for phenotyping purposes. Specific subtypes of each CPS reportedly have strong associations with candidate cytokines when compared with healthy controls.34 35 For example, one report noted MCP-1 was not elevated in IBS patients compared with controls, but levels were significantly higher in IBS patients with metabolic syndrome than controls.35 Our study sample may have had an under-representation of these phenotypes and an over-representation of other phenotypes unassociated with our selected cytokines. This could potentially explain opposing concentration trends and overlapping cytokine levels seen in CPS patients when compared with healthy controls in current literature and should be further explored in future studies.

We recognise the limitations in our study. First, common dynamic QST modalities were not included in our protocol, and larger QST sample size was restricted to static heat and mechanical modalities. Compared with minimal detectable effect sizes of 0.163–0.186 at 80.0% power (1-β) in heat tests, we were only powered to detect effect sizes of 0.802–0.941 in the dynamic thermal burn tests (online supplemental table S7). This was unavoidable as a secondary data analysis. Our conclusions, therefore, draw primarily from heat and mechanical static tests; other QST results must be interpreted more cautiously. Further population studies with dynamic QST modalities and increased sample size may indicate stronger associations in CPS.

We recognise that there is a large variation in the number for participants completing each QST. As a secondary data analysis, we were unable improve sample sizes. Our sample was primarily restricted to women, with men being excluded from inflammatory marker analyses, meaning results cannot be generalised to men.

Not all participants received QST and had serum samples collected on the same day. Given the agreement between results of the main analyses and temporally restricted sensitivity analyses, we believe comparison between the QST and inflammation analyses are viable.

Perhaps our greatest limitation is that participants with common painful conditions, beyond CPS (ie, osteoarthritis), were not excluded from the analysis. While CPS case status was determined through validated diagnostic questionnaires, comorbid pain conditions have the potential to influence pain sensitivity and inflammation levels.

Our findings have several implications. We found no association between single QST and CPS in a large cross-sectional analysis of over 3000 adult volunteers. Despite using a highly sensitive proteomic assay from Olink, we did not detect association between individual circulating inflammatory markers and CPS. The lack of associations demonstrates limitations of both approaches in CPS.

supplementary material

10.1136/bmjopen-2024-085814 online supplemental file 1

Data availability statement

Data are available upon reasonable request at https://twinsuk.ac.uk/resources-for-researchers/access-our-data/.

Review Process File
3 9 2024

Funding: This research was funded in whole, or in part, by the CORE Wellcome grant [WT 212904/Z/18/Z]. For the purpose of Open Access, the author has applied a CC BY public copyright license to any Author Accepted Manuscript version arising from this submission.

Prepublication history and additional supplemental material for this paper are available online. To view these files, please visit the journal online (https://doi.org/10.1136/bmjopen-2024-085814).

Provenance and peer review: Not commissioned; externally peer-reviewed.

Patient consent for publication: Not applicable.

Ethics approval: This study was carried out under TwinsUK BioBank ethics, approved by North West—Liverpool Central Research Ethics Committee (REC reference 19/NW/0187), IRAS ID 258513. Participants gave informed consent to participate in the study before taking part.

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.
==== Refs
References

1 GBD 2019 Diseases and Injuries Collaborators Global burden of 369 diseases and injuries in 204 countries and territories, 1990-2019: a systematic analysis for the Global Burden of Disease Study 2019 Lancet 2020 396 1204 22 10.1016/S0140-6736(20)30925-9 33069326
2 Edwards RR Dworkin RH Turk DC et al Patient phenotyping in clinical trials of chronic pain treatments: IMMPACT recommendations Pain 2016 157 1851 71 10.1097/j.pain.0000000000000602 27152687
3 Vehof J Zavos HMS Lachance G et al Shared genetic factors underlie chronic pain syndromes Pain 2014 155 1562 8 10.1016/j.pain.2014.05.002 24879916
4 Mansfield KE Sim J Jordan JL et al A systematic review and meta-analysis of the prevalence of chronic widespread pain in the general population Pain 2016 157 55 64 10.1097/j.pain.0000000000000314 26270591
5 Canavan C West J Card T The epidemiology of irritable bowel syndrome Clin Epidemiol 2014 6 71 80 10.2147/CLEP.S40245 24523597
6 Stapleton F Alves M Bunya VY et al TFOS DEWS II Epidemiology Report Ocul Surf 2017 15 334 65 10.1016/j.jtos.2017.05.003 28736337
7 Yang TY Chen CS Lin CL et al Risk for Irritable Bowel Syndrome in Fibromyalgia Patients: A National Database Study Medicine (Baltimore) 2017 96 e6657 10.1097/MD.0000000000006657 28383443
8 Lee CJ Levitt RC Felix ER et al Evidence that dry eye is a comorbid pain condition in a U.S. veteran population Pain Rep 2017 2 e629 10.1097/PR9.0000000000000629 29392243
9 Asproudis I Tsoumani AT Katsanos KH et al Irritable bowel syndrome might be associated with dry eye disease Ann Gastroenterol 2016 29 487 91 10.20524/aog.2016.0064 27708515
10 Fitzcharles M-A Cohen SP Clauw DJ et al Nociplastic pain: towards an understanding of prevalent pain conditions Lancet 2021 397 2098 110 10.1016/S0140-6736(21)00392-5 34062144
11 Arendt-Nielsen L Yarnitsky D Experimental and clinical applications of quantitative sensory testing applied to skin, muscles and viscera J Pain 2009 10 556 72 10.1016/j.jpain.2009.02.002 19380256
12 Desmeules JA Cedraschi C Rapiti E et al Neurophysiologic evidence for a central sensitization in patients with fibromyalgia Arthritis Rheum 2003 48 1420 9 10.1002/art.10893 12746916
13 Baron R Binder A Wasner G Neuropathic pain: diagnosis, pathophysiological mechanisms, and treatment Lancet Neurol 2010 9 807 19 10.1016/S1474-4422(10)70143-5 20650402
14 Rolke R Baron R Maier C et al Quantitative sensory testing in the German Research Network on Neuropathic Pain (DFNS): standardized protocol and reference values Pain 2006 123 231 43 10.1016/j.pain.2006.01.041 16697110
15 Cruz-Almeida Y Fillingim RB Can quantitative sensory testing move us closer to mechanism-based pain management? Pain Med 2014 15 61 72 10.1111/pme.12230 24010588
16 Nijs J Lahousse A Kapreli E et al Nociplastic Pain Criteria or Recognition of Central Sensitization? Pain Phenotyping in the Past, Present and Future J Clin Med 2021 10 3203 10.3390/jcm10153203 34361986
17 Treede R-D Jensen TS Campbell JN et al Neuropathic pain Neurol (ECronicon) 2008 70 1630 5 10.1212/01.wnl.0000282763.29778.59
18 IASP Classification of chronic pain 2nd edn Seattle, WA IASP Press 1994
19 Baron R Maier C Attal N et al Peripheral neuropathic pain: a mechanism-related organizing principle based on sensory profiles Pain 2017 158 261 72 10.1097/j.pain.0000000000000753 27893485
20 Schmelz M What can we learn from the failure of quantitative sensory testing? Pain 2021 162 663 4 10.1097/j.pain.0000000000002059 32868753
21 Vollert J Forstenpointner J Enax-Krumova EK et al The need for previous knowledge does not render quantitative sensory testing a “failure” but part of a larger picture of the relationship between nociception and pain Pain 2021 162 1273 4 10.1097/j.pain.0000000000002138 33730009
22 Kosek E Clauw D Nijs J et al Chronic nociplastic pain affecting the musculoskeletal system: clinical criteria and grading system Pain 2021 162 2629 34 10.1097/j.pain.0000000000002324 33974577
23 Cruccu G Anand P Attal N et al EFNS guidelines on neuropathic pain assessment Eur J Neurol 2004 11 153 62 10.1111/j.1468-1331.2004.00791.x 15009162
24 Mobasheri A Saarakkala S Finnilä M et al Recent advances in understanding the phenotypes of osteoarthritis F1000Res 2019 8 F1000 10.12688/f1000research.20575.1
25 Bäckryd E Tanum L Lind A-L et al Evidence of both systemic inflammation and neuroinflammation in fibromyalgia patients, as assessed by a multiplex protein panel applied to the cerebrospinal fluid and to plasma J Pain Res 2017 10 515 25 10.2147/JPR.S128508 28424559
26 Theoharides TC Tsilioni I Bawazeer M Mast Cells, Neuroinflammation and Pain in Fibromyalgia Syndrome Front Cell Neurosci 2019 13 353 10.3389/fncel.2019.00353 31427928
27 Roda M Corazza I Bacchi Reggiani ML et al Dry Eye Disease and Tear Cytokine Levels-A Meta-Analysis Int J Mol Sci 2020 21 3111 10.3390/ijms21093111 32354090
28 Bashashati M Moradi M Sarosiek I Interleukin-6 in irritable bowel syndrome: A systematic review and meta-analysis of IL-6 (-G174C) and circulating IL-6 levels Cytokine 2017 99 132 8 10.1016/j.cyto.2017.08.017 28886490
29 Verdi S Abbasian G Bowyer RCE et al TwinsUK: The UK Adult Twin Registry Update Twin Res Hum Genet 2019 22 523 9 10.1017/thg.2019.65 31526404
30 Williams FMK Scollen S Cao D et al Genes contributing to pain sensitivity in the normal population: an exome sequencing study PLoS Genet 2012 8 e1003095 10.1371/journal.pgen.1003095 23284290
31 Andrew T Hart DJ Snieder H et al Are twins and singletons comparable? A study of disease-related and lifestyle characteristics in adult women Twin Res 2001 4 464 77 10.1375/1369052012803 11780939
32 Norbury TA MacGregor AJ Urwin J et al Heritability of responses to painful stimuli in women: a classical twin study Brain (Bacau) 2007 130 3041 9 10.1093/brain/awm233
33 Ang DC Moore MN Hilligoss J et al MCP-1 and IL-8 as pain biomarkers in fibromyalgia: a pilot study Pain Med 2011 12 1154 61 10.1111/j.1526-4637.2011.01179.x 21714842
34 Bashashati M Rezaei N Shafieyoun A et al Cytokine imbalance in irritable bowel syndrome: a systematic review and meta-analysis Neurogastroenterol Motil 2014 26 1036 48 10.1111/nmo.12358 24796536
35 Chira A Braicu C Budisan L et al Monocyte chemotactic protein-1 and nitrotyrosine in irritable bowel syndrome Exp Ther Med 2020 20 24 30 10.3892/etm.2020.8665
36 Na K-S Mok J-W Kim JY et al Correlations between tear cytokines, chemokines, and soluble receptors and clinical severity of dry eye disease Invest Ophthalmol Vis Sci 2012 53 5443 50 10.1167/iovs.11-9417 22789923
37 Uçeyler N Häuser W Sommer C Systematic review with meta-analysis: cytokines in fibromyalgia syndrome BMC Musculoskelet Disord 2011 12 245 10.1186/1471-2474-12-245 22034969
38 Kumar S Shukla R Ranjan P et al Interleukin-10: A Compelling Therapeutic Target in Patients With Irritable Bowel Syndrome Clin Ther 2017 39 632 43 10.1016/j.clinthera.2017.01.030 28237672
39 Olink Olink target 96 inflammation panels n.d. Available https://www.olink.com/products-services/target/inflammation/
40 Olink What is NPX? n.d. Available https://www.olink.com/faq/what-is-npx/
41 Olink Data normalization n.d. Available https://olink.com/faq/data-normalization/
42 Cederroth CR Hong M-G Freydin MB et al Screening for Circulating Inflammatory Proteins Does Not Reveal Plasma Biomarkers of Constant Tinnitus J Assoc Res Otolaryngol 2023 24 593 606 10.1007/s10162-023-00920-3 38079022
43 White KP Harth M Speechley M et al Testing an instrument to screen for fibromyalgia syndrome in general population studies: the London Fibromyalgia Epidemiology Study Screening Questionnaire J Rheumatol 1999 26 880 4 10229410
44 Gulati A Sullivan R Buring JE et al Validation and repeatability of a short questionnaire for dry eye syndrome Am J Ophthalmol 2006 142 125 31 10.1016/j.ajo.2006.02.038 16815260
45 Bates D Mächler M Bolker B et al Fitting Linear Mixed-Effects Models Using lme4 J Stat Softw 2015 67 1 48 10.18637/jss.v067.i01
46 Zierle-Ghosh A Jan A Physiology, body mass index StatPearls [Internet] Treasure Island (FL) StatPearls Publishing 2023 Available https://www.ncbi.nlm.nih.gov/books/NBK535456/
47 Therneau T Grambsch P Modeling Survival Data: Extending the Cox Model Springer 2000
48 Faul F Erdfelder E Buchner A et al Statistical power analyses using G*Power 3.1: tests for correlation and regression analyses Behav Res Methods 2009 41 1149 60 10.3758/BRM.41.4.1149 19897823
49 Pavlaković G Petzke F The role of quantitative sensory testing in the evaluation of musculoskeletal pain conditions Curr Rheumatol Rep 2010 12 455 61 10.1007/s11926-010-0131-0 20857243
50 Zhu GC Böttger K Slater H et al Concurrent validity of a low-cost and time-efficient clinical sensory test battery to evaluate somatosensory dysfunction Eur J Pain 2019 23 1826 38 10.1002/ejp.1456 31325385
51 Krimmel SR Keaser ML Speis D et al Migraine disability, pain catastrophizing, and headache severity are associated with evoked pain and targeted by mind-body therapy Pain 2022 163 e1030 7 10.1097/j.pain.0000000000002578 35297801
52 Nahman-Averbuch H Shefi T Schneider VJ et al Quantitative sensory testing in patients with migraine: a systematic review and meta-analysis Pain 2018 159 1202 23 10.1097/j.pain.0000000000001231 29781957
53 Vehof J Kozareva D Hysi PG et al Relationship between dry eye symptoms and pain sensitivity JAMA Ophthalmol 2013 131 1304 8 10.1001/jamaophthalmol.2013.4399 23907167
54 Fagnani C Bellani M Soares JC et al Discordant twins as a tool to unravel the aetiology of bipolar disorder Epidemiol Psychiatr Sci 2014 23 137 40 10.1017/S2045796014000055 24513394
55 O’Mahony LF Srivastava A Mehta P et al Is fibromyalgia associated with a unique cytokine profile? A systematic review and meta-analysis Rheumatology (Oxford) 2021 60 2602 14 10.1093/rheumatology/keab146
56 Gür A Karakoç M Nas K et al Cytokines and depression in cases with fibromyalgia J Rheumatol 2002 29 358 61 11838856
57 Bote ME García JJ Hinchado MD et al Inflammatory/stress feedback dysregulation in women with fibromyalgia Neuroimmunomodulation 2012 19 343 51 10.1159/000341664 22986514
58 Bote ME Garcia JJ Hinchado MD et al Fibromyalgia: anti-inflammatory and stress responses after acute moderate exercise PLoS One 2013 8 e74524 10.1371/journal.pone.0074524 24023948
59 Mendieta D De la Cruz-Aguilera DL Barrera-Villalpando MI et al IL-8 and IL-6 primarily mediate the inflammatory response in fibromyalgia patients J Neuroimmunol 2016 290 22 5 10.1016/j.jneuroim.2015.11.011
60 Öhman L Törnblom H Simrén M Crosstalk at the mucosal border: importance of the gut microenvironment in IBS Nat Rev Gastroenterol Hepatol 2015 12 36 49 10.1038/nrgastro.2014.200 25446728
