The properties of diffusion processes are drastically affected by heterogeneities of the medium that can induce non-Gaussian behavior of the propagator in contrast with the idealized realm of Brownian motion. In this paper we analyze the diffusion propagator when distinct origins of heterogeneity (e.g. time-fractional diffusion, diffusing diffusivity, distributed diffusivity across a population) are combined. These combinations allow one to describe new classes of strongly heterogeneous processes relevant to biology. Based on a combined subordination technique, we obtain the exact propagator for different instructive examples. This approach is then used to calculate analytically the first-passage time statistics (on half-real line and in any bounded domain) for a particle undergoing non-Gaussian diffusion of mixed origins.