En cierto sentido, la respuesta es sí, existe una relación, pasando por Teorema de Riemann-Roch (dos corolarios de la misma, en realidad). De hecho, se puede demostrar Gauss-Bonnet (al menos para las superficies de Riemann $^1$ ) con Cauchy y RR.
Pruebas: En una superficie riemanniana siempre podemos encontrar coordenadas locales, en las que el tensor métrico es conforme y la curvatura gaussiana toma la forma particularmente fácil
$$\mathbf{g}=\frac12\lambda^2(z,\overline{z}) (dz\otimes d\overline{z}+d\overline{z}\otimes dz)=\lambda^2(dx\otimes dx+dy\otimes dy)\\ K=-\frac1{\lambda^2}\Delta\log(\lambda)$$
Por partición de la unidad, podemos unirlos a una construcción global. La forma de volumen $dS$ adopta la forma $dS=\lambda^2 dx\wedge dy$ por lo que obtenemos
$$KdS=-\Delta \log(\lambda)dx\wedge dy$$
Utilizando Operadores de Dolbeaut podemos escribirlo como $$KdS=2id(\partial \log(\lambda))$$ Ahora, dada una diferencial meromorfa $1-$ formulario $\omega$ (que existe gracias a RR), se puede demostrar que, escribiendo localmente $\omega=fdz$ la expresión $$\varphi=\frac{\lambda}{|f|}$$ define una función meromorfa. Dado que en particular $\ln(|f|)$ es armónico, podemos escribir $$KdS=2id(\partial\log(\varphi))$$
Pasemos ahora al teorema de Gauss-Bonnet: escribiendo $X_\varepsilon$ como la superficie de riemann menos pequeños "discos" $D_{k,\varepsilon}$ alrededor de las singularidades $z_k$ de $\varphi$ tenemos
$$\int_X KdS=\lim_{\varepsilon\to 0}\int_{X_\varepsilon}2id(\partial \log(\varphi))=\lim_{\varepsilon\to 0} 2i\int_{\partial X_\varepsilon}\partial \log(\varphi)=\lim_{\varepsilon\to 0}\sum _k\int_{\partial D_{k,\varepsilon}}\partial \log(\varphi) $$
Ahora, cerca de las singularidades de $\varphi$ es decir, los polos o ceros de $f$ podemos escribir $\varphi$ como $\frac{\psi}{|z|^m}$ . Así, por la fórmula integral de Cauchy $$\lim_{\varepsilon\to 0}\int_{\partial D_{k,\varepsilon}}\partial \log(\varphi)=\lim_{\varepsilon\to 0}-m_k\int_{\partial D_{k,\varepsilon}}\partial \log(|z|)=\lim_{\varepsilon\to 0}-\frac {m_k}{2}\int_{\partial D_{k,\varepsilon}}\frac{dz}{z}=-\pi im_k$$
Dado que, por RR, el grado del divisor canónico es $2g-2$ y $\sum m_k=-\text{deg}(\omega)$ Así, obtenemos
$$\int_X KdS=-2\pi \text{deg}(\omega)=2\pi (2-2g)=2\pi \chi(X)$$
$^1:$ En realidad se puede demostrar, aunque es mucho más difícil, que toda superficie orientable admite un sistema de coordenadas local de este tipo, llamado coordenadas isotérmicas . Las pruebas habituales se basan en la solución de la Ecuación de Beltrami