Let Φ : 𝔖(M1) → 𝔖(M2) be a bijection (not assumed affine nor continuous) between the sets of normal states of two quantum systems, modelled on the self-adjoint parts of von Neumann algebras M1 and M2, respectively. This paper concerns with the situation when Φ preserves (or partially preserves) one of the following three notions of “transition probability” on the normal state spaces: the transition probability PU introduced by Uhlmann [Rep. Math. Phys. 9, 273-279 (1976)], the transition probability PR introduced by Raggio [Lett. Math. Phys. 6, 233-236 (1982)], and an “asymmetric transition probability” P0 (as introduced in this article). It is shown that the two systems are isomorphic, i.e., M1 and M2 are Jordan ∗-isomorphic, if Φ preserves all pairs with zero Uhlmann (respectively, Raggio or asymmetric) transition probability, in the sense that for any normal states μ and ν, we have PΦ(μ),Φ(ν)=0 if and only if P(μ, ν) = 0, where P stands for PU (respectively, PR or P0). Furthermore, as an extension of Wigner’s theorem, it is shown that there is a Jordan ∗-isomorphism Θ : M2 → M1 satisfying Φ = Θ∗|𝔖(M1) if and only if Φ preserves the “asymmetric transition probability.” This is also equivalent to Φ preserving the Raggio transition probability. Consequently, if Φ preserves the Raggio transition probability, it will preserve the Uhlmann transition probability as well. As another application, the sets of normal states equipped with either the usual metric, the Bures metric or “the metric induced by the self-dual cone,” are complete Jordan ∗-invariants for the underlying von Neumann algebras.