We prove the existence of infinitely many nonnegative solutions to the following nonlocal elliptic partial differential equation involving singularitieswhere Ω ⊂ R N , N ≥ 2 is a smooth, bounded domain, λ > 0, s ∈ (0, 1), γ(x) ∈ (0, 1) for all x ∈ Ω, N > sp(x, y) for all (x, y) ∈ Ω × Ω and (−∆) s p(•) is the fractional p(•)-Laplacian operator with variable exponent. The nonlinear function f satisfies certain growth conditions. Moreover, we establish a uniform L ∞ ( Ω) estimate of the solution(s) by the Moser iteration technique.