Multiple sclerosis (MS) is an inflammatory demyelinating disease with heterogeneous clinical presentations and variable long-term disability accumulation. There are currently no standard criteria to accurately predict disease outcomes. In this study we investigated the cross-sectional relationship between disease phenotype and immune-modulating cytokines and chemokines in cerebrospinal fluid (CSF). We analyzed CSF from 20 DMT-naïve MS patients using Olink Proteomics’ Target 96 Inflammation panel and correlated the resulting analytes with respect to (1) disease subtype, (2) patient age and sex, (3) extent of clinical disability, and (4) MRI segmental brain volumes. We found that intrathecal IL-4 correlated with higher Expanded Disability Status Scale (EDSS) scores and longer 25-foot walk times, and CD8A correlated with decreased thalamic volumes and longer 9-hole peg test times. Male sex was associated with higher FGF-19 expression, and Tumefactive MS with elevated CCL4. Several inflammatory markers were correlated with older age at the time of LP. Finally, higher intrathecal IL-33 correlated with increased MS lesion burden and multi-compartment brain atrophy. This study confirms immune heterogeneity underlying CSF profiles in MS, but also identifies several inflammatory protein biomarkers that may be of use for predicting clinical outcomes in future algorithms.