Hadronic spectral functions of $\tau$ decays have been used in the past to provide an alternative determination of the LO Hadronic Vacuum Polarization relevant for the (g-2) of the muon. Following recent developments and results in Lattice QCD+QED calculations, we explore the possibility of studying the isospin breaking corrections of $\tau$ spectral functions for this prediction. We present preliminary results at physical pion mass based on Domain Wall Fermion ensembles generated by the RBC/UKQCD collaboration.