The following is only a partial answer, but I believe it is relevant to the question.
For $\,n \ge 13\,$ the result follows from a far more general, albeit lesser known, irreducibility criterion of Pólya from Verschiedene Bemerkungen zur Zahlentheorie (1919). Below is the English translation of Prasolov's Russian translation in his book Polynomials.
2.2.3 Irreducibility of polynomials attaining small values
Theorem 2.2.8 (Pólya). Let $\,f\,$ be a polynomial of degree $\,n\,$ with integer coefficients and define $\,m = \left\lfloor\frac{n+1}{2}\right\rfloor\,$. Suppose that, for $\,n\,$ different integers $\,a_1,\ldots,a_n\,$ we have $\,|f(a_i)| \lt 2^{-m}m!\,$ and the numbers $\,a_1,\ldots,a_n\,$ are not roots of $\,f\,$. Then $\,f\,$ is irreducible.
With $\,a_i=i\,$, the criterion implies that all polynomials $\,(x-1)(x-2)\ldots(x-n) \pm k\,$ are irreducible for $\,k = 1, 2, \ldots, \left\lfloor\dfrac{m!}{2^m}\right\rfloor\,$ where $\,m = \left\lfloor\dfrac{n+1}{2}\right\rfloor\,$. For $\,n \ge 13\,$ it is straightforward to show that $\,n \lt \left\lfloor\dfrac{m!}{2^m}\right\rfloor\,$, so the conclusion in OP's question follows.
I can't read German, but the very specific example $\,x(x-1)\ldots(x-n+1) + \lambda n\,$ appears to be mentioned in Pólya's paper as a concrete solution $\lambda=\pm 1$ to Hilbert's existential irreducibility theorem in this case.

An English translation of Pólya's proof (adapted to the present case):
Assume contrariwise that $P_n(x)=g(x) h(x)$ for some $g,h\in\Bbb{Q}[x]$. By Gauss' Lemma we can w.l.o.g. assume that $g,h\in\Bbb{Z}[x]$. One of the factors, say $g(x)$, has degree at least $m\ge\lfloor\dfrac{n+1}2\rfloor$.
Observe that $m+1=\deg g(x)+1\le n$ as otherwise the other factor $h(x)$ is a constant.
By Lagrange interpolation formula we can write
$$
g(x)=\sum_{i=1}^{m+1} g(i) \Delta_i(x),
$$
where for all $i, 0< i\le m+1$,
$$\Delta_i(x)=\frac{\prod_{0< j\le m+1, j\neq i}(x-j)}{\prod_{0< j\le m+1, j\neq i}(i-j)}.$$
The leading coefficient $c_i$ of $\Delta_i(x)$ has absolute value
$$|c_i|=\frac1{i!(m-i)!}=\frac1{m!}\binom m i.$$
Furthermore, $g(i)$ must be a factor of $P_n(i)=n$, so $|g(i)|\le n$.
Putting all these pieces together, the triangle inequality says that the leading coefficient $c=1$ of $g(x)$ has an upper bound
$$
|c|\le\sum_{i=1}^{m+1} |g(i)|\cdot |c_i|\le \frac{n}{m!}\sum_{i=0}^m\binom m i=\frac{n\cdot 2^m}{m!}.
$$
As $c=1$, this is a contradiction whenever $n\cdot 2^m<m!$.