Extensive synteny conservation of holocentric chromosomes in Lepidoptera despite high rates of local genome rearrangements

The recent assembly of the silkworm Bombyx mori genome with 432 Mb on 28 holocentric chromosomes has become a reference in the genomic analysis of the very diverse Order of Lepidoptera. We sequenced BACs from two major pests, the noctuid moths Helicoverpa armigera and Spodoptera frugiperda, corresponding to 15 regions distributed on 11 B. mori chromosomes, each BAC/region being anchored by known orthologous gene(s) to analyze syntenic relationships and genome rearrangements among the three species. Nearly 300 genes and numerous transposable elements were identified, with long interspersed nuclear elements and terminal inverted repeats the most abundant transposable element classes. There was a high degree of synteny conservation between B. mori and the two noctuid species. Conserved syntenic blocks of identified genes were very small, however, approximately 1.3 genes per block between B. mori and the two noctuid species and 2.0 genes per block between S. frugiperda and H. armigera. This corresponds to approximately two chromosome breaks per Mb DNA per My. This is a much higher evolution rate than among species of the Drosophila genus and may be related to the holocentric nature of the lepidopteran genomes. We report a large cluster of eight members of the aminopeptidase N gene family that we estimate to have been present since the Jurassic. In contrast, several clusters of cytochrome P450 genes showed multiple lineage-specific duplication events, in particular in the lepidopteran CYP9A subfamily. Our study highlights the value of the silkworm genome as a reference in lepidopteran comparative genomics.