We present a general theory to calculate the steady-state heat and electronic currents for nonlinear systems using a perturbative expansion in the system-bath coupling. We explicitly demonstrate that using the truncated Dyson-series leads to divergences in the steady-state limit, thus making it impossible to be used for actual applications. In order to resolve the divergences we propose a unique choice of initial-condition for the reduced density matrix, which removes the divergences at each order. Our approach not only allows us to use the truncated Dyson-series, with a reasonable choice of initial condition, but also gives the expected result that the steady-state solutions should be independent of initial preparations. Using our improved Dyson series we evaluate the heat and electronic currents upto fourth-order in system-bath coupling, a considerable improvement over the standard quantum master equation techniques. We then numerically corroborate our theory for archetypal settings of linear systems using the exact nonequilibrium Green's function approach. Lastly, to demonstrate the advantage of our approach we deal with the nonlinear spin-boson model to evaluate heat current upto fourth-order and find signatures of cotunnelling process. ; Comment: 14 pages, 5 figures. This version is accepted in the Journal of Chemical Physics