Introduction {#S0001}
============

Post-stroke depression (PSD) refers to self-accusation and self-guilt, lack of pleasure, and loss of interest in patients after stroke. Symptoms can last for more than 2 weeks.[@CIT0001] PSD significantly reduces quality of life, delays recovery of neurological function, increases cost of treatment, and increases the risk of death by 3.4 times within 10 years.[@CIT0002] However, the incidence of PSD is mostly concentrated from the second month to the sixth month after stroke, so how we know the patient's condition depends on the follow-up after discharge.[@CIT0003],[@CIT0004] Untimely follow-up care, poor objectivity, and lack of awareness of the illness by patients and their families may affect the diagnosis of PSD. An objective, accurate prediction of PSD coupled with early intervention during hospitalization may be a promising approach to understand the pathogenesis of PSD and reduce its burden.

Proton magnetic resonance spectroscopy (^1^H-MRS) is a non-invasive imaging modality that observes metabolic changes in damaged brain regions, and is used in a variety of neurological conditions. In recent years, research on depression, schizophrenia, and other mental diseases, especially the metabolic changes in the brain areas related to PSD, has become a hot topic. MRS reflects the changes of neurobiochemical metabolites in the brain. The commonly used neurotransmitters related to depression include choline (Cho), N-acetylaspartic acid (NAA) and creatine (Cr). These neurobiochemical metabolites are closely related to the structural and functional integrity of nerve cells, membrane phospholipid metabolism, neurotransmitter metabolism, second messenger metabolism and high-energy phosphate metabolism. International[@CIT0005] and domestic[@CIT0006] studies have suggested that ^1^H-MRS may be used to obtain a wide range of information related to the occurrence, progression, and prognosis of PSD. ^1^H-MRS can predict the occurrence and guide the treatment of PSD and has the potential to be a novel imaging biomarker. However, there is a dearth of reported research on the establishment of a PSD prediction model utilizing ^1^H-MRS.

A ﻿nomogram prediction model is simple, intuitive, and practical. With the programming advantage of R language, it can be easily drawn with only a series of codes, no matter its subtype or continuous variable.[@CIT0007] This study aimed to use LASSO regression to screen relevant indicators and build a PSD risk prediction model, to explore neurological metabolic mechanisms related to the onset of PSD, and to evaluate the probability of a patient developing PSD, to allow for timely intervention for PSD patients at an early stage.

Patients and Methods {#S0002}
====================

Selection and Grouping of Subjects {#S0002-S2001}
----------------------------------

Individuals who presented at the Department of Neurology of the First Affiliated Hospital of Jinzhou Medical University and hospitalized within 7 days after a stroke from November 2017 to March 2019 were selected as subjects for the study. After routine clinical information collection and medical treatment, when the patients were in stable condition, proton MR spectroscopy (^1^H-MRS) examinations were performed to measure the ratio of choline to creatine (Cho/Cr) and N-acetylaspartate to creatine (NAA/Cr) in brain regions related to emotion. All stroke patients were evaluated using the Hamilton Depression Scale (HAMD) from the second month to the sixth month once a month after the stroke onset. According to the scores, the stroke patients were divided into a post-stroke depression (PSD) group and a post-stroke non-depression (N-PSD) group.

This study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethics Committee of the First Hospital of Jinzhou Medical University. All subjects signed the medical informed consent documents.

Inclusion Criteria {#S0002-S2002}
------------------

### Diagnosis {#S0002-S2002-S3001}

At onset, patients were diagnosed with first-time ischemic stroke confined by the head computed tomography (CT) in accordance with the diagnostic criteria of stroke of the Fourth National Conference on Cerebrovascular Diseases.[@CIT0008]

### Evaluation of Depression {#S0002-S2002-S3002}

Depression was evaluated from the second month to the sixth month once a month after onset by the HAMD scale in accordance with the diagnostic criteria of the third edition of the Chinese Diagnostic and Classification Standard of Mental Disorders.[@CIT0009]

Exclusion Criteria {#S0002-S2003}
------------------

The exclusion criteria were individuals older than 80 years of age; left-handed patients; history of cerebrovascular disease (cerebral infarction or cerebral hemorrhage), head trauma, or brain organic disease; intellectual disability, severe visual acuity, hearing impairment, aphasia, or inability to complete the dementia scale; a personal or family history of depression and anxiety or other mental disorders; suspected history of alcohol or drug abuse; history of taking psychoactive drugs; use of steroid drugs during treatment.

Collection of Clinical Data {#S0002-S2004}
---------------------------

At the time of the first visit, age, sex, level of education, family income, and basic medical history regarded as the risk of stroke (hypertension, diabetes, hyperlipidemia, and heart disease) for the subjects were recorded. Computed tomography (CT), electrocardiography, echocardiography, carotid artery ultrasonography, and transcranial Doppler were performed.

Evaluation Standards {#S0002-S2005}
--------------------

### PSD Diagnosis {#S0002-S2005-S3001}

Patients were evaluated with the HAMD17, which was scored independently by two evaluators; the average of 4-month values was calculated (non-depression, \<7 points; depression, ≥7 points).

### Family Care Function Assessment {#S0002-S2005-S3002}

The Family APGAR ("Adaptability," "Partnership," "Growth," "Affection," and "Resolve") index was used (good family function, 7--10 points; family dysfunction, \<7 points).

### Personality {#S0002-S2005-S3003}

The Eysenck Personality Questionnaire was used to evaluate patient personality as extroverted or introverted.

### Activities of Daily Living (ADL) Scale {#S0002-S2005-S3004}

The ADL scale was adopted to evaluate capacity for normal independent function. A score of 100 points denotes that the patients are capable of self-care; ≥60 points signify the capacity for basic self-care; 60--40 points signify that the patient requires assistance; 40--20 points signify that the patient requires considerable assistance; and \<20 points denote the patient is completely dependent.

### National Institutes of Health Stroke Scale (NIHSS) {#S0002-S2005-S3005}

The NIHSS was used to evaluate the degree of neurological impairment at stroke onset and the 14th day after onset. The scores ranged from 0 to 42 points. The higher the score, the more severe the neurological impairment.

MR Scanning and Analysis {#S0002-S2006}
------------------------

When the patients were in stable condition, mostly on the 14th day after onset, they were scanned using MRI.

### Morphological Information Collection by Cerebral MRI {#S0002-S2006-S3001}

Morphological brain information was collected with a Siemens Magentom Verio 3.0T scanner. The scanning parameters were axis T1 weighted imaging (T1WI) DRIVE: TR=2000 ms, TE=20 ms, IR=800 ms, thickness=6 mm, gap=1 mm, NSA=1, matrix=400×351, TSE Factor=7; T2 weighted imaging (T2WI) DRIVE: TR=3000 ms, TE=80 ms, NSA=1, matrix=484×377, thickness=6 mm, gap=1 mm, TSE Factor=15.

### Proton Magnetic Resonance Spectroscopy (^1^H-MRS) {#S0002-S2006-S3002}

Multivoxel spectral scanning was performed, and the left frontal lobe, left temporal lobe, left anterior cingulate gyrus, left posterior cingulate gyrus, left precuneus, and the contralateral middle cerebellar peduncles of stroke lesion[@CIT0010] were selected as regions of interest from conventional MRI T1WI or T2WI images.[@CIT0011]--[@CIT0013] The voxel size was 1.0 cm x 1.0 cm x 1.5 cm. The structures of the skull, fat, air cavity, and cerebrospinal fluid should be avoided to prevent affecting the examination results. The scanning parameters were TR=1500 ms and TE=35 ms and the average scanning time for each part was 4min12s. The transmitting and receiving adjustment magnetic field and anhydrous suppression scanning were completed automatically, and the full width at half height of the water peak was less than 3--5 Hz.

Changes in metabolites were observed.^1^H-MRS of normal brain tissues showed different types of metabolites at different echo times (TE). In this study, the metabolites detected at long TE were N-acetylaspartic acid (NAA), choline (Cho), and creatine (Cr) and their chemical shifts were 2.02, 3.22, and 3.02 ppm, respectively. The relative concentrations of NAA and Cho were recorded, and NAA/Cr and Cho/Cr were calculated ([Figures 1](#F0001){ref-type="fig"} and [2](#F0002){ref-type="fig"}).Figure 1MR spectroscopic images of the post-stroke group. The MR spectroscopy image is shown on the left side. The right indicates the location of the left prefrontal lobe on MR images.**Abbreviations:** MR, magnetic resonance; NAA, N-acetylaspartic acid; Cr, creatine; Cho, choline; ppm, parts per million.Figure 2MR spectroscopic images of the post-stroke depression group. The MR spectroscopy image is shown on the left side. The right indicates the location of the left temporal lobe on MR images.**Abbreviations:** MR, magnetic resonance; NAA, N-acetylaspartic acid; Cr, creatine; Cho, choline; ppm, parts per million.

Statistical Calculation {#S0002-S2007}
-----------------------

### Screening of Predictors {#S0002-S2007-S3001}

In this study, whether patients suffered from PSD during follow-up was a binary classification problem. The risk factors related to PSD were screened by LASSO (least absolute shrinkage and selection operator) regression through R Programming Language (version 3.6.1, package "glmnet"). For an excellent model fit, a 10-fold cross validation method was used and we finally selected the parameters with modeling ability.

### Drawing of the Forest Plots {#S0002-S2007-S3002}

Stata (Version: 13.1) was used to analyze the selected factors, and the odds ratio (OR) of each risk factor for PSD was calculated, which was expressed as 95% confidence interval (95% CI). An intuitive and easy to understand forest plot was drawn.

### Establishment and Evaluation of the Nomogram Risk Model {#S0002-S2007-S3003}

Based on the final multivariate analysis results, the R Programming Language (package "rms") was used to establish the risk prediction model for PSD.

In this study, Bootstrap method was used for repeated sampling for further internal verification. In order to better verify the accuracy of the model, we promoted the established risk prediction model to other neurology departments of various hospitals in Liaoning Province and asked the local neurologists to help collecting the relevant medical data of stroke patients to facilitate our external validation. Finally, we performed external validation to the model on 382 stroke patients (133 PSD and 249 N-PSD) who were hospitalized from September 2019 to January 2020. (The patients who participated the external validation set were 148 from Shengjing Hospital of China Medical University, 110 from the First Affiliated Hospital of Dalian Medical University, 27 from People's Hospital of Dalian Lyshunkou District, 59 from the First Affiliated Hospital of Jinzhou Medical University and 38 from the Third Affiliated Hospital of Jinzhou Medical University.)

The performance of nomogram can be evaluated by a calibration plot. The calibration plot compares predicted probability and actual occurrence rate by drawing a scatter plot. Each scatter represents the deviation of the predicted probability from the actual occurrence rate. The scatter points of the ideal ﻿nomogram prediction model should all fall near the 45° diagonal line. The area under the curve (AUC) of the receiver operating characteristic (ROC) can be used to evaluate the prediction accuracy of the model. The AUC value ranges from 0.5 to 1.0 and a larger value indicates a stronger discrimination ability. In this study, the calibration plot and ROC curve were drawn using the R statistical programming language (packages "rms" and "ROCR").

Results {#S0003}
=======

Analyses of Data Based on Clinical and ^1^H-MRS Imaging Metabolic Indexes {#S0003-S2001}
-------------------------------------------------------------------------

Sixty-one subjects were included in this study (PSD: n=30; N-PSD: n=31;). The average age, education level, monthly income, APGAR score, ADL score, initial NIHSS score, NIHSS score 14 days after admission, HAMD score (average of 5 months), Cho/Cr, and NAA/Cr of the two groups are shown in [Table 1](#T0001){ref-type="table"}.Table 1Clinical Data and MRS Characteristics of the Internal and External Validation SetsInternal Validation SetExternal Validation SetPSD (n=30)N-PSD (n=31)PSD (n=133)N-PSD (n=249)Age (years)62.33 ± 4.7660.00 ± 4.5461.81 ± 6.6760.60 ± 3.38Sex (male)40%42.86%48.18%45.24%Education time (years)10.40 ± 2.5611.24 ± 3.009.33 ± 1.7210.76 ± 3.14Monthly per capita income (yuan)3786.67 ± 1356.924353.33 ± 1885.994022.38 ± 1905.624203.36 ± 1220.05Basic diseases of stroke (n%)80%80.95%71.96%78.01%Personality (introverted)50.00%52.38%53.11%48.63%APGAR questionnaire8.73 ± 0.708.07 ± 1.039.15 ± 0.298.56 ± 1.12ADL score82.33 ± 8.0095.00 ± 4.3189.52 ± 6.6897.35 ± 2.05Initial NIHSS score4.47 ± 1.063.13 ± 1.853.50 ± 0.213.25 ± 0.73NIHSS score 14 days after admission2.20 ± 0.861.40 ± 0.852.80 ± 1.161.52 ± 0.87HAMD score (average of 5 months)27.33 ± 6.524.00 ± 1.8924.16 ± 4.363.28 ± 1.08Frontal Cho/Cr1.36 ± 0.151.28 ± 0.131.32 ± 0.331.30 ± 0.18Frontal NAA/Cr1.33 ± 0.091.67 ± 0.081.41 ± 0.161.55 ± 0.47Temporal Cho/Cr1.05 ± 0.130.90 ± 0.061.18 ± 0.140.98 ± 0.39Temporal NAA/Cr1.07 ± 0.111.30 ± 0.051.16 ± 0.141.22 ± 0.05Cerebellar Cho/Cr0.70 ± 0.060.65 ± 0.050.69 ± 0.070.62 ± 0.03Cerebellar NAA/Cr1.08 ± 0.051.10 ± 0.051.10 ± 0.051.10 ± 0.06ACC Cho/Cr1.19 ± 0.070.87 ± 0.050.98 ± 0.050.92 ± 0.11ACC NAA/Cr1.60 ± 0.051.59 ± 0.061.60 ± 0.281.62 ± 0.04PCC Cho/Cr0.8 ± 0.030.79 ± 0.040.79 ± 0.050.77 ± 0.17PCC NAA/Cr1.65 ± 0.021.66 ± 0.061.66 ± 0.041.65 ± 0.10Precuneus Cho/Cr0.79 ± 0.040.75 ± 0.050.78 ± 0.080.75 ± 0.10Precuneus NAA/Cr1.64 ± 0.031.64 ± 0.021.63 ± 0.311.67 ± 0.03[^1]

LASSO for Characteristic Parameter Selection {#S0003-S2002}
--------------------------------------------

In order to determine the risk factors related to PSD, we used LASSO regression to screen variables from high-dimensional clinical statistics and ^1^H-MRS imaging metabolic indexes. As shown in [Figure 3](#F0003){ref-type="fig"}, after the 10-fold cross validation, when λ took the minimum value, five factors (ADL score, initial NIHSS score, frontal-lobe Cho/Cr, temporal-lobe Cho/Cr, and anterior cingulate-gyrus Cho/Cr) were finally screened to be related to PSD. [Figure 4](#F0004){ref-type="fig"} shows the regression coefficient plot of 22 factors by LASSO. [Figure 5](#F0005){ref-type="fig"} shows the forest plot after analyzing the selected factors by Stata. ADL score (OR = 0.83; 0.73\~0.92), initial NIHSS score (OR = 3.87; 2.24\~8.64), frontal lobe Cho/Cr (OR = 1.00; 0.65\~1.54), temporal lobe Cho/Cr (OR = 4.12; 2.22\~7.64), anterior cingulate gyrus Cho/Cr (OR = 2.23; 1.83\~5.97) were independent risk factors for PSD ([Figure 5](#F0005){ref-type="fig"}). The ADL score was on the left side of the invalid line, which indicates that it is negatively correlated with PSD. The initial NIHSS score, temporal lobe Cho/Cr, anterior cingulate gyrus Cho/Cr were positively correlated with PSD.Figure 3Elucidation of LASSO coefficient profiles for selected factors. After the 10-fold cross validation, when λ took the minimum value, five factors were finally screened to be related to PSD.**Abbreviation:** LASSO, the least absolute shrinkage and selection operator.Figure 4Parameter selection was performed through LASSO regression. The figure shows the regression coefficient plot with 22 factors by LASSO, with five factors screened (ADL score, initial NIHSS score, frontal-lobe Cho/Cr, temporal-lobe Cho/Cr, and anterior cingulate-gyrus Cho/Cr).**Abbreviation:** LASSO, the least absolute shrinkage and selection operator.Figure 5The forest plot of the odds ratios (ORs) and 95% confidence intervals (CIs) obtained from Stata's analysis of the five independent risk factors. The ADL score is on the left side of the invalid line, which indicates a negative correlation with PSD. The initial NIHSS score, temporal lobe Cho/Cr, anterior cingulate gyrus Cho/Cr are positively correlated with PSD.

Establishment of Nomogram {#S0003-S2003}
-------------------------

A nomogram was established in combination with the above independent prediction factors ([Figure 6](#F0006){ref-type="fig"}). According to the value of each independent risk factor corresponding to the upper scale (0--100), we could determine the corresponding score of this risk factor. The total score was obtained by adding the scores of each factor. Projecting downward from the total score, the corresponding PSD risk prediction probability value could be obtained.Figure 6Post-stroke depression (PSD) risk prediction nomogram. According to the value of each independent risk factor corresponding to the upper scale (0--100), we could get the corresponding score of this risk factor. The total score was obtained by adding the scores of each factor. Projecting downward from the total score, the corresponding PSD risk prediction probability value could be obtained.**Abbreviations:** ACC, anterior cingulate cortex; ADL, activities of daily living; Cho, choline; Cr, creatine; NIHSS, National Institutes of Health Stroke Scale; PSD, post-stroke depression.

Validation of the Prediction Model {#S0003-S2004}
----------------------------------

The consistency test between the predicted value and the actual value through the calibration plot found that the predicted probability of the ﻿nomogram for PSD correlates well with the actual occurrence rate. The AUC values for internal and external validation were 0.8635 (95% CI 0.823 to 0.910) and 0.8851 (95% CI 0.831 to 0.908), respectively, indicating that the model had a better ability to identify and predict the risk of PSD ([Figures 7](#F0007){ref-type="fig"}--[10](#F0010){ref-type="fig"}).Figure 7Scatter diagram of the calibration plot for internal verification of the nomogram. Bootstrap method was used for repeated sampling and internal verification. X axis: predicted prevalence of PSD; Y axis: the actually diagnosed PSD; diagonal dotted line: the perfect prediction of an ideal model; solid line: the performance of the line diagram. The closer the scatter points to the diagonal dashed line, the better the prediction efficiency of the nomogram.**Abbreviation:** PSD, post-stroke depression.Figure 8Scatter diagram of the calibration plot for external verification of the nomogram. Three hundred and eighty-two stroke patients (133PSD, 249N-PSD) admitted to the Department of Neurology from September 2019 to January 2020 in Liaoning province were analyzed for external verification.Figure 9The ROC curve for internal validation of the prediction model. The AUC was 0.8635 (n=61).**Abbreviations:** AUC, area under the ROC curve; PSD, post-stroke depression; ROC, receiver operating characteristic.Figure 10The ROC curve for external validation of the prediction model. The AUC was 0.8851 (n=382).**Abbreviations:** AUC, area under the ROC curve; PSD, post-stroke depression; ROC, receiver operating characteristic.

Discussion {#S0004}
==========

The incidence of PSD varies from 18% to 79%[@CIT0014] and is attributed to the lack of attention to PSD onset during hospitalization, untimely follow-up after discharge, and different diagnostic criteria of PSD. In this study, ^1^H-MRS technology was used to detect the metabolic indexes of emotion﻿-related brain areas. Utilizing the advantages of nomogram such as simple operation, easy understanding, and high accuracy, we hoped to carry out clinical intervention in the early stage of PSD.

Clinical Data Factors {#S0004-S2001}
---------------------

The results of the clinical data in this study showed that the ADL score and the initial NIHSS score were the strongest predictors of PSD and had high clinical relevance in improving the quality of life of stroke patients. It has been well documented that the degree of neurological impairment and the ability to perform ADL at the early stages of stroke are significantly correlated with the occurrence of PSD.[@CIT0015],[@CIT0016] Specifically, the more severe the degree of neurological impairment and the lower the ADL ability, the greater the risk of PSD. Therefore, for patients with cerebrovascular diseases, we should actively promote the recovery of limb function, help patients carry out systematic functional training early on, improve the ability of ADL, and encourage the rehabilitation of physical and psychological health, which may additionally help to avoid secondary PSD and to improve the prognostic quality of life for patients.[@CIT0017]

MRS Imaging Factors {#S0004-S2002}
-------------------

With regard to the study of ^1^H-MRS imaging, major factors affecting the PSD model were the involvement of the left frontal lobe and left temporal lobe Cho/Cr, which were considered to be significantly associated with the occurrence of PSD in previous studies. Most scholars have agreed that depression is accompanied by regional Cho abnormalities in the brain,[@CIT0018] and research has reported Cho increase in depressed patients, suggesting that Cho complexes may have a depressive effect on the central nervous system. Cho, reflecting the total choline compounds reserve in the brain, is the precursor of acetylcholine, a neurotransmitter closely related to emotion and cognition, previously confirmed to be a pathophysiological foundation of affective disorder. In this study, the ratio of Cho/Cr in the brain of PSD patients was significantly higher than that in the brains of patients without PSD. These results suggested that the left prefrontal lobe and temporal hippocampus may show accelerated phospholipid metabolism and abnormal choline changes at the early stage of stroke, which also confirmed the predictive validity of this model. Another key indicator of ^1^H-MRS, the NAA/Cr of the frontal lobe and hippocampus, has been previously reported to be reduced in the brain of PSD patients[@CIT0005] revealing that NAA decrease may be involved in the occurrence of PSD. NAA is mainly distributed in neurons and their axons, is produced by the neuronal mitochondria, and is considered to be a crucial metabolic marker of neuronal integrity. Nevertheless, in this study, the NAA/Cr index of affective brain regions was not selected as a risk factor by the LASSO regression, probably because various nervous system diseases might lead to NAA/Cr decrease. That is, NAA/Cr is not as specific as Cho/Cr as a risk factor for PSD.

Besides the results of our study that the frontal-lobe and temporal-lobe (hippocampal) Cho/Cr can be used as predictors of PSD, it is worth mentioning that in this study, the Cho/Cr of the anterior cingulate gyrus, another structure closely related to emotion, was also entered in the LASSO regression as a risk factor. Unfortunately, ^1^H-MRS study on limbic system components such as the anterior cingulate gyrus and PSD has been rarely reported. Most scholars have supported that the gray matter volume of the anterior cingulate gyrus is decreased in patients with depression.[@CIT0019] It is also well reported that among the ^1^H-MRS imaging indexes of the anterior cingulate gyrus, the change in Glu/Cr is significantly correlated with a decrease in the HAMD score in patients with depression.[@CIT0020] From the perspective of anatomy, the cingulate gyrus is the main brain area for behavioral adjustment and adaptation to the environment, which is the common boundary between cognition and emotion. James Papez[@CIT0021] claimed that emotional regulation was concentrated in the deep structures of our brain and that emotional experience was determined by the close functional relationship between the limbic system, such as the cingulate gyrus, and brain regions such as the frontal lobe, amygdala, and temporal lobe. Based on the LASSO regression and correlation results, this prediction model selected the anterior cingulate-gyrus Cho/Cr as a risk factor of PSD, indicating that the anterior cingulate gyrus played a crucial role in the occurrence and development of neurological dysfunction, such as PSD, and also reflected the rigorous reliability of the model in screening for related factors.

Verification of PSD Prediction Model {#S0004-S2003}
------------------------------------

The evaluation indexes of clinical prediction model mainly include model calibration and model discrimination.

Calibration focuses on the accuracy of the risk prediction value of the model, that is, the consistency between the probability of the event predicted by the model and the actual occurrence rate of the event. The consistency test between the predicted value and the actual value through the calibration plot found that the predicted probability of the ﻿nomogram for PSD correlates well with the actual occurrence rate.

The discrimination focuses on the model's differential ability, which refers to the ability of the model to distinguish the objects with or without terminal events. As a well-differentiated prediction model, the risk score of the objects with terminal events and the risk score of the non-terminal-event objects will be as scattered as possible, but the model with low differentiation does not have such ability. AUC is one of the most commonly used statistical methods to measure discrimination. The AUC values for internal and external validation were 0.8635 and 0.8851, respectively, indicating that the model had good calibration and identification capabilities.

This study has a few limitations such as the small number of patients included, therefore, a larger sample size is required for further study. This study has not yet been terminated, and we will include more research objects that meet the criterion to strengthen the external verification of this risk prediction model. Furthermore, we will also use DTI, rfMRI and other multimodality MRI indicators to establish PSD risk prediction model. We will compare and improve the prediction models of multimodality MRI from the perspectives of convenience, accuracy and price.

Preventive intervention for PSD is a novel concept. The establishment of a PSD prediction model may be conducive to more effective early treatment of PSD, halting its progression and benefiting stroke patients.

Conclusion {#S0005}
==========

The PSD prediction model based on the ADL score, initial NIHSS score, frontal-lobe Cho/Cr, temporal-lobe Cho/Cr, and anterior cingulate-gyrus Cho/Cr has clinical predictive value.

The study was financially supported by the National Natural Science Foundation of China (Grant No 81371461), Science-Technology Foundation of Education Bureau of Liaoning Province, China (Grant No JYTJCZR201903) and Youth Foundation of Education Bureau of Liaoning Province, China (Grant No JYTQN201707). The funding source had no other involvement. We would also like to thank Editage for English language editing.

Disclosure {#S0006}
==========

The authors report no conflicts of interest in this work.

[^1]: **Abbreviations:** APGAR, "Adaptability," "Partnership," "Growth," "Affection," and "Resolve"; ADL, ﻿activities of daily living; HAMD, Hamilton Depression Scale; NIHSS, National Institutes of Health Stroke scale; Cho, choline; Cr, creatine; NAA, N-acetylaspartic acid; ACC, anterior cingulate cortex; PCC, posterior cingulate cortex.
