In this paper, we consider a class of normalized harmonic functions in the unit disk satisfying a third-order differential inequality and we investigate several properties of this class such as close to convexity, coefficient bounds, growth estimates, sufficient coefficient condition and convolution. Moreover, as an application, we construct harmonic polynomials involve Gaussian hypergeometric function which belong to the considered class. We also provide examples illustrating graphically with the help of Maple.