It has been a longstanding question how stars are packed or distributed in the universe. Here we show that the spatial distribution of stars follows some geometric rules although there may be local inhomogeneities. The topological and metric properties of their Voronoi tessellations exhibit a lognormal distribution. Moreover, these distributions can be connected with those established for the packing of spheres, corresponding to the loosest packing state. The findings provide a new angle to understand the packing structure of stars in the universe. (C) 2020 Elsevier B.V. All rights reserved.