This paper studies an adaptive neural control for nonlinear multiple-input multiple-output systems in interconnected form. The studied systems are composed of N subsystems in pure feedback structure and the interconnection terms are contained in every equation of each subsystem. Moreover, the studied systems consider the effects of Prandtl-Ishlinskii (PI) hysteresis model. It is for the first time to study the control problem for such a class of systems. In addition, the proposed scheme removes an important assumption imposed on the previous works that the bounds of the parameters in PI hysteresis are known. The radial basis functions neural networks are employed to approximate unknown functions. The adaptation laws and the controllers are designed by employing the backstepping technique. The closed-loop system can be proven to be stable by using Lyapunov theorem. A simulation example is studied to validate the effectiveness of the scheme. © 2015 IEEE.