We find the asymptotic representation of the solution of the variable-order fractional diffusion equation, which remains unsolved since it was proposed in [Checkin et. al., J. Phys. A, 2005]. We identify a new advection term that causes ultra-slow spatial aggregation of subdiffusive particles due to dominance over the standard advection and diffusion terms, in the long-time limit. This uncovers the anomalous mechanism by which non-uniform distributions can occur. We perform Monte Carlo simulations of the underlying anomalous random walk and find good agreement with the asymptotic solution.Anomalous diffusion has attracted immense interest in the past due to many physical, chemical and biological processes characterized by the mean square displacement (MSD) involving the fractional exponent µ:Anomalous diffusion is observed also in many other areas, for instance, in finance and economics [7]. An influential paper by Metzler and Klafter [2] reviews anomalous diffusion in the scope of a constant exponent µ. However, anomalous transport in realistic inhomogeneous and complex environments [8], such as lipid granules [9], porous media [10] and entangled polymer liquids [11], requires a multi-fractional approach involving the space-dependent variable-order fractional exponent [12][13][14][15][16][17][18]. Important examples of anomalous transport involving multi-fractional exponents are lateral diffusion of proteins on crowded lipid membranes [19], intracellular subdiffusion of proteins [20], mRNA [21] and organelles [22] due in part to inhomogeneous crowding [23] and weak interactions between components in the cell [24]. Recent observations show that lysosomes, which are key organelles for cellular metabolism, predominantly move subdiffusively and maintain a non-uniform spatial distribution in the cell [24]. The majority of these organelles are concentrated in the perinuclear area. A fundamental unresolved question is how lysosomes are self-organized spatially to coordinate their roles [24]. In this Letter, we propose a new anomalous mechanism by which nonuniform distribution of subdiffusing organelles can occur.A generic model for anomalous diffusion in inhomogeneous media is the space-dependent variable-order fractional diffusion equation [12][13][14][15][16] ∂p(x, t) ∂twhere p(x, t) is the probability density function (PDF) of a particle at position x and time t. This function can be also interpreted as the mean number density of subdiffusive particles. In Eq. (1), D µ(x) = a 2 /2τis the fractional diffusion coefficient with the microscopic time scale τ 0 , length scale a, and space-dependent fractional exponent µ(x) ∈ (0, 1). The Riemann-Liouville deriva- (2) also involves spatial dependence. Equation (1) was first derived by Chechkin, Gorenflo and Sokolov [12], and since then, many attempts have been made to find a solution through composite regions with constant anomalous exponents and numerically [12,14,25]. However, Eq. (1) remains unsolved for the general case of a spacedependent anomalous exponent µ(x).In thi...