As verified structures for landslide mitigation, stabilizing piles are often adopted to treat local failure zones of width-limited soil slopes. To achieve a balance between slope safety and construction cost, an optimized piling scheme for treating width-limited soil slopes should be obtained through three-dimension (3D) stability analysis. In this regard, this paper presents a novel calculation framework for the multi-objective optimization (MOO) design of stabilizing piles. It is based on the first-order reliability method (FORM) and considers a 3D width-limited slope failure with geological uncertainties. The study first develops a deterministic 3D stability model of pile-reinforced slopes using limit analysis. Accounting for soil shear strength uncertainties, reliability analyses of 3D reinforced slopes are conducted based on the prescribed pile-reinforcement patterns. Then, a multi-objective probabilistic design procedure combining the Pareto front and reliability analysis results is proposed. The effectiveness and significance of the proposed MOO design framework are demonstrated through two illustrative examples: one involves designing stabilizing piles in a homogenous slope, and the other involves designing for an inhomogeneous earth slope with depth-dependent soil cohesion. To gain better understanding of the probabilistic impact of uncertain pile design parameters on reinforced slope stability, comprehensive parametric studies are conducted.