Preservers of Gershgorin Set of Jordan Product of Matrices

For $A,B{\in}M_2(\mathbb{C})$, let the Jordan product be AB + BA and G(A) the eigenvalue inclusion set, the Gershgorin set of A. Characterization is obtained for maps ${\phi}:M_2(\mathbb{C}){\rightarrow}M_2(\mathbb{C})$ satisfying $$G[{\phi}(A){\phi}(B)+{\phi}(B){\phi}(A)]=G(AB+BA)$$ for all matrices A and B. In fact, it is shown that such a map has the form ${\phi}(A)={\pm}(PD)A(PD)^{-1}$, where P is a permutation matrix and D is a unitary diagonal matrix in $M_2(\mathbb{C})$.