Molecular traces of alternative social organization in a termite genome

Genome-Wide Profiling of P450 Gene Expression Reveals Caste-Specific and Developmental Patterns in Solenopsis invicta

P450 enzymes are integral to insect physiology, metabolism, hormone regulation, and adaptation to environmental challenges. By leveraging transcriptomic and genomic data, this study characterized the expression of 68 unique P450 genes across developmental stages and castes in the red imported fire ant (Solenopsis invicta), uncovering stage- and caste-specific differential expression patterns. Genes from the CYP4, CYP6, and CYP9 families, known for metabolizing exogenous and endogenous compounds, were highly expressed in early larval stages and minim workers, underscoring their roles in supporting rapid growth, hormone metabolism, colony maintenance, and brood care. The overexpression of CYP4AA1-linked to pheromone production-in queens, female alates, and female alate pupae highlights its critical functions in reproductive dominance, social structure maintenance, and colony dynamics. Here, juvenile hormone biosynthesis genes, including CYP305A1 and CYP315A1, exhibited significant overexpression in later instar larvae and larger workers, emphasizing their roles in development and in fulfilling colony-wide physiological demands. The "Halloween genes" (CYP302A1, CYP306A1, CYP315A1, CYP307A1, and CYP314A1) and CYP18A1 demonstrated dynamic regulation across developmental stages and castes, reflecting their essential contributions to hormonal production and balance throughout S. invicta's lifecycle. These findings offer valuable insights into the molecular and biological mechanisms driving S. invicta's social organization, developmental transitions, physiological adaptations, and evolutionary success. They also provide a foundation for future research into the regulatory pathways governing P450 gene expression and function.

Comparison Between Worker and Soldier Transcriptomes of Termite Neotermes binovatus Reveals Caste Specialization of Host-Flagellate Symbiotic System

Termites are eusocial insects with functionally specialized workers and soldiers, both sharing the same genotype. Additionally, lower termites host flagellates in their hindguts that assist in wood digestion. However, worker-biased and soldier-biased gene expression patterns of the host-flagellate symbiotic system remain underexplored in most taxonomic groups. In this study, we sequenced high-depth transcriptomes from the workers and soldiers of a lower termite, Neotermes binovatus (Kalotermitidae), to investigate the differentially expressed termite transcripts, flagellate transcript abundance, and co-expression patterns of the host-flagellate transcript pairs in both castes. The worker-biased transcripts were enriched in functions related to cuticle development, nervous system regulation, pheromone biosynthesis, and metabolism, whereas the soldier-biased transcripts were predominantly involved in muscle development and kinesis, body morphogenesis, protein modification, and aggression. Flagellate transcripts from the orders Cristamonadida, Trichomonadida, Tritrichomonadida, and Oxymonadida were identified in both workers and soldiers, with the abundance of most flagellate transcripts tending to be higher in workers than in soldiers. Furthermore, we observed a much larger number of strong co-expression correlations between the termite and flagellate transcripts in workers than in soldiers, suggesting the possibility that soldiers depend more on food processed by worker holobionts than on their own symbiotic system. This research provides insights into the functional specialization of the host-flagellate symbiotic system in the worker and soldier castes of termites, supporting the workers' roles in nest maintenance, preliminary food processing, and communication, while emphasizing the defensive role of soldiers. Additionally, it offers new perspectives on the potential termite-flagellate interactions and underscores the need for whole-genome data of termite flagellates in further studies.

Genome Report: Improved chromosome-level genome assembly of the American cockroach, Periplaneta americana

The American cockroach, Periplaneta americana, is a cosmopolitan insect notorious for thriving among humans undeterred by attempts to eliminate it. The traits that contribute to its ubiquity as an opportunistic pest, such as long lifespan, expansive neurosensory capacity, and nutritional flexibility, also make P. americana an excellent invertebrate model organism with a long history in neuroscience and physiological research. Current genetic resources available for P. americana highlight its large, complex genome and richly diverse transcriptional capabilities, but fall short of producing a complete, chromosome-level genome. Here, we present a high-quality de novo genome assembly of a laboratory-raised adult female P. americana using a combination of high fidelity PacBio long reads and Hi-C sequencing. The final 3.23 Gb genome was assembled with chromosomal resolution into 17 scaffolds, consistent with previous karyotype analysis, and has a scaffold N50 of 188.1 Mb and genome BUSCO score of 99.7%. This assembly includes a chromosome that was missing from the previous reference genome for this species. Protein prediction and annotation were performed via the NCBI Eukaryotic Genome Annotation Pipeline, which identified 16,780 protein-coding genes and generated an annotation BUSCO score of 97.8%. Ortholog comparisons with available Blattodea assemblies highlight the expanded chemosensory and immune capabilities of P. americana compared to termite relatives. This genome assembly is a valuable tool for facilitating future research on the biology and evolution of this remarkable insect.

Identification of an additional periplanone receptor family gene preferentially expressed in the male antennae of the American cockroach

Periplaneta cockroaches use periplanone analogs as female sex pheromones to attract males. We previously identified two periplanone receptor genes, PameOR1 and PameOR2, in the American cockroach Periplaneta americana. Here, we report the identification of PameOR1-like, an additional olfactory receptor resembling PameOR1 in P. americana. PameOR1-like showed high-level sequence similarity to PameOR1, and is preferentially expressed in the male antennae. Quantitative expression analysis revealed that, in the adult male antennae, the expression level of PameOR1 is 2.4-fold higher than that of PameOR1-like. Fluorescent in situ hybridization revealed that PameOR1 and PameOR1-like are co-expressed in the periplanone-A-responsive sensory neurons within the single walled-B sensilla on the male adult antennae. These data support the idea that PameOR1 and PameOR1-like are generated by a recent gene duplication event and play a redundant function in sex pheromone reception in P. americana.

Lipophorin receptor knockdown reduces hatchability of kissing bug Rhodnius prolixus eggs

Lipophorin is the primary lipoprotein present in the hemolymph of insects, responsible for the lipids' transport between organs. It interacts with specific sites on cell membranes in an essential process for transferring lipids. The lipophorin receptor is the protein responsible for the interaction between lipophorin and cell membranes. In the kissing bug Rhodnius prolixus, much information on the interaction of lipophorin with organs is available. However, molecular data on the lipophorin receptor and its functions is still needed. Here, we explored lipophorin receptor gene expression and functions using a functional genomics approach. The R. prolixus genome encodes seven genes from the low-density lipoprotein receptor family, including a single ortholog of the lipophorin receptor. All organs analyzed (anterior and posterior midguts, fat body, ovaries, and flight muscle) expressed this gene. In the fat body, blood-feeding strongly reduced lipophorin receptor gene expression. Lipophorin receptor knockdown by RNA interference delayed egg laying and reduced the triacylglycerol in laid eggs without altering lipid stores in the fat body or lipid levels in the hemolymph. In the ovaries, lipophorin receptor knockdown reduces the expression of acetyl-CoA carboxylase and a fatty acid synthase while altered the gene expression profile in the fat body, causing an increase in the expression of carnitine palmitoyltransferase 1 and a reduction in Brummer lipase and vitellogenin 2. RNA interference treatment reduced the hatching of the eggs, causing the collapse and darkening of the laid eggs, in addition to the hatching of deformed first-stage nymphs. Furthermore, the structure of the chorion showed distortions in patterns and cracks and reduced hydrocarbon levels. These results show that the lipophorin receptor alone is not essential for lipid physiology in R. prolixus. However, this protein plays a fundamental role in the viability of eggs and, consequently, in insect reproduction.

Evidence for a Novel X Chromosome in Termites

Termites, together with cockroaches, belong to the Blattodea. They possess an XX/XY sex determination system which has evolved from an XX/X0 system present in other Blattodean species, such as cockroaches and wood roaches. Little is currently known about the sex chromosomes of termites, their gene content, or their evolution. We here investigate the X chromosome of multiple termite species and compare them with the X chromosome of cockroaches using genomic and transcriptomic data. We find that the X chromosome of the termite Macrotermes natalensis is large and differentiated showing hall marks of sex chromosome evolution such as dosage compensation, while this does not seem to be the case in the other two termite species investigated here where sex chromosomes may be evolutionary younger. Furthermore, the X chromosome in M. natalensis is different from the X chromosome found in the cockroach Blattella germanica indicating that sex chromosome turn-over events may have happened during termite evolution.

Termites and subsocial roaches inherited many bacterial-borne carbohydrate-active enzymes (CAZymes) from their common ancestor

Termites digest wood using Carbohydrate-Active Enzymes (CAZymes) produced by gut bacteria with whom they have cospeciated at geological timescales. Whether CAZymes were encoded in the genomes of their ancestor's gut bacteria and transmitted to modern termites or acquired more recently from bacteria not associated with termites is unclear. We used gut metagenomes from 195 termites and one Cryptocercus, the sister group of termites, to investigate the evolution of termite gut bacterial CAZymes. We found 420 termite-specific clusters in 81 bacterial CAZyme gene trees, including 404 clusters showing strong cophylogenetic patterns with termites. Of the 420 clusters, 131 included at least one bacterial CAZyme sequence associated with Cryptocercus or Mastotermes, the sister group of all other termites. Our results suggest many bacterial CAZymes have been encoded in the genomes of termite gut bacteria since termite origin, indicating termites rely upon many bacterial CAZymes endemic to their guts to digest wood.

Genomic and transcriptomic analyses of a social hemipteran provide new insights into insect sociality

The origin of sociality represents one of the most important evolutionary transitions. Insect sociality evolved in some hemipteran aphids, which can produce soldiers and normal nymphs with distinct morphology and behaviour through parthenogenesis. The lack of genomic data resources has hindered the investigations into molecular mechanisms underlying their social evolution. Herein, we generated the first chromosomal-level genome of a social hemipteran (Pseudoregma bambucicola) with highly specialized soldiers and performed comparative genomic and transcriptomic analyses to elucidate the molecular signatures and regulatory mechanisms of caste differentiation. P. bambucicola has a larger known aphid genome of 582.2 Mb with an N50 length of 11.24 Mb, and about 99.6% of the assembly was anchored to six chromosomes with a scaffold N50 of 98.27 Mb. A total of 14,027 protein-coding genes were predicted and 37.33% of the assembly were identified as repeat sequences. The social evolution is accompanied by a variety of changes in genome organization, including expansion of gene families related to transcription factors, transposable elements, as well as species-specific expansions of certain sugar transporters and UGPases involved in carbohydrate metabolism. We also characterized large candidate gene sets linked to caste differentiation and found evidence of expression regulation and positive selection acting on energy metabolism and muscle structure, explaining the soldier-specific traits including morphological and behavioural specialization, developmental arrest and infertility. Overall, this study offers new insights into the molecular basis of social aphids and the evolution of insect sociality and also provides valuable data resources for further comparative and functional studies.

Orthoptera-specific target enrichment (OR-TE) probes resolve relationships over broad phylogenetic scales

Phylogenomic data are revolutionizing the field of insect phylogenetics. One of the most tenable and cost-effective methods of generating phylogenomic data is target enrichment, which has resulted in novel phylogenetic hypotheses and revealed new insights into insect evolution. Orthoptera is the most diverse insect order within polyneoptera and includes many evolutionarily and ecologically interesting species. Still, the order as a whole has lagged behind other major insect orders in terms of transitioning to phylogenomics. In this study, we developed an Orthoptera-specific target enrichment (OR-TE) probe set from 80 transcriptomes across Orthoptera. The probe set targets 1828 loci from genes exhibiting a wide range of evolutionary rates. The utility of this new probe set was validated by generating phylogenomic data from 36 orthopteran species that had not previously been subjected to phylogenomic studies. The OR-TE probe set captured an average of 1037 loci across the tested taxa, resolving relationships across broad phylogenetic scales. Our detailed documentation of the probe design and bioinformatics process is intended to facilitate the widespread adoption of this tool.

Colony identity clues for Syntermes grandis (Blattodea: Termitidae) individuals using near-infrared spectroscopy and PLS-DA approach

Termites are social insects with high species diversity in tropical ecosystems. Multivariate analysis with near-infrared spectroscopy (NIRS) and data interpretation can separate social insects belonging to different colonies of the same species. The objective of this study was to propose the use of discriminant analysis by partial least squares (PLS-DA) combined with NIRS to identify the colonial origin of the Syntermes grandis (Rambur, 1842) (Blattodea: Termitidae) in 2 castes. Six ground S. grandis colonies were identified and mapped; 30 workers and 30 soldier termites in each colony were submitted to spectral measurement with NIRS. PLS-DA applied to the termites' spectral absorbance was used to detect a spectral pattern per S. grandis colony by caste. PLS-DA regression with NIRS proved to be an approach with 99.9% accuracy for identifying the colonial origin of S. grandis workers and 98.3% for soldiers. The methodology showed the importance of qualitatively characterizing the colonial phenotypic response of this species. NIRS is a high-precision approach to identifying the colony origin of S. grandis workers and soldiers. The PLS-DA can be used to design ecological field studies to identify colony territorial competition and foraging behavior of subterranean termite species.

Genomic data provide insights into the classification of extant termites

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.

Evolution of the sex-determination gene Doublesex within the termite lineage

The molecular mechanism of sex determination has long been considered conserved in insects. However, recent studies of hemimetabolous insects have challenged this notion. One notable example is termites. In Reticulitermes speratus, a homolog of sex determination gene, Doublesex (RsDsx), exhibits characteristics that are distinct from those of other insects, including sister-group cockroaches. It comprises a single exon, contains only doublesex/mab-3 DNA-binding domain (DM) but lacks a conserved oligomerization domain (OD), and exhibits transcriptional activity only in males. To investigate whether these characteristics are widespread within the termite lineage, we identified Dsx homologs in three different families. The absence of the conserved OD sequences was observed in all termite species examined, whereas the number of exons and expression patterns between sexes varied among families. Particularly, distinctive differences in Dsx were found in species from the Archotermopsidae and Kalotermitidae, both of which have a linear caste developmental pathway. Our findings indicate that diversification of Dsx structure and expression patterns may have contributed to ecological diversification, such as caste developmental pathways, within the termite lineage.

Major changes in domain arrangements are associated with the evolution of termites

Domains as functional protein units and their rearrangements along the phylogeny can shed light on the functional changes of proteomes associated with the evolution of complex traits like eusociality. This complex trait is associated with sterile soldiers and workers, and long-lived, highly fecund reproductives. Unlike in Hymenoptera (ants, bees, and wasps), the evolution of eusociality within Blattodea, where termites evolved from within cockroaches, was accompanied by a reduction in proteome size, raising the question of whether functional novelty was achieved with existing rather than novel proteins. To address this, we investigated the role of domain rearrangements during the evolution of termite eusociality. Analysing domain rearrangements in the proteomes of three solitary cockroaches and five eusocial termites, we inferred more than 5,000 rearrangements over the phylogeny of Blattodea. The 90 novel domain arrangements that emerged at the origin of termites were enriched for several functions related to longevity, such as protein homeostasis, DNA repair, mitochondrial activity, and nutrient sensing. Many domain rearrangements were related to changes in developmental pathways, important for the emergence of novel castes. Along with the elaboration of social complexity, including permanently sterile workers and larger, foraging colonies, we found 110 further domain arrangements with functions related to protein glycosylation and ion transport. We found an enrichment of caste-biased expression and splicing within rearranged genes, highlighting their importance for the evolution of castes. Furthermore, we found increased levels of DNA methylation among rearranged compared to non-rearranged genes suggesting fundamental differences in their regulation. Our findings indicate the importance of domain rearrangements in the generation of functional novelty necessary for termite eusociality to evolve.

Genome reduction and horizontal gene transfer in the evolution of Endomicrobia-rise and fall of an intracellular symbiosis with termite gut flagellates

Bacterial endosymbionts of eukaryotic hosts typically experience massive genome reduction, but the underlying evolutionary processes are often obscured by the lack of free-living relatives. Endomicrobia, a family-level lineage of host-associated bacteria in the phylum Elusimicrobiota that comprises both free-living representatives and endosymbionts of termite gut flagellates, are an excellent model to study evolution of intracellular symbionts. We reconstructed 67 metagenome-assembled genomes (MAGs) of Endomicrobiaceae among more than 1,700 MAGs from the gut microbiota of a wide range of termites. Phylogenomic analysis confirmed a sister position of representatives from termites and ruminants, and allowed to propose eight new genera in the radiation of Endomicrobiaceae. Comparative genome analysis documented progressive genome erosion in the new genus Endomicrobiellum, which comprises all flagellate endosymbionts characterized to date. Massive gene losses were accompanied by the acquisition of new functions by horizontal gene transfer, which led to a shift from a glucose-based energy metabolism to one based on sugar phosphates. The breakdown of glycolysis and many anabolic pathways for amino acids and cofactors in several subgroups was compensated by the independent acquisition of new uptake systems, including an ATP/ADP antiporter, from other gut microbiota. The putative donors are mostly flagellate endosymbionts from other bacterial phyla, including several, hitherto unknown lineages of uncultured Alphaproteobacteria, documenting the importance of horizontal gene transfer in the convergent evolution of these intracellular symbioses. The loss of almost all biosynthetic capacities in some lineages of Endomicrobiellum suggests that their originally mutualistic relationship with flagellates is on its decline.IMPORTANCEUnicellular eukaryotes are frequently colonized by bacterial and archaeal symbionts. A prominent example are the cellulolytic gut flagellates of termites, which harbor diverse but host-specific bacterial symbionts that occur exclusively in termite guts. One of these lineages, the so-called Endomicrobia, comprises both free-living and endosymbiotic representatives, which offers the unique opportunity to study the evolutionary processes underpinning the transition from a free-living to an intracellular lifestyle. Our results revealed a progressive gene loss in energy metabolism and biosynthetic pathways, compensated by the acquisition of new functions via horizontal gene transfer from other gut bacteria, and suggest the eventual breakdown of an initially mutualistic symbiosis. Evidence for convergent evolution of unrelated endosymbionts reflects adaptations to the intracellular environment of termite gut flagellates.

Body part-specific development in termite caste differentiation: crosstalk between hormonal actions and developmental toolkit genes

In social insects, interactions among colony members trigger caste differentiation with morphological modifications. During caste differentiation in termites, body parts and caste-specific morphologies are modified during postembryonic development under endocrine controls such as juvenile hormone (JH) and ecdysone. In addition to endocrine factors, developmental toolkit genes such as Hox- and appendage-patterning genes also contribute to the caste-specific body part modifications. These toolkits are thought to provide spatial information for specific morphogenesis. During social evolution, the complex crosstalks between physiological and developmental mechanisms should be established, leading to the sophisticated caste systems. This article reviews recent studies on these mechanisms underlying the termite caste differentiation and addresses implications for the evolution of caste systems in termites.

The Function of Termicin from Odontotermes formosanus (Shiraki) in the Defense against Bacillus thuringiensis (Bt) and Beauveria bassiana (Bb) Infection

Odontotermes formosanus (Shiraki) is a subterranean termite species known for causing severe damage to trees and structures such as dams. During the synergistic evolution of O. formosanus with pathogenic bacteria, the termite has developed a robust innate immunity. Termicin is a crucial antimicrobial peptide in termites, significantly contributing to the defense against external infections. Building upon the successful construction and expression of the dsRNA-HT115 engineering strains of dsOftermicin1 and dsOftermicin2 in our laboratory, this work employs the ultrasonic breaking method to establish an inactivated dsOftermicins-HT115 technological system capable of producing a substantial quantity of dsRNA. This approach also addresses the limitation of transgenic strains which cannot be directly applied. Treatment of O. formosanus with dsOftermicins produced by this method could enhance the virulence of both Bt and Bb to the termites. This study laid the theoretical groundwork for the development of novel termite immunosuppressants and for the advancement and application of termite biological control strategies.

Pervasive relaxed selection in termite genomes

Genetic changes that enabled the evolution of eusociality have long captivated biologists. More recently, attention has focussed on the consequences of eusociality on genome evolution. Studies have reported higher molecular evolutionary rates in eusocial hymenopteran insects compared with their solitary relatives. To investigate the genomic consequences of eusociality in termites, we analysed nine genomes, including newly sequenced genomes from three non-eusocial cockroaches. Using a phylogenomic approach, we found that termite genomes have experienced lower rates of synonymous substitutions than those of cockroaches, possibly as a result of longer generation times. We identified higher rates of non-synonymous substitutions in termite genomes than in cockroach genomes, and identified pervasive relaxed selection in the former (24-31% of the genes analysed) compared with the latter (2-4%). We infer that this is due to reductions in effective population size, rather than gene-specific effects (e.g. indirect selection of caste-biased genes). We found no obvious signature of increased genetic load in termites, and postulate efficient purging of deleterious alleles at the colony level. Additionally, we identified genomic adaptations that may underpin caste differentiation, such as genes involved in post-translational modifications. Our results provide insights into the evolution of termites and the genomic consequences of eusociality more broadly.

Untangling the gordian knot: The intertwining interactions between developmental hormone signaling and epigenetic mechanisms in insects

Hormones control developmental and physiological processes, often by regulating the expression of multiple genes simultaneously or sequentially. Crosstalk between hormones and epigenetics is pivotal to dynamically coordinate this process. Hormonal signals can guide the addition and removal of epigenetic marks, steering gene expression. Conversely, DNA methylation, histone modifications and non-coding RNAs can modulate regional chromatin structure and accessibility and regulate the expression of numerous (hormone-related) genes. Here, we provide a review of the interplay between the classical insect hormones, ecdysteroids and juvenile hormones, and epigenetics. We summarize the mode-of-action and roles of these hormones in post-embryonic development, and provide a general overview of epigenetic mechanisms. We then highlight recent advances on the interactions between these hormonal pathways and epigenetics, and their involvement in development. Furthermore, we give an overview of several 'omics techniques employed in the field. Finally, we discuss which questions remain unanswered and possible avenues for future research.

Oleic acid triggers burial behavior in a termite population through an odorant binding protein

Social insects maintain hygienic conditions through their social immunity behaviors. Among these behaviors, burial behavior of termites is central for protecting healthy individuals from corpses. Many factors trigger burial behavior, and it is generally believed that chemicals released by corpses, such as oleic acid, are the most important cues for triggering burial behavior in termites. However, the contribution of the olfactory system to this behavior remains unclear. Here we report an odorant binding protein (OBP) that transports oleic acid and triggers burial behavior in Coptotermes formosanus Shiraki. We demonstrated that CforOBP7 is highly expressed in the antennae of workers. Fluorescent competition binding experiments exhibited that CforOBP7 has a strong affinity for oleic acid. Furthermore, the antennal response to oleic acid was significantly reduced, and oleic acid-triggered burial behavior was also inhibited in CforOBP7-silenced termites. We conclude that CforOBP7 governs the burial behavior of C. formosanus triggered by oleic acid.

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.

The evolution and functional divergence of 10 Apolipoprotein D-like genes in Nilaparvata lugens

Apolipoprotein D (ApoD), a member of the lipocalin superfamily of proteins, is involved in lipid transport and stress resistance. Whereas only a single copy of the ApoD gene is found in humans and some other vertebrates, there are typically several ApoD-like genes in insects. To date, there have been relatively few studies that have examined the evolution and functional differentiation of ApoD-like genes in insects, particularly hemi-metabolous insects. In this study, we identified 10 ApoD-like genes (NlApoD1-10) with distinct spatiotemporal expression patterns in Nilaparvata lugens (BPH), which is an important pest of rice. NlApoD1-10 were found to be distributed on 3 chromosomes in a tandem array of NlApoD1/2, NlApoD3-5, and NlApoD7/8, and show sequence and gene structural divergence in the coding regions, indicating that multiple gene duplication events occurred during evolution. Phylogenetic analysis revealed that NlApoD1-10 can be clustered into 5 clades, with NlApoD3-5 and NlApoD7/8 potentially evolving exclusively in the Delphacidae family. Functional screening using an RNA interference approach revealed that only NlApoD2 was essential for BPH development and survival, whereas NlApoD4/5 are highly expressed in testes, and might play roles in reproduction. Moreover, stress response analysis revealed that NlApoD3-5/9, NlApoD3-5, and NlApoD9 were up-regulated after treatment with lipopolysaccharide, H2 O2 , and ultraviolet-C, respectively, indicating their potential roles in stress resistance.

Genomic signatures of eusocial evolution in insects

The genomes of eusocial insects allow the production and regulation of highly distinct phenotypes, largely independent of genotype. Although rare, eusociality has evolved convergently in at least three insect orders (Hymenoptera, Blattodea and Coleoptera). Despite such disparate origins, eusocial phenotypes show remarkable similarity, exhibiting long-lived reproductives and short-lived sterile workers and soldiers. In this article, we review current knowledge on genomic signatures of eusocial evolution. We confirm that especially an increased regulatory complexity and the adaptive evolution of chemical communication are common to several origins of eusociality. Furthermore, colony life itself can shape genomes of divergent taxa in a similar manner. Future research should be geared towards generating more high-quality genomic resources, especially in hitherto understudied clades, such as ambrosia beetles and termites. The application of more sophisticated tools such as machine learning techniques may allow the detection of more subtle convergent genomic footprints of eusociality.

A genetic toolkit underlying the queen phenotype in termites with totipotent workers

Social insect castes (e.g., queens, workers) are prime examples of phenotypic plasticity (i.e., different phenotypes arising from the same genotype). Yet, the mechanisms that give rise to highly fertile, long-lived queens versus non-reproducing, short-lived workers are not well understood. Recently, a module of co-expressed genes has been identified that characterizes queens compared to workers of the termite Cryptotermes secundus (Kalotermitidae): the Queen Central Module (QCM). We tested whether the QCM is shared in termite species, in which queens gradually develop via early larval and late larval instars, the latter functioning as totipotent workers (linear development). Similar as in C. secundus, gene expression profiles revealed an enrichment of QCM genes in Zootermopsis angusticollis queens, a species from another termite family (Archotermopsidae). The expression of these QCM genes became gradually enriched during development from early larval instars via workers to queens. Thus, our results support the hypothesis of a conserved genetic toolkit that characterizes termite queens with gradual linear development. Our data also imply a strong caste-specific tissue specificity with the QCM signal being restricted to head-prothorax tissues in termite queens. This tissue-specific expression of key aging-related genes might have facilitated the evolution of a long lifespan in termite queens.

Neuropeptide regulations on behavioral plasticity in social insects

Social insects demonstrate remarkable behavioral flexibility in response to complex external and social environments. One of the most striking examples of this adaptability is the context-dependent division of labor among workers of bees and ants. Neuropeptides, the brain's most diverse group of messenger molecules, play an essential role in modulating this phenotypic plasticity related to labor division in social insects. Integrated omics research and mass spectrometry imaging technology have greatly accelerated the identification and spatiotemporal analysis of neuropeptides. Moreover, key roles of several neuropeptides in age- and caste-dependent behavioral plasticity have been uncovered. This review summarizes recent advances in the characterization, expression, distribution, and functions of neuropeptides in controlling behavioral plasticity in social insects, particularly bees and ants. The article concludes with a discussion of future directions and challenges in understanding the regulation of social behavior by neuropeptides.

Nutrition- and hormone-controlled developmental plasticity in Blattodea

Blattodea, which includes cockroaches and termites, possesses high developmental plasticity that is mainly controlled by nutritional conditions and insect hormones. Insulin/insulin-like growth factor signaling (IIS), target of rapamycin complex 1 (TORC1), and adenosine monophosphate-activated protein complex are the three primary nutrition-responsive signals. Juvenile hormone (JH) and 20-hydroxyecdysone (20E) constitute the two most vital insect hormones that might interact with each other through the Met, Kr-h1, E93 (MEKRE93) pathway. Nutritional and hormonal signals interconnect to create a complex regulatory network. Here we summarize recent progress in our understanding of how nutritional and hormonal signals coordinately control the developmental plasticity of metamorphosis, reproduction, and appendage regeneration in cockroaches as well as caste differentiation in termites. We also highlight several perspectives that should be further emphasized in the studies of developmental plasticity in Blattodea. This review provides a general landscape in the field of nutrition- and hormone-controlled developmental plasticity in insects.

Uncovering the Chemosensory System of a Subterranean Termite, Odontotermes formosanus (Shiraki) (Isoptera: Termitidae): Revealing the Chemosensory Genes and Gene Expression Patterns

Termites are eusocial insects. Chemical signals between colony members are crucial to the smooth running of colony operations, but little is known about their olfactory system and the roles played by various chemosensory genes in this process. Chemosensory genes are involved in basic olfactory perception in insects. Odontotermes formosanus (Shiraki) is one of the most damaging pests to agricultural crops, forests, and human-made structures. To better understand the olfactory system and the genes involved in olfactory processing in O. formosanus, we produced a transcriptome of worker termites. In this study, we identified 13 OforOBPs, 1 OforCSP, 15 OforORs, 9 OforGRs, and 4 OforSNMPs. Multiple sequence alignments were used in the phylogenetic study, which included data from other termite species and a wide variety of insect species. Moreover, we also investigated the mRNA expression levels using qRT-PCR. The significantly high expression levels of OforCSP1, OforOBP2, OforOR1, and OforSNMP1 suggest that these genes may play important roles in olfactory processing in termite social behavior, including caste differentiation, nestmate and non-nestmate discrimination, and the performance of colony operations among members. Our research establishes a foundation for future molecular-level functional studies of chemosensory genes in O. formosanus, which might lead to the identification of novel targets for termite integrated pest management.

Live-bearing cockroach genome reveals convergent evolutionary mechanisms linked to viviparity in insects and beyond

Live birth (viviparity) has arisen repeatedly and independently among animals. We sequenced the genome and transcriptome of the viviparous Pacific beetle-mimic cockroach and performed comparative analyses with two other viviparous insect lineages, tsetse flies and aphids, to unravel the basis underlying the transition to viviparity in insects. We identified pathways undergoing adaptive evolution for insects, involved in urogenital remodeling, tracheal system, heart development, and nutrient metabolism. Transcriptomic analysis of cockroach and tsetse flies revealed that uterine remodeling and nutrient production are increased and the immune response is altered during pregnancy, facilitating structural and physiological changes to accommodate and nourish the progeny. These patterns of convergent evolution of viviparity among insects, together with similar adaptive mechanisms identified among vertebrates, highlight that the transition to viviparity requires changes in urogenital remodeling, enhanced tracheal and heart development (corresponding to angiogenesis in vertebrates), altered nutrient metabolism, and shifted immunity in animal systems.

Genome and transcriptome of Ips nitidus provide insights into high-altitude hypoxia adaptation and symbiosis

Ips nitidus is a well-known conifer pest that has contributed significantly to spruce forest disturbance in the Qinghai-Tibet Plateau and seriously threatens the ecological balance of these areas. We report a chromosome-level genome of I. nitidus determined by PacBio and Hi-C technology. Phylogenetic inference showed that it diverged from the common ancestor of I. typographus ∼2.27 mya. Gene family expansion in I. nitidus was characterized by DNA damage repair and energy metabolism, which may facilitate adaptation to high-altitude hypoxia. Interestingly, differential gene expression analysis revealed upregulated genes associated with high-altitude hypoxia adaptation and downregulated genes associated with detoxification after feeding and tunneling in fungal symbiont Ophiostoma bicolor-colonized substrates. Our findings provide evidence of the potential adaptability of I. nitidus to conifer host, high-altitude hypoxia and insight into how fungal symbiont assist in this process. This study enhances our understanding of insect adaptation, symbiosis, and pest management.

Epigenetic weapons in plant-herbivore interactions: Sulforaphane disrupts histone deacetylases, gene expression, and larval development in Spodoptera exigua while the specialist feeder Trichoplusia ni is largely resistant to these effects

Cruciferous plants produce sulforaphane (SFN), an inhibitor of nuclear histone deacetylases (HDACs). In humans and other mammals, the consumption of SFN alters enzyme activities, DNA-histone binding, and gene expression within minutes. However, the ability of SFN to act as an HDAC inhibitor in nature, disrupting the epigenetic machinery of insects feeding on these plants, has not been explored. Here, we demonstrate that SFN consumed in the diet inhibits the activity of HDAC enzymes and slows the development of the generalist grazer Spodoptera exigua, in a dose-dependent fashion. After consuming SFN for seven days, the activities of HDAC enzymes in S. exigua were reduced by 50%. Similarly, larval mass was reduced by 50% and pupation was delayed by 2-5 days, with no additional mortality. Similar results were obtained when SFN was applied topically to eggs. RNA-seq analyses confirm that SFN altered the expression of thousands of genes in S. exigua. Genes associated with energy conversion pathways were significantly downregulated while those encoding for ribosomal proteins were dramatically upregulated in response to the consumption of SFN. In contrast, the co-evolved specialist feeder Trichoplusia ni was not negatively impacted by SFN, whether it was consumed in their diet at natural concentrations or applied topically to eggs. The activities of HDAC enzymes were not inhibited and development was not disrupted. In fact, SFN exposure sometimes accelerated T. ni development. RNA-seq analyses revealed that the consumption of SFN alters gene expression in T. ni in similar ways, but to a lesser degree, compared to S. exigua. This apparent resistance of T. ni can be overwhelmed by unnaturally high levels of SFN or by exposure to more powerful pharmaceutical HDAC inhibitors. These results demonstrate that dietary SFN interferes with the epigenetic machinery of insects, supporting the hypothesis that plant-derived HDAC inhibitors serve as "epigenetic weapons" against herbivores.

iORbase: A database for the prediction of the structures and functions of insect olfactory receptors

Insect olfactory receptors (iORs) with atypical 7-transmembrane domains, unlike Chordata olfactory receptors, are not in the GPCR protein family. iORs selectively bind to volatile ligands in the environment and affect essential insect behaviors. In this study, we constructed a new platform (iORbase, https://www.iorbase.com) for the structural and functional analysis of iORs based on a combined algorithm for gene annotation and protein structure prediction. Moreover, it provides the option to calculate the binding affinities and binding residues between iORs and pheromone molecules by virtual screening of docking. Furthermore, iORbase supports the automatic structural and functional prediction of user-submitted iORs or pheromones. iORbase contains the well-analyzed results of approximately 6 000 iORs and their 3D protein structures identified from 59 insect species and 2 077 insect pheromones from the literature, as well as approximately 12 million pairs of simulated interactions between functional iORs and pheromones. We also built 4 online modules, iORPDB, iInteraction, iModelTM, and iOdorTool to easily retrieve and visualize the 3D structures and interactions. iORbase can help greatly improve the experimental efficiency and success rate, identify new insecticide targets, or develop electronic nose technology. This study will shed light on the olfactory recognition mechanism and evolutionary characteristics from the perspectives of omics and macroevolution.

Epigenetic regulation and division of labor in social insects

Analogous to multicellular organisms, social insect colonies are characterized by division of labor with queens and workers reflecting germline and soma, respectively. In multicellular organisms, such division is achieved through epigenetic factors regulating cell differentiation during development. Analogously, epigenetic regulation is postulated to regulate caste differences in social insects. We summarize recent findings about the role of epigenetics in social insects, focusing on DNA methylation and histone modifications. We specifically address (i) queen versus worker caste differentiation, (ii) queen versus worker caste differences, and (iii) division of labor among workers. Our review provides an overview of an exciting and controversially discussed field in developmental and molecular biology. It shows that our current understanding about the role of epigenetics in regulating division of labor in social insects is still fragmentary but that refined methods with well-replicated samples and targeted questions offer promising insights into this emerging field of socio-epigenomics.

Gene identification and RNAi-silencing of p62/SQSTM1 in the vector Rhodnius prolixus reveals a high degree of sequence conservation but no apparent deficiency-related phenotypes in vitellogenic females

Autophagy and the ubiquitin-proteasome system (UPS) are important cellular mechanisms that coordinate protein degradation essential for proteostasis. P62/SQSTM1 is a receptor cargo protein able to deliver ubiquitinated targets to the proteasome proteolytic complex and/or to the autophagosome. In the insect vector of Chagas disease, Rhodnius prolixus, previous works have shown that the knockdown of different autophagy-related genes (ATGs) and ubiquitin-conjugating enzymes resulted in abnormal oogenesis phenotypes and embryo lethality. Here, we investigate the role of the autophagy/UPS adaptor protein p62 during the oogenesis and reproduction of this vector. We found that R. prolixus presents one isoform of p62 encoded by a non-annotated gene. The predicted protein presents the domain architecture anticipated for p62: PB1 (N-term), ZZ-finger, and UBA (C-term) domains, and phylogenetic analysis showed that this pattern is highly conserved within insects. Using parental RNAi, we found that although p62 is expressed in the ovary, midgut, and fat body of adult females, systemic silencing of this gene did not result in any apparent phenotypes under in-house conditions. The insects' overall levels of blood meal digestion, lifespan, yolk protein production, oviposition, and embryo viability were not altered when compared to controls. Because it is known that autophagy and UPS can undergo compensatory mechanisms, we asked whether the silencing of p62 was triggering adaptative changes in the expression of genes of the autophagy, UPS, and the unfolded protein response (UPR) and found that only ATG1 was slightly up regulated in the ovaries of silenced females. In addition, experiments to further investigate the role of p62 in insects previously silenced for the E1-conjugating enzyme (a condition known to trigger the upregulation of p62), also did not result in any apparent phenotypes in vitellogenic females.

Knockdown of carnitine palmitoyltransferase I (CPT1) reduces fat body lipid mobilization and resistance to starvation in the insect vector Rhodnius prolixus

The energy stored in fatty acids is essential for several critical activities of insects, such as embryogenesis, oviposition, and flight. Rhodnius prolixus is an obligatory hematophagous hemipteran and vector of Chagas disease, and it feeds infrequently on very large blood meals. As digestion slowly occurs, lipids are synthesized and accumulate in the fat body, mainly as triacylglycerol, in lipid droplets. Between feeding bouts, proper mobilization and oxidation of stored lipids are crucial for survival, and released fatty acids are oxidized by mitochondrial β-oxidation. Carnitine palmitoyl transferase I (CPT1) is the enzyme that catalyzes the first reaction of the carnitine shuttle, where the activated fatty acid, acyl-CoA, is converted to acyl-carnitine to be transported into the mitochondria. Here, we investigated the role of CPT1 in lipid metabolism and in resistance to starvation in Rhodnius prolixus. The expression of the CPT1 gene (RhoprCpt1) was determined in the organs of adult females on the fourth day after a blood meal, and the flight muscle showed higher expression levels than the ovary, fat body, and anterior and posterior midgut. RhoprCpt1 expression in the fat body dramatically decreased after feeding, and started to increase again 10 days later, but no changes were observed in the flight muscle. β-oxidation rates were determined in flight muscle and fat body homogenates with the use of ³H-palmitate, and in unfed females, they were higher in the flight muscle. In the fat body, lipid oxidation activity did not show any variation before or at different days after feeding, and was not affected by the presence of etomoxir or malonyl-CoA. We used RNAi and generated RhoprCPT1-deficient insects, which surprisingly did not show a decrease in measured ³H-palmitate oxidation rates. However, the RNAi-knockdown females presented increased amounts of triacylglycerol and larger lipid droplets in the fat body, but not in the flight muscle. When subjected to starvation, these insects had a shorter lifespan. These results indicated that the inhibition of RhoprCpt1 expression compromised lipid mobilization and affected resistance to starvation.

Evidence of cospeciation between termites and their gut bacteria on a geological time scale

Termites host diverse communities of gut microbes, including many bacterial lineages only found in this habitat. The bacteria endemic to termite guts are transmitted via two routes: a vertical route from parent colonies to daughter colonies and a horizontal route between colonies sometimes belonging to different termite species. The relative importance of both transmission routes in shaping the gut microbiota of termites remains unknown. Using bacterial marker genes derived from the gut metagenomes of 197 termites and one Cryptocercus cockroach, we show that bacteria endemic to termite guts are mostly transferred vertically. We identified 18 lineages of gut bacteria showing cophylogenetic patterns with termites over tens of millions of years. Horizontal transfer rates estimated for 16 bacterial lineages were within the range of those estimated for 15 mitochondrial genes, suggesting that horizontal transfers are uncommon and vertical transfers are the dominant transmission route in these lineages. Some of these associations probably date back more than 150 million years and are an order of magnitude older than the cophylogenetic patterns between mammalian hosts and their gut bacteria. Our results suggest that termites have cospeciated with their gut bacteria since first appearing in the geological record.

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.

Evolution of a Cockroach Allergen into the Major Protein of Termite Royal Jelly

Termites live in colonies, and their members belong to different castes that each have their specific role within the termite society. In well-established colonies of higher termites, the only food the founding female, the queen, receives is saliva from workers; such queens can live for many years and produce up to 10,000 eggs per day. In higher termites, worker saliva must thus constitute a complete diet and therein resembles royal jelly produced by the hypopharyngeal glands of honeybee workers that serves as food for their queens; indeed, it might as well be called termite royal jelly. However, whereas the composition of honeybee royal jelly is well established, that of worker termite saliva in higher termites remains largely unknown. In lower termites, cellulose-digesting enzymes constitute the major proteins in worker saliva, but these enzymes are absent in higher termites. Others identified a partial protein sequence of the major saliva protein of a higher termite and identified it as a homolog of a cockroach allergen. Publicly available genome and transcriptome sequences from termites make it possible to study this protein in more detail. The gene coding the termite ortholog was duplicated, and the new paralog was preferentially expressed in the salivary gland. The amino acid sequence of the original allergen lacks the essential amino acids methionine, cysteine and tryptophan, but the salivary paralog incorporated these amino acids, thus allowing it to become more nutritionally balanced. The gene is found in both lower and higher termites, but it is in the latter that the salivary paralog gene got reamplified, facilitating an even higher expression of the allergen. This protein is not expressed in soldiers, and, like the major royal jelly proteins in honeybees, it is expressed in young but not old workers.

Silencing of the 20S proteasomal subunit-α6 triggers full oogenesis arrest and increased mRNA levels of the selective autophagy adaptor protein p62/SQSTM1 in the ovary of the vector Rhodnius prolixus

The high reproductive rates of insects contribute significantly to their ability to act as vectors of a variety of vector-borne diseases. Therefore, it is strategically critical to find molecular targets with biotechnological potential through the functional study of genes essential for insect reproduction. The ubiquitin-proteasome system is a vital degradative pathway that contributes to the maintenance of regular eukaryotic cell proteostasis. This mechanism involves the action of enzymes to covalently link ubiquitin to proteins that are meant to be delivered to the 26S proteasome and broken down. The 26S proteasome is a large protease complex (including the 20S and 19S subcomplexes) that binds, deubiquitylates, unfolds, and degrades its substrates. Here, we used bioinformatics to identify the genes that encode the seven α and β subunits of the 20S proteasome in the genome of R. prolixus and learned that those transcripts are accumulated into mature oocytes. To access proteasome function during oogenesis, we conducted RNAi functional tests employing one of the 20S proteasome subunits (Prosα6) as a tool to suppress 20S proteasomal activity. We found that Prosα6 silencing resulted in no changes in TAG buildup in the fat body and unaffected availability of yolk proteins in the hemolymph of vitellogenic females. Despite this, the silencing of Prosα6 culminated in the impairment of oocyte maturation at the early stages of oogenesis. Overall, we discovered that proteasome activity is especially important for the signals that initiate oogenesis in R. prolixus and discuss in what manner further investigations on the regulation of proteasome assembly and activity might contribute to the unraveling of oogenesis molecular mechanisms and oocyte maturation in this vector.

Caste-specific expressions and diverse roles of takeout genes in the termite Reticulitermes speratus

Acquisition of novel functions caused by gene duplication may be important for termite social evolution. To clarify this possibility, additional evidence is needed. An important example is takeout, encoding juvenile hormone binding protein. We identified 25 takeouts in the termite Reticulitermes speratus genome. RNA-seq revealed that many genes were highly expressed in specific castes. Two novel paralogs (RsTO1, RsTO2) were tandemly aligned in the same scaffold. Real-time qPCR indicated that RsTO1 and RsTO2 were highly expressed in queens and soldiers, respectively. Moreover, the highest RsTO1 expression was observed in alates during queen formation. These patterns were different from vitellogenins, encoding egg-yolk precursors, which were highly expressed in queens than alates. In situ hybridization showed that RsTO1 mRNA was localized in the alate-frontal gland, indicating that RsTO1 binds with secretions probably used for the defence during swarming flight. In contrast, increased RsTO2 expression was observed approximately 1 week after soldier differentiation. Expression patterns of geranylgeranyl diphosphate synthase, whose product functions in the terpenoid synthesis, were similar to RsTO2 expression. In situ hybridization indicated RsTO2-specific mRNA signals in the soldier-frontal gland. RsTO2 may interact with terpenoids, with a soldier-specific defensive function. It may provide additional evidence for functionalization after gene duplication in termites.

ATP synthase affects lipid metabolism in the kissing bug Rhodnius prolixus beyond its role in energy metabolism

ATP synthase plays an essential role in mitochondrial metabolism, being responsible for the production of ATP in oxidative phosphorylation. However, recent results have shown that it may also be present in the cell membrane, involved in lipophorin binding to its receptors. Here, we used a functional genetics approach to investigate the roles of ATP synthase in lipid metabolism in the kissing bug Rhodnius prolixus. The genome of R. prolixus encodes five nucleotide-binding domain genes of the ATP synthase alpha and beta family, including the alpha and beta subunits of ATP synthase (RpATPSynA and RpATPSynB), and the catalytic and non-catalytic subunits of the vacuolar ATPase (RpVha68 and RpVha55). These genes were expressed in all analyzed organs, being their expression highest in the ovaries, fat body and flight muscle. Feeding did not regulate the expression of ATP synthases in the posterior midgut or fat body. Furthermore, ATP synthase is present in the fat body's mitochondrial and membrane fractions. RpATPSynB knockdown by RNAi impaired ovarian development and reduced egg-laying by approximately 85%. Furthermore, the lack of RpATPSynB increased the amount of triacylglycerol in the fat body due to increased de novo fatty acid synthesis and reduced transfer of lipids to lipophorin. RpATPSynA knockdown had similar effects, with altered ovarian development, reduced oviposition, and triacylglycerol accumulation in the fat body. However, ATP synthases knockdown had only a slight effect on the amount of ATP in the fat body. These results support the hypothesis that ATP synthase has a direct role in lipid metabolism and lipophorin physiology, which are not directly due to changes in energy metabolism.

Efficient RNA interference method during caste differentiation with hormone treatment in the termite Reticulitermes speratus (Isoptera: Rhinotermitidae)

Unveiling the proximate mechanism of caste differentiation is crucial for understanding insect social evolution, and gene function analysis is an important tool in this endeavor. The RNA interference (RNAi) technique is useful in termites, but its knockdown effects may differ among species. One of the most important model species in the field of termite sociogenomics is Reticulitermes speratus Kolbe (Isoptera: Rhinotermitidae). Presoldier and worker differentiation of this species can be artificially induced by juvenile hormone and 20-hydroxyecdysone application, respectively. However, appropriate RNAi technique of genes expressed during caste differentiation has never been considered. To clarify this issue, first, we injected nine different volumes of nuclease-free water (NFW, 0-404.8 nL) into workers and found that survival and caste differentiation rates were strongly reduced by the application of the top three largest volumes. Second, we injected double-stranded (ds) RNA of ecdysone receptor homolog (RsEcR) (2.0 µg/151.8 nL NFW) into workers with hormone treatments. The expression levels of RsEcR were significantly reduced at 9 days after dsRNA injection. RsEcR RNAi strongly affected both molting events during presoldier and worker differentiation induced by hormone treatments. The present results highlight the need for caution regarding injection volumes for RNAi experiments using hormone treatments. We suggest that the injection of dsRNA solution (2 µg; approximately 100-200 nL) is suitable for RNAi experiments during caste differentiation induced by hormone application in R. speratus.

Gene expression profiles of chemosensory genes of termite soldier and worker antennae

Termites have an elaborate social system that involves cooperation and division of labour among colony members. Although this social system is regulated by chemical signals produced in the colony, it remains unclear how these signals are perceived by other members. Signal transduction is well known to be triggered by the reception of odorant molecules by some binding proteins in the antennae, after which, a signal is transmitted to chemosensory receptors. However, there is insufficient information on the role of chemosensory genes involved in signal transduction in termites. Here, we identified the genes involved in chemosensory reception in the termite Reticulitermes speratus and performed a genome-wide comparative transcriptome analysis of worker and soldier antennae. First, we identified 31 odorant-binding proteins (OBPs), and three chemosensory protein A (CheA) from the genome data. Thereafter, we performed RNA sequencing to compare the expression levels of OBPs, CheAs, and previously identified chemosensory receptor genes between worker and soldier antennae. There were no receptor genes with significant differences in expression between castes. However, the expression levels of three non-receptor odorant-detection/binding proteins (OBP, CheA, and Sensory neuron membrane protein) were significantly different between castes. Real-time qPCR (RT-qPCR) analysis using antennae and other head parts confirmed that these genes were highly expressed in soldier antennae. Finally, independent RT-qPCR analysis showed that the expression patterns of these genes were altered in soldiers from different social contexts. Present results suggest that gene expression levels of some non-receptors are affected by both castes and behavioural interactions among colony members in termites.

Wood feeding and social living: Draft genome of the subterranean termite Reticulitermes lucifugus (Blattodea; Termitoidae)

Termites (Insecta, Blattodea, Termitoidae) are a widespread and diverse group of eusocial insects known for their ability to digest wood matter. Herein, we report the draft genome of the subterranean termite Reticulitermes lucifugus, an economically important species and among the most studied taxa with respect to eusocial organization and mating system. The final assembly (~813 Mb) covered up to 88% of the estimated genome size and, in agreement with the Asexual Queen Succession Mating System, it was found completely homozygous. We predicted 16,349 highly supported gene models and 42% of repetitive DNA content. Transposable elements of R. lucifugus show similar evolutionary dynamics compared to that of other termites, with two main peaks of activity localized at 25% and 8% of Kimura divergence driven by DNA, LINE and SINE elements. Gene family turnover analyses identified multiple instances of gene duplication associated with R. lucifugus diversification, with significant lineage-specific gene family expansions related to development, perception and nutrient metabolism pathways. Finally, we analysed P450 and odourant receptor gene repertoires in detail, highlighting the large diversity and dynamical evolutionary history of these proteins in the R. lucifugus genome. This newly assembled genome will provide a valuable resource for further understanding the molecular basis of termites biology as well as for pest control.

Two Antenna-Enriched Carboxylesterases Mediate Olfactory Responses and Degradation of Ester Volatiles in the German Cockroach Blattella germanica

Insects have evolved an extremely sensitive olfactory system that is essential for a series of physiological and behavioral activities. Some carboxylesterases (CCEs) comprise a major subfamily of odorant-degrading enzymes (ODEs) playing a crucial role in odorant signal inactivation to maintain the odorant receptor sensitivity. In this study, 93 CCEs were annotated in the genome of the German cockroach Blattella germanica, a serious urban pest. Phylogenetic and digital tissue expression pattern analyses identified two antenna-enriched CCEs, BgerCCE021e3 and BgerCCE021d1, as candidate ODEs. RNA interference (RNAi)-mediated knockdown of BgerCCE021e3 and BgerCCE021d1 resulted in partial anosmia with experimental insects exhibiting reduced attraction to ester volatile resources and slower olfactory responses than controls. Furthermore, enzymatic conversion of geranyl acetate by crude male antennal extracts from BgerCCE021e3 and BgerCCE021d1 RNAi insects was also significantly reduced. Our results provide evidence for CCE function in German cockroach olfaction and provide a basis for further exploring behavioral inhibitors that target olfactory-related CCEs.

A holobiont approach towards polysaccharide degradation by the highly compartmentalised gut system of the soil-feeding higher termite Labiotermes labralis

BACKGROUND

Termites are among the most successful insects on Earth and can feed on a broad range of organic matter at various stages of decomposition. The termite gut system is often referred to as a micro-reactor and is a complex structure consisting of several components. It includes the host, its gut microbiome and fungal gardens, in the case of fungi-growing higher termites. The digestive tract of soil-feeding higher termites is characterised by radial and axial gradients of physicochemical parameters (e.g. pH, O₂ and H₂ partial pressure), and also differs in the density and structure of residing microbial communities. Although soil-feeding termites account for 60% of the known termite species, their biomass degradation strategies are far less known compared to their wood-feeding counterparts.

RESULTS

In this work, we applied an integrative multi-omics approach for the first time at the holobiont level to study the highly compartmentalised gut system of the soil-feeding higher termite Labiotermes labralis. We relied on 16S rRNA gene community profiling, metagenomics and (meta)transcriptomics to uncover the distribution of functional roles, in particular those related to carbohydrate hydrolysis, across different gut compartments and among the members of the bacterial community and the host itself. We showed that the Labiotermes gut was dominated by members of the Firmicutes phylum, whose abundance gradually decreased towards the posterior segments of the hindgut, in favour of Bacteroidetes, Proteobacteria and Verrucomicrobia. Contrary to expectations, we observed that L. labralis gut microbes expressed a high diversity of carbohydrate active enzymes involved in cellulose and hemicelluloses degradation, making the soil-feeding termite gut a unique reservoir of lignocellulolytic enzymes with considerable biotechnological potential. We also evidenced that the host cellulases have different phylogenetic origins and structures, which is possibly translated into their different specificities towards cellulose. From an ecological perspective, we could speculate that the capacity to feed on distinct polymorphs of cellulose retained in soil might have enabled this termite species to widely colonise the different habitats of the Amazon basin.

CONCLUSIONS

Our study provides interesting insights into the distribution of the hydrolytic potential of the highly compartmentalised higher termite gut. The large number of expressed enzymes targeting the different lignocellulose components make the Labiotermes worker gut a relevant lignocellulose-valorising model to mimic by biomass conversion industries.

Enhancement of Pathogen Toxicity by Feeding Reticulitermes chinensis Snyder Sonicated Bacteria Expressing Double-Stranded RNA That Interferes with Olfaction

Reticulitermes chinensis Snyder is a serious pest in China, and the odorant receptor co-receptor gene RcOrco plays a crucial role in olfaction. However, the function of RcOrco in the resistance of termites to entomopathogens has not been reported. We constructed dsRcOrco-HT115 engineered bacteria based on the RcOrco sequence from the full-length transcriptome data of R. chinensis. The engineered bacteria expressed dsRNA of RcOrco. Sonication was used to inactivate the dsRNA-HT115 strain and obtain a large amount of dsRcOrco. The dsRcOrco produced using this method overcame the problem that genetically engineered bacteria could not be applied directly and improved its effectiveness against termites. Bioassays using the dsRcOrco generated using this method showed that dsRcOrco significantly increased the toxicity of the bacterial and fungal pathogens to R. chinensis. The present study showed, for the first time, the function of Orco in termite resistance to pathogens, and the results provide a theoretical basis for the development and application of termite RNA biopesticides.

Conserved orthology in termite chemosensory gene families

Termites are eusocial insects known to use a variety of pheromones in tasks necessary for maintenance of their societies. As such, olfaction and pheromone communication in termites has been an object of intense study; trail-following pheromones (TFPs) and sex-pairing pheromones (SPPs), for example, have been identified in many termite species. In contrast, the molecular basis of olfactory detection is understudied in the group. Here, we present chemosensory genes of three species of termites belonging to three distinct lineages, Neotermes cubanus (Kalotermitidae), Prorhinotermes simplex (Rhinotermitidae), and Inquilinitermes inquilinus (Termitidae). Using antennal transcriptome screening of termite workers, we identified the chemosensory genes, which allowed us to perform phylogenetic analysis. We found a comparatively large repertoires of odorant receptors (ORs), gustatory receptors (GRs), ionotropic receptors (IRs), odorant binding proteins (OBPs), chemosensory proteins (CSPs), and sensory neuron membrane proteins (SNMPs). The evolutionary analysis of termite chemosensory genes revealed Isoptera-specific expansions with a 1:1 orthologous pattern, indicating the existence of conserved olfactory functions. Our findings on basal eusocial insects will further enhance our understanding of the molecular underpinnings of eusociality and the evolution of olfactory communication in termites.

Epigenetic regulations as drivers of insecticide resistance and resilience to climate change in arthropod pests

Arthropod pests are remarkably capable of rapidly adapting to novel forms of environmental stress, including insecticides and climate change. The dynamic interplay between epigenetics and genetics explains the largely unexplored reality underlying rapid climatic adaptation and the development of insecticide resistance in insects. Epigenetic regulation modulates gene expression by methylating DNA and acetylating histones that play an essential role in governing insecticide resistance and adaptation to climate change. This review summarises and discusses the significance of recent advances in epigenetic regulation that facilitate phenotypic plasticity in insects and their symbiotic microbes to cope with selection pressure implied by extensive insecticide applications and climate change. We also discuss how epigenetic changes are passed on to multiple generations through sexual recombination, which remains enigmatic. Finally, we explain how these epigenetic signatures can be utilized to manage insecticide resistance and pest resilience to climate change in Anthropocene.

Differential expression of some termite neuropeptides and insulin/IGF-related hormones and their plausible functions in growth, reproduction and caste determination

Background

Insulin-like growth factor (IGF) and other insulin-like peptides (ilps) are important hormones regulating growth and development in animals. Whereas most animals have a single female and male adult phenotype, in some insect species the same genome may lead to different final forms. Perhaps the best known example is the honeybee where females can either develop into queens or workers. More extreme forms of such polyphenism occur in termites, where queens, kings, workers and soldiers coexist. Both juvenile hormone and insulin-like peptides are known to regulate growth and reproduction as well as polyphenism. In termites the role of juvenile hormone in reproduction and the induction of the soldier caste is well known, but the role of IGF and other ilps in these processes remains largely unknown. Here the various termite ilps are identified and hypotheses regarding their functions suggested.

Methods

Genome assemblies and transcriptome short read archives (SRAs) were used to identify insulin-like peptides and neuropeptides in termites and to determine their expression in different species, tissues and castes.

Results and Discussion

Termites have seven different ilps, i.e. gonadulin, IGF and an ortholog of Drosophila insulin-like peptide 7 (dilp7), which are commonly present in insects, and four smaller peptides, that have collectively been called short IGF-related peptides (sirps) and individually atirpin, birpin, cirpin and brovirpin. Gonadulin is lost from the higher termites which have however amplified the brovirpin gene, of which they often have two or three paralogs. Based on differential expression of these genes it seems likely that IGF is a growth hormone and atirpin an autocrine tissue factor that is released when a tissue faces metabolic stress. Birpin seems to be responsible for growth and in the absence of juvenile hormone this may lead to reproductive adults or, when juvenile hormone is present, to soldiers. Brovirpin is expressed both by the brain and the ovary and likely stimulates vitellogenesis, while the function of cirpin is less clear.

Evolution and functionalization of vitellogenin genes in the termite Reticulitermes speratus

Eusociality has been commonly observed in distinct animal lineages. The reproductive division of labor is a particular feature, achieved by the coordination between fertile and sterile castes within the same nest. The sociogenomic approach in social hymenopteran insects indicates that vitellogenin (Vg) has undergone neo-functionalization in sterile castes. Here, to know whether Vgs have distinct roles in nonreproductive castes in termites, we investigated the unique characteristics of Vgs in the rhinotermitid termite Reticulitermes speratus. The four Vgs were identified from R. speratus (RsVg1-4), and RsVg3 sequences were newly identified using the RACE method. Molecular phylogenetic analysis supported the monophyly of the four termite Vgs. Moreover, the termites Vg1-3 and Vg4 were positioned in two different clades. The  dN/dS ratios indicated that the branch leading to the common ancestor of termite Vg4 was under weak purifying selection. Expression analyses among castes (reproductives, workers, and soldiers) and females (nymphs, winged alates, and queens) showed that RsVg1-3 was highly expressed in fertile queens. In contrast, RsVg4 was highly expressed in workers and female nonreproductives (nymphs and winged adults). Localization of RsVg4 messenger RNA was confirmed in the fat body of worker heads and abdomens. These results suggest that Vg genes are functionalized after gene duplication during termite eusocial transition and that Vg4 is involved in nonreproductive roles in termites.

The antennal transcriptome of Triatoma infestans reveals substantial expression changes triggered by a blood meal

BACKGROUND

Triatoma infestans is the main vector of Chagas disease in the Americas, currently transmitting it in Argentina, Paraguay, and Bolivia. Many T. infestans populations present insecticide resistance, reducing the efficiency of control campaigns. Alternative vector control methods are needed, and molecular targets mediating fundamental physiological processes can be a promising option to manipulate kissing bug behavior. Therefore, it is necessary to characterize the main sensory targets, as well as to determine whether they are modulated by physiological factors. In order to identify gene candidates potentially mediating host cue detection, the antennal transcripts of T. infestans fifth instar larvae were sequenced and assembled. Besides, we evaluated whether a blood meal had an effect on transcriptional profiles, as responsiveness to host-emitted sensory cues depends on bug starvation.

RESULTS

The sensory-related gene families of T. infestans were annotated (127 odorant receptors, 38 ionotropic receptors, 11 gustatory receptors, 41 odorant binding proteins, and 25 chemosensory proteins, among others) and compared to those of several other hemipterans, including four triatomine species. Several triatomine-specific lineages representing sensory adaptations developed through the evolution of these blood-feeding heteropterans were identified. As well, we report here various conserved sensory gene orthogroups shared by heteropterans. The absence of the thermosensor pyrexia, of pickpocket receptor subfamilies IV and VII, together with clearly expanded takeout repertoires, are revealed features of the molecular bases of heteropteran antennal physiology. Finally, out of 2,122 genes whose antennal expression was significantly altered by the ingestion of a blood meal, a set of 41 T. infestans sensory-related genes (9 up-regulated; 32 down-regulated) was detected.

CONCLUSIONS

We propose that the set of genes presenting nutritionally-triggered modulation on their expression represent candidates to mediate triatomine host-seeking behavior. Besides, the triatomine-specific gene lineages found represent molecular adaptations to their risky natural history that involves stealing blood from an enormously diverse set of vertebrates. Heteropteran gene orthogroups identified may represent unknown features of the sensory specificities of this largest group of hemipteroids. Our work is the first molecular characterization of the peripheral modulation of sensory processes in a non-dipteran vector of human disease.