First draft genome assembly of the desert locust, Schistocerca gregaria

Chromosome-level genome assembly of the morabine grasshopper Vandiemenella viatica19

Morabine grasshoppers in the Vandiemenella viatica species group, which show karyotype diversity, have been studied for their ecological distribution and speciation in relation to their genetic and chromosomal diversity. They are good models for studying sex chromosome evolution as "old" and newly emerged sex chromosomes co-exist within the group. Here we present a reference genome for the viatica19 chromosomal race, that possesses the ancestral karyotype within the group. Using PacBio HiFi and Hi-C sequencing, we generated a chromosome-level assembly of 4.09 Gb in span, scaffold N50 of 429 Mb, and complete BUSCO score of 98.1%, containing 10 pseudo-chromosomes. We provide Illumina datasets of males and females, used to identify the X chromosome. The assembly contains 19,034 predicted protein-coding genes, and a total of 75.21% of repetitive DNA sequences. By leveraging HiFi reads, we mapped the genome-wide distribution of methylated bases (5mC and 6 mA). This comprehensive assembly offers a robust reference for morabine grasshoppers and supports further research into speciation and sex chromosome diversification within the group and its related species.

The grasshopper genome reveals long-term gene content conservation of the X Chromosome and temporal variation in X Chromosome evolution

We present the first chromosome-level genome assembly of the grasshopper, Locusta migratoria, one of the largest insect genomes. We use coverage differences between females (XX) and males (X0) to identify the X Chromosome gene content, and find that the X Chromosome shows both complete dosage compensation in somatic tissues and an underrepresentation of testis-expressed genes. X-linked gene content from L. migratoria is highly conserved across seven insect orders, namely Orthoptera, Odonata, Phasmatodea, Hemiptera, Neuroptera, Coleoptera, and Diptera, and the 800 Mb grasshopper X Chromosome is homologous to the fly ancestral X Chromosome despite 400 million years of divergence, suggesting either repeated origin of sex chromosomes with highly similar gene content, or long-term conservation of the X Chromosome. We use this broad conservation of the X Chromosome to test for temporal dynamics to Fast-X evolution, and find evidence of a recent burst evolution for new X-linked genes in contrast to slow evolution of X-conserved genes.

Genome size and chromosome number are critical metrics for accurate genome assembly assessment in Eukaryota

The number of genome assemblies has rapidly increased in recent history, with NCBI databases reaching over 41,000 eukaryotic genome assemblies across about 2,300 species. Increases in read length and improvements in assembly algorithms have led to increased contiguity and larger genome assemblies. While this number of assemblies is impressive, only about a third of these assemblies have corresponding genome size estimations for their respective species on publicly available databases. In this paper, genome assemblies are assessed regarding their total size compared to their respective publicly available genome size estimations. These deviations in size are assessed related to genome size, kingdom, sequencing platform, and standard assembly metrics, such as N50 and BUSCO values. A large proportion of assemblies deviate from their estimated genome size by more than 10%, with increasing deviations in size with increased genome size, suggesting nonprotein coding and structural DNA may be to blame. Modest differences in performance of sequencing platforms are noted as well. While standard metrics of genome assessment are more likely to indicate an assembly approaching the estimated genome size, much of the variation in this deviation in size is not explained with these raw metrics. A new, proportional N50 metric is proposed, in which N50 values are made relative to the average chromosome size of each species. This new metric has a stronger relationship with complete genome assemblies and, due to its proportional nature, allows for a more direct comparison across assemblies for genomes with variation in sizes and architectures.

Chromosomal-Level Reference Genome for the Chinese Endemic Pygmy Grasshopper, Zhengitettix transpicula, Sheds Light on Tetrigidae Evolution and Advancing Conservation Efforts

The pygmy grasshopper, Zhengitettix transpicula, is a Chinese endemic species with an exceedingly limited distribution and fragile population structure, rendering it vulnerable to extinction. We present a high-continuity, chromosome-scale reference genome assembly to elucidate this species' distinctive biology and inform conservation. Employing an integrated sequencing approach, we achieved a 970.40 Mb assembly with 96.32% coverage across seven pseudo-chromosomes and impressive continuity (N50 > 220 Mb). Genome annotation achieves identification with 99.2% BUSCO completeness, supporting quality. Comparative analyses with 14 genomes from Orthoptera-facilitated phylogenomics and revealed 549 significantly expanded gene families in Z. transpicula associated with metabolism, stress response, and development. However, genomic analysis exposed remarkably low heterozygosity (0.02%), implying a severe genetic bottleneck from small, fragmented populations, characteristic of species vulnerable to extinction from environmental disruptions. Elucidating the genetic basis of population dynamics and specialization provides an imperative guideline for habitat conservation and restoration of this rare organism. Moreover, divergent evolution analysis of the CYP305m2 gene regulating locust aggregation highlighted potential structural and hence functional variations between Acrididae and Tetrigidae. Our chromosomal genomic characterization of Z. transpicula advances Orthopteran resources, establishing a framework for evolutionary developmental explorations and applied conservation genomics, reversing the trajectory of this unique grasshopper lineage towards oblivion.

Neuromodulation and the toolkit for behavioural evolution: can ecdysis shed light on an old problem?

The geneticist Thomas Dobzhansky famously declared: 'Nothing in biology makes sense except in the light of evolution'. A key evolutionary adaptation of Metazoa is directed movement, which has been elaborated into a spectacularly varied number of behaviours in animal clades. The mechanisms by which animal behaviours have evolved, however, remain unresolved. This is due, in part, to the indirect control of behaviour by the genome, which provides the components for both building and operating the brain circuits that generate behaviour. These brain circuits are adapted to respond flexibly to environmental contingencies and physiological needs and can change as a function of experience. The resulting plasticity of behavioural expression makes it difficult to characterize homologous elements of behaviour and to track their evolution. Here, we evaluate progress in identifying the genetic substrates of behavioural evolution and suggest that examining adaptive changes in neuromodulatory signalling may be a particularly productive focus for future studies. We propose that the behavioural sequences used by ecdysozoans to moult are an attractive model for studying the role of neuromodulation in behavioural evolution.

Sizing Up the Onychophoran Genome: Repeats, Introns, and Gene Family Expansion Contribute to Genome Gigantism in Epiperipatus broadwayi

Genome assemblies are growing at an exponential rate and have proved indispensable for studying evolution but the effort has been biased toward vertebrates and arthropods with a particular focus on insects. Onychophora or velvet worms are an ancient group of cryptic, soil dwelling worms noted for their unique mode of prey capture, biogeographic patterns, and diversity of reproductive strategies. They constitute a poorly understood phylum of exclusively terrestrial animals that is sister group to arthropods. Due to this phylogenetic position, they are crucial in understanding the origin of the largest phylum of animals. Despite their significance, there is a paucity of genomic resources for the phylum with only one highly fragmented and incomplete genome publicly available. Initial attempts at sequencing an onychophoran genome proved difficult due to its large genome size and high repeat content. However, leveraging recent advances in long-read sequencing technology, we present here the first annotated draft genome for the phylum. With a total size of 5.6Gb, the gigantism of the Epiperipatus broadwayi genome arises from having high repeat content, intron size inflation, and extensive gene family expansion. Additionally, we report a previously unknown diversity of onychophoran hemocyanins that suggests the diversification of copper-mediated oxygen carriers occurred independently in Onychophora after its split from Arthropoda, parallel to the independent diversification of hemocyanins in each of the main arthropod lineages.

Neuroendocrinal and molecular basis of flight performance in locusts

Insect flight is a complex physiological process that involves sensory and neuroendocrinal control, efficient energy metabolism, rhythmic muscle contraction, and coordinated wing movement. As a classical study model for insect flight, locusts have attracted much attention from physiologists, behaviorists, and neuroendocrinologists over the past decades. In earlier research, scientists made extensive efforts to explore the hormone regulation of metabolism related to locust flight; however, this work was hindered by the absence of molecular and genetic tools. Recently, the rapid development of molecular and genetic tools as well as multi-omics has greatly advanced our understanding of the metabolic, molecular, and neuroendocrinal basis of long-term flight in locusts. Novel neural and molecular factors modulating locust flight and their regulatory mechanisms have been explored. Moreover, the molecular mechanisms underlying phase-dependent differences in locust flight have also been revealed. Here, we provide a systematic review of locust flight physiology, with emphasis on recent advances in the neuroendocrinal, genetic, and molecular basis. Future research directions and potential challenges are also addressed.

Studying phenotypic variation and DNA methylation across development, ecology and evolution in the clonal marbled crayfish: a paradigm for investigating epigenotype-phenotype relationships in macro-invertebrates

Animals can produce different phenotypes from the same genome during development, environmental adaptation and evolution, which is mediated by epigenetic mechanisms including DNA methylation. The obligatory parthenogenetic marbled crayfish, Procambarus virginalis, whose genome and methylome are fully established, proved very suitable to study this issue in detail. Comparison between developmental stages and DNA methylation revealed low expression of Dnmt methylation and Tet demethylation enzymes from the spawned oocyte to the 256 cell embryo and considerably increased expression thereafter. The global 5-methylcytosine level was 2.78% at mid-embryonic development and decreased slightly to 2.41% in 2-year-old adults. Genetically identical clutch-mates raised in the same uniform laboratory setting showed broad variation in morphological, behavioural and life history traits and differences in DNA methylation. The invasion of diverse habitats in tropical to cold-temperate biomes in the last 20 years by the marbled crayfish was associated with the expression of significantly different phenotypic traits and DNA methylation patterns, despite extremely low genetic variation on the whole genome scale, suggesting the establishment of epigenetic ecotypes. The evolution of marbled crayfish from its parent species Procambarus fallax by autotriploidy a few decades ago was accompanied by a significant increase in body size, fertility and life span, a 20% reduction of global DNA methylation and alteration of methylation in hundreds of genes, suggesting that epigenetic mechanisms were involved in speciation and fitness enhancement. The combined analysis of phenotypic traits and DNA methylation across multiple biological contexts in the laboratory and field in marbled crayfish may serve as a blueprint for uncovering the role of epigenetic mechanisms in shaping of phenotypes in macro-invertebrates.

Prothoracicostatic Activity of the Ecdysis-Regulating Neuropeptide Crustacean Cardioactive Peptide (CCAP) in the Desert Locust

Accurate control of innate behaviors associated with developmental transitions requires functional integration of hormonal and neural signals. Insect molting is regulated by a set of neuropeptides, which trigger periodic pulses in ecdysteroid hormone titers and coordinate shedding of the old cuticle during ecdysis. In the current study, we demonstrate that crustacean cardioactive peptide (CCAP), a structurally conserved neuropeptide described to induce the ecdysis motor program, also exhibits a previously unknown prothoracicostatic activity to regulate ecdysteroid production in the desert locust, Schistocerca gregaria. We identified the locust genes encoding the CCAP precursor and three G protein-coupled receptors that are activated by CCAP with EC₅₀ values in the (sub)nanomolar range. Spatiotemporal expression profiles of the receptors revealed expression in the prothoracic glands, the endocrine organs where ecdysteroidogenesis occurs. RNAi-mediated knockdown of CCAP precursor or receptors resulted in significantly elevated transcript levels of several Halloween genes, which encode ecdysteroid biosynthesis enzymes, and in elevated ecdysteroid levels one day prior to ecdysis. Moreover, prothoracic gland explants exhibited decreased secretion of ecdysteroids in the presence of CCAP. Our results unequivocally identify CCAP as the first prothoracicostatic peptide discovered in a hemimetabolan species and reveal the existence of an intricate interplay between CCAP signaling and ecdysteroidogenesis.

Comparative Analysis of Transposable Elements in Genus Calliptamus Grasshoppers Revealed That Satellite DNA Contributes to Genome Size Variation

Simple Summary

Calliptamus is a genus of grasshoppers belonging to the family Acrididae. The genus Calliptamus includes approximately 17 recognized species. Calliptamus abbreviatus, Calliptamus italicus, and Calliptamus barbarus are three species that are widely found in northern China. These species are polyphagous, feeding on a variety of wild plants as well as crops, particularly legumes. The genome sizes, phylogenetic position, and transcriptome analysis of the genus Calliptamus were already known previous to this research. The repeatome analysis of these species was missing, which is directly linked to the larger genome sizes of the grasshoppers. Here, we classified repetitive DNA sequences at the level of superfamilies and sub-families, and found that LINE, TcMar-Tc1 and Ty3-gypsy LTR retrotransposons dominated the repeatomes of all genomes, accounting for 16–34% of the total genomes of these species. Satellite DNA dynamic evolutionary changes in all three genomes played a role in genome size evolution. This study would be a valuable source for future genome assemblies.

Abstract

Transposable elements (TEs) play a significant role in both eukaryotes and prokaryotes genome size evolution, structural changes, duplication, and functional variabilities. However, the large number of different repetitive DNA has hindered the process of assembling reference genomes, and the genus level TEs diversification of the grasshopper massive genomes is still under investigation. The genus Calliptamus diverged from Peripolus around 17 mya and its species divergence dated back about 8.5 mya, but their genome size shows rather large differences. Here, we used low-coverage Illumina unassembled short reads to investigate the effects of evolutionary dynamics of satDNAs and TEs on genome size variations. The Repeatexplorer2 analysis with 0.5X data resulted in 52%, 56%, and 55% as repetitive elements in the genomes of Calliptamus barbarus, Calliptamus italicus, and Calliptamus abbreviatus, respectively. The LINE and Ty3-gypsy LTR retrotransposons and TcMar-Tc1 dominated the repeatomes of all genomes, accounting for 16–35% of the total genomes of these species. Comparative analysis unveiled that most of the transposable elements (TEs) except satDNAs were highly conserved across three genomes in the genus Calliptamus grasshoppers. Out of a total of 20 satDNA families, 17 satDNA families were commonly shared with minor variations in abundance and divergence between three genomes, and 3 were Calliptamus barbarus specific. Our findings suggest that there is a significant amplification or contraction of satDNAs at genus phylogeny which is the main cause that made genome size different.

Genome Size Estimation and Full-Length Transcriptome of Sphingonotus tsinlingensis: Genetic Background of a Drought-Adapted Grasshopper

Sphingonotus Fieber, 1852 (Orthoptera: Acrididae), is a grasshopper genus comprising approximately 170 species, all of which prefer dry environments such as deserts, steppes, and stony benchlands. In this study, we aimed to examine the adaptation of grasshopper species to arid environments. The genome size of Sphingonotus tsinlingensis was estimated using flow cytometry, and the first high-quality full-length transcriptome of this species was produced. The genome size of S. tsinlingensis is approximately 12.8 Gb. Based on 146.98 Gb of PacBio sequencing data, 221.47 Mb full-length transcripts were assembled. Among these, 88,693 non-redundant isoforms were identified with an N50 value of 2,726 bp, which was markedly longer than previous grasshopper transcriptome assemblies. In total, 48,502 protein-coding sequences were identified, and 37,569 were annotated using public gene function databases. Moreover, 36,488 simple tandem repeats, 12,765 long non-coding RNAs, and 414 transcription factors were identified. According to gene functions, 61 cytochrome P450 (CYP450) and 66 heat shock protein (HSP) genes, which may be associated with drought adaptation of S. tsinlingensis, were identified. We compared the transcriptomes of S. tsinlingensis and two other grasshopper species which were less tolerant to drought, namely Mongolotettix japonicus and Gomphocerus licenti. We observed the expression of CYP450 and HSP genes in S. tsinlingensis were higher. We produced the first full-length transcriptome of a Sphingonotus species that has an ultra-large genome. The assembly characteristics were better than those of all known grasshopper transcriptomes. This full-length transcriptome may thus be used to understand the genetic background and evolution of grasshoppers.

The Evolutionary Patterns of Genome Size in Ensifera (Insecta: Orthoptera)

Genomic size variation has long been a focus for biologists. However, due to the lack of genome size data, the mechanisms behind this variation and the biological significance of insect genome size are rarely studied systematically. The detailed taxonomy and phylogeny of the Ensifera, as well as the extensive documentation concerning their morphological, ecological, behavioral, and distributional characteristics, make them a strong model for studying the important scientific problem of genome size variation. However, data on the genome size of Ensifera are rather sparse. In our study, we used flow cytometry to determine the genome size of 32 species of Ensifera, the smallest one being only 1C = 0.952 pg with the largest species up to 1C = 19.135 pg, representing a 20-fold range. This provides a broader blueprint for the genome size variation of Orthoptera than was previously available. We also completed the assembly of nine mitochondrial genomes and combined mitochondrial genome data from public databases to construct phylogenetic trees containing 32 species of Ensifera and three outgroups. Based on these inferred phylogenetic trees, we detected the phylogenetic signal of genome size variation in Ensifera and found that it was strong in both males and females. Phylogenetic comparative analyses revealed that there were no correlations between genome size and body size or flight ability in Tettigoniidae. Reconstruction of ancestral genome size revealed that the genome size of Ensifera evolved in a complex pattern, in which the genome size of the grylloid clade tended to decrease while that of the non-grylloid clade expanded significantly albeit with fluctuations. However, the evolutionary mechanisms underlying variation of genome size in Ensifera are still unknown.

Cholinergic calcium responses in cultured antennal lobe neurons of the migratory locust

Locusts are advantageous organisms to elucidate mechanisms of olfactory coding at the systems level. Sensory input is provided by the olfactory receptor neurons of the antenna, which send their axons into the antennal lobe. So far, cellular properties of neurons isolated from the circuitry of the olfactory system, such as transmitter-induced calcium responses, have not been studied. Biochemical and immunocytochemical investigations have provided evidence for acetylcholine as classical transmitter of olfactory receptor neurons. Here, we characterize cell cultured projection and local interneurons of the antennal lobe by cytosolic calcium imaging to cholinergic stimulation. We bulk loaded the indicator dye Cal-520 AM in dissociated culture and recorded calcium transients after applying cholinergic agonists and antagonists. The majority of projection and local neurons respond with increases in calcium levels to activation of both nicotinic and muscarinic receptors. In local interneurons, we reveal interactions lasting over minutes between intracellular signaling pathways, mediated by muscarinic and nicotinic receptor stimulation. The present investigation is pioneer in showing that Cal-520 AM readily loads Locusta migratoria neurons, making it a valuable tool for future research in locust neurophysiology, neuropharmacology, and neurodevelopment.

Gene transcription changes in a locust model of noise-induced deafness

Locusts have auditory structures called Müller’s organs attached to tympanic membranes on either side of the abdomen. We measured the normalized abundances of 500 different mRNA transcripts in 320 Müller’s organs obtained from 160 locusts (Schistocerca gregaria) that had been subjected to a loud continuous 3-kHz tone for 24 h. Abundance ratios were then measured relative to transcripts from 360 control organs. A histogram of the number of observed transcripts versus their abundance ratios (noise exposed/control) was well fitted by a Cauchy distribution with median value near one. Transcripts below 5% and above 95% of the cumulative distribution function of the fitted Cauchy distribution were selected as putatively different from the expected values of an untreated preparation. This yielded eight transcripts with ratios increased by noise exposure (ratios 1.689–3.038) and 18 transcripts with reduced ratios (0.069–0.457). Most of the transcripts with increased abundance represented genes responsible for cuticular construction, suggesting extensive remodeling of some or all the cuticular components of the auditory structure, whereas the reduced abundance transcripts were mostly involved in lipid and protein storage and metabolism, suggesting a profound reduction in metabolic activity in response to the overstimulation.

NEW & NOTEWORTHY Locust ears have functional and genetic similarities to human ears, including loss of hearing from age or noise exposure. We measured transcript abundances in transcriptomes of noise-exposed and control locust ears. The data indicate remodeling of the ear tympanum and profound reductions in metabolism that may explain reduced sound transduction. These findings advance our understanding of this useful model and suggest further experiments to elucidate mechanisms that ears use to cope with excessive stimulation.

Crucial Role of Juvenile Hormone Receptor Components Methoprene-Tolerant and Taiman in Sexual Maturation of Adult Male Desert Locusts

Currently (2020), Africa and Asia are experiencing the worst desert locust (Schistocerca gregaria) plague in decades. Exceptionally high rainfall in different regions caused favorable environmental conditions for very successful reproduction and population growth. To better understand the molecular mechanisms responsible for this remarkable reproductive capacity, as well as to fill existing knowledge gaps regarding the regulation of male reproductive physiology, we investigated the role of methoprene-tolerant (Scg-Met) and Taiman (Scg-Tai), responsible for transducing the juvenile hormone (JH) signal, in adult male locusts. We demonstrated that knockdown of these components by RNA interference strongly inhibits male sexual maturation, severely disrupting reproduction. This was evidenced by the inability to show mating behavior, the absence of a yellow-colored cuticle, the reduction of relative testes weight, and the drastically reduced phenylacetonitrile (PAN) pheromone levels of the treated males. We also observed a reduced relative weight, as well as relative protein content, of the male accessory glands in Scg-Met knockdown locusts. Interestingly, in these animals the size of the corpora allata (CA), the endocrine glands where JH is synthesized, was significantly increased, as well as the transcript level of JH acid methyltransferase (JHAMT), a rate-limiting enzyme in the JH biosynthesis pathway. Moreover, other endocrine pathways appeared to be affected by the knockdown, as evidenced by changes in the expression levels of the insulin-related peptide and two neuroparsins in the fat body. Our results demonstrate that JH signaling pathway components play a crucial role in male reproductive physiology, illustrating their potential as molecular targets for pest control.

Identification of Gonadulin and Insulin-Like Growth Factor From Migratory Locusts and Their Importance in Reproduction in Locusta migratoria

Many insect species have several genes coding for insulin-related peptides (IRPs), but so far only a single IRP gene has been identified in migratory locusts. Here, we report and characterize two other genes coding for peptides that are related to insulin, namely gonadulin and arthropod insulin-like growth factor (aIGF); peptides postulated to be orthologs of Drosophila melanogaster insulin-like peptides 8 and 6 respectively. In Locusta migratoria the aIGF transcript is expressed in multiple tissues as was previously reported for IRP in both L. migratoria and Schistocerca gregaria, but there are significant differences in expression patterns between the two species. The gonadulin transcript, however, seems specific to the ovary, whereas its putative receptor transcript is expressed most abundantly in the ovary, fat body and the central nervous system. Since the central nervous system-fat body-ovary axis is essential for successful reproduction, we studied the influence of gonadulin on vitellogenesis and oocyte growth. A reduction in the gonadulin transcript (via RNA interference) led to a significant reduction in vitellogenin mRNA levels in the fat body and a strong oocyte growth inhibition, thus suggesting an important role for gonadulin in reproduction in this species.

Comparative analysis of morabine grasshopper genomes reveals highly abundant transposable elements and rapidly proliferating satellite DNA repeats

BACKGROUND

Repetitive DNA sequences, including transposable elements (TEs) and tandemly repeated satellite DNA (satDNAs), collectively called the "repeatome", are found in high proportion in organisms across the Tree of Life. Grasshoppers have large genomes, averaging 9 Gb, that contain a high proportion of repetitive DNA, which has hampered progress in assembling reference genomes. Here we combined linked-read genomics with transcriptomics to assemble, characterize, and compare the structure of repetitive DNA sequences in four chromosomal races of the morabine grasshopper Vandiemenella viatica species complex and determine their contribution to genome evolution.

RESULTS

We obtained linked-read genome assemblies of 2.73-3.27 Gb from estimated genome sizes of 4.26-5.07 Gb DNA per haploid genome of the four chromosomal races of V. viatica. These constitute the third largest insect genomes assembled so far. Combining complementary annotation tools and manual curation, we found a large diversity of TEs and satDNAs, constituting 66 to 75% per genome assembly. A comparison of sequence divergence within the TE classes revealed massive accumulation of recent TEs in all four races (314-463 Mb per assembly), indicating that their large genome sizes are likely due to similar rates of TE accumulation. Transcriptome sequencing showed more biased TE expression in reproductive tissues than somatic tissues, implying permissive transcription in gametogenesis. Out of 129 satDNA families, 102 satDNA families were shared among the four chromosomal races, which likely represent a diversity of satDNA families in the ancestor of the V. viatica chromosomal races. Notably, 50 of these shared satDNA families underwent differential proliferation since the recent diversification of the V. viatica species complex.

CONCLUSION

This in-depth annotation of the repeatome in morabine grasshoppers provided new insights into the genome evolution of Orthoptera. Our TEs analysis revealed a massive recent accumulation of TEs equivalent to the size of entire Drosophila genomes, which likely explains the large genome sizes in grasshoppers. Despite an overall high similarity of the TE and satDNA diversity between races, the patterns of TE expression and satDNA proliferation suggest rapid evolution of grasshopper genomes on recent timescales.

First draft genome assembly of the desert locust, Schistocerca gregaria

Background: At the time of publication, the most devastating desert locust crisis in decades is affecting East Africa, the Arabian Peninsula and South-West Asia. The situation is extremely alarming in East Africa, where Kenya, Ethiopia and Somalia face an unprecedented threat to food security and livelihoods. Most of the time, however, locusts do not occur in swarms, but live as relatively harmless solitary insects. The phenotypically distinct solitarious and gregarious locust phases differ markedly in many aspects of behaviour, physiology and morphology, making them an excellent model to study how environmental factors shape behaviour and development. A better understanding of the extreme phenotypic plasticity in desert locusts will offer new, more environmentally sustainable ways of fighting devastating swarms. Methods: High molecular weight DNA derived from two adult males was used for Mate Pair and Paired End Illumina sequencing and PacBio sequencing. A reliable reference genome of Schistocerca gregaria was assembled using the ABySS pipeline, scaffolding was improved using LINKS. Results: In total, 1,316 Gb Illumina reads and 112 Gb PacBio reads were produced and assembled. The resulting draft genome consists of 8,817,834,205 bp organised in 955,015 scaffolds with an N50 of 157,705 bp, making the desert locust genome the largest insect genome sequenced and assembled to date. In total, 18,815 protein-encoding genes are predicted in the desert locust genome, of which 13,646 (72.53%) obtained at least one functional assignment based on similarity to known proteins. Conclusions: The desert locust genome data will contribute greatly to studies of phenotypic plasticity, physiology, neurobiology, molecular ecology, evolutionary genetics and comparative genomics, and will promote the desert locust's use as a model system. The data will also facilitate the development of novel, more sustainable strategies for preventing or combating swarms of these infamous insects.

First draft genome assembly of the desert locust, Schistocerca gregaria

Background: At the time of publication, the most devastating desert locust crisis in decades is affecting East Africa, the Arabian Peninsula and South-West Asia. The situation is extremely alarming in East Africa, where Kenya, Ethiopia and Somalia face an unprecedented threat to food security and livelihoods. Most of the time, however, locusts do not occur in swarms, but live as relatively harmless solitary insects. The phenotypically distinct solitarious and gregarious locust phases differ markedly in many aspects of behaviour, physiology and morphology, making them an excellent model to study how environmental factors shape behaviour and development. A better understanding of the extreme phenotypic plasticity in desert locusts will offer new, more environmentally sustainable ways of fighting devastating swarms. Methods: High molecular weight DNA derived from two adult males was used for Mate Pair and Paired End Illumina sequencing and PacBio sequencing. A reliable reference genome of Schistocerca gregaria was assembled using the ABySS pipeline, scaffolding was improved using LINKS. Results: In total, 1,316 Gb Illumina reads and 112 Gb PacBio reads were produced and assembled. The resulting draft genome consists of 8,817,834,205 bp organised in 955,015 scaffolds with an N50 of 157,705 bp, making the desert locust genome the largest insect genome sequenced and assembled to date. In total, 18,815 protein-encoding genes are predicted in the desert locust genome, of which 13,646 (72.53%) obtained at least one functional assignment based on similarity to known proteins. Conclusions: The desert locust genome data will contribute greatly to studies of phenotypic plasticity, physiology, neurobiology, molecular ecology, evolutionary genetics and comparative genomics, and will promote the desert locust's use as a model system. The data will also facilitate the development of novel, more sustainable strategies for preventing or combating swarms of these infamous insects.