A bstract In view of both the latest LHCb measurement of $$ {R}_{K^{\left(\ast \right)}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>R</mml:mi> <mml:msup> <mml:mi>K</mml:mi> <mml:mfenced> <mml:mo>∗</mml:mo> </mml:mfenced> </mml:msup> </mml:msub> </mml:math> and the new 2 . 7 σ deviation reported by Belle II on $$ {B}^{+}\to {K}^{+}\nu \overline{\nu} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msup> <mml:mi>B</mml:mi> <mml:mo>+</mml:mo> </mml:msup> <mml:mo>→</mml:mo> <mml:msup> <mml:mi>K</mml:mi> <mml:mo>+</mml:mo> </mml:msup> <mml:mi>ν</mml:mi> <mml:mover> <mml:mi>ν</mml:mi> <mml:mo>¯</mml:mo> </mml:mover> </mml:math> decays, we present a fit to the B meson anomalies for various one and two dimensional hypothesis including complex Wilson coefficients. We show in a model-independent way that the generic non-universal U(1) ′ extensions of the SM, without flavour violation, fail to simultaneously fit those observables and corroborate that they can modify $$ \textrm{BR}\left({B}^{+}\to {K}^{+}\nu \overline{\nu}\right) $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>BR</mml:mi> <mml:mfenced> <mml:mrow> <mml:msup> <mml:mi>B</mml:mi> <mml:mo>+</mml:mo> </mml:msup> <mml:mo>→</mml:mo> <mml:msup> <mml:mi>K</mml:mi> <mml:mo>+</mml:mo> </mml:msup> <mml:mi>ν</mml:mi> <mml:mover> <mml:mi>ν</mml:mi> <mml:mo>¯</mml:mo> </mml:mover> </mml:mrow> </mml:mfenced> </mml:math> up to only a 10%. In view of this deficit, we propose a new way in which those models can accommodate the data at tree level by introducing lepton flavour violating couplings and non-diagonal elements of the charged lepton mixing matrix, with implications in future charged lepton flavour violation searches.