The microstructural and crystallographic aspects, reflected at the macroscopic scale on yield surface and its subsequent evolution, are reappraised by application of crystal plasticity simulations. Strain hardening rule in the slip system is coupled to cope with latent hardening and Bauschinger effect. Uniaxial tension simulation on an isotropic polycrystalline aggregate leads to anisotropic strain hardening. Typical elements of phenomenological plastic anisotropy and hardening rules such as expansion, kinematic shift and distortion of the yield surface, are shown to be featured in crystal plasticity by tuning the slip system hardening rules appropriately.