Combination of 15 lipid metabolites and motilin to diagnose spleen-deficiency FD

Jiaqi Zhang, Xue Wang, Xiaoshuang Shi, Jingyi Xie, Min Zhang, Jinxin Ma, Fengyun Wang, Xudong Tang, Jiaqi Zhang, Xue Wang, Xiaoshuang Shi, Jingyi Xie, Min Zhang, Jinxin Ma, Fengyun Wang, Xudong Tang

Abstract

Background: This study aims to assess clinical characteristics in FD with spleen deficiency syndrome and metabolic perturbations involved in FD progress. We combined metabolic biomarkers and clinical features into a better prediction for FD with Spleen Deficiency syndrome.

Methods: A total of 276 people were recruited, including 215 FD patients and 61 healthy control group (HC). The clinical characteristics and gastric emptying rate were compared between spleen deficiency-FD group and non-spleen deficiency-FD. The serum lipids metabonomics analysis was performed to determine the metabolic differences in spleen deficiency-FD group and HC.

Results: The symptoms of postprandial discomfort in Spleen Deficiency group were more severe (P < 0.05), and delayed gastric emptying was more pronounced (P < 0.05) vs. non-Spleen deficiency. Decreased motilin (OR = 0.990, 95% confidence interval (CI) 0.982-0.997) was independent risk factor related to Spleen Deficiency group. We identified 15 metabolites for spleen deficiency group vs. HC, majority of those biomarkers belonged to the glycerophospholipid metabolic pathway. The combination of 15 metabolics could diagnose spleen deficiency-FD, with the AUC of 0.9943, 95% CI 0.9854-1.0000), and the combination of 15 metabolics and motilin could diagnose spleen deficiency-FD, with the AUC of 0.9615, 95% CI 0.9264-9967).

Conclusions: This study provides supportive evidence that Spleen deficiency syndrome was associated with delayed gastric emptying and the glycerophospholipid metabolic pathway was perturbed in FD patients. The combination of metabolic biomarkers and clinical features provided us with new ideas for multidimensional diagnosis of FD.Trial registration http://www.chictr.org.cn, no: ChiCTR-TRC-13003200. clinicaltrials.gov, no: NCT02762136.

Keywords: Biomarker; Functional dyspepsia; Metabolomics; Spleen deficiency syndrome.

Conflict of interest statement

The authors declare that they have no competing interests.

Figures

Fig. 1
Fig. 1
The flow chart of the selection process
Fig. 2
Fig. 2
Postprandial discomfort severity scale (*p 

Fig. 3

Gastric emptying rates of proximal…

Fig. 3

Gastric emptying rates of proximal and distal stomach (*p

Fig. 3
Gastric emptying rates of proximal and distal stomach (*p 

Fig. 4

Risk factors for the FD…

Fig. 4

Risk factors for the FD patients with spleen deficiency syndrome. The results of…

Fig. 4
Risk factors for the FD patients with spleen deficiency syndrome. The results of multivariate logistic regression shown that decreased motilin was an independent risk factor related to spleen deficiency in patients with FD

Fig. 5

The PCA scores of spleen…

Fig. 5

The PCA scores of spleen deficiency-FD group and healthy control group. ( a…

Fig. 5
The PCA scores of spleen deficiency-FD group and healthy control group. (a All spleen deficiency-FD patients and healthy control group; b spleen deficiency-FD patients of XiYuan Hospital and healthy control group; c spleen deficiency-FD patients of GuangDong Province Traditional Chinese Medical Hospital and healthy control group; d spleen deficiency-FD patients of WuHan Integrated TCM &Western Medicine Hospital and healthy control group)

Fig. 6

The OPLS-DA scores of spleen…

Fig. 6

The OPLS-DA scores of spleen deficiency-FD group and healthy control group. ( a…

Fig. 6
The OPLS-DA scores of spleen deficiency-FD group and healthy control group. (a spleen deficiency-FD patients of XiYuan Hospital and healthy control group; b spleen deficiency-FD patients of GuangDong Province Traditional Chinese Medical Hospital and healthy control group; c spleen deficiency-FD patients of WuHan Integrated TCM &Western Medicine Hospital and healthy control group)

Fig. 7

The S-Plot of spleen deficiency-FD…

Fig. 7

The S-Plot of spleen deficiency-FD group and healthy control group. ( a spleen…

Fig. 7
The S-Plot of spleen deficiency-FD group and healthy control group. (a spleen deficiency-FD patients of XiYuan Hospital and healthy control group; b spleen deficiency-FD patients of GuangDong Province Traditional Chinese Medical Hospital and healthy control group; c spleen deficiency-FD patients of WuHan Integrated TCM &Western Medicine Hospital and healthy control group)

Fig. 8

Venn plot for all metabonomics…

Fig. 8

Venn plot for all metabonomics biomarkers among three groups. There are 15 overlapping…

Fig. 8
Venn plot for all metabonomics biomarkers among three groups. There are 15 overlapping metabolites: LysoPC (18:0), PC (22:6/0:0), LysoPC (17:0), LysoPC (18:1), LysoPC (20:4), PC (17:2/0:0), PC (17:1/0:0), LysoPC (16:1), LysoPE (18:2/0:0), PC (15:0/0:0), PC (16:0/0-1:0), LysoPE (0:0/16:0), LysoPE (20:4/0:0), LysoPE (0:0/22:6), LysoPC (20:3)

Fig. 9

Metabolism map of FD with…

Fig. 9

Metabolism map of FD with spleen deficiency syndrome. (1) glycerol phospholipid metabolism; (2)…

Fig. 9
Metabolism map of FD with spleen deficiency syndrome. (1) glycerol phospholipid metabolism; (2) linoleic acid metabolism; (3) sphingolipid metabolism; (4) alpha-linolenic acid metabolism; (5) glyceride metabolism; (6) arachidonic acid metabolism

Fig. 10

ROC curves of biomarkers for…

Fig. 10

ROC curves of biomarkers for Spleen-deficiency in prediction. The prediction of spleen deficiency-FD…

Fig. 10
ROC curves of biomarkers for Spleen-deficiency in prediction. The prediction of spleen deficiency-FD vs. HC was robust, and the area under the curve (AUC) of 15 biomarkers for spleen deficiency-FD vs. HC ranged from 0.6790 to 0.9926. The combination of 15 biomarkers could offer a highly accurate diagnosis of spleen deficiency in FD with an AUC of 0.9943
All figures (10)
Similar articles
Cited by
References
    1. Tack J, Talley NJ. Functional dyspepsia—symptoms, definitions and validity of the Rome III criteria. Nat Rev Gastroenterol Hepatol. 2013;10(03):134–141. doi: 10.1038/nrgastro.2013.14. - DOI - PubMed
    1. Ford AC, Luthra P, Tack J, et al. Efficacy of psychotropic drugs in functional dyspepsia: systematic review and meta-analysis. Gut. 2017;66(03):411–420. doi: 10.1136/gutjnl-2015-310721. - DOI - PubMed
    1. Wang F-Y, Ma X-X, et al. Efficacy and safety of Xiangsha Liujunzi granules for functional dyspepsia: a multi-center randomized double-blind placebo-controlled clinical study. World J Gastroenterol. 2017;23(30):5589–5601. doi: 10.3748/wjg.v23.i30.5589. - DOI - PMC - PubMed
    1. Lin L, Wang FY, Tang XD, et al. Effect of Pixu I recipe on cytochrome c oxidase subunit va expression in gastric tissues of FD rats with spleen deficiency. Chin J Exp Tradit Med Formulae. 2017;23(07):102–108.
    1. Zhong ZS, Zhang HY, Zhang W, et al. Effect of Si Junzitang on CaM-MLCK pathway in rats with spleen deficiency and gastrointestinal motility disorder. Chin J Exp Tradit Med Formulae. 2018;24(05):95–99.
Show all 43 references
Associated data
Related information
[x]
Cite
Copy Download .nbib
Format: AMA APA MLA NLM
Fig. 3
Fig. 3
Gastric emptying rates of proximal and distal stomach (*p 

Fig. 4

Risk factors for the FD…

Fig. 4

Risk factors for the FD patients with spleen deficiency syndrome. The results of…

Fig. 4
Risk factors for the FD patients with spleen deficiency syndrome. The results of multivariate logistic regression shown that decreased motilin was an independent risk factor related to spleen deficiency in patients with FD

Fig. 5

The PCA scores of spleen…

Fig. 5

The PCA scores of spleen deficiency-FD group and healthy control group. ( a…

Fig. 5
The PCA scores of spleen deficiency-FD group and healthy control group. (a All spleen deficiency-FD patients and healthy control group; b spleen deficiency-FD patients of XiYuan Hospital and healthy control group; c spleen deficiency-FD patients of GuangDong Province Traditional Chinese Medical Hospital and healthy control group; d spleen deficiency-FD patients of WuHan Integrated TCM &Western Medicine Hospital and healthy control group)

Fig. 6

The OPLS-DA scores of spleen…

Fig. 6

The OPLS-DA scores of spleen deficiency-FD group and healthy control group. ( a…

Fig. 6
The OPLS-DA scores of spleen deficiency-FD group and healthy control group. (a spleen deficiency-FD patients of XiYuan Hospital and healthy control group; b spleen deficiency-FD patients of GuangDong Province Traditional Chinese Medical Hospital and healthy control group; c spleen deficiency-FD patients of WuHan Integrated TCM &Western Medicine Hospital and healthy control group)

Fig. 7

The S-Plot of spleen deficiency-FD…

Fig. 7

The S-Plot of spleen deficiency-FD group and healthy control group. ( a spleen…

Fig. 7
The S-Plot of spleen deficiency-FD group and healthy control group. (a spleen deficiency-FD patients of XiYuan Hospital and healthy control group; b spleen deficiency-FD patients of GuangDong Province Traditional Chinese Medical Hospital and healthy control group; c spleen deficiency-FD patients of WuHan Integrated TCM &Western Medicine Hospital and healthy control group)

Fig. 8

Venn plot for all metabonomics…

Fig. 8

Venn plot for all metabonomics biomarkers among three groups. There are 15 overlapping…

Fig. 8
Venn plot for all metabonomics biomarkers among three groups. There are 15 overlapping metabolites: LysoPC (18:0), PC (22:6/0:0), LysoPC (17:0), LysoPC (18:1), LysoPC (20:4), PC (17:2/0:0), PC (17:1/0:0), LysoPC (16:1), LysoPE (18:2/0:0), PC (15:0/0:0), PC (16:0/0-1:0), LysoPE (0:0/16:0), LysoPE (20:4/0:0), LysoPE (0:0/22:6), LysoPC (20:3)

Fig. 9

Metabolism map of FD with…

Fig. 9

Metabolism map of FD with spleen deficiency syndrome. (1) glycerol phospholipid metabolism; (2)…

Fig. 9
Metabolism map of FD with spleen deficiency syndrome. (1) glycerol phospholipid metabolism; (2) linoleic acid metabolism; (3) sphingolipid metabolism; (4) alpha-linolenic acid metabolism; (5) glyceride metabolism; (6) arachidonic acid metabolism

Fig. 10

ROC curves of biomarkers for…

Fig. 10

ROC curves of biomarkers for Spleen-deficiency in prediction. The prediction of spleen deficiency-FD…

Fig. 10
ROC curves of biomarkers for Spleen-deficiency in prediction. The prediction of spleen deficiency-FD vs. HC was robust, and the area under the curve (AUC) of 15 biomarkers for spleen deficiency-FD vs. HC ranged from 0.6790 to 0.9926. The combination of 15 biomarkers could offer a highly accurate diagnosis of spleen deficiency in FD with an AUC of 0.9943
All figures (10)
Fig. 4
Fig. 4
Risk factors for the FD patients with spleen deficiency syndrome. The results of multivariate logistic regression shown that decreased motilin was an independent risk factor related to spleen deficiency in patients with FD
Fig. 5
Fig. 5
The PCA scores of spleen deficiency-FD group and healthy control group. (a All spleen deficiency-FD patients and healthy control group; b spleen deficiency-FD patients of XiYuan Hospital and healthy control group; c spleen deficiency-FD patients of GuangDong Province Traditional Chinese Medical Hospital and healthy control group; d spleen deficiency-FD patients of WuHan Integrated TCM &Western Medicine Hospital and healthy control group)
Fig. 6
Fig. 6
The OPLS-DA scores of spleen deficiency-FD group and healthy control group. (a spleen deficiency-FD patients of XiYuan Hospital and healthy control group; b spleen deficiency-FD patients of GuangDong Province Traditional Chinese Medical Hospital and healthy control group; c spleen deficiency-FD patients of WuHan Integrated TCM &Western Medicine Hospital and healthy control group)
Fig. 7
Fig. 7
The S-Plot of spleen deficiency-FD group and healthy control group. (a spleen deficiency-FD patients of XiYuan Hospital and healthy control group; b spleen deficiency-FD patients of GuangDong Province Traditional Chinese Medical Hospital and healthy control group; c spleen deficiency-FD patients of WuHan Integrated TCM &Western Medicine Hospital and healthy control group)
Fig. 8
Fig. 8
Venn plot for all metabonomics biomarkers among three groups. There are 15 overlapping metabolites: LysoPC (18:0), PC (22:6/0:0), LysoPC (17:0), LysoPC (18:1), LysoPC (20:4), PC (17:2/0:0), PC (17:1/0:0), LysoPC (16:1), LysoPE (18:2/0:0), PC (15:0/0:0), PC (16:0/0-1:0), LysoPE (0:0/16:0), LysoPE (20:4/0:0), LysoPE (0:0/22:6), LysoPC (20:3)
Fig. 9
Fig. 9
Metabolism map of FD with spleen deficiency syndrome. (1) glycerol phospholipid metabolism; (2) linoleic acid metabolism; (3) sphingolipid metabolism; (4) alpha-linolenic acid metabolism; (5) glyceride metabolism; (6) arachidonic acid metabolism
Fig. 10
Fig. 10
ROC curves of biomarkers for Spleen-deficiency in prediction. The prediction of spleen deficiency-FD vs. HC was robust, and the area under the curve (AUC) of 15 biomarkers for spleen deficiency-FD vs. HC ranged from 0.6790 to 0.9926. The combination of 15 biomarkers could offer a highly accurate diagnosis of spleen deficiency in FD with an AUC of 0.9943

References

    1. Tack J, Talley NJ. Functional dyspepsia—symptoms, definitions and validity of the Rome III criteria. Nat Rev Gastroenterol Hepatol. 2013;10(03):134–141. doi: 10.1038/nrgastro.2013.14.
    1. Ford AC, Luthra P, Tack J, et al. Efficacy of psychotropic drugs in functional dyspepsia: systematic review and meta-analysis. Gut. 2017;66(03):411–420. doi: 10.1136/gutjnl-2015-310721.
    1. Wang F-Y, Ma X-X, et al. Efficacy and safety of Xiangsha Liujunzi granules for functional dyspepsia: a multi-center randomized double-blind placebo-controlled clinical study. World J Gastroenterol. 2017;23(30):5589–5601. doi: 10.3748/wjg.v23.i30.5589.
    1. Lin L, Wang FY, Tang XD, et al. Effect of Pixu I recipe on cytochrome c oxidase subunit va expression in gastric tissues of FD rats with spleen deficiency. Chin J Exp Tradit Med Formulae. 2017;23(07):102–108.
    1. Zhong ZS, Zhang HY, Zhang W, et al. Effect of Si Junzitang on CaM-MLCK pathway in rats with spleen deficiency and gastrointestinal motility disorder. Chin J Exp Tradit Med Formulae. 2018;24(05):95–99.
    1. Xin X, Lin H, Ren L. Metabonomics research of functional dysplasia of hepatic stagnation and spleen deficiency type. World Chin J Digestol. 2017;25(29):2591–2596. doi: 10.11569/wcjd.v25.i29.2591.
    1. Mahadeva S, Goh KL. Epidemiology of functional dyspepsia: a global perspective[J] World J Gastroenterol. 2006;12(17):2660–2666. doi: 10.3748/wjg.v12.i17.2661.
    1. Tack J, Carbone F. Functional dyspepsia and gastroparesis. Curr Opin Gastroenterol. 2017;33(00):1–9.
    1. Sarnelli G, Caenepeel P, Geypens B, et al. Symptoms associated with impaired gastric emptying of solids and liquids in functional dyspepsia. Am J Gastroenterol. 2003;98(4):783–788. doi: 10.1111/j.1572-0241.2003.07389.x.
    1. Talley NJ, Walker MM, Holtmann G. Functional dyspepsia. Curr Opin Gastroenterol. 2016;32:467–473. doi: 10.1097/MOG.0000000000000306.
    1. Enck P, Azpiroz F, Boeckxstaens G, et al. Functional dyspepsia. Nat Rev Dis Primers. 2017;3:17081–17099. doi: 10.1038/nrdp.2017.81.
    1. Stanghellini V, Cogliandro R. Review article: adherence to Rome criteria in therapeutic trials in functional dyspepsia. Aliment Pharmacol Ther. 2014;40(5):435–466.
    1. Stanghellini V, Tack J. Gastroparesis: separate entity or just a part of dyspepsia? Gut. 2014;63(12):1972–1978. doi: 10.1136/gutjnl-2013-306084.
    1. Kumar A, Pate J, Sawant P. Epidemiology of functional dyspepsia. J Assoc Phys India. 2012;60(3):9–12.
    1. Stanghellini V, Chan FK, Hasler WL, et al. Gastroduodenal disorders. Gastroenterology. 2016;150(6):1380–1392. doi: 10.1053/j.gastro.2016.02.011.
    1. Wang HM, Zhang LP, Chen LR, et al. Gastrointestinal motility and spleen deficiency in gastroesophageal reflux disease. Chin J Integr Tradit West Med Digest. 2015;23(3):196–198.
    1. Chen YL, Fu R, Liu XJ, et al. Literature research of spleen asthenic syndrome symptoms regularity. China J Chin Med. 2013;28(1):150–152.
    1. Wei-Hua XU, Ni-Jiao LI, Zhang YL. Study of radionuclide gastric emptying on different TCM syndrome types in patients with functional dyspepsia. Glob Trad Chin Med. 2013;6(5):321–324.
    1. Zhao-hong SHI, Li-he FU, Lei ZHAO, et al. Clinical study on differences of gastric receptivity and emptying function of functional dyspepsia patients with different syndromes. J Nanjing Univ Tradit Chin Med. 2017;33(6):566–569.
    1. Deng JM, Wang HB, Xiang LI, et al. A clinical trial on functional dyspepsia with spleen deficiency and qi stagnation treated by the decoction of invigorating spleen and regulating qi and its effect on radionuclide gastric emptying. Chin J Integr Tradit West Med Digest. 2012;20(5):193–196.
    1. Wei-hua XU, Wei WANG, Ni-jiao LI, et al. Effects of modified Xiangsha Liujunzi Decoction combined with Zhizhu Pills in the treatment of functional dyspepsia with spleen deficiency and qi stagnation and its influence on radio nuclide gastric emptying. China J Tradit Chin Med Pharm. 2017;32(3):1025–1028.
    1. Khoo J, Rayner CK, Feinlebisset C, et al. Gastrointestinal hormonal dysfunction in gastroparesis and functional dyspepsia. Neurogastroenterol Motil. 2010;22(12):1270–1278. doi: 10.1111/j.1365-2982.2010.01609.x.
    1. Ogiso K, Asakawa A, Amitani H, et al. Ghrelin: a gut hormonal basis of motility regulation and functional dyspepsia. J Gastroenterol Hepatol. 2011;26(s3):67–72. doi: 10.1111/j.1440-1746.2011.06630.x.
    1. Monkemuller K, Neumann HM, Vieth M, et al. Serum gastrin and pepsinogens do not correlate with the different grades of severity of gastro-oesophageal reflux disease: a matched case-control study. Aliment Pharmacol Ther. 2010;28(4):491–496. doi: 10.1111/j.1365-2036.2008.03769.x.
    1. Garashi M, Nagano J, Tsuda A, et al. Correlation between the serum pepsinogen I level and the symptom degree in proton pump inhibitor-users administered with a probiotic. Pharmaceuticals. 2014;7(7):754–764. doi: 10.3390/ph7070754.
    1. Germaná B, Di MF, Cavallaro LG, et al. Clinical usefulness of serum pepsinogens I and II, gastrin-17 and anti-Helicobacter pylori antibodies in the management of dyspeptic patients in primary care. Digest Liver Dis. 2005;37(7):501–508. doi: 10.1016/j.dld.2005.01.016.
    1. Väänänen H, Vauhkonen M, Helske T, et al. Non-endoscopic diagnosis of atrophic gastritis with a blood test. Correlation between gastric histology and serum levels of gastrin-17 and pepsinogen I: a multicentre study. Eur J Gastroenterol Hepatol. 2003;15(8):885–891. doi: 10.1097/00042737-200308000-00009.
    1. Sanger GJ, Lee K. Hormones of the gut-brain axis as targets for the treatment of upper gastrointestinal disorders. Nat Rev Drug Discov. 2008;7(3):241–254. doi: 10.1038/nrd2444.
    1. Sanger GJ. Motilin, ghrelin and related neuropeptides as targets for the treatment of GI diseases. Drug Discov Today. 2008;13(5):234–239. doi: 10.1016/j.drudis.2007.10.024.
    1. Wu Q, Zhang Q, Sun B, et al. 1H NMR-based metabonomic study on the metabolic changes in the plasma of patients with functional dyspepsia and the effect of acupuncture. J Pharm Biomed Anal. 2010;51(3):698–704. doi: 10.1016/j.jpba.2009.09.042.
    1. Pilichiewicz AN, Horowitz M, Holtmann GJ, Talley NJ, Feinle-Bisset C. Relationship between symptoms and dietary patterns in patients with functional dyspepsia. Clin Gastroenterol Hepatol. 2009;7:317–322. doi: 10.1016/j.cgh.2008.09.007.
    1. Feinle C, Meier B, Fried M. Role of cognitive factors in symptom induction in functional dyspepsia (FD) after high and low fat meals. Gastroenterology. 2003;124(4):A253. doi: 10.1016/S0016-5085(03)81267-2.
    1. Hishikawa D, Hashidate T, Shimizu T, et al. Thematic Review Series: lysophospholipids and their receptors: diversity and function of membrane glycerophospholipids generated by the remodeling pathway in mammalian cells. J Lipid Res. 2014;55(5):799–807. doi: 10.1194/jlr.R046094.
    1. Robichaud JC, Jn VDV, Yao Z, et al. Hepatic uptake and metabolism of phosphatidylcholine associated with high density lipoproteins. Biochim Et Biophys Acta Gener Subj. 2009;1790(6):538–551. doi: 10.1016/j.bbagen.2009.02.010.
    1. Wu MK, Cohen DE. Phosphatidylcholine transfer protein regulates size and hepatic uptake of high-density lipoproteins. Am J Physiol Gastrointest Liver Physiol. 2005;289(6):G1067–G1074. doi: 10.1152/ajpgi.00194.2005.
    1. Abstract A. Phosphatidylinositol transfer protein controls the phosphatidylcholine content of yeast Golgi membranes. J Cell Biol. 1994;124(3):273–287. doi: 10.1083/jcb.124.3.273.
    1. Abdel-Latif AA. Calcium-mobilizing receptors, polyphosphoinositides, generation of second messengers and contraction in the mammalian iris smooth muscle: historical perspectives and current status. Life Sci. 1989;45(9):757–786. doi: 10.1016/0024-3205(89)90170-7.
    1. Rosel P, Arranz B, Oros M, et al. Different regional distribution of the 5-HT reuptake complex, the 5-HT2A receptors and their second messenger IP3 in human brain. Neuroscience Research Communications. 2015;24(2):107–115. doi: 10.1002/(SICI)1520-6769(199903/04)24:2<107::AID-NRC6>;2-Y.
    1. Nair R, Manikkath J, Hegde AR, et al. Liposome-encapsulated diacyl glycerol and inositol triphosphate-induced delayed oocyte activation and poor development of parthenotes. J Turk German Gynecol Assoc. 2017;18(3):102–109.
    1. Adibhatla RM, Hatcher JF. Cytidine 5′-diphosphocholine (CDP-choline) in stroke and other CNS disorders. Neurochem Res. 2005;30(1):15–23. doi: 10.1007/s11064-004-9681-8.
    1. Fagone P, Jackowski S. Phosphatidylcholine and the CDP-choline cycle. Biochim Biophys Acta. 2013;1831(3):523–532. doi: 10.1016/j.bbalip.2012.09.009.
    1. Nardone G, Compare D. The psyche and gastric functions. Dig Dis. 2014;32(3):206–212. doi: 10.1159/000357851.
    1. Xiao LY, Wang LY, Zhou KM, et al. Effects of herbal cake-separated moxibustion on related receptors and hormones of HPA axis in rat model of functional dyspepsia with liver stagnation and spleen deficiency syndrome. China J Tradit Chin Med Pharm. 2017;32(1):238–241.

Source: PubMed

3
Sottoscrivi