Closed-form Inference and Prediction in Gaussian Process State-Space Models