Assembling large genomes: analysis of the stick insect (Clitarchus hookeri) genome reveals a high repeat content and sex-biased genes associated with reproduction

The chromosome-level genome of Chinese praying mantis Tenodera sinensis (Mantodea: Mantidae) reveals its biology as a predator

BACKGROUND

The Chinese praying mantis, Tenodera sinensis (Saussure), is a carnivorous insect that preys on a variety of arthropods and small vertebrates, including pest species. Several studies have been conducted to understand its behavior and physiology. However, there is limited knowledge about the genetic information underlying its genome evolution, digestive demands, and predatory behaviors.

FINDINGS

Here we have assembled the chromosome-level genome of T. sinensis, representing the first sequenced genome of the family Mantidae, with a genome size of 2.54 Gb and scaffold N50 of 174.78 Mb. Our analyses revealed that 98.6% of BUSCO genes are present, resulting in a well-annotated assembly compared to other insect genomes, containing 25,022 genes. The reconstructed phylogenetic analysis showed the expected topology placing the praying mantis in an appropriate position. Analysis of transposon elements suggested the Gypsy/Dirs family, which belongs to long terminal repeat (LTR) transposons, may be a key factor resulting in the larger genome size. The genome shows expansions in several digestion and detoxification associated gene families, including trypsin and glycosyl hydrolase (GH) genes, ATP-binding cassette (ABC) transporter, and carboxylesterase (CarE), reflecting the possible genomic basis of digestive demands. Furthermore, we have found 1 ultraviolet-sensitive opsin and 2 long-wavelength-sensitive (LWS) opsins, emphasizing the core role of LWS opsins in regulating predatory behaviors.

CONCLUSIONS

The high-quality genome assembly of the praying mantis provides a valuable repository for studying the evolutionary patterns of the mantis genomes and the gene expression profiles of insect predators.

Genome assembly and annotation of the European earwig Forficula auricularia (subspecies B)

The European earwig Forficula auricularia is an important model for studies of maternal care, sexual selection, sociality, and host–parasite interactions. However, detailed genetic investigations of this species are hindered by a lack of genomic resources. Here, we present a high-quality hybrid genome assembly for Forficula auricularia using Nanopore long-reads and 10× linked-reads. The final assembly is 1.06 Gb in length with 31.03% GC content. It consists of 919 scaffolds with an N50 of 12.55 Mb. Half of the genome is present in only 20 scaffolds. Benchmarking Universal Single-Copy Orthologs scores are ∼90% from 3 sets of single-copy orthologs (eukaryotic, insect, and arthropod). The total repeat elements in the genome are 64.62%. The MAKER2 pipeline annotated 12,876 protein-coding genes and 21,031 mRNAs. Phylogenetic analysis revealed the assembled genome as that of species B, one of the 2 known genetic subspecies of Forficula auricularia. The genome assembly, annotation, and associated resources will be of high value to a large and diverse group of researchers working on dermapterans.

Characterization of transposable elements within the Bemisia tabaci species complex

Background

Whiteflies are agricultural pests that cause negative impacts globally to crop yields resulting at times in severe economic losses and food insecurity. The Bemisia tabaci whitefly species complex is the most damaging in terms of its broad crop host range and its ability to serve as vector for over 400 plant viruses. Genomes of whiteflies belonging to this species complex have provided valuable genomic data; however, transposable elements (TEs) within these genomes remain unexplored. This study provides the first accurate characterization of TE content within the B. tabaci species complex.

Results

This study identified that an average of 40.61% of the genomes of three whitefly species (MEAM1, MEDQ, and SSA-ECA) consists of TEs. The majority of the TEs identified were DNA transposons (22.85% average) while SINEs (0.14% average) were the least represented. This study also compared the TE content of the three whitefly genomes with three other hemipteran genomes and found significantly more DNA transposons and less LINEs in the whitefly genomes. A total of 63 TE superfamilies were identified to be present across the three whitefly species (39 DNA transposons, six LTR, 16 LINE, and two SINE). The sequences of the identified TEs were clustered which generated 5766 TE clusters. A total of 2707 clusters were identified as uniquely found within the whitefly genomes while none of the generated clusters were from both whitefly and non-whitefly TE sequences.

This study is the first to characterize TEs found within different B. tabaci species and has created a standardized annotation workflow that could be used to analyze future whitefly genomes.

Conclusion

This study is the first to characterize the landscape of TEs within the B. tabaci whitefly species complex. The characterization of these elements within the three whitefly genomes shows that TEs occupy significant portions of B. tabaci genomes, with DNA transposons representing the vast majority. This study also identified TE superfamilies and clusters of TE sequences of potential interest, providing essential information, and a framework for future TE studies within this species complex.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13100-022-00270-6.

Testes morphology and the identification of transcripts of the hormonal pathways of the velvetbean caterpillar Anticarsia gemmatalis Hübner, 1818 (Lepidoptera: Erebidae)

Anticarsia gemmatalis is one of the main defoliating pests of soybeans in Brazil. In the current study, we characterized the histomorphology of the testes and the spermatogenesis process in A. gemmatalis. We also identified transcripts involved in the biosynthesis, metabolism, and signaling of juvenile and ecdysteroid hormones, in order to provide information about potential mechanisms of regulation of hormonal pathways in this species. Our analyses revealed that the A. gemmatalis larvae have a pair of kidney-shaped testicles. These are divided into four testicular follicles, where there are germ cell cysts at different stages of development. In the pupal stage, the testicles are fused, so adults have a single spherical testis, with a variable number of follicles. The A. gemmatalis has centripetal spermatogenesis and exhibits spermatic dimorphism. We identified 31 transcripts that encode proteins involved in juvenile hormone and ecdysteroid pathways, such as mevalonate kinase, CYP14A1, ecdysone receptor, among others. Our results on the morphology of the testes and spermatogenesis process, as well as identification of the genes involved in hormonal pathways in A. gemmatalis, provide important data for understanding the biology of this agricultural pest, which can be used as a basis for further research in other economically important lepidopterans.

A molecular toolkit for superorganisms

Social insects, such as ants, bees, wasps, and termites, draw biologists' attention due to their distinctive lifestyles. As experimental systems, they provide unique opportunities to study organismal differentiation, division of labor, longevity, and the evolution of development. Ants are particularly attractive because several ant species can be propagated in the laboratory. However, the same lifestyle that makes social insects interesting also hampers the use of molecular genetic techniques. Here, we summarize the efforts of the ant research community to surmount these hurdles and obtain novel mechanistic insight into the biology of social insects. We review current approaches and propose novel ones involving genomics, transcriptomics, chromatin and DNA methylation profiling, RNA interference (RNAi), and genome editing in ants and discuss future experimental strategies.

Divergent Gene Expression Following Duplication of Meiotic Genes in the Stick Insect Clitarchus hookeri

Some animal groups, such as stick insects (Phasmatodea), have repeatedly evolved alternative reproductive strategies, including parthenogenesis. Genomic studies have found modification of the genes underlying meiosis exists in some of these animals. Here we examine the evolution of copy number, evolutionary rate, and gene expression in candidate meiotic genes of the New Zealand geographic parthenogenetic stick insect Clitarchus hookeri. We characterized 101 genes from a de novo transcriptome assembly from female and male gonads that have homology with meiotic genes from other arthropods. For each gene we determined copy number, the pattern of gene duplication relative to other arthropod orthologs, and the potential for meiosis-specific expression. There are five genes duplicated in C. hookeri, including one also duplicated in the stick insect Timema cristinae, that are not or are uncommonly duplicated in other arthropods. These included two sister chromatid cohesion associated genes (SA2 and SCC2), a recombination gene (HOP1), an RNA-silencing gene (AGO2) and a cell-cycle regulation gene (WEE1). Interestingly, WEE1 and SA2 are also duplicated in the cyclical parthenogenetic aphid Acyrthosiphon pisum and Daphnia duplex, respectively, indicating possible roles in the evolution of reproductive mode. Three of these genes (SA2, SCC2, and WEE1) have one copy displaying gonad-specific expression. All genes, with the exception of WEE1, have significantly different nonsynonymous/synonymous ratios between the gene duplicates, indicative of a shift in evolutionary constraints following duplication. These results suggest that stick insects may have evolved genes with novel functions in gamete production by gene duplication.

Genome Size Estimation of Callipogon relictus Semenov (Coleoptera: Cerambycidae), an Endangered Species and a Korea Natural Monument

Simple Summary

The longhorned beetle Calipogon relictus has been considered as a class I endangered species since 2012 in Korea. In an attempt towards beetle conservation, we estimated its genome size at 1.8 ± 0.2 Gb, representing one of the largest cerambycid genomes. This study provides useful insight at the genome level and facilitates the development of an effective conservation strategy.

Abstract

We estimated the genome size of a relict longhorn beetle, Callipogon relictus Semenov (Cerambycidae: Prioninae)—the Korean natural monument no. 218 and a Class I endangered species—using a combination of flow cytometry and k-mer analysis. The two independent methods enabled accurate estimation of the genome size in Cerambycidae for the first time. The genome size of C. relictus was 1.8 ± 0.2 Gb, representing one of the largest cerambycid genomes studied to date. An accurate estimation of genome size of a critically endangered longhorned beetle is a major milestone in our understanding and characterization of the C. relictus genome. Ultimately, the findings provide useful insight into insect genomics and genome size evolution, particularly among beetles.

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.

A draft genome sequence of the elusive giant squid, Architeuthis dux

BACKGROUND

The giant squid (Architeuthis dux; Steenstrup, 1857) is an enigmatic giant mollusc with a circumglobal distribution in the deep ocean, except in the high Arctic and Antarctic waters. The elusiveness of the species makes it difficult to study. Thus, having a genome assembled for this deep-sea-dwelling species will allow several pending evolutionary questions to be unlocked.

FINDINGS

We present a draft genome assembly that includes 200 Gb of Illumina reads, 4 Gb of Moleculo synthetic long reads, and 108 Gb of Chicago libraries, with a final size matching the estimated genome size of 2.7 Gb, and a scaffold N50 of 4.8 Mb. We also present an alternative assembly including 27 Gb raw reads generated using the Pacific Biosciences platform. In addition, we sequenced the proteome of the same individual and RNA from 3 different tissue types from 3 other species of squid (Onychoteuthis banksii, Dosidicus gigas, and Sthenoteuthis oualaniensis) to assist genome annotation. We annotated 33,406 protein-coding genes supported by evidence, and the genome completeness estimated by BUSCO reached 92%. Repetitive regions cover 49.17% of the genome.

CONCLUSIONS

This annotated draft genome of A. dux provides a critical resource to investigate the unique traits of this species, including its gigantism and key adaptations to deep-sea environments.

Transcriptome Analysis of Female and Male Conopomorpha sinensis (Lepidoptera: Gracilariidae) Adults With a Focus on Hormone and Reproduction

Conopomorpha sinensis Bradley is the dominant borer pest of litchi and longan in the Asian-pacific area. Reduction or interference of reproduction and mating of adult moths is one of the most used strategies to control C. sinensis. Insect reproduction is a critical biological process closely related to endocrine control. Conopomorpha sinensis genome and transcriptome information is limited, hampering both our understanding of the molecular mechanisms underlying hormone activity and reproduction and the development of control strategies for this borer pest. To explore the sex differences in gene expression profiles influencing these biological processes, de novo transcriptomes were assembled from female and male adult C. sinensis specimens. This analysis yielded 184,422 unigenes with an average length of 903 bp and 405,961 transcripts after sequencing and assembly. About 45.06, 22.41, 19.53, 34.05, 35.82, 36.42, and 19.85% of the unigenes had significant matches in seven public databases. Subsequently, gene ontology (GO) and kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis revealed comprehensive information about the function of each gene and identified enriched categories and pathways that were associated with the 2,890 female-biased genes and 2,964 male-biased genes. In addition, we identified some important unigenes related to hormone activity and reproduction among the sex-differentially expressed genes (DEGs), including unigenes coding for ecdysone-induced protein 78C, juvenile hormone (JH)-regulated gene fatty acyl-CoA reductase, vitellogenin, etc. Our findings provide a more comprehensive portrait of the sex differences involved in the relationship of two important physiological features-hormone activity and reproduction in C. sinensis and members of the family Gracillariidae.

Incomer, a DD36E family of Tc1/mariner transposons newly discovered in animals

BACKGROUND

The Tc1/mariner superfamily might represent the most diverse and widely distributed group of DNA transposons. Several families have been identified; however, exploring the diversity of this superfamily and updating its classification is still ongoing in the life sciences.

RESULTS

Here we identified a new family of Tc1/mariner transposons, named Incomer (IC), which is close to, but distinct from the known family DD34E/Tc1. ICs have a total length of about 1.2 kb, and harbor a single open reading frame encoding a ~ 346 amino acid transposase with a DD36E motif and flanked by short terminal inverted repeats (TIRs) (22-32 base pairs, bp). This family is absent from prokaryotes, and is mainly distributed among vertebrates (141 species of four classes), including Agnatha (one species of jawless fish), Actinopterygii (132 species of ray-finned fish), Amphibia (four species of frogs), and Mammalia (four species of bats), but have a restricted distribution in invertebrates (four species in Insecta and nine in Arachnida). All ICs in bats (Myotis lucifugus, Eptesicus fuscus, Myotis davidii, and Myotis brandtii) are present as truncated copies in these genomes, and most of them are flanked by relatively long TIRs (51-126 bp). High copy numbers of miniature inverted-repeat transposable elements (MITEs) derived from ICs were also identified in bat genomes. Phylogenetic analysis revealed that ICs are more closely related to DD34E/Tc1 than to other families of Tc1/mariner (e.g., DD34D/mariner and DD × D/pogo), and can be classified into four distinct clusters. The host and IC phylogenies and pairwise distance comparisons between RAG1 genes and all consensus sequences of ICs support the idea that multiple episodes of horizontal transfer (HT) of ICs have occurred in vertebrates. In addition, the discovery of intact transposases, perfect TIRs and target site duplications of ICs suggests that this family may still be active in Insecta, Arachnida, frogs, and fish.

CONCLUSIONS

Exploring the diversity of Tc1/mariner transposons and revealing their evolutionary profiles will help provide a better understanding of the evolution of DNA transposons and their impact on genomic evolution. Here, a newly discovered family (DD36E/Incomer) of Tc1/mariner transposons is described in animals. It displays a similar structural organization and close relationship with the known DD34E/Tc1 elements, but has a relatively narrow distribution, indicating that DD36E/IC might have originated from the DD34E/Tc1 family. Our data also support the hypothesis of horizontal transfer of IC in vertebrates, even invading one lineage of mammals (bats). This study expands our understanding of the diversity of Tc1/mariner transposons and updates the classification of this superfamily.

Loss and gain of sexual reproduction in the same stick insect

The outcome of competition between different reproductive strategies within a single species can be used to infer selective advantage of the winning strategy. Where multiple populations have independently lost or gained sexual reproduction it is possible to investigate whether the advantage is contingent on local conditions. In the New Zealand stick insect Clitarchus hookeri, three populations are distinguished by recent change in reproductive strategy and we determine their likely origins. One parthenogenetic population has established in the United Kingdom and we provide evidence that sexual reproduction has been lost in this population. We identify the sexual population from which the parthenogenetic population was derived, but show that the UK females have a post‐mating barrier to fertilisation. We also demonstrate that two sexual populations have recently arisen in New Zealand within the natural range of the mtDNA lineage that otherwise characterizes parthenogenesis in this species. We infer independent origins of males at these two locations using microsatellite genotypes. In one population, a mixture of local and nonlocal alleles suggested males were the result of invasion. Males in another population were most probably the result of loss of an X chromosome that produced a male phenotype in situ. Two successful switches in reproductive strategy suggest local competitive advantage for outcrossing over parthenogenetic reproduction. Clitarchus hookeri provides remarkable evidence of repeated and rapid changes in reproductive strategy, with competitive outcomes dependent on local conditions.

Two Lys-vasopressin-like peptides, EFLamide, and other phasmid neuropeptides

Phasmid neuropeptide genes were identified in the genomes of two phasmids, Timema cristinae and Clitarchus hookeri. The two species belong to two sisters groups, the Timematodea and Euphasmatodea respectively. Neuropeptide genes were identified using the BLAST+ program on the genome assemblies and the absence of some neuropeptides was confirmed by the concomitant absence of their G-protein coupled receptors. Both genomes were assembled using short reads and the average coverage of the genome is more than 166 times for both species. This makes it virtually impossible that there would not be a single short read for at least one of the conserved transmembrane regions of a GPCR coded by such a genome. Hence, when not a single read can be found for a specific GPCR, it can be concluded that the particular gene is absent from that species. Most previously identified insect neuropeptides are used by these two species. Of the three arthropod allatostatin C related peptides, only allatostatins CC and CCC are present. Both species lack leucokinin, while sulfakinin and dilp8 signaling is absent from Clitarchus, but present in Timema. Interestingly, whereas Timema has lost a vasopressin-related peptide, the gene coding such a peptide is amplified in the Clitarchus genome. Furthermore, while Clitarchus has a specific tryptopyrokinin gene, Timema does not and in this species tryptopyrokinin is coded only by the pyrokinin and periviscerokinin genes. Finally, both species have genes coding EFLamide and its GPCR; in phasmids these genes codes for one (Clitarchus) or two (Timema) EFLamide paracopies.

Sex loss in insects: causes of asexuality and consequences for genomes

Boasting a staggering diversity of reproductive strategies, insects provide attractive models for the comparative study of the causes and consequences of transitions to asexuality. We provide an overview of some contemporary studies of reproductive systems in insects and compile an initial database of asexual insect genome resources. Insect systems have already yielded some important insights into various mechanisms by which sex is lost, including genetic, endosymbiont-mediated, and hybridization. Studies of mutation and substitution after loss of sex provide the strongest empirical support for hypothesized effects of asexuality, whereas there is mixed evidence for ecological hypotheses such as increased parasite load and altered niche breadth in asexuals. Most hypotheses have been explored in a select few taxa (e.g. stick insects, aphids), such that much of the great taxonomic breadth of insects remain understudied. Given the variation in the proximate causes of asexuality in insects, we argue for expanding the taxonomic breadth of study systems. Despite some challenges for investigating sex in insects, the increasing cost-effectiveness of genomic sequencing makes data generation for closely-related asexual and sexual lineages increasingly feasible.

The Draft Genome of the Invasive Walking Stick, Medauroidea extradendata, Reveals Extensive Lineage-Specific Gene Family Expansions of Cell Wall Degrading Enzymes in Phasmatodea

Plant cell wall components are the most abundant macromolecules on Earth. The study of the breakdown of these molecules is thus a central question in biology. Surprisingly, plant cell wall breakdown by herbivores is relatively poorly understood, as nearly all early work focused on the mechanisms used by symbiotic microbes to breakdown plant cell walls in insects such as termites. Recently, however, it has been shown that many organisms make endogenous cellulases. Insects, and other arthropods, in particular have been shown to express a variety of plant cell wall degrading enzymes in many gene families with the ability to break down all the major components of the plant cell wall. Here we report the genome of a walking stick, Medauroidea extradentata, an obligate herbivore that makes uses of endogenously produced plant cell wall degrading enzymes. We present a draft of the 3.3Gbp genome along with an official gene set that contains a diversity of plant cell wall degrading enzymes. We show that at least one of the major families of plant cell wall degrading enzymes, the pectinases, have undergone a striking lineage-specific gene family expansion in the Phasmatodea. This genome will be a useful resource for comparative evolutionary studies with herbivores in many other clades and will help elucidate the mechanisms by which metazoans breakdown plant cell wall components.