You can use Ramanujan's master formula to evaluate the Mellin transform of $$ {}_pF_{q} (a_{1},a_{2}, \ldots, a_{p};b_{1},b_{2}, \ldots, b_{q}; -x) = \sum_{k=0}^{\infty} \frac{\Gamma(a_{1}+k) \Gamma(a_{2}+k) \cdots \Gamma(a_{p}+k) \Gamma(b_{1}) \Gamma(b_{2}) \cdots \Gamma(b_{q})}{\Gamma(a_{1}) \Gamma(a_{2}) \cdots \Gamma(a_{p}) \Gamma (b_{1}+k) \Gamma(b_{2}+k) \cdots \Gamma(b_{q}+k)} \frac{(-x)^{k}}{k!}$$
Specifically, $$ \int_{0}^{\infty} x^{s-1} {}_pF_{q} (a_{1},a_{2}, \ldots, a_{p};b_{1},b_{2}, \ldots, b_{q}; -x) \ dx = \Gamma(s) \ \frac{\Gamma(a_{1}-s) \Gamma(a_{2}-s) \cdots \Gamma(a_{p}-s) \Gamma(b_{1}) \Gamma(b_{2}) \cdots \Gamma(b_{q})}{\Gamma(a_{1}) \Gamma(a_{2}) \cdots \Gamma(a_{p}) \Gamma (b_{1}-s) \Gamma(b_{2}-s) \cdots \Gamma(b_{q}-s)}$$
I've used the formula to evaluate the Mellin transforms of functions with simple hypergeometric representations of that form.