The product ϕ_{λ}^{(α,β)}(t_{1})ϕ_{λ}^{(α,β)}(t_{2}) of two Jacobi functions is expressed as an integral in terms of ϕ_{λ}^{(α,β)}(t_{3}) with explicit non-negative kernel, when α≧β≧−1/2. The resulting convolution structure for Jacobi function expansions is studied. For special values of α and β the results are known from the theory of symmetric spaces.