While surrogate-based optimization has encountered a growing success in engineering design, the development of stochastic metamodels, i.e. capable of representing the complete random responses with respect to random inputs, is still an open issue, although they could be fruitfully used in optimization under uncertainty, both with single and multiple objectives. Therefore, the contribution of the paper is twofold. First, hierarchical stochastic metamodels based on moving least squares and spectral decomposition (by polynomial chaos expansion) are proposed in order to obtain a complete description of the random responses with respect to the deterministic and random input parameters. Then, these metamodels are incorporated into a novel multiobjective reliability-based formulation leaning on the concept of probabilistic nondominance. The whole procedure is applied to an analytical test case as well as to the design optimization of space truss structures, demonstrating the ability of the proposed method to provide accurate solutions at an affordable computational time.