Herba Rhodiolae alleviates depression via the BDNF/TrkB-GSK-3β signaling pathway
Introduction
Depression is one of the most common psychiatric disorders worldwide. The illness is clinically characterized by various symptoms including persistent sadness, loss of interest and enjoyment, feelings of guilt, anxiety, fatigue, low self-worth, poor concentration, and recurrent suicidal thoughts (1). Depression threatens the physical and psychological health of patients and greatly affects their work and personal lives, contributing to the increasing rate of suicide of this century (2,3). Previous findings have demonstrated a complex of factors contributing to the development of depression, including genetic, environmental stress, endocrine, and immunologic factors. Some molecular pathways are associated with the development of depression, including reductions of brain monoaminergic transmission, function of neurotransmitter receptors and neurotrophic factors, increase of proinflammatory cytokines, abnormality of nitric oxide, oxidative disturbances, and dysregulation of the hypothalamic-pituitary-adrenal axis (4). Nonetheless, the pathogenesis of depression has not been well clarified and effective antidepressant agents require further exploration. Physiological mechanisms of an anti-depressive effect of physical exercise in major depressive disorder has been believed to be closely linked with brain-derived neurotrophic factor (BDNF) (5-7), and BDNF/TrkB/GSK-3β pathway has been proved to be closely related with the development of depression (5-7). GSK-3β is an enzyme that phosphorylates glycogen synthase, which in turn inhibits glycogen biosynthesis. Moreover, TrkB/GSK-3β is now believed to play an important role in the pathophysiology of depression and is implicated to be a drug target for the treatment of depression (7).
There are 3 types of treatment targeting depression, namely, electroconvulsive therapy, psychotherapy, and pharmacotherapy. Among them, electroconvulsive therapy is the oldest and most controversial treatment, given the cognitive side effects, and its application has dwindled despite it having been proven effective in some cases (8). The psychotherapeutic techniques are composed of cognitive behavioral therapy, interpersonal psychotherapy, problem-solving therapy, non-directive supportive therapy, behavioral activation therapy, and exercise. Psychological treatment is effective for moderate depression, but its efficacy appears to be limited for chronic and severe depression (9). Antidepressant medications have become increasingly prevalent in recent decades, and the drugs mainly include tricyclic antidepressants, monoamine oxidase inhibitors, selective serotonin reuptake inhibitors, and serotonin norepinephrine reuptake inhibitors (10). However, such drugs may only be effective for specific patients and long-term utilization can lead to treatment resistance (11). Besides, side effects and pathogenesis have been reported as associated with most pharmacotherapies (12).
Traditional Chinese medicine (TCM) has been increasingly applied to treat depression in China (13). Alongside the associated safety advantages, TCM herbal formulae could be more effective and less prone to treatment resistance, given its multiple targets, compared to conventional chemical drugs (14). Rhodiola rosea root, an herbal medicine known as “golden root” or “arctic root”, has been used as an antidepressant for years, and clinical trials have shown that it is effective (15). Previous findings have indicated that Rhodiola rosea contains hundreds of bioactive compounds, for example, phenylpropanoids, phenylethanol derivatives, flavanoids, monoterpernes, triterpenes, and phenolic acids (16), which are supposed to be involved in neuroprotective, cardio-protective, anti-fatigue, anti-stress, anti-depressive, and anxiolytic activities, and immunity improvement (17). Nonetheless, the molecular mechanism of Rhodiola rosea in the treatment of depression is still unknown.
We present the following article in accordance with the ARRIVE reporting checklist (available at https://dx.doi.org/10.21037/atm-21-5849).
Methods
Animals
Adult male Sprague-Dawley (SD) rats (8 weeks, 220–240 g) were purchased from Hangzhou Medical College, Hangzhou, China. The rats were separately housed in the animal facility under the 12 h light-dark cycle at 23–25 °C. The humidity was around 50%, and the rats were provided free access to food and water. After 7 days of acclimatization, 30 rats were randomly divided into 5 groups. The rats in group I (Sham group) were treated orally with physiological saline (1 mL per day). Rats in all other groups were exposed to chronic unpredictable mild stress (CUMS), combined with physiological saline (Oral, 1 mL per day) in group II, Herba Rhodiolae suspension (Oral, 0.04 g/kg, once a day) in group III, Herba Rhodiolae (Oral, 0.04 g/kg, once a day) and glycogen synthase kinase 3 beta (GSK-3β) inhibitor (Intravenous injection, 10 µM, 1 mL, once a day) group IV, and venlafaxine suspension (Oral, 0.035 g/kg, once a day) in group V. Venlafaxine, the common used anti-depression drug in the clinic, was chose in this study. Herba Rhodiolae was purchased from Beijing Tongrentang (Beijing, China), and the dose of Herba Rhodiolae was determined through our pre-experiment and dose conversion from human to rat. Animal experiments were performed under a project license (#2021-016) granted by the Ethics Committee of The Second Affiliated Hospital of Fujian Traditional Chinese Medical University, in compliance with the institutional guidelines for the care and use of animals.
CUMS protocol
The rats in group II-V were exposed to CUMS to induce depression, according to a modified protocol described by Dubey et al. (2015) and Lu et al. (2017) (12,18). The stressors used in the present study consisted of water deprivation for 24 h, food deprivation for 24 h, soiled bedding for 24 h, cage tilting (45°) for 24 h, inversion of light-dark cycle for 24 h, social isolation for 12 h, restraint stress for 2 h, hanging by the tail for 30 min, and swimming in cold water (4 °C) for 5 min. Each rat was randomly subjected to 1 stressor every day, with each stressor applied to each rat at least 1 time but less than 3 times. The CUMS stimulation was performed for 14 consecutive days.
Behavioral tests
The fluid consumption test was performed to evaluate potential anhedonia. Prior to the test, the rats were supplied with 2 bottles of 1% sucrose solution on day 1, and 1 bottle of 1% sucrose solution and 1 bottle of pure water on day 2. On day 3, the rats were also provided 1 bottle of sucrose solution and 1 bottle of pure water, but the position of the 2 bottles were switched every 6 h. After 3 days of acclimatization, the rats were subjected to food and water deprivation for 12 h, and then was given 1 bottle of sucrose solution and 1 bottle of pure water. After 12 h, the sucrose preference was calculated as follows: sucrose preference (%) = sucrose solution consumption/(sucrose solution consumption + water consumption) ×100%.
To evaluate locomotor activity, the open field test was conducted in a dark box with a square floor (80 cm × 80 cm) and walls 45 cm in height. The floor was further subdivided into 16 equal squares (20 cm × 20 cm) with white stripes. The rat was initially placed in the center of floor, and its behavior was recorded for 5 min with a camera, after environmental adaptation for 2 min. The score of locomotor activity was estimated by the number of crossing behaviors, with 3 paws crossing the line accounting for 1 score. Each rearing accounted for 1 score. The box was cleaned using 70% ethanol in water after each test.
To evaluate susceptibility to negative mood, the forced swimming test was conducted in a cylindrical glass (diameter 14 cm, height 20 cm) filled with water at 23 °C with a depth of 10 cm. After swimming training (10 min) at 24 h before the test, each rat was placed individually in the glass and video recorded for 5 min. The immobility duration was calculated for each rat.
Tissue sampling and brain morphology
After behavioral tests, the rats were sacrificed with a lethal dose of pentobarbital. Hippocampi were subsequently separated from brains of 3 individuals in each group, and hippocampus samples were immediately frozen in liquid nitrogen and stored at −80 °C. For the other 3 animals in each group, brain tissues including the hippocampus regions were fixed overnight with the mixture of 2.5% glutaraldehyde and 2.5% paraformaldehyde at 4 °C. The fixed samples were washed with 0.1 M cold phosphate-buffered saline (PBS) and post-fixed with 1% osmium in 0.1 M PBS for 1 h at 4 °C. The samples were dehydrated with a gradient of acetone solution (30–100%) and the dehydrated samples were embedded in epoxy resin overnight. The prepared samples were cut into semi-thin sections (100–120 nm), stained with 1% toluidine blue for 5 min, and blocked with rubber. The semi-thin sections were further cut into ultra-thin sections (40–60 nm) and stained with uranyl acetate for 30 min. The ultra-thin sections were observed by JEM-1200EX transmission electron microscopy (TEM, JEOL, Tokyo, Japan).
High performance liquid chromatography
The glutamate level in the hippocampus was determined by high performance liquid chromatography (HPLC; Agilent 1260, Thermo Fisher Scientific, Waltham, MA, USA), according to the protocol described by Cui et al. (2017) (19). The hippocampus sample was added to homogenization buffer and centrifuged to obtain supernatant, which was suspended with O-phthaldialdehyde. The derivatization mixture was eluted in an Easy C18 analytical column (Agilent, 150 mm × 75 µm, 3 µm particle size), using solvent A (50 mmol/L sodium acetate buffer) and solvent B (methanol). The gradient elution was run at a flow rate of 0.1 mL/min, during which solvent B concentration was initially at 30% for 5 min, increased to 37% for 15 min, and returned to 30% for 20 min. The column temperature was maintained at 35 °C. The glutamate was detected using a fluorescence detector, of which the excitation and emission wavelengths were set at 340 and 440 nm, respectively.
Western blot
The proteins were isolated from hippocampus samples using radio-immunoprecipitation assay (RIPA) lysis buffer, including a cocktail of protease, phosphatase inhibitors, and phenylmethylsulfonyl fluoride (PMSF; Beyotime, Shanghai, China). The isolated proteins were separated on a 10% sodium dodecyl-sulfate polyacrylamide gel electrophoresis (SDS-PAGE) gel and transferred to polyvinylidene difluoride (PVDF) membranes (Bio-Rad, Hercules, CA, USA). After blocking with 5% non-fat milk in Tris-buffered saline containing 0.1% Tween 20 (TBST), the membrane was incubated overnight at 4 °C with monoclonal rabbit primary antibodies specific to tropomyosin receptor kinase B (TrkB), brain derived neurotrophic factor (BDNF), P-GSK3-β, GSK3-β, LC3I, LC3II, Beclin 1, and glyceraldehyde 3-phosphate dehydrogenase (GAPDH; acting as internal control). After washing with TBST, the membrane was incubated for 1 h at room temperature with horseradish peroxidase (HRP)-conjugated secondary antibodies. The protein band was visualized using the Immun-StarTM HRP Chemiluminescence Kit (Bio-Rad) with ImageQuant LAS 4000 system (GE Healthcare, Hino, Japan). The intensities of protein bands were estimated using Image J analysis software (https://imagej.nih.gov/ij/).
Statistical analysis
One-way analysis of variance (ANOVA) test was applied for comparison among groups, in combination with Turkey’s post-hoc test for multiple comparisons. Prior to one-way ANOVA, the normality of data was verified by the Shapiro-Wilk test. All statistical analyses were performed on Graph Pad Prism 5.0 (Graph Pad Software, San Diego, CA, USA). A P value <0.05 indicated a significant difference.
Results
Behavioral performance
In the sucrose preference test, the percentage of sucrose consumption in group I (51.59%±3.85%) was significantly greater than that in the other 4 groups [F(4,28)=80.92, P<0.05]. The sucrose consumption of group III (29.68%±3.85%) was significantly greater than that in group IV (24.87%±3.76%) and group II (29.68%±3.85%). The percentage of sucrose consumption in group II was significantly smaller than that in group V (29.68%±2.80%) (Figure 1A). In the open field test, significant differences were observed in all comparisons between groups [F(4,28)=144.50, P<0.05]. The average movement score was 106.80 (SD 3.95) in group I, 91.86 (SD 3.54) in group IV, 81.63 (SD 4.67) in group V, 68.09 (SD 1.92) in group III, and 55.05 (SD 5.99) in group II (Figure 1B). In the forced swimming test, significant differences were observed in all comparisons between different groups [F(4,28)=201.10, P<0.05]. The average of immobility time was 239.30 seconds (SD 10.33) in group II, 201.30 (SD 11.07) in group III, 178.60 (SD 11.75) in group V, 151.10 (SD 10.98) in group IV, and 79.50 (SD 6.89) in group I (Figure 1C).
Protein expression
The expression level of TrkB was significantly decreased in group II-V (0.53±0.03 for group II, 0.66±0.06 for group III, 0.76±0.04 for group IV, and 0.64±0.01 for group V), compared with group I (1.00±0.01). The level of TrkB in group II was significantly lower than that in group III-V (F(4,10)=68.70, P<0.05). Significant differences were observed between group IV and group V (Figure 2). The expression of BDNF was significantly decreased in group II-V (0.47±0.04 for group II, 0.68±0.10 for group III, 0.65±0.09 for group IV, and 0.59±0.04 for group V), compared with group I (1.00±0.03). The level of BDNF was significantly lower compared with group III and IV [F(4,10)=27.16, P<0.05, Figure 2]. The ratio of LC3II and LC3I was significantly increased in group II-V (1.96±0.07 for group II, 1.53±0.15 for group III, 1.54±0.04 for group IV, and 1.23±0.12 for group V), compared with group I (1.00±0.05). The expression ratio of LC3II and LC3I in group II was significantly higher than that of group III-V [F(4,10)=43.13, P<0.05, Figure 2]. The expression level of P-GSK/GSK was significantly increased in group II-V (4.20±0.07 for group II, 2.58±0.27 for group III, 1.70±0.06 for group IV, and 1.89±0.16 for group V), compared with group I (1.00±0.08). Significant differences were observed in Beclin 1 expression between groups (1.00±0.06 for group I, 1.93±0.07 for group II, 1.73±0.01 for group III, and 1.40±0.05 for group IV), except for group I and group V (1.11±0.13) [F(4,10)=85.13, P<0.05, Figure 2].
TEM and glutamate concentration
The glutamate concentration in group I (589.7± 41.79) was significantly lower than that in group II-V (848.00±29.31 for group II, 750.30±5.69 for group III, 678.70±22.50 for group IV, and 713.30±8.51 for group V) (Figure 3). The glutamate concentration in group II was significantly greater than that in group III-V. The glutamate concentration in group III was significantly higher than that in group IV. No differences were observed in comparisons between group III and group V, group IV and group V [F(4,10)=41.94, P<0.05] (Figure 3). No auto-phagosome was observed in brain neurons of rats in the group I and group IV. Significant increases of auto-phagosomes number in brain neurons, especially in the hippocampus region of animals were observed in group II compared with group III and V. A larger number of auto-phagosomes were found in group V than in group III (Figure 4).
Discussion
Compared with the sham group, the sucrose consumption and locomotor activity of CUMS rats were significantly decreased, and the immobility time was contrarily increased, in agreement with previous findings (20). The decrease of sucrose consumption and locomotor activity, anhedonia, loss of interesting, and despair, respectively, are the core symptoms of human depression (21,22). The depressive-like behaviors were probably related to the impairments of the central nervous system in CUMS rats, given the abnormalities in brain morphology. In morphology, a large number of auto-phagosomes were induced by CUMS stimulation in neurons, especially in the hippocampus regions. Auto-phagosome is the unique organelle involving in autophagy, which is the intracellular degradation system transporting cytoplasmic constituents to the lysosome. The autophagy activated in the brain was in consistent with the increased levels of the autophagy-related proteins Beclin-1 and LC3II/I.
Abnormality of the central nervous system was observed in physiological examinations. Specifically, glutamate level in the brain was elevated in CUMS rats compared to the sham group, which was in agreement with findings in humans (23). Glutamate is a major excitatory neurotransmitter in the central nervous system, and it is involved in synaptic plasticity and synaptogenesis (24). Glutamate released from presynaptic neurons could bind to glutamate receptors of postsynaptic neurons. Glutamate receptors mainly include α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid (AMPA) receptor, N-methyl-D-aspartate (NMDA) receptor, and metabotropic glutamate receptors (mGluRs). The excessive glutamate likely resulted in neurotoxic-like effects in the central nervous system (25). Intracellular BDNF/TrkB signaling, which closely interacts with the glutamate system in neurons, has been demonstrated to be involved in pathological mechanism of depression. The BDNF/TrkB signaling has been shown to activate the PI3K/mTOR and Ras/ERK pathways and play an important role in the maintenance and survival of neurons and synaptic plasticity (26,27). In agreement with previous findings (28), the levels of BDNF and TrkB were significantly reduced by CUMS stimulation in the present study, which were associated with the development of depression. In concert with the BDNF/TrkB signaling (28,29), the expression of GSK-3β, a component of PI3K/mTOR and Ras/ERK pathways, was contrarily increased in CUMS rats.
The depressive-like behaviors induced by CUMS stimulation were alleviated by the treatments of Herba Rhodiolae, combination of Herba Rhodiolae and GSK-3β inhibitor, and venlafaxine. Among these treatments, the combination of Herba Rhodiolae and GSK-3β inhibitor appeared to be most effective, incorporating the benefits of both Herba Rhodiolae and GSK-3β inhibitor. The discrepancy in treatment effect likely resulted from different effects on pathological processes. Specifically, Herba Rhodiolae, GSK-3β inhibitor, and venlafaxine were understood to affect the auto-phagosome number and glutamate level stimulated by CUMS in the central nervous system. Herba Rhodiolae was more effective in rehabilitating BDNF/TrkB signaling than venlafaxine. The rehabilitation of BDNF/TrkB signaling was not observed after treatment with GSK-3β inhibitor. Nonetheless, GSK-3β inhibitor was the most effective in down-regulation of GSK-3β expression compared to Herba Rhodiolae and venlafaxine. These findings indicated that the combined treatments of Herba Rhodiolae and GSK-3β inhibitor rehabilitate the BDNF/TrkB-GSK-3β signaling pathway.
Acknowledgments
Funding: This study was supported by the Fujian Provincial Health Commission (2019-ZQN-78) and Natural Science Foundation of the Fujian Province (2020J01247)
Footnote
Reporting Checklist: The authors have completed the ARRIVE reporting checklist. Available at https://dx.doi.org/10.21037/atm-21-5849
Data Sharing Statement: Available at https://dx.doi.org/10.21037/atm-21-5849
Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://dx.doi.org/10.21037/atm-21-5849). The authors have no conflicts of interest to declare.
Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. Animal experiments were performed under a project license (#2021-016) granted by the Ethics Committee of The Second Affiliated Hospital of Fujian Traditional Chinese Medical University, in compliance with the institutional guidelines for the care and use of animals.
Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.
References
- Nestler EJ, Barrot M, DiLeone RJ, et al. Neurobiology of depression. Neuron 2002;34:13-25. [Crossref] [PubMed]
- Kessler RC, Berglund P, Demler O, et al. Lifetime prevalence and age-of-onset distributions of DSM-IV disorders in the National Comorbidity Survey Replication. Arch Gen Psychiatry 2005;62:593-602. [Crossref] [PubMed]
- Hawton K, Casañas I, Comabella C, Haw C, et al. Risk factors for suicide in individuals with depression: a systematic review. J Affect Disord 2013;147:17-28. [Crossref] [PubMed]
- Chopra K, Kumar B, Kuhad A. Pathobiological targets of depression. Expert Opin Ther Targets 2011;15:379-400. [Crossref] [PubMed]
- Kallies G, Rapp MA, Fydrich T, et al. Serum brain-derived neurotrophic factor (BDNF) at rest and after acute aerobic exercise in major depressive disorder. Psychoneuroendocrinology 2019;102:212-5. [Crossref] [PubMed]
- Amidfar M, Kim YK, Wiborg O. Effectiveness of memantine on depression-like behavior, memory deficits and brain mRNA levels of BDNF and TrkB in rats subjected to repeated unpredictable stress. Pharmacol Rep 2018;70:600-6. [Crossref] [PubMed]
- Chiu CH, Chyau CC, Chen CC, et al. Erinacine A-Enriched Hericium erinaceus Mycelium Produces Antidepressant-Like Effects through Modulating BDNF/PI3K/Akt/GSK-3beta Signaling in Mice. Int J Mol Sci 2018;19:341. [Crossref] [PubMed]
- Barbui C, Butler R, Cipriani A, et al. Depression in adults: drug and physical treatments. BMJ Clin Evid 2007.
- Cuijpers P, Karyotaki E. The effects of psychological treatment of perinatal depression: an overview. Arch Womens Ment Health 2021;24:801-6. [Crossref] [PubMed]
- Holtzheimer PE 3rd, Nemeroff CB. Advances in the treatment of depression. NeuroRx 2006;3:42-56. [Crossref] [PubMed]
- Little A. Treatment-resistant depression. Am Fam Physician 2009;80:167-72. [PubMed]
- Dubey VK, Ansari F, Vohora D, et al. Possible involvement of corticosterone and serotonin in antidepressant and antianxiety effects of chromium picolinate in chronic unpredictable mild stress induced depression and anxiety in rats. J Trace Elem Med Biol 2015;29:222-6. [Crossref] [PubMed]
- Pan YJ, Cheng IC, Yeh LL, et al. Utilization of traditional Chinese medicine in patients treated for depression: a population-based study in Taiwan. Complement Ther Med 2013;21:215-23. [Crossref] [PubMed]
- Li C, Huang J, Cheng YC, et al. Traditional Chinese Medicine in Depression Treatment: From Molecules to Systems. Front Pharmacol 2020;11:586. [Crossref] [PubMed]
- Gao L, Wu C, Liao Y, et al. Antidepressants effects of Rhodiola capsule combined with sertraline for major depressive disorder: A randomized double-blind placebo-controlled clinical trial. J Affect Disord 2020;265:99-103. [Crossref] [PubMed]
- Ming DS, Hillhouse BJ, Guns ES, et al. Bioactive compounds from Rhodiola rosea (Crassulaceae). Phytother Res 2005;19:740-3. [Crossref] [PubMed]
- Kamel Ismail ZM, Morcos MA, Eldin Mohammad MD, et al. Enhancement of Neural Stem Cells after Induction of Depression in Male Albino Rats (A histological & Immunohistochemical Study). Int J Stem Cells 2014;7:70-8. [Crossref] [PubMed]
- Lu J, Shao RH, Jin SY, et al. Acupuncture ameliorates inflammatory re-sponse in a chronic unpredictable stress rat model of depression. Brain Res Bull 2017;128:106-12. [Crossref] [PubMed]
- Cui T, Qiu HM, Huang D, et al. Abnormal levels of seven amino neuro-transmitters in depressed rat brain and determination by HPLC-FLD. Biomed Chromatogr 2017;31:e3937. [Crossref]
- Wang Q, Timberlake MA 2nd, Prall K, et al. The recent progress in animal models of depression. Prog Neuropsychopharmacol Biol Psychiatry 2017;77:99-109. [Crossref] [PubMed]
- Bai M, Zhu X, Zhang Y, et al. Abnormal hippocampal BDNF and miR-16 expression is associated with depression-like behaviors induced by stress during early life. PLoS One 2012;7:e46921. [Crossref] [PubMed]
- Anisman H, Matheson K. Stress, depression, and anhedonia: caveats concerning animal models. Neurosci Biobehav Rev 2005;29:525-46. [Crossref] [PubMed]
- Hashimoto K. Role of the mTOR signaling pathway in the rapid antidepressant action of ketamine. Expert Rev Neurother 2011;11:33-6. [Crossref] [PubMed]
- Hashimoto K. Emerging role of glutamate in the pathophysiology of major depressive disorder. Brain Res Rev 2009;61:105-23. [Crossref] [PubMed]
- Murrough JW, Abdallah CG, Mathew SJ. Targeting glutamate signalling in depression: progress and prospects. Nat Rev Drug Discov 2017;16:472-86. [Crossref] [PubMed]
- Dwivedi Y. Brain-derived neurotrophic factor: role in depression and suicide. Neuropsychiatr Dis Treat 2009;5:433-49. [Crossref] [PubMed]
- Leal G, Afonso PM, Salazar IL, et al. Regulation of hippocampal synaptic plasticity by BDNF. Brain Res 2015;1621:82-101. [Crossref] [PubMed]
- Zhang SS, Tian YH, Jin SJ, et al. Isoflurane produces antidepressant effects inducing BDNF-TrkB signaling in CUMS mice. Psychopharmacology (Berl) 2019;236:3301-15. [Crossref] [PubMed]
- Duda P, Hajka D, Wojcicka O, et al. GSK3beta: A Master Player in Depressive Disorder Pathogenesis and Treatment Responsiveness. Cells 2020;9:727. [Crossref] [PubMed]