For delayed reaction-diffusion Schnakenberg systems with Neumann boundary conditions, critical conditions for Turing instability are derived, which are necessary and sufficient. And existence conditions for Turing, Hopf and Turing-Hopf bifurcations are established. Normal forms truncated to order 3 at Turing-Hopf singularity of codimension 2, are derived. By investigating Turing-Hopf bifurcation, the parameter regions for the stability of a periodic solution, a pair of spatially inhomogeneous steady states and a pair of spatially inhomogeneous periodic solutions, are derived in (τ, ε) parameter plane (τ for time delay, ε for diffusion rate). It is revealed that joint effects of diffusion and delay can lead to the occurrence of mixed spatial and temporal patterns. Moreover, it is also demonstrated that various spatially inhomogeneous patterns with different spatial frequencies can be achieved via changing the diffusion rate. And, the phenomenon that time delay may induce a failure of Turing instability observed by Gaffney and Monk (2006) are theoretically explained.