Enhanced genome assembly and a new official gene set for Tribolium castaneum

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.

BeetleAtlas: An Ontogenetic and Tissue-specific Transcriptomic Atlas of the Red Flour Beetle Tribolium castaneum

The red flour beetle Tribolium castaneum has emerged as a powerful model in insect functional genomics. However, a major limitation in the field is the lack of a detailed spatio-temporal view of the genetic signatures underpinning the function of distinct tissues and life stages. Here, we present an ontogenetic and tissue-specific web-based resource for Tribolium transcriptomics: BeetleAtlas (https://www.beetleatlas.org). This web application provides access to a database populated with quantitative expression data for nine adult and seven larval tissues, as well as for four embryonic stages of Tribolium. BeetleAtlas allows one to search for individual Tribolium genes to obtain values of both total gene expression and enrichment in different tissues, together with data for individual isoforms. To facilitate cross-species studies, one can also use Drosophila melanogaster gene identifiers to search for related Tribolium genes. For retrieved genes there are options to identify and display the tissue expression of related Tribolium genes or homologous Drosophila genes. Five additional search modes are available to find genes conforming to any of the following criteria: exhibiting high expression in a particular tissue; showing significant differences in expression between larva and adult; having a peak of expression at a specific stage of embryonic development; belonging to a particular functional category; and displaying a pattern of tissue expression similar to that of a query gene. We illustrate how the different feaures of BeetleAtlas can be used to illuminate our understanding of the genetic mechanisms underpinning the biology of what is the largest animal group on earth.

Resources modulate developmental shifts but not infection tolerance upon coinfection in an insect system

Energetic resources fuel immune responses and parasite growth within organisms, but it is unclear whether energy allocation is sufficient to explain changes in infection outcomes under the threat of multiple parasites. We manipulated diet in flour beetles (Tribolium confusum) infected with two natural parasites to investigate the role of resources in shifting metabolic and immune responses after single and co-infection. Our results suggest that gregarine parasites alter the within-host energetic environment, and by extension juvenile development time, in a diet-dependent manner. Gregarines do not affect host resistance to acute bacterial infection but do stimulate the expression of an alternative set of immune genes and promote damage to the gut, ultimately contributing to reduced survival regardless of diet. Thus, energy allocation is not sufficient to explain the immunological contribution to coinfection outcomes, emphasizing the importance of mechanistic insight for predicting the impact of coinfection across levels of biological organization.

A chromosome-level genome assembly of Prosopocoilus inquinatus Westwood, 1848 (Coleoptera: Lucanidae)

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.

Tissue-Level Integration Overrides Gradations of Differentiating Cell Identity in Beetle Extraembryonic Tissue

During animal embryogenesis, one of the earliest specification events distinguishes extraembryonic (EE) from embryonic tissue fates: the serosa in the case of the insects. While it is well established that the homeodomain transcription factor Zen1 is the critical determinant of the serosa, the subsequent realization of this tissue's identity has not been investigated. Here, we examine serosal differentiation in the beetle Tribolium castaneum based on the quantification of morphological and morphogenetic features, comparing embryos from a Tc-zen1 RNAi dilution series, where complete knockdown results in amnion-only EE tissue identity. We assess features including cell density, tissue boundary morphology, and nuclear size as dynamic readouts for progressive tissue maturation. While some features exhibit an all-or-nothing outcome, other key features show dose-dependent phenotypic responses with trait-specific thresholds. Collectively, these findings provide nuance beyond the known status of Tc-Zen1 as a selector gene for serosal tissue patterning. Overall, our approach illustrates how the analysis of tissue maturation dynamics from live imaging extends but also challenges interpretations based on gene expression data, refining our understanding of tissue identity and when it is achieved.

PhyloAln: A Convenient Reference-Based Tool to Align Sequences and High-Throughput Reads for Phylogeny and Evolution in the Omic Era

The current trend in phylogenetic and evolutionary analyses predominantly relies on omic data. However, prior to core analyses, traditional methods typically involve intricate and time-consuming procedures, including assembly from high-throughput reads, decontamination, gene prediction, homology search, orthology assignment, multiple sequence alignment, and matrix trimming. Such processes significantly impede the efficiency of research when dealing with extensive data sets. In this study, we develop PhyloAln, a convenient reference-based tool capable of directly aligning high-throughput reads or complete sequences with existing alignments as a reference for phylogenetic and evolutionary analyses. Through testing with simulated data sets of species spanning the tree of life, PhyloAln demonstrates consistently robust performance compared with other reference-based tools across different data types, sequencing technologies, coverages, and species, with percent completeness and identity at least 50 percentage points higher in the alignments. Additionally, we validate the efficacy of PhyloAln in removing a minimum of 90% foreign and 70% cross-contamination issues, which are prevalent in sequencing data but often overlooked by other tools. Moreover, we showcase the broad applicability of PhyloAln by generating alignments (completeness mostly larger than 80%, identity larger than 90%) and reconstructing robust phylogenies using real data sets of transcriptomes of ladybird beetles, plastid genes of peppers, or ultraconserved elements of turtles. With these advantages, PhyloAln is expected to facilitate phylogenetic and evolutionary analyses in the omic era. The tool is accessible at https://github.com/huangyh45/PhyloAln.

De novo assembly of transcriptomes and differential gene expression analysis using short-read data from emerging model organisms - a brief guide

Many questions in biology benefit greatly from the use of a variety of model systems. High-throughput sequencing methods have been a triumph in the democratization of diverse model systems. They allow for the economical sequencing of an entire genome or transcriptome of interest, and with technical variations can even provide insight into genome organization and the expression and regulation of genes. The analysis and biological interpretation of such large datasets can present significant challenges that depend on the 'scientific status' of the model system. While high-quality genome and transcriptome references are readily available for well-established model systems, the establishment of such references for an emerging model system often requires extensive resources such as finances, expertise and computation capabilities. The de novo assembly of a transcriptome represents an excellent entry point for genetic and molecular studies in emerging model systems as it can efficiently assess gene content while also serving as a reference for differential gene expression studies. However, the process of de novo transcriptome assembly is non-trivial, and as a rule must be empirically optimized for every dataset. For the researcher working with an emerging model system, and with little to no experience with assembling and quantifying short-read data from the Illumina platform, these processes can be daunting. In this guide we outline the major challenges faced when establishing a reference transcriptome de novo and we provide advice on how to approach such an endeavor. We describe the major experimental and bioinformatic steps, provide some broad recommendations and cautions for the newcomer to de novo transcriptome assembly and differential gene expression analyses. Moreover, we provide an initial selection of tools that can assist in the journey from raw short-read data to assembled transcriptome and lists of differentially expressed genes.

The Genome Organization of 5S rRNA Genes in the Model Organism Tribolium castaneum and Its Sibling Species Tribolium freemani

5S ribosomal DNAs (rDNAs) are arranged in tandem and are often under-represented in genome assemblies. In the present study, we performed a global and in-depth analysis of the 5S rDNAs in the model insect Tribolium castaneum and its closely related species Tribolium freemani. To accomplish this goal, we used our recently published genome assemblies based on Nanopore and PacBio long-read sequencing. Although these closely related species share the 5S rRNA gene sequence with high homology, they show a different organization of the 5S rDNA locus. Analysis of 5S rDNA arrays in T. castaneum revealed a typical tandemly repeated organization characterized by repeat units consisting of the 121 bp long 5S rRNA gene and the 71 bp long nontranscribed spacer (NTS). In contrast, T. freemani showed a much more complex organization of 5S rDNA arrays characterized by two patterns. The first is based on the association of 5S rRNA gene with arrays of a satellite DNA, representing the NTS sequence of the 5S rDNA genes in T. freemani. The second, more complex type is characterized by a somewhat less frequent occurrence of the 5S rRNA gene and its association with longer satellite DNA arrays that are regularly interrupted by Jockey-like retrotransposons. This organization, in which the ribosomal gene is associated with two completely different repetitive elements such as satellite DNAs and retrotransposons, suggests that the 5S rRNA gene, regardless of its crucial function in the genome, could be a subject of extremely dynamic genomic rearrangements.

A near-complete genome assembly of Monochamus alternatus a major vector beetle of pinewood nematode

The Japanese sawyer beetle, Monochamus alternatus, is not only one of the most important wood boring pest itself, but also a major vector of the invasive pinewood nematode (PWN), which is the causal agent of the devastative pine wilt disease (PWD) and threats the global pine forest. Here, we present a near-complete genome of M. alternatus at the chromosome level. The assembled genome was 792.05 Mb with contig N50 length of 55.99 Mb, which is the largest N50 size among the sequenced Coleoptera insects currently. 99.57% of sequence was anchored onto ten pseudochromosomes (one X-chromosome and nine autosomes), and the final genome harbored only 13 gaps. BUSCO evaluation revealed the presence of 99.0% of complete core genes. Thus, our genome assembly represented the highest-contiguity genome assembly as well as high completeness in insects so far. We identified 20,471 protein-coding genes, of which 20,070 (98.04%) were functionally annotated. The genome assembly of M. alternatus provides a valuable resource for exploring the evolution of the symbiosis between PWN and the vector insects.

Genomic insights into the chromosomal elongation in a family of Collembola

Collembola is a highly diverse and abundant group of soil arthropods with chromosome numbers ranging from 5 to 11. Previous karyotype studies indicated that the Tomoceridae family possesses an exceptionally long chromosome. To better understand chromosome size evolution in Collembola, we obtained a chromosome-level genome of Yoshiicerus persimilis with a size of 334.44 Mb and BUSCO completeness of 97.0% (n = 1013). Both genomes of Y. persimilis and Tomocerus qinae (recently published) have an exceptionally large chromosome (ElChr greater than 100 Mb), accounting for nearly one-third of the genome. Comparative genomic analyses suggest that chromosomal elongation occurred independently in the two species approximately 10 million years ago, rather than in the ancestor of the Tomoceridae family. The ElChr elongation was caused by large tandem and segmental duplications, as well as transposon proliferation, with genes in these regions experiencing weaker purifying selection (higher dN/dS) than conserved regions. Moreover, inter-genomic synteny analyses indicated that chromosomal fission/fusion events played a crucial role in the evolution of chromosome numbers (ranging from 5 to 7) within Entomobryomorpha. This study provides a valuable resource for investigating the chromosome evolution of Collembola.

scRNA-seq Reveals Novel Genetic Pathways and Sex Chromosome Regulation in Tribolium Spermatogenesis

Spermatogenesis is critical to sexual reproduction yet evolves rapidly in many organisms. High-throughput single-cell transcriptomics promises unparalleled insight into this important process but understanding can be impeded in nonmodel systems by a lack of known genes that can reliably demarcate biologically meaningful cell populations. Tribolium castaneum, the red flour beetle, lacks known markers for spermatogenesis found in insect species like Drosophila melanogaster. Using single-cell sequencing data collected from adult beetle testes, we implement a strategy for elucidating biologically meaningful cell populations by using transient expression stage identification markers, weighted principal component clustering, and SNP-based haploid/diploid phasing. We identify populations that correspond to observable points in sperm differentiation and find species specific markers for each stage. Our results indicate that molecular pathways underlying spermatogenesis in Coleoptera are substantially diverged from those in Diptera. We also show that most genes on the X chromosome experience meiotic sex chromosome inactivation. Temporal expression of Drosophila MSL complex homologs coupled with spatial analysis of potential chromatin entry sites further suggests that the dosage compensation machinery may mediate escape from meiotic sex chromosome inactivation and postmeiotic reactivation of the X chromosome.

Genome evolution and divergence in cis-regulatory architecture is associated with condition-responsive development in horned dung beetles

Phenotypic plasticity is thought to be an important driver of diversification and adaptation to environmental variation, yet the genomic mechanisms mediating plastic trait development and evolution remain poorly understood. The Scarabaeinae, or true dung beetles, are a species-rich clade of insects recognized for their highly diversified nutrition-responsive development including that of cephalic horns-evolutionarily novel, secondary sexual weapons that exhibit remarkable intra- and interspecific variation. Here, we investigate the evolutionary basis for horns as well as other key dung beetle traits via comparative genomic and developmental assays. We begin by presenting chromosome-level genome assemblies of three dung beetle species in the tribe Onthophagini (> 2500 extant species) including Onthophagus taurus, O. sagittarius, and Digitonthophagus gazella. Comparing these assemblies to those of seven other species across the order Coleoptera identifies evolutionary changes in coding sequence associated with metabolic regulation of plasticity and metamorphosis. We then contrast chromatin accessibility in developing head horn tissues of high- and low-nutrition O. taurus males and females and identify distinct cis-regulatory architectures underlying nutrition- compared to sex-responsive development, including a large proportion of recently evolved regulatory elements sensitive to horn morph determination. Binding motifs of known and new candidate transcription factors are enriched in these nutrition-responsive open chromatin regions. Our work highlights the importance of chromatin state regulation in mediating the development and evolution of plastic traits, demonstrates gene networks are highly evolvable transducers of environmental and genetic signals, and provides new reference-quality genomes for three species that will bolster future developmental, ecological, and evolutionary studies of this insect group.

A chromosome-level assembly of the seed beetle Callosobruchus maculatus genome with annotation of its repetitive elements

Callosobruchus maculatus is a major agricultural pest of legume crops worldwide and an established model system in ecology and evolution. Yet, current molecular biological resources for this species are limited. Here, we employ Hi-C sequencing to generate a greatly improved genome assembly and we annotate its repetitive elements in a dedicated in-depth effort where we manually curate and classify the most abundant unclassified repeat subfamilies. We present a scaffolded chromosome-level assembly, which is 1.01 Gb in total length with 86% being contained within the 9 autosomes and the X chromosome. Repetitive sequences accounted for 70% of the total assembly. DNA transposons covered 18% of the genome, with the most abundant superfamily being Tc1-Mariner (9.75% of the genome). This new chromosome-level genome assembly of C. maculatus will enable future genetic and evolutionary studies not only of this important species but of beetles more generally.

Mapping the functional form of the trade-off between infection resistance and reproductive fitness under dysregulated immune signaling

Immune responses benefit organismal fitness by clearing parasites but also exact costs associated with immunopathology and energetic investment. Hosts manage these costs by tightly regulating the induction of immune signaling to curtail excessive responses and restore homeostasis. Despite the theoretical importance of turning off the immune response to mitigate these costs, experimentally connecting variation in the negative regulation of immune responses to organismal fitness remains a frontier in evolutionary immunology. In this study, we used a dose-response approach to manipulate the RNAi-mediated knockdown efficiency of cactus (IκBα), a central regulator of Toll pathway signal transduction in flour beetles (Tribolium castaneum). By titrating cactus activity across four distinct levels, we derived the shape of the relationship between immune response investment and traits associated with host fitness, including infection susceptibility, lifespan, fecundity, body mass, and gut homeostasis. Cactus knock-down increased the overall magnitude of inducible immune responses and delayed their resolution in a dsRNA dose-dependent manner, promoting survival and resistance following bacterial infection. However, these benefits were counterbalanced by dsRNA dose-dependent costs to lifespan, fecundity, body mass, and gut integrity. Our results allowed us to move beyond the qualitative identification of a trade-off between immune investment and fitness to actually derive its functional form. This approach paves the way to quantitatively compare the evolution and impact of distinct regulatory elements on life-history trade-offs and fitness, filling a crucial gap in our conceptual and theoretical models of immune signaling network evolution and the maintenance of natural variation in immune systems.

Multiple Origins of Bioluminescence in Beetles and Evolution of Luciferase Function

Bioluminescence in beetles has long fascinated biologists, with diverse applications in biotechnology. To date, however, our understanding of its evolutionary origin and functional variation mechanisms remains poor. To address these questions, we obtained high-quality reference genomes of luminous and nonluminous beetles in 6 Elateroidea families. We then reconstructed a robust phylogenetic relationship for all luminous families and related nonluminous families. Comparative genomic analyses and biochemical functional experiments suggested that gene evolution within Elateroidea played a crucial role in the origin of bioluminescence, with multiple parallel origins observed in the luminous beetle families. While most luciferase-like proteins exhibited a conserved nonluminous amino acid pattern (TLA346 to 348) in the luciferin-binding sites, luciferases in the different luminous beetle families showed divergent luminous patterns at these sites (TSA/CCA/CSA/LVA). Comparisons of the structural and enzymatic properties of ancestral, extant, and site-directed mutant luciferases further reinforced the important role of these sites in the trade-off between acyl-CoA synthetase and luciferase activities. Furthermore, the evolution of bioluminescent color demonstrated a tendency toward hypsochromic shifts and variations among the luminous families. Taken together, our results revealed multiple parallel origins of bioluminescence and functional divergence within the beetle bioluminescent system.

A life-history allele of large effect shortens developmental time in a wild insect population

Developmental time is a key life-history trait with large effects on Darwinian fitness. In many insects, developmental time is currently under strong selection to minimize ecological mismatches in seasonal timing induced by climate change. The genetic basis of responses to such selection, however, is poorly understood. To address this problem, we set up a long-term evolve-and-resequence experiment in the beetle Tribolium castaneum and selected replicate, outbred populations for fast or slow embryonic development. The response to this selection was substantial and embryonic developmental timing of the selection lines started to diverge during dorsal closure. Pooled whole-genome resequencing, gene expression analysis and an RNAi screen pinpoint a 222 bp deletion containing binding sites for Broad and Tramtrack upstream of the ecdysone degrading enzyme Cyp18a1 as a main target of selection. Using CRISPR/Cas9 to reconstruct this allele in the homogenous genetic background of a laboratory strain, we unravel how this single deletion advances the embryonic ecdysone peak inducing dorsal closure and show that this allele accelerates larval development but causes a trade-off with fecundity. Our study uncovers a life-history allele of large effect and reveals the evolvability of developmental time in a natural insect population.

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

Advances in genome sequencing reveal changes in gene content that contribute to arthropod macroevolution

Current sequencing technology allows for the relatively affordable generation of highly contiguous genomes. Technological advances have made it possible for researchers to investigate the consequences of diverse sorts of genomic variants, such as gene gain and loss. With the extraordinary number of high-quality genomes now available, we take stock of how these genomic variants impact phenotypic evolution. We take care to point out that the identification of genomic variants of interest is only the first step in understanding their impact. Painstaking lab or fieldwork is still required to establish causal relationships between genomic variants and phenotypic evolution. We focus mostly on arthropod research, as this phylum has an impressive degree of phenotypic diversity and is also the subject of much evolutionary genetics research. This article is intended to both highlight recent advances in the field and also to be a primer for learning about evolutionary genetics and genomics.

Next generation marker-based vector concepts for rapid and unambiguous identification of single and double homozygous transgenic organisms

For diploid model organisms, the actual transgenesis processes require subsequent periods of transgene management, which are challenging in emerging model organisms due to the lack of suitable methodology. We used the red flour beetle Tribolium castaneum, a stored-grain pest, to perform a comprehensive functional evaluation of our AClashOfStrings (ACOS) and the combined AGameOfClones/AClashOfStrings (AGOC/ACOS) vector concepts, which use four clearly distinguishable markers to provide full visual control over up to two independent transgenes. We achieved comprehensive statistical validation of our approach by systematically creating seventeen novel single and double homozygous sublines intended for fluorescence live imaging, including several sublines in which the microtubule cytoskeleton is labeled. During the mating procedures, we genotyped more than 20,000 individuals in less than 80 working hours, which corresponds to about 10 to 15 s per individual. We also confirm the functionality of our combined concept in two double transgene special cases, i.e. integration of both transgenes in close proximity on the same chromosome and integration of one transgene on the X allosome. Finally, we discuss our vector concepts regarding performance, genotyping accuracy, throughput, resource saving potential, fluorescent protein choice, modularity, adaptation to other diploid model organisms and expansion capability.

Mapping the functional form of the trade-off between infection resistance and reproductive fitness under dysregulated immune signaling

Immune responses benefit organismal fitness by clearing parasites but also exact costs associated with immunopathology and energetic investment. Hosts manage these costs by tightly regulating the induction of immune signaling to curtail excessive responses and restore homeostasis. Despite the theoretical importance of turning off the immune response to mitigate these costs, experimentally connecting variation in the negative regulation of immune responses to organismal fitness remains a frontier in evolutionary immunology. In this study, we used a dose-response approach to manipulate the RNAi-mediated knockdown efficiency of cactus (IκBα), a central regulator of Toll pathway signal transduction in flour beetles (Tribolium castaneum). By titrating cactus activity along a continuous gradient, we derived the shape of the relationship between immune response investment and traits associated with host fitness, including infection susceptibility, lifespan, fecundity, body mass, and gut homeostasis. Cactus knock-down increased the overall magintude of inducible immune responses and delayed their resolution in a dsRNA dose-dependent manner, promoting survival and resistance following bacterial infection. However, these benefits were counterbalanced by dsRNA dose-dependent costs to lifespan, fecundity, body mass, and gut integrity. Our results allowed us to move beyond the qualitative identification of a trade-off between immune investment and fitness to actually derive its functional form. This approach paves the way to quantitatively compare the evolution and impact of distinct regulatory elements on life-history trade-offs and fitness, filling a crucial gap in our conceptual and theoretical models of immune signaling network evolution and the maintenance of natural variation in immune systems.

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.

Pest status, molecular evolution, and epigenetic factors derived from the genome assembly of Frankliniella fusca, a thysanopteran phytovirus vector

BACKGROUND

The tobacco thrips (Frankliniella fusca Hinds; family Thripidae; order Thysanoptera) is an important pest that can transmit viruses such as the tomato spotted wilt orthotospovirus to numerous economically important agricultural row crops and vegetables. The structural and functional genomics within the order Thysanoptera has only begun to be explored. Within the > 7000 known thysanopteran species, the melon thrips (Thrips palmi Karny) and the western flower thrips (Frankliniella occidentalis Pergrande) are the only two thysanopteran species with assembled genomes.

RESULTS

A genome of F. fusca was assembled by long-read sequencing of DNA from an inbred line. The final assembly size was 370 Mb with a single copy ortholog completeness of ~ 99% with respect to Insecta. The annotated genome of F. fusca was compared with the genome of its congener, F. occidentalis. Results revealed many instances of lineage-specific differences in gene content. Analyses of sequence divergence between the two Frankliniella species' genomes revealed substitution patterns consistent with positive selection in ~ 5% of the protein-coding genes with 1:1 orthologs. Further, gene content related to its pest status, such as xenobiotic detoxification and response to an ambisense-tripartite RNA virus (orthotospovirus) infection was compared with F. occidentalis. Several F. fusca genes related to virus infection possessed signatures of positive selection. Estimation of CpG depletion, a mutational consequence of DNA methylation, revealed that F. fusca genes that were downregulated and alternatively spliced in response to virus infection were preferentially targeted by DNA methylation. As in many other insects, DNA methylation was enriched in exons in Frankliniella, but gene copies with homology to DNA methyltransferase 3 were numerous and fragmented. This phenomenon seems to be relatively unique to thrips among other insect groups.

CONCLUSIONS

The F. fusca genome assembly provides an important resource for comparative genomic analyses of thysanopterans. This genomic foundation allows for insights into molecular evolution, gene regulation, and loci important to agricultural pest status.

Whole genome assemblies of Zophobas morio and Tenebrio molitor

Zophobas morio (=Zophobas atratus) and Tenebrio molitor are darkling beetles with industrial importance due to their use as feeder insects and their apparent ability to biodegrade plastics. High quality genome assemblies were recently reported for both species. Here, we report additional independent Z. morio and T. molitor genome assemblies generated from Nanopore and Illumina data. Following scaffolding against the published genomes, haploid assemblies of 462 Mb (scaffold N90 of 16.8 Mb) and 258 Mb (scaffold N90 of 5.9 Mb) were produced for Z. morio and T. molitor, respectively. Gene prediction led to the prediction of 28,544 and 19,830 genes for Z. morio and T. molitor, respectively. Benchmarking Universal Single Copy Orthologs (BUSCO) analyses suggested that both assemblies have a high level of completeness; 91.5 and 89.0% of the BUSCO endopterygota marker genes were complete in the Z. morio assembly and proteome, respectively, while 99.1 and 92.8% were complete in the T. molitor assembly and proteome, respectively. Phylogenomic analyses of four genera from the family Tenebrionidae yielded phylogenies consistent with those previously constructed based on mitochondrial genomes. Synteny analyses revealed large stretches of macrosynteny across the family Tenebrionidae, as well as numerous within-chromosome rearrangements. Finally, orthogroup analysis identified ∼28,000 gene families across the family Tenebrionidae, of which 8,185 were identified in all five of the analyzed species, and 10,837 were conserved between Z. morio and T. molitor. We expect that the availability of multiple whole genome sequences for Z. morio and T. molitor will facilitate population genetics studies to identify genetic variation associated with industrially relevant phenotypes.

The draft genome sequence of the Japanese rhinoceros beetle Trypoxylus dichotomus septentrionalis towards an understanding of horn formation

The Japanese rhinoceros beetle Trypoxylus dichotomus is a giant beetle with distinctive exaggerated horns present on the head and prothoracic regions of the male. T. dichotomus has been used as a research model in various fields such as evolutionary developmental biology, ecology, ethology, biomimetics, and drug discovery. In this study, de novo assembly of 615 Mb, representing 80% of the genome estimated by flow cytometry, was obtained using the 10 × Chromium platform. The scaffold N50 length of the genome assembly was 8.02 Mb, with repetitive elements predicted to comprise 49.5% of the assembly. In total, 23,987 protein-coding genes were predicted in the genome. In addition, de novo assembly of the mitochondrial genome yielded a contig of 20,217 bp. We also analyzed the transcriptome by generating 16 RNA-seq libraries from a variety of tissues of both sexes and developmental stages, which allowed us to identify 13 co-expressed gene modules. We focused on the genes related to horn formation and obtained new insights into the evolution of the gene repertoire and sexual dimorphism as exemplified by the sex-specific splicing pattern of the doublesex gene. This genomic information will be an excellent resource for further functional and evolutionary analyses, including the evolutionary origin and genetic regulation of beetle horns and the molecular mechanisms underlying sexual dimorphism.

The Low-Copy-Number Satellite DNAs of the Model Beetle Tribolium castaneum

The red flour beetle Tribolium castaneum is an important pest of stored agricultural products and the first beetle whose genome was sequenced. So far, one high-copy-number and ten moderate-copy-number satellite DNAs (satDNAs) have been described in the assembled part of its genome. In this work, we aimed to catalog the entire collection of T. castaneum satDNAs. We resequenced the genome using Illumina technology and predicted potential satDNAs via graph-based sequence clustering. In this way, we discovered 46 novel satDNAs that occupied a total of 2.1% of the genome and were, therefore, considered low-copy-number satellites. Their repeat units, preferentially 140-180 bp and 300-340 bp long, showed a high A + T composition ranging from 59.2 to 80.1%. In the current assembly, we annotated the majority of the low-copy-number satDNAs on one or a few chromosomes, discovering mainly transposable elements in their vicinity. The current assembly also revealed that many of the in silico predicted satDNAs were organized into short arrays not much longer than five consecutive repeats, and some of them also had numerous repeat units scattered throughout the genome. Although 20% of the unassembled genome sequence masked the genuine state, the predominance of scattered repeats for some low-copy satDNAs raises the question of whether these are essentially interspersed repeats that occur in tandem only sporadically, with the potential to be satDNA "seeds".

Analysis of the Toll and Spaetzle Genes Involved in Toll Pathway-Dependent Antimicrobial Gene Induction in the Red Flour Beetle, Tribolium castaneum (Coleoptera; Tenebrionidae)

Insects rely only on their innate immune system to protect themselves from pathogens. Antimicrobial peptide (AMP) production is the main immune reaction in insects. In Drosophila melanogaster, the reaction is regulated mainly by the Toll and immune deficiency (IMD) pathways. Spaetzle proteins, activated by immune signals from upstream components, bind to Toll proteins, thus, activating the Toll pathway, which in turn, induces AMP genes. Previous studies have shown the difference in immune systems related to Toll and IMD pathways between D. melanogaster and Tribolium castaneum. In T. castaneum, nine Toll and seven spaetzle (spz) genes were identified. To extend our understanding of AMP production by T. castaneum, we conducted functional assays of Toll and spaetzle genes related to Toll-pathway-dependent AMP gene expression in T. castaneum under challenge with bacteria or budding yeast. The results revealed that Toll3 and Toll4 double-knockdown and spz7 knockdown strongly and moderately reduced the Toll-pathway-dependent expression of AMP genes, respectively. Moreover, Toll3 and Toll4 double-knockdown pupae more rapidly succumbed to entomopathogenic bacteria than the control pupae, but spz7 knockdown pupae did not. The results suggest that Toll3 and Toll4 play a large role in Toll-pathway-dependent immune reactions, whereas spz7 plays a small part.

Global Gene Expression Analysis Reveals Complex Cuticle Organization of the Tribolium Compound Eye

The red flour beetle Tribolium castaneum is a resource-rich model for genomic and developmental studies. To extend previous studies on Tribolium eye development, we produced transcriptomes for normal-eyed and eye-depleted heads of pupae and adults to identify differentially transcript-enriched (DE) genes in the visual system. Unexpectedly, cuticle-related genes were the largest functional class in the pupal compound eye DE gene population, indicating differential enrichment in three distinct cuticle components: clear lens facet cuticle, highly melanized cuticle of the ocular diaphragm, which surrounds the Tribolium compound eye for internal fortification, and newly identified facet margins of the tanned cuticle, possibly enhancing external fortification. Phylogenetic, linkage, and high-throughput gene knockdown data suggest that most cuticle proteins (CPs) expressed in the Tribolium compound eye stem from the deployment of ancient CP genes. Consistent with this, TcasCPR15, which we identified as the major lens CP gene in Tribolium, is a beetle-specific but pleiotropic paralog of the ancient CPR RR-2 CP gene family. The less abundant yet most likely even more lens-specific TcasCP63 is a member of a sprawling family of noncanonical CP genes, documenting a role of local gene family expansions in the emergence of the Tribolium compound eye CP repertoire. Comparisons with Drosophila and the mosquito Anopheles gambiae reveal a steady turnover of lens-enriched CP genes during insect evolution.

The apple 14-3-3 gene MdGRF6 negatively regulates salt tolerance

The 14-3-3 (GRF, general regulatory factor) regulatory proteins are highly conserved and are widely distributed throughout the eukaryotes. They are involved in the growth and development of organisms via target protein interactions. Although many plant 14-3-3 proteins were identified in response to stresses, little is known about their involvement in salt tolerance in apples. In our study, nineteen apple 14-3-3 proteins were cloned and identified. The transcript levels of Md14-3-3 genes were either up or down-regulated in response to salinity treatments. Specifically, the transcript level of MdGRF6 (a member of the Md14-3-3 genes family) decreased due to salt stress treatment. The phenotypes of transgenic tobacco lines and wild-type (WT) did not affect plant growth under normal conditions. However, the germination rate and salt tolerance of transgenic tobacco was lower compared to the WT. Transgenic tobacco demonstrated decreased salt tolerance. The transgenic apple calli overexpressing MdGRF6 exhibited greater sensitivity to salt stress compared to the WT plants, whereas the MdGRF6-RNAi transgenic apple calli improved salt stress tolerance. Moreover, the salt stress-related genes (MdSOS2, MdSOS3, MdNHX1, MdATK2/3, MdCBL-1, MdMYB46, MdWRKY30, and MdHB-7) were more strongly down-regulated in MdGRF6-OE transgenic apple calli lines than in the WT when subjected to salt stress treatment. Taken together, these results provide new insights into the roles of 14-3-3 protein MdGRF6 in modulating salt responses in plants.

Chromosome-scale assembly of the yellow mealworm genome

Background: The yellow mealworm beetle, Tenebrio molitor, is a promising alternative protein source for animal and human nutrition and its farming involves relatively low environmental costs. For these reasons, its industrial scale production started this century. However, to optimize and breed sustainable new T. molitor lines, the access to its genome remains essential. Methods: By combining Oxford Nanopore and Illumina Hi-C data, we constructed a high-quality chromosome-scale assembly of T. molitor. Then, we combined RNA-seq data and available coleoptera proteomes for gene prediction with GMOVE. Results: We produced a high-quality genome with a N50 = 21.9Mb with a completeness of 99.5% and predicted 21,435 genes with a median size of 1,780 bp. Gene orthology between T. molitor and Tribolium castaneum showed a highly conserved synteny between the two coleoptera and paralogs search revealed an expansion of histones in the T. molitor genome. Conclusions: The present genome will greatly help fundamental and applied research such as genetic breeding and will contribute to the sustainable production of the yellow mealworm.

Persistent Parental RNAi in the Beetle Tribolium castaneum Involves Maternal Transmission of Long Double-Stranded RNA

Parental RNA interference (pRNAi) is a powerful and widely used method for gene-specific knockdown. Yet in insects its efficacy varies between species, and how the systemic response is transmitted from mother to offspring remains elusive. Using the beetle Tribolium castaneum, an RT-qPCR strategy to distinguish the presence of double-stranded RNA (dsRNA) from endogenous mRNA is reported. It is found that injected dsRNA is directly transmitted into the egg and persists throughout embryogenesis. Despite this depletion of dsRNA from the mother, it is shown that strong pRNAi can persist for months before waning at strain-specific rates. In seeking the receptor proteins for cellular uptake of long dsRNA into the egg, a phylogenomics profiling approach of candidate proteins is also presented. A visualization strategy based on taxonomically hierarchical assessment of orthology clustering data to rapidly assess gene age and copy number changes, refined by sequence-based evidence, is demonstrated. Repeated losses of SID-1-like channel proteins in the arthropods, including wholesale loss in the Heteroptera (true bugs), which are nonetheless highly sensitive to pRNAi, are thereby documented. Overall, practical considerations for insect pRNAi against a backdrop of outstanding questions on the molecular mechanism of dsRNA transmission for long-term, systemic knockdown are elucidated.

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.

Phenotypic screen and transcriptomics approach complement each other in functional genomics of defensive stink gland physiology

Background

Functional genomics uses unbiased systematic genome-wide gene disruption or analyzes natural variations such as gene expression profiles of different tissues from multicellular organisms to link gene functions to particular phenotypes. Functional genomics approaches are of particular importance to identify large sets of genes that are specifically important for a particular biological process beyond known candidate genes, or when the process has not been studied with genetic methods before.

Results

Here, we present a large set of genes whose disruption interferes with the function of the odoriferous defensive stink glands of the red flour beetle Tribolium castaneum. This gene set is the result of a large-scale systematic phenotypic screen using RNA interference applied in a genome-wide forward genetics manner. In this first-pass screen, 130 genes were identified, of which 69 genes could be confirmed to cause phenotypic changes in the glands upon knock-down, which vary from necrotic tissue and irregular reservoir size to irregular color or separation of the secreted gland compounds. Gene ontology analysis revealed that many of those genes are encoding enzymes (peptidases and cytochromes P450) as well as proteins involved in membrane trafficking with an enrichment in lysosome and mineral absorption pathways. The knock-down of 13 genes caused specifically a strong reduction of para-benzoquinones in the gland reservoirs, suggesting a specific function in the synthesis of these toxic compounds. Only 14 of the 69 confirmed gland genes are differentially overexpressed in stink gland tissue and thus could have been detected in a transcriptome-based analysis. However, only one out of eight genes identified by a transcriptomics approach known to cause phenotypic changes of the glands upon knock-down was recognized by this phenotypic screen, indicating the limitation of such a non-redundant first-pass screen.

Conclusion

Our results indicate the importance of combining diverse and independent methodologies to identify genes necessary for the function of a certain biological tissue, as the different approaches do not deliver redundant results but rather complement each other. The presented phenotypic screen together with a transcriptomics approach are now providing a set of close to hundred genes important for odoriferous defensive stink gland physiology in beetles.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08822-z.

A high-quality genome of the dobsonfly Neoneuromus ignobilis reveals molecular convergences in aquatic insects

Neoneuromus ignobilis is an archaic holometabolous aquatic predatory insect. However, a lack of genomic resources hinders the use of whole genome sequencing to explore their genetic basis and molecular mechanisms for adaptive evolution. Here, we provided a high-contiguity, chromosome-level genome assembly of N. ignobilis using high coverage Nanopore and PacBio reads with the Hi-C technique. The final assembly is 480.67 MB in size, containing 12 telomere-ended pseudochromosomes with only 17 gaps. We compared 42 hexapod species genomes including six independent lineages comprising 11 aquatic insects, and found convergent expansions of long wavelength-sensitive and blue-sensitive opsins, thermal stress response TRP channels, and sulfotransferases in aquatic insects, which may be related to their aquatic adaptation. We also detected strong nonrandom signals of convergent amino acid substitutions in aquatic insects. Collectively, our comparative genomic analysis revealed the evidence of molecular convergences in aquatic insects during both gene family evolution and convergent amino acid substitutions.

The red flour beetle T. castaneum: elaborate genetic toolkit and unbiased large scale RNAi screening to study insect biology and evolution

The red flour beetle Tribolium castaneum has emerged as an important insect model system for a variety of topics. With respect to studying gene function, it is second only to the vinegar fly D. melanogaster. The RNAi response in T. castaneum is exceptionally strong and systemic, and it appears to target all cell types and processes. Uniquely for emerging model organisms, T. castaneum offers the opportunity of performing time- and cost-efficient large-scale RNAi screening, based on commercially available dsRNAs targeting all genes, which are simply injected into the body cavity. Well established transgenic and genome editing approaches are met by ease of husbandry and a relatively short generation time. Consequently, a number of transgenic tools like UAS/Gal4, Cre/Lox, imaging lines and enhancer trap lines are already available. T. castaneum has been a genetic experimental system for decades and now has become a workhorse for molecular and reverse genetics as well as in vivo imaging. Many aspects of development and general biology are more insect-typical in this beetle compared to D. melanogaster. Thus, studying beetle orthologs of well-described fly genes has allowed macro-evolutionary comparisons in developmental processes such as axis formation, body segmentation, and appendage, head and brain development. Transgenic approaches have opened new ways for in vivo imaging. Moreover, this emerging model system is the first choice for research on processes that are not represented in the fly, or are difficult to study there, e.g. extraembryonic tissues, cryptonephridial organs, stink gland function, or dsRNA-based pesticides.

A high-quality de novo genome assembly based on nanopore sequencing of a wild-caught coconut rhinoceros beetle (Oryctes rhinoceros)

BACKGROUND

An optimal starting point for relating genome function to organismal biology is a high-quality nuclear genome assembly, and long-read sequencing is revolutionizing the production of this genomic resource in insects. Despite this, nuclear genome assemblies have been under-represented for agricultural insect pests, particularly from the order Coleoptera. Here we present a de novo genome assembly and structural annotation for the coconut rhinoceros beetle, Oryctes rhinoceros (Coleoptera: Scarabaeidae), based on Oxford Nanopore Technologies (ONT) long-read data generated from a wild-caught female, as well as the assembly process that also led to the recovery of the complete circular genome assemblies of the beetle's mitochondrial genome and that of the biocontrol agent, Oryctes rhinoceros nudivirus (OrNV). As an invasive pest of palm trees, O. rhinoceros is undergoing an expansion in its range across the Pacific Islands, requiring new approaches to management that may include strategies facilitated by genome assembly and annotation.

RESULTS

High-quality DNA isolated from an adult female was used to create four ONT libraries that were sequenced using four MinION flow cells, producing a total of 27.2 Gb of high-quality long-read sequences. We employed an iterative assembly process and polishing with one lane of high-accuracy Illumina reads, obtaining a final size of the assembly of 377.36 Mb that had high contiguity (fragment N50 length = 12 Mb) and accuracy, as evidenced by the exceptionally high completeness of the benchmarked set of conserved single-copy orthologous genes (BUSCO completeness = 99.1%). These quality metrics place our assembly ahead of the published Coleopteran genomes, including that of an insect model, the red flour beetle (Tribolium castaneum). The structural annotation of the nuclear genome assembly contained a highly-accurate set of 16,371 protein-coding genes, with only 2.8% missing BUSCOs, and the expected number of non-coding RNAs. The number and structure of paralogous genes in a gene family like Sigma GST is lower than in another scarab beetle (Onthophagus taurus), but higher than in the red flour beetle (Tribolium castaneum), which suggests expansion of this GST class in Scarabaeidae. The quality of our gene models was also confirmed with the correct placement of O. rhinoceros among other members of the rhinoceros beetles (subfamily Dynastinae) in a phylogeny based on the sequences of 95 protein-coding genes in 373 beetle species from all major lineages of Coleoptera. Finally, we provide a list of 30 candidate dsRNA targets whose orthologs have been experimentally validated as highly effective targets for RNAi-based control of several beetles.

CONCLUSIONS

The genomic resources produced in this study form a foundation for further functional genetic research and management programs that may inform the control and surveillance of O. rhinoceros populations, and we demonstrate the efficacy of de novo genome assembly using long-read ONT data from a single field-caught insect.

Exploring new genomic territories with emerging model insects

Improvements in reference genome generation for insects and across the tree of life are extending the concept and utility of model organisms beyond traditional laboratory-tractable supermodels. Species or groups of species with comprehensive genome resources can be developed into model systems for studying a large variety of biological phenomena. Advances in sequencing and assembly technologies are supporting these emerging genome-enabled model systems by producing resources that are increasingly accurate and complete. Nevertheless, quality controls including assessing gene content completeness are required to ensure that these data can be included in expanding catalogues of high-quality references that will greatly advance understanding of insect biology and evolution.

Reference-Guided De Novo Genome Assembly of the Flour Beetle Tribolium freemani

The flour beetle Tribolium freemani is a sibling species of the model organism and important pest Tribolium castaneum. The two species are so closely related that they can produce hybrid progeny, but the genetic basis of their differences has not been revealed. In this work, we sequenced the T. freemani genome by applying PacBio HiFi technology. Using the well-assembled T. castaneum genome as a reference, we assembled 262 Mb of the T. freemani genomic sequence and anchored it in 10 linkage groups corresponding to nine autosomes and sex chromosome X. The assembly showed 99.8% completeness of conserved insect genes, indicating a high-quality reference genome. Comparison with the T. castaneum assembly revealed that the main differences in genomic sequence between the two sibling species come from repetitive DNA, including interspersed and tandem repeats. In this work, we also provided the complete assembled mitochondrial genome of T. freemani. Although the genome assembly needs to be ameliorated in tandemly repeated regions, the first version of the T. freemani reference genome and the complete mitogenome presented here represent useful resources for comparative evolutionary studies of related species and for further basic and applied research on different biological aspects of economically important pests.

Superoxide dismutase 6 is required during metamorphosis for the development of properly movable legs in Tribolium castaneum

The body form of holometabolous insects dramatically transforms from larval to adult stages during metamorphosis that occurs in the pupal stage. The larval disorganization and then new adult tissues are built up at this time. In motoneuron, larval neuronal cells degenerate, and new adult neurons are remodeled. Finally, adult neurons reconnect to new adult muscles. However, the factors that control metamorphosis have not yet been fully elucidated. Here, we show that an antioxidant enzyme, Tribolium castaneum superoxide dismutase 6 (TcSOD6), is secreted into the haemolymph and is required for proper movable legs during metamorphosis. TcSOD6 has a unique domain architecture and is mainly expressed in the pupal stage. The depletion of TcSOD6 expression in the pupa inhibits normal axon development and results in adults that display dysfunctional leg motions, suggesting that SOD6 expression is required for the development of properly movable legs. Therefore, we speculate that TcSOD6 might participate in some of the processes for larval neurons to be remodelled to new adult functions in the legs during metamorphosis, providing new insight into the evolution of SOD functions.

Screens in fly and beetle reveal vastly divergent gene sets required for developmental processes

Background

Most of the known genes required for developmental processes have been identified by genetic screens in a few well-studied model organisms, which have been considered representative of related species, and informative—to some degree—for human biology. The fruit fly Drosophila melanogaster is a prime model for insect genetics, and while conservation of many gene functions has been observed among bilaterian animals, a plethora of data show evolutionary divergence of gene function among more closely-related groups, such as within the insects. A quantification of conservation versus divergence of gene functions has been missing, without which it is unclear how representative data from model systems actually are.

Results

Here, we systematically compare the gene sets required for a number of homologous but divergent developmental processes between fly and beetle in order to quantify the difference of the gene sets. To that end, we expanded our RNAi screen in the red flour beetle Tribolium castaneum to cover more than half of the protein-coding genes. Then we compared the gene sets required for four different developmental processes between beetle and fly. We found that around 50% of the gene functions were identified in the screens of both species while for the rest, phenotypes were revealed only in fly (~ 10%) or beetle (~ 40%) reflecting both technical and biological differences. Accordingly, we were able to annotate novel developmental GO terms for 96 genes studied in this work. With this work, we publish the final dataset for the pupal injection screen of the iBeetle screen reaching a coverage of 87% (13,020 genes).

Conclusions

We conclude that the gene sets required for a homologous process diverge more than widely believed. Hence, the insights gained in flies may be less representative for insects or protostomes than previously thought, and work in complementary model systems is required to gain a comprehensive picture. The RNAi screening resources developed in this project, the expanding transgenic toolkit, and our large-scale functional data make T. castaneum an excellent model system in that endeavor.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01231-4.

Deeply conserved synteny and the evolution of metazoan chromosomes

Animal genomes show networks of deeply conserved gene linkages whose phylogenetic scope and chromosomal context remain unclear. Here, we report chromosome-scale conservation of synteny among bilaterians, cnidarians, and sponges and use comparative analysis to reconstruct ancestral chromosomes across major animal groups. Comparisons among diverse metazoans reveal the processes of chromosome evolution that produced contemporary karyotypes from their Precambrian progenitors. On the basis of these findings, we introduce a simple algebraic representation of chromosomal change and use it to establish a unified systematic framework for metazoan chromosome evolution. We find that fusion-with-mixing, a previously unappreciated mode of chromosome change, has played a central role. We find that relicts of several metazoan chromosomal units are preserved in unicellular eukaryotes. These conserved pre-metazoan linkages include the chromosomal unit that encodes the most diverse set of metazoan homeobox genes, suggesting a candidate genomic context for the early diversification of this key gene family.

Tribolium castaneum: A Model Insect for Fundamental and Applied Research

Tribolium castaneum has a long history as a model species in many distinct subject areas, but improved connections among the genetics, genomics, behavioral, ecological, and pest management fields are needed to fully realize this species' potential as a model. Tribolium castaneum was the first beetle whose genome was sequenced, and a new genome assembly and enhanced annotation, combined with readily available genomic research tools, have facilitated its increased use in a wide range of functional genomics research. Research into T. castaneum's sensory systems, response to pheromones and kairomones, and patterns of movement and landscape utilization has improved our understanding of behavioral and ecological processes. Tribolium castaneum has also been a model in the development of pest monitoring and management tactics, including evaluation of insecticide resistance mechanisms. Application of functional genomics approaches to behavioral, ecological, and pest management research is in its infancy but offers a powerful tool that can link mechanism with function and facilitate exploitation of these relationships to better manage this important food pest.

BUSCO: Assessing Genomic Data Quality and Beyond

Evaluation of the quality of genomic "data products" such as genome assemblies or gene sets is of critical importance in order to recognize possible issues and correct them during the generation of new data. It is equally essential to guide subsequent or comparative analyses with existing data, as the correct interpretation of the results necessarily requires knowledge about the quality level and reliability of the inputs. Using datasets of near universal single-copy orthologs derived from OrthoDB, BUSCO can estimate the completeness and redundancy of genomic data by providing biologically meaningful metrics based on expected gene content. These can complement technical metrics such as contiguity measures (e.g., number of contigs/scaffolds, and N50 values). Here, we describe the use of the BUSCO tool suite to assess different data types that can range from genome assemblies of single isolates and assembled transcriptomes and annotated gene sets to metagenome-assembled genomes where the taxonomic origin of the species is unknown. BUSCO is the only tool capable of assessing all these types of sequences from both eukaryotic and prokaryotic species. The protocols detail the various BUSCO running modes and the novel workflows introduced in versions 4 and 5, including the batch analysis on multiple inputs, the auto-lineage workflow to run assessments without specifying a dataset, and a workflow for the evaluation of (large) eukaryotic genomes. The protocols further cover the BUSCO setup, guidelines to interpret the results, and BUSCO "plugin" workflows for performing common operations in genomics using BUSCO results, such as building phylogenomic trees and visualizing syntenies. © 2021 The Authors. Current Protocols published by Wiley Periodicals LLC. Basic Protocol 1: Assessing an input sequence with a BUSCO dataset specified manually Basic Protocol 2: Assessing an input sequence with a dataset automatically selected by BUSCO Basic Protocol 3: Assessing multiple inputs Alternate Protocol: Decreasing analysis runtime when assessing a large number of small genomes with BUSCO auto-lineage workflow and Snakemake Support Protocol 1: BUSCO setup Support Protocol 2: Visualizing BUSCO results Support Protocol 3: Building phylogenomic trees.

The T-box gene optomotor-blind organizes proximodistal leg patterning in the beetle Tribolium castaneum by repressing dorsal Dpp pathway activity

Leg axis formation in Drosophila is organized by Wingless (Wg) and Decapentaplegic (Dpp) that control a number of downstream factors to pattern the dorsoventral (DV) and proximodistal (PD) axis. The T-box genes are important downstream factors mainly involved in dorsoventral leg axis formation. The ventral side is specified by H15 and midline, whereas optomotor-blind (omb) and Dorsocross (Doc1) are factors to specify dorsal cell fates. We show here that omb also organizes PD leg axis patterning in the beetle Tribolium castaneum. In the legs, Tc-omb is expressed along the dorsal side and represses ventral factors like wg and H15. Intriguingly, removing Tc-omb function leads to the activation of the Dpp pathway along the dorsal side of the legs, thus mimicking normal dpp expression in Drosophila. Dpp activity along the dorsal side leads to altered expression of proximal-distal patterning genes such as Distal-less (Dll) and dachshund (dac). Our results indicate a cell-autonomous activation of Dll and repression of dac by dpp. These findings are compatible with the cross-regulatory "cascade model" of proximal-distal leg imaginal disc patterning of Drosophila.

The transposable element-rich genome of the cereal pest Sitophilus oryzae

BACKGROUND

The rice weevil Sitophilus oryzae is one of the most important agricultural pests, causing extensive damage to cereal in fields and to stored grains. S. oryzae has an intracellular symbiotic relationship (endosymbiosis) with the Gram-negative bacterium Sodalis pierantonius and is a valuable model to decipher host-symbiont molecular interactions.

RESULTS

We sequenced the Sitophilus oryzae genome using a combination of short and long reads to produce the best assembly for a Curculionidae species to date. We show that S. oryzae has undergone successive bursts of transposable element (TE) amplification, representing 72% of the genome. In addition, we show that many TE families are transcriptionally active, and changes in their expression are associated with insect endosymbiotic state. S. oryzae has undergone a high gene expansion rate, when compared to other beetles. Reconstruction of host-symbiont metabolic networks revealed that, despite its recent association with cereal weevils (30 kyear), S. pierantonius relies on the host for several amino acids and nucleotides to survive and to produce vitamins and essential amino acids required for insect development and cuticle biosynthesis.

CONCLUSIONS

Here we present the genome of an agricultural pest beetle, which may act as a foundation for pest control. In addition, S. oryzae may be a useful model for endosymbiosis, and studying TE evolution and regulation, along with the impact of TEs on eukaryotic genomes.

Deep Conservation of Hid-Like RHG Gene Family Homologs in Winged Insects Revealed by "Taxon Hopping" BLAST

Together with sickle (skl), the Drosophila paralogs reaper (rpr), head involution defective (hid), and grim (RHG) control a critical switch in the induction of programmed cell death. RHG homologs have been identified in other dipteran and lepidopteran species but not beyond. Revisiting this issue with a "taxon hopping" BLAST search strategy in current genome and transcriptome resources, I detected high confidence RHG homologs in Coleoptera, Hymenoptera, Hemiptera, and Dictyoptera. Analyses of gene structure and protein sequence conservation revealed aconserved splicing pattern and highly conserved amino acid residues at both the N- and C-terminal ends that identify hid as the most ancestrally organized RHG gene family member in Drosophila. hid-like RHG homologs were also detected in mosquitoes, redefining their michelob_x (mx) genes as an expansion of derived RHG homologs. Only singleton homologs were detected in the large majority of other insect clades. Lepidopteran RHG homologs, however, stand out by producing an evolutionarily-derived splice isoform, identified in previous work, in addition to the newly detected hid-like isoform. Exceptional sequence diversification of select RHG homologs at the family- and genus-level explain their previous elusiveness in important insect genome model species like the red flour beetle Tribolium castaneum and the pea aphid Acyrthosiphon pisum. Combined, these findings expand the minimal age of the RHG gene family by about 100 million years and open new avenues for molecular cell death studies in insects.

Mutually Exclusive Expression of Closely Related Odorant-Binding Proteins 9A and 9B in the Antenna of the Red Flour Beetle Tribolium castaneum

Olfaction is crucial for insects to find food sources, mates, and oviposition sites. One of the initial steps in olfaction is facilitated by odorant-binding proteins (OBPs) that translocate hydrophobic odorants through the aqueous olfactory sensilla lymph to the odorant receptor complexes embedded in the dendritic membrane of olfactory sensory neurons. The Tribolium castaneum (Coleoptera, Tenebrionidae) OBPs encoded by the gene pair TcasOBP9A and TcasOBP9B represent the closest homologs to the well-studied Drosophila melanogaster OBP Lush (DmelOBP76a), which mediates pheromone reception. By an electroantennographic analysis, we can show that these two OBPs are not pheromone-specific but rather enhance the detection of a broad spectrum of organic volatiles. Both OBPs are expressed in the antenna but in a mutually exclusive pattern, despite their homology and gene pair character by chromosomal location. A phylogenetic analysis indicates that this gene pair arose at the base of the Cucujiformia, which dates the gene duplication event to about 200 Mio years ago. Therefore, this gene pair is not the result of a recent gene duplication event and the high sequence conservation in spite of their expression in different sensilla is potentially the result of a common function as co-OBPs.

Whole-Genome Sequence of the Trypoxylus dichotomus Japanese rhinoceros beetle

The draft whole-genome sequence of the Japanese rhinoceros beetle, Trypoxylus dichotomus was obtained using long-read PacBio sequence technology. The final assembled genome consisted of 739 Mbp in 2,347 contigs, with 24.5× mean coverage and a G+C content of 35.99%.

Shaking hands is a homeodomain transcription factor that controls axon outgrowth of central complex neurons in the insect model Tribolium

Gene regulatory mechanisms that specify subtype identity of central complex (CX) neurons are the subject of intense investigation. The CX is a compartment within the brain common to all insect species and functions as a ‘command center’ that directs motor actions. It is made up of several thousand neurons, with more than 60 morphologically distinct identities. Accordingly, transcriptional programs must effect the specification of at least as many neuronal subtypes. We demonstrate a role for the transcription factor Shaking hands (Skh) in the specification of embryonic CX neurons in Tribolium. The developmental dynamics of skh expression are characteristic of terminal selectors of subtype identity. In the embryonic brain, skh expression is restricted to a subset of neurons, many of which survive to adulthood and contribute to the mature CX. skh expression is maintained throughout the lifetime in at least some CX neurons. skh knockdown results in axon outgrowth defects, thus preventing the formation of an embryonic CX primordium. The previously unstudied Drosophila skh shows a similar embryonic expression pattern, suggesting that subtype specification of CX neurons may be conserved.

Summary: A detailed examination of the developmental expression of the homeodomain transcription factor Shaking hands in Tribolium reveals a role in the formation of the central complex primordium.

The β2Tubulin, Rad50-ATPase and enolase cis-regulatory regions mediate male germline expression in Tribolium castaneum

Genetics-based pest management processes, including the sterile insect technique, are an effective method for the control of some pest insects. However, current SIT methods are not directly transferable to many important pest insect species due to the lack of genetic sexing strains. Genome editing is revolutionizing the way we conduct genetics in insects, including in Tribolium castaneum, an important genetic model and agricultural pest. We identified orthologues of β₂Tubulin, Rad50-ATPase and enolase in T. castaneum. Using RT-PCR, we confirmed that these genes are predominantly expressed in the testis. PiggyBac-based transformation of T. castaneum cis-regulatory regions derived from Tc-β₂t, Tc-rad50 or Tc-eno resulted in EGFP expression specifically in the T. castaneum testis. Additionally, we determined that each of these regulatory regions regulates EGFP expression in different cell types of the male gonad. Cis-regulatory regions from Tc-β₂t produced EGFP expression throughout spermatogenesis and also in mature sperms; Tc-rad50 resulted in expression only in the haploid spermatid, while Tc-eno expressed EGFP in late spermatogenesis. In summary, the regulatory cis-regions characterized in this study are not only suited to study male gonadal function but could be used for development of transgenic sexing strains that produce one sex in pest control strategies.

A highly-contiguous genome assembly of the Eurasian spruce bark beetle, Ips typographus, provides insight into a major forest pest

Conifer-feeding bark beetles are important herbivores and decomposers in forest ecosystems. These species complete their life cycle in nutritionally poor substrates and some can kill enormous numbers of trees during population outbreaks. The Eurasian spruce bark beetle (Ips typographus) can destroy >100 million m3 of spruce in a single year. We report a 236.8 Mb I. typographus genome assembly using PacBio long-read sequencing. The final phased assembly has a contig N50 of 6.65 Mb in 272 contigs and is predicted to contain 23,923 protein-coding genes. We reveal expanded gene families associated with plant cell wall degradation, including pectinases, aspartyl proteases, and glycosyl hydrolases. This genome sequence from the genus Ips provides timely resources to address questions about the evolutionary biology of the true weevils (Curculionidae), one of the most species-rich animal families. In forests of today, increasingly stressed by global warming, this draft genome may assist in developing pest control strategies to mitigate outbreaks.