I already proved that $\zeta(z)=\frac{1}{\Gamma(z)}\int_0^\infty\frac{t^{z-1}}{e^t-1}dt=\frac{\Gamma(z-1)}{2\pi i}\int_{-\infty}^0\frac{t^{z-1}}{e^{-t}-1}dt$
Now the Benoulli numbers are defined by $\frac{1}{e^t-1}=\sum_{m=0}^{\infty}B_m\frac{t^{m-1}}{m!}$ where $B_0=1, B_1=1/2, B_{2m+1}=0$
How can I use these things to get an expression for $\zeta(-n), n=0,1,2,3...$ in terms of $B_n$
I assume you mean $$\zeta(z) = \frac{\Gamma(1-z)}{2 \pi i} \int_{C} \frac{t^{z-1} }{e^{-t}-1} \, \mathrm dt = -\frac{\Gamma(1-z)}{2 \pi i} \int_{C} \frac{t^{z-1}e^{t} }{e^{t}-1} \, \mathrm dt,$$
where $C$ is a contour that starts at $- \infty$ below the branch cut on the negative real axis, goes around the origin (without enclosing any of the points $z= \pm 2 \pi i, \pm 4 \pi i, \ldots$), and then goes to back to $-\infty$ above the branch cut.
The above integral representation of the Riemann zeta function is valid for all complex values of $z$ except positive integers.
If we let $z= - n$, where $n \in \mathbb{N}_{\ge 0}$, we get
$$ \zeta(-n) = -\frac{\Gamma(n+1)}{2 \pi i} \int_{C} \frac{t^{-n-1}e^{t}}{e^{t}-1} \, \mathrm dt.$$
But since $z$ is now an integer, the integral above and below the branch cut cancel each other, and all we're left with is the circle around the origin.
Therefore,
$$ \begin{align} \zeta(-n) &= -\frac{\Gamma(n+ 1)}{2 \pi i } \, 2 \pi i \, \operatorname{Res}_{t=0} \left(\frac{t^{-n-1} e^{t}}{e^{t}-1} \right) \\ &= - n! \, \operatorname{Res}_{t=0} \left( t^{-n-1} \sum_{m=0}^{\infty} \frac{B_{m}(1)}{m!} t^{m-1} \right) \tag{1}\\ & = - n! \, \frac{B_{n+1}(1)}{(n+1)!} \\ &= - \frac{B_{n+1}(1)}{n+1} \\&= (-1)^{n} \frac{B_{n+1}}{n+1}. \end{align}$$
$(1)$ https://en.wikipedia.org/wiki/Bernoulli_polynomials#Generating_functions