Observational results from strong lensing and dynamical modeling indicate that the total density profiles of early-type galaxies are close to isothermal, i.e. ρtot ∝ rγ with γ ≈ −2. To understand the origin of this universal slope we study a set of simulated spheroids formed in cosmological hydrodynamical zoom-in simulations (see Oser et al. 2010 for more details). We find that the total stellar plus dark matter density profiles of all our simulations on average can be described by a power law with a slope of γ ≈ −2.1, with a tendency towards steeper slopes for more compact, lower mass ellipticals, while the total intrinsic velocity dispersion is flat for all simulations, independent of the values of γ. Our results are in good agreement with observations of Coma cluster ellipticals (Thomas et al. 2007) and results from strong lensing (Sonnenfeld et al. 2012). We find that for z ≳ 2 the majority of the stellar build-up occurs through in-situ star formation, i.e. the gas falls to the center of the galaxy and forms stars, causing the galaxy to be more compact and thus the stellar component to be more dominant. As a result, the total density slopes at z ≈ 2 are generally steeper (around γ ≈ −3). Between z = 2 and z = 0 galaxies grow mostly through dry merging, with each merging event shifting the slope more towards γ ≈ −2. We conclude from our simulations that the steepness of the slope of present day galaxies is a signature of the importance of mostly dry mergers in the formation of an elliptical, and suggest that all elliptical galaxies will with time end up in a configuration with a density slope of γ ≈ −2. For a more detailed analysis with a larger sample of simulations see Remus et al. (2013).