By combining analytical and numerical calculations, we investigate the minimal-energy shape of short DNA loops of approximately 100 base pairs (bp). We show that in these loops the excess twist density oscillates as a response to an imposed bending stress, as recently found in DNA minicircles and observed in nucleosomal DNA. These twist oscillations, here referred to as twist waves, are due to the coupling between twist and bending deformations, which in turn originates from the asymmetry between DNA major and minor grooves. We introduce a simple analytical variational shape that reproduces the exact loop energy up to the fourth significant digit and is in very good agreement with shapes obtained from coarse-grained simulations. We, finally...