  We initiate the study of spectral zeta functions $\zeta_{X}$ for finite and
infinite graphs $X$, instead of the Ihara zeta function, with a perspective
towards zeta functions from number theory and connections to hypergeometric
functions. The Riemann hypothesis is shown to be equivalent to an approximate
functional equation of graph zeta functions. The latter holds at all points
where Riemann's zeta function $\zeta(s)$ is non-zero. This connection arises
via a detailed study of the asymptotics of the spectral zeta functions of
finite torus graphs in the critcal strip and estimates on the real part of the
logarithmic derivative of $\zeta(s)$. We relate $\zeta_{\mathbb{Z}}$ to Euler's
beta integral and show how to complete it giving the functional equation
$\xi_{\mathbb{Z}}(1-s)=\xi_{\mathbb{Z}}(s)$. This function appears in the
theory of Eisenstein series although presumably with this spectral
intepretation unrecognized. In higher dimensions $d$ we provide a meromorphic
continuation of $\zeta_{\mathbb{Z}^{d}}(s)$ to the whole plane and identify the
poles. From our aymptotics several known special values of $\zeta(s)$ are
derived as well as its non-vanishing on the line $Re(s)=1$. We determine the
spectral zeta functions of regular trees and show it to be equal to a
specialization of Appell's hypergeometric function $F_{1}$ via an Euler-type
integral formula due to Picard.
