Let Sym denote the algebra of symmetric functions and P µ ( · ; q, t) and Q µ ( · ; q, t) be the Macdonald symmetric functions (recall that they differ by scalar factors only). The (q, t)-Cauchy identityexpresses the fact that the P µ ( · ; q, t)'s form an orthogonal basis in Sym with respect to a special scalar product · , · q,t . The present paper deals with the inhomogeneous interpolation Macdonald symmetric functionsThese functions come from the N -variate interpolation Macdonald polynomials, extensively studied in the 90's by Knop, Okounkov, and Sahi. The goal of the paper is to construct symmetric functions H µ ( · ; q, t) with the biorthogonality property I µ ( · ; q, t), H ν ( · ; q, t) q,t = δ µν .These new functions live in a natural completion Sym ⊃ Sym. As a corollary one obtains a new Cauchy-type identity in which the interpolation Macdonald polynomials are paired with certain multivariate rational symmetric functions. The degeneration of this identity in the Jack limit (q, t) = (q, q k ) → (1, 1) is also described.