We study the effect of primordial isocurvature perturbations on non-Gaussian properties of cosmic microwave background (CMB) temperature anisotropies. We consider generic forms of the non-linearity of isocurvature perturbations which can be applied to a wide range of theoretical models. We derive analytical expressions for the bispectrum and the Minkowski Functionals for CMB temperature fluctuations to describe the non-Gaussianity from isocurvature perturbations. We find that the isocurvature non-Gaussianity in the quadratic isocurvature model, where the isocurvature perturbation S is written as a quadratic function of the Gaussian variable σ, S=σ2−〈σ2〉, can give the same signal-to-noise ratio as fNL= 30 even if we impose the current observational limit on the fraction of isocurvature perturbations contained in the primordial power spectrum . We give constraints on isocurvature non-Gaussianity from Minkowski Functionals using the Wilkinson Microwave Anisotropy Probe (WMAP) 5-year data. We do not find a significant signal of isocurvature non-Gaussianity. For the quadratic isocurvature model, we obtain a stringent upper limit on the isocurvature fraction < 0.070 (95 per cent CL) for a scale-invariant spectrum which is comparable to the limit obtained from the power spectrum.