A genomic comparison of two termites with different social complexity

No single rescue recipe: genome complexities modulate insect response to climate change

Declines in insect populations have gained formidable attention. Given their crucial role in the ecosystem, the causes of declining insect populations must be investigated. However, the insect clade has been associated with low extinction and high diversification rates. It is unlikely that insects underwent mass extinctions in the past. However, the pace of current climate change could make insect populations vulnerable to extinction. We propose genome size (GS) and transposable elements (TEs) to be rough estimates to assess extinction risk. Larger GS and/or proliferating TEs have been associated with adaptation in rapid climate change scenarios. We speculate that unstable, stressful environmental conditions are strongly associated with GS and TE expansion, which could be further correlated with adaptations. Alternately, stressful conditions trigger TE bursts that are not purged in smaller populations. GS and TE loads could be indicators of small effective populations in the wild, likely experiencing bottlenecks or drastic climatic perturbations, which calls for an urgent assessment of extinction risk.

Chemical and transcriptomic diversity do not correlate with ascending levels of social complexity in the insect order Blattodea

Eusocial insects, such as ants and termites, are characterized by high levels of coordinated social organization. This is contrasted by solitary insects that display more limited forms of collective behavior. It has been hypothesized that this gradient in sociobehavioral sophistication is positively correlated with chemical profile complexity, due to a potentially increased demand for diversity in chemical communication mechanisms in insects with higher levels of social complexity. However, this claim has rarely been assessed empirically. Here, we compare different levels of chemical and transcriptomic complexity in selected species of the order Blattodea that represent different levels of social organization, from solitary to eusocial. We primarily focus on cuticular hydrocarbon (CHC) complexity, since it has repeatedly been demonstrated that CHCs are key signaling molecules conveying a wide variety of chemical information in solitary as well as eusocial insects. We assessed CHC complexity and divergence between our studied taxa of different social complexity levels as well as the differentiation of their respective repertoires of CHC biosynthesis gene transcripts. Surprisingly, we did not find any consistent pattern of chemical complexity correlating with social complexity, nor did the overall chemical divergence or transcriptomic repertoire of CHC biosynthesis genes reflect on the levels of social organization. Our results challenge the assumption that increasing social complexity is generally reflected in more complex chemical profiles and point toward the need for a more cautious and differentiated view on correlating complexity on a chemical, genetic, and social level.

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.

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.

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.

More effective transposon regulation in fertile, long-lived termite queens than in sterile workers

Transposable elements (TEs) are mobile genetic sequences, which can cause the accumulation of genomic damage in the lifetime of an organism. The regulation of TEs, for instance via the piRNA-pathway, is an important mechanism to protect the integrity of genomes, especially in the germ-line where mutations can be transmitted to offspring. In eusocial insects, soma and germ-line are divided among worker and reproductive castes, so one may expect caste-specific differences in TE regulation to exist. To test this, we compared whole-genome levels of repeat element transcription in the fat body of female workers, kings and five different queen stages of the higher termite, Macrotermes natalensis. In this species, queens can live over 20 years, maintaining near maximum reproductive output, while sterile workers only live weeks. We found a strong, positive correlation between TE expression and the expression of neighbouring genes in all castes. However, we found substantially higher TE activity in workers than in reproductives. Furthermore, TE expression did not increase with age in queens, despite a sevenfold increase in overall gene expression, due to a significant upregulation of the piRNA-pathway in 20-year-old queens. Our results suggest a caste- and age-specific regulation of the piRNA-pathway has evolved in higher termites that is analogous to germ-line-specific activity in solitary organisms. In the fat body of these termite queens, an important metabolic tissue for maintaining their extreme longevity and reproductive output, an efficient regulation of TEs likely protects genome integrity, thus further promoting reproductive fitness even at high age.

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.

Complex regulatory role of DNA methylation in caste- and age-specific expression of a termite

The reproductive castes of eusocial insects are often characterized by extreme lifespans and reproductive output, indicating an absence of the fecundity/longevity trade-off. The role of DNA methylation in the regulation of caste- and age-specific gene expression in eusocial insects is controversial. While some studies find a clear link to caste formation in honeybees and ants, others find no correlation when replication is increased across independent colonies. Although recent studies have identified transcription patterns involved in the maintenance of high reproduction throughout the long lives of queens, the role of DNA methylation in the regulation of these genes is unknown. We carried out a comparative analysis of DNA methylation in the regulation of caste-specific transcription and its importance for the regulation of fertility and longevity in queens of the higher termite Macrotermes natalensis. We found evidence for significant, well-regulated changes in DNA methylation in mature compared to young queens, especially in several genes related to ageing and fecundity in mature queens. We also found a strong link between methylation and caste-specific alternative splicing. This study reveals a complex regulatory role of fat body DNA methylation both in the division of labour in termites, and during the reproductive maturation of queens.

Termite sociogenomics: evolution and regulation of caste-specific expressed genes

Termite genomes have been sequenced in at least five species from four different families. Genome-based transcriptome analyses have identified large numbers of protein-coding genes with caste-specific expression patterns. These genes include those involved in caste-specific morphologies and roles, for example high fecundity and longevity in reproductives. Some caste-specific expressed genes belong to multi-gene families, and their genetic architecture and expression profiles indicate they have evolved via tandem gene duplication. Candidate regulatory mechanisms of caste-specific expression include epigenetic regulation (e.g. histone modification and non-coding RNA) and diversification of transcription factors and cis-regulatory elements. We review current knowledge in the area of termite sociogenomics, focussing on the evolution and regulation of caste-specific expressed genes, and discuss future research directions.

Satellitome Analysis and Transposable Elements Comparison in Geographically Distant Populations of Spodoptera frugiperda

Spodoptera frugiperda (fall armyworm) is a member of the superfamily Noctuoidea that accounts for more than a third of all Lepidoptera and includes a considerable number of agricultural and forest pest species. Spodoptera frugiperda is a polyphagous species that is a significant agricultural pest worldwide, emphasizing its economic importance. Spodoptera frugiperda’s genome size, assembly, phylogenetic classification, and transcriptome analysis have all been previously described. However, the different studies reported different compositions of repeated DNA sequences that occupied the whole assembled genome, and the Spodoptera frugiperda genome also lacks the comprehensive study of dynamic satellite DNA. We conducted a comparative analysis of repetitive DNA across geographically distant populations of Spodoptera frugiperda, particularly satellite DNA, using publicly accessible raw genome data from eight different geographical regions. Our results showed that most transposable elements (TEs) were commonly shared across all geographically distant samples, except for the Maverick and PIF/Harbinger elements, which have divergent repeat copies. The TEs age analysis revealed that most TEs families consist of young copies 1–15 million years old; however, PIF/Harbinger has some older/degenerated copies of 30–35 million years old. A total of seven satellite DNA families were discovered, accounting for approximately 0.65% of the entire genome of the Spodoptera frugiperda fall armyworm. The repeat profiling analysis of satellite DNA families revealed differential read depth coverage or copy numbers. The satellite DNA families range in size from the lowest 108 bp SfrSat06-108 families to the largest (1824 bp) SfrSat07-1824 family. We did not observe a statistically significant correlation between monomer length and K2P divergence, copy number, or abundance of each satellite family. Our findings suggest that the satellite DNA families identified in Spodoptera frugiperda account for a considerable proportion of the genome’s repetitive fraction. The satellite DNA families’ repeat profiling revealed a point mutation along the reference sequences. Limited TEs differentiation exists among geographically distant populations of Spodoptera frugiperda.

Genomic and transcriptomic analyses of the subterranean termite Reticulitermes speratus: Gene duplication facilitates social evolution

Gene duplication is a major source of evolutionary innovation and is associated with the increases in biological complexity and adaptive radiation. Termites are model social organisms characterized by a sophisticated caste system. We analyzed the genome of the Japanese subterranean termite, an ecologically and economically important insect acting as a destructive pest. The analyses revealed the significance of gene duplication in social evolution. Gene duplication associated with caste-biased gene expression was prevalent in the termite genome. Many of the duplicated genes were related to social functions, such as chemical communication, social immunity, and defense, and they were often expressed in caste-specific organs. We propose that gene duplication facilitates social evolution through regulatory diversification leading to caste-biased expression and functional specialization.

Termites are model social organisms characterized by a polyphenic caste system. Subterranean termites (Rhinotermitidae) are ecologically and economically important species, including acting as destructive pests. Rhinotermitidae occupies an important evolutionary position within the clade representing a transitional taxon between the higher (Termitidae) and lower (other families) termites. Here, we report the genome, transcriptome, and methylome of the Japanese subterranean termite Reticulitermes speratus. Our analyses highlight the significance of gene duplication in social evolution in this termite. Gene duplication associated with caste-biased gene expression was prevalent in the R. speratus genome. The duplicated genes comprised diverse categories related to social functions, including lipocalins (chemical communication), cellulases (wood digestion and social interaction), lysozymes (social immunity), geranylgeranyl diphosphate synthase (social defense), and a novel class of termite lineage–specific genes with unknown functions. Paralogous genes were often observed in tandem in the genome, but their expression patterns were highly variable, exhibiting caste biases. Some of the assayed duplicated genes were expressed in caste-specific organs, such as the accessory glands of the queen ovary and the frontal glands of soldier heads. We propose that gene duplication facilitates social evolution through regulatory diversification, leading to caste-biased expression and subfunctionalization and/or neofunctionalization conferring caste-specialized functions.

Lifespan prolonging mechanisms and insulin upregulation without fat accumulation in long-lived reproductives of a higher termite

Kings and queens of eusocial termites can live for decades, while queens sustain a nearly maximal fertility. To investigate the molecular mechanisms underlying their long lifespan, we carried out transcriptomics, lipidomics and metabolomics in Macrotermes natalensis on sterile short-lived workers, long-lived kings and five stages spanning twenty years of adult queen maturation. Reproductives share gene expression differences from workers in agreement with a reduction of several aging-related processes, involving upregulation of DNA damage repair and mitochondrial functions. Anti-oxidant gene expression is downregulated, while peroxidability of membranes in queens decreases. Against expectations, we observed an upregulated gene expression in fat bodies of reproductives of several components of the IIS pathway, including an insulin-like peptide, Ilp9. This pattern does not lead to deleterious fat storage in physogastric queens, while simple sugars dominate in their hemolymph and large amounts of resources are allocated towards oogenesis. Our findings support the notion that all processes causing aging need to be addressed simultaneously in order to prevent it.

Eusociality in snapping shrimps is associated with larger genomes and an accumulation of transposable elements

Despite progress uncovering the genomic underpinnings of sociality, much less is known about how social living affects the genome. In different insect lineages, for example, eusocial species show both positive and negative associations between genome size and structure, highlighting the dynamic nature of the genome. Here, we explore the relationship between sociality and genome architecture in Synalpheus snapping shrimps that exhibit multiple origins of eusociality and extreme interspecific variation in genome size. Our goal is to determine whether eusociality leads to an accumulation of repetitive elements and an increase in genome size, presumably due to reduced effective population sizes resulting from a reproductive division of labor, or whether an initial accumulation of repetitive elements leads to larger genomes and independently promotes the evolution of eusociality through adaptive evolution. Using phylogenetically informed analyses, we find that eusocial species have larger genomes with more transposable elements (TEs) and microsatellite repeats than noneusocial species. Interestingly, different TE subclasses contribute to the accumulation in different species. Phylogenetic path analysis testing alternative causal relationships between sociality and genome architecture is most consistent with the hypothesis that TEs modulate the relationship between sociality and genome architecture. Although eusociality appears to influence TE accumulation, ancestral state reconstruction suggests moderate TE abundances in ancestral species could have fueled the initial transitions to eusociality. Ultimately, we highlight a complex and dynamic relationship between genome and social evolution, demonstrating that sociality can influence the evolution of the genome, likely through changes in demography related to patterns of reproductive skew.

Why and how do termite kings and queens live so long?

Lifespan varies greatly across the tree of life. Of the various explanations for this phenomenon, those that involve trade-offs between reproduction and longevity have gained considerable support. There is an important exception: social insect reproductives (queens and in termites, also kings) exhibit both high reproductive outputs and extraordinarily long lives. As both the ultimate and proximate mechanisms underlying the absence of the fecundity/longevity trade-off could shed light on the unexpected dynamics and molecular mechanisms of extended longevity, reproductives of social insects have attracted much attention in the field of ageing research. Here, we highlight current ecological and physiological studies on ageing and discuss the various possible evolutionary and molecular explanations of the extended lifespans of termite reproductives. We integrate these findings into a coherent framework revealing the evolution of longevity in these reproductives. Studies on termites may explain why and how ageing is shaped by natural selection. This article is part of the theme issue 'Ageing and sociality: why, when and how does sociality change ageing patterns?'

Comparative transcriptomic analysis of the mechanisms underpinning ageing and fecundity in social insects

The exceptional longevity of social insect queens despite their lifelong high fecundity remains poorly understood in ageing biology. To gain insights into the mechanisms that might underlie ageing in social insects, we compared gene expression patterns between young and old castes (both queens and workers) across different lineages of social insects (two termite, two bee and two ant species). After global analyses, we paid particular attention to genes of the insulin/insulin-like growth factor 1 signalling (IIS)/target of rapamycin (TOR)/juvenile hormone (JH) network, which is well known to regulate lifespan and the trade-off between reproduction and somatic maintenance in solitary insects. Our results reveal a major role of the downstream components and target genes of this network (e.g. JH signalling, vitellogenins, major royal jelly proteins and immune genes) in affecting ageing and the caste-specific physiology of social insects, but an apparently lesser role of the upstream IIS/TOR signalling components. Together with a growing appreciation of the importance of such downstream targets, this leads us to propose the TI-J-LiFe (TOR/IIS-JH-Lifespan and Fecundity) network as a conceptual framework for understanding the mechanisms of ageing and fecundity in social insects and beyond. This article is part of the theme issue 'Ageing and sociality: why, when and how does sociality change ageing patterns?'

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.

Natural variation in colony inbreeding does not influence susceptibility to a fungal pathogen in a termite

Reduced genetic diversity through inbreeding can negatively affect pathogen resistance. This relationship becomes more complicated in social species, such as social insects, since the chance of disease transmission increases with the frequency of interactions among individuals. However, social insects may benefit from social immunity, whereby individual physiological defenses may be bolstered by collective-level immune responses, such as grooming or sharing of antimicrobial substance through trophallaxis. We set out to determine whether differences in genetic diversity between colonies of the subterranean termite, Reticulitermes flavipes, accounts for colony survival against pathogens. We sampled colonies throughout the United States (Texas, North Carolina, Maryland, and Massachusetts) and determined the level of inbreeding of each colony. To assess whether genetically diverse colonies were better able to survive exposure to diverse pathogens, we challenged groups of termite workers with two strains of a pathogenic fungus, one local strain present in the soil surrounding sampled colonies and another naïve strain, collected outside the range of this species. We found natural variation in the level of inbreeding between colonies, but this variation did not explain differences in susceptibility to either pathogen. Although the naïve strain was found to be more hazardous than the local strain, colony resistance was correlated between two strains, meaning that colonies had either relatively high or low susceptibility to both strains regardless of their inbreeding coefficient. Overall, our findings may reflect differential virulence between the strains, immune priming of the colonies via prior exposure to the local strain, or a coevolved resistance toward this strain. They also suggest that colony survival may rely more upon additional factors, such as different behavioral response thresholds or the influence of a specific genetic background, rather than the overall genetic diversity of the colony.

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.

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.

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.

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.

Evolutionary relationship of fat body endoreduplication and queen fecundity in termites

Endoreduplication or nuclear genome replication without cell division is widely observed in the metabolically active tissues of plants and animals. The fat body cells of adult female insects produce abundant yolk proteins and become polyploid, which is assumed to accelerate egg production. Recently, it was reported that in termites, endopolyploidy in the fat body occurs only in queens but not in the other females; however, the relationship between the fecundity and ploidy level in the fat body remains unclear. Termite queens exhibit a huge variation in their egg-producing capacity among different species; queens in the species with a foraging lifestyle, in which workers leave the nest to forage outside, are much more fecund than those in the species living in a single piece of wood. In this study, we conducted ploidy analyses on three foragings and three wood-dwelling termites via flow cytometry. In all the species, the fat body of queens contained significantly more polyploid cells than that of other nonreproductive females, considering their body size effect. However, the male fat body, which is not involved in yolk production, did not show consistency in polyploid cell numbers among the species studied. Moreover, highly fecund queens in foraging termites exhibit higher levels of endopolyploidy in their fat body than those with less fecundity in wood-dwelling termites. These results suggest that endopolyploidy in the fat body of termite queens can boost their egg production, and the level of endopolyploidy in their fat body is linked to their fecundity. Our study provides a novel insight into the evolutionary relationship between endoreduplication and caste specialization in social insects.

Comparative genomics and transcriptomics of host-pathogen interactions in insects: evolutionary insights and future directions

Classical evolutionary studies of protein-coding genes have established that genes in the canonical immune system are often among the most rapidly evolving within and between species. As more genomes and transcriptomes across insects are sequenced, it is becoming clear that duplications and losses of immune genes are also a likely consequence of host-pathogen interactions. Furthermore, particular species respond to diverse pathogenic challenges with a wide range of challenge-specific responses that are still poorly understood. Transcriptional studies, using RNA-seq to characterize the infection-regulated transcriptome of diverse insects, are crucial for additional progress in understanding the ecology and evolution of the full complexity of the host response.

The evolution of sociality in termites from cockroaches: A taxonomic and phylogenetic perspective

Despite multiple studies and advances, sociality still puzzles evolutionary biologists in numerous ways, which might be partly addressed with the advent of sociogenomics. In insects, the majority of sociogenomic studies deal with Hymenoptera, one of the two groups that evolved eusociality with termites. But, to fully grasp the evolution of sociality, studies must obviously not restrict to eusocial lineages. Multiple kinds of social system transitions have been recorded and they all bring complementary insights. For instance, cockroaches, the closest relatives to termites, display a wide range of social interactions and evolved convergently subsocial behaviors (i.e., brood care). In this context, we emphasize the need for natural history, taxonomic, and phylogenetic studies. Natural history studies provide the foundations on which building hypotheses, whereas taxonomy provides the taxa to sample to test these hypotheses, and phylogenetics brings the historical framework necessary to test evolutionary scenarios of sociality evolution.

Chemical Fertility Signaling in Termites: Idiosyncrasies and Commonalities in Comparison with Ants

Termites evolved eusociality independently from social Hymenoptera. As a common trait, reproductive monopoly is maintained through chemical communication. The queen (and in termites also a king) prevents workers from reproduction by conveying their reproductive status. In termites all soldiers are sterile, but workers' potential to reproduce differs between species. It ranges from totipotency in wood-dwelling lower termites where workers are a transient stage from which all other castes develop, to sterile workers in some higher termites. Intermediate are species in which workers can develop into replacement sexuals within the nest but not into winged sexuals. I summarize the patchy picture about fertility signaling that we currently have for termites, pointing also to potential conflicts over reproduction that differ from those in social Hymenoptera. Recent findings imply that, similar to many social Hymenoptera, wood-dwelling termites that live in confined nests use long-chain cuticular hydrocarbons (CHCs) as fertility signals. Yet other compounds are important as well, comprising proteinaceous secretions and especially volatiles. For a subterranean termite, two volatiles have been identified as primer pheromones that prevent reproductive differentiation of workers. It requires more data to test whether wood-dwelling termites use CHCs, while species with larger colonies and less confined nests use volatiles, or whether all species rely on multicomponent signals. Ultimately, we need more effort to model and test potential conflicts over reproduction between queens, kings and workers. Here results from social Hymenoptera cannot be transferred to termites as the latter are diploid and commonly inbred. This review illustrates promising future research avenues.

The Complete Mitochondrial Genome of Coptotermes 'suzhouensis' (syn. Coptotermes formosanus) (Isoptera: Rhinotermitidae) and Molecular Phylogeny Analysis

Coptotermes suzhouensis (Isoptera: Rhinotermitidae) is a significant subterranean termite pest of wooden structures and is widely distributed in southeastern China. The complete mitochondrial DNA sequence of C. suzhouensis was analyzed in this study. The mitogenome was a circular molecule of 15,764 bp in length, which contained 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and an A+T-rich region with a gene arrangement typical of Isoptera mitogenomes. All PCGs were initiated by ATN codons and terminated by complete termination codons (TAA), except COX2, ND5, and Cytb, which ended with an incomplete termination codon T. All tRNAs displayed a typical clover-leaf structure, except for tRNASer(AGN), which did not contain the stem-loop structure in the DHU arm. The A+T content (69.23%) of the A+T-rich region (949 bp) was higher than that of the entire mitogenome (65.60%), and two different sets of repeat units (A+B) were distributed in this region. Comparison of complete mitogenome sequences with those of Coptotermes formosanus indicated that the two taxa have very high genetic similarity. Forty-one representative termite species were used to construct phylogenetic trees by maximum likelihood, maximum parsimony, and Bayesian inference methods. The phylogenetic analyses also strongly supported (BPP, MLBP, and MPBP = 100%) that all C. suzhouensis and C. formosanus samples gathered into one clade with genetic distances between 0.000 and 0.002. This study provides molecular evidence for a more robust phylogenetic position of C. suzhouensis and inferrs that C. suzhouensis was the synonymy of C. formosanus.

Hemimetabolous genomes reveal molecular basis of termite eusociality

Around 150 million years ago, eusocial termites evolved from within the cockroaches, 50 million years before eusocial Hymenoptera, such as bees and ants, appeared. Here, we report the 2-Gb genome of the German cockroach, Blattella germanica, and the 1.3-Gb genome of the drywood termite Cryptotermes secundus. We show evolutionary signatures of termite eusociality by comparing the genomes and transcriptomes of three termites and the cockroach against the background of 16 other eusocial and non-eusocial insects. Dramatic adaptive changes in genes underlying the production and perception of pheromones confirm the importance of chemical communication in the termites. These are accompanied by major changes in gene regulation and the molecular evolution of caste determination. Many of these results parallel molecular mechanisms of eusocial evolution in Hymenoptera. However, the specific solutions are remarkably different, thus revealing a striking case of convergence in one of the major evolutionary transitions in biological complexity.

Epigenetics and developmental plasticity in orthopteroid insects

Developmental plasticity is a key driver of the extraordinary ecological success of insects. Epigenetic mechanisms provide an important link between the external stimuli that initiate polyphenisms, and the stable changes in gene expression that govern alternative insect morphs. We review the epigenetics of orthopteroid insects, focussing on recent research on locusts and termites, two groups which display high levels of phenotypic plasticity, and for which genome sequences have become available in recent years. We examine research on the potential role of DNA methylation, histone modifications, and non-coding RNAs in the regulation of gene expression in these insects. DNA methylation patterns in orthopteroids share a number of characteristics with those of hymenopteran insects, although methylation levels are much higher, and extend to introns and repeat elements. Future examinations of epigenetic mechanisms in these insects will benefit from comparison of tissues from aged-matched individuals from alternative morphs, and adequate biological replication.

An alternative pathway to eusociality: Exploring the molecular and functional basis of fortress defense

Some animals express a form of eusociality known as "fortress defense," in which defense rather than brood care is the primary social act. Aphids are small plant-feeding insects, but like termites, some species express division of labor and castes of aggressive juvenile "soldiers." What is the functional basis of fortress defense eusociality in aphids? Previous work showed that the acquisition of venoms might be a key innovation in aphid social evolution. We show that the lethality of aphid soldiers derives in part from the induction of exaggerated immune responses in insects they attack. Comparisons between closely related social and nonsocial species identified a number of secreted effector molecules that are candidates for immune modulation, including a convergently recruited protease described in unrelated aphid species with venom-like functions. These results suggest that aphids are capable of antagonizing conserved features of the insect immune response, and provide new insights into the mechanisms underlying the evolution of fortress defense eusociality in aphids.

Why do social insect queens live so long? Approaches to unravel the sociality-aging puzzle

Social insects are characterized by an apparent reshaping of the fecundity/longevity trade-off with sociality. Currently, we have only sketchy information about the potential underlying causes and mechanisms of aging and senescence which in addition are restricted to few model insect organisms (mainly the fruit fly Drosophila melanogaster and the honey bee Apis mellifera). How can we gain a more thorough understanding how sociality shapes senescence and the fecundity/longevity trade-off? By reviewing available literature, I propose a comparative approach that offers the opportunity to gain fundamental insights into uncovering the basis for this life history trade-off and its reshaping with sociality.

Chapter 6. Comparative genomic approaches to investigate molecular traits specific to social insects

Ageing is a feature of nearly all known organisms and, by its connection to survival, appears to trade off with fecundity. However, in some organisms such as in queens of social insects, this negative relation appears reversed and individuals live long and reproduce much. Since new experimental techniques, transcriptomes and genomes of many social insects have recently become available, a comparison of these data in a phylogenetic framework becomes feasible. This allows the study of general trends, species specific oddities and evolutionary dynamics of the molecular properties and changes which underlie ageing, fecundity and the reversal of this negative association. In the framework of social insect evolution, we review the most important recent insights, computational methods, their applications and data resources which are available.

The benefits of expanding studies of trait exaggeration to hemimetabolous insects and beyond morphology

Trait exaggeration, well known to naturalists and evolutionary biologists, has recently become a prominent research subject in the modern field of Evolutionary Developmental Biology. A large number of traits that can be considered as cases of exaggeration exist in nature. Yet, the field has almost exclusively focused on the study of growth-related exaggerated traits in a selection of holometabolous insects. The absence of the hemimetabola from studies of exaggeration leaves a significant gap in our understanding of the development and evolution of such traits. Here we argue that efforts to understand the mechanisms of trait exaggeration would benefit from expanding the study subjects to include other kinds of exaggeration and other model species.

Genomic sources of phenotypic novelty in the evolution of eusociality in insects

Genomic resources are now available for closely related species that vary in social behavior, providing insight on the genomics of social evolution. Changes in the architecture of gene regulatory networks likely influence the evolutionary trajectory of social traits. Evolutionarily novel genes are likely important in the evolution of social diversity among insects, but it is unclear whether new genes played a driving role in the advent or elaboration of eusociality or if they were instead a result of other genomic features of eusociality. The worker phenotype appears to be the center of genetic novelty, but the mechanisms for this remain unresolved. Future studies are needed to understand how genetic novelty arises, becomes incorporated into existing gene regulatory networks, and the effects this has on the evolution of social traits in closely related social and solitary species.

Comparative Analysis of Transcriptomes from Secondary Reproductives of Three Reticulitermes Termite Species

Termites are eusocial insects related to cockroaches that feed on lignocellulose. These insects are key species in ecosystems since they recycle a large amount of nutrients but also are pests, exerting major economic impacts. Knowledge on the molecular pathways underlying reproduction, caste differentiation or lignocellulose digestion would largely benefit from additional transcriptomic data. This study focused on transcriptomes of secondary reproductive females (nymphoid neotenics). Thirteen transcriptomes were used: 10 of Reticulitermes flavipes and R. grassei sequenced from a previous study, and two transcriptomes of R. lucifugus sequenced for the present study. After transcriptome assembly and read mapping, we examined interspecific variations of genes expressed by termites or gut microorganisms. A total of 18,323 orthologous gene clusters were detected. Functional annotation and taxonomic assignment were performed on a total of 41,287 predicted contigs in the three termite species. Between the termite species studied, functional categories of genes were comparable. Gene ontology (GO) terms analysis allowed the discovery of 9 cellulases and a total of 79 contigs potentially involved in 11 enzymatic activities used in wood metabolism. Altogether, results of this study illustrate the strong potential for the use of comparative interspecific transcriptomes, representing a complete resource for future studies including differentially expressed genes between castes or SNP analysis for population genetics.