Understanding catalysts strain dynamic behaviours is crucial for the development of cost-effective, efficient, stable and long-lasting catalysts. Here, we reveal in situ three-dimensional strain evolution of single gold nanocrystals during a catalytic CO oxidation reaction under operando conditions with coherent X-ray diffractive imaging. We report direct observation of anisotropic strain dynamics at the nanoscale, where identically crystallographically-oriented facets are qualitatively differently affected by strain leading to preferential active sites formation. Interestingly, the single nanoparticle elastic energy landscape, which we map with attojoule precision, depends on heating versus cooling cycles. The hysteresis observed at the single particle level is following the normal/inverse hysteresis loops of the catalytic performances. This approach opens a powerful avenue for studying, at the single particle level, catalytic nanomaterials and deactivation processes under operando conditions that will enable profound insights into nanoscale catalytic mechanisms.