Because the polynomial is homogeneous in degree (every term is 4th degree), the expression may reduced to two variables by "factoring out" a 4th power of one of the variables.
If we set $X=x/z$ and $Y=y/z$, then the original expression is $z^4$ times the following:
$$ X^3 - X^3Y - Y^3 + Y + XY^3 - X $$
Collecting terms in descending powers of one variable (say $X$), we have:
$$ (1-Y)X^3 + (Y^3 - 1)X - Y^3 + Y $$
By inspection the leading coefficient is a factor:
$$ (1-Y)(X^3 - (Y^2 + Y + 1)X + (Y + 1)Y) $$
Since the second factor is now monic in $X$ and polynomials in $X$ over $\mathbb{Z}[Y]$ form a UFD, we check for factors of the form $X-d$ where root $d$ is a divisor of $(Y+1)Y$. Equivalently we can substitute $X=d$ rather than perform long division, which amounts to synthetic division/Horner's method if you practice it. In any case one gets the factor $X-Y$ rather easily:
$$ (1-Y)(X-Y)(X^2 + YX - (Y+1)) $$
Finally the third factor yields either to inspection, to root testing, or to the quadratic formula:
$$ (1-Y)(X-Y)(X - 1)(X + Y + 1) $$
Multiplying by $z^4$ in such a way that each factor above gets one copy of $z$ converts back to the original variables:
$$ (z-y)(x-y)(x-z)(x + y + z)$$
And we are done!