In this paper, we make a detailed discussion on the η and η′-meson leading-twist light-cone distribution amplitude ϕ2;η(′)(u,μ) by using QCD sum rules approach under the background field theory. Taking both the non-perturbative condensates up to dimension-six and NLO QCD corrections to the perturbative part, its first three moments ⟨ξ2;η(′)n⟩∣μ0 with n=(2,4,6) at initial scale μ0=1 GeV can be determined. e.g. ⟨ξ2;η2⟩∣μ0=0.231−0.013+0.010, ⟨ξ2;η4⟩∣μ0=0.109−0.007+0.007, and ⟨ξ2;η6⟩∣μ0=0.066−0.006+0.006 for η-meson, ⟨ξ2;η′2⟩∣μ0=0.211−0.017+0.015, ⟨ξ2;η′4⟩∣μ0=0.093−0.009+0.009, and ⟨ξ2;η′6⟩∣μ0=0.054−0.008+0.008 for η′-meson. Next, we calculate Ds→η(′) TFFs f+η(′)(q2) within QCD light-cone sum rules approach up to NLO level. The values at large recoil region are f+η(0)=0.476−0.036+0.040 and f+η′(0)=0.544−0.042+0.046. After extrapolating TFFs to the allowable physical regions within the series expansion, we obtain the branching fractions of the semi-leptonic decay, i.e. Ds+→η(′)ℓ+νℓ, i.e. B(Ds+→η(′)e+νe)=2.346−0.331+0.418(0.792−0.118+0.141)×10−2 and B(Ds+→η(′)μ+νμ)=2.320−0.327+0.413(0.773−0.115+0.138)×10−2 for ℓ=(e,μ) channels respectively. And in addition to that, the mixing angle for η−η′ with φ and ratio for the different decay channels Rη′/ηℓ are given, which show good agreement with the recent BESIII measurements.