The Topological Sensitivity Analysis (TSA) is represented by a scalar function, called Topological Derivative (TD), that gives for each point of the domain the sensitivity of a given cost function when an infinitesimal hole is created. Applications to the Laplace, Poisson, Helmoltz, Navier, Stokes and Navier-Stokes equations can be found in the literature. In the present work, an approximated TD expression applied to nonlinear hyperelasticity using the two parameter Mooney-Rivlin constitutive model is obtained by a numerical asymptotic analysis. The cost function is the total potential energy functional. The weak form of state equation is the constraint and the total Lagrangian formulation used. Numerical results of the presented approach are considered for hyperelastic plane problems.