The nested Kingman coalescent describes the dynamics of particles (called genes) contained in larger components (called species), where pairs of species coalesce at constant rate and pairs of genes coalesce at constant rate provided they lie within the same species. We prove that starting from rn species, the empirical distribution of species masses (numbers of genes/n) at time t/n converges as n → ∞ to a solution of the deterministic coagulation-transport equation ∂ t d = ∂ x (ψd) + a(t) (d d − d) , where ψ(x) = cx 2 , denotes convolution and a(t) = 1/(t + δ) with δ = 2/r. The most interesting case when δ = 0 corresponds to an infinite initial number of species. This equation describes the evolution of the distribution of species of mass x, where pairs of species can coalesce and each species' mass evolves likeẋ = −ψ(x). We provide two natural probabilistic solutions of the latter IPDE and address in detail the case when δ = 0. The first solution is expressed in terms of a branching particle system where particles carry masses behaving as independent continuous-state branching processes. The second one is the law of the solution to the following McKean-Vlasov equation dx t = −ψ(x t) dt + v t ∆J t where J is an inhomogeneous Poisson process with rate 1/(t + δ) and (v t ; t ≥ 0) is a sequence of independent random variables such that L(v t) = L(x t). We show that there is a unique solution to this equation and we construct this solution with the help of a marked Brownian coalescent point process. When ψ(x) = x γ , we show the existence of a self-similar solution for the PDE which relates when γ = 2 to the speed of coming down from infinity of the nested Kingman coalescent.