I present some recent progress in the calculation of quarkonium inclusive and electromagnetic decay widths in the framework of the NRQCD/pNRQCD factorization formalism. I will express, under general dynamical assumptions, the 40 NRQCD matrix elements, which describe up to order mv5 all the ratios of hadronic and electromagnetic inclusive decay widths for bottomonium and charmonium states below threshold, in terms of only 6 non-perturbative universal parameters. New predictions for charmonium and bottomonium decay widths, following from this drastic reduction, will be presented.