Hemimetabolous genomes reveal molecular basis of termite eusociality

Differential Expression of Hormone-Related Genes in the Heads of Adult and Nymphal Woodroaches (Cryptocercus)

Termites are eusocial cockroaches, but the crucial distinctions in gene expression during the evolution of eusociality remain unclear. One reason for the lack of this information is that comparative transcriptome analysis of termites with their sister group, the cockroach genus Cryptocercus, has not been conducted. We identified genes associated with three vital hormones (juvenile hormone [JH], 20-hydoroxyecdysone [20E], and insulin) from the genome sequence of Cryptocercus punctulatus and conducted RNA-seq analysis using the heads of female/male adults and nymphs to elucidate their expression levels. The comprehensive gene expression analysis revealed a multitude of genes exhibiting differences in expression between developmental stages rather than between sexes. Subsequently, we compared the differences in expression patterns of each hormone-related gene by combining the results of a previous RNA-seq study conducted on the heads of castes (reproductives, workers, and soldiers) in the termite Reticulitermes speratus. The results indicated that genes with expression differences among castes in R. speratus, particularly those related to JH and 20E, were significantly more abundant compared to genes with expression differences between adults and nymphs in C. punctulatus. While no significant difference was observed in the number of genes within the insulin signaling pathway, a trend of homologs highly expressed in adult woodroaches but not in adult termites was observed, and the expression patterns of positive and negative regulators in the pathway differed significantly between adults and nymphs. The differences in the expression patterns between Cryptocercus and termites are believed to reflect variations in hormone levels and signaling activities between adults and juveniles, the latter encompassing workers and soldiers in the case of termites.

Unexpectedly low recombination rates and presence of hotspots in termite genomes

Meiotic recombination is a fundamental evolutionary process that facilitates adaptation and the removal of deleterious genetic variation. Social Hymenoptera exhibit some of the highest recombination rates among metazoans, whereas high recombination rates have not been found among nonsocial species from this insect order. It is unknown whether elevated recombination rates are a ubiquitous feature of all social insects. In many metazoan taxa, recombination is mainly restricted to hotspots a few kilobases in length. However, little is known about the prevalence of recombination hotspots in insect genomes. Here we infer recombination rate and its fine-scale variation across the genomes of two social species from the insect order Blattodea: the termites Macrotermes bellicosus and Cryptotermes secundus We used linkage disequilibrium-based methods to infer recombination rate. We infer that recombination rates are close to 1 cM/Mb in both species, similar to the average metazoan rate. We also observe a highly punctate distribution of recombination in both termite genomes, indicative of the presence of recombination hotspots. We infer the presence of full-length PRDM9 genes in the genomes of both species, which suggests recombination hotspots in termites might be determined by PRDM9, as they are in mammals. We also find that recombination rates in genes are correlated with inferred levels of germline DNA methylation. The finding of low recombination rates in termites indicates that eusociality is not universally connected to elevated recombination rate. We speculate that the elevated recombination rates in social Hymenoptera are instead promoted by intense selection among haploid males.

Repeated behavioural evolution is associated with convergence of gene expression in cavity-nesting songbirds

Uncovering the genomic bases of phenotypic adaptation is a major goal in biology, but this has been hard to achieve for complex behavioural traits. Here we leverage the repeated, independent evolution of obligate cavity nesting in birds to test the hypothesis that pressure to compete for a limited breeding resource has facilitated convergent evolution in behaviour, hormones and gene expression. We used an integrative approach, combining aggression assays in the field, testosterone measures and transcriptome-wide analyses of the brain in wild-captured females and males. Our experimental design compared species pairs across five avian families, each including one obligate cavity-nesting species and a related species with a more flexible nest strategy. We find behavioural convergence, with higher levels of territorial aggression in obligate cavity nesters, particularly among females. Across species, levels of testosterone in circulation were not associated with nest strategy nor aggression. Phylogenetic analyses of individual genes and co-regulated gene networks revealed more shared patterns of brain gene expression than expected by drift, although the scope of convergent gene expression evolution was limited to a small percentage of the genome. When comparing our results to other studies that did not use phylogenetic methods, we suggest that accounting for shared evolutionary history may reduce the number of genes inferred as convergently evolving. Altogether, we find that behavioural convergence in response to shared ecological pressures is associated with largely independent evolution of gene expression across different avian families, punctuated by a narrow set of convergently evolving genes.

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.

Molecular basis of E93-dependent tissue morphogenesis and histolysis during insect metamorphosis

The evolution of insect metamorphosis has profoundly influenced their successful adaptation and diversification. Two key physiological processes during insect metamorphosis are notable: wing maturation and prothoracic gland (PG) histolysis. The ecdysone-induced protein 93 (E93) is a transcription factor indispensable for metamorphosis. While it has been established that both wing maturation and PG histolysis are dependent on E93, the molecular mechanisms through which E93 regulates these seemingly 'opposing' events remain poorly understood. In this study, time-course transcriptome profiles were generated for wing pads and PGs during metamorphosis in Blattella germanica, a hemimetabolous model insect. Comparative transcriptomic analyses demonstrated that E93 exerts predominant control over extensive gene transcription during wing morphogenesis and PG histolysis. During wing morphogenesis, E93 selectively enhances the expression of genes associated with cell proliferation, energy supply, signal transduction, actin cytoskeleton organization, and cell adhesion, etc. Additionally, E93 activates the transcription of the majority of genes within the wing gene network that are crucial for wing development in B. germanica. During PG histolysis, E93 preferentially promotes the expression of genes related to endocytosis, focal adhesion, the AMPK signaling pathway, adipocytokine signaling pathway, Toll and Imd signaling pathways, and autophagy, etc. The key genes involved in the aforementioned pathways were subsequently confirmed to contribute to the E93-dependent degeneration of the PG in B. germanica. In summary, our results reveal that E93 functions as a master transcriptional regulator orchestrating both tissue morphogenesis and histolysis during insect metamorphosis. These findings contribute to a deeper understanding of the genetic underpinnings of insect metamorphosis.

Synthesis and biological evaluation of rationally designed pyrazoles as insecticidal agents

The current research focused on the synthesis of two series of pyrazole derivatives and evaluation of their insecticidal effectiveness. In the first series, seven pyrazole Schiff bases 3a-g were successfully synthesized with yields (79-95%) by condensing phenylfuran-2-carbaldehyde with substituted pyrazole rings. In the second series, eleven amino acid-pyrazole conjugates 6a-k were synthesized utilizing acetic acid, sulfuric acid, morpholine, and EDC. HCl achieving yields of 59% to 94%. The synthesized compounds were assessed for their chemotherapeutic efficacy against locusts and termites by calculating LC50 values, thereby determining their potential as anti-termite and anti-locust agents. Among the eighteen synthesized pyrazole compounds, the Schiff base pyrazole molecules 3f (LC50 = 0.001 μg/mL) and 3d (LC50 = 0.006 μg/mL) demonstrated excellent anti-termite activity compared to the reference drug fipronil (LC50 = 0.038 μg/mL). Pyrazole derivative 6 h with LC50 = 47.68 μg/mL exhibited superior anti-locust activity than the reference drug fipronil (LC50 = 63.09 μg/mL). Additionally, compound 3b, containing NO2 functionality at the meta position, exhibited notable and significant anti-locust activity with an LC50 values of 100.00 μg/mL. However, the highest mortality was caused by the glycine conjugate of fipronil 6 h of the 2nd series with an LC50 value of 47.68 μg/mL, which also proved to be a potent anti-locust agent. This study explores the efficacy of biologically active pyrazole structures as potential insecticidal agents through a combination of virtual molecular docking analysis and biological experimental investigations. The results demonstrate a strong correlation between the computational predictions and experimental outcomes, suggesting that the pyrazole derivatives exhibit significant insecticidal properties. The findings highlight the potential of these compounds in the development of innovative insecticides, paving the way for future research in pest control strategies.

Binding properties of olfactory proteins to host volatiles, free fatty acids and cuticular hydrocarbons in the termite Reticulitermes aculabialis

As eusocial insects prevalent in tropical and subtropical regions, termites are characterized by highly organized behaviors and exceptional adaptability, rooted in caste differentiation and chemical communication. These traits make them excellent models for studying insect social structures and ecological interactions. Investigating how termites use chemical signals to perceive and respond to their environment provides insights into their coordination and adaptation within complex ecosystems. This study delved into the chemosensory mechanisms of Reticulitermes aculabialis, examining the interactions of four olfactory proteins with 70 ligands, including host volatiles, cuticular hydrocarbons (CHCs), and free fatty acids (FFAs). Molecular docking simulations revealed varied affinities of the olfactory proteins for long-chain hydrocarbons (n-C23 to n-C28), suggesting a nuanced chemical communication system through specific hydrocarbon detection. RacuCSP1 and RacuCSP2 exhibited specific binding to linoleic acid and undecanoic acid, respectively, highlighting the significance of FFAs in the physiological and behavioral processes of termites. The four olfactory proteins showed a strong affinity for longifolene in fluorescence competitive binding experiments. Notably, RacuOBPs exhibited unique affinities for terpenoid volatiles such as β-lonone and neocembrene, while RacuCSPs specifically bound with terpenoids like 3-carene, myrtenol, α-pinene oxide and β-pinene indicating their critical roles in host detection. Behavioral observations following gene silencing revealed that RacuOBP5 was essential for recognizing longifolene and α-lonone recognition, while RacuCSP1 was key for detecting α-pinene in termites. These findings enhance our understanding of the termite chemosensory system and offer insights for developing precise pest management strategies.

Non-gonadal expression of piRNAs is widespread across Arthropoda

PIWI-interacting RNAs (piRNAs) were discovered in the early 2000s and became known for their role in protecting the germline genome against mobile genetic elements. Successively, piRNAs were also detected in the somatic cells of gonads in multiple animal species. In recent years, piRNAs have been reported in non-gonadal tissues in various arthropods, contrary to the initial assumptions of piRNAs being exclusive to gonads. Here, we performed an extensive literature review, which revealed that reports on non-gonadal somatic piRNA expression are not limited to a few specific species. Instead, when multiple studies are considered collectively, it appears to be a widespread phenomenon across arthropods. Furthermore, we systematically analyzed 168 publicly available small RNA-seq datasets from diverse tissues in 17 species, which further supported the bibliographic reports that piRNAs are expressed across tissues and species in Arthropoda.

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.

Genetic Underpinnings of Cuticular Hydrocarbon Biosynthesis in the German Cockroach, Blattella germanica (L.): Progress and Perspectives

Insect cuticular hydrocarbons (CHCs) serve as important waterproofing barriers and as signals and cues in chemical communication. Over the past 30 years, numerous studies on CHCs have been conducted in the German cockroach, Blattella germanica, leading to substantial progress in the field. However, there has not been a systematic review of CHC studies in this species in recent years. This review aims to provide a concise overview of the chemical composition, storage, transport, and physical properties of different CHCs in B. germanica. Additionally, we focus on the biosynthetic pathway and the genetic regulation of HC biosynthesis in this species. A considerable amount of biochemical evidence regarding the biosynthetic pathway of insect CHCs has been gathered from studies conducted in B. germanica. In recent years, there has also been an improved understanding of the molecular mechanisms that underlie CHC production in this insect. In this article, we summarize the biosynthesis of different classes of CHCs in B. germanica. Then, we review CHCs reaction to various environmental conditions and stressors and internal physiological states. Additionally, we review a body of work showing that in B. germanica, CHC profiles exhibit significant sexual dimorphism, specific CHCs act as essential precursors for female contact sex pheromone components, and we summarize the molecular regulatory mechanisms that underlie sexual dimorphism of CHC profiles. Finally, we highlight future directions and challenges in research on the biosynthesis and regulatory mechanisms of CHCs in B. germanica, and also identify potential applications of CHC studies in the pest control.

Histone Deacetylase 3 Is Involved in Maintaining Queen Hallmarks of a Termite

The role of epigenetics in regulating caste polyphenism in social insects has been debated. Here, we tested the importance of histone de/acetylation processes for the maintenance of queen hallmarks like a high fecundity and a long lifespan. To this end, we performed RNA interference experiments against histone deacetylase 3 (HDAC3) in the termite Cryptotermes secundus. Fat body transcriptomes and chemical communication profiles revealed that silencing of HDAC3 leads to signals indicative of queen hallmarks. This includes fostering of queen signalling, defence against ageing and a reduction of life-shortening IIS (insulin/insulin-like growth factor signalling) and endocrine JH (juvenile hormone) signalling via Kr-h1 (Krüppel-homologue 1). These observed patterns were similar to those of a protein-enriched diet, which might imply that histone acetylation conveys nutritional effects. Strikingly, in contrast to solitary insects, reduced endocrine JH signalling had no negative effect on fecundity-related vitellogenesis in the fat bodies. This suggests an uncoupling of longevity pathways from fecundity in fat bodies, which can help explain queens' extraordinary lifespans combined with high fecundity.

E93 is indispensable for reproduction in ametabolous and hemimetabolous insects

Ecdysone-induced protein 93 (E93), known as the 'adult-specifier' transcription factor in insects, triggers metamorphosis in both hemimetabolous and holometabolous insects. Although E93 is conserved in ametabolous insects, its spatiotemporal expression and physiological function remain poorly understood. In this study, we first discover that, in the ametabolous firebrat Thermobia domestica, the previtellogenic ovary exhibits cyclically high E93 expression, and E93 mRNA is broadly distributed in previtellogenic ovarioles. E93 homozygous mutant females of T. domestica exhibit severe fecundity deficiency due to impaired previtellogenic development of the ovarian follicles, likely because E93 induces the expression of genes involved in ECM (extracellular matrix)-receptor interactions during previtellogenesis. Moreover, we reveal that in the hemimetabolous cockroach Blattella germanica, E93 similarly promotes previtellogenic ovarian development. In addition, E93 is also essential for vitellogenesis that is necessary to guarantee ovarian maturation and promotes the vitellogenesis-previtellogenesis switch in the fat body of adult female cockroaches. Our findings deepen the understanding of the roles of E93 in controlling reproduction in insects, and of E93 expression and functional evolution, which are proposed to have made crucial contributions to the origin of insect metamorphosis.

An integrative genomic toolkit for studying the genetic, evolutionary, and molecular underpinnings of eusociality in insects

While genomic resources for social insects have vastly increased over the past two decades, we are still far from understanding the genetic and molecular basis of eusociality. Here, we briefly review three scientific advancements that, when integrated, can be highly synergistic for advancing our knowledge of the genetics and evolution of eusocial traits. Population genomics provides a natural way to quantify the strength of natural selection on coding and regulatory sequences, highlighting genes that have undergone adaptive evolution during the evolution or maintenance of eusociality. Genome-wide association studies (GWAS) can be used to characterize the complex genetic architecture underlying eusocial traits and identify candidate causal variants. Concurrently, CRISPR/Cas9 enables the precise manipulation of gene function to both validate genotype-phenotype associations and study the molecular biology underlying interesting traits. While each approach has its own advantages and disadvantages, which we discuss herein, we argue that their combination will ultimately help us better understand the genetics and evolution of eusocial behavior. Specifically, by triangulating across these three different approaches, researchers can directly identify and study loci that have a causal association with key phenotypes and have evidence of positive selection over the relevant timescales associated with the evolution and maintenance of eusociality in insects.

Functional redundancy of the three insulin receptors of cockroaches

Gene duplication is a fundamental evolutionary process which provides opportunities to acquire new gene functions. In the case of the insulin receptors (InRs) in cockroaches and close-related insects, two successive duplications determined the occurrence of three InR genes: InR2, InR1 and InR3, the last two forming a sister cluster to InR2. The biological role of each of the gene duplicates and whether they resulted from neofunctionalization or subfunctionalization is still unclear. The analysis of the sequences from different lineages did not detect positive selection as driving the divergence of InR1 and InR3, discarding neofunctionalization, and suggesting that there is no functional divergence between both gene copies. Using the cockroach Blattella germanica as a model, we have determined that BgInR2 is the gene with the highest expression levels in all the tissues analyzed, both in adult females and males, as well as in nymphs and embryos. BgInR3 is second in expression levels while BgInR1 is expressed at lower levels and only in some tissues. The selective depletion by RNAi of each of the three InRs, analyzed in terms of phenotype and fat body transcriptomic profiles, resulted in essentially redundant effects, with a magnitude approximately proportional to the level of expression of the respective InR. Therefore, the results indicate that the InR duplicates likely experienced a subfunctionalization process, by which the three InRs maintained similar functions but contributing to those functions proportionally to their expression levels.

Social aphids: emerging model for studying insect sociality

Sociality is also evolved in parthenogenetic herbivorous hemipteran aphids, encompassing species with complex life history traits and significant social diversity. Owing to their interesting biological characteristics comparing to other social insect groups, social aphids can be a good model for studying insect sociality. Here, we review the species, behavior, and trait diversity of social aphids, and present recent findings on environmental, physiological, and molecular regulations of caste differentiation and behavior in social aphids. We propose the unique value of social aphids in investigating the evolution and mechanisms of insect sociality as well as future research directions using the social aphid model, including social evolution, caste differentiation, behavioral polymorphism, morphological plasticity, physical mechanics, and interspecific interactions.

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.

Glucose aversion: a behavioral resistance mechanism in the German cockroach

The German cockroach is a valuable model for research on indoor pest management strategies and for understanding mechanisms of adaptive evolution under intense anthropogenic selection. Under the selection pressure of toxic baits, populations of the German cockroach have evolved a variety of physiological and behavioral resistance mechanisms. In this review, we focus on glucose aversion, an adaptive trait that underlies a behavioral resistance to baits. Taste polymorphism, a change in taste quality of glucose from sweet to bitter, causes cockroaches to avoid glucose-containing baits. We summarize recent findings, including the contribution of glucose aversion to olfactory learning-based avoidance of baits, aversion to other sugars, and assortative mating under sexual selection, which underscores the behavioral phenotype to all oligosaccharides that contain glucose. It is a remarkable example of how anthropogenic selection drove the evolution of an altered gustatory trait that reshapes the foraging ecology and sexual communication.

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.

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.

Comparative Transcriptomics of Fat Bodies between Symbiotic and Quasi-Aposymbiotic Adult Females of Blattella germanica with Emphasis on the Metabolic Integration with Its Endosymbiont Blattabacterium and Its Immune System

We explored the metabolic integration of Blattella germanica and its obligate endosymbiont Blattabacterium cuenoti by the transcriptomic analysis of the fat body of quasi-aposymbiotic cockroaches, where the endosymbionts were almost entirely removed with rifampicin. Fat bodies from quasi-aposymbiotic insects displayed large differences in gene expression compared to controls. In quasi-aposymbionts, the metabolism of phenylalanine and tyrosine involved in cuticle sclerotization and pigmentation increased drastically to compensate for the deficiency in the biosynthesis of these amino acids by the endosymbionts. On the other hand, the uricolytic pathway and the biosynthesis of uric acid were severely decreased, probably because the reduced population of endosymbionts was unable to metabolize urea to ammonia. Metabolite transporters that could be involved in the endosymbiosis process were identified. Immune system and antimicrobial peptide (AMP) gene expression was also reduced in quasi-aposymbionts, genes encoding peptidoglycan-recognition proteins, which may provide clues for the maintenance of the symbiotic relationship, as well as three AMP genes whose involvement in the symbiotic relationship will require additional analysis. Finally, a search for AMP-like factors that could be involved in controlling the endosymbiont identified two orphan genes encoding proteins smaller than 200 amino acids underexpressed in quasi-aposymbionts, suggesting a role in the host-endosymbiont relationship.

RNAi of yellow-y, required for normal cuticle pigmentation, impairs courtship behavior and oviposition in the German cockroach (Blattella germanica)

The insect cuticle plays a key role in maintaining the insect's physiological function and behavior. Herein, the yellow-y protein is required to produce black melanin, and is expressed in a pattern that correlates with the distribution of this pigment. However, yellow-y can also have other functions, for instance, in insect behavior, but not much is known. In this study, we have studied the yellow-y gene in one important model and pest species, namely the German cockroach (Blattella germanica), which is to our knowledge the first time reported. In essence, we identified the yellow-y gene (BgY-y) and characterized its function by using RNA interference (RNAi). Silencing of BgY-y gene led to different developmental abnormalities (body weight and wings) in both genders. Specifically, there was an abundant decrease in melanin, turning the body color in pale yellow and the cuticle softer and more transparent. Interestingly, we also observed that the knockdown of BgY-y impaired the male cockroaches to display a weaker response to female-emitted contact sex pheromones, and also that the oviposition ability was weakened in the RNAi females. This study comprehensively analyzed the biological functions of the yellow-y gene in German cockroaches from the perspectives of development, body color, courtship behavior and oviposition, and as a consequence, this may opens new avenues to explore it as a novel pest control gene.

The melanin pigment gene black mediates body pigmentation and courtship behaviour in the German cockroach Blattella germanica

Genes involved in melanin production directly impact insect pigmentation and can affect diverse physiology and behaviours. The role these genes have on sex behaviour, however, is unclear. In the present study, the crucial melanin pigment gene black was functionally characterised in an urban pest, the German cockroach, Blattella germanica. RNAi knockdown of B. germanica black (Bgblack) had no effect on survival, but did result in black pigmentation of the thoraxes, abdomens, heads, wings, legs, antennae, and cerci due to cuticular accumulation of melanin. Sex-specific variation in the pigmentation pattern was apparent, with females exhibiting darker coloration on the abdomen and thorax than males. Bgblack knockdown also resulted in wing deformation and negatively impacted the contact sex pheromone-based courtship behaviour of males. This study provides evidence for black function in multiple aspects of B. germanica biology and opens new avenues of exploration for novel pest control strategies.

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.

Gene regulation supporting sociality shared across lineages and variation in complexity

Across evolutionary lineages, insects vary in social complexity, from those that exhibit extended parental care to those with elaborate divisions of labor. Here, we synthesize the sociogenomic resources from hundreds of species to describe common gene regulatory mechanisms in insects that regulate social organization across phylogeny and levels of social complexity. Different social phenotypes expressed by insects can be linked to the organization of co-expressing gene networks and features of the epigenetic landscape. Insect sociality also stems from processes like the emergence of parental care and the decoupling of ancestral genetic programs. One underexplored avenue is how variation in a group's social environment affects the gene expression of individuals. Additionally, an experimental reduction of gene expression would demonstrate how the activity of specific genes contributes to insect social phenotypes. While tissue specificity provides greater localization of the gene expression underlying social complexity, emerging transcriptomic analysis of insect brains at the cellular level provides even greater resolution to understand the molecular basis of social insect evolution.

Neuroregulation of foraging behavior mediated by the olfactory co-receptor Orco in termites

Olfaction is critical for survival because it allows animals to look for food and detect pheromonal cues. Neuropeptides modulate olfaction and behaviors in insects. While how the neuroregulation of olfactory recognition affects foraging behavior in termites is still unclear. Here, we analyzed the change after silencing the olfactory co-receptor gene (Orco) and the neuropeptide Y gene (NPY), and then investigated the impact of olfactory recognition on foraging behavior in Odontotermes formosanus under different predation pressures. The knockdown of Orco resulted in the reduced Orco protein expression in antennae and the decreased EAG response to trail pheromones. In addition, NPY silencing led to the damaged ability of olfactory response through downregulating Orco expression. Both dsOrco- and dsNPY-injected worker termites showed significantly reduced walking activity and foraging success. Additionally, we found that 0.1 pg/cm trail pheromone and nestmate soldiers could provide social buffering to relieve the adverse effect of predator ants on foraging behavior in worker termites with the normal ability of olfactory recognition. Our orthogonal experiments further verified that Orco/NPY genes are essential in manipulating termite olfactory recognition during foraging under different predation pressures, suggesting that the neuroregulation of olfactory recognition plays a crucial role in regulating termite foraging behavior.

The influence of first desaturase subfamily genes on fatty acid synthesis, desiccation tolerance and inter-caste nutrient transfer in the termite Coptotermes formosanus

Desaturase enzymes play an essential role in the biosynthesis of unsaturated fatty acids (UFAs). In this study, we identified seven "first desaturase" subfamily genes (Cfor-desatA1, Cfor-desatA2-a, Cfor-desatA2-b, Cfor-desatB-a, Cfor-desatB-b, Cfor-desatD and Cfor-desatE) from the Formosan subterranean termite Coptotermes formosanus. These desaturases were highly expressed in the cuticle and fat body of C. formosanus. Inhibition of either the Cfor-desatA2-a or Cfor-desatA2-b gene resulted in a significant decrease in the contents of fatty acids (C16:0, C18:0, C18:1 and C18:2) in worker castes. Moreover, we observed that inhibition of most of desaturase genes identified in this study had a negative impact on the survival rate and desiccation tolerance of workers. Interestingly, when normal soldiers were reared together with dsCfor-desatA2-b-treated workers, they exhibited higher mortality, suggesting that desaturase had an impact on trophallaxis among C. formosanus castes. Our findings shed light on the novel roles of desaturase family genes in the eusocial termite C. formosanus.

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.

Dynamic Evolution of Repetitive Elements and Chromatin States in Apis mellifera Subspecies

In this study, we elucidate the contribution of repetitive DNA sequences to the establishment of social structures in honeybees (Apis mellifera). Despite recent advancements in understanding the molecular mechanisms underlying the formation of honeybee castes, primarily associated with Notch signaling, the comprehensive identification of specific genomic cis-regulatory sequences remains elusive. Our objective is to characterize the repetitive landscape within the genomes of two honeybee subspecies, namely A. m. mellifera and A. m. ligustica. An observed recent burst of repeats in A. m. mellifera highlights a notable distinction between the two subspecies. After that, we transitioned to identifying differentially expressed DNA elements that may function as cis-regulatory elements. Nevertheless, the expression of these sequences showed minimal disparity in the transcriptome during caste differentiation, a pivotal process in honeybee eusocial organization. Despite this, chromatin segmentation, facilitated by ATAC-seq, ChIP-seq, and RNA-seq data, revealed a distinct chromatin state associated with repeats. Lastly, an analysis of sequence divergence among elements indicates successive changes in repeat states, correlating with their respective time of origin. Collectively, these findings propose a potential role of repeats in acquiring novel regulatory functions.

Comparative Evolutionary Genomics in Insects

Genome sequencing quality, in terms of both read length and accuracy, is constantly improving. By combining long-read sequencing technologies with various scaffolding techniques, chromosome-level genome assemblies are now achievable at an affordable price for non-model organisms. Insects represent an exciting taxon for studying the genomic underpinnings of evolutionary innovations, due to ancient origins, immense species-richness, and broad phenotypic diversity. Here we summarize some of the most important methods for carrying out a comparative genomics study on insects. We describe available tools and offer concrete tips on all stages of such an endeavor from DNA extraction through genome sequencing, annotation, and several evolutionary analyses. Along the way we describe important insect-specific aspects, such as DNA extraction difficulties or gene families that are particularly difficult to annotate, and offer solutions. We describe results from several examples of comparative genomics analyses on insects to illustrate the fascinating questions that can now be addressed in this new age of genomics research.

Dual function of the transcription factor Ftz-f1 on oviposition in the cockroach Blattella germanica

The transcription factor Ftz-f1 has multiple functions in insect development in a spatial-temporal line. One of these roles is in the insect ovaries, specifically in the regulation of steroidogenic enzymes production. We studied the function of F in Blattella germanica oogenesis, as it shows two moments of high expression in ovaries: before the imaginal moult, and just before ovulation in the adult. Injecting dsftz-f1 into adult females, either just after the imaginal moult or just prior to choriogenesis, prevented oviposition, with differences between the two approaches. In 3-day-old adult females treated with dsftz-f1 just after the emergence, the expression of ftz-f1 was not modified, but the steroidogenic genes increased their expression. ftz-f1 transcript levels in the ovaries of 5-day-old dsftz-f1-treated females were significantly depleted, and the expression levels of the same steroidogenic genes began to decrease. These results suggest that Ftz-f1 regulates the expression of steroidogenic genes in B. germanica, with phm possibly being a key target. Ftz-f1 has a different temporal function in the cytoskeleton of follicular cells of the basal ovarian follicles. Early in the gonadotrophic cycle, Ftz-f1 promotes the expression of genes related to the cytoskeleton and muscle proteins, while at the end of the cycle it maintains the expression levels of these genes, thus ensuring correct ovulation.

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.

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.

Genome-wide identification of candidate chemosensory receptors in the bean bug Riptortus pedestris (Hemiptera: Alydidae) and the functional verification of its odorant receptor co-receptor (Orco) in recognizing aggregation pheromone

A sophisticated and sensitive olfactory system plays a vital role in the survival and reproduction of insects. Chemosensory receptors are indispensable for the molecular recognition and discrimination of semiochemicals. Riptortus pedestris is a notorious pest of legume plants, resulting in yield losses and quality decreases in soybeans. It is well accepted that R. pedestris highly relies on its olfactory system in detecting aggregation pheromones, host volatiles, and pesticides; however, little research focused on its chemosensory receptors. In the present study, we identified 237 odorant receptors (ORs), 42 gustatory receptors (GRs), and 31 ionotropic receptors (IRs) from the reported genome of R. pedestris, and analyzed their phylogenetic relationship with other hemipteran species. Through the results of RNA-seq and real-time quantitative PCR (qRT-PCR), we found that RpedORs displayed different expression levels in the antennae of R. pedestris at different development stages. To further verify the function of odorant receptor co-receptor (Orco), an obligate and unique insect OR, we silenced RpedOrco by RNA interference (RNAi) method. The results showed that silencing RpedOrco could significantly impair the response to aggregation pheromone in R. pedestris, indicating that RpedOrco plays an essential role in odorant detection. Our results can provide the theoretical foundations for revealing the olfactory recognition mechanism of R. pedestris and help explore and develop novel olfactory-based agents against this pest.

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.

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.

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.

Meta-Analysis of Public RNA Sequencing Data Revealed Potential Key Genes Associated with Reproductive Division of Labor in Social Hymenoptera and Termites

Eusociality in insects has evolved independently many times. One of the most notable characteristics of eusociality is the reproductive division of labor. In social insects, the reproductive division of labor is accomplished by queens and workers. Transcriptome analyses of queens and workers have been conducted for various eusocial species. However, the genes that regulate the reproductive division of labor across all or multiple eusocial species have not yet been fully elucidated. Therefore, we conducted a meta-analysis using publicly available RNA-sequencing data from four major groups of social insects. In this meta-analysis, we collected 258 pairs (queen vs. worker) of RNA-sequencing data from 34 eusocial species. The meta-analysis identified a total of 20 genes that were differentially expressed in queens or workers. Out of these, 12 genes have not previously been reported to be involved in the reproductive division of labor. Functional annotation of these 20 genes in other organisms revealed that they could be regulators of behaviors and physiological states related to the reproductive division of labor. These 20 genes, revealed using massive datasets of numerous eusocial insects, may be key regulators of the reproductive division of labor.

The remoulding of dietary effects on the fecundity / longevity trade-off in a social insect

BACKGROUND

In many organisms increased reproductive effort is associated with a shortened life span. This trade-off is reflected in conserved molecular pathways that link nutrient-sensing with fecundity and longevity. Social insect queens apparently defy the fecundity / longevity trade-off as they are both, extremely long-lived and highly fecund. Here, we have examined the effects of a protein-enriched diet on these life-history traits and on tissue-specific gene expression in a termite species of low social complexity.

RESULTS

On a colony level, we did not observe reduced lifespan and increased fecundity, effects typically seen in solitary model organisms, after protein enrichment. Instead, on the individual level mortality was reduced in queens that consumed more of the protein-enriched diet - and partially also in workers - while fecundity seemed unaffected. Our transcriptome analyses supported our life-history results. Consistent with life span extension, the expression of IIS (insulin/insulin-like growth factor 1 signalling) components was reduced in fat bodies after protein enrichment. Interestingly, however, genes involved in reproductive physiology (e.g., vitellogenin) were largely unaffected in fat body and head transcriptomes.

CONCLUSION

These results suggest that IIS is decoupled from downstream fecundity-associated pathways, which can contribute to the remoulding of the fecundity/longevity trade-off in termites as compared to solitary insects.

Tracking Adaptive Pathways of Invasive Insects: Novel Insight from Genomics

Despite the huge human and economic costs of invasive insects, which are the main group of invasive species, their environmental impacts through various mechanisms remain inadequately explained in databases and much of the invasion biology literature. High-throughput sequencing technology, especially whole-genome sequencing, has been used as a powerful method to study the mechanisms through which insects achieve invasion. In this study, we reviewed whole-genome sequencing-based advances in revealing several important invasion mechanisms of invasive insects, including (1) the rapid genetic variation and evolution of invasive populations, (2) invasion history and dispersal paths, (3) rapid adaptation to different host plant ranges, (4) strong environmental adaptation, (5) the development of insecticide resistance, and (6) the synergistic damage caused by invasive insects and endosymbiotic bacteria. We also discussed prevention and control technologies based on whole-genome sequencing and their prospects.

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.