Fickian yet non-Gaussian diffusion has been observed in several biological and soft matter systems, but the underlying reasons behind the emergence of non-Gaussianity while simultaneously retaining the linear nature of the mean square displacement remain speculative. Here, we perform a set of controlled experiments that quantitatively explore the effect of spatial heterogeneities on the appearance of non-Gaussianity in Fickian diffusion. We study the diffusion of fluorescent colloidal particles in a matrix of micropillars having a range of structural configurations: from completely ordered to completely random. Structural randomness and density are found to be the two most important factors in making diffusion non-Gaussian. We show that non-Gaussianity emerges as a direct consequence of two coupled factors. First, individual particle diffusivities become spatially dependent in a heterogeneous environment. Second, the spatial distribution of the particles varies significantly in heterogeneous environments, which further influences the diffusivity of a single particle. As a result, we find that considerable non-Gaussianity appears even for weak disorder in the arrangement of the micropillars. A simple simulation validates our hypothesis that non-Gaussian yet Fickian diffusion in our system arises from the superstatistical behavior of the ensemble in a structurally heterogeneous environment. The two mechanisms identified here are relevant for many systems of crowded heterogeneous environments where non-Gaussian diffusion is frequently observed, for example in biological systems, polymers, gels and porous materials. PACS numbers: 82.70.Dd, 05.40.-a, 47.56.+rEinsteins theory of Brownian motion shows that for colloidal particles diffusing in two-dimensions in a homogenous, Newtonian fluid, the mean square displacement (MSD) is given by MSD = 4Dτ where D is the diffusion coefficient and τ is the lag time. The solution of the diffusion equation is obtained as a Gaussian probability distribution of displacements (G(∆x)) as is expected for random, independent displacements. However, in many systems including granular materials [1], turbulent flow [2], active gels [3], glassy materials [4], porous materials [5,6], nanoparticles diffusing in polymer melt [7], log-return of stock prices [8] and biological systems [9-13], a non-Gaussian G(∆x) has been observed. In many of these systems the non-Gaussian nature of G(∆x) is related to a process of anomalous diffusion. In such cases the MSD itself is non-linear in time and given by MSD = 4Dτ n , n being the diffusion exponent with a value < 1 for subdiffusion and > 1 for superdiffusion. However, in several very surprising cases, diffusion has been observed to be Fickian, that is the MSD remains linear in time, but counter-intuitively G(∆x) is non-Gaussian. This peculiar behavior has been referred to as 'Anomalous, yet Brownian' or 'Fickian yet non-Gaussian' diffusion (FNG), and has been observed in a wide variety of systems ranging from tracer colloids diffusing in suspensions of ...