Hemimetabolous genomes reveal molecular basis of termite eusociality

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.

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.

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.

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.

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.

DNA sequencing in the classroom: complete genome sequence of two earwig (Dermaptera; Insecta) species

BACKGROUND

Despite representing the largest fraction of animal life, the number of insect species whose genome has been sequenced is barely in the hundreds. The order Dermaptera (the earwigs) suffers from a lack of genomic information despite its unique position as one of the basally derived insect groups and its importance in agroecosystems. As part of a national educational and outreach program in genomics, a plan was formulated to engage the participation of high school students in a genome sequencing project. Students from twelve schools across Chile were instructed to capture earwig specimens in their geographical area, to identify them and to provide material for genome sequencing to be carried out by themselves in their schools.

RESULTS

The school students collected specimens from two cosmopolitan earwig species: Euborellia annulipes (Fam. Anisolabididae) and Forficula auricularia (Fam. Forficulidae). Genomic DNA was extracted and, with the help of scientific teams that traveled to the schools, was sequenced using nanopore sequencers. The sequence data obtained for both species was assembled and annotated. We obtained genome sizes of 1.18 Gb (F. auricularia) and 0.94 Gb (E. annulipes) with the number of predicted protein coding genes being 31,800 and 40,000, respectively. Our analysis showed that we were able to capture a high percentage (≥ 93%) of conserved proteins indicating genomes that are useful for comparative and functional analysis. We were also able to characterize structural elements such as repetitive sequences and non-coding RNA genes. Finally, functional categories of genes that are overrepresented in each species suggest important differences in the process underlying the formation of germ cells, and modes of reproduction between them, features that are one of the distinguishing biological properties that characterize these two distant families of Dermaptera.

CONCLUSIONS

This work represents an unprecedented instance where the scientific and lay community have come together to collaborate in a genome sequencing project. The versatility and accessibility of nanopore sequencers was key to the success of the initiative. We were able to obtain full genome sequences of two important and widely distributed species of insects which had not been analyzed at this level previously. The data made available by the project should illuminate future studies on the Dermaptera.

What Are the Functional Roles of Piwi Proteins and piRNAs in Insects?

Research on Piwi proteins and piRNAs in insects has focused on three experimental models: oogenesis and spermatogenesis in Drosophila melanogaster, the antiviral response in Aedes mosquitoes and the molecular analysis of primary and secondary piRNA biogenesis in Bombyx mori-derived BmN4 cells. Significant unique and complementary information has been acquired and has led to a greater appreciation of the complexity of piRNA biogenesis and Piwi protein function. Studies performed in other insect species are emerging and promise to add to the current state of the art on the roles of piRNAs and Piwi proteins. Although the primary role of the piRNA pathway is genome defense against transposons, particularly in the germline, recent findings also indicate an expansion of its functions. In this review, an extensive overview is presented of the knowledge of the piRNA pathway that so far has accumulated in insects. Following a presentation of the three major models, data from other insects were also discussed. Finally, the mechanisms for the expansion of the function of the piRNA pathway from transposon control to gene regulation were considered.

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.

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 chromosome-level genome of Chinese praying mantis Tenodera sinensis (Mantodea: Mantidae) reveals its biology as a predator

BACKGROUND

The Chinese praying mantis, Tenodera sinensis (Saussure), is a carnivorous insect that preys on a variety of arthropods and small vertebrates, including pest species. Several studies have been conducted to understand its behavior and physiology. However, there is limited knowledge about the genetic information underlying its genome evolution, digestive demands, and predatory behaviors.

FINDINGS

Here we have assembled the chromosome-level genome of T. sinensis, representing the first sequenced genome of the family Mantidae, with a genome size of 2.54 Gb and scaffold N50 of 174.78 Mb. Our analyses revealed that 98.6% of BUSCO genes are present, resulting in a well-annotated assembly compared to other insect genomes, containing 25,022 genes. The reconstructed phylogenetic analysis showed the expected topology placing the praying mantis in an appropriate position. Analysis of transposon elements suggested the Gypsy/Dirs family, which belongs to long terminal repeat (LTR) transposons, may be a key factor resulting in the larger genome size. The genome shows expansions in several digestion and detoxification associated gene families, including trypsin and glycosyl hydrolase (GH) genes, ATP-binding cassette (ABC) transporter, and carboxylesterase (CarE), reflecting the possible genomic basis of digestive demands. Furthermore, we have found 1 ultraviolet-sensitive opsin and 2 long-wavelength-sensitive (LWS) opsins, emphasizing the core role of LWS opsins in regulating predatory behaviors.

CONCLUSIONS

The high-quality genome assembly of the praying mantis provides a valuable repository for studying the evolutionary patterns of the mantis genomes and the gene expression profiles of insect predators.

Egfr signaling promotes juvenile hormone biosynthesis in the German cockroach

BACKGROUND

In insects, an interplay between the activities of distinct hormones, such as juvenile hormone (JH) and 20-hydroxyecdysone (20E), regulates the progression through numerous life history hallmarks. As a crucial endocrine factor, JH is mainly synthesized in the corpora allata (CA) to regulate multiple physiological and developmental processes, including molting, metamorphosis, and reproduction. During the last century, significant progress has been achieved in elucidating the JH signal transduction pathway, while less progress has been made in dissecting the regulatory mechanism of JH biosynthesis. Previous work has shown that receptor tyrosine kinase (RTK) signaling regulates hormone biosynthesis in both insects and mammals. Here, we performed a systematic RNA interference (RNAi) screening to identify RTKs involved in regulating JH biosynthesis in the CA of adult Blattella germanica females.

RESULTS

We found that the epidermal growth factor receptor (Egfr) is required for promoting JH biosynthesis in the CA of adult females. The Egf ligands Vein and Spitz activate Egfr, followed by Ras/Raf/ERK signaling, and finally activation of the downstream transcription factor Pointed (Pnt). Importantly, Pnt induces the transcriptional expression of two key enzyme-encoding genes in the JH biosynthesis pathway: juvenile hormone acid methyltransferase (JHAMT) and methyl farnesoate epoxidase (CYP15A1). Dual-luciferase reporter assay shows that Pnt is able to activate a promoter region of Jhamt. In addition, electrophoretic mobility shift assay confirms that Pnt directly binds to the - 941~ - 886 nt region of the Jhamt promoter.

CONCLUSIONS

This study reveals the detailed molecular mechanism of Egfr signaling in promoting JH biosynthesis in the German cockroach, shedding light on the intricate regulation of JH biosynthesis during insect development.

Melanin synthesis genes BgTH and BgDdc affect body color and cuticle permeability in Blattella germanica

Melanin is involved in cuticle pigmentation and sclerotization of insects, which is critical for maintaining structural integrity and functional completeness of insect cuticle. The 2 key enzymes of tyrosine hydroxylase (TH) and dopa decarboxylase (DDC) predicted in melanin biosynthesis are usually conserved in insects. However, it is unclear whether their function is related to epidermal permeability. In this study, we identified and cloned the gene sequences of BgTH and BgDdc from Blattella germanica, and revealed that they both showed a high expression at the molting, and BgTH was abundant in the head and integument while BgDdc was expressed highest in the fat body. Using RNA interference (RNAi), we found that knockdown of BgTH caused molting obstacles in some cockroaches, with the survivors showing pale color and softer integuments, while knockdown of BgDdc was viable and generated an abnormal light brown body color. Desiccation assay showed that the dsBgTH-injected adults died earlier than control groups under a dry atmosphere, but dsBgDdc-injected cockroaches did not. In contrast, when dsRNA-treated cockroaches were reared under a high humidity condition, almost no cockroaches died in all treatments. Furthermore, with eosin Y staining assay, we found that BgTH-RNAi resulted in a higher cuticular permeability, and BgDdc-RNAi also caused slight dye penetration. These results demonstrate that BgTH and BgDdc function in body pigmentation and affect the waterproofing ability of the cuticle, and the reduction of cuticular permeability may be achieved through cuticle melanization.

Evolutionary endocrinology and the problem of Darwin's tangled bank

Like Darwin's tangled bank of biodiversity, the endocrine mechanisms that give rise to phenotypic diversity also exhibit nearly endless forms. This tangled bank of mechanistic diversity can prove problematic as we seek general principles on the role of endocrine mechanisms in phenotypic evolution. A key unresolved question is therefore: to what degree are specific endocrine mechanisms re-used to bring about replicated phenotypic evolution? Related areas of inquiry are booming in molecular ecology, but behavioral traits are underrepresented in this literature. Here, I leverage the rich comparative tradition in evolutionary endocrinology to evaluate whether and how certain mechanisms may be repeated hotspots of behavioral evolutionary change. At one extreme, mechanisms may be parallel, such that evolution repeatedly uses the same gene or pathway to arrive at multiple independent (or, convergent) origins of a particular behavioral trait. At the other extreme, the building blocks of behavior may be unique, such that outwardly similar phenotypes are generated via lineage-specific mechanisms. This review synthesizes existing case studies, phylogenetic analyses, and experimental evolutionary research on mechanistic parallelism in animal behavior. These examples show that the endocrine building blocks of behavior have some elements of parallelism across replicated evolutionary events. However, support for parallelism is variable among studies, at least some of which relates to the level of complexity at which we consider sameness (i.e. pathway vs. gene level). Moving forward, we need continued experimentation and better testing of neutral models to understand whether, how - and critically, why - mechanism A is used in one lineage and mechanism B is used in another. We also need continued growth of large-scale comparative analyses, especially those that can evaluate which endocrine parameters are more or less likely to undergo parallel evolution alongside specific behavioral traits. These efforts will ultimately deepen understanding of how and why hormone-mediated behaviors are constructed the way that they are.

Eusocial Transition in Blattodea: Transposable Elements and Shifts of Gene Expression

(1) Unravelling the molecular basis underlying major evolutionary transitions can shed light on how complex phenotypes arise. The evolution of eusociality, a major evolutionary transition, has been demonstrated to be accompanied by enhanced gene regulation. Numerous pieces of evidence suggest the major impact of transposon insertion on gene regulation and its role in adaptive evolution. Transposons have been shown to be play a role in gene duplication involved in the eusocial transition in termites. However, evidence of the molecular basis underlying the eusocial transition in Blattodea remains scarce. Could transposons have facilitated the eusocial transition in termites through shifts of gene expression? (2) Using available cockroach and termite genomes and transcriptomes, we investigated if transposons insert more frequently in genes with differential expression in queens and workers and if those genes could be linked to specific functions essential for eusocial transition. (3) The insertion rate of transposons differs among differentially expressed genes and displays opposite trends between termites and cockroaches. The functions of termite transposon-rich queen- and worker-biased genes are related to reproduction and ageing and behaviour and gene expression, respectively. (4) Our study provides further evidence on the role of transposons in the evolution of eusociality, potentially through shifts in gene expression.

Coevolution of Metabolic Pathways in Blattodea and Their Blattabacterium Endosymbionts, and Comparisons with Other Insect-Bacteria Symbioses

Many insects harbor bacterial endosymbionts that supply essential nutrients and enable their hosts to thrive on a nutritionally unbalanced diet. Comparisons of the genomes of endosymbionts and their insect hosts have revealed multiple cases of mutually-dependent metabolic pathways that require enzymes encoded in 2 genomes. Complementation of metabolic reactions at the pathway level has been described for hosts feeding on unbalanced diets, such as plant sap. However, the level of collaboration between symbionts and hosts that feed on more variable diets is largely unknown. In this study, we investigated amino acid and vitamin/cofactor biosynthetic pathways in Blattodea, which comprises cockroaches and termites, and their obligate endosymbiont Blattabacterium cuenoti (hereafter Blattabacterium). In contrast to other obligate symbiotic systems, we found no clear evidence of "collaborative pathways" for amino acid biosynthesis in the genomes of these taxa, with the exception of collaborative arginine biosynthesis in 2 taxa, Cryptocercus punctulatus and Mastotermes darwiniensis. Nevertheless, we found that several gaps specific to Blattabacterium in the folate biosynthetic pathway are likely to be complemented by their host. Comparisons with other insects revealed that, with the exception of the arginine biosynthetic pathway, collaborative pathways for essential amino acids are only observed in phloem-sap feeders. These results suggest that the host diet is an important driving factor of metabolic pathway evolution in obligate symbiotic systems. IMPORTANCE The long-term coevolution between insects and their obligate endosymbionts is accompanied by increasing levels of genome integration, sometimes to the point that metabolic pathways require enzymes encoded in two genomes, which we refer to as "collaborative pathways". To date, collaborative pathways have only been reported from sap-feeding insects. Here, we examined metabolic interactions between cockroaches, a group of detritivorous insects, and their obligate endosymbiont, Blattabacterium, and only found evidence of collaborative pathways for arginine biosynthesis. The rarity of collaborative pathways in cockroaches and Blattabacterium contrasts with their prevalence in insect hosts feeding on phloem-sap. Our results suggest that host diet is a factor affecting metabolic integration in obligate symbiotic systems.

An aggregation model of cockroaches with fast-or-slow motion dichotomy

We propose a mathematical model, namely a reaction-diffusion system, to describe social behaviour of cockroaches. An essential new aspect in our model is that the dispersion behaviour due to overcrowding effect is taken into account as a counterpart to commonly studied aggregation. This consideration leads to an intriguing new phenomenon which has not been observed in the literature. Namely, due to the competition between aggregation towards areas of higher concentration of pheromone and dispersion avoiding overcrowded areas, the cockroaches aggregate more at the transition area of pheromone. Moreover, we also consider the fast reaction limit where the switching rate between active and inactive subpopulations tends to infinity. By utilising improved duality and energy methods, together with the regularisation of heat operator, we prove that the weak solution of the reaction-diffusion system converges to that of a reaction-cross-diffusion system.