Non-radial solutions to a bi-harmonic equation with negative exponent

We prove the existence of non-radial entire solution to $$\Delta^2 u+u^{-q}=0\quad\text{in }\mathbb{R}^3,\quad u>0,$$ for $q>1$. This answers an open question raised by P. J. McKenna and W. Reichel (E. J. D. E. \textbf{37} (2003) 1-13).