We present observations of massive star-forming regions selected from the IRAS Point Source Catalog. The observations were made with the Very Large Array and the Large Millimeter Telescope to search for Class I methanol masers. We made interferometric observations of 125 massive star-forming regions in the 44 GHz methanol maser transition; 53 of the 125 fields showed emission. The data allow us to demonstrate associations, at arcsecond precision, of the Class I maser emission with outflows, HII regions and shocks traced by 4.5 μm emission. We made single-dish observations toward 38 of the 53 regions with 44 GHz masers detected to search for the methanol transitions at 84.5, 95.1, 96.7, 107.0, and 108.8 GHz. We find detection rates of 74, 55, 100, 3, and 45%, respectively. We used a wide-band receiver which revealed many other spectral lines that are common in star-forming regions.