Statistical analyses
Statistical analyses were conducted in R (version 3.5.2). Variables with skewed distribution are presented as median and interquartile range (75 percentile—25 percentile). Response to treatment was defined as ≥ 50% reduction in IDS-C30 scores at week 12, relative to baseline. Differences between responders and non-responders were assessed by exact Wilcoxon–Mann–Whitney test. Partial Least-Squares Discriminant Analysis (PLS-DA) of the change in plasma lipid mediators was performed on normalized and scaled data using MetaboAnalyst 5.0 (