Divergence along the benthic-pelagic axis is one of the most widespread and repeated patterns of morphological variation in fishes, producing body shape diversity associated with ecology and swimming mechanics. This ecological shift is also the first stage of the explosive adaptive radiation of cichlid fishes in the East African Rift Lakes. We use two hybrid crosses of cichlids (Metriaclima sp. x Aulonocara sp. and Labidochromis sp. x Labeotropheus sp., >975 animals total) along the benthic-pelagic ecomorphological axis to determine the genetic basis of body shape diversification. Using a series of both linear and geometric shape measurements, we identify 55 quantitative trait loci (QTL) that underlie various aspects of body shape variation associated with benthic-pelagic divergence. These QTL are spread throughout the genome, each explain 3.0-7.2% of phenotypic variation, and are largely modular. Further, QTL are distinct both between these two crosses of Lake Malawi cichlids and compared to previously identified QTL for body shape in fishes such as sticklebacks. We find that body shape is controlled by many genes of small effects. In all, we find that convergent benthic and pelagic body phenotypes commonly observed across fish clades are most likely due to distinct genetic and molecular mechanisms.