Among over 2,000 species of mealybugs (Hemiptera: Pseudococcidae), only 13 genomes have been published so far, seriously limiting the researches on the phylogeny and adaptive evolution of this group. The continuous publication of mealybug genomes will significantly facilitate our exploration of the biological characteristics, detrimental attributes, and control strategies of the Pseudococcidae family. Jack Beardsley mealybug (Pseudococcus jackbeardsleyi) as one of the hazardous invasive pests, it could cause enormous losses to the fruit and vegetable industries worldwide. Herein, we combined Nanopore long-read, short-read Illumina and Hi-C sequencing, generating a high-quality chromosome-level genome assembly of P. jackbeardsleyi. The genome size was determined to be 334.818 Mb, which was assembled into 5 linkage groups with a N50 of 67.233 Mb. The BUSCO analysis demonstrated the completeness of the genome assembly and annotation are 95.7% and 92.8%, respectively. The developed high-quality genome will serve as an asset for delving into the genetic mechanisms underlying the invasiveness of P. jackbeardsleyi, thereby offering a crucial theoretical foundation for the prevention and management of Pseudococcidae pests.

The Stag beetle (Coleoptera: Lucanidae) is a fascinating group, often considered one of the most primitive within the Scarabaeoidea. They are valuable models for studying beetle evolution. However, the lack of high-quality genomes hinders our understanding of the evolution and ecology of Lucanidae. In this study, we present a chromosome-level genome of Serrognathus titanus by combining PacBio HiFi long reads, Illumina short reads, and Hi-C data. The genome spans 384.07 Mb, with a scaffold N50 size of 75.81 Mb, and most contigs (97.45%, 374.30 Mb) were anchored into six chromosomes. Our BUSCO analysis of the assembly indicates a completeness of 97.6% (n = 1,367), with 92.8% single-copy BUSCOs and 4.8% duplicated BUSCOs identified. Additionally, we found that the genome contains 43.87% (168.50 Mb) repeat elements and identified 14,263 predicted protein-coding genes. The high-quality genome of S. titanus provides valuable genomic information for comprehending the evolution and ecology of Lucanidae.

The higher classification of termites requires substantial revision as the Neoisoptera, the most diverse termite lineage, comprise many paraphyletic and polyphyletic higher taxa. Here, we produce an updated termite classification using genomic-scale analyses. We reconstruct phylogenies under diverse substitution models with ultraconserved elements analyzed as concatenated matrices or within the multi-species coalescence framework. Our classification is further supported by analyses controlling for rogue loci and taxa, and topological tests. We show that the Neoisoptera are composed of seven family-level monophyletic lineages, including the Heterotermitidae Froggatt, Psammotermitidae Holmgren, and Termitogetonidae Holmgren, raised from subfamilial rank. The species-rich Termitidae are composed of 18 subfamily-level monophyletic lineages, including the new subfamilies Crepititermitinae, Cylindrotermitinae, Forficulitermitinae, Neocapritermitinae, Protohamitermitinae, and Promirotermitinae; and the revived Amitermitinae Kemner, Microcerotermitinae Holmgren, and Mirocapritermitinae Kemner. Building an updated taxonomic classification on the foundation of unambiguously supported monophyletic lineages makes it highly resilient to potential destabilization caused by the future availability of novel phylogenetic markers and methods. The taxonomic stability is further guaranteed by the modularity of the new termite classification, designed to accommodate as-yet undescribed species with uncertain affinities to the herein delimited monophyletic lineages in the form of new families or subfamilies.

Episyrphus balteatus can provide dual ecosystem services including pest control and pollination, which the larvae are excellent predators of aphid pest whereas adults are efficient pollinator. In this study, we assembled a high-quality genome of E. balteatus from northern China geographical population at the chromosome level by using Illumina, PacBio long reads, and Hi-C technologies. The 467.42 Mb genome was obtained from 723 contigs, with a contig N50 of 9.16 Mb and Scaffold N50 of 118.85 Mb, and 90.25% (431.75 Mb) of the assembly was anchored to 4 pseudo-autosomes and one pseudo-heterosome. In total, 14,848 protein-coding genes were annotated, and 95.14% of genes were fully represented in NR, GO, KEGG databases. Besides, we also obtained the mitochondrial genome of E. balteatus of 16, 837 bp in length with 37 typical mitochondrial genes. Overall, this high-quality genome is valuable for evolutionary and genetic studies of E. balteatus and other Syrphidae hoverfly species.

To reveal the molecular function of elongation family of very long chain fatty acids(ELO) protein in Cyrtotrachelus buqueti, we have identified 15 ELO proteins from C.buqueti genome. 15 CbuELO proteins were located on four chromosomes. Their isoelectric points ranged from 9.22 to 9.68, and they were alkaline. These CbuELO proteins were stable and hydrophobic. CbuELO proteins had transmembrane movement, and had multiple phosphorylation sites. The secondary structure of CbuELO proteins was mainly α-helix. A total of 10 conserved motifs were identified in CbuELO protein family. Phylogenetic analysis showed that molecular evolutionary relationships of ELO protein family between C. buqueti and Tribolium castaneum was the closest. Developmental transcriptome analysis indicated that CbuELO10, CbuELO13 and CbuELO02 genes were key enzyme genes that determine the synthesis of very long chain fatty acids in pupae and eggs, CbuELO6 and CbuELO7 were that in the male, and CbuELO8 and CbuELO11 were that in the larva. Transcriptome analysis under different temperature conditions indicated that CbuELO1, CbuELO5, CbuELO12 and CbuELO14 participated in regulating temperature stress responses. Transcriptome analysis at different feeding times showed CbuELO12 gene expression level in all feeding time periods was significant downregulation. The qRT-PCR experiment verified expression level changes of CbuELO gene family under different temperature and feeding time conditions. Protein-protein interaction analysis showed that 9 CbuELO proteins were related to each other, CbuELO1, CbuELO4 and CbuELO12 had more than one interaction relationship. These results lay a theoretical foundation for further studying its molecular function during growth and development of C. buqueti.

Sexual dimorphism arises because of divergent fitness optima between the sexes. Phenotypic divergence between sexes can range from mild to extreme. Fireflies, bioluminescent beetles, present various degrees of sexual dimorphism, with species showing very mild sexual dimorphism to species presenting female-specific neoteny, posing a unique framework to investigate the evolution of sexually dimorphic traits across species. In this work, we present novel assembled genomes of two firefly species, Lamprohiza splendidula and Luciola italica, species with different degrees of sexual dimorphism. We uncover high synteny conservation of the X-chromosome across ~ 180 Mya and find full X-chromosome dosage compensation in our two fireflies, hinting at common mechanism upregulating the single male X-chromosome. Different degrees of sex-biased expressed genes were found across two body parts showing different proportions of expression conservation between species. Interestingly, we do not find X-chromosome enrichment of sex-biased genes, but retrieve autosomal enrichment of sex-biased genes. We further uncover higher nucleotide diversity in the intronic regions of sex-biased genes, hinting at a maintenance of heterozygosity through sexual selection. We identify different levels of sex-biased gene expression divergence including a set of genes showing conserved sex-biased gene expression between species. Divergent and conserved sex-biased genes are good candidates to test their role in the maintenance of sexually dimorphic traits.

Chromosome rearrangements may distort 3D chromatin architectures and thus change gene regulation, yet how 3D chromatin structures evolve in insects is largely unknown. Here, we obtain chromosome-level genomes for four butterfly species, Graphium cloanthus, Graphium sarpedon, Graphium eurypylus with 2n = 30, 40, and 60, respectively, and Papilio bianor with 2n = 60. Together with large-scale Hi-C data, we find that inter-chromosome rearrangements very rarely disrupted the pre-existing 3D chromatin structure of ancestral chromosomes. However, some intra-chromosome rearrangements changed 3D chromatin structures compared to the ancestral configuration. We find that new TADs and subTADs have emerged across the rearrangement sites where their adjacent compartments exhibit uniform types. Two intra-chromosome rearrangements altered Rel and lft regulation, potentially contributing to wing patterning differentiation and host plant choice. Notably, butterflies exhibited chromatin loops between Hox gene cluster ANT-C and BX-C, unlike Drosophila. Our CRISPR-Cas9 experiments in butterflies confirm that knocking out the CTCF binding site of the loops in BX-C affected the phenotypes regulated by Antp in ANT-C, resulting in legless larva. Our results reveal evolutionary patterns of insect 3D chromatin structures and provide evidence that 3D chromatin structure changes can play important roles in the evolution of traits.

The pink stem borer, Sesamia inferens Walker (Lepidoptera: Noctuidae), is one of the most notorious pest insects of rice and maize crops in the world. Here, we generated a high-quality chromosome-level genome assembly of S. inferens, using a combination of Illumina, PacBio HiFi and Hi-C technologies. The total assembly size was 973.18 Mb with a contig N50 of 33.39 Mb, anchored to 31 chromosomes, revealing a karyotype of 30 + Z. The BUSCO analysis indicated a high completeness of 98.90% (n = 5286), including 5172 (97.8%) single-copy BUSCOs and 58 (1.1%) duplicated BUSCOs. The genome contains 58.59% (564.58 Mb) repeat elements and 26628 predicted protein-coding genes. The chromosome-level genome assembly of S. inferens provides in-depth knowledge and will be a helpful resource for the Lepidoptera and pest control research communities.

Lucanidae (Coleoptera: Scarabaeidae) are fascinating beetles exhibiting significant dimorphism and are widely used as beetle evolutionary study models. However, lacking high-quality genomes prohibits our understanding of Lucanidae. Herein, we proposed a chromosome-level genome assembly of a widespread species, Prosopocoilus inquinatus, combining PacBio HiFi, Illumina, and Hi-C data. The genome size reaches 649.73 Mb, having the scaffold N50 size of 59.50 Mb, and 99.6% (647.13 Mb) of the assembly successfully anchored on 12 chromosomes. The BUSCO analysis of the genome exhibits a completeness of 99.6% (n = 1,367), including 1,362 (98.5%) single-copy BUSCOs and 15 (1.1%) duplicated BUSCOs. The genome annotation identifies that the genome contains 61.41% repeat elements and 13,452 predicted protein-coding genes. This high-quality Lucanidae genome provides treasured genomic information to our knowledge of stag beetles.

Slavum lentiscoides and Chaetogeoica ovagalla are two aphid species from the subtribe Fordina of Fordini within the subfamily Eriosomatinae, and they produce galls on their primary host plants Pistacia. We assembled chromosome-level genomes of these two species using Nanopore long-read sequencing and Hi-C technology. A 332 Mb genome assembly of S. lentiscoides with a scaffold N50 of 19.77 Mb, including 11,747 genes, and a 289 Mb genome assembly of C. ovagalla with a scaffold N50 of 11.85 Mb, containing 14,492 genes, were obtained. The Benchmarking Universal Single-Copy Orthologs (BUSCO) benchmark of the two genome assemblies reached 93.7% (91.9% single-copy) and 97.0% (95.3% single-copy), respectively. The high-quality genome assemblies in our study provide valuable resources for future genomic research of galling aphids.