Efficient Two-Stage Gaussian Process Regression Via Automatic Kernel Search and Subsampling