For considered model existence and uniqueness of the global positive solution is proven. Then, some other properties, such as boundedness, moment and pathwise estimation, extinction, non-persistence in time average and weak persistence, are established. Considered features are natural requirements from the biological point of view. Finally, examples with numerical simulations are given to illustrate our results.