In this talk, we will review some recent applications of operad theory to singular SPDEs. It is an essential tool for characterising the chain rule symmetry in the full subcritical regime which leads to renormalising quasilinear SPDEs with local counterterms. Also, it provides negative results about the existence of another combinatorial set lying between multi-indices and decorated trees for encoding expansions of solutions of singular SPDEs.