General characteristics of methanol (CH3OH) maser emission are summarized. It is shown that methanol maser sources are concentrated in the spiral arms. Most of the methanol maser sources from the Perseus arm are associated with embedded stellar clusters and a considerable portion is situated close to compact HII regions. Almost 1/3 of the Perseus Arm sources lie at the edges of optically identified HII regions which means that massive star formation in the Perseus Arm is to a great extent triggered by local phenomena. A multiline analysis of the methanol masers allows us to determine the physical parameters in the regions of maser formation. Maser modelling shows that class II methanol masers can be pumped by the radiation of the warm dust as well as by free-free emission of a hypercompact region hcHII with a turnover frequency exceeding 100 GHz. Methanol masers of both classes can reside in the vicinity of hcHIIs. Modelling shows that periodic changes of maser fluxes can be reproduced by variations of the dust temperature by a few percent which may be caused by variations in the brightness of the central young stellar object reflecting the character of the accretion process. Sensitive observations have shown that the masers with low flux densities can still have considerable amplification factors. The analysis of class I maser surveys allows us to identify four distinct regimes that differ by the series of their brightest lines.

Methanol masers which are traditionally divided into two classes provide possibility to study important parts of the star forming regions: Class~II masers trace vicinities of the massive YSOs while class~I masers are likely to trace more distant parts of the outflows where newer stars can form. There are many methanol transitions which produce observed masers. This allows to use pumping analysis for estimation of the physical parameters in the maser formation regions and its environment, for the study of their evolution. Extensive surveys in different masing transitions allow to conclude on the values of the temperatures, densities, dust properties, etc. in the bulk of masing regions. Variability of the brightest masers is monitored during several years. In some cases it is probably caused by the changes of the dust temperature which follow variations in the brightness of the central YSO reflecting the character of the accretion process. A unified catalogue of the class II methanol masers consisting of more than 500 objects is compiled. Analysis of the data shows that: physical conditions within the usual maser source vary considerably; maser brightness is determined by parameters of some distinguished part of the object - maser formation region; class II methanol masers are formed not within the outflows but in the regions affected by their propagation. It is shown that the "near" solutions for the kinematic distances to the sources can be used for statistical analysis. The luminosity function of the 6.7 GHz methanol masers is constructed. It is shown that improvement of the sensitivity of surveys can increase number of detected maser sources considerably.