We present a generic approach for treating the effect of nuclear motion in high-order harmonic generation from polyatomic molecules. Our procedure relies on a separation of nuclear and electron dynamics where we account for the electronic part using the Lewenstein model and nuclear motion enters as a nuclear correlation function. We express the nuclear correlation function in terms of Franck-Condon factors, which allows us to decompose nuclear motion into modes and identify the modes that are dominant in the high-order harmonic generation process. We show results for the isotopes CH 4 and CD 4 and thereby provide direct theoretical support for a recent experiment [S. Baker et al., Science 312, 424 (2006)] that uses high-order harmonic generation to probe the ultrafast structural nuclear rearrangement of ionized methane.