This paper presents a novel spectral algorithm for the numerical solution of multi-dimensional fractional-order telegraph equations, a critical model used to capture the combined effects of diffusion and wave propagation. The core innovation of this work is the application of Jacobi-Romanovski polynomials as the basis functions for spectral discretization. These polynomials offer unique advantages, including the ability to handle nonstandard domains and boundary conditions, making them particularly suitable for partial differential equation (PDE) applications. A comprehensive error analysis is conducted, providing deep insights into the convergence rates and factors affecting the accuracy of the numerical solutions. Extensive numerical experiments further demonstrate the superior performance of the proposed spectral algorithm in solving a wide range of multi-dimensional fractional-order telegraph equation models. The results show a significant improvement in accuracy and computational efficiency compared to traditional numerical methods, such as finite difference or finite element techniques. This research advances the field of computational science by offering a robust, efficient, and versatile numerical framework for the precise solution of complex multi-dimensional PDEs.