This paper focuses on the stability and convergence analysis of the first-order Euler implicit/explicit scheme based on mixed finite element approximation for three-dimensional (3D) time-dependent MHD equations. Firstly, for initial data u0,B0∈Hαu0,B0∈Hα with α=1,2α=1,2, the regularity results of the continuous solution (u,p,B)(u,p,B) and the spatial semi-discretization solution (uh,ph,Bh)(uh,ph,Bh) are obtained, and L2L2-error estimate of (uh,Bh)(uh,Bh) is deduced by using the negative norm technique. Next, through the use of mathematic induction, the H2H2-stability of the fully discrete first-order scheme is proved under the stability condition depending on the smoothness of initial data. Here, the stability condition is Δt≤C0Δt≤C0 for α=2α=2 and Δth−12≤C0Δth−12≤C0 for α=1α=1 where C0C0 is some positive constant. Then, under the stability condition, the optimal H1H1-L2L2 error estimate of the fully discrete solution (unh,Bnh)(unh,Bnh) and optimal L2L2-error estimate of the fully discrete solution pnhpnh are established by using the parabolic dual argument.