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 $\mu$m emission. We made single-dish observations toward 38 of 53 regions with 44 GHz masers detected to search for 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.

We present a survey of 56 massive star-forming regions in the 44 GHz methanol maser transition made with the Karl G. Jansky Very Large Array (VLA); 24 of the 56 fields showed maser 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 micron emission. We find a total of 83 maser components with linewidths ranging from 0.17 to 3.3 km s$^{-1}$ with a nearly flat distribution and a median value of 1.1 km s$^{-1}$. The relative velocities of the masers with respect to the systemic velocity of the host clouds range from $-$2.5 to 3.1 km s$^{-1}$ with a distribution peaking near zero. We also study the correlation between the masers and the so-called extended green objects (EGOs) from the GLIMPSE survey. Multiple sources in each field are revealed from IR images as well as from centimeter continuum emission from VLA archival data; in the majority of cases the 44 GHz masers are positionally correlated with EGOs which seem to trace the younger sources in the fields. We report a possible instance of a 44 GHz maser associated with a low-mass protostar. If confirmed, this region will be the fifth known star-forming region that hosts Class I masers associated with low-mass protostars. We discuss three plausible cases of maser variability.

Sensitive and high angular resolution ($\sim$ $0\rlap.{"}7$) (sub)millimeter line and continuum observations of the massive star forming region W3(OH) made with the Submillimeter Array are presented. We report the first detection of two bipolar outflows emanating from the young and massive "Turner-Welch" [TW] protobinary system detected by the emission of the carbon monoxide. The outflows are massive (10 M$_\odot$), highly-collimated (10$^\circ$), and seem to be the extended molecular component of the strong radio jets and a 22 GHz maser water outflow energized also by the stars in the W3(OH)TW system. Observations of the 890 $\mu$m continuum emission and the thermal emission of the CH$_3$OH might suggest the presence of two rotating circumstellar disk-like structures associated with the binary system. The disks-like structures have sizes of about 1500 AU, masses of a few M$_\odot$ and appear to energize the molecular outflows and radio jets. We estimate that the young stars feeding the outflows and that are surrounded by the massive disk-like structures maybe are B-type.