Hemimetabolous genomes reveal molecular basis of termite eusociality

Long-lived termite kings and queens activate telomerase in somatic organs

Kings and queens of termites, like queens of other advanced eusocial insects, are endowed with admirable longevity, which dramatically exceeds the life expectancies of their non-reproducing nest-mates and related solitary insects. In the quest to find the mechanisms underlying the longevity of termite reproductives, we focused on somatic maintenance mediated by telomerase. This ribonucleoprotein is well established for pro-longevity functions in vertebrates, thanks primarily to its ability of telomere extension. However, its participation in lifespan regulation of insects, including the eusocial taxa, remains understudied. Here, we report a conspicuous increase of telomerase abundance and catalytic activity in the somatic organs of primary and secondary reproductives of the termite Prorhinotermes simplex and confirm a similar pattern in two other termite species. These observations stand in contrast with the telomerase downregulation characteristic for most adult somatic tissues in vertebrates and also in solitary insects and non-reproducing castes of termites. At the same time, we did not observe caste-specific differences in telomere lengths that might explain the differential longevity of termite castes. We conclude that although the telomerase activation in termite reproductives is in line with the broadly assumed association between telomerase and longevity, its direct phenotypic impact remains to be elucidated.

Eusociality Shapes Convergent Patterns of Molecular Evolution across Mitochondrial Genomes of Snapping Shrimps

Eusociality is a highly conspicuous and ecologically impactful behavioral syndrome that has evolved independently across multiple animal lineages. So far, comparative genomic analyses of advanced sociality have been mostly limited to insects. Here, we study the only clade of animals known to exhibit eusociality in the marine realm-lineages of socially diverse snapping shrimps in the genus Synalpheus. To investigate the molecular impact of sociality, we assembled the mitochondrial genomes of eight Synalpheus species that represent three independent origins of eusociality and analyzed patterns of molecular evolution in protein-coding genes. Synonymous substitution rates are lower and potential signals of relaxed purifying selection are higher in eusocial relative to noneusocial taxa. Our results suggest that mitochondrial genome evolution was shaped by eusociality-linked traits-extended generation times and reduced effective population sizes that are hallmarks of advanced animal societies. This is the first direct evidence of eusociality impacting genome evolution in marine taxa. Our results also strongly support the idea that eusociality can shape genome evolution through profound changes in life history and demography.

Expression profiles of neotropical termites reveal microbiota-associated, caste-biased genes and biotechnological targets

Termites are well recognized by their complex development trajectories, involving dynamic differentiation process between non-reproductive castes, workers and soldiers. These insects are associated with endosymbiotic microorganisms, which help in lignocellulose digestion and nitrogen metabolism. Aiming to identify genes harbouring biotechnological potential, we analyzed workers and soldiers RNA-Seq data of three neotropical termites: Heterotermes tenuis (Isoptera: Rhinotermitidae), Velocitermes heteropterus (Isoptera: Termitidae) and Cornitermes cumulans (Isoptera: Termitidae). We observed differences in the microbiota associated with each termite family, and found protists' genes in both Termitidae species. We found an opposite pattern of caste-biased gene expression between H. tenuis and the termitids studied. Moreover, the two termitids are considerably different concerning the number of differentially expressed genes (DEGs). Functional annotation indicated considerable differences in caste-biased gene content between V. heteropterus and C. cumulans, even though they share similar diet and biological niche. Among the most DEGs, we highlighted those involved in caste differentiation and cellulose digestion, which are attractive targets for studying more efficient technologies for termite control, biomass digestion and other biotechnological applications.

Genetic basis of chemical communication in eusocial insects

Social behavior is one of the most fascinating and complex behaviors in humans and animals. A fundamental process of social behavior is communication among individuals. It relies on the capability of the nervous system to sense, process, and interpret various signals (e.g., pheromones) and respond with appropriate decisions and actions. Eusocial insects, including ants, some bees, some wasps, and termites, display intriguing cooperative social behavior. Recent advances in genetic and genomic studies have revealed key genes that are involved in pheromone synthesis, chemosensory perception, and physiological and behavioral responses to varied pheromones. In this review, we highlight the genes and pathways that regulate queen pheromone-mediated social communication, discuss the evolutionary changes in genetic systems, and outline prospects of functional studies in sociobiology.

Phylogenetic and sequence analyses of insect transferrins suggest that only transferrin 1 has a role in iron homeostasis

Iron is essential to life, but surprisingly little is known about how iron is managed in nonvertebrate animals. In mammals, the well-characterized transferrins bind iron and are involved in iron transport or immunity, whereas other members of the transferrin family do not have a role in iron homeostasis. In insects, the functions of transferrins are still poorly understood. The goals of this project were to identify the transferrin genes in a diverse set of insect species, resolve the evolutionary relationships among these genes, and predict which of the transferrins are likely to have a role in iron homeostasis. Our phylogenetic analysis of transferrins from 16 orders of insects and two orders of noninsect hexapods demonstrated that there are four orthologous groups of insect transferrins. Our analysis suggests that transferrin 2 arose prior to the origin of insects, and transferrins 1, 3, and 4 arose early in insect evolution. Primary sequence analysis of each of the insect transferrins was used to predict signal peptides, carboxyl-terminal transmembrane regions, GPI-anchors, and iron binding. Based on this analysis, we suggest that transferrins 2, 3, and 4 are unlikely to play a major role in iron homeostasis. In contrast, the transferrin 1 orthologs are predicted to be secreted, soluble, iron-binding proteins. We conclude that transferrin 1 orthologs are the most likely to play an important role in iron homeostasis. Interestingly, it appears that the louse, aphid, and thrips lineages have lost the transferrin 1 gene and, thus, have evolved to manage iron without transferrins.

Complex Evolution of Insect Insulin Receptors and Homologous Decoy Receptors, and Functional Significance of Their Multiplicity

Evidence accumulates that the functional plasticity of insulin and insulin-like growth factor signaling in insects could spring, among others, from the multiplicity of insulin receptors (InRs). Their multiple variants may be implemented in the control of insect polyphenism, such as wing or caste polyphenism. Here, we present a comprehensive phylogenetic analysis of insect InR sequences in 118 species from 23 orders and investigate the role of three InRs identified in the linden bug, Pyrrhocoris apterus, in wing polymorphism control. We identified two gene clusters (Clusters I and II) resulting from an ancestral duplication in a late ancestor of winged insects, which remained conserved in most lineages, only in some of them being subject to further duplications or losses. One remarkable yet neglected feature of InR evolution is the loss of the tyrosine kinase catalytic domain, giving rise to decoys of InR in both clusters. Within the Cluster I, we confirmed the presence of the secreted decoy of insulin receptor in all studied Muscomorpha. More importantly, we described a new tyrosine kinase-less gene (DR2) in the Cluster II, conserved in apical Holometabola for ∼300 My. We differentially silenced the three P. apterus InRs and confirmed their participation in wing polymorphism control. We observed a pattern of Cluster I and Cluster II InRs impact on wing development, which differed from that postulated in planthoppers, suggesting an independent establishment of insulin/insulin-like growth factor signaling control over wing development, leading to idiosyncrasies in the co-option of multiple InRs in polyphenism control in different taxa.

Silencing the odorant coreceptor (Orco) disrupts sex pheromonal communication and feeding responses in Blattella germanica: toward an alternative target for controlling insect-transmitted human diseases

BACKGROUND

The German cockroach, Blattella germanica, is one of the most severe pests of urban and rural areas. High-throughput genetic screening approaches indicate that the olfactory system of this pest is extremely powerful because it has an extensive array of olfactory receptor genes compared with many other insect species. Several of these genes have been identified previously, but their functions have not yet been characterized.

RESULTS

This study describes the sequence of five transcriptomes of B. germanica adult male antennae, female antennae, maxillary palps, legs, and fifth-instar nymph antennae to investigate expression patterns of odorant receptors (ORs). Approximately 90% of ORs were found to be the most highly expressed genes in adult or nymph antennae. Additionally, every OR requires an odorant co-receptor (Orco) to become fully functional, and this was selected and successfully inhibited by injection of the corresponding double-stranded (ds)RNA targeting the Orco. A strong RNA interference (RNAi) effect was observed in which > 75% of Orco messenger RNA (mRNA) was clearly suppressed after 72 h of treatment. Olfactory behavioral assays showed that Orco-impaired B. germanica respond more slowly and show less attraction to one volatile sex pheromone and food resources compared with a control group.

CONCLUSION

The results show that Orco plays a pivotal role in both sex pheromone and food-seeking olfactory processes, and provide an alternative genetic technique for controlling this urban pest species by olfactory disruption. © 2020 Society of Chemical Industry.

Transcriptomic analyses of the termite, Cryptotermes secundus, reveal a gene network underlying a long lifespan and high fecundity

Organisms are typically characterized by a trade-off between fecundity and longevity. Notable exceptions are social insects. In insect colonies, the reproducing caste (queens) outlive their non-reproducing nestmate workers by orders of magnitude and realize fecundities and lifespans unparalleled among insects. How this is achieved is not understood. Here, we identified a single module of co-expressed genes that characterized queens in the termite species Cryptotermes secundus. It encompassed genes from all essential pathways known to be involved in life-history regulation in solitary model organisms. By manipulating its endocrine component, we tested the recent hypothesis that re-wiring along the nutrient-sensing/endocrine/fecundity axis can account for the reversal of the fecundity/longevity trade-off in social insect queens. Our data from termites do not support this hypothesis. However, they revealed striking links to social communication that offer new avenues to understand the re-modelling of the fecundity/longevity trade-off in social insects.

How Can Termites Achieve Their Unparalleled Postembryonic Developmental Plasticity? A Test for the Role of Intermolt-Specific High Juvenile Hormone Titers

Termites are “social cockroaches” and amongst the most phenotypically plastic insects. The different castes (i.e., two types of reproductives, workers, and soldiers) within termite societies are all encoded by a single genome and present the result of differential postembryonic development. Besides the default progressive development into winged sexuals of solitary hemimetabolous insects, termites have two postembryonic, non-terminal molts (stationary and regressive; i.e., molts associated, respectively, with no change or reduction of size/morphological differentiation) which allow them to retain workers, and two terminal developmental types to become soldiers and replacement reproductives. Despite this unique plasticity, especially the mechanisms underlying the non-terminal development are poorly understood. In 1982, Nijhout and Wheeler proposed a model how this diversity might have evolved. They proposed that varying juvenile hormone (JH) titers at the start, mid-phase, and end of each intermolt period account for the developmental diversity. We tested this rarely addressed model in the lower termite Cryptotermes secundus using phase-specific pharmacological manipulations of JH titers. Our results partially support the Nijhout and Wheeler model. These data are supplemented with gene expression studies of JH-related genes that characterize different postembryonic developmental trajectories. Our study provides new insights into the evolution of the unique postembryonic developmental plasticity of termites that constitutes the foundation of their social life.

The Structure and Function of Ionotropic Receptors in Drosophila

Ionotropic receptors (IRs) are a highly divergent subfamily of ionotropic glutamate receptors (iGluR) and are conserved across Protostomia, a major branch of the animal kingdom that encompasses both Ecdysozoa and Lophothrochozoa. They are broadly expressed in peripheral sensory systems, concentrated in sensory dendrites, and function in chemosensation, thermosensation, and hygrosensation. As iGluRs, four IR subunits form a functional ion channel to detect environmental stimuli. Most IR receptors comprise individual stimulus-specific tuning receptors and one or two broadly expressed coreceptors. This review summarizes the discoveries of the structure of IR complexes and the expression and function of each IR, as well as discusses the future direction for IR studies.

A Review of Termite Pheromones: Multifaceted, Context-Dependent, and Rational Chemical Communications

Termite colonies, composed of large numbers of siblings, develop an important caste-based division of labor; individuals in these societies interact via intra- or intercaste chemical communications. For more than 50 years, termites have been known to use a variety of pheromones to perform tasks necessary for maintenance of their societies, similar to eusocial hymenopterans. Although trail-following pheromones have been chemically identified in various termites, other types of pheromones have not been elucidated chemically or functionally. In the past decade, however, chemical compositions and biological functions have been successfully identified for several types of termite pheromones; accordingly, the details of the underlying pheromone communications have been gradually revealed. In this review, we summarize both the functions of all termite pheromones identified so far and the chemical interactions among termites and other organisms. Subsequently, we argue how termites developed their sophisticated pheromone communication. We hypothesize that termites have diverted defensive and antimicrobial substances to pheromones associated in caste recognition and caste-specific roles. Furthermore, termites have repeatedly used a pre-existing pheromone or have added supplementary compounds to it in accordance with the social context, leading to multifunctionalization of pre-existing pheromones and emergence of new pheromones. These two mechanisms may enable termites to transmit various context-dependent information with a small number of chemicals, thus resulting in formation of coordinated, complex, and rational chemical communication systems.

Transposable Elements and the Evolution of Insects

Insects are major contributors to our understanding of the interaction between transposable elements (TEs) and their hosts, owing to seminal discoveries, as well as to the growing number of sequenced insect genomes and population genomics and functional studies. Insect TE landscapes are highly variable both within and across insect orders, although phylogenetic relatedness appears to correlate with similarity in insect TE content. This correlation is unlikely to be solely due to inheritance of TEs from shared ancestors and may partly reflect preferential horizontal transfer of TEs between closely related species. The influence of insect traits on TE landscapes, however, remains unclear. Recent findings indicate that, in addition to being involved in insect adaptations and aging, TEs are seemingly at the cornerstone of insect antiviral immunity. Thus, TEs are emerging as essential insect symbionts that may have deleterious or beneficial consequences on their hosts, depending on context.

A Century of Synergy in Termite Symbiosis Research: Linking the Past with New Genomic Insights

Termites have long been studied for their symbiotic associations with gut microbes. In the late nineteenth century, this relationship was poorly understood and captured the interest of parasitologists such as Joseph Leidy; this research led to that of twentieth-century biologists and entomologists including Cleveland, Hungate, Trager, and Lüscher. Early insights came via microscopy, organismal, and defaunation studies, which led to descriptions of microbes present, descriptions of the roles of symbionts in lignocellulose digestion, and early insights into energy gas utilization by the host termite. Focus then progressed to culture-dependent microbiology and biochemical studies of host-symbiont complementarity, which revealed specific microhabitat requirements for symbionts and noncellulosic mechanisms of symbiosis (e.g., N₂ fixation). Today, knowledge on termite symbiosis has accrued exponentially thanks to omic technologies that reveal symbiont identities, functions, and interdependence, as well as intricacies of host-symbiont complementarity. Moving forward, the merging of classical twentieth-century approaches with evolving omic tools should provide even deeper insights into host-symbiont interplay.

No obvious transcriptome-wide signature of indirect selection in termites

The evolution of sterile helper castes in social insects implies selection on genes that underlie variation in this nonreproductive phenotype. These focal genes confer no direct fitness and are presumed to evolve through indirect fitness effects on the helper's reproducing relatives. This separation of a gene's phenotypic effect on one caste and its fitness effect on another suggests that genes for this and other forms of reproductive altruism are buffered from selection and will thus evolve closer to the neutral rate than genes directly selected for selfish reproduction. We test this hypothesis by comparing the strength of selection at loci associated in their expression with reproductive versus sterile castes in termites. Specifically, we gather caste-biased gene expression data from four termite transcriptomes and measure the global dN/dS ratio across gene sets and phylogenetic lineages. We find that the majority of examined orthologous gene groups show patterns of nucleotide substitution that are consistent with strong purifying selection and display little evidence for distinct signatures of direct versus indirect selection in reproductive and sterile castes. For one particular species (Reticulitermes flavipes), the strength of purifying selection is relaxed in a reproductive nymph-biased gene set, which opposes the nearly neutral idea. In other species, the synonymous rate (dS) alone was often found to be the highest in the sterile worker caste, suggesting a more subtle signature of indirect selection or an altogether different relationship between caste-biased expression and rates of molecular evolution.

DNMT1 Promotes Genome Methylation and Early Embryo Development in Cockroaches

The influence of DNA methylation on gene behavior and its consequent phenotypic effects appear to be very important, but the details are not well understood. Insects offer a diversity of DNA methylation modes, making them an excellent lineage for comparative analyses. However, functional studies have tended to focus on quite specialized holometabolan species, such as wasps, bees, beetles, and flies. Here, we have studied DNA methylation in the hemimetabolan insect Blattella germanica. In this cockroach, a gene involved in DNA methylation, DNA methyltransferase 1 (DNMT1), is expressed in early embryogenesis. In our experiments, RNAi of DNMT1 reduces DNA methylation and impairs blastoderm formation. Using reduced representation bisulfite sequencing and transcriptome analyses, we observed that methylated genes are associated with metabolism and are highly expressed, whereas unmethylated genes are related to signaling and show low expression. Moreover, methylated genes show greater expression change and less expression variability than unmethylated genes.

Blattella germanica displays a large arsenal of antimicrobial peptide genes

Defence systems against microbial pathogens are present in most living beings. The German cockroach Blattella germanica requires these systems to adapt to unhealthy environments with abundance of pathogenic microbes, in addition to potentially control its symbiotic systems. To handle this situation, four antimicrobial gene families (defensins, termicins, drosomycins and attacins) were expanded in its genome. Remarkably, a new gene family (blattellicins) emerged recently after duplication and fast evolution of an attacin gene, which is now encoding larger proteins with the presence of a long stretch of glutamines and glutamic acids. Phylogenetic reconstruction, within Blattellinae, suggests that this duplication took place before the divergence of Blattella and Episymploce genera. The latter harbours a long attacin gene (pre-blattellicin), but the absence of the encoded Glx-region suggests that this element evolved recently in the Blattella lineage. A screening of AMP gene expression in available transcriptomic SR projects of B. germanica showed that, while some AMPs are expressed during almost the whole development, others are restricted to shorter periods. Blattellicins are highly expressed only in adult females. None of the available SR tissue projects could be associated with blattellicins’ expression, suggesting that it takes place in other tissues, maybe the gut.

Expression, Affinity, and Functional Characterization of the Specific Binding of Two Putative Pheromone-Binding Proteins in the Omnivorous German Cockroach Blattella germanica

The German cockroach Blattella germanica (L.) is an important pest in medical, veterinary, and public health. Studies on the olfaction mechanism of hemimetabolous insects have rarely been reported, especially in cockroaches. Pheromone-binding proteins (PBPs) play a vital role in insect sex pheromone recognition, which solubilize and carry the hydrophobic pheromonal compounds through the antennal lymph to receptors. In this study, two potential PBPs (BgerOBP26 and BgerOBP40) were identified on the basis of their biased expression in male antennae using tissue transcriptome data and verified by the quantitative real-time polymerase chain reaction approach. We then expressed and purified the two identified odorant-binding proteins (OBPs) using the Escherichia coli expression system and affinity purification. In vitro binding studies showed that the two OBPs display stronger binding affinities to the female volatile sex pheromone blattellaquinone than to its analogues and contact sex pheromone components. Finally, three-dimensional modeling of the two OBPs and dock conformation with sex pheromone molecules showed BgerOBP26 has a larger odorant cavity and more conservative active amino acid residues than BgerOBP40. These results illuminated the binding characteristics of potential PBPs of B. germanica, which could lay the groundwork for improved understanding of many aspects of the chemical ecology of B. germanica. Moreover, this information complements the understanding of the olfactory molecular mechanism in cockroaches and provides potential gene targets for B. germanica control.

On the roles of AA15 lytic polysaccharide monooxygenases derived from the termite Coptotermes gestroi

Lytic polysaccharide monooxygenases (LPMOs) are copper-dependent enzymes which catalyze the oxidative cleavage of polysaccharides. LPMOs belonging to family 15 in the Auxiliary Activity (AA) class from the Carbohydrate-Active Enzyme database are found widespread across the Tree of Life, including viruses, algae, oomycetes and animals. Recently, two AA15s from the firebrat Thermobia domestica were reported to have oxidative activity, one towards cellulose or chitin and the other towards chitin, signalling that AA15 LPMOs from insects potentially have different biochemical functions. Herein, we report the identification and characterization of two family AA15 members from the lower termite Coptotermes gestroi. Addition of Cu(II) to CgAA15a or CgAA15b had a thermostabilizing effect on both. Using ascorbate and O₂ as co-substrates, CgAA15a and CgAA15b were able to oxidize chitin, but showed no activity on celluloses, xylan, xyloglucan and starch. Structural models indicate that the LPMOs from C. gestroi (CgAA15a/CgAA15b) have a similar fold but exhibit key differences in the catalytic site residues when compared to the cellulose/chitin-active LPMO from T. domestica (TdAA15a), especially the presence of a non-coordinating phenylalanine nearby the Cu ion in CgAA15a/b, which appears as a tyrosine in the active site of TdAA15a. Despite the overall similarity in protein folds, however, mutation of the active site phenylalanine in CgAA15a to a tyrosine did not expanded the enzymatic specificity from chitin to cellulose. Our data show that CgAA15a/b enzymes are likely not involved in lignocellulose digestion but might play a role in termite developmental processes as well as on chitin and nitrogen metabolisms.

R2 and Non-Site-Specific R2-Like Retrotransposons of the German Cockroach, Blattella germanica

The structural and functional organization of the ribosomal RNA gene cluster and the full-length R2 non-LTR retrotransposon (integrated into a specific site of 28S ribosomal RNA genes) of the German cockroach, Blattella germanica, is described. A partial sequence of the R2 retrotransposon of the cockroach Rhyparobia maderae is also analyzed. The analysis of previously published next-generation sequencing data from the B. germanica genome reveals a new type of retrotransposon closely related to R2 retrotransposons but with a random distribution in the genome. Phylogenetic analysis reveals that these newly described retrotransposons form a separate clade. It is shown that proteins corresponding to the open reading frames of newly described retrotransposons exhibit unequal structural domains. Within these retrotransposons, a recombination event is described. New mechanism of transposition activity is discussed. The essential structural features of R2 retrotransposons are conserved in cockroaches and are typical of previously described R2 retrotransposons. However, the investigation of the number and frequency of 5′-truncated R2 retrotransposon insertion variants in eight B. germanica populations suggests recent mobile element activity. It is shown that the pattern of 5′-truncated R2 retrotransposon copies can be an informative molecular genetic marker for revealing genetic distances between insect populations.

SMRT sequencing of the full-length transcriptome of Odontotermes formosanus (Shiraki) under Serratia marcescens treatment

Odontotermes formosanus (Shiraki) is an important pest in the world. Serratia marcescens have a high lethal effect on O. formosanus, but the specific insecticidal mechanisms of S. marcescens on O. formosanus are unclear, and the immune responses of O. formosanus to S. marcescens have not been clarified. At present, genetic database resources of O. formosanus are extremely scarce. Therefore, using O. formosanus workers infected by S. marcescens and the control as experimental materials, a full-length transcriptome was sequenced using the PacBio Sequel sequencing platform. A total of 10,364 isoforms were obtained as the final transcriptome. The unigenes were further annotated with the Nr, Swiss-Prot, EuKaryotic Orthologous Groups (KOG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) Ortholog public databases. In a comparison between the control group and a Serratia marcescens-infected group, a total of 259 differentially expressed genes (DEGs) were identified, including 132 upregulated and 127 downregulated genes. Pathway enrichment analysis indicated that the expression of the mitogen-activated protein kinase (MAPK) pathway, oxidative stress genes and the AMP-activated protein kinase (AMPK) pathway in O. formosanus may be associated with S. marcescens treatment. This research intensively studied O. formosanus at the high-throughput full-length transcriptome level, laying a foundation for further development of molecular markers and mining of target genes in this species and thereby promoting the biological control of O. formosanus. Furthermore, these results will be helpful to clarify the action mechanisms of S. marcescens on O. formosanus, and also explore the relationship between O. formosanus and S. marcescens. In addition, this study will identify the immune response of O. formosanus to S. marcescens, which will provide a theoretical foundation for the development of new immunosuppressants for O. formosanus.

gNOMO: a multi-omics pipeline for integrated host and microbiome analysis of non-model organisms

The study of bacterial symbioses has grown exponentially in the recent past. However, existing bioinformatic workflows of microbiome data analysis do commonly not integrate multiple meta-omics levels and are mainly geared toward human microbiomes. Microbiota are better understood when analyzed in their biological context; that is together with their host or environment. Nevertheless, this is a limitation when studying non-model organisms mainly due to the lack of well-annotated sequence references. Here, we present gNOMO, a bioinformatic pipeline that is specifically designed to process and analyze non-model organism samples of up to three meta-omics levels: metagenomics, metatranscriptomics and metaproteomics in an integrative manner. The pipeline has been developed using the workflow management framework Snakemake in order to obtain an automated and reproducible pipeline. Using experimental datasets of the German cockroach Blattella germanica, a non-model organism with very complex gut microbiome, we show the capabilities of gNOMO with regard to meta-omics data integration, expression ratio comparison, taxonomic and functional analysis as well as intuitive output visualization. In conclusion, gNOMO is a bioinformatic pipeline that can easily be configured, for integrating and analyzing multiple meta-omics data types and for producing output visualizations, specifically designed for integrating paired-end sequencing data with mass spectrometry from non-model organisms.

Arthropod IGF, relaxin and gonadulin, putative orthologs of Drosophila insulin-like peptides 6, 7 and 8, likely originated from an ancient gene triplication

BACKGROUND

Insects have several genes coding for insulin-like peptides and they have been particularly well studied in Drosophila. Some of these hormones function as growth hormones and are produced by the fat body and the brain. These act through a typical insulin receptor tyrosine kinase. Two other Drosophila insulin-like hormones are either known or suspected to act through a G-protein coupled receptor. Although insulin-related peptides are known from other insect species, Drosophila insulin-like peptide 8, one that uses a G-protein coupled receptor, has so far only been identified from Drosophila and other flies. However, its receptor is widespread within arthropods and hence it should have orthologs. Such putative orthologs were recently identified in decapods and have been called gonadulins.

METHODOLOGY

In an effort to identify gonadulins in other arthropods public genome assemblies and short-read archives from insects and other arthropods were explored for the presence of genes and transcripts coding insulin-like peptides and their putative receptors.

RESULTS

Gonadulins were detected in a number of arthropods. In those species for which transcriptome data from the gonads is available insect gonadulin genes are expressed in the ovaries and at least in some species also in the testes. In some insects differences in gonadulin expression in the ovary between actively reproducing and non-reproducing females differs more than 100-fold. Putative orthologs of Drosophila ilp 6 were also identified. In several non-Dipteran insects these peptides have C-terminally extensions that are alternatively spliced. The predicted peptides have been called arthropod insulin-like growth factors. In cockroaches, termites and stick insects genes coding for the arthropod insulin-like growth factors, gonadulin and relaxin, a third insulin-like peptide, are encoded by genes that are next to one another suggesting that they are the result of a local gene triplication. Such a close chromosomal association was also found for the arthropod insulin-like growth factor and gonadulin genes in spiders. Phylogenetic tree analysis of the typical insulin receptor tyrosine kinases from insects, decapods and chelicerates shows that the insulin signaling pathway evolved differently in these three groups. The G-protein coupled receptors that are related to the Drosophila ilp 8 receptor similarly show significant differences between those groups.

CONCLUSION

A local gene triplication in an early ancestor likely yielded three genes coding gonadulin, arthropod insulin-like growth factor and relaxin. Orthologs of these genes are now commonly present in arthropods and almost certainly include the Drosophila insulin-like peptides 6, 7 and 8.

Nutritional symbiosis and ecology of host-gut microbe systems in the Blattodea

Cockroaches and termites (Order: Blattodea) have been the subject of substantial research attention for over a century due, in part, to a subset of them having a strong propensity to cohabitate with humans and their structures. Recent research has led to numerous insights into their behavior, physiology, and ecology, as well as their ability to harbor taxonomically diverse microbial communities within their digestive systems, which include taxa that contribute to host growth and development. Further, recent investigations into the physiological and behavioral adaptations that enable recalcitrant polysaccharide digestion and the maintenance of microbial symbionts in cockroaches and termites suggests that symbionts contribute significantly to nutrient provisioning and processing.

Functional insights from the GC-poor genomes of two aphid parasitoids, Aphidius ervi and Lysiphlebus fabarum

BACKGROUND

Parasitoid wasps have fascinating life cycles and play an important role in trophic networks, yet little is known about their genome content and function. Parasitoids that infect aphids are an important group with the potential for biological control. Their success depends on adapting to develop inside aphids and overcoming both host aphid defenses and their protective endosymbionts.

RESULTS

We present the de novo genome assemblies, detailed annotation, and comparative analysis of two closely related parasitoid wasps that target pest aphids: Aphidius ervi and Lysiphlebus fabarum (Hymenoptera: Braconidae: Aphidiinae). The genomes are small (139 and 141 Mbp) and the most AT-rich reported thus far for any arthropod (GC content: 25.8 and 23.8%). This nucleotide bias is accompanied by skewed codon usage and is stronger in genes with adult-biased expression. AT-richness may be the consequence of reduced genome size, a near absence of DNA methylation, and energy efficiency. We identify missing desaturase genes, whose absence may underlie mimicry in the cuticular hydrocarbon profile of L. fabarum. We highlight key gene groups including those underlying venom composition, chemosensory perception, and sex determination, as well as potential losses in immune pathway genes.

CONCLUSIONS

These findings are of fundamental interest for insect evolution and biological control applications. They provide a strong foundation for further functional studies into coevolution between parasitoids and their hosts. Both genomes are available at https://bipaa.genouest.org.

Multi-faceted analysis provides little evidence for recurrent whole-genome duplications during hexapod evolution

BACKGROUND

Gene duplication events play an important role in the evolution and adaptation of organisms. Duplicated genes can arise through different mechanisms, including whole-genome duplications (WGDs). Recently, WGD was suggested to be an important driver of evolution, also in hexapod animals.

RESULTS

Here, we analyzed 20 high-quality hexapod genomes using whole-paranome distributions of estimated synonymous distances (KS), patterns of within-genome co-linearity, and phylogenomic gene tree-species tree reconciliation methods. We observe an abundance of gene duplicates in the majority of these hexapod genomes, yet we find little evidence for WGD. The majority of gene duplicates seem to have originated through small-scale gene duplication processes. We did detect segmental duplications in six genomes, but these lacked the within-genome co-linearity signature typically associated with WGD, and the age of these duplications did not coincide with particular peaks in KS distributions. Furthermore, statistical gene tree-species tree reconciliation failed to support all but one of the previously hypothesized WGDs.

CONCLUSIONS

Our analyses therefore provide very limited evidence for WGD having played a significant role in the evolution of hexapods and suggest that alternative mechanisms drive gene duplication events in this group of animals. For instance, we propose that, along with small-scale gene duplication events, episodes of increased transposable element activity could have been an important source for gene duplicates in hexapods.

Genomic adaptations to aquatic and aerial life in mayflies and the origin of insect wings

The evolution of winged insects revolutionized terrestrial ecosystems and led to the largest animal radiation on Earth. However, we still have an incomplete picture of the genomic changes that underlay this diversification. Mayflies, as one of the sister groups of all other winged insects, are key to understanding this radiation. Here, we describe the genome of the mayfly Cloeon dipterum and its gene expression throughout its aquatic and aerial life cycle and specific organs. We discover an expansion of odorant-binding-protein genes, some expressed specifically in breathing gills of aquatic nymphs, suggesting a novel sensory role for this organ. In contrast, flying adults use an enlarged opsin set in a sexually dimorphic manner, with some expressed only in males. Finally, we identify a set of wing-associated genes deeply conserved in the pterygote insects and find transcriptomic similarities between gills and wings, suggesting a common genetic program. Globally, this comprehensive genomic and transcriptomic study uncovers the genetic basis of key evolutionary adaptations in mayflies and winged insects.

Molecular mechanisms underlying milk production and viviparity in the cockroach, Diploptera punctata

Viviparous reproduction is characterized by maternal retention of developing offspring within the reproductive tract during gestation, culminating in live birth. In some cases, a mother will provide nutrition beyond that present in the yolk; this is known as matrotrophic viviparity. While this phenomenon is best associated with mammals, it is observed in insects such as the viviparous cockroach, Diploptera punctata. Female D. punctata carry developing embryos in the brood sac, a reproductive organ that acts as both a uterus and a placenta by protecting and providing a nutritive secretion to the intrauterine developing progeny. While the basic physiology of D. punctata pregnancy has been characterized, little is known about the molecular mechanisms underlying this phenomenon. This study combined RNA-seq analysis, RNA interference, and other assays to characterize molecular and physiological changes associated with D. punctata reproduction. A comparison of four stages of the female reproductive cycle and males revealed unique gene expression profiles corresponding to each stage and between sexes. Differentially regulated transcripts of interest include the previously identified family of milk proteins and transcripts associated with juvenile hormone metabolism. RNA interference and methoprene application experiments established the potential impacts of bothbreakdown and synthesis reduction of juvenile hormone in maintaining pregnancy in D. punctata. These studies provide the comprehensive molecular mechanisms associated with cockroach viviparity, which will be a critical resource for comparative purposes among viviparity in insect systems.

The making of the defensive caste: Physiology, development, and evolution of the soldier differentiation in termites

Termites (Blattodea, Termitoidea, or Isoptera) constitute one of the major lineages of eusocial insects. In termite societies, multiple types of functional individuals, that is, castes, perform divisions of labors to coordinate social behaviors. Among other castes, the soldier caste is distinctive since it is sterile and exclusively specialized into defensive behavior with largely modified morphological features. Therefore, many of the previous studies have been focused on soldiers, in terms of ecology, behavior, and evolution as well as developmental and physiological mechanisms. This article overviews the accumulation of studies especially focusing on the developmental and physiological mechanisms underlying the soldier differentiation in termites. Furthermore, the evolutionary trajectories that have led the acquisition of soldier caste and have diversified the soldier characteristics in association with the social evolution are discussed.

Transcriptomic analysis of epigenetic modification genes in the termite Reticulitermes speratus

Eusocial insects display a caste system in which different castes are morphologically and physiologically specialized for different tasks. Recent studies have revealed that epigenetic modifications, including DNA methylation and histone modification, mediate caste determination and differentiation, longevity, and polyethism in eusocial insects. Although there has been a growing interest in the relationship between epigenetic mechanisms and phenotypic plasticity in termites, there is little information about differential expression levels among castes and expression sites for these genes in termites. Here we show royal-tissue-specific expression of epigenetic modification genes in the termite Reticulitermes speratus. Using RNA-seq, we identified 74 genes, including three DNA methyltransferases, seven sirtuins, 48 Trithorax group proteins, and 16 Polycomb group proteins. Among these genes, 15 showed king-specific expression, and 52 showed age-dependent differential expression in kings and queens. Quantitative real-time PCR revealed that DNA methyltransferase 3 is expressed specifically in the king's testis and fat body, whereas some histone modification genes are remarkably expressed in the king's testis and queen's ovary. These findings imply that epigenetic modification plays important roles in the gamete production process in termite kings and queens.

Comparative Genomics Identifies Putative Signatures of Sociality in Spiders

Comparative genomics has begun to elucidate the genomic basis of social life in insects, but insight into the genomic basis of spider sociality has lagged behind. To begin, to characterize genomic signatures associated with the evolution of social life in spiders, we performed one of the first spider comparative genomics studies including five solitary species and two social species, representing two independent origins of sociality in the genus Stegodyphus. We found that the two social spider species had a large expansion of gene families associated with transport and metabolic processes and an elevated genome-wide rate of molecular evolution compared with the five solitary spider species. Genes that were rapidly evolving in the two social species relative to the five solitary species were enriched for transport, behavior, and immune functions, whereas genes that were rapidly evolving in the solitary species were enriched for energy metabolism processes. Most rapidly evolving genes in the social species Stegodyphus dumicola were broadly expressed across four tissues and enriched for transport functions, but 12 rapidly evolving genes showed brain-specific expression and were enriched for social behavioral processes. Altogether, our study identifies putative genomic signatures and potential candidate genes associated with spider sociality. These results indicate that future spider comparative genomic studies, including broader sampling and additional independent origins of sociality, can further clarify the genomic causes and consequences of social life.

Transcriptome analysis of life stages of the house cricket, Acheta domesticus, to improve insect crop production

To develop genetic resources for the improvement of insects as food, we sequenced transcripts from embryos, one-day hatchlings, three nymphal stages, and male and female adults of the house cricket, Acheta domesticus. A draft transcriptome was assembled from more than 138 million sequences combined from all life stages and sexes. The draft transcriptome assembly contained 45,866 contigs, and more than half were similar to sequences at NCBI (e value < e⁻³). The highest sequence identity was found in sequences from the termites Cryptotermes secundus and Zootermopsis nevadensis. Sequences with identity to Gregarina niphandrodes suggest that these crickets carry the parasite. Among all life stages, there were 5,042 genes with differential expression between life stages (significant at p < 0.05). An enrichment analysis of gene ontology terms from each life stage or sex highlighted genes that were important to biological processes in cricket development. We further characterized genes that may be important in future studies of genetically modified crickets for improved food production, including those involved in RNA interference, and those encoding prolixicin and hexamerins. The data represent an important first step in our efforts to provide genetically improved crickets for human consumption and livestock feed.

The modular nature of protein evolution: domain rearrangement rates across eukaryotic life

BACKGROUND

Modularity is important for evolutionary innovation. The recombination of existing units to form larger complexes with new functionalities spares the need to create novel elements from scratch. In proteins, this principle can be observed at the level of protein domains, functional subunits which are regularly rearranged to acquire new functions.

RESULTS

In this study we analyse the mechanisms leading to new domain arrangements in five major eukaryotic clades (vertebrates, insects, fungi, monocots and eudicots) at unprecedented depth and breadth. This allows, for the first time, to directly compare rates of rearrangements between different clades and identify both lineage specific and general patterns of evolution in the context of domain rearrangements. We analyse arrangement changes along phylogenetic trees by reconstructing ancestral domain content in combination with feasible single step events, such as fusion or fission. Using this approach we explain up to 70% of all rearrangements by tracing them back to their precursors. We find that rates in general and the ratio between these rates for a given clade in particular, are highly consistent across all clades. In agreement with previous studies, fusions are the most frequent event leading to new domain arrangements. A lineage specific pattern in fungi reveals exceptionally high loss rates compared to other clades, supporting recent studies highlighting the importance of loss for evolutionary innovation. Furthermore, our methodology allows us to link domain emergences at specific nodes in the phylogenetic tree to important functional developments, such as the origin of hair in mammals.

CONCLUSIONS

Our results demonstrate that domain rearrangements are based on a canonical set of mutational events with rates which lie within a relatively narrow and consistent range. In addition, gained knowledge about these rates provides a basis for advanced domain-based methodologies for phylogenetics and homology analysis which complement current sequence-based methods.

Gene content evolution in the arthropods

BACKGROUND

Arthropods comprise the largest and most diverse phylum on Earth and play vital roles in nearly every ecosystem. Their diversity stems in part from variations on a conserved body plan, resulting from and recorded in adaptive changes in the genome. Dissection of the genomic record of sequence change enables broad questions regarding genome evolution to be addressed, even across hyper-diverse taxa within arthropods.

RESULTS

Using 76 whole genome sequences representing 21 orders spanning more than 500 million years of arthropod evolution, we document changes in gene and protein domain content and provide temporal and phylogenetic context for interpreting these innovations. We identify many novel gene families that arose early in the evolution of arthropods and during the diversification of insects into modern orders. We reveal unexpected variation in patterns of DNA methylation across arthropods and examples of gene family and protein domain evolution coincident with the appearance of notable phenotypic and physiological adaptations such as flight, metamorphosis, sociality, and chemoperception.

CONCLUSIONS

These analyses demonstrate how large-scale comparative genomics can provide broad new insights into the genotype to phenotype map and generate testable hypotheses about the evolution of animal diversity.

Molecular characterization and distribution of the voltage-gated sodium channel, Para, in the brain of the grasshopper and vinegar fly

Voltage-gated sodium (NaV) channels, encoded by the gene para, play a critical role in the rapid processing and propagation of visual information related to collision avoidance behaviors. We investigated their localization by immunostaining the optic lobes and central brain of the grasshopper Schistocerca americana and the vinegar fly Drosophila melanogaster with an antibody that recognizes the channel peptide domain responsible for fast inactivation gating. NaV channels were detected at high density at all stages of development. In the optic lobe, they revealed stereotypically repeating fascicles consistent with the regular structure of the eye. In the central brain, major axonal tracts were strongly labeled, particularly in the grasshopper olfactory system. We used the NaV channel sequence of Drosophila to identify an ortholog in the transcriptome of Schistocerca. The grasshopper, vinegar fly, and human NaV channels exhibit a high degree of conservation at gating and ion selectivity domains. Comparison with three species evolutionarily close to Schistocerca identified splice variants of Para and their relation to those of Drosophila. The anatomical distribution of NaV channels molecularly analogous to those of humans in grasshoppers and vinegar flies provides a substrate for rapid signal propagation and visual processing in the context of visually-guided collision avoidance.

Involvement of integument-rich CYP4G19 in hydrocarbon biosynthesis and cuticular penetration resistance in Blattella germanica (L.)

BACKGROUND

Cuticle penetration plays an important role as a mechanism of insecticide resistance, but the underlying molecular mechanism remains poorly understood. In Blattella germanica, the cytochrome P450 gene, CYP4G19, is overexpressed in a pyrethroid-resistant strain. Here, we investigated whether CYP4G19 is involved in the biosynthesis of hydrocarbons and further contributes to cuticular penetration resistance in B. germanica.

RESULTS

Compared with the susceptible strain, pyrethroid-resistant cockroaches showed lower cuticular permeability with Eosin Y staining. Removal of epicuticular lipids, mainly nonpolar hydrocarbons, with a hexane wash intensified the cuticular permeability and decreased the resistance index of the resistant strain. CYP4G19 was predominately expressed in the abdominal integument and could be upregulated by desiccation stress or short exposure to beta-cypermethrin. Overexpression of CYP4G19 in the resistant strain was positively correlated with a higher level of cuticular hydrocarbons (CHCs). RNAi-mediated knockdown of CYP4G19 significantly decreased its expression and caused a reduction in CHCs. Meanwhile, CYP4G19 suppression resulted in a non-uniform array of the lipid layer, enhanced cuticle permeability, and compromised insecticide tolerance.

CONCLUSION

Our findings confirm that CYP4G19 is involved in hydrocarbon production and appears to contribute to hydrocarbon-based penetration resistance in B. germanica. This study highlights the lipid-based penetration resistance, advancing our understanding of the molecular mechanism underlying P450-mediated cuticular penetration resistance in insects. © 2019 Society of Chemical Industry.

The de novo transcriptome of workers head of the higher group termite Globitermes sulphureus Haviland (Blattodea: Termitidae)

The subterranean termite Globitermus sulphureus is an important Southeast Asian pest with limited genomic resources that causes damages to agriculture crops and building structures. Therefore, the main goal of this study was to survey the G. sulphureus transcriptome composition. Here, we performed de novo transcriptome for G. sulphureus workers’ heads using Illumina HiSeq paired-end sequencing technology. A total of 88, 639, 408 clean reads were collected and assembled into 243, 057 transcripts and 193, 344 putative genes. The transcripts were annotated with the Trinotate pipeline. In total, 27, 061 transcripts were successfully annotated using BLASTX against the SwissProt database and 17, 816 genes were assigned to 47, 598 GO terms. We classified 14, 223 transcripts into COG classification, resulting in 25 groups of functional annotations. Next, a total of 12, 194 genes were matched in the KEGG pathway and 392 metabolic pathways were predicted based on the annotation. Moreover, we detected two endogenous cellulases in the sequences. The RT-qPCR analysis showed that there were significant differences in the expression levels of two genes β-glucosidase and endo-β-1,4-glucanase between worker and soldier heads of G. sulphureus. This is the first study to characterize the complete head transcriptome of a higher termite G. sulphureus using a high-throughput sequencing. Our study may provide an overview and comprehensive molecular resource for comparative studies of the transcriptomics and genomics of termites.

The X chromosome of the German cockroach, Blattella germanica, is homologous to a fly X chromosome despite 400 million years divergence

Background

Sex chromosome evolution is a dynamic process that can proceed at varying rates across lineages. For example, different chromosomes can be sex-linked between closely related species, whereas other sex chromosomes have been conserved for > 100 million years. Cases of long-term sex chromosome conservation could be informative of factors that constrain sex chromosome evolution. Cytological similarities between the X chromosomes of the German cockroach (Blattella germanica) and most flies suggest that they may be homologous—possibly representing an extreme case of long-term conservation.

Results

To test the hypothesis that the cockroach and fly X chromosomes are homologous, we analyzed whole-genome sequence data from cockroaches. We found evidence in both sequencing coverage and heterozygosity that a significant excess of the same genes are on both the cockroach and fly X chromosomes. We also present evidence that the candidate X-linked cockroach genes may be dosage compensated in hemizygous males. Consistent with this hypothesis, three regulators of transcription and chromatin on the fly X chromosome are conserved in the cockroach genome.

Conclusions

Our results support our hypothesis that the German cockroach shares the same X chromosome as most flies. This may represent the convergent evolution of the X chromosome in the lineages leading to cockroaches and flies. Alternatively, the common ancestor of most insects may have had an X chromosome that resembled the extant cockroach and fly X. Cockroaches and flies diverged ∼ 400 million years ago, which would be the longest documented conservation of a sex chromosome. Cockroaches and flies have different mechanisms of sex determination, raising the possibility that the X chromosome was conserved despite the evolution of the sex determination pathway.

Genome-wide identification of chemosensory receptor genes in the small brown planthopper, Laodelphax striatellus

The small brown planthopper (SBPH), Laodelphax striatellus is one of the major insect pests of rice, but little is known about the molecular-level means by which it locates its hosts. SBPH host-seeking behavior heavily relies on chemosensory receptors (CRs). In this study, we utilized genome analysis of the SBPH to identify 169 CRs, including: 133 odorant receptors (ORs), 13 gustatory receptors (GRs) and 23 ionotropic receptors (IRs). The phylogenetic relationships of OR genes from three rice planthoppers and other insect species revealed that the odorant co-receptor (Orco) clade is the most conserved group. Among the candidate GRs, two sugar receptors and five fructose receptors have been identified but no carbon dioxide receptors investigated. Furthermore, we identified homologs of the three highly conserved IR co-receptors. The obtained results will provide us with precious information needed to better understand the interaction between insect pests and crop plants required for effective crop protection.

Keeping Allergen Names Clear and Defined

The World Health Organization/International Union of Immunological Societies (WHO/IUIS) Allergen Nomenclature Sub-Committee was established in 1986 by leading allergists to standardize names given to proteins that cause IgE-mediated reactions in humans. The Sub-Committee's objective is to assign unique names to allergens based on a critical analysis of confidentially submitted biochemical and clinical data from researchers, often prior to publication to preserve consistency. The Sub-Committee maintains and revises the database as the understanding of allergens evolves. This report summarizes recent developments that led to updates in classification of cockroach group 1 and 5 allergens to animal as well as environmental and occupational allergens. Interestingly, routes, doses, and frequency of exposure often affects allergenicity as does the biochemical properties of the proteins and similarity to self and other proteins. Information required by the Sub-Committee now is more extensive than previously as technology has improved. Identification of new allergens requires identification of the amino acid sequence and physical characteristics of the protein as well as demonstration of IgE binding from subjects verified by described clinical histories, proof of the presence of the protein in relevant exposure substances, and demonstration of biological activity (skin prick tests, activation of basophils, or mast cells). Names are assigned based on taxonomy with the abbreviation of genus and species and assignment of a number, which reflects the priority of discovery, but more often now, the relationships with homologous proteins in related species.

Ecological specificity of the metagenome in a set of lower termite species supports contribution of the microbiome to adaptation of the host

Background

Elucidating the interplay between hosts and their microbiomes in ecological adaptation has become a central theme in evolutionary biology. A textbook example of microbiome-mediated adaptation is the adaptation of lower termites to a wood-based diet, as they depend on their gut microbiome to digest wood. Lower termites have further adapted to different life types. Termites of the wood-dwelling life type never leave their nests and feed on a uniform diet. Termites of the foraging life type forage for food outside the nest and have access to other nutrients. Here we sought to investigate whether the microbiome that is involved in food substrate breakdown and nutrient acquisition might contribute to adaptation to these dietary differences. We reasoned that this should leave ecological imprints on the microbiome.

Results

We investigated the protist and bacterial microbiomes of a total of 29 replicate colonies from five termite species, covering both life types, using metagenomic shotgun sequencing. The microbiome of wood-dwelling species with a uniform wood diet was enriched for genes involved in lignocellulose degradation. Furthermore, metagenomic patterns suggest that the microbiome of wood-dwelling species relied primarily on direct fixation of atmospheric nitrogen, while the microbiome of foraging species entailed the necessary pathways to utilize nitrogen in the form of nitrate for example from soil.

Conclusion

Our findings are consistent with the notion that the microbiome of wood-dwelling species bears an imprint of its specialization on degrading a uniform wood diet, while the microbiome of the foraging species might reflect its adaption to access growth limiting nutrients from more diverse sources. This supports the idea that specific subsets of functions encoded by the microbiome can contribute to host adaptation.

Transcriptomic and Functional Analyses of Phenotypic Plasticity in a Higher Termite, Macrotermes barneyi Light

Eusocial termites have a complex caste system, which leads to the division of labor. Previous studies offered some insight into the caste differentiation in lower termites; however, few studies were focusing on the molecular mechanisms of higher termites with sophisticated societies. Comparative transcriptomic analyses of five immature castes of a higher termite, Macrotermes barneyi Light, suggest that phenotypic plasticity is modulated by an array of transcriptional changes, including differentially expressed genes (e.g., caste-biased genes Vtg and TnC), co-expression networks (e.g., genes associated with nymph reproduction), and alternative splicing (e.g., events related to muscle development in presoldiers). Transcriptional (RT-PCR and RT-qPCR) and functional (in vivo RNAi) validation studies reveal multiple molecular mechanisms contributing to the phenotypic plasticity in eusocial termites. Molecular mechanisms governing the phenotypic plasticity in M. barneyi could be a rule rather than an exception in the evolution of sociality.

Transcriptome changes reveal the genetic mechanisms of the reproductive plasticity of workers in lower termites

Background

The reproductive plasticity of termite workers provides colonies with tremendous flexibility to respond to environmental changes, which is the basis for evolutionary and ecological success. Although it is known that all colony members share the same genetic background and that differences in castes are caused by differences in gene expression, the pattern of the specific expression of genes involved in the differentiation of workers into reproductives remains unclear. In this study, the isolated workers of Reticulitermes labralis developed into reproductives, and then comparative transcriptomes were used for the first time to reveal the molecular mechanisms underlying the reproductive plasticity of workers.

Results

We identified 38,070 differentially expressed genes and found a pattern of gene expression involved in the differentiation of the workers into reproductives. 12, 543 genes were specifically upregulated in the isolated workers. Twenty-five signal transduction pathways classified into environmental information processing were related to the differentiation of workers into reproductives. Ras functions as a signalling switch regulates the reproductive plasticity of workers. The catalase gene which is related to longevity was up-regulated in reproductives.

Conclusion

We demonstrate that workers leaving the natal colony can induce the expression of stage-specific genes in the workers, which leads to the differentiation of workers into reproductives and suggests that the signal transduction along the Ras-MAPK pathway crucially controls the reproductive plasticity of the workers. This study also provides an important model for revealing the molecular mechanism of longevity changes.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6037-y) contains supplementary material, which is available to authorized users.

2018 Highlights of Urban Entomology

The field of urban entomology is primarily associated with the study of ants, bed bugs, cockroaches, termites, and other occasional invader pests that are found within or near human-made structures. A wide array of peer-reviewed studies were published in 2018. The topics of these articles ranged from genomes and basic biology of urban insects to various applied aspects of pest management. Key findings of these papers are presented and discussed from the perspective of the contributions they make to the discipline of urban entomology. Additionally, potential future research opportunities that are evident from these publications have been outlined.

BgFas1: A fatty acid synthase gene required for both hydrocarbon and cuticular fatty acid biosynthesis in the German cockroach, Blattella germanica (L.)

Insect cuticular hydrocarbons (CHCs), the evolutionary products of aquatic hexapod ancestors expanding to terrestrial environment, are deposited on the surface of insect integument and originally functioned primarily as waterproofing agents. CHCs are derived from the conserved fatty acid synthesis pathway in insects. However, the pivotal fatty acid synthase (FAS) involved in hydrocarbon (HC) biosynthesis remains unknown in many insect orders including the primitive Blattodea. Here, we investigated functional FAS genes that modulate cuticular lipid biogenesis in the German cockroach, Blattella germanica (L.). Based on our full-length transcriptomic data and the available genomic data, seven FAS genes (BgFas1-7) were identified from B. germanica. Tissue-specific expression analysis revealed that BgFas1, BgFas3, BgFas4 and BgFas7 were highly expressed in the integument, whereas BgFas2 was dominantly expressed in the fat body. BgFas5/6 mRNA was almost negligible in the tested tissues. Systemic RNAi screen was performed against BgFas1-7, we found that only RNAi knockdown of BgFas1 caused a dramatic reduction of methyl-branched HCs (mbHCs) and a slight decrease of straight-chain HCs (scHCs) for both internal and external HCs. Significant reduction of cuticular free fatty acids (cFFAs) was also detected within BgFas1-repressed cockroaches, while repression of CYP4G19 resulted in dramatic increase of cFFAs. Moreover, we found that BgFas1 mRNA levels were correlated with insect molting cycles, and could be induced by long-term mild dryness treatment. Furthermore, desiccation assay revealed that BgFas1 suppression accelerated water loss and led to early death of cockroaches under desiccation. Our results indicate that BgFas1 is necessary for both HC and cFFA biosynthesis in B. germanica. In addition, our study also confirms that cuticular lipids, particularly mbCHCs, are critical for desiccation resistance in B. germanica.

Leveraging technological innovations to investigate evolutionary transitions to eusociality

The study of the major transition to eusociality presents several challenges to researchers, largely resulting from the importance of complex behavioral phenotypes and the shift from individual to group level selection. These challenges are being met with corresponding technological improvements. Advances in resource development for non-model taxa, behavioral tracking, nucleic acid sequencing, and reverse genetics are facilitating studies of hypotheses that were previously intractable. These innovations are resulting in the development of new model systems tailored to the exploration of specific behavioral phenotypes and the querying of underlying molecular mechanisms that drive eusocial behaviors. Here, we present a brief overview of how methodological innovations are advancing our understanding of the evolution of eusociality.