Integrally convex functions constitute a fundamental function class in discrete convex analysis. This paper shows that an integer-valued integrally convex function admits an integral subgradient and that the integral biconjugate of an integer-valued integrally convex function coincides with itself. The proof is based on the Fourier-Motzkin elimination. The latter result provides a unified proof of integral biconjugacy for various classes of integer-valued discrete convex functions, including L-convex, M-convex, L 2 -convex, M 2 -convex, BS-convex, and UJ-convex functions as well as multimodular functions. Our results of integral subdifferentiability and integral biconjugacy make it possible to extend the theory of discrete DC (difference of convex) functions developed for L-and M-convex functions to that for integrally convex functions, including an analogue of the Toland-Singer duality for integrally convex functions. ♮ 2 -convex, and M ♮ 2 -convex functions are known to be integrally convex [11]. Multimodular functions [4] are also integrally convex, as pointed out in [13]. Moreover, BS-convex and UJ-convex functions [3] are integrally convex.The concept of integral convexity is used in formulating discrete fixed point theorems and found applications in economics and game theory [6,14,19]. A proximity theorem for integrally convex functions has recently been established in [9] together with a proximityscaling algorithm for minimization. Fundamental operations for integrally convex functions such as projection and convolution are investigated in [8].