In this article we obtain estimates of Neumann eigenvalues of p-Laplace operators in a large class of space domains satisfying quasihyperbolic boundary conditions. The suggested method is based on composition operators generated by quasiconformal mappings and their applications to Sobolev-Poincaré-inequalities. By using a sharp version of the inverse Hölder inequality we refine our estimates for quasi-balls, that is, images of balls under quasiconformal mappings of the whole space.