In this study, the connection between certain subfamilies of harmonic univalent functions is established by utilizing a convolution operator involving the Mittag–Leffler function. The investigation reveals inclusion relations concerning harmonic γ-uniformly starlike mappings in the open unit disc, harmonic starlike functions and harmonic convex functions, highlighting the improvements given by the results presented here on previously published outcomes.