I wanted to know if there is a theorem that in writing a Lagrangian if one missed out a term which preserves the (Lie?) symmetry of the other terms and is also marginal then that will necessarily be generated by RG flow to the IR.
(..like a massless scalar field Lagrangian in 3+1 will generate the mass term as a counter term..though there are cases like Gross-Neveu (GN) where where the discrete symmetry is broken dynamically and not in perturbative RG and so there are 2 in equivalent QFTs corresponding to the GN Lagrangian but the relationship between the two non-perturbative vacua re-realizes the discrete symmetry..)
Can one make a statement along these lines in Ward identities? Like what can be in general be said about the fate in the IR of the Ward identities corresponding to the UV symmetries?
Also if nothing else prohibits is RG flow towards the IR also guaranteed to produce all those terms which are higher dimensional and preserve the symmetries of the UV theory?
Will any such terms necessarily come in the form $O/\Lambda^n$ where $\Lambda$ is the UV cut-off and $O$ a higher dimension operator and "n" is such that the whole term has the marginal dimension. (...and then I can interpret that these effective terms vanish in the UV when I take the cut-off to infinity..)
And how much of the answers to the above remain if one is doing non-relativistic theory?
For starters, I don't think this has anything to do with relativity. The same principles should hold in 3d non-reltivistic systems and 4d lorentz invariant theories.
I don't know if my answer will satisfy you and I'm not fully sure about the details, but I'll think out loud.
Your picture of the terms being of the form $O/\Lambda^n$ seems right to me.
Marginal couplings might be a little bit more subtle, since classically they should not be renormalized. But quantum (loop) corrections will make the coupling run logarithmically with the scale. This is called dimensional transmutation.
I would like to mention a related example of custodial symmetry. In condensed matter systems, the UV physics is typically described by a theory on a lattice. Then, the UV theory only has a discrete subgroup of the rotation group as a symmetry. But in the IR, we get a nice effective field theory with full rotational symmetry! To understand this, let's look at the picture in momentum space. For the UV theory, the (pseudo)momentum can take values within, say a cube (assuming our lattice to be cubic). With RG flow to the IR, all operators of dimension $p^4$ or higher are irrelevant and become negligible. So, only the operator with two derivatives ($p^2$) survives in the IR, giving us full rotational symmetry. However, if we didn't have the cubic lattice in the UV, we might have had terms breaking rotational symmetry, which become relevant in the IR. So, a (discrete) symmetry in the UV protects an ehanced symmetry in the IR. This phenomenon is aptly called custodial symmetry. The enhanced IR symmetry may not be exact, but terms breaking that symmetry are suppressed.
Note that the above example is in line with the principles raised in your question and doesn't say anything different. I just found it to be a very cool example, so I hope you indulge me.