Analysis of ME/CFS Patient Clusters in Vaes Data Reveals Distinct Subgroups

1. Introduction

Dr. Anouk W. Vaes and her collaborators at CIRO performed one of the largest systematic surveys of ME/CFS symptom patterns, applying clustering methods to patient-reported outcomes to identify recurring constellations of symptoms [1]. They used symptom surveys based on the DePaul Symptom Questionnaire version 2 (DSQ-2) from 337 patients to identify a total of 45 patient clusters of which 13 were of size ≥ 10. We used the publicly available data of their final 13 clusters for analysis [2].

While Vaes described each cluster, and noted some differences among symptom severity, an overall organizing principle was not proposed. Other methods have been used to group ME/CFS phenotypes, including by severity alone, or a combination of severity and symptoms but no consensus has emerged on how to best group patients.

Our analysis suggests that the Vaes ME/CFS clusters can be organized into two overarching families defined by their relationship to PEM. This reduces the complexity of the cluster symptom intensities and highlights a fibromyalgia-like subgroup within ME/CFS that has not been explicitly recognized in prior work.

Our use of the terms “high” and “low intensity” does not reproduce the clusters described by Vaes or others but reflects a new organization based on how symptom domains scale with PEM. Vaes used the DePaul symptom groupings which we keep as-is.

2. Symptom Domains vs. PEM

The strong weight of PEM symptoms makes the Vaes dataset tricky to interpret if you only look at raw intensity. We compared overall symptom levels to PEM severity and grouped clusters by how those two measures relate.

In the following sections we will start by showing how the groups are most similar, then how they differ. We add an “all” category to gauge overall symptom severity to the original symptom domains.

Notes: The DSQ-2 uses a single symptom for fatigue. Also, for the “Overall” chart we used the arithmetic mean instead of geometric mean for the Y axis. All other graphs in this study use geometric means exclusively.

2.1. Similarities

In terms of fatigue, the two groups overlap completely, as ME/CFS patient clusters are expected to.

While fatigue (Figure 1) shows complete overlap in our groups the next charts (Figure 2) show how they stand out. In the first chart (top left) we compare the average of all symptoms to PEM. In this chart one can easily discern two tiers of clusters. One has a higher overall symptom burden at any given PEM level compared to the other and this relationship is where we derive the names for our groups: high and low intensity.

Autonomic, Neuroendocrine and Other also follow a similar pattern: Parallel but lower than the high intensity group.

We note that the limited range of PEM severity in the low intensity group (~2.2 to 2.9) may have otherwise caused it to remain undetected.

Groupings: From visual inspection of the charts above we group clusters C9, C19, C26, C28, C31, C36, C37, C40 as high intensity, C2, C4, C7, C11, C24 as low intensity. We’ll keep these groupings for all plots that follow.

The last chart (Other vs. PEM) shows a markedly lower offset and is perhaps the most visibly obvious difference between the two groups. The two symptoms in `other’ are sensitivity to mold and vibration.

While the groups maintain integrity, the parallelism observed in the charts above does not hold in the next charts and suggests a more complex physiological cause for the differences.

2.2. Increasing Amplification vs. PEM

Pain and neurocognitive domains show that while the low intensity group remains offset, symptom domain amplification is increased. As PEM increases to 3, C7 comes very close to the high intensity line.

Figure 3. Pain and Neurocognitive vs. PEM

Figure 3. Pain and Neurocognitive vs. PEM

The two charts above, while not severe enough to claim this is fibromyalgia with PEM, begin to suggest it.

2.3. Different Correlations to PEM

In all of the previous charts we’ve shown that both groups maintain a strong but distinct relationship to PEM. By contrast, the immune and sleep domains show little to no correlation with PEM severity in the low intensity group. Notably, despite this lack of correlation, sleep disturbance can be pronounced in the low intensity group.

Figure 4. Immune and Sleep vs. PEM

Figure 4. Immune and Sleep vs. PEM

3. Summary of Differences

High Intensity: Has a consistent amplification pattern: as PEM increases, all other symptom domains rise together.

Low Intensity: Presents a more selective profile. Pain, neurocognitive, and neuroendocrine symptoms escalate more rapidly with PEM, while sleep and immune symptoms show little or no correlation. The elevated neurocognitive and pain symptoms, together with variable sleep symptom intensity makes the low-intensity group more suggestive of fibromyalgia than the high-intensity group. Additionally, the low intensity group is notably nestled in the middle of the overall PEM range which could allow it to hide among the other clusters if not actively searched for.

We summarize these differences in the table below.

Table 1. Summary of Group Differences.

Feature	High group	Low group
Overall vs. PEM	Higher	Lower - esp. “other”
Fatigue	Same	Same
PEM range	~1.3 – 3.4	Narrower ~2.2 – 2.9
Neurocognitive	Tracks PEM	Elevated amplification
Pain	Tracks PEM	Elevated amplification
Immune	Tracks PEM	Flat at ~1
Sleep	Tracks PEM	~1.4 – 2.5, but uncorrelated [1]

Table 1. Summary of Group Differences.

Feature	High group	Low group
Overall vs. PEM	Higher	Lower - esp. “other”
Fatigue	Same	Same
PEM range	~1.3 – 3.4	Narrower ~2.2 – 2.9
Neurocognitive	Tracks PEM	Elevated amplification
Pain	Tracks PEM	Elevated amplification
Immune	Tracks PEM	Flat at ~1
Sleep	Tracks PEM	~1.4 – 2.5, but uncorrelated [1]

Symptom Severity: In terms of overall symptom range the two groups largely overlap. The high intensity group, by symptom intensity alone, could be considered a superset of the low intensity group. It is the differing relationships to PEM and the selective amplification of certain symptom domains in the low intensity group that set the two groups apart.

4. Discussion

We have shown that the Vaes clusters can be organized into two subgroups defined by symptom domains relative to PEM, and that these subgroups remain consistent across most symptom domains.

ME/CFS is widely recognized as a heterogeneous condition with complex symptom presentations. Grouping the Vaes clusters by their relationship to PEM simplifies this complexity at the domain and cluster levels, but it only partially reduces heterogeneity at the individual-symptom level and may not translate directly to grouping individual patients. Vaes [1] provided cluster-specific summaries of standout symptoms; some of those align with our domain-based charts while others show a different emphasis, underscoring the challenge of reconciling symptom-level observations with domain-level analyses.

Although these patterns may guide future phenotyping, several cautions are warranted. This is a secondary analysis of patient clusters defined by symptom domains with small subgroup sizes (8 and 5 clusters). The original clustering methods of Vaes et al. may have obscured or exaggerated some relationships we show here. Moreover, visual inspection of the plots may suggest directional trends where none truly exist.

Despite these limitations, we believe at least two further steps should be taken. Using the original patient scores could validate, refine, or repudiate these subgroups and possibly produce a short list of differentiating symptoms. If this proves useful then examining biological markers across these subgroups could help to reveal underlying pathophysiological differences.

5. Conclusions

Our secondary analysis shows that the many heterogeneous Vaes clusters can be organized into two broader structures defined by overall symptom intensity relative to PEM severity. These two groups remain coherent across symptom domains and display distinct, largely linear relationships with PEM. Although exploratory and limited to cluster-level data, these findings offer a reproducible framework that may help clarify apparent heterogeneity in ME/CFS.

Appendix A.1. Cluster Summary Data

We use the Vaes cluster spreadsheet [1] exclusively as our data source. We use R for significant cleaning and transformation from Excel to CSV. After removing everything but the mean severity of the symptoms in the 13 clusters we:

• Group each cluster by symptom domain and calculate geometric means for each symptom domain within each cluster.
• Rotate the table
• Add a new column, “all_mean” which is the arithmetic mean of all symptoms within each cluster.
• Save this file as cluster_grouped_tidy.csv

This file is what we then used for our analysis, and the plots. We present the final data used below. Values were rounded and abbreviations used for formatting. See the [sec-code]Code Availability section for links to the code used to generate this file.

Table A1. Cluster Summary Data with Geometric Means by Symptom Domain

cluster	all_mean	PEM	fatigue	Auto	Immune	NC	NE	pain	other
C2	1.21	2.30	2.79	0.76	0.53	1.15	0.67	0.95	0.27
C4	1.55	2.40	3.13	1.24	0.87	1.66	0.96	1.19	0.46
C7	1.92	2.85	3.20	1.48	1.03	2.32	1.05	1.94	0.86
C9	2.06	2.64	2.85	2.16	1.40	2.24	1.52	1.71	1.55
C11	1.42	2.23	2.94	1.02	1.08	1.11	0.64	1.00	0.33
C19	1.37	1.97	2.53	1.12	0.57	1.81	0.61	1.06	0.60
C24	1.63	2.50	2.95	1.44	0.82	1.52	0.86	1.53	0.56
C26	2.14	2.80	3.28	1.87	1.76	2.06	1.52	2.34	1.77
C28	1.03	1.62	2.50	0.83	0.56	0.78	0.60	0.88	0.30
C31	1.47	2.15	2.73	1.31	0.77	1.31	1.04	1.54	1.41
C36	2.74	3.35	3.50	2.59	2.30	3.08	2.04	2.73	2.04
C37	0.67	1.30	2.39	0.41	0.16	0.43	0.39	0.51	0.19
C40	1.21	1.76	2.80	0.74	0.52	0.95	0.82	1.19	0.67

Table A1. Cluster Summary Data with Geometric Means by Symptom Domain

cluster	all_mean	PEM	fatigue	Auto	Immune	NC	NE	pain	other
C2	1.21	2.30	2.79	0.76	0.53	1.15	0.67	0.95	0.27
C4	1.55	2.40	3.13	1.24	0.87	1.66	0.96	1.19	0.46
C7	1.92	2.85	3.20	1.48	1.03	2.32	1.05	1.94	0.86
C9	2.06	2.64	2.85	2.16	1.40	2.24	1.52	1.71	1.55
C11	1.42	2.23	2.94	1.02	1.08	1.11	0.64	1.00	0.33
C19	1.37	1.97	2.53	1.12	0.57	1.81	0.61	1.06	0.60
C24	1.63	2.50	2.95	1.44	0.82	1.52	0.86	1.53	0.56
C26	2.14	2.80	3.28	1.87	1.76	2.06	1.52	2.34	1.77
C28	1.03	1.62	2.50	0.83	0.56	0.78	0.60	0.88	0.30
C31	1.47	2.15	2.73	1.31	0.77	1.31	1.04	1.54	1.41
C36	2.74	3.35	3.50	2.59	2.30	3.08	2.04	2.73	2.04
C37	0.67	1.30	2.39	0.41	0.16	0.43	0.39	0.51	0.19
C40	1.21	1.76	2.80	0.74	0.52	0.95	0.82	1.19	0.67