Summary: In this paper we study general nonlinear stochastic differential equations, where the usual Brownian motion is replaced by a Lévy process. We also suppose that the coefficient multiplying the increments of this process is merely Lipschitz continuous and not necessarily linear in the time-marginals of the solution as is the case in the classical McKean-Vlasov model. We first study existence, uniqueness and particle approximations for these stochastic differential equations. When the driving process is a pure jump Lévy process with a smooth but unbounded Lévy measure, we develop a stochastic calculus of variations to prove that the time-marginals of the solutions are absolutely continuous with respect to the Lebesgue measure. In the case of a symmetric stable driving process, we deduce the existence of a function solution to a nonlinear integro-differential equation involving the fractional Laplacian.