How can I insert message bits into the frequency spectrum of an audio signal using Matlab?
The FFT is just an implementation of the DFT. There's, mathematically, no difference.

The DFT is a linear operation. So it makes absolutely zero difference whether you add up two signals in time or frequency domain.