Recall the archetypal Frullani integral
$$\int_{0}^{\infty} \frac{e^{-u} - e^{-nu}}{u} \, du = \log n. $$
Since the integrand is non-negative for all $n \geq 1$, when $x \in [0, 1)$ we can apply the Tonelli's theorem to interchange the sum and integral unconditionally to get
\begin{align*}
\sum_{n=1}^{\infty} x^n \log n
&= \int_{0}^{\infty} \left( \sum_{n=1}^{\infty} x^n \cdot \frac{e^{-u} - e^{-nu}}{u} \right) \, du \\
&= \frac{x^2}{1-x} \int_{0}^{\infty} \frac{e^{-u}(1 - e^{-u})}{u(1-xe^{-u})} \, du. \tag{1}
\end{align*}
Notice that the last integral converges absolutely. So this computation can be fed back to the Fubini's theorem, showing that exactly the same computation can be carried out to prove $\text{(1)}$ for all $|x| < 1$.
Now we would like to take limit as $x \to -1^{+}$. When $x \in (-1, 0]$, the integrand of the last integral of $\text{(1)}$ is uniformly bounded by the integrable function $u^{-1}e^{-u}(1-e^{-u})$. Therefore by the dominated convergence theorem, as $x \to -1^{+}$ we have
$$
\lim_{x\to -1^+} \sum_{n=1}^{\infty} x^n \log n
= \frac{1}{2} \int_{0}^{\infty} \frac{e^{-u}(1 - e^{-u})}{u(1+e^{-u})} \, du.
$$
This already proves that the limit exists, but it even tells more that the limit is indeed $\eta'(0)$. To this end, we first perform integration by parts to remove the pesky factor $u$ in the denominator. Then the right-hand side becomes
$$ \frac{1}{2} \int_{0}^{\infty} \frac{e^{-u}(1 - e^{-u})}{u(1+e^{-u})} \, du
= \int_{0}^{\infty} \left( \frac{e^u}{(e^u + 1)^2} - \frac{e^{-u}}{2} \right) \log u \, du$$
In order to compute this integral, it suffices to prove the following claim.
Claim. We have
$$ \int_{0}^{\infty} e^{-u} \log u = -\gamma, \qquad \int_{0}^{\infty} \frac{e^u \log u}{(e^u + 1)^2} \, du = -\frac{\gamma}{2} + \eta'(0). $$
Notice that the first claim is an immediate consequence of the identity $\psi(1) = -\gamma$, where $\psi$ is the digamma function. Next, term-wise integration gives
$$ \int_{0}^{\infty} \frac{u^{s-1}}{e^{\alpha u} + 1} \, du = \alpha^{-s} \Gamma(s)\eta(s) $$
where $\alpha > 0$ and $s$ is initially assumed to satisfy $\Re(s) > 1$ (so that interchanging the summation and integration works smoothly). Differentiating both sides w.r.t. $\alpha$ gives
$$ \int_{0}^{\infty} \frac{u^s e^{\alpha u}}{(e^{\alpha u} + 1)^2} \, du = \alpha^{-s-1} \Gamma(s+1)\eta(s). $$
Although we initially assumed $\Re(s) > 1$, now both sides define an analytic function for $\Re(s) > -1$, hence by the principle of analytic continuation this identity extends to this region as well. Now plugging $\alpha = 1$ and differentiating both sides w.r.t. $s$, we get
$$ \int_{0}^{\infty} \frac{u^s e^u \log u}{(e^u + 1)^2} \, du = \Gamma(s+1)\psi(s+1)\eta(s) + \Gamma(s+1)\eta'(s). $$
Plugging $s = 0$ and using known values $\psi(1) = -\gamma$ and $\eta(0) = \frac{1}{2}$, this yields
$$ \int_{0}^{\infty} \frac{e^u \log u}{(e^u + 1)^2} \, du = -\frac{\gamma}{2} + \eta'(0). $$