Symmetry-protected topological phases (SPTs) are phases of matter lying outside the Landau paradigm: they do not break any symmetries and hence cannot be detected by any local order parameters. Although SPTs are characterized by their support of gapless edge modes, it is natural to expect that information about which SPT phase a wavefunction belongs to should also be encoded in the bulk. How to extract this information for SPTs in 2D has, to date, been largely unclear. In this talk, I will close this gap and propose two non-local order parameters which can detect and distinguish all bosonic SPTs in 2D protected by Abelian, internal symmetries. Our order parameters take as input the reduced density matrix on a disk and extract the corresponding topological invariants by simulating the SPT path integral on a non-trivial spacetime manifold. Remarkably, these order parameters correspond to symmetry-twisted entanglement measures and imply 2D SPTs are characterized by four-party and six-party entanglement, suggesting symmetry-protected multipartite entanglement is a defining feature of SPT order.