The respiratory tract is protected by mucus, a complex fluid transported along the epithelial surface by the coordinated beating of millions of microscopic cilia, hence the name of mucociliary clearance. Its impairment is a strong marker of severe chronic respiratory diseases. Yet, the relationship between ciliary density and the spatial scale of mucus transport, as well as the mechanisms that drive ciliary-beat orientations during ciliogenesis are much debated. Here, we show on polarized human bronchial epithelia that mucus swirls and circular orientational order of the underlying ciliary beats emerge and grow during ciliogenesis, until a macroscopic mucus transport is achieved for physiological ciliary densities. By establishing that the macroscopic ciliary-beat order is lost and recovered by removing and adding mucus respectively, we demonstrate that cilia/mucus hydrodynamic interactions govern the collective dynamics of ciliary-beat directions. We propose a twodimensional model that predicts a phase diagram of mucus transport in accordance with the experiments. It paves the way to a predictive in-silico modeling of bronchial mucus transport in health and disease.