Algebraic multilevel preconditioners for algebraic problems arising from the discretization of a class of systems of coupled elliptic partial differential equations (PDEs) are presented. These preconditioners are based on modifications of Schwarz methods and of the smoothed aggregation technique, where the coarsening strategy and the restriction and prolongation operators are defined using a point-based approach with a primary matrix corresponding to a single PDE. The preconditioners are implemented in a parallel computing framework and are tested on two representative PDE systems. The results of the numerical experiments show the effectiveness and the scalability of the proposed methods. A convergence theory for the twolevel case is presented