This article is a brief introduction to a modern polarization theory and first-principles calculations of maximally localized Wannier functions. Practical calculations for these quantities are presented as well as a useful application to the downfolding problem with the ab initio Wannier functions.