Comparative analysis of adipokinetic hormones and their receptors in Blattodea reveals novel patterns of gene evolution.

Adipokinetic hormone (AKH) is a neuropeptide produced in the insect corpora cardiaca that plays an essential role in mobilising carbohydrates and lipids from the fat body to the haemolymph. AKH acts by binding to a rhodopsin-like G protein-coupled receptor (GPCR), the adipokinetic hormone receptor (AKHR). In this study, we tackle AKH ligand and receptor gene evolution as well as the evolutionary origins of AKH gene paralogues from the order Blattodea (termites and cockroaches). Phylogenetic analyses of AKH precursor sequences point to an ancient AKH gene duplication event in the common ancestor of Blaberoidea, yielding a new group of putative decapeptides. In total, 16 different AKH peptides from 90 species were obtained. Two octapeptides and seven putatively novel decapeptides are predicted for the first time. AKH receptor sequences from 18 species, spanning solitary cockroaches and subsocial wood roaches as well as lower and higher termites, were subsequently acquired using classical molecular methods and in silico approaches employing transcriptomic data. Aligned AKHR open reading frames revealed 7 highly conserved transmembrane regions, a typical arrangement for GPCRs. Phylogenetic analyses based on AKHR sequences support accepted relationships among termite, subsocial (Cryptocercus spp.) and solitary cockroach lineages to a large extent, while putative post-translational modification sites do not greatly differ between solitary and subsocial roaches and social termites. Our study provides important information not only for AKH and AKHR functional research but also for further analyses interested in their development as potential candidates for biorational pest control agents against invasive termites and cockroaches.

Evidence for reduced immune gene diversity and activity during the evolution of termites.

The evolution of biological complexity is associated with the emergence of bespoke immune systems that maintain and protect organism integrity. Unlike the well-studied immune systems of cells and individuals, little is known about the origins of immunity during the transition to eusociality, a major evolutionary transition comparable to the evolution of multicellular organisms from single-celled ancestors. We aimed to tackle this by characterizing the immune gene repertoire of 18 cockroach and termite species, spanning the spectrum of solitary, subsocial and eusocial lifestyles. We find that key transitions in termite sociality are correlated with immune gene family contractions. In cross-species comparisons of immune gene expression, we find evidence for a caste-specific social defence system in termites, which appears to operate at the expense of individual immune protection. Our study indicates that a major transition in organismal complexity may have entailed a fundamental reshaping of the immune system optimized for group over individual defence.

Unity in defence: honeybee workers exhibit conserved molecular responses to diverse pathogens.

BACKGROUND: Organisms typically face infection by diverse pathogens, and hosts are thought to have developed specific responses to each type of pathogen they encounter. The advent of transcriptomics now makes it possible to test this hypothesis and compare host gene expression responses to multiple pathogens at a genome-wide scale. Here, we performed a meta-analysis of multiple published and new transcriptomes using a newly developed bioinformatics approach that filters genes based on their expression profile across datasets. Thereby, we identified common and unique molecular responses of a model host species, the honey bee (Apis mellifera), to its major pathogens and parasites: the Microsporidia Nosema apis and Nosema ceranae, RNA viruses, and the ectoparasitic mite Varroa destructor, which transmits viruses. RESULTS: We identified a common suite of genes and conserved molecular pathways that respond to all investigated pathogens, a result that suggests a commonality in response mechanisms to diverse pathogens. We found that genes differentially expressed after infection exhibit a higher evolutionary rate than non-differentially expressed genes. Using our new bioinformatics approach, we unveiled additional pathogen-specific responses of honey bees; we found that apoptosis appeared to be an important response following microsporidian infection, while genes from the immune signalling pathways, Toll and Imd, were differentially expressed after Varroa/virus infection. Finally, we applied our bioinformatics approach and generated a gene co-expression network to identify highly connected (hub) genes that may represent important mediators and regulators of anti-pathogen responses. CONCLUSIONS: Our meta-analysis generated a comprehensive overview of the host metabolic and other biological processes that mediate interactions between insects and their pathogens. We identified key host genes and pathways that respond to phylogenetically diverse pathogens, representing an important source for future functional studies as well as offering new routes to identify or generate pathogen resilient honey bee stocks. The statistical and bioinformatics approaches that were developed for this study are broadly applicable to synthesize information across transcriptomic datasets. These approaches will likely have utility in addressing a variety of biological questions.

Elevated virulence of an emerging viral genotype as a driver of honeybee loss.

Emerging infectious diseases (EIDs) have contributed significantly to the current biodiversity crisis, leading to widespread epidemics and population loss. Owing to genetic variation in pathogen virulence, a complete understanding of species decline requires the accurate identification and characterization of EIDs. We explore this issue in the Western honeybee, where increasing mortality of populations in the Northern Hemisphere has caused major concern. Specifically, we investigate the importance of genetic identity of the main suspect in mortality, deformed wing virus (DWV), in driving honeybee loss. Using laboratory experiments and a systematic field survey, we demonstrate that an emerging DWV genotype (DWV-B) is more virulent than the established DWV genotype (DWV-A) and is widespread in the landscape. Furthermore, we show in a simple model that colonies infected with DWV-B collapse sooner than colonies infected with DWV-A. We also identify potential for rapid DWV evolution by revealing extensive genome-wide recombination in vivo The emergence of DWV-B in naive honeybee populations, including via recombination with DWV-A, could be of significant ecological and economic importance. Our findings emphasize that knowledge of pathogen genetic identity and diversity is critical to understanding drivers of species decline.

Interspecific competition in honeybee intracellular gut parasites is asymmetric and favours the spread of an emerging infectious disease.

There is increasing appreciation that hosts in natural populations are subject to infection by multiple parasite species. Yet the epidemiological and ecological processes determining the outcome of mixed infections are poorly understood. Here, we use two intracellular gut parasites (Microsporidia), one exotic and one co-evolved in the western honeybee (Apis mellifera), in an experiment in which either one or both parasites were administered either simultaneously or sequentially. We provide clear evidence of within-host competition; order of infection was an important determinant of the competitive outcome between parasites, with the first parasite significantly inhibiting the growth of the second, regardless of species. However, the strength of this 'priority effect' was highly asymmetric, with the exotic Nosema ceranae exhibiting stronger inhibition of Nosema apis than vice versa. Our results reveal an unusual asymmetry in parasite competition that is dependent on order of infection. When incorporated into a mathematical model of disease prevalence, we find asymmetric competition to be an important predictor of the patterns of parasite prevalence found in nature. Our findings demonstrate the wider significance of complex multi-host-multi-parasite interactions as drivers of host-pathogen community structure.

A sting in the spit: widespread cross-infection of multiple RNA viruses across wild and managed bees.

Declining populations of bee pollinators are a cause of concern, with major repercussions for biodiversity loss and food security. RNA viruses associated with honeybees represent a potential threat to other insect pollinators, but the extent of this threat is poorly understood. This study aims to attain a detailed understanding of the current and ongoing risk of emerging infectious disease (EID) transmission between managed and wild pollinator species across a wide range of RNA viruses. Within a structured large-scale national survey across 26 independent sites, we quantify the prevalence and pathogen loads of multiple RNA viruses in co-occurring managed honeybee (Apis mellifera) and wild bumblebee (Bombus spp.) populations. We then construct models that compare virus prevalence between wild and managed pollinators. Multiple RNA viruses associated with honeybees are widespread in sympatric wild bumblebee populations. Virus prevalence in honeybees is a significant predictor of virus prevalence in bumblebees, but we remain cautious in speculating over the principle direction of pathogen transmission. We demonstrate species-specific differences in prevalence, indicating significant variation in disease susceptibility or tolerance. Pathogen loads within individual bumblebees may be high and in the case of at least one RNA virus, prevalence is higher in wild bumblebees than in managed honeybee populations. Our findings indicate widespread transmission of RNA viruses between managed and wild bee pollinators, pointing to an interconnected network of potential disease pressures within and among pollinator species. In the context of the biodiversity crisis, our study emphasizes the importance of targeting a wide range of pathogens and defining host associations when considering potential drivers of population decline.

Draft Genome of the Entomopathogenic Fungus Metarhizium robertsii DSM 1490.

Metarhizium robertsii DSM 1490 is a generalist entomopathogenic fungus. The mechanisms of pathogenesis of such fungi in insects like termites are not completely understood. Here, we report the draft genome sequence, as sequenced on the Oxford Nanopore platform. The genome has a GC% of 47.82 and a size of 45,688,865 bp.

The ultrastructure of the rostral gland in soldiers of Verrucositermes tuberosus (Blattodea: Termitidae: Nasutitermitinae).

The soil-feeding habit is an evolutionary novelty found in some advanced groups of termites. The study of such groups is important to revealing interesting adaptations to this way-of-life. The genus Verrucositermes is one such example, characterized by peculiar outgrowths on the head capsule, antennae and maxillary palps, which are not found in any other termite. These structures have been hypothesized to be linked to the presence of a new exocrine organ, the rostral gland, whose structure has remained unexplored. We have thus studied the ultrastructure of the epidermal layer of the head capsule of Verrucositermes tuberosus soldiers. We describe the ultrastructure of the rostral gland, which consists of class 3 secretory cells only. The dominant secretory organelles comprise rough endoplasmic reticulum and Golgi apparatus, which provide secretions delivered to the surface of the head, likely made of peptide-derived components of unclear function. We discuss a possible role of the rostral gland of soldiers as an adaptation to the frequent encounter with soil pathogens during search for new food resources.

Defensive behavior is linked to altered surface chemistry following infection in a termite society.

The care-kill response determines whether a sick individual will be treated or eliminated from an insect society, but little is known about the physiological underpinnings of this process. We exploited the stepwise infection dynamics of an entomopathogenic fungus in a termite to explore how care-kill transitions occur, and identify the chemical cues behind these shifts. We found collective responses towards pathogen-injected individuals to vary according to severity and timing of pathogen challenge, with elimination, via cannibalism, occurring sooner in response to a severe active infection. However, injection with inactivated fungal blastospores also resulted in increased albeit delayed cannibalism, even though it did not universally cause host death. This indicates that the decision to eliminate an individual is triggered before pathogen viability or terminal disease status has been established. We then compared the surface chemistry of differently challenged individuals, finding increased amounts of long-chained methyl-branched alkanes with similar branching patterns in individuals injected with both dead and viable fungal blastospores, with the latter showing the largest increase. This coincided with the highest amounts of observed cannibalism as well as signs of severe moribundity. Our study provides new mechanistic insight into the emergent collective behaviors involved in the disease defense of a termite society.

Tripartite interactions: how immunity, microbiota and pathogens interact and affect pathogen virulence evolution.

The bipartite interactions between insect hosts and their bacterial gut microbiota, or their bacterial pathogens, are empirically and theoretically well-explored. However, direct, and indirect tripartite interactions will also likely occur inside a host. These interactions will almost certainly affect the trajectory of pathogen virulence evolution, an area that is currently under researched. The interactions within tripartite associations can be competitive, that is, exploitative-competition, interference-competition or apparent-competition. Competitive interactions will be significantly influenced by non-competitive effects, for example, immunopathology, immunosuppression, and microbiota-mediated tolerance. Considering a combination of these interactions and effects, will enable an increased understanding of the evolution of pathogen virulence. This new perspective allows us to identify several novel research questions, which we hope will be a useful framework for future research.

Eating in a losing cause: limited benefit of modified macronutrient consumption following infection in the oriental cockroach Blatta orientalis.

BACKGROUND: Host-pathogen interactions can lead to dramatic changes in host feeding behaviour. One aspect of this includes self-medication, where infected individuals consume substances such as toxins or alter their macronutrient consumption to enhance immune competence. Another widely adopted animal response to infection is illness-induced anorexia, which is thought to assist host immunity directly or by limiting the nutritional resources available to pathogens. Here, we recorded macronutrient preferences of the global pest cockroach, Blatta orientalis to investigate how shifts in host macronutrient dietary preference and quantity of carbohydrate (C) and protein (P) interact with immunity following bacterial infection. RESULTS: We find that B. orientalis avoids diets enriched for P under normal conditions, and that high P diets reduce cockroach survival in the long term. However, following bacterial challenge, cockroaches significantly reduced their overall nutrient intake, particularly of carbohydrates, and increased the relative ratio of protein (P:C) consumed. Surprisingly, these behavioural shifts had a limited effect on cockroach immunity and survival, with minor changes to immune protein abundance and antimicrobial activity between individuals placed on different diets, regardless of infection status. CONCLUSIONS: We show that cockroach feeding behaviour can be modulated by a pathogen, resulting in an illness-induced anorexia-like feeding response and a shift from a C-enriched to a more P:C equal diet. However, our results also indicate that such responses do not provide significant immune protection in B. orientalis, suggesting that the host's dietary shift might also result from random rather than directed behaviour. The lack of an apparent benefit of the shift in feeding behaviour highlights a possible reduced importance of diet in immune regulation in these invasive animals, although further investigations employing pathogens with alternative infection strategies are warranted.

Fast-forwarding evolution-Accelerated adaptation in a proofreading-deficient hypermutator herpesvirus.

Evolution relies on the availability of genetic diversity for fitness-based selection. However, most deoxyribonucleic acid (DNA) viruses employ DNA polymerases (Pol) capable of exonucleolytic proofreading to limit mutation rates during DNA replication. The relative genetic stability produced by high-fidelity genome replication can make studying DNA virus adaptation and evolution an intensive endeavor, especially in slowly replicating viruses. Here, we present a proofreading-impaired Pol mutant (Y547S) of Marek's disease virus that exhibits a hypermutator phenotype while maintaining unimpaired growth in vitro and wild-type (WT)-like pathogenicity in vivo. At the same time, mutation frequencies observed in Y547S virus populations are 2-5-fold higher compared to the parental WT virus. We find that Y547S adapts faster to growth in originally non-permissive cells, evades pressure conferred by antiviral inhibitors more efficiently, and is more easily attenuated by serial passage in cultured cells compared to WT. Our results suggest that hypermutator viruses can serve as a tool to accelerate evolutionary processes and help identify key genetic changes required for adaptation to novel host cells and resistance to antiviral therapy. Similarly, the rapid attenuation achieved through adaptation of hypermutators to growth in cell culture enables identification of genetic changes underlying attenuation and virulence, knowledge that could practically exploited, e.g. in the rational design of vaccines.

Cryptic diversity and female host specificity in a parasitoid where the sexes utilize hosts from separate orders.

Investigating complex parasitic life cycles is important for understanding the major fitness components that drive the evolution of host-parasite systems. The rare condition of heterotrophic heteronomy, in which the sexes utilize disparate host taxa, is a poorly understood complex parasitic lifestyle. One of only two known examples occurs in the Myrmecolacidae, an unusual family of the parasitoid order Strepsiptera (Insecta), in which males parasitize ants while females parasitize grasshoppers, crickets, and praying mantids. Here, we reconstruct the evolutionary pattern and timescale of host-use in a set of morphologically cryptic myrmecolacid taxa currently identified as Caenocholax fenyesi. We find that (i) C. fenyesi contains at least ten cryptic lineages consistent with separate species; (ii) Fossil evidence suggests a very low molecular clock rate and an ancient origin for cryptic lineages; (iii) Diversity among Caenocholax species is partitioned by geography and host association of the female; and (iv) Switches in host usage are uncoupled between the sexes, with changes in female host preference accompanying diversification. This study represents the first phylogeographical analysis of any strepsipteran, and the first molecular examination of host-use for a heterotrophic heteronomous taxon. Our results have implications for the understanding of evolution, host usage and estimated species richness in parasitic taxa.

High-Throughput Cultivation for the Selective Isolation of Acidobacteria From Termite Nests.

Microbial communities in the immediate environment of socialized invertebrates can help to suppress pathogens, in part by synthesizing bioactive natural products. Here we characterized the core microbiomes of three termite species (genus Coptotermes) and their nest material to gain more insight into the diversity of termite-associated bacteria. Sampling a healthy termite colony over time implicated a consolidated and highly stable microbiome, pointing toward the fact that beneficial bacterial phyla play a major role in termite fitness. In contrast, there was a significant shift in the composition of the core microbiome in one nest during a fungal infection, affecting the abundance of well-characterized Streptomyces species (phylum Actinobacteria) as well as less-studied bacterial phyla such as Acidobacteria. High-throughput cultivation in microplates was implemented to isolate and identify these less-studied bacterial phylogenetic group. Amplicon sequencing confirmed that our method maintained the bacterial diversity of the environmental samples, enabling the isolation of novel Acidobacteriaceae and expanding the list of cultivated species to include two strains that may define new species within the genera Terracidiphilus and Acidobacterium.

The mitochondrial genome of the 'twisted-wing parasite' Mengenilla australiensis (Insecta, Strepsiptera): a comparative study.

BACKGROUND: Strepsiptera are an unusual group of sexually dimorphic, entomophagous parasitoids whose evolutionary origins remain elusive. The lineage leading to Mengenilla australiensis (Family Mengenillidae) is the sister group to all remaining extant strepsipterans. It is unique in that members of this family have retained a less derived condition, where females are free-living from pupation onwards, and are structurally much less simplified. We sequenced almost the entire mitochondrial genome of M. australiensis as an important comparative data point to the already available genome of its distant relative Xenos vesparum (Family Xenidae). This study represents the first in-depth comparative mitochondrial genomic analysis of Strepsiptera. RESULTS: The partial genome of M. australiensis is presented as a 13421 bp fragment, across which all 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes and 18 transfer RNA (tRNA) sequences are identified. Two tRNA translocations disrupt an otherwise ancestral insect mitochondrial genome order. A+T content is measured at 84.3%, C-content is also very skewed. Compared with M. australiensis, codon bias in X. vesparum is more balanced. Interestingly, the size of the protein coding genome is truncated in both strepsipterans, especially in X. vesparum which, uniquely, has 4.3% fewer amino acids than the average holometabolan complement. A revised assessment of mitochondrial rRNA secondary structure based on comparative structural considerations is presented for M. australiensis and X. vesparum. CONCLUSIONS: The mitochondrial genome of X. vesparum has undergone a series of alterations which are probably related to an extremely derived lifestyle. Although M. australiensis shares some of these attributes; it has retained greater signal from the hypothetical most recent common ancestor (MRCA) of Strepsiptera, inviting the possibility that a shift in the mitochondrial selective environment might be related to the specialization accompanying the evolution of a small, morphologically simplified completely host-dependent lifestyle. These results provide useful insights into the nature of the evolutionary transitions that accompanied the emergence of Strepsiptera, but we emphasize the need for adequate sampling across the order in future investigations concerning the extraordinary developmental and evolutionary origins of this group.

Exposure of Larvae of the Solitary Bee Osmia bicornis to the Honey Bee Pathogen Nosema ceranae Affects Life History.

Wild bees are important pollinators of wild plants and agricultural crops and they are threatened by several environmental stressors including emerging pathogens. Honey bees have been suggested as a potential source of pathogen spillover. One prevalent pathogen that has recently emerged as a honey bee disease is the microsporidian Nosema ceranae. While the impacts of N. ceranae in honey bees are well documented, virtually nothing is known about its effects in solitary wild bees. The solitary mason bee Osmia bicornis is a common pollinator in orchards and amenable to commercial management. Here, we experimentally exposed larvae of O. bicornis to food contaminated with N. ceranae and document spore presence during larval development. We measured mortality, growth parameters, and timing of pupation in a semi-field experiment. Hatched individuals were assessed for physiological state including fat body mass, wing muscle mass, and body size. We recorded higher mortality in the viable-spore-exposed group but could only detect a low number of spores among the individuals of this treatment. Viable-spore-treated individuals with higher head capsule width had a delayed pupation start. No impact on the physiological status could be detected in hatched imagines. Although we did not find overt evidence of O. bicornis infection, our findings indicate that exposure of larvae to viable N. ceranae spores could affect bee development.

Evolution of Termite Symbiosis Informed by Transcriptome-Based Phylogenies.

Termitidae comprises ∼80% of all termite species [1] that play dominant decomposer roles in tropical ecosystems [2, 3]. Two major events during termite evolution were the loss of cellulolytic gut protozoans in the ancestor of Termitidae and the subsequent gain in the termitid subfamily Macrotermitinae of fungal symbionts cultivated externally in "combs" constructed within the nest [4, 5]. How these symbiotic transitions occurred remains unresolved. Phylogenetic analyses of mitochondrial data previously suggested that Macrotermitinae is the earliest branching termitid lineage, followed soon after by Sphaerotermitinae [6], which cultivates bacterial symbionts on combs inside its nests [7]. This has led to the hypothesis that comb building was an important evolutionary step in the loss of gut protozoa in ancestral termitids [8]. We sequenced genomes and transcriptomes of 55 termite species and reconstructed phylogenetic trees from up to 4,065 orthologous genes of 68 species. We found strong support for a novel sister-group relationship between the bacterial comb-building Sphaerotermitinae and fungus comb-building Macrotermitinae. This key finding indicates that comb building is a derived trait within Termitidae and that the creation of a comb-like "external rumen" involving bacteria or fungi may not have driven the loss of protozoa from ancestral termitids, as previously hypothesized. Instead, associations with gut prokaryotic symbionts, combined with dietary shifts from wood to other plant-based substrates, may have played a more important role in this symbiotic transition. Our phylogenetic tree provides a platform for future studies of comparative termite evolution and the evolution of symbiosis in this taxon.

A proofreading-impaired herpesvirus generates populations with quasispecies-like structure.

RNA virus populations are composed of highly diverse individuals that form a cloud of related sequences commonly referred to as a 'quasispecies'1-3. This diversity arises as a consequence of low-fidelity genome replication4,5. By contrast, DNA virus populations contain more uniform individuals with similar fitness6. Genome diversity is often correlated with increased fitness in RNA viruses, while DNA viruses are thought to require more faithful genome replication. During DNA replication, erroneously incorporated bases are removed by a 3'-5' exonuclease, a highly conserved enzymatic function of replicative DNA but not RNA polymerases. This proofreading process enhances replication fidelity and ensures the genome integrity of DNA organisms, including large DNA viruses7. Here, we show that a herpesvirus can tolerate impaired exonucleolytic proofreading, resulting in DNA virus populations, which, as in RNA viruses8, are composed of highly diverse genotypes of variable individual fitness. This indicates that herpesvirus mutant diversity may compensate for individual fitness loss. Notably, in vivo infection with diverse virus populations results in a marked increase in virulence compared to genetically homogenous parental virus. While we cannot exclude that the increase in virulence is caused by selection of and/or interactions between individual genotypes, our findings are consistent with quasispecies dynamics. Our results contrast with traditional views of DNA virus replication and evolution, and indicate that a substantial increase in population diversity can lead to higher virulence.