If I have some time-dependent interaction Hamiltonian $H_I(\tau)$, with the interaction beginning in the asymptotic past ($\tau$$\rightarrow$$-\infty$) and ending in the asymptotic future ($\tau$$\rightarrow$$\infty$), then the first two terms in the Dyson series will read:
\begin{equation} U=U^{(1)}+U^{(2)}+\cdots \end{equation}
where
\begin{equation} U^{(1)}=-i\int_{-\infty}^{\infty}d\tau H_I(\tau) \end{equation} \begin{equation} U^{(2)}=-\int_{-\infty}^{\infty}d\tau\int_{-\infty}^{\infty}d\tau'\Theta(\tau-\tau')H_I(\tau)H_I(\tau'). \end{equation}
The interaction Hamiltonian is Hermitian $H_I^{\dagger}(\tau)=H_I(\tau)$. I want to calculate the first two terms of the series expansion of $U^{\dagger}$. I know that
\begin{equation} U^{(1)\dagger}=i\int_{-\infty}^{\infty}d\tau H_I(\tau). \end{equation}
My question is, when I write down $U^{(2)\dagger}$, do I reverse the time ordering (i.e, replace $\Theta(\tau-\tau')$ with $\Theta(\tau'-\tau)$), and write down
\begin{equation} U^{(2)\dagger}=-\int_{-\infty}^{\infty}d\tau\int_{-\infty}^{\infty}d\tau'\Theta(\tau'-\tau)H_I(\tau)H_I(\tau'), \end{equation}
or do I simply keep the term the same (since there is no factor of $i$. and the Hamiltonians are already Hermitian)?