This study investigates the effect of shear band evolution on the nucleation of Goss {110}<001> texture during the primary recrystallization of 3.24 wt% Si grain-oriented electrical steel. Nucleation at the early stage of primary recrystallization of the steel is explored both experimentally and numerically. The experimental approach involves cold rolling the steel specimens to obtain a thickness reduction ratio of 76 % and then applying heat treatment to them at 600 °C for less than 1 min. The numerical simulation employes crystal plasticity (CP) finite element model (FEM) to simulate the plastic deformation induced by the dislocation slips on predefined slip systems and non-crystallographic shear bands during cold rolling. Based on the CPFEM results, the generalized strain energy release maximization (GSERM) model is used to predict the preferential orientation probability of recrystallized nuclei for the steel by considering shear band formation. Subsequently, the microstructure evolution during the early stage of primary recrystallization of the steel is simulated using the phase field model (PFM). The developed CP model successfully predicted shear band activation and evolution in the γ-fibers centered on the {111}<112> texture component. The model also demonstrated that shear bands would be the preferred nucleation sites at the early stage of primary recrystallization because of their high stored energy. Moreover, by coupling with the GSERM model, the PFM could reproduce the nucleation of Goss grains at the beginning of primary recrystallization in shear bands.