Abstract. We construct a commutative algebra A z , generated by d algebraically independent q-difference operators acting on variables z 1 , z 2 , . . . , z d , which is diagonalized by the multivariable Askey-Wilson polynomials P n (z) considered by Gasper and Rahman (2005). Iterating Sears' 4 φ 3 transformation formula, we show that the polynomials P n (z) possess a certain duality between z and n. Analytic continuation allows us to obtain another commutative algebra A n , generated by d algebraically independent difference operators acting on the discrete variables n 1 , n 2 , . . . , n d , which is also diagonalized by P n (z). This leads to a multivariable q-Askey-scheme of bispectral orthogonal polynomials which parallels the theory of symmetric functions.