The Easter Egg Weevil (Pachyrhynchus) genome reveals syntenic patterns in Coleoptera across 200 million years of evolution

Genomic signatures associated with the evolutionary loss of egg yolk in parasitoid wasps

Trait regression and loss have occurred repeatedly in numerous lineages in response to environmental changes. In parasitoid wasps, a megadiverse group of hymenopteran insects, yolk protein reduction or loss has been observed in many species, likely linked to the transition from ectoparasitism to endoparasitism. However, the genetic basis of this trait and the impact of its loss on genome evolution remain poorly understood. Here, we performed a comparative genomic analysis of 64 hymenopteran insects. The conserved insect yolk protein gene vitellogenin (Vg) underwent five independent loss events in four families, involving 23 of the analyzed endoparasitoid species. Whole-genome alignment suggested that Vg loss occurred during genome rearrangement events. Analysis of Vg receptor gene (VgR) loss, selection, and structural variation in lineages lacking Vg demonstrated functional biases in the patterns of gene loss. The ectoparasitism to endoparasitism transition did not appear to be the primary driver of Vg loss or the subsequent VgR evolution. Together, these findings reveal the genomic changes underlying a unique trait loss in parasitoid wasps. More broadly, this study enhances our understanding of yolk protein loss evolution outside the class Mammalia, highlighting a potential evolutionary trend arising from the availability of an alternative nutrient source for embryonic development.

Cytogenetics of insects in the era of chromosome-level genome assemblies

Over the past few years, a revolution has occurred in cytogenetics, driven by the emergence and spread of methods for obtaining high-quality chromosome-level genome assemblies. In fact, this has led to a new tool for studying chromosomes and chromosomal rearrangements, and this tool is thousands of times more powerful than light microscopy. This tool has revolutionized the cytogenetics of many groups of insects for which previously karyotype information, if available at all, was limited to the chromosome number. Even more impressive are the achievements of the genomic approach for studying the general patterns of chromosome organization and evolution in insects. Thus, it has been shown that rapid transformations of chromosomal numbers, which are often found in the order Lepidoptera, are most often carried out in the most parsimonious way, as a result of simple fusions and fissions of chromosomes. It has been established that these fusions and fissions are not random and occur independently in different phylogenetic lineages due to the reuse of the same ancestral chromosomal breakpoints. It has been shown that the tendency for chromosome fissions is correlated with the presence in chromosomes of the so-called interstitial telomeres, i. e. telomere-like structures located not at the ends of chromosomes, but inside them. It has been revealed that, in most insects, telomeric DNA is not just a set of short repeats, but a very long sequence consisting of (TTAGG)n (or other telomeric motifs), regularly and specifically interrupted by retrotransposons, and the telomeric motifs are diverse in terms of their length and nucleotide composition. The number of high-quality chromosome-level genome assemblies available for insects in the GenBank database is growing exponentially and now exceeds a thousand species. Therefore, the exceptional prospects for using genomic data for karyotype analysis are beyond doubt.

Genome Assembly of Pyrocephalus nanus: A Step Toward the Genetic Conservation of the Endangered Little Vermilion Flycatcher of the Galapagos Islands

Incredibly powerful whole genome studies of conservation genetics, evolution, and biogeography become possible for non-model organisms when reference genomes are available. Here, we report the sequence and assembly of the whole genome of the little vermilion flycatcher (Pyrocephalus nanus; family Tyrannidae), which is an endemic, endangered, and declining species of the Galapagos Islands. Using PacBio HiFi reads to assemble long contigs and Hi-C reads for scaffolding, we assembled a genome of 1.07 Gb comprising 267 contigs in 152 scaffolds, scaffold N50 74 M, contig N50 17.8 M, with 98.9% assigned to candidate chromosomal sequences and 99.72% of the BUSCO passeriformes 10,844 single-copy orthologs present. In addition, we used the novel HiFiMiTie pipeline to fully assemble and verify all portions of the mitochondrial genome from HiFi reads, obtaining a mitogenome of 17,151 bases, containing 13 protein-coding genes, 22 tRNAs, 2 rRNAs, two control regions, and a unique structure of control region duplication and repeats. These genomes will be a critical tool for much-needed studies of phylogenetics, population genetics, biogeography, and conservation genetics of Pyrocephalus and related genera. This genome and other studies that use it will be able to provide recommendations for conservation management, taxonomic improvement, and to understand the evolution and diversification of this genus within the Galapagos Islands.

The genome of the invasive and broadly polyphagous Diaprepes root weevil, Diaprepes abbreviatus (Coleoptera), reveals an arsenal of putative polysaccharide-degrading enzymes

The Diaprepes root weevil (DRW), Diaprepes abbreviatus, is a broadly polyphagous invasive pest of agriculture in the southern United States and the Caribbean. Its genome was sequenced, assembled, and annotated to study genomic correlates of specialized plant-feeding and invasiveness and to facilitate the development of new methods for DRW control. The 1.69 Gb D. abbreviatus genome assembly was distributed across 653 contigs, with an N50 of 7.8 Mb and the largest contig of 62 Mb. Most of the genome was comprised of repetitive sequences, with 66.17% in transposable elements, 5.75% in macrosatellites, and 2.06% in microsatellites. Most expected orthologous genes were present and fully assembled, with 99.5% of BUSCO genes present and 1.5% duplicated. One hundred and nine contigs (27.19 Mb) were identified as putative fragments of the X and Y sex chromosomes, and homology assessment with other beetle X chromosomes indicated a possible sex chromosome turnover event. Genome annotation identified 18,412 genes, including 43 putative horizontally transferred (HT) loci. Notably, 258 genes were identified from gene families known to encode plant cell wall degrading enzymes and invertases, including carbohydrate esterases, polysaccharide lyases, and glycoside hydrolases (GH). GH genes were unusually numerous, with 239 putative genes representing 19 GH families. Interestingly, several other beetle species with large numbers of GH genes are (like D. abbreviatus) successful invasive pests of agriculture or forestry.

The de novo genome of the Black-necked Snakefly (Venustoraphidia nigricollis Albarda, 1891): A resource to study the evolution of living fossils

Snakeflies (Raphidioptera) are the smallest order of holometabolous insects that have kept their distinct and name-giving appearance since the Mesozoic, probably since the Jurassic, and possibly even since their emergence in the Carboniferous, more than 300 million years ago. Despite their interesting nature and numerous publications on their morphology, taxonomy, systematics, and biogeography, snakeflies have never received much attention from the general public, and only a few studies were devoted to their molecular biology. Due to this lack of molecular data, it is therefore unknown, if the conserved morphological nature of these living fossils translates to conserved genomic structures. Here, we present the first genome of the species and of the entire order of Raphidioptera. The final genome assembly has a total length of 669 Mbp and reached a high continuity with an N50 of 5.07 Mbp. Further quality controls also indicate a high completeness and no meaningful contamination. The newly generated data was used in a large-scaled phylogenetic analysis of snakeflies using shared orthologous sequences. Quartet score and gene concordance analyses revealed high amounts of conflicting signals within this group that might speak for substantial incomplete lineage sorting and introgression after their presumed re-radiation after the asteroid impact 66 million years ago. Overall, this reference genome will be a door-opening dataset for many future research applications, and we demonstrated its utility in a phylogenetic analysis that provides new insights into the evolution of this group of living fossils.

Metagenomic binning of PacBio HiFi data prior to assembly reveals a complete genome of Cosmopolites sordidus (Germar) (Coleopterea: Curculionidae, Dryophthorinae) the most damaging arthropod pest of bananas and plantains

PacBio HiFi sequencing was employed in combination with metagenomic binning to produce a high-quality reference genome of Cosmopolites sordidus. We compared k-mer and alignment reference based pre-binning and post-binning approaches to remove contamination. We were also interested to know if the post-binning approach had interspersed bacterial contamination within intragenic regions of Arthropoda binned contigs. Our analyses identified 3,433 genes that were composed with reads identified as of putative bacterial origins. The pre-binning approach yielded a C. sordidus genome of 1.07 Gb genome composed of 3,089 contigs with 98.6% and 97.1% complete and single copy genome and protein BUSCO scores respectively. In this article we demonstrate that in this case the pre-binning approach does not sacrifice assembly quality for more stringent metagenomic filtering. We also determine post-binning allows for increased intragenic contamination increased with increasing coverage, but the frequency of gene contamination increased with lower coverage. Future work should focus on developing reference free pre-binning approaches for HiFi reads produced from eukaryotic based metagenomic samples.

Unraveling the complex evolutionary history of lepidopteran chromosomes through ancestral chromosome reconstruction and novel chromosome nomenclature

BACKGROUND

Lepidoptera is one of the most species-rich animal groups, with substantial karyotype variations among species due to chromosomal rearrangements. Knowledge of the evolutionary patterns of lepidopteran chromosomes still needs to be improved.

RESULTS

Here, we used chromosome-level genome assemblies of 185 lepidopteran insects to reconstruct an ancestral reference genome and proposed a new chromosome nomenclature. Thus, we renamed over 5000 extant chromosomes with this system, revealing the historical events of chromosomal rearrangements and their features. Additionally, our findings indicate that, compared with autosomes, the Z chromosome in Lepidoptera underwent a fast loss of conserved genes, rapid acquisition of lineage-specific genes, and a low rate of gene duplication. Moreover, we presented evidence that all available 67 W chromosomes originated from a common ancestor chromosome, with four neo-W chromosomes identified, including one generated by fusion with an autosome and three derived through horizontal gene transfer. We also detected nearly 4000 inter-chromosomal gene movement events. Notably, Geminin is transferred from the autosome to the Z chromosome. When located on the autosome, Geminin shows female-biased expression, but on the Z chromosome, it exhibits male-biased expression. This contributes to the sexual dimorphism of body size in silkworms.

CONCLUSIONS

Our study sheds light on the complex evolutionary history of lepidopteran chromosomes based on ancestral chromosome reconstruction and novel chromosome nomenclature.

Genome sequencing provides insights into the evolution of gene families encoding plant cell wall-degrading enzymes in longhorned beetles

With more than 36,000 species, the longhorned beetles (family Cerambycidae) are a mega-diverse lineage of mostly xylophagous insects, all of which are represented by the sole sequenced genome of the Asian longhorned beetle (Anoplophora glabripennis; Lamiinae). Their successful radiation has been linked to their ability to degrade plant cell wall components using a range of so-called plant cell wall-degrading enzymes (PCWDEs). Our previous analysis of larval gut transcriptomes demonstrated that cerambycid beetles horizontally acquired genes encoding PCWDEs from various microbial donors; these genes evolved through multiple duplication events to form gene families. To gain further insights into the evolution of these gene families during the Cerambycidae radiation, we assembled draft genomes for four beetle species belonging to three subfamilies using long-read nanopore sequencing. All the PCWDE-encoding genes we annotated from the corresponding larval gut transcriptomes were present in these draft genomes. We confirmed that the newly discovered horizontally acquired glycoside hydrolase family 7 (GH7), subfamily 26 of GH43 (GH43_26), and GH53 (all of which are absent from the A. glabripennis genome) were indeed encoded by these beetles' genome. Most of the PCWDE-encoding genes of bacterial origin gained introns after their transfer into the beetle genome. Altogether, we show that draft genome assemblies generated from nanopore long-reads offer meaningful information to the study of the evolution of gene families in insects. We anticipate that our data will support studies aiming to better understand the biology of the Cerambycidae and other beetles in general.

Chromosomal-level Genome Assembly of the Coffee Bee Hawk Moth Reveals the Evolution of Chromosomes and the Molecular Basis of Distinct Phenotypes

Cephonodes hylas, the coffee bee hawk moth is a hawk moth species with unique characteristics, such as larvae feeding on gardenia, overcoming the toxicity of its iridoid glycosides, diurnal adults, and transparent wings. Although C. hylas is a fascinating model for molecular biological research, genome sequence analysis-based genetic approaches to elucidate these peculiarities have not yet been undertaken. We successfully achieved de novo genome assembly at the chromosome level of C. hylas comparable to the Lepidoptera model organism, silkworm. Additionally, 16,854 protein-coding genes were annotated, and the constructed genome sequence and annotated genes were of the highest quality BUSCO completion compared to closely related species. Comparative genome analysis revealed the process of chromosomal evolution from the Bombycoidea ancestral (n = 31) genome and changes in turnover at the chromosome level associated with chromosomal fusion events, such as the rate of repetitive sequence insertion. These analyses were only possible because the genome was constructed at the chromosome level. Additionally, increased the nonsynonymous/synonymous rate (dN/dS) ratios were observed in multiple photoreceptor-related genes that were strongly associated with the acquisition of diurnal activity. Furthermore, tandemly duplicated expanded genes containing many digestive and other enzymes and larval midgut-specific expression were also confirmed. These genes may be involved in the metabolism of genipin, a toxin found in gardenias. Using the genome sequence of C. hylas determined at the chromosome level, we have successfully identified new insights into the chromosomal evolution of Bombycoidea, as well as the relationship between the genome sequence and its characteristic traits.

A chromosome-scale genome assembly and evaluation of mtDNA variation in the willow leaf beetle Chrysomela aeneicollis

The leaf beetle Chrysomela aeneicollis has a broad geographic range across Western North America but is restricted to cool habitats at high elevations along the west coast. Central California populations occur only at high altitudes (2,700-3,500 m) where they are limited by reduced oxygen supply and recent drought conditions that are associated with climate change. Here, we report a chromosome-scale genome assembly alongside a complete mitochondrial genome and characterize differences among mitochondrial genomes along a latitudinal gradient over which beetles show substantial population structure and adaptation to fluctuating temperatures. Our scaffolded genome assembly consists of 21 linkage groups; one of which we identified as the X chromosome based on female/male whole genome sequencing coverage and orthology with Tribolium castaneum. We identified repetitive sequences in the genome and found them to be broadly distributed across all linkage groups. Using a reference transcriptome, we annotated a total of 12,586 protein-coding genes. We also describe differences in putative secondary structures of mitochondrial RNA molecules, which may generate functional differences important in adaptation to harsh abiotic conditions. We document substitutions at mitochondrial tRNA molecules and substitutions and insertions in the 16S rRNA region that could affect intermolecular interactions with products from the nuclear genome. This first chromosome-level reference genome will enable genomic research in this important model organism for understanding the biological impacts of climate change on montane insects.

Insight into weevil biology from a reference quality genome of the boll weevil, Anthonomus grandis grandis Boheman (Coleoptera: Curculionidae)

The boll weevil, Anthonomus grandis grandis Boheman, is one of the most historically impactful insects due to its near destruction of the US cotton industry in the early 20th century. Contemporary efforts to manage this insect primarily use pheromone baited traps for detection and organophosphate insecticides for control, but this strategy is not sustainable due to financial and environmental costs. We present a high-quality boll weevil genome assembly, consisting of 306 scaffolds with approximately 24,000 annotated genes, as a first step in the identification of gene targets for novel pest control. Gene content and transposable element distribution are similar to those found in other Curculionidae genomes; however, this is the most contiguous and only assembly reported to date for a member in the species-rich genus Anthonomus. Transcriptome profiles across larval, pupal, and adult life stages led to identification of several genes and gene families that could present targets for novel control strategies.

Chromosome-level genome assembly and population genomic analyses provide insights into adaptive evolution of the red turpentine beetle, Dendroctonus valens

BACKGROUND

Biological invasions are responsible for substantial environmental and economic losses. The red turpentine beetle (RTB), Dendroctonus valens LeConte, is an important invasive bark beetle from North America that has caused substantial tree mortality in China. The lack of a high-quality reference genome seriously limits deciphering the extent to which genetic adaptions resulted in a secondary pest becoming so destructive in its invaded area.

RESULTS

Here, we present a 322.41 Mb chromosome-scale reference genome of RTB, of which 98% of assembled sequences are anchored onto fourteen linkage groups including the X chromosome with a N50 size of 24.36 Mb, which is significantly greater than other Coleoptera species. Repetitive sequences make up 45.22% of the genome, which is higher than four other Coleoptera species, i.e., Mountain pine beetle Dendroctonus ponderosae, red flour beetle Tribolium castaneum, blister beetle Hycleus cichorii, and Colorado potato beetle Leptinotarsa decemlineata. We identify rapidly expanded gene families and positively selected genes in RTB, which may be responsible for its rapid environmental adaptation. Population genetic structure of RTB was revealed by genome resequencing of geographic populations in native and invaded regions, suggesting substantial divergence of the North American population and illustrates the possible invasion and spread route in China. Selective sweep analysis highlighted the enhanced ability of Chinese populations in environmental adaptation.

CONCLUSIONS

Overall, our high-quality reference genome represents an important resource for genomics study of invasive bark beetles, which will facilitate the functional study and decipher mechanism underlying invasion success of RTB by integrating the Pinus tabuliformis genome.

The genome of the forest insect pest Pissodes strobi reveals genome expansion and evidence of a Wolbachia endosymbiont

The highly diverse insect family of true weevils, Curculionidae, includes many agricultural and forest pests. Pissodes strobi, commonly known as the spruce weevil or white pine weevil, is a major pest of spruce and pine forests in North America. Pissodes strobi larvae feed on the apical shoots of young trees, causing stunted growth and can destroy regenerating spruce or pine forests. Here, we describe the nuclear and mitochondrial Pissodes strobi genomes and their annotations, as well as the genome of an apparent Wolbachia endosymbiont. We report a substantial expansion of the weevil nuclear genome, relative to other Curculionidae species, possibly driven by an abundance of class II DNA transposons. The endosymbiont observed belongs to a group (supergroup A) of Wolbachia species that generally form parasitic relationships with their arthropod host.