We characterize the continuous $q$-ultraspherical polynomials in terms of the special form of the coefficients in the expansion $\mathcal{D}_q P_n(x)$ in the basis $\{P_n(x)\}$, $\mathcal{D}_q$ being the Askey--Wilson divided difference operator. The polynomials are assumed to be symmetric, and the connection coefficients are multiples of the reciprocal of the square of the $L^2$ norm of the polynomials. A similar characterization is given for the discrete $q$-ultraspherical polynomials. A new proof of the evaluation of the connection coefficients for big $q$-Jacobi polynomials is given.

We define the higher order Riesz transforms and the Littlewood--Paley $g$-function associated to the differential operator $L_\l f(\theta)=-f''(\theta)-2\l\cot\theta f'(\theta)+\l^2f(\theta)$. We prove that these operators are Calder\'{o}n--Zygmund operators in the homogeneous type space $((0,\pi),(\sin t)^{2\l}\,dt)$. Consequently, $L^p$ weighted, $H^1-L^1$ and $L^\infty-BMO$ inequalities are obtained.

We develop a method for deriving integral representations of certain orthogonal polynomials as moments. These moment representations are applied to find linear and multilinear generating functions for $q$-orthogonal polynomials. As a byproduct we establish new transformation formulas for combinations of basic hypergeometric functions, including a new representation of the $q$-exponential function $\mathcal{E}_q$.

We show that the Meixner, Pollaczek, Meixner-Pollaczek, the continuous $q$-ultraspherical polynomials and Al-Salam-Chihara polynomials, in certain normalization, are moments of probability measures. We use this fact to derive bilinear and multilinear generating functions for some of these polynomials. We also comment on the corresponding formulas for the Charlier, Hermite and Laguerre polynomials.