In this paper, we define a differential operator on an open unit disk Δ by using the novel Borel distribution (BD) operator and means of convolution. This operator is adopted to introduce new subclasses of p-valent functions through the principle of differential subordination, and we focus on some interesting inclusion relations of these classes. Additionally, some inclusion relations are derived by using the Bernardi integral operator. Moreover, relevant convolution results are established for a class of analytic functions on Δ, and other results of analytic univalent functions are derived in detail. This study provides a new perspective for developing p-univalent functions with BD and offers valuable understanding for further research in complex analysis.