By using a large enough number of particles and implementing a parallel algorithm on the CUDA platform, we have performed brute-force molecular dynamics simulations to study the template-induced heterogeneous crystallization in charged colloids. Six kinds of templates, whose patterns include the planes of fcc(100), fcc(110), fcc(111), bcc(100), bcc(110) and bcc(111), have been implanted into the middle of the simulation box. Except the fcc(111) template, whose structure benefits not only fcc but also hcp crystals resulting in a similar behavior to homogeneous crystallization, bcc-type templates favor the formation of bcc crystals and bcc-like precursors while fcc-type templates favor the formation of fcc crystals and fcc-like precursors. Therefore, for fcc(100) and fcc(110) templates, heterogeneous crystallization will definitely result in a fcc crystallite. However, the results of heterogeneous crystallization that are induced by bcc-type templates are subtly different at different state points. At the state points where the interaction strength of charged colloids is weak and the fcc phase is thermodynamically stable, the bcc crystals formed with the promotion of bcc-type templates are not stable so as to tend to transform into fcc or hcp crystals. When the interaction strength of charged colloids is high, the predominant bcc crystals formed with the promotion of bcc-type templates can always persist within the time scale of simulation although not bcc but fcc crystals are thermodynamically stable.