In this paper, we consider the following Schrödinger systems involving pseudo-differential operator in R n (−∆) α 2 u(x) = u β 1 (x)v τ 1 (x), in R n , (−∆) γ 2 v(x) = u β 2 (x)v τ 2 (x), in R n , (1) where α and γ are any number between 0 and 2, α does not identically equal to γ. We employ a direct method of moving planes to partial differential equations (PDEs) (1). Instead of using the Caffarelli-Silvestre's extension method and the method of moving planes in integral forms, we directly apply the method of moving planes to the nonlocal fractional order pseudo-differential system. We obtained radial symmetry in the critical case and non-existence in the subcritical case for positive solutions. In the proof, combining a new approach and the integral definition of the fractional Laplacian, we derive the key tools, which are needed in the method of moving planes, such as, narrow region principle, decay at infinity. The new idea may hopefully be applied to many other problems.