ObjectiveTo identify candidate biomarkers associated with neuromyelitis optica spectrum disorder (NMOSD) using high-throughput technologies that broadly assay the concentrations of serum analytes and frequencies of immune cell subsets.MethodsSera, peripheral blood mononuclear cells (PBMCs), and matched clinical data from participants with NMOSD and healthy controls (HCs) were obtained from the Collaborative International Research in Clinical and Longitudinal Experience Study NMOSD biorepository. Flow cytometry panels were used to measure the frequencies of 39 T-cell, B-cell, regulatory T-cell, monocyte, natural killer (NK) cell, and dendritic cell subsets in unstimulated PBMCs. In parallel, multiplex proteomics assays were used to measure 46 serum cytokines and chemokines in 2 independent NMOSD and HC cohorts. Multivariable regression models were used to assess molecular and cellular profiles in NMOSD compared with HC.ResultsNMOSD samples had a lower frequency of CD16+CD56+ NK cells. Both serum cohorts and multivariable logistic regression revealed increased levels of B-cell activating factor associated with NMOSD. Interleukin 6, CCL22, and CCL3 were also elevated in 1 NMOSD cohort of the 2 analyzed. Multivariable linear regression of serum analyte levels revealed a correlation between CX3CL1 (fractalkine) levels and the number of days since most recent disease relapse.ConclusionsIntegrative analyses of cytokines, chemokines, and immune cells in participants with NMOSD and HCs provide congruence with previously identified biomarkers of NMOSD and highlight CD16+CD56+ NK cells and CX3CL1 as potential novel biomarker candidates.