Computationally Efficient Bayesian Learning of Gaussian Process State Space Models