Relationship between major depressive disorder and [ME/CFS]: a two-sample mendelian randomization study analysis, 2025, Zhu et al

On the other hand the authors haven't published before on ME/CFS or mendelian randomisation from what I can see while the senior author has published a lot on acupuncture.

I did find that the lead author, Wenjing Song, has done at least two mendelian randomization studies before.

For Daytime napping and the incidence of Parkinson’s disease: a prospective cohort study with Mendelian randomization, I looked at the SNPs in the paper, and compared them to the SNPs in the 23andMe paper they cite on daytime napping, and most of them match.

For the other, I didn't see any SNPs listed.
 
Last edited:
Seriously impressive digging again on here.

I’m afraid I can’t help, to pass on a couple of things Chris has mentioned to me in general. First, because UK buying bank data is so accessible, it’s used in a lot of GWAS Studies, and some of these are of questionable quality. Second, Mendelian randomisation studies are sophisticated and hard to get right.

These are general comments and don’t tell us anything about this study, but the apparent errors are not a great sign. Hopefully, the authors will clarify things
 
Will send a message to the authors asking for a clarification.
Did they ever respond? If not, it might be worthwhile to post something to PubPeer. I drafted this up and could post it:

There appears to be an issue with the SNPs reported in this paper in that they do not match the SNPs reported in the referenced MDD study, and curiously, every single SNP in this paper has been associated with BMI or obesity in other papers.

The paper says:
> Associations between single-nucleotide polymorphisms (SNPs) and MDD were estimated using data from the largest published GWAS meta-analysis of European ancestry to date[8]. [...] Following a meticulous reconciliation of exposure and outcome data, a subset of 57 SNPs, which exhibited the strongest associations with MDD and had the highest levels of statistical confidence, were selected as IVs for the subsequent MR analysis.

[8] is Howard, D. M. et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. _Nat. Neurosci. vol_. **22** (3), 343–352. [https://doi.org/10.1038/s41593-018-0326-7](https://doi.org/10.1038/s41593-018-0326-7) (2019).

Supplementary Material 1 has a list of 57 SNPs, which the paper says were chosen based on strong associations with MDD in Howard et al. Supplementary Table 1 from the Howard et al paper lists 102 variants that were found to be associated with MDD. There is no overlap between the SNPs reported in these papers.

All SNPs listed in Song et al have previously been reported to be associated with BMI or obesity, and can all be found mentioned as such among the following papers:

- Li, J., Tian, A., Zhu, H., Chen, L., Wen, J., Liu, W., & Chen, P. (2022). Mendelian Randomization Analysis Reveals No Causal Relationship Between Nonalcoholic Fatty Liver Disease and Severe COVID-19. _Clinical gastroenterology and hepatology : the official clinical practice journal of the American Gastroenterological Association_, _20_(7), 1553–1560.e78. https://doi.org/10.1016/j.cgh.2022.01.045
- Garfield, V., Llewellyn, C. H., Wichstrøm, L., & Steinsbekk, S. (2021). Shared genetic architecture underlying sleep and weight in children. _Sleep medicine_, _83_, 40–44. https://doi.org/10.1016/j.sleep.2021.04.021
- Shimaoka, I., Kamide, K., Ohishi, M., Katsuya, T., Akasaka, H., Saitoh, S., Sugimoto, K., Oguro, R., Congrains, A., Fujisawa, T., Shimamoto, K., Ogihara, T., & Rakugi, H. (2010). Association of gene polymorphism of the fat-mass and obesity-associated gene with insulin resistance in Japanese. _Hypertension research : official journal of the Japanese Society of Hypertension_, _33_(3), 214–218. https://doi.org/10.1038/hr.2009.215
- Fall, T., Hägg, S., Ploner, A., Mägi, R., Fischer, K., Draisma, H. H., Sarin, A. P., Benyamin, B., Ladenvall, C., Åkerlund, M., Kals, M., Esko, T., Nelson, C. P., Kaakinen, M., Huikari, V., Mangino, M., Meirhaeghe, A., Kristiansson, K., Nuotio, M. L., Kobl, M., … ENGAGE Consortium (2015). Age- and sex-specific causal effects of adiposity on cardiovascular risk factors. _Diabetes_, _64_(5), 1841–1852. https://doi.org/10.2337/db14-0988
- Nikpay, M., Turner, A. W., & McPherson, R. (2018). Partitioning the Pleiotropy Between Coronary Artery Disease and Body Mass Index Reveals the Importance of Low Frequency Variants and Central Nervous System-Specific Functional Elements. _Circulation. Genomic and precision medicine_, _11_(2), e002050. https://doi.org/10.1161/CIRCGEN.117.002050
- He, Y., Zheng, C., He, M. H., & Huang, J. R. (2021). The Causal Relationship Between Body Mass Index and the Risk of Osteoarthritis. _International journal of general medicine_, _14_, 2227–2237. https://doi.org/10.2147/IJGM.S314180
- Beckers, S., Peeters, A., Zegers, D., Mertens, I., Van Gaal, L., & Van Hul, W. (2008). Association of the BDNF Val66Met variation with obesity in women. _Molecular genetics and metabolism_, _95_(1-2), 110–112. https://doi.org/10.1016/j.ymgme.2008.06.008
- Young, A. I., Wauthier, F. L., & Donnelly, P. (2018). Identifying loci affecting trait variability and detecting interactions in genome-wide association studies. _Nature genetics_, _50_(11), 1608–1614. https://doi.org/10.1038/s41588-018-0225-6

Edit: Or feel free to post something yourself since it was your discovery.
 
Last edited:
Did they ever respond?
Yes but not sure if I can copy their reply so I'll try to summarise.

They basically said that they did not use the 102 significant SNPs highlighted in the publication by Howard et al. (2019). Instead they screen the full dataset where summary data for all tested SNPs is listed. They don't say this explicitly but I believe this info is available in the text file named ‘PGC_UKB_depression_genome-wide.txt’ on this webpage, which the Howard 2019 publication refers to:
https://datashare.ed.ac.uk/handle/10283/3203

But even then, the data don't seem to match from what I can tell. The SNP listed in the supplementary material still do not seem to be associated with depression and the p-values don't match. To take one example: the first SNP in the supplementary material is rs1000940 for which Song et al. report a p-value (pval.exposure) of 1.63E-08. But in the text file listed by Howard et al. 2019 I found a p-value of 0.5313.
 
Yes but not sure if I can copy their reply so I'll try to summarise.

They basically said that they did not use the 102 significant SNPs highlighted in the publication by Howard et al. (2019). Instead they screen the full dataset where summary data for all tested SNPs is listed. They don't say this explicitly but I believe this info is available in the text file named ‘PGC_UKB_depression_genome-wide.txt’ on this webpage, which the Howard 2019 publication refers to:
https://datashare.ed.ac.uk/handle/10283/3203

But even then, the data don't seem to match from what I can tell. The SNP listed in the supplementary material still do not seem to be associated with depression and the p-values don't match. To take one example: the first SNP in the supplementary material is rs1000940 for which Song et al. report a p-value (pval.exposure) of 1.63E-08. But in the text file listed by Howard et al. 2019 I found a p-value of 0.5313.
I'm concerned by their response. It sounds like they're saying they independently found 57 new SNPs associated with MDD from a GWAS already analyzed by others, which sounds like a study worth a paper of its own, but its just a small section of another study that then uses those SNPs in mendelian randomization.

I ran a script to get all 57 rows from that PGC text file that match Song et al's SNPs:
MarkerName P
rs3101336 6.242E-18
rs2365389 0.00000007333
rs17724992 0.0005739
rs12986742 0.0009833
rs3736485 0.0009973
rs3817334 0.001145
rs2049045 0.001492
rs663129 0.001656
rs13078960 0.004567
rs1167827 0.009531
rs6804842 0.01044
rs7903146 0.0117
rs29941 0.0138
rs1528435 0.01902
rs1558902 0.01961
rs1928295 0.02317
rs879620 0.0233
rs3845344 0.02551
rs16951275 0.03372
rs12940622 0.03775
rs9926784 0.03867
rs4740619 0.06251
rs758747 0.06652
rs10182181 0.08368
rs2820292 0.09436
rs657452 0.1396
rs7138803 0.1658
rs891389 0.2129
rs10938397 0.2183
rs10733682 0.2633
rs2303223 0.2978
rs7599312 0.3169
rs1412235 0.3248
rs2112347 0.3515
rs7144011 0.3609
rs9400239 0.3628
rs2033732 0.3777
rs13021737 0.386
rs12286929 0.3861
rs2207139 0.3925
rs11165643 0.431
rs6457796 0.4417
rs11672660 0.4682
rs9579083 0.4935
rs1000940 0.5313
rs2245368 0.5447
rs10132280 0.5518
rs17094222 0.5602
rs1016287 0.6559
rs6477694 0.7286
rs10840100 0.738
rs3849570 0.7559
rs9304665 0.7952
rs543874 0.8204
rs17405819 0.8423
rs3888190 0.9293
rs2176598 0.9732
 
Back
Top Bottom