A rotation of the anisotropic soil fabric pattern is commonly observed in natural slopes with a tilted stratification. This study investigates the rotated anisotropy effects on slope reliability considering spatially varied soils. Karhunen–Loève expansion is used to generate the random fields of the soil shear strength properties (i.e., cohesion and friction angle). The presented probabilistic analyses are based on a meta-model combining Sparse Polynomial Chaos Expansion (SPCE) and Global Sensitivity Analysis (GSA). This method allows the number of involved random variables to be reduced and then the computational efficiency to be improved. Two kinds of deterministic models, namely a discretization kinematic approach and a finite element limit analysis, are considered. A variety of valuable results (i.e., failure probability, probability density function, statistical moments of model response, and sensitivity indices of input variables) can be effectively provided. Moreover, the influences of the rotated anisotropy, autocorrelation length, coefficient of variation and cross-correlation between the cohesion and friction angle on the probabilistic analysis results are discussed. The rotation of the anisotropic soil stratification has a significant effect on the slope stability, particularly for the cases with large values of autocorrelation length, coefficient of variation, and cross-correlation coefficient.