The talk reviews recent results on statistics of wave functions in disordered systems.  The main emphasis is put on multifractality of critical wave functions at localization transitions.  The analytical and numerical studies of multifractality for various types of critical points are presented: Anderson transition for different spatial dimensionalities, power-law random banded matrix model, Quantum Hall transitions in normal and superconducting systems.