At transonic condition, the design of natural laminar flow (NLF) wing is challenging, as the extension of the laminar flow needs to be finely balanced with the potential wave drag increase. To achieve this balance, it is proposed to unlock the wing sweep and introduce three-dimensional (3-D) contour shock-control bump (SCB) in the optimization of NLF infinite swept wing aiming at total drag reduction. A 3-D Reynolds-Averaged Navier-Stokes flow solver is extended to incorporate laminar-turbulent transition prediction due to streamwise and crossflow instabilities. The flow solver is integrated in a gradient-based optimization framework. The transition criteria including streamwise and crossflow instabilities are coupled in the sensitivity calculation using a discrete adjoint solver. Transonic design optimization of the sectional profile and wing sweep angle is first conducted at a Mach number of 0.78. The optimization managed to alleviate the shock wave, reducing the pressure drag while it failed to extend the laminar flow. Further, a combined optimization of the wing with a parameterized SCB at the same condition is carried out. The optimized wing with SCB features a long favorable pressure gradient region and a low sweep angle which allows for a large proportion of laminar flow without significant pressure drag. The shock wave is controlled by the 3-D bump and the SCB has little effect on the upstream flow.