Studies of the large $N$ behaviour of the topological properties of gauge theories typically focused on the large $N$ scaling of the topological susceptibility. A much more difficult task is the study of the behaviour of higher cumulants of the topological charge in the large $N$ limit, which up to now remained elusive. We will present first results confirming the expected large $N$ scaling of the coefficient commonly denoted by $b_2$, related to the kurtosis of the topological charge.