Abstract Invariant ensembles, which are characterised by the joint distribution of eigenvalues P ( λ 1 , …, λ N ), play a central role in random matrix theory. We consider the truncated linear statistics L K = ∑ n = 1 K f ( λ n ) with 1 ⩽ K ⩽ N , λ 1 > λ 2 > ⋯ > λ N and f a given function. This quantity has been studied recently in the case where the function f is monotonous. Here, we consider the general case, where this function can be non-monotonous. Motivated by the physics of cold atoms, we study the example f ( λ ) = λ 2 in the Gaussian ensembles of random matrix theory. Using the Coulomb gas method, we obtain the distribution of the truncated linear statistics, in the limit N → ∞ and K → ∞, with κ = K / N fixed. We show that the distribution presents two essential singularities, which arise from two infinite order phase transitions for the underlying Coulomb gas. We further argue that this mechanism is universal, as it depends neither on the choice of the ensemble, nor on the function f .