Cell migration through 3D extracellular matrices is critical to the normal development of tissues and organs and in disease processes, yet adequate analytical tools to characterize 3D migration are lacking. Here, we quantified the migration patterns of individual fibrosarcoma cells on 2D substrates and in 3D collagen matrices and found that 3D migration does not follow a random walk. Both 2D and 3D migration features a non-Gaussian, exponential mean cell velocity distribution, which we show is primarily a result of cellto-cell variations. Unlike in the 2D case, 3D cell migration is anisotropic: velocity profiles display different speed and selfcorrelation processes in different directions, rendering the classical persistent random walk (PRW) model of cell migration inadequate. By incorporating cell heterogeneity and local anisotropy to the PRW model, we predict 3D cell motility over a wide range of matrix densities, which identifies density-independent emerging migratory properties. This analysis also reveals the unexpected robust relation between cell speed and persistence of migration over a wide range of matrix densities.theory | 3D motility | cancer R andom walks are ubiquitous in biology (1). In particular, the motility of bacteria and eukaryotic cells in the absence of symmetry-breaking gradients has long been described in terms of random walk statistics. Eukaryotic cell migration is a complex process that is a tightly regulated and critical to the normal development of organs and tissues (2-4). Cell migration is activated in a wide range of human diseases, including cancer metastasis (5, 6), immunological responses (7), and wound healing (8). Most of what we know about eukaryotic cell migration at a mechanistic molecular level has stemmed from well-controlled studies of cell migration on flat dishes (i.e., 2D environment). However, cell migration in vivo often forces cells to remodel, exert pulling forces on, and move through a 3D collagen I-rich matrix. Recent work has demonstrated that mechanisms of 3D migration are often different from their 2D counterparts (9-15). Migration on 2D dishes, which induces a basal-apical polarization of the cell, is driven by actomyosin contractility of stress fibers between large focal adhesions and the formation of a wide lamellipodium terminated by thin filopodial protrusions at the leading cellular edge (4, 16). The same cells in collagen-rich 3D matrix do not display a lamellipodium or filopodia. Instead, they display highly dendritic pseudopodial protrusions controlled by distinct proteins that rely both on acto-myosin contractility and microtubule assembly/disassembly dynamics (11, 17). 3D cell migration depends on the expression of metalloproteinases (MMPs), which are dispensable in 2D migration, and physical properties of the 3D matrix (18), such as pore size (6). Recent work has also shown how cancer cells in 3D can alternate between a mesenchymal and an amoeboid migratory phenotype depending on the physical properties of the matrix (19,20) and MMP inhibition (17...