Making use of the generalized $(p, q)$-Bernardi integral operator, we introduce and study a new class $\mathcal{F J}_{p, q}^m(\alpha, \delta, \lambda, \gamma)$ of multivalent analytic functions with negative coefficients in the open unit disk $E$. Several geometric characteristics are obtained, like, coefficient estimate, radii of convexity, close-to-convexity and starlikeness, closure theorems, extreme points, integral means inequalities, neighborhood property and convolution properties for functions belonging to the class $\mathcal{F} \mathcal{J}_{p, q}^m(\alpha, \delta, \lambda, \gamma)$.