This text deals with multidimensional Borg-Levinson inverse theory. Its main purpose is to establish that the Dirichlet eigenvalues and Neumann boundary data of the Dirichlet Laplacian acting in a bounded domain of dimension 2 or greater, uniquely determine the real-valued bounded potential. We first address the case of incomplete spectral data, where finitely many boundary spectral eigen-pairs remain unknown. Under suitable summability condition on the Neumann data, we also consider the case where only the asymptotic behavior of the eigenvalues is known. Finally, we use the multidimensional Borg-Levinson theory for solving parabolic inverse coefficient problems.