1. Participant recruitment
1) Recruit adolescents aged 12–17 meeting DSM-5 criteria for MDD.
2) Confirm diagnosis using K-SADS-PL by trained psychologists.
3) Apply exclusion criteria (psychosis, developmental disorders, prior antidepressant use, etc.).
4) Obtain written informed consent.
2. SSRI treatment protocol (8 weeks)
1) Administer escitalopram 5 mg/day during week 1.
2) Increase to 10 mg/day from week 2.
3) Adjust dose at weeks 4 and 6 based on clinical judgment (max 30 mg/day).
4) Assess depressive symptoms using CDRS-R at weeks 0, 2, 4, 6, and 8.
5) Monitor adherence via pill counts; exclude participants <60% adherence.
3. MRI acquisition (approx. 10 min per scan)
1) Scanner: Siemens Tim Trio 3T
2) Sequence: T2*-weighted echo-planar imaging
3) Parameters: TR = 3000 ms; TE = 40 ms; flip angle = 90°; 35 axial slices; voxel size = 3.4 × 3.4 × 4.0 mm; FOV = 240 mm
4) Instruct participants to keep eyes closed, remain awake, and avoid structured thought.
4. Preprocessing (AFNI + FreeSurfer, 4–6 hours per subject)
1) Discard first 4 volumes for signal stabilization.
2) Perform spike correction, slice-timing correction, and motion correction.
3) Align functional data to base volume with minimal outliers.
4) Co-register to structural T1 image.
5) Nonlinear warp to MNI152_T1_2009c template (AFNI).
6) Apply 4 mm FWHM Gaussian smoothing.
7) Regress out white matter and CSF signals using FreeSurfer-derived masks (eroded).
8) Include 6 motion parameters as nuisance regressors.
9) Censor volumes with >0.2 mm motion or >10% outlier voxels.
10) Apply bandpass filter (0.01–0.1 Hz).
5. ROI definition (10 min)
1) Define 10-mm spherical masks for the following DMN seeds:
2) Extract mean time series for each seed.
6. Functional connectivity analysis (3dRegAna, 1–2 hours per subject)
1) Compute voxel-wise correlation maps with each seed’s mean time series.
2) Convert r-values to z-scores via Fisher’s transformation.
3) Run regression analyses predicting change in CDRS-R (week 8 − baseline).
4) Include age and sex as covariates.
5) Correct for multiple comparisons using 3dClustSim (10,000 Monte Carlo iterations).
6) Apply voxel-wise p < 0.005, cluster-level FWE α < 0.017 (Bonferroni correction).
7. Statistical analysis and visualization (SPSS v29, GraphPad Prism)
1) Extract z-scores from significant clusters.
2) Replicate regression analysis in SPSS to confirm direction and significance.
3) Plot scatterplots of baseline rsFC vs. ΔCDRS-R.