In quite general N-component ferromagnetic spin systems, it is proved that an arbitrary correlation function is bounded by the corresponding correlation function of a Gaussian model. The bound is useful for the analysis of high-temperature behavior of the system. Similar bounds for truncated correlation functions are also obtained for a class of single-component spin systems.