Feature selection is crucial to develop a brain computer interface (BCI) system which has high classification accuracy and less computational complexity in especially a large feature space. Feature selection (FS) problem has been solved by many various methods. Among these methods, especially evolutionary computation (EC) techniques have gained a lot of attention in recent years. However, there are very few studies in the literature that consider FS problem as a multi-objective problem to find the optimal trade-off between classification ac-curacy and the number of selected features. Therefore, in this paper, a non-dominated sorting multi-objective symbiotic organism search (NSMOSOS) algorithm is proposed to generate the optimal feature subset in BCI. The efficiency and robustness of the proposed algorithm as a feature selection method is investigated in two datasets based on motor imagery. The highest classification accuracies of NSMOSOS for dataset 1 and dataset 2 are obtained 97.86% with 11 features and 96.57% with average 19 features, respectively. The obtained results demonstrate that the proposed method achieves satisfying results with regard to both the classification accuracy improvement and feature reduction rates for both datasets. The superiority of the proposed method is verified compared with the existing methods for both datasets. Besides, three different versions of symbiotic search or-ganism (SOS) algorithm are improved, and pros and cons of these algorithms are evaluated compared with each other. In conclusion, the paper indicates that the proposed NSMOSOS algorithm is an efficient and practicable technique for FS problem and could be helpful in developing the BCI applications.