An operator performing the mapping is $$O= D/(e^D-1)=e^{B.(0)D},$$

where $D=d/dx $ and $(B.(0))^n=B_n(x)|_{x=0}$, since the Bernoulli polynomials are an Appell sequence.

Edit (6/20/2017):

This operator is essentially the Todd operator. See the discussions on pg. 30 and Appendix B of "[Permutohedra, associahedra, and beyond][1]" by Postnikov of the Todd operator as a transform of the homogeneous volume polynomials for classes  of polytopes into a generalized Ehrhart polynomial coding the number of lattice points in the polytopes.


  [1]: https://arxiv.org/abs/math/0507163