We consider random orthonormal polynomials
F
n
(
x
)
=
∑
i
=
0
n
ξ
i
p
i
(
x
)
,
\begin{equation*} F_{n}(x)=\sum _{i=0}^{n}\xi _{i}p_{i}(x), \end{equation*}
where
ξ
0
\xi _{0}
, …,
ξ
n
\xi _{n}
are independent random variables with zero mean, unit variance and uniformly bounded
(
2
+
ε
)
(2+\varepsilon )
moments, and
(
p
n
)
n
=
0
∞
(p_n)_{n=0}^{\infty }
is the system of orthonormal polynomials with respect to a fixed compactly supported measure on the real line.
Under mild technical assumptions satisfied by many classes of classical polynomial systems, we establish universality for the leading asymptotics of the average number of real roots of
F
n
F_n
, both globally and locally.
Prior to this paper, these results were known only for random orthonormal polynomials with Gaussian coefficients (see D. D. Lubinsky, I. E. Pritsker, and X. Xie [Proc. Amer. Math. Soc. 144 (2016), pp. 1631–1642]) using the Kac-Rice formula, a method that does not extend to the generality of our paper.