By using Malliavin calculus, Bismut type formulas are established for the Lions derivative of P t f (µ) := Ef (X µ t ), where t > 0, f is a bounded measurable function, and X µ t solves a distribution dependent SDE with initial distribution µ. As applications, explicit estimates are derived for the Lions derivative and the total variational distance between distributions of solutions with different initial data. Both degenerate and nondegenerate situations are considered. Due to the lack of the semigroup property and the invalidity of the formula P t f (µ) = P t f (x)µ(dx), essential difficulties are overcome in the study.