Two coupled spatial birth-and-death Markov evolutions on $\mathbb{R}^d$ are obtained as unique weak solutions to the associated Fokker-Planck equations. Such solutions are constructed by its associated sequence of correlation functions satisfying the so-called Ruelle-bound. Using the general scheme of Vlasov scaling we are able to derive a system of non-linear, non-local mesoscopic equations describing the effective density of the particle system. The results are applied to several models of ecology and biology.