For visualisation, the following code was run in Mathematica.
ContourPlot[{Abs[x + I y - 1/(x + I y)] == 2}, {x, -2.42, 2.42}, {y, -2.42, 2.42}, Axes -> True]
Here is the locus plot generated by the above code.
I do not see a nice way of showing the locus is what is is apart from brute force algebraic factorisation, which basically requires knowing the form of the locus a priori.
A less cumbersome approach is desired, if such a method exists.

We can bring a little geometry into the picture.
First render
$|z+1|^2|z-1|^2=4|z|^2$
by algebraic methods as explained in other answers. Now draw the quadrilateral whose vertices are $\pm 1, \pm z$. This is a parallelogram, meaning the sum of its squared diagonal equals the sum of all four squared sides. Thereby
$4|z|^2+4=2(|z+1|^2+|z-1|^2)$
Eliminating $|z|^2$ between these two equations and collecting all terms to one side then leads to:
$|z+1|^2|z-1|^2-2(|z+1|^2+|z-1|^2)+4=0$
which is factorable!
$\color{blue}{(|z+1|^2-2)(|z-1|^2-2)=0}$
The union of circles then follows.