$\scriptstyle \color{grey}{\text{(eliding the convergence problems)}}$
The Eiseinstein series $$G_{2k}(\tau) = \sum_{(n,m) \in \mathbb{Z}^2 \setminus (0,0)} \frac{1}{(m+n\tau)^{2k}} \tag{1}$$ Has the Fourier series representation $$G_{2k}(\tau) = 2 \zeta(2k) + 2 \frac{(2i\pi)^{2k}}{(2k-1)!} \sum_{n=1}^\infty \sigma_{2k-1}(n) e^{2i \pi n \tau}$$
Using $(2\pi n)^{-s} \Gamma(s) = \int_0^\infty x^{s-1} e^{-2 \pi nx}dx$
we have that the Mellin transform of $G_{2k}(ix)-2\zeta(2k)$ is
$$F(s) = \int_0^\infty x^{s-1} (G_{2k}(ix)-2\zeta(k))dx = 2\frac{(2i\pi)^{2k}}{(2k-1)!} \int_0^\infty x^{s-1} \sum_{n=1}^\infty \sigma_{2k-1}(n) e^{-2 \pi n x}dx$$
$$ = 2\frac{(2i\pi)^{2k}}{(2k-1)!}\sum_{n=1}^\infty \sigma_{2k-1}(n)\int_0^\infty x^{s-1}e^{-2 \pi n x}dx = 2\frac{(2i\pi)^{2k}}{(2k-1)!} \Gamma(s) (2\pi)^{-s}\sum_{n =1}^\infty \sigma_{2k-1}(n) n^{-s}$$
$$ = \boxed{\frac{2(-1)^{k}}{(2k-1)!} (2i\pi)^{2k-s}\Gamma(s) \zeta(s) \zeta(s-2k+1)}$$
Since $\zeta(s) = \sum_{n=1}^\infty n^{-s}$ and $\sigma_{2k-1}(n) = \sum_{d | n} d^{2k-1}$ we have $$\zeta(s) \zeta(s-2k+1) = (\sum_{l=1}^\infty l^{-s})(\sum_{m=1}^\infty m^{2k-1} m^{-s})= \sum_{n=1}^\infty n^{-s} \sum_{d | n} d^{2k-1} = \sum_{n=1}^\infty n^{-s} \sigma_{2k-1}(n)$$
Using $i \pi - 2 i \pi \sum_{d=0}^\infty e^{2 i \pi d \tau} = i\pi - 2 i \pi \frac{1}{1-e^{2i \pi \tau}} = i \pi \frac{e^{2i \pi \tau}+1}{e^{2i \pi \tau}-1}=\pi \cot \pi \tau = \frac{1}{\tau} + \sum_{m=1}^\infty \frac{1}{\tau-m}+\frac{1}{\tau-m}$ and differentiating $2k-1$ times we have $$\frac{(- 2 i \pi)^{2k}}{(2k-1)!} \sum_{d=1}^\infty d^{2k-1} e^{2 i \pi d \tau} = \sum_{m=-\infty}^\infty \frac{1}{(m+\tau)^{2k}}$$
So that
$$G_{2k}(\tau) = \sum_{m \in \mathbb{Z}^*} \frac{1}{m^{2k}} + \sum_{n=1}^\infty \sum_{m=-\infty}^\infty \frac{1}{(m+n\tau)^{2k}}$$
$$ = 2\zeta(2k) + \sum_{n=1}^\infty \frac{(- 2 i \pi)^{2k}}{(2k-1)!} \sum_{d=1}^\infty d^{2k-1} e^{2 i \pi d n \tau}$$
$$ = 2\zeta(2k)+\frac{(- 2 i \pi)^{2k}}{(2k-1)!} \sum_{n=1}^\infty \sigma_{2k-1}(n) e^{2 i \pi n \tau}$$
Of course, the functional equation for $\zeta(s)$ proves that $G_{2k}(\tau)$ is a modular form, but $(1)$ proves directly that $$G_{2k}(-1/\tau) = \tau^{2k} G_{2k}(\tau)$$
i.e. $G_{2k}(i/x)-2\zeta(2k) = (ix)^{2k} G_{2k}(ix) -2\zeta(2k)$
and
$$F(s) = \int_0^\infty x^{s-1} (G_{2k}(ix)-2\zeta(2k))dx \underset{y = 1/x}{=} \int_0^\infty y^{-s-1} (G_{2k}(i/y)-2\zeta(2k))dy $$
$$ = \int_0^\infty y^{-s-1} ((iy)^{2k} G_{2k}(iy) -2\zeta(2k))dy = (-1)^{k} F(2k-s)$$