The rapid evolution of repetitive DNA sequences, including heterochromatic regions, satellite DNA, tandem duplications, and transposable elements, can underlie phenotypic evolution and contribute to hybrid incompatibilities between species. However, repetitive genomic regions are fragmented in most contemporary genome assemblies. We generated highly contiguous de novo assemblies for the Drosophila simulans species complex ( D. simulans, D. mauritiana, and D. sechellia ), which speciated~250,000 years ago. These species diverged from their common ancestor with D. melanogaster 3 million years ago. Our assemblies are comparable in contiguity and accuracy to the current D. melanogaster genome, allowing us to directly compare repetitive regions in genomes across different evolutionary times. We find a rapid turnover of satellite DNA and extensive structural variation in heterochromatic regions, while the euchromatic gene content is mostly conserved. Despite the overall preservation of synteny, euchromatin of each species has been sculpted by clade and species-specific inversions, transposable elements (TE), satellite and tRNA tandem arrays, and gene duplications. We also find Y-linked genes rapidly diverging, in terms of copy number and recent duplications from the autosomes. Our assemblies provide a valuable resource for studying genome evolution and its consequences for phenotypic evolution in these genetic model species.The group of four fruit fly species composed of D. melanogaster, D. simulans, D.sechellia and D. mauritiana is collectively known as the D. melanogaster species complex (or mel-complex for short) (Hey and Kliman 1993) and serves as a model system for studying speciation (Tao et al.