In ΛCDM with simple single-field inflation, the primordial density field is expected to be extremely close to Gaussian, so non-Gaussian signatures in large-scale structure (LSS) should arise mainly from late-time gravitational evolution and galaxy bias, and be modest and well modeled (Desjacques & Seljak 2010; Verde 2010). However, LSS surveys reveal a mix of higher-order clustering statistics, bispectrum shapes, and scale-dependent biases that are sometimes difficult to reconcile cleanly with a purely Gaussian, single-parameter fNL description plus standard nonlinear growth, leaving room for either more complex early-universe physics or unaccounted-for structure formation effects beyond the minimal ΛCDM framework (Matarrese & Verde 2008; Assassi et al. 2015).