The supersymmetric approach is used to analyse a class of two-dimensional quantum systems with periodic potentials. In particular, the method of SUSY-separation of variables allowed us to find a part of the energy spectra and the corresponding wave functions (partial solvability) for several models. These models are not amenable to conventional separation of variables, and they can be considered as two-dimensional generalizations of Lamé, associated Lamé, and trigonometric Razavy potentials. All these models have the symmetry operators of fourth order in momenta, and one of them (the Lamé potential) obeys the property of self-isospectrality. © 2007 Elsevier B.V. All rights reserved.