Identification of faulty variables is an important component of multivariate statistical process monitoring (MSPM); it provides crucial information for further analysis of the root cause of the detected fault. The main challenge is the large number of combinations of process variables under consideration, usually resulting in a combinatorial optimization problem. This paper develops a generic reconstruction based multivariate contribution analysis (RBMCA) framework to identify the variables that are the most responsible for the fault. A branch and bound (BAB) algorithm is proposed to efficiently solve the combinatorial optimization problem. The formulation of the RBMCA does not depend on a specific model, which allows it to be applicable to any MSPM model. We demonstrate the application of the RBMCA to a specific model: the mixture of probabilistic principal component analysis (PPCA mixture) model. Finally, we illustrate the effectiveness and computational efficiency of the proposed methodology through a numerical example and the benchmark simulation of the Tennessee Eastman process.