This paper is concerned with monotonic and load reversal deformation of four dual-phase (DP) 590, 780, 980, 1180 and one martensitic (MS) 1700 steel sheets. While the monotonic data are presented for all steels, the load reversal data are provided for DP 590, DP 780, and DP 1180. Particularities pertaining to the reversal deformation including the decreasing hardening rate during forward tension, a linear and then a non-linear unloading, followed by the Bauschinger effect, and a shift in the hardening rate during continuous straining were quantified and discussed as a function of loading history. Moreover, parameters such as the reloading stress differential, reloading softening stress, ratcheting strain, and unloading deviation stress were determined and discussed as a function of martensite fraction. The data were interpreted and predicted using an elasto-plastic self-consistent (EPSC) crystal plasticity model incorporating anisotropic elasticity, a dislocation density based hardening law, and a slip system backstress law. The model parameters associated with the slip strengths of ferrite and martensite and backstress were established. This work demonstrated the ability of crystal plasticity modeling to account for the co-dependent nature of crystallographic slip in ferrite and martensite and the sources of hardening caused by history-dependent dislocation density evolution and backstress to predict not only monotonic but also hysteresis in plastic response during the forward-reversal cycles. The combination of comprehensive experimental data and modeling results allowed us to infer that the tradeoff between the magnitude of backstress per phase and the volume fraction of ferrite versus martensite per steel governs the unloading and subsequent yielding per steel, while the dissolution of dislocations facilitates capturing the hardening rates during load reversal deformation.