The equation of state of nuclear matter is of great importance in the interpretation of heavy-ion reactions and a crucial ingredient in many astrophysical models and simulations. One of the challenges of nuclear theory is a microscopic derivation of the equation of state starting from the bare nucleon-nucleon interaction in the vacuum. In order to achieve that, one has to properly take into account the modifications of the single-particle properties in the medium and the correlations induced by the nucleon-nucleon interaction. After introducing the general aspects of the problem and the state of the art of microscopic calculations of the equation of state, I will show results obtained with the self-consistent Green’s function method implemented with two- and three-nucleon forces. The effect of three-body forces on thermodynamic observables and single-particle properties will be analyzed. Finally, the connection with approaches using a low-momentum potential and energy density functionals will be addressed.