Abstract-A novel algorithm is proposed for automatic modulation classification in multiple-input multiple-output spatial multiplexing systems, which employs fourth-order cumulants of the estimated transmit signal streams as discriminating features and a likelihood ratio test (LRT) for decision making. The asymptotic likelihood function of the estimated feature vector is analytically derived and used with the LRT. Hence, the algorithm can be considered as asymptotically optimal for the employed feature vector when the channel matrix and noise variance are known. Both the case with perfect channel knowledge and the practically more relevant case with blind channel estimation are considered. The results show that the proposed algorithm provides a good classification performance while exhibiting a significantly lower computational complexity when compared with conventional algorithms.