Abstract Results of the long-term monitoring observations by the Hitachi 32 m radio telescope of the 6.7 GHz Class II methanol masers associated with four high-mass star-forming regions are presented. We detected periodic flux variability in G06.795−0.257, G10.472+0.027, G12.209−0.102, and G13.657−0.599 with the periods of 968, 1624, 1272, and 1266 d, respectively, although the detected period is tentative due to the short monitoring term relative to the estimated period. The facts that the flux variation patterns show the symmetric sine curves and that the luminosities of the central protostar and periods of maser flux variation are consistent with the expected period–luminosity (PL) relation suggest that the mechanisms of maser flux variability of G10.472+0.027 and G12.209−0.102 can be explained by protostellar pulsation instability. From the PL relation, the central stars of these two sources are expected to be very high-mass protostars with a mass of $\sim 40\, M_{\odot }$ and to have a mass accretion rate of $\sim 2 \times 10^{-2}\, M_{\odot }\:$yr−1. On the other hand, G06.795−0.257 and G13.657−0.599 have intermittent variation patterns and have luminosities that are an order of magnitude smaller than those expected from the PL relation, suggesting that the variation mechanisms of these sources originated from a binary system. Since almost all the maser features vary with the same period regardless of the geometry, periodic accretion models may be appropriate mechanisms for flux variability in G06.795−0.257 and G13.657−0.599.