chillier
Senior Member (Voting Rights)
This figure is the bit that really matters for sure. They have taken muscle biopsies from these 10 controls and 14 ME/CFS patients (from the intramural study i think so should be diagnosed well) and done western blots on them for WASF3, MTCO1, COX17, PERK and BiP. Here's the rest of the figure showing the western blots themselves and the quantified bar plots (just quantifying how large the blobs are from the western blot).
View attachment 20151
They see elevated WASF3, reduced ComplexIV proteins MTCO1 and COX17, increased endoplasmic reticulum (ER) stress marker PERK, and reduced protein-folding chaperone BiP. I really don't like when this type of data is represented in bar plots without showing the individual datapoints as a strip style scatter plot. We have no idea how much the groups overlap and they also don't tell us the effect size. Our only indication of variance is the error bar which as you point out represents the standard error which is the wrong metric to use - it should be standard deviation. The standard error of the mean is always smaller than the standard deviation so it looks better too. They do at least show all of the blots for all these proteins for the controls and patients (some hidden in the supplementary) so theoretically we could make these plots properly ourselves.
They go on to see the same set of proteins showing the same phenotype when they chemically induce ER stress in a separate experiment which is potentially interesting.
Their overall model for what is happening is that ER stress leads to an increase in PERK expression (a marker for ER stress), a reduction in BiP expression (an ER protein which helps proteins fold as they're being made) and in someway therefore an increase in WASF3 expression (being translated directly into the ER co-translationally I think). WASF3 translocates over to the mitochondria and binds to oxidative phosphorylation complex 3 preventing it from forming a complex with complex 4, thereby inhibiting oxidative phosphorylation activity.
According to the text S1 had all sorts of different diseases and that she had chronic fatigue and exercise intolerance - though no specific mention of ME/CFS diagnosis.
So yeah as always it needs to be replicated in a much bigger cohort with sedentary controls carefully selected and with disease controls - possible it could be a broadly applicable phenotype not specific to ME/CFS.
I've done the densitometry myself and quantified relative to GAPDH as they would presumably have done (a "housekeeping" gene there to normalise for total protein levels). Bear in mind I'm working with the snippets of the blot membranes they are showing in the paper and not the membrane itself, and it is split across different figures. Hopefully the contrast has been maintained between the different crops they've done. This also made it tricky to get a good read of the background signal which you need to normalise properly. That being said:

P-values (two-tailed unpaired student's t test) and effect sizes (cohen's D):
Gene: WASF3
P-value: 0.0859249
Cohen's d: 0.7596534
Gene: MTCO1
P-value: 0.0007785032
Cohen's d: 2.143763
Gene: COX17
P-value: 0.00237593
Cohen's d: 1.633564
Gene: PERK
P-value: 0.02902655
Cohen's d: 0.9279234
Gene: BiP
P-value: 0.004006675
Cohen's d: 1.298791
Weirdly WASF3 doesn't actually make the 0.05 threshold in my attempt, but it is close so I can see how it might have squeezed past the threshold when they did it. Some of the others though scream through with tiny p values and enormous effect sizes.