MCMC for variationally sparse Gaussian processes