Bayesian, Gradient-Free, and Multi-Fidelity Supervised Dimension Reduction Methods for Surrogate Modeling of Expensive Analyses with High-Dimensional Inputs