Here goes (my solution using symmetric polynomials and residues):
First, observe that by an arithmetic argument,
$$P_k(s) = [u^k] \prod_p \left(1 + \sum_{r \geq 1} p^{-rs} u^r\right) = \prod_p \left(1-\frac{u}{p^s}\right)^{-1},$$ where setting $s \mapsto 1$ gives the usual (non-split) Euler product for $\zeta(s)$. We can interpret the RHS product in terms of the complete homogeneous symmetric polynomials for $x_j := p_j^{-s}$:
$$P_k(s) = h_k(x_1,x_2,\ldots).$$ Now by generating function expansions for these polynomials, we have that
$$h_k = [z^k] \exp\left(-\sum_{k \geq 1} \frac{p_k}{k} z^k\right),$$ where $$p_k = \sum_p p^{-ks} = P(sk).$$ Here, $P(s)$ is the prime zeta function.
So by generating function arithmetic, we have that
$$\zeta(s)_N = [z^N] \frac{1}{1-z} \exp\left(-\sum_{k \geq 1} \frac{p_k}{k} z^k\right).$$ The RHS has poles at $s := 1/k$ for $k \geq 1$ by the nature of singularities of $P(s)$. From an inverse Z-transform, we can recover these coefficients as
$$\begin{align}
\zeta(s)_N & = \frac{1}{2\pi\imath} \oint_C \frac{s^{N-1}}{1-s} \exp\left(-\sum_{k \geq 1} \frac{P(sk)}{k} s^k\right) ds \\
& = \sum_{k \geq 1} \operatorname{Res}_{s=\frac{1}{k}}\left[\frac{s^{N-1}}{1-s} \exp\left(-\sum_{k \geq 1} \frac{P(sk)}{k} s^k\right)\right].
\end{align}$$
Since we can expand $P(s)$ near $s \rightarrow 1^{+}$ by
$$P(1+\varepsilon) \sim -\log(\varepsilon) + B + O(\varepsilon),$$ my back-of-the-envelope calculations suggest that these residues correspond to (NOTE: I haven't checked these out numerically)
$$R_k = \begin{cases}
-e^{-B} \exp\left(-\sum\limits_{j \geq 2} \frac{P(j)}{j}\right), & k = 1 \\
\frac{1}{k^N (k-1)} \exp(-B / k^{k+1}) \exp\left(-\sum\limits_{\substack{j \geq 1 \\ j \neq k}} \frac{P(kj)}{j \cdot k^j}\right), & k \geq 2.
\end{cases}$$ The constant $B$ in the previous equations is defined by
$$B = \sum_{n \geq 2} \frac{\mu(n)}{n} \log \zeta(n) \approx -0.315718452.$$
Now you should be able to proceed with evaluating zeros for the function.
N.b.: With respect to there not being much in the way of literature on $P(s)$, my friend has pointed out to me offline that the prime zeta function is essentially meaningless due to the essential and messy, complicated structure of its poles. It cannot be analytically continued past zero, and hence, is not nicely modular, nor does it have a nice functional equation like most zeta functions do. It's easier just not to work with the function unless you are forced.
$\sum_{n\le x, \Omega(n) = k} 1 \approx \frac{x}{\log x} \frac{(\log \log x)^k}{k!} \approx \frac1k \sum_{p \le x} P_{k-1}(x/p) \approx \frac1k\sum_{2\le n\le x} P_{k-1}(x/n)\frac{1}{\log n}$.
– reuns Mar 10 '19 at 23:53