I think the path-integral is a complete red herring here! I'll try to convince you that Wick rotation yields completely equivalent way of writing the Lagrangian in classical field theory.
Consider a classical action
$$S[x] = \int L[x(t)] dt$$
where $x:\mathbb{R} \to \mathcal{M}$ for some target manifold $\mathcal{M}$. The Lagrangian is schematically given by
$$L[x(t)] = \left(\left.\frac{dx(s)}{ds}\right|_{s=t}\right)^2-V(x(t))$$
where $V(x(t))$ is some polynomial in $x(t)$ which (critically) involves no derivatives.
Now analytically continue $x$ to a function $\tilde{x}: \mathbb{C}\to\mathcal{M}$ by defining $\tilde{x}(c) = x(|c|)$ which is obviously analytic. Relabel $\tilde{x}$ as $x$ for simplicity. Define a new variable
$$\tau = it$$
inside the integral, and substitute. (Warning: there are mathematical subtleties about complex substitutions, which should be dealt with using Jordan's lemma). Ignoring the subtleties, the resulting integral is
$$S[x] = \int L[x(-i\tau)](-i)d\tau$$
Now let's examine what happens to the Lagrangian more carefully. Looking at the first term we have
$$\left(\left.\frac{dx(s)}{ds}\right|_{s=-i\tau}\right)^2$$
Change the differentiation variable to $u = is$ and this term becomes
$$\left(i\left.\frac{dx(u)}{du}\right|_{u=\tau}\right)^2$$
Relabelling $u\to s$ we see that the first term has the same form as originally, but with $t$ replaced by $\tau$ and an extra minus sign, viz.
$$-\left(\left.\frac{dx(s)}{ds}\right|_{s=\tau}\right)^2$$
Now on to the potential term. This is much simpler because $x(-i\tau)=x(|-i\tau|)=x(\tau)$ by definition so the potential term is just
$$V(x(\tau))$$
which is of exactly the same form as originally. Now we define a Euclidean Lagrangian
$$L_E(x(\tau)) = \left(\left.\frac{dx(s)}{ds}\right|_{s=\tau}\right)^2+V(x(\tau))$$
Putting it all together we find
$$S[x] = i\int L_E[x(\tau)]d\tau$$
Finally defining
$$S_E[x] = \int L[x(t)]dt$$
we see that it's mathematically equivalent to calculate the path integral as
$$\int \mathcal{D}x\exp(iS[x]) \textrm{ or } \int \mathcal{D}x\exp(-S_E[x])$$
The final condition for a realistic theory in physics is that its predictions (of scattering amplitudes or correlators etc.) have to agree with observations. This implies that the predictions have to be consistent and obey some general consistency conditions (unitarity, non-negativity of probabilities, some symmetries, locality or approximate locality, and so on).
For theoretical theories, it's just the general consistencies that hold. Now, the task is to classify all possible theories and learn how to calculate with them.
It turns out that the Euclidean spacetime or world sheet is simply a simpler, more straightforward, more free-of-subtleties approach to produce a machine that calculates some scattering amplitudes or other observables.
At least formally, the Euclidean theories may be continued to analytic ones and vice versa. For nontrivial spacetime topologies, the Euclidean objects are likely to be more manageable. For example, the world sheets in string theory (think about a torus or pants diagrams etc.) are much more well-behaved in the Euclidean signature so we may consider this approach "primary". Covariant calculations in string theory are almost always done with the Euclidean ones. The result may be continued to the Minkowski momenta etc. and some of the consistency conditions above are still guaranteed to hold because of some properties of the complex calculus.
In the light-cone gauge, we may work directly with the Minkowski-signature world sheets. But we pay the price that the interaction points where strings split or join are singular and the direction of the "future time" is ambiguous. We must also include contact terms, higher-order interaction terms, to deal with some divergences caused by the singular world sheets, but when these things are summed over, we may prove that the resulting amplitudes agree with the covariantly computed ones (in the Euclidean signature).
Gravity in $d\geq 4$ (and maybe 3) suffers from the "negative norm conformal factor". The Euclideanized Einstein-Hilbert action $\int R \sqrt{g}$ is no longer positively definite. In particular, if you consider scalar waves that scale the metric by an overall number, $g_{\mu\nu}=e^F\eta_{\mu\nu}$, and derive the kinetic term for $F$, it will have the opposite sign than the kinetic term for other components of the metric tensor (the physical polarizations of the gravitational waves, like $g_{xy}$).
It follows that the action will be bounded neither from below nor from above, and $\exp(-S_E)$ in the Euclidean path integral will diverge in some region of the configuration space. In this sense, people believe that the Minkowskian path integral must be the "more kosher one" for higher-dimensional gravity. But this is a bit empty statement because at the quantum level, higher-dimensional gravity obtained as a direct quantization of Einstein's equations is inconsistent, anyway. And string theory which is consistent and contains gravity doesn't give us any tool to directly rewrite the path integral in terms of spacetime fields including the metric; it is not a field theory in the ordinary sense. So the preference for the "Minkowski signature" is a bit vacuous. After all, the Minkowskian action isn't bounded from either side, either. This is considered "not to be a problem" because the integrand is $\exp(iS)$ which still has the absolute value equal to one, so it doesn't diverge. But I would personally say that the unboundedness of the Euclidean action is the "same" problem for the Minkowskian path integral.
Quite generally, the Wick rotation is extremely important in quantum field theory and it is actually even more important in quantum gravity or places with many spacetime (or world sheet) topologies, i.e. in situations where one has many different "time variables" in which we might try to expand things. One shouldn't be afraid but at the end, whatever theory he deals with, he gets some amplitudes whose self-consistency (and/or consistency with observations) must be verified. With some "good rules of behavior" while Wick-rotating, one may be "pretty sure" that some tests will be passed.
Best Answer
1st comment:
It's worth thinking for a second about where Wick rotation comes from. You can do this in the context of the quantum mechanics of a free particle. In QFT, all of the details are more complicated, but the basic idea is the same.
In free particle, QM, we get the path integral by inserting sums over intermediate states at various times. The need for Wick rotation arises as soon as you do this just once.
$\langle q' | e^{- \frac{ iP^2 t}{2m\hbar}}|q\rangle = \int_{-\infty}^\infty \langle q'| p \rangle \langle p |e^{- \frac{ iP^2 t}{2m\hbar}}|q\rangle dp = \frac{1}{2\pi\hbar} \int_{-\infty}^\infty e^{\frac{-i t }{2m\hbar} p^2 + i\frac{q' - q}{\hbar} p} dp$.
This is an oscillatory integral. The integrand has norm 1 because the argument of the exponential is purely imaginary. Such integrals don't converge absolutely, so the right hand side of this equation is not obviously well-defined. It's not Lebesgue integralable, although it is convergent as a Riemann integral, thanks to some rather delicate cancellations. To make the integral well defined -- equivalently to see how these cancellations occur -- we need to supply some additional information.
Wick rotation provides a way of doing this. You observe that the left hand side is analytic in $t$, and that the right hand side is well-defined if $Im(t) < 0$. Then you can define the integral for real $t$ by saying that it's analytic continued from complex $t$ with negative imaginary part.
2nd comment:
As V. Moretti pointed out, in QFT, it's in some sense backwards to think of analytically continuing from Minkowski signature to Euclidean signature. Rather, one finds something in Euclidean signature which has nice properties and then analytically continues from Euclidean to Minkowski. However, one can often begin this process by taking a Minkowski action and finding its Euclidean version, and then trying to build up a QFT from there. There's no guarantee that this will work though. Spinor fields may have reality conditions that depend on the signature of spacetime. Or the Euclidean action you derive may be badly behaved. This is famously the case for Einstein's gravity; the Euclidean action is not bounded below, so one does not get a sensible Euclidean theory.