I'm not entirely sure how to show $PQ + I$ is invertible in infinite dimensions (Edit: a proof of this fact following @MartinArgerami 's arguments has been added to the end of this answer), but assuming that isn't an issue, the following provides a proof that your proposition is correct:
Since $Q$ is strictly positive, $Q^{-1/2}$ is a densely-defined, closed operator. Then $PQ + I = Q^{-1/2}Q^{1/2}PQ^{1/2}Q^{1/2} + Q^{-1/2}Q^{1/2} = Q^{-1/2}(Q^{1/2}PQ^{1/2} + I)Q^{1/2}$. Observe that $Q^{1/2}PQ^{1/2}$ is positive, hence $Q^{1/2}PQ^{1/2} + I$ is invertible. As such, we may consider the operator $Q^{-1/2}(Q^{1/2}PQ^{1/2} + I)^{-1}Q^{1/2}$. We see that it is closed and defined on the range of $Q^{-1/2}(Q^{1/2}PQ^{1/2} + I)Q^{1/2} = PQ + I$, which, as $PQ + I$ is invertible, is the entire $l^2$. Hence, it is bounded and it is then easy to see it is exactly the inverse of $Q^{-1/2}(Q^{1/2}PQ^{1/2} + I)Q^{1/2} = PQ + I$. (In finite dimensions this produces a proof that $PQ + I$ is invertible. Domain issues in infinite dimensions seem to complicate things.) Hence, $Q(PQ + I)^{-1} = QQ^{-1/2}(Q^{1/2}PQ^{1/2} + I)^{-1}Q^{1/2} = Q^{1/2}(Q^{1/2}PQ^{1/2} + I)^{-1}Q^{1/2}$. The same, of course, holds for $Q_k$ and $P_k$, so we need to show $\lim_{k \rightarrow \infty} \mathrm{Tr}(Q_k^{1/2}(Q_k^{1/2}P_kQ_k^{1/2} + I)^{-1}Q_k^{1/2}) = \mathrm{Tr}(Q^{1/2}(Q^{1/2}PQ^{1/2} + I)^{-1}Q^{1/2})$.
From this point onwards the best way I know is through functional calculus. We first note the following: For any fixed polynomial $f$ with $f(0) = 0$, since both $P$ and $Q$ are compact, we see that $Q_k^{1/2}P_kQ_k^{1/2} \rightarrow Q^{1/2}PQ^{1/2}$ in norm, hence $f(Q_k^{1/2}P_kQ_k^{1/2}) \rightarrow f(Q^{1/2}PQ^{1/2})$ in norm. Now, we note that $Q_k^{1/2}P_kQ_k^{1/2}$ are uniformly bounded in norm, so there exists a sequence of polynomials $g_n$ which approximates $(\cdot + I)^{-1}$ in a uniform way, i.e., for each $\epsilon > 0$, there exists an $N$ s.t. for all $n \geq N$, $||g_n(Q_k^{1/2}P_kQ_k^{1/2}) - (Q_k^{1/2}P_kQ_k^{1/2} + I)^{-1}|| \leq \epsilon$ for all $k$ and furthermore $||g_n(Q^{1/2}PQ^{1/2}) - (Q^{1/2}PQ^{1/2} + I)^{-1}|| \leq \epsilon$. Since $(0 + 1)^{-1} = 1$, we may choose $g_n$ s.t. $g_n(0) = 1$. Furthermore, such $g_n$ may be chosen so that they are uniformly bounded on the spectra of all $Q_k^{1/2}P_kQ_k^{1/2}$ as well as $Q^{1/2}PQ^{1/2}$, i.e., there exists $C' > 0$ s.t. $||g_n(Q_k^{1/2}P_kQ_k^{1/2})||, ||g_n(Q^{1/2}PQ^{1/2})|| \leq C'$ for all $n, k$. Let $f_n = g_n - 1$, which are then polynomials satisfying $f_n(0) = 0$. Let $C = C' + 1$, we observe that $||f_n(Q_k^{1/2}P_kQ_k^{1/2})||, ||f_n(Q^{1/2}PQ^{1/2})|| \leq C$ for all $n, k$.
Choose a large enough $N$ satisfying the condition in the preceding paragraph. Now, for any $k$,
$$|\mathrm{Tr}(Q_k^{1/2}(Q_k^{1/2}P_kQ_k^{1/2} + I)^{-1}Q_k^{1/2}) - (\mathrm{Tr}(Q_kf_N(Q_k^{1/2}P_kQ_k^{1/2})) + \mathrm{Tr}(Q_k))|$$
$$= |\mathrm{Tr}(Q_k(Q_k^{1/2}P_kQ_k^{1/2} + I)^{-1}) - \mathrm{Tr}(Q_k(f_N(Q_k^{1/2}P_kQ_k^{1/2}) + I))|$$
$$= |\mathrm{Tr}(Q_k(Q_kP_kQ_k^{1/2} + I)^{-1}) - \mathrm{Tr}(Q_kg_N(Q_k^{1/2}P_kQ_k^{1/2}))|$$
$$= |\mathrm{Tr}(Q_k((Q_kP_kQ_k^{1/2} + I)^{-1} - g_N(Q_k^{1/2}P_kQ_k^{1/2})))|$$
$$\leq \mathrm{Tr}(Q_k)||(Q_kP_kQ_k^{1/2} + I)^{-1} - g_N(Q_k^{1/2}P_kQ_k^{1/2})||$$
$$\leq \mathrm{Tr}(Q)\epsilon$$
Where in the final step we used $\mathrm{Tr}(Q_k) \leq \mathrm{Tr}(Q)$. Similarly,
$$|\mathrm{Tr}(Q^{1/2}(Q^{1/2}PQ^{1/2} + I)^{-1}Q^{1/2}) - (\mathrm{Tr}(Qf_N(Q^{1/2}PQ^{1/2})) + \mathrm{Tr}(Q))| \leq \mathrm{Tr}(Q)\epsilon$$
Since $f_N$ is a polynomial satisfying $f_N(0) = 0$, by what we have proved before we have $f_N(Q_k^{1/2}P_kQ_k^{1/2}) \rightarrow f_N(Q^{1/2}PQ^{1/2})$ in norm. We also have $\mathrm{Tr}(Q_k) \rightarrow \mathrm{Tr}(Q)$, so for large $k$ we have $||f_N(Q_k^{1/2}P_kQ_k^{1/2}) - f_N(Q^{1/2}PQ^{1/2})|| \leq \epsilon$ and $|\mathrm{Tr}(Q_k) - \mathrm{Tr}(Q)| \leq \epsilon$. Thus, for large $k$,
$$|\mathrm{Tr}(Q_kf_N(Q_k^{1/2}P_kQ_k^{1/2})) - \mathrm{Tr}(Qf_N(Q^{1/2}PQ^{1/2}))|$$
$$\leq |\mathrm{Tr}(Q_k(f_N(Q_k^{1/2}P_kQ_k^{1/2}) - f_N(Q^{1/2}PQ^{1/2})))| + |\mathrm{Tr}((Q_k - Q)f_N(Q^{1/2}PQ^{1/2}))|$$
$$\leq \mathrm{Tr}(Q_k)||f_N(Q_k^{1/2}P_kQ_k^{1/2}) - f_N(Q^{1/2}PQ^{1/2})|| + |\mathrm{Tr}(Q_k) - \mathrm{Tr}(Q)|||f_N(Q^{1/2}PQ^{1/2})||$$
$$\leq \mathrm{Tr}(Q)\epsilon + C\epsilon$$
Hence, putting everything together, we have, for large $k$,
$$|\mathrm{Tr}(Q_k^{1/2}(Q_k^{1/2}P_kQ_k^{1/2} + I)^{-1}Q_k^{1/2}) - \mathrm{Tr}(Q^{1/2}(Q^{1/2}PQ^{1/2} + I)^{-1}Q^{1/2})|$$
$$\leq |\mathrm{Tr}(Q_k^{1/2}(Q_k^{1/2}P_kQ_k^{1/2} + I)^{-1}Q_k^{1/2}) - (\mathrm{Tr}(Q_kf_N(Q_k^{1/2}P_kQ_k^{1/2})) + \mathrm{Tr}(Q_k))| + |\mathrm{Tr}(Q^{1/2}(Q^{1/2}PQ^{1/2} + I)^{-1}Q^{1/2}) - (\mathrm{Tr}(Qf_N(Q^{1/2}PQ^{1/2})) + \mathrm{Tr}(Q))| + |\mathrm{Tr}(Q_kf_N(Q_k^{1/2}P_kQ_k^{1/2})) - \mathrm{Tr}(Qf_N(Q^{1/2}PQ^{1/2}))| + |\mathrm{Tr}(Q_k) - \mathrm{Tr}(Q)|$$
$$\leq \mathrm{Tr}(Q)\epsilon + \mathrm{Tr}(Q)\epsilon + \mathrm{Tr}(Q)\epsilon + C\epsilon + \epsilon$$
$$= (3\mathrm{Tr}(Q) + C + 1)\epsilon$$
Letting $\epsilon \rightarrow 0$ concludes the proof.
A proof that $PQ + I$ is invertible in infinite dimensions following @MartinArgerami 's arguments:
We show that $\sigma(PQ) \subseteq [0, \infty)$, which is sufficient since then $-1 \notin \sigma(PQ)$, so by definition of the spectrum we have $PQ + I$ is invertible.
We note that $PQ = P^{1/2}(P^{1/2}Q)$. We note the following general fact: In a unital algebra (for example $\mathbb{B}(l^2)$), for any $a, b$ in the algebra, $\sigma(ab)\setminus\{0\} = \sigma(ba)\setminus\{0\}$. Indeed, by scaling and symmetry, it is sufficient to show that $1 - ab$ is invertible implies $1 - ba$ is invertible. A direct calculation shows that if $c$ is the inverse of $1 - ab$, then $1 + bca$ is the inverse of $1 - ba$. (This proof can be found in Remark 1.2.1 in $C^*$ -algebras and operator theory by Gerard Murphy.) Applying this fact to $a = P^{1/2}, b = P^{1/2}Q$, we see that $\sigma(PQ) \subseteq \{0\} \cup \sigma(P^{1/2}QP^{1/2})$. Since $P^{1/2}QP^{1/2} \geq 0$, $\sigma(P^{1/2}QP^{1/2}) \subseteq [0, \infty)$. This proves the claim.