DSolve[Derivative[1][y][x] == (-x + x*(1 + y[x]^(-1)) + x^2*_F1[x*(1 + y[x]^(-1))] - x^2*(1 + y[x]^(-1))*_F1[x*(1 + y[x]^(-1))])^(-1),y[x],x]