We study the Berezin transform in the context of solvable groups AN (acting on homogeneous cones and Siegel domains) and determine its spectral decomposition, using an explicit integral kernel representation for the associated ‘eigen-operators’ in terms of multivariable hypergeometric functions.