the corrections to neutrino sector go as $M_R/M_{\rm Pl}$ since $M_R \propto u$. In fact the usual Weinberg operator $(L L H H)/M_{\rm Pl}$ can give a correction like $\frac{M_R}{M_{\rm Pl}} M_{\nu}$ where $M_{\nu}$ is the leading order type I seesaw contribution to light neutrino mass. For example, if the scale of $A_4$ breaking at renormalisable level is $u \sim 10^{16}$ GeV, similar to typical grand unified theory (GUT) scale, one can have a slightly lower scale $M_R \sim 10^{14}$ GeV and hence the scale $a, b$ in the right handed neutrino mass matrix given in , by appropriate tuning of relevant Yukawa couplings. Such $M_R$ can generate sub-eV scale light neutrino mass matrix at leading order for order one Dirac Yukawa couplings. If we consider the corrections to Dirac neutrino mass matrix, such a scale would correspond to $y_2, y_3 \leq 10^{-3}-10^{-2}$ which will correspond to the lowermost regions (along the y-axis) of the parameter space shown in figure \[fig:nhy1y2\] and figure \[fig:ihy1y2\] described before. Similar estimates about these couplings can be made for different symmetry breaking scale. However, in order to have sizeable corrections to the TBM mixing and to avoid unnaturally fine-tuned as well as non-perturbative Yukawa 