A series study of $Y(4260)$ are performed with the world's largest data sample around $\sqrt{s} =$ 4.26 GeV at BESIII, similar structures around 4.22 GeV, denoted as $Y(4220)$, are observed in the cross section distributions for $e^+ e^- \to \pi^+ \pi^- J/\psi$, $\pi^+\pi^- h_c$, $\pi^+ \pi^- \psi(3686)$ and $\pi^+ D^0 D^{*-}$. The parameters are consist with each other within errors.
Enhancements around similar positions are observed for $e^+ e^- \to K^+ K^- J/\psi$, $\eta h_c$ and $\gamma \eta_c$. Several light hadrons decay channels are investigated to search for $Y(4220)$ and no clear signal is observed.