The Ehrhart polynomial of a matroid specializes to the beta invariant

We show that the linear coefficient of the Ehrhart polynomial of a matroid base polytope evaluated at $t-1$ is equal to, up to normalization, the $\beta$-invariant of the matroid. This yields a lattice-point counting formula for the $\beta$-invariant and establishes a new and unexpected positivity property of Ehrhart polynomials of matroid polytopes.