In this article, using methods from geometric analysis and theory of heat kernels, we derive qualitative estimates of automorphic cusp forms defined over quaternion algebras. Using which, we prove an average version of the holomorphic QUE conjecture. We then derive quantitative estimates of classical Hilbert modular cusp forms. This is a generalization of the results from [3] and [9] to higher dimensions.The Bergman kernel, which is the reproducing kernel for L 2 -holomorphic functions defined on a domain in C n has been extensively studied in complex analysis. The generalization of the Bergman kernel to complex manifolds as the reproducing kernel for the space of global holomorphic sections of a vector bundle carries the information on the algebraic and geometric structures of the underlying manifolds.Automorphic forms defined over quaternion algebras are global sections of a holomorphic vector bundle. In this article, we derive asymptotic estimates of the Bergman kernel associated to cuspidal holomorphic vector-valued modular forms of large weight defined over quaternion algebras. We prove an average version of the holomorphic QUE conjecture, when the associated complex analytic space is compact. We also derive quantitative estimates of the Bergman kernel associated to classical Hilbert modular cusp forms.In the first half of this article, utilizing certain results from geometric analysis, we derive asymptotic estimates of the Bergman kernel associated to automorphic cusp forms defined over quaternion algebras. Using which, we prove an average version of the holomorphic QUE conjecture.Jorgenson and Kramer have derived optimal estimates of heat kernels defined over hyperbolic Riemann surfaces of finite volume, and used the estimates to derive optimal estimates of the Bergman kernel associated to cusp forms. In the second half of the article, we extend the heat kernel estimates from [9] to the setting of Hilbert modular varieties, to derive quantitative estimates of the Bergman kernel associated to classical Hilbert modular cusp forms. Our estimates for the Bergman kernel are optimal, and complement the asymptotic estimates derived in the first half of the article.Some of the results of this article are known to experts, and the others are very much expected. So this article serves also as a comprehensive collection of estimates for the Bergman kernel associated to automorphic forms defined over quaternion algebras.