Multitudes of motile cilia covering epithelial surfaces in the mammalian airways, fallopian tubes, and the brain synchronize their beat to enable efficient fluid transport. The nature of the mechanisms leading to (large scale) ciliary coordination remains unclear. I will review our understanding of this process based on physics-based mathematical models where hydrodynamic forces between near beating cilia are sufficient to lead to cilia coordination. Further, I will demonstrate that it is possible to reach, and transition between,multiple synchronization states by varying the intrinsic activity of the cilium and the strength of hydrodynamic coupling between neighboring cilia. I will conclude by commenting on the implications of these findings to our understanding of cilia coordination in epithelial tissues under healthy and diseased conditions as well as in unicellular organisms where transitions between distinct modes of cilia coordination are crucial for the cell behavior and survival.