Locusts are one of the world’s most destructive agricultural pests and represent a useful model system in entomology. Here we present a draft 6.5 Gb genome sequence of Locusta migratoria, which is the largest animal genome sequenced so far. Our findings indicate that the large genome size of L. migratoria is likely to be because of transposable element proliferation combined with slow rates of loss for these elements. Methylome and transcriptome analyses reveal complex regulatory mechanisms involved in microtubule dynamic-mediated synapse plasticity during phase change. We find significant expansion of gene families associated with energy consumption and detoxification, consistent with long-distance flight capacity and phytophagy. We report hundreds of potential insecticide target genes, including cys-loop ligand-gated ion channels, G-protein-coupled receptors and lethal genes. The L. migratoria genome sequence offers new insights into the biology and sustainable management of this pest species, and will promote its wide use as a model system.

The migratory locust, Locusta migratoria, is the most widely distributed grasshopper species in the world. However, its global genetic structure and phylogeographic relationships have not been investigated. In this study, we explored the worldwide genetic structure and phylogeography of the locust populations based on the sequence information of 65 complete mitochondrial genomes and three mitochondrial genes of 263 individuals from 53 sampling sites. Although this locust can migrate over long distances, our results revealed high genetic differentiation among the geographic populations. The populations can be divided into two different lineages: the Northern lineage, which includes individuals from the temperate regions of the Eurasian continent, and the Southern lineage, which includes individuals from Africa, southern Europe, the Arabian region, India, southern China, South-east Asia and Australia. An analysis of population genetic diversity indicated that the locust species originated from Africa. Ancestral populations likely separated into Northern and Southern lineages 895 000 years ago by vicariance events associated with Pleistocene glaciations. These two lineages evolved in allopatry and occupied their current distributions in the world via distinct southern and northern dispersal routes. Genetic differences, caused by the long-term independent diversification of the two lineages, along with other factors, such as geographic barriers and temperature limitations, may play important roles in maintaining the present phylogeographic patterns. Our phylogeographic evidence challenged the long-held view of multiple subspecies in the locust species and tentatively divided it into two subspecies, L. m. migratoria and L. m. migratorioides.

BackgroundLepidoptera encompasses more than 160,000 described species that have been classified into 45–48 superfamilies. The previously determined Lepidoptera mitochondrial genomes (mitogenomes) are limited to six superfamilies of the lineage Ditrysia. Compared with the ancestral insect gene order, these mitogenomes all contain a tRNA rearrangement. To gain new insights into Lepidoptera mitogenome evolution, we sequenced the mitogenomes of two ghost moths that belong to the non-ditrysian lineage Hepialoidea and conducted a comparative mitogenomic analysis across Lepidoptera.ResultsThe mitogenomes of Thitarodes renzhiensis and T. yunnanensis are 16,173 bp and 15,816 bp long with an A + T content of 81.28 % and 82.34 %, respectively. Both mitogenomes include 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and the A + T-rich region. Different tandem repeats in the A + T-rich region mainly account for the size difference between the two mitogenomes. All the protein-coding genes start with typical mitochondrial initiation codons, except for cox1 (CGA) and nad1 (TTG) in both mitogenomes. The anticodon of trnS(AGN) in T. renzhiensis and T. yunnanensis is UCU instead of the mostly used GCU in other sequenced Lepidoptera mitogenomes. The 1,584-bp sequence from rrnS to nad2 was also determined for an unspecified ghost moth (Thitarodes sp.), which has no repetitive sequence in the A + T-rich region. All three Thitarodes species possess the ancestral gene order with trnI-trnQ-trnM located between the A + T-rich region and nad2, which is different from the gene order trnM-trnI-trnQ in all previously sequenced Lepidoptera species. The formerly identified conserved elements of Lepidoptera mitogenomes (i.e. the motif ‘ATAGA’ and poly-T stretch in the A + T-rich region and the long intergenic spacer upstream of nad2) are absent in the Thitarodes mitogenomes.ConclusionThe mitogenomes of T. renzhiensis and T. yunnanensis exhibit unusual features compared with the previously determined Lepidoptera mitogenomes. Their ancestral gene order indicates that the tRNA rearrangement event(s) likely occurred after Hepialoidea diverged from other lepidopteran lineages. Characterization of the two ghost moth mitogenomes has enriched our knowledge of Lepidoptera mitogenomes and contributed to our understanding of the mechanisms underlying mitogenome evolution, especially gene rearrangements.