A common approach to calculate the solution of a scalar advection-diffusion equation is by a Feynman-Kac representation which averages over stochastic Lagrangian trajectories going backward in time to the initial conditions and boundary data. The trajectories are obtained by solving SDE's with the advecting velocity as drift and a backward Itō term representing the scalar diffusivity.  In this framework, we present an exact formula for scalar dissipation in terms of the variance of the scalar values acquired along each random trajectory. As an important application, we study the connection between anomalous scalar dissipation in turbulent flows for large Reynolds and Péclet numbers and the spontaneous stochasticity of the Lagrangian particle trajectories. The latter property corresponds to the Lagrangian trajectories remaining random in the limit \(Re,Pe \to \infty\), when the backward Itō term formally vanishes but the advecting velocity field becomes non-Lipschitz. For flows on domains without boundaries and for wall-bounded flows with no-flux Neumann conditions for the scalar, we prove that spontaneous stochasticity is necessary and sufficient for anomalous scalar dissipation. The fluctuation-dissipation relation provides a Lagrangian representation of scalar dissipation also in turbulent flows where present experiments suggest that dissipation is tending to zero as \(Re,Pe \to \infty\). We discuss an illustrative example of Rayleigh-Bénard convection with imposed heat-flux at the top and bottom plates.  This talk presents joint work with Gregory Eyink.