We study the limiting spectral distribution of sample covariance matrices XX T , where X are p × n random matrices with correlated entries and p/n → y ∈ [0, ∞). If y > 0, we obtain the Marčenko-Pastur distribution and in the case y = 0 the semicircle distribution after appropriate rescaling. The entries we consider are Curie-Weiss spins, which are correlated random signs, where the degree of the correlation is governed by an inverse temperature β > 0. The model exhibits a phase transition at β = 1. The correlation between any two entries is of order O((np) −1 ) for β ∈ (0, 1), O((np) −1/2 ) for β = 1, and for β > 1 the correlation does not vanish in the limit. In our proofs we use Stieltjes transforms and concentration of random quadratic forms.