A topological shape optimization method is developed using a dual evolution scheme. Design sensitivity at fixed grids is derived using isogeometric boundary integral equation. Transformation of implicit geometry from level sets into parametric NURBS curves. NURBS geometry is used for shape variations and the level sets for topological ones. Minimization of energy functional to resolve mismatch between level sets and NURBS.