Consider a time homogeneous {0, 1}-valued m-dependent Markov chain \(\{X_{- m + 1 + n}, n \geqslant 0\}\). In this paper, we study the joint probability distribution of number of 0-runs of length \(k_{0} (k_{0} \geqslant m)\) and number of 1-runs of length \(k_{1} (k_{1} \geqslant m)\) in n trials. We study the joint distributions based on five popular counting schemes of runs. The main tool used to obtain the probability generating function of the joint distribution is the conditional probability generating function method. Further a compact method for the evaluation of exact joint distribution is developed. For higher-order two-state Markov chain, these joint distributions are new in the literature of distributions of run statistics. We use these distributions to derive some waiting time distributions.