Skip to main content

Advertisement

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

Article metrics

Abstract

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.

Background

The importance of microbes for the evolution of higher organisms is starting to be realized [1, 2]. Metazoan evolution is not only driven by pathogenic microbes, as reflected by fast evolution of immune genes [3]. Rather, microbes often are facilitators of metabolic and environmental adaptations [2, 4, 5]. For instance, the gut microbial communities of wood-feeding roaches and termites facilitate thriving on a wood diet that is difficult to digest and poor in nitrogen. Nitrogen fixation and the digestion of wood depend on the termite gut microbiome [2, 6, 7]. In lower termites, lignocellulose degradation was initially mainly attributed to unicellular eukaryotes (protists) in the gut [8]. Recently, it has become evident that lignocellulose degradation is a synergistic effort of the termite, its associated protists, and bacteria [9,10,11]. In addition to their role in lignocellulose degradation, bacteria are also essential for the assimilation of nitrogen taken up from the environment. Nitrogen can be acquired from the environment either via fixation from the atmosphere [12, 13], or via nitrate reduction [14]. Also, nitrogen can be recycled from the metabolic waste product uric acid [15, 16]. By using genome sequencing and pathway reconstruction, these processes have been assigned to four major bacterial phyla in the termite gut: Proteobacteria (Desulfovibrio [17]), Spirochetes (Treponema [18, 19]), Bacteroidetes (Azobacteroides) [16], and Elusimicrobia (Endomicrobium [20, 21]).

Many bacteria in the termite gut live in tight association with protists, where they sit on the surface [22, 23], in invaginations of the cell membrane [17], or even inside the protist cells [24]. Such tight associations lead to frequent vertical transmission of bacteria between protist generations. In return, protists and bacteria are vertically transmitted between termite generations via proctodeal trophallaxis during colony foundation [25]. Vertical transmission has led to co-speciation between bacteria and their protist hosts, and sometimes even the termite hosts [26,27,28,29]. Evidence for horizontal transfer of protists between termite species, so called transfaunations, is limited to a few exceptions [30]. Hence, the termite host species association is rather strict, leading to strong phylogenetic imprints on protist community structure [31,32,33]. In comparison, the bacterial microbiome is more flexible, frequently transferred between termite host species [34], and affected by diet [33, 35,36,37,38,39,40,41].

There is evidence that the gut microbiome of termites has contributed to the adaptation of different termite species to their specific ecologies [33, 36, 42,43,44]. There are pronounced ecological differences between the so-called termite life types [45, 46]. Termite species of the wood-dwelling life type never leave their nest except for the mating flight. They feed on a relatively uniform bonanza resource, that is the piece of wood they built their nest in [47, 48]. On the other hand, foraging species leave their nest to forage for food and have access to additional nutrients [47, 49]. This likely imposes different selection pressures on the termite holobiont, in particular with regard to nutrient uptake. Because the microbiome is directly involved in nutrient uptake, it seems reasonable to hypothesize that it may also play a role in adaptation to life type related dietary differences. In this scenario, one would expect the life types to leave an imprint on microbiome structure and function. As such, searching for microbial imprints of a given life type can possibly provide us with a lead for microbiome-mediated adaptation.

One potential pitfall of such an endeavor is that microbiomes may bear imprints from transient microbes that were ingested from the environment. Transient microbes rarely form evolutionary relevant relationships with the host [50, 51]. Instead, they reflect short-term associations with microbes from the local environment the termites were collected from. For example, in previous studies [34, 36,37,38,39,40,41] microbiomes were profiled from termites that were collected from a variety of substrates and locations, sometimes including natural and laboratory environments in the same study. In such an experimental setup, microbiomes can be similar because termites ingested transient microbes from similar environments. For instance, microbiomes from different termite species may become similar because the termites ingested soil that harbors a discrete microbiome. In such a case it is challenging to disentangle effects stemming from ingested, transient soil microbes from those that resulted from potential adaptations of the holobiont to a soil diet. Critically, such similarities may then erroneously interpreted as involved in the adaptation to a soil diet. Similarly, rearing termites in different laboratories without acclimation to common conditions in the same laboratory could lead to artificial clustering of microbiomes. Strong laboratory effects have been shown for the microbiome of Drosophila [52]. Therefore, it is essential to reduce the impact of transient microbes in experimental setups to test evolutionary hypotheses.

In order to explore potential ecological imprints on the microbiome, we focused on an evolutionary switch between wood-dwelling and foraging life types in the Rhinotermitidae (Fig. 1). Reticulitermes species are of the foraging life type, while Prorhinotermes simplex is wood-dwelling. If the microbiome was affected by life type specific ecology, we would expect that the microbiome of Prorhinotermes simplex was similar to that of the other wood-dwelling species (Cryptotermes) although these are from a different family (Kalotermitidae). At the same time, the microbiome of the foraging Reticulitermes species should bear distinct features. Alternatively, if there was no ecological imprint, we would expect the microbiome to follow a phylogenetic pattern, with the Rhinotermitidae Prorhinotermes and Reticulitermes forming a cluster and the Cryptotermes species (Kalotermitidae) forming a second cluster. Using this experimental setup, we recently showed that protist community composition aligned with phylogeny, but bacterial communities aligned more strongly with wood-dwelling and foraging life types [33].

Fig. 1
figure1

Schematic phylogeny of the five lower termite species used in this study from [33]. Branch length not drawn to scale. Colored boxes indicate the life type

To further explore this, we investigated whether changes in microbiome composition are also reflected by changes in microbiome function, as would be expected if the microbiome played a role in adaptation. For instance, we would expect dietary adaptations to be reflected by changes to pathways involved in substrate breakdown and effective provisioning of limiting nutrients such as nitrogen. In order to test whether and which changes in the functional repertoire align with life type and could be involved in potential adaptation to different diets, we characterized the metagenome of two foraging species; Reticulitermes flavipes and Reticulitermes grassei. We compared their functional repertoire to that of three wood-dwelling species Prorhinotermes simplex, Cryptotermes secundus, and Cryptotermes domesticus. Because there can be substantial variation in microbial communities between colonies [53,54,55,56], we analyzed five C. domesticus, eight C. secundus, seven P. simplex, five R. flavipes, and four R. grassei replicate colonies. We focused on the persistent, long-term differences between microbiomes by controlling short-term effects caused by the influx of transient microbes. This was achieved by feeding a common diet of sterile Pinus wood for several weeks prior to sample collection.

Results

We analyzed a total of ~ 440 million metagenomic shotgun sequences. Between 974,176 and 8,949,734 sequences per sample were of microbial origin (Additional file 2: Table S1). Sequences were subsampled (rarefied) to 1,386,882 bacterial and 2781 protist annotated sequences per sample. For annotation, sequences were aligned to a reference database of clusters of orthologous groups of genes (COGs) with known function. These COGs represent the lowest level of the evolutionary genealogy of genes - non-supervised orthologous groups (eggNOG) hierarchical annotation. At the next higher level, the COGs are grouped into pathways (Additional file 1: Fig. S2, Fig. S3), and at the third and highest level, the pathways are grouped into three categories; “information storage and processing”, “cellular process and signaling”, and “metabolism”. We adhere to this definition of the eggNOG hierarchical terms throughout the study.

“information storage and processing” differentiates the protist metagenomes of wood-dwelling and foraging lower termite species

In our previous study [33] on the identical samples, the protist communities of the Rhinotermitidae Prorhinotermes and Reticulitermes clustered together, supporting a phylogenetic imprint on community composition. Here, we tested whether this pattern was also reflected by the functions encoded by the protist metagenome. Therefore, we annotated metagenome encoded functions in the shotgun sequences and compared the functional metagenome profiles across host species, using Bray-Curtis-Dissimilarity [57]. This index considers abundance of functional categories, thus avoiding arbitrary coverage cutoffs.

The protist functional repertoire clustered according to host family and genus (Fig. 2a), thus showing a dominant phylogenetic imprint. Family-wise clustering was supported by Redundancy Analysis (RDA): the model including host family explained more variance in the functional repertoire and produced lower Akaike’s Information Criteria (AICs) than the model based on life type (Table 1). For a more detailed view, we analyzed the three categories at the highest level in the eggNOG hierarchical annotation (Additional file 1: Fig. S2) separately. Cluster analysis of the categories “cellular process and signaling “and “metabolism” supported the notion that phylogenetic relatedness is an important factor for functional similarity (Additional file 1: Figs. S4B and D). In contrast, the portion of the metagenome assigned to “information storage and processing” (Fig. 2b) clustered primarily by life type. The stronger effect of life type than phylogeny on this functional category was also supported by higher explanatory power and lower AICs in RDA (Table 1).

Fig. 2
figure2

Cluster dendrograms of the functional profiles of the protist and bacterial community. Community distances are based on Bray-Curtis Dissimilarities of a) all functions of the protist community (25,795 sequences), b) The category “information storage and processing” of the protist community (4527 sequences), c) all functions of the bacterial community (21,215,480 sequences), and d) the category “metabolism” of the bacterial community (10,586,058 sequences). Cd (red) = C. domesticus colonies; Cs (orange) = C. secundus colonies; Ps (green) = P. simplex colonies; Rf (blue) = R. flavipes colonies; Rg (lightblue) = R. grassei colonies. Green background = wood-dwelling life types; orange background = foraging life type. For protist functions involved in “information storage and processing” in the protist community, samples clustered according to life type. Similarly, the bacterial metabolic metagenomes clustered according to life type. Cluster dendrograms of all functional categories of protist and bacterial communities can be found in Additional file 1: Fig. S4 and S7

Table 1 Models of the effects of life type and host family (phylogeny) on functional community profiles

Identification of the functions that differentiate the protist metagenomes of the wood-dwelling and foraging species can hold clues with regard to the nature of potentially adaptive phenotypes in the protist metagenome. In order to do so, we performed a linear discriminant analysis (LEfSe: [58]). This analysis identified 22 over-represented COGs in foraging and 14 in wood-dwelling species (Fig. 3a, Additional file 3: Table S2, p < 0.05, q < 0.05, LDA score > 2, Additional file 1: Figure S6).

Fig. 3
figure3

Differences in the functional content of the protist metagenomes of wood-dwelling and foraging species. a Circular dendrogram/hierarchy of all over-represented COGs in the category “information storage and processing” in wood-dwelling species (green) or foraging species (orange). Circle size at the edges scales with abundance of the COG. Colored branches indicate over-represented pathways. Over-representation was detected with LEfSe [58] (p < 0.05, q < 0.05, LDA > 2). A Venn diagram visualizing the total number, and the differentially abundant number of functions in each of the five pathways that constitute the category “information storage and processing” can be found in Additional file 1: Fig. S6. b Sequence coverage of wood-dwelling (green) and foraging (orange) species of examples of over-represented COGs mentioned in the text. Error bars represent 95% confidence intervals across replicate colonies

The pathway “replication, recombination, and repair” was over-represented in the foraging species (Fig. 3a, Additional file 3: Table S2, p = 0.0001, q = 0.002). The over-represented COGs in this pathway included a DNA dependent DNA-polymerase (COG0470) and five helicases (COG0514, COG0553, COG1199, COG1204, ENOG410XNUT, see Fig. 3b for grouped analysis and Additional file 3: Table S2 for individual COG p- and q-values). In wood-dwelling species, the pathway “transcription” was over-represented (p = 0.0004, q = 0.003). The over-represented COGs in this pathway contained DNA binding domains and were supposedly involved in transcriptional regulation (COG5147, ENOG4111SAB).

The bacterial metabolic metagenome aligns with host ecology

In our previous study [33], the bacterial community composition of termite hosts clustered primarily by life type, which is consistent with ecology related differences between microbiomes. Following the rationale above, we tested whether this pattern was also reflected by the functions encoded by the metagenome.

Against the expectation from our previous study, the functional bacterial profiles showed no life type, but a phylogenetic imprint, which is in line with the protist functional profiles. Most samples clustered according to host family (Fig. 2c). Analyzing the three high-level eggNOG functional categories separately provided more detailed insight. The categories “cellular process and signaling “and “information storage and processing” supported the notion of strong phylogenetic effects on metagenome function (Additional file 1: Figs. S7B and C). In contrast, the metabolic metagenomes (Fig. 2d) clustered primarily according to host life type. Host life type was also a better predictor for metabolic functions than host family in RDA (Table 1).

Aside from these general patterns, several samples stood out. Samples Rg2 and Rg4 of R. grassei were on long branches in the dendrograms (Fig. 2 and Additional file 1: Figure S7), suggesting unusual functional profiles. Notably, these samples already stood out in our previous study [33] because of their unusual abundance of microbial taxa potentially due to infection with pathogens. This unusual composition was confirmed by taxonomic annotation in this study (see Additional file 1: Figure S8). Sample Cs7 (C. secundus) also clustered separately from the other samples. This was mainly driven by abundant transposases in this sample (53.1% of sequences) (for example COG1662, COG3385, or ENOG410XT1T, see Additional file 4: Table S3), accompanied by an increase in the frequency of Bacteroides (Additional file 1: Fig. S8) that are rich in conjugative transposons [59, 60]. We performed all analyses with and without these samples and found no qualitative differences (data not shown).

Bacterial metabolic functions that differentiated wood-dwelling from foraging species were identified using linear discriminant analysis (LDA). 105 metabolic COGs were over-represented in the wood-dwelling species, while 151 were over-represented in the foraging species (Additional file 5: Table S4, p < 0.05, q < 0.05, LDA score > 2, Additional file 1: Figure S11). All COGs described as over-represented or enriched in the following were subject to these p-value, q-value and LDA cutoffs. Because of their specialized diet, genes involved in nitrogen metabolism and lignocellulose break down like glycoside hydrolases (GH) are of particular interest, when focusing on metabolic differences among gut microbiomes in wood-feeding termites with different diets. In fact, among the genes involved in ‘carbohydrate transport and metabolism’ that were enriched in the microbiome of wood-dwelling termites, GHs were over-represented (43.3% of enriched genes versus 12% expected, exact binomial test, p = 2.12e-05, Additional file 5: Table S4, Additional file 6: S5). In the foraging termite species, only one gene with putative lignocellulolytic activity was over-represented (COG3858), suggesting that the wood-dwelling species have a higher potential for complex carbohydrate degradation. To further investigate differences in GH abundance between the microbiomes of wood-dwelling and foraging species, we performed a detailed pathway analysis using the carbohydrate-active enzymes (CAZy) database ([61], Fig. 4). All GHs acting in hemicellulose break-down were more abundant in the wood-dwelling species (Fig. 4b). Among the cellulolytic enzymes, ß-glucosidases were significantly more abundant in the wood-dwelling species. The other two enzymes involved (cellulase (endo-ß-1.4-glucanase), cellobiohydrolase) showed a trend into the same direction. All of the genes with cellulolytic or hemicelllulolytic activity were affiliated with Bacteroidetes (mostly members of the genus Bacteroides) or the genus Treponema. Additional support for the increased importance of hemicellulose utilization in the wood-dwelling species was provided by the over-representation of twelve COGs annotated as TonB-dependent receptors (ENOG410XNNV, ENOG410XNPQ or COG4206, see Additional file 5: Table S4). Apart from other substrates, these receptors are important for the uptake of plant-derived hemicellulose [62, 63]. All functions annotated as TonB-dependent receptors (or TonB-dependent associated receptor plugs) were affiliated with the genus Bacteroides (see Additional file 5: Table S4).

Fig. 4
figure4

Differences in the functional content of the bacterial metagenomes of wood-dwelling and foraging species. a Circular dendrogram/hierarchy of all COGs in the category “metabolism” over-represented in wood-dwelling species (green) or foraging species (orange). Circle size at the leafs scales with the abundance of the COG. Over-representation was detected with LefSe [58] (p < 0.05, q < 0.05, LDA > 2). A Venn diagram visualizing the total number, and the differentially abundant number of functions in each of the five pathways that constitute the category “metbolism” can be found in Additional file 1: Fig. S11. b Pathway analysis of cellulose and hemicellulose degradation. Colored boxes of cellulolytic or hemicellulolytic genes indicate proportion of relative abundance of sequences affiliated with wood-dwelling (green) or foraging (orange) species. c Pathway analysis of nitrogen metabolism. Boxes for genes with functions in nitrogen metabolism indicate relative abundance in the two life types. d Pathway analysis of glutamate synthesis. Boxes in c) and d) show relative abundance in the two life types of genes with functions in nitrogen/glutamate metabolism. Pie charts show taxonomic association of the gene. All hemicellulolytic genes were over-represented in the wood-dwelling species. Also, a nitrogenase was enriched in the wood-dwelling species, while in the foraging species, genes involved in dissimilatory nitrate reduction were over-represented

Because wood is poor in nitrogen, termites depend on an efficient system for conserving and upgrading nitrogen [6]. In the wood-dwelling species, a potential nitrogenase (nifH, COG1348) was over-represented (Fig. 4c, Additional file 5: Table S4). Nitrogenases are key enzymes in the fixation of atmospheric nitrogen and downstream ammonia synthesis. Nitrogenase activity was mainly affiliated with members of the genus Treponema (Fig. 4c). In contrast, in the foraging species, COGs involved in dissimilatory nitrate reduction (COG1251, COG5013, COG2181, COG0243, Fig. 4c, Additional file 5: Table S4) were over-represented. They were affiliated with a variety of different genera ranging from Desulfovibrio and Gordonibacter to Stenoxybacter, Enterobacter and Serratia. Serratia and Enterobacter are potential insect pathogens and contributed to the prevalence of one of the three nitrate reductases, narG (COG5013). Closer inspection of the source of these bacteria revealed that they mainly stemmed from the abnormal samples Rg2 and Rg4 that we suspected to carry a potential pathogenic infection. When we remove these samples from the analysis the increase of narG in foragers remains significant (p = 0.034).

For living on a nitrogen poor substrate, it can also be adaptive to effectively recycle nitrogen from the main waste product of the host’s amino acid metabolism, uric acid. Uric acid can be recycled through anaerobic ammonia production and downstream glutamate synthesis [6, 15, 20, 64]. In the wood-dwelling species a putative glutamate dehydrogenase (COG0334), involved in glutamate synthesis by ammonia assimilation, was over-represented. This glutamate dehydrogenase gene was mainly affiliated with members of the genera Bacteroides, Treponema and Desulfovibrio. In the foraging species, COGs with putative glutamine (COG0174) and glutamate synthase (COG0067, COG0069) function, were enriched (Fig. 4d). These COGs were affiliated with Desulfovibrio, Treponema, Pseudomonas and Acetobacterium.

Discussion

In this study, we assessed functional differences of termite metagenomes that underwent an evolutionary switch from wood-dwelling to foraging to identify putative contributions of the microbiome to ecological niche adaptation. To do this, we chose a set of five termite species (two foraging, three wood-dwelling species) and determined whether the functional profiles of the termite gut microbiome followed phylogeny of the host or aligned with host ecology. We hypothesized that alignment of microbiome function with termite life type is consistent with a contribution of the microbiome to termite holobiont adaptation to different diets. By comparing the functional content of microbiomes of different host species we focused on long-term evolutionary processes.

A potential pitfall of such an approach is that an alignment of the termite microbiome with life type-related ecology could also be caused by short-term differences between microbiomes that are merely transient. For example, microbes in the environment might differ between collection sites for the different host species. Further, ingestion of environmental microbes might lead to an association between microbiome and ecology. Similarly, differences in local food supply can lead to transient, short-term effects on the termite microbiome [56]. Consequently, such short-term differences reflect environmental differences at termite collection sites, rather than potentially adaptive, evolved differences between host-lineage-specific microbiomes.

For this reason, we chose to follow an approach where we control for environmental and dietary differences by acclimating all termites on the same (sterile) food source and to identical environmental conditions. We consider metagenomic patterns that persist under such highly controlled experimental conditions as robust and indicative of long-term, evolutionary acquired differences, rather than short-term imprints originating from differences in the environment or food source. It should be noted that the experimental setup poses a restriction to the number of sampled host species [33].

Increased potential for replication in the protists of foraging termite species

In the protist metagenome of foraging species, genes involved in replication were more abundant. High replication rates are expected to be more frequently under positive selection during recolonization of the gut with protists, when the gut environment has not yet reached carrying capacity [65]. Therefore, we would like to speculate that this difference is related to the fact that Reticulitermes guts have to be recolonized more frequently because they molt more frequently; the intermolt periods in Reticulitermes are about two weeks long [49], while they average almost two months in Cryptotermes [48]. During molting the protists are lost and the guts have to be recolonized through proctodeal trophallaxis from nest mates [66]. However, we are aware that differences in the relative abundance of housekeeping genes like those required for replication between protist microbiomes can not be clearly disentangled from differences in average protist genome size and therefore should be interpreted with caution.

Enrichment of genes for lignocellulose degradation in the microbiome of wood-dwelling termite species

While genes involved in replication differentiated the protist metagenomes of wood-dwelling and foraging species in our study, metabolic genes differentiated the bacterial metagenomes. Consistent with differences in their respective diets, the metagenomes of foraging and wood-dwelling species in our study differed by their potential for cellulose and hemicellulose utilization. Several GHs that have cellulolytic and hemicellulolytic function were over-represented in the metagenomes of wood-dwelling species (GH families 2, 3, 16, 43, mannosidases, xylosidases, glucanases, xylanases, Fig. 4b, Additional file 5: Table S4). A more detailed pathway analysis confirmed that hemicellulases are more abundant in the wood-dwelling species. This suggests a more pronounced role for lignocellulose degradation in the metabolism of the wood-dwelling species in our study. Accordingly, TonB dependent transporters were enriched in the microbiome of wood-dwellers. These transporters can shuttle hemicellulose and its building-blocks, in particular xylans and xylose through bacterial membranes [67, 68]. A large fraction of cellulases, hemicellulases, and putative TonB transporters were attributed to the genus Bacteroides. In Bacteroides, TonB dependent transporters are often co-localized and co-regulated with enzymes for polysaccharide degradation like hemicellulases [60, 69]. This suggests a partnership of enzymes and transporters in polysaccharide degradation. Bacteroides species from the human gut are also hemicellulose degraders [70], suggesting a distinctive role for the genus in hemicellulose degradation in termites as well.

The above-identified differences in functional potential between the wood-dwelling and foraging species in our study are suggestive of adaptations to utilize diets that differ in hemicellulose content. Hemicellulose content differs between wood species [71, 72]. The wood-dwelling Cryptotermes species in our study are mostly found in hardwood mangroves [73] where they can thrive on a bonanza food resource. The other wood-dwelling genus in our study, Prorhinotermes, lives in similar coastal habitats with a similar arboreal flora [74]. Hardwood is richer in hemicelluloses and the potential to use hemicelluloses is larger in the microbiome of species living on hardwood. On the other hand, Reticulitermes species originated in inland habitats [75], prefer soft woods like pine [76, 77] with lower hemicellulose levels, and accordingly, hemicellulolytic pathways are depleted.

Termites with different life-types rely on different forms of nitrogen uptake and recycling

Nitrogen is scarce in a wood-based diet. As a consequence, termites need to acquire additional nitrogen from the environment. The microbiome is essential for this process. In the microbiome of wood-dwelling species, which feed on a uniform lignocellulose diet, a potential nitrogenase gene was enriched (nifH, COG1348). Nitrogenases are the key enzymes in the fixation of atmospheric nitrogen and downstream ammonia synthesis. This nifH was mainly affiliated with treponemes that have been shown to play an important role in nitrogen fixation before [12, 18, 19]. In contrast, the microbiome of the foraging species in our study has a higher potential to provide nitrogen to the termite holobiont by dissimilatory reduction of nitrate (Fig. 4c). Nitrogen in the form of nitrate naturally occurs in soil. R. flavipes has been shown to acquire micro nutrients from soil [78] and to actively balance mineral uptake by food choice [79]. Therefore, it seems reasonable to assume that the microbiome of Reticulitermes relies on nitrogen from soil in the form of nitrate to balance the low nitrogen content of wood. The necessary nitrate reductases were found primarily in Desulfovibrio, Gordonibacter and Stenoxybacter that were found in association with Reticulitermes before and are shared between a wide range of termites [33, 80, 81].

Aside from obtaining nitrogen from the environment (atmosphere, soil), bacteria can also recycle uric acid nitrogen. All of these processes result in ammonia synthesis, the central metabolite of nitrogen metabolism. Ammonia is then further assimilated to glutamate. In the wood-dwelling species a glutamate dehydrogenase (COG0334) was over-represented. It was mainly affiliated with members of the Bacteroides, Desulfovibrio and treponemes. The foraging species seem to rely on another glutamate synthesis pathway, including glutamine (COG0174) and glutamate synthases (COG0067, COG0069). Accordingly, they were associated with a different set of bacteria including Pseudomonas, Acetobacterium, Desulfovibrio, and treponemes (Fig. 4d).

Phylogeny and ecology align with metagenome-encoded functions

Differences in the propensity for nitrogen uptake and recycling are likely to reflect differences in diet of the termite host species. Given differences in diet between the species that represent the different life types, it seems also reasonable to suggest that the changes in the repertoire of hemicellulases reflects adaptations of the microbiome to diets with different hemicellulose content. The finding that this manifested specifically in the metabolic functional repertoire, may suggest that potential selection acts in particular on metabolic functions.

Metabolic microbiome mediated adaptation to different diets can happen in two ways. First, acquisition of new microbes with adaptive functions could lead to adaptive changes of the microbiome. Second, genome evolution of microbes that are already associated with the host could lead to adaptation. Microbes that were already present before the onset of lineage specific adaptation are likely to be shared among host species. By contrast, newly acquired microbes are expected to be host lineage specific. We found that the bacterial groups that contributed most to the differentiation of metabolic functions are shared among all five host species (Treponema, Bacteroides, Desulfovibrio, Dysgomonas, Gordonibacter, Pseudomonas, Additional file 5: Table S4, Additional file 1: Figure S8). This supports that genome evolution of microbes that were already associated with the host contributed to potential adaptation in our model system.

Conclusion

We applied metagenomic sequencing of gut microbiomes from a controlled experimental setup to assess a putative contribution of the microbiome to host dietary adaptation that accompanies the evolutionary switch from wood-dwelling to foraging life types. We found that the overall pattern of microbiome variation reflected a phylogenetic signal. Interestingly, specific functions of the microbiome aligned with the underlying host ecology. The specific ecology related differences in microbiome function led us to hypothesize that the microbiome contributed to dietary adaptations, namely different hemicellulose and nitrogen contents. This hypothesis can now be tested, assessing host fitness under different dietary conditions. Such experiments will be crucial to disentangle adaptive functional changes from selectively neutral functional turnover or side effects from other adaptations.

Methods

Termite samples

All termites were collected from typical natural habitats (see [33]). They were kept under constant conditions (27 °C, 70% humidity) on autoclaved Pinus radiata wood from the same source for at least six weeks prior to the experiment. The feeding of Pinus represents a natural or near natural treatment; Pinus is a natural food source of P. simplex and Reticulitermes. Cryptotermes growth and behavior on Pinus recapitulates that on natural substrate [73]. The time of the acclimation period was chosen to lie well beyond the gut passage time of 24 h in lower termites [82, 83] and following Huang et al. [84], who showed that six weeks are sufficient for the microbiota to adjust to a new diet. That way, all excretable material like remaining food, transient microbes taken up from the environment that have no mechanisms to persist in the gut, and microbial DNA taken up before the experiment was made sure to be excreted. The samples were identical to those analyzed in our previous study, [33] where detailed information about animal collection, keeping, and cytochrome oxidase II based species identification and a phylogeny can be found.

DNA extraction and shotgun sequencing

DNA was extracted from a pool of three worker guts per colony using bead beating, chloroform extraction and isopropanol precipitation (see Additional file 1: methods section S13). Each of the 29 colony samples went through independent metagenomic shotgun library preparation and sequencing on an Illumina HiSeq platform (150 bp paired end reads).

Analysis

We employed a double filtering strategy to remove host DNA from our analysis. First, sequences were removed that mapped to available host genomes from C. secundus [85] and transcriptomes from P. simplex [86] and R. flavipes, provided by the 1KITE consortium (www.1kite.org, BioSample SAMN04005235) using BBMap [87] (for detailed workflow and more detailed information about used genomes and transcriptomes see Additional file 1: Fig. S14 and section S15). Of note, the sequences were not assembled, but individual reads were directly annotated. In a second step we used taxonomic and functional annotations with Megan6 [88] to retrieve only sequences that could be unambiguously assigned to either bacteria or protists. In order to compare the bacterial and protist data sets of all samples, they were rarefied to the number of sequences in the sample with lowest coverage, resulting in 1,386,882 and 2781 sequences per sample, respectively. Sample Cs4 was excluded from the analysis for insufficient sequence coverage (974,176 sequences), so was Cs5 from the protist data. Sample Ps5 did not pass the analysis pipeline and was also excluded.

Functional annotation with the eggNOG database resulted in the highest number of annotated sequences (21,215,480 annotated sequences in total) and was chosen for further functional analysis. Bray-Curtis Distances of functional abundances were clustered with the pvClust package in R [89]. Multivariate modeling was performed via RDA (Redundancy Analysis) and AICs as well as values for the proportion of variance explained were derived with the model selection tool ordistep and ordiR2step, as implemented in the R vegan package [90]. Models were compared to the null-model via ANOVA. To identify over-represented functions associated with the two termite life types, a Linear Discriminant Analysis (LDR) was performed using LEfSe [58] and visualized using GraPhlAn [91]. Pathway analysis of CAZy GHs was performed by blasting bacterial reads of all samples against the full CAZy protein database, using Diamond [92]. GH abundance was estimated by counting reads with matches on proteins with cellulolytic and hemicellulolytic functions [93]. Pathway analysis of the nitrogen metabolism was performed by searching COG IDs corresponding to the KEGG IDs among the over-represented COGs from the LEfSe analysis. A detailed workflow for full reproducibility, including scripts can be found in the supplementary pdf file in Additional file 1: Figure S14 and in sections S13 and S15 in the same file.

Availability of data and materials

The raw data has been uploaded to the NCBI short-read archive (BioProject ID PRJNA509211, Accession: SAMN10573992 – SAMN10574019). Supporting information and analysis workflows are included in the Additional file in this article. For the analysis workflow see Additional file 1: Fig. S14 and sections S13 and S15 in the supplementary pdf.

Abbreviations

AIC:

Akaike’s Information Criterion

CAZy:

Carbohydrate-Active enZymes

COGs:

Clusters of Orthologous Groups of genes

eggNOG:

Evolutionary Genealogy of Genes: Non-supervised Orthologous Groups

GH:

Glycoside Hydrolase

LDA:

Linear Discriminant Analysis

RDA:

Redundancy Analysis

References

  1. 1.

    McFall-Ngai M, Hadfield MG, Bosch TCG, Carey HV, Domazet-Lošo T, Douglas AE, et al. Animals in a bacterial world, a new imperative for the life sciences. PNAS. 2013;110:3229–36.

  2. 2.

    Bang C, Dagan T, Deines P, Dubilier N, Duschl WJ, Fraune S, et al. Metaorganisms in extreme environments: do microbes play a role in organismal adaptation? Zoology. 2018;127:1–19.

  3. 3.

    Drosophila 12 Genomes Consortium, Clark AG, Eisen MB, Smith DR, Bergman CM, Oliver B, et al. Evolution of genes and genomes on the Drosophila phylogeny. Nature. 2007;450:203.

  4. 4.

    Jaenike J, Unckless R, Cockburn SN, Boelio LM, Perlman SJ. Adaptation via symbiosis: recent spread of a Drosophila defensive symbiont. Science. 2010;329:212–5.

  5. 5.

    Douglas AE. Multiorganismal insects: diversity and function of resident microorganisms. Annu Rev Entomol. 2015;60:17–34.

  6. 6.

    Brune A. Symbiotic digestion of lignocellulose in termite guts. Nat Rev Microbiol. 2014;12:168–80.

  7. 7.

    Brune A, Dietrich C. The gut microbiota of termites: digesting the diversity in the light of ecology and evolution. Annu Rev Microbiol. 2015;69:145–66.

  8. 8.

    Cleveland LR. Symbiosis between termites and their intestinal protozoa. Proc Natl Acad Sci U S A. 1923;9:424–8.

  9. 9.

    Scharf ME, Karl ZJ, Sethi A, Boucias DG. Multiple levels of synergistic collaboration in termite lignocellulose digestion. PLoS One. 2011;6. https://doi.org/10.1371/journal.pone.0021709.

  10. 10.

    Peterson BF, Stewart HL, Scharf ME. Quantification of symbiotic contributions to lower termite lignocellulose digestion using antimicrobial treatments. Insect Biochem Mol Biol. 2015;59:80–8.

  11. 11.

    Peterson BF, Scharf ME. Metatranscriptome analysis reveals bacterial symbiont contributions to lower termite physiology and potential immune functions. BMC Genomics. 2016;17. https://doi.org/10.1186/s12864-016-3126-z.

  12. 12.

    Lilburn TG, Kim KS, Ostrom NE, Byzek KR, Leadbetter JR, Breznak JA. Nitrogen fixation by symbiotic and free-living spirochetes. Science. 2001;292:2495–8.

  13. 13.

    Desai MS, Brune A. Bacteroidales ectosymbionts of gut flagellates shape the nitrogen-fixing community in dry-wood termites. ISME J. 2012;6:1302–13.

  14. 14.

    Tiedje J. Ecology of denitrification and dissimilatory nitrate reduction to ammonium. In: Methods of soil analysis. Part 2. Chemical and microbiological properties. 1988. p. 179–244.

  15. 15.

    Potrikus CJ, Breznak JA. Gut bacteria recycle uric acid nitrogen in termites: a strategy for nutrient conservation. Proc Natl Acad Sci U S A. 1981;78:4601–5.

  16. 16.

    Hongoh Y, Sharma VK, Prakash T, Noda S, Toh H, Taylor TD, et al. Genome of an endosymbiont coupling N2 fixation to cellulolysis within protist cells in termite gut. Science. 2008;322 November:1108–9.

  17. 17.

    Kuwahara H, Yuki M, Izawa K, Ohkuma M, Hongoh Y. Genome of “Ca. Desulfovibrio trichonymphae”, an H2-oxidizing bacterium in a tripartite symbiotic system within a protist cell in the termite gut. ISME J. 2017;11:766–76.

  18. 18.

    Warnecke F, Luginbühl P, Ivanova N, Ghassemian M, Richardson TH, Stege JT, et al. Metagenomic and functional analysis of hindgut microbiota of a wood-feeding higher termite. Nature. 2007;450:560–5.

  19. 19.

    Ohkuma M, Noda S, Hattori S, Iida T, Yuki M, Starns D, et al. Acetogenesis from H2 plus CO2 and nitrogen fixation by an endosymbiotic spirochete of a termite-gut cellulolytic protist. Proc Natl Acad Sci U S A . 2015;112:10224–10230.

  20. 20.

    Hongoh Y, Sharma VK, Prakash T, Noda S, Taylor TD, Kudo T, et al. Complete genome of the uncultured termite group 1 bacteria in a single host protist cell. Proc Natl Acad Sci U S A . 2008;105:5555–5560.

  21. 21.

    Zheng H, Dietrich C, Radek R, Brune A. Endomicrobium proavitum, the first isolate of Endomicrobia class. Nov. (phylum Elusimicrobia) - an ultramicrobacterium with an unusual cell cycle that fixes nitrogen with a group IV nitrogenase. Environ Microbiol. 2016;18:191–204.

  22. 22.

    Hongoh Y, Sato T, Noda S, Ui S, Kudo T, Ohkuma M. Candidatus Symbiothrix dinenymphae: bristle-like Bacteroidales ectosymbionts of termite gut protists. Environ Microbiol . 2007;9:2631–2635.

  23. 23.

    Yuki M, Kuwahara H, Shintani M, Izawa K, Sato T, Starns D, et al. Dominant ectosymbiotic bacteria of cellulolytic protists in the termite gut also have the potential to digest lignocellulose. Environ Microbiol. 2015;17:4942–53.

  24. 24.

    Brune A. Endomicrobia: intracellular symbionts of termite gut flagellates. Endocytobiosis Cell Res. 2012;23:11–5.

  25. 25.

    Shimada K, Lo N, Kitade O, Wakui A, Maekawa K. Cellulolytic protist numbers rise and fall dramatically in termite queens and kings during colony foundation. Eukaryot Cell. 2013;12:545–50.

  26. 26.

    Noda S, Kitade O, Inoue T, Kawai M, Kanuka M, Hiroshima K, et al. Cospeciation in the triplex symbiosis of termite gut protists (Pseudotrichonympha spp.), their hosts, and their bacterial endosymbionts. Mol Ecol. 2007;16:1257–66.

  27. 27.

    Ohkuma M. Symbioses of flagellates and prokaryotes in the gut of lower termites. Trends Microbiol. 2008;16:345–52.

  28. 28.

    Ikeda-Ohtsubo W, Brune A. Cospeciation of termite gut flagellates and their bacterial endosymbionts: Trichonympha species and “Candidatus Endomicrobium trichonymphae.” Mol Ecol 2009;18:332–342.

  29. 29.

    Desai MS, Strassert JFH, Meuser K, Hertel H, Ikeda-Ohtsubo W, Radek R, et al. Strict cospeciation of devescovinid flagellates and Bacteroidales ectosymbionts in the gut of dry-wood termites (Kalotermitidae). Environ Microbiol. 2010;12:2120–32.

  30. 30.

    Radek R, Meuser K, Strassert JFH, Arslan O, Teßmer A, Šobotník J, et al. Exclusive gut flagellates of Serritermitidae suggest a major transfaunation event in lower termites: description of Heliconympha glossotermitis gen. Nov. spec. Nov. J Eukaryot Microbiol. 2018;65:77–92.

  31. 31.

    Yamin MA. Flagellates of the orders Trichomonadida Kirby, Oxymonadida Grassé, and Hypermastigida Grassi & Foà reported from lower termites (Isoptera families Mastotermitidae, Kalotermitidae, Hodotermitidae, Termopsidae, Rhinotermitidae, and Serritermitidae) and from the wood-feeding roach Cryptocercus (Dictyoptera: Cryptocercidae). Sociobiology. 1979;4:1–120.

  32. 32.

    Kitade O. Comparison of symbiotic flagellate faunae between termites and a wood-feeding cockroach of the genus Cryptocercus. Microbes Environ. 2004;19:215–20.

  33. 33.

    Waidele L, Korb J, Voolstra CR, Künzel S, Dedeine F, Staubach F. Differential ecological specificity of protist and bacterial microbiomes across a set of termite species. Front Microbiol. 2017;8. https://doi.org/10.3389/fmicb.2017.02518.

  34. 34.

    Bourguignon T, Lo N, Dietrich C, Šobotník J, Sidek S, Roisin Y, et al. Rampant host switching shaped the termite gut microbiome. Curr Biol. 2018;28:649–654.e2.

  35. 35.

    Boucias DG, Cai Y, Sun Y, Lietze V-U, Sen R, Raychoudhury R, et al. The hindgut lumen prokaryotic microbiota of the termite Reticulitermes flavipes and its responses to dietary lignocellulose composition. Mol Ecol . 2013;22:1836–1853.

  36. 36.

    He S, Ivanova N, Kirton E, Allgaier M, Bergin C, Scheffrahn RH, et al. Comparative metagenomic and metatranscriptomic analysis of hindgut paunch microbiota in wood- and dung-feeding higher termites. PLoS One. 2013;8:e61126.

  37. 37.

    Dietrich C, Köhler T, Brune A, Kohler T, Brune A. The cockroach origin of the termite gut microbiota: patterns in bacterial community structure reflect major evolutionary events. Appl Environ Microbiol. 2014;80:2261–9.

  38. 38.

    Mikaelyan A, Dietrich C, Köhler T, Poulsen M, Sillam-Dussès D, Brune A, et al. Diet is the primary determinant of bacterial community structure in the guts of higher termites. Mol Ecol . 2015;24:5284–5295.

  39. 39.

    Abdul Rahman N, Parks DH, Willner DL, Engelbrektson AL, Goffredi SK, Warnecke F, et al. A molecular survey of Australian and north American termite genera indicates that vertical inheritance is the primary force shaping termite gut microbiomes. Microbiome. 2015;3:5.

  40. 40.

    Rossmassler K, Dietrich C, Thompson C, Mikaelyan A, Nonoh JO, Scheffrahn RH, et al. Metagenomic analysis of the microbiota in the highly compartmented hindguts of six wood- or soil-feeding higher termites. Microbiome. 2015;3:56.

  41. 41.

    Tai V, James ER, Nalepa CA, Scheffrahn RH, Perlman SJ, Keeling PJ, et al. The role of host phylogeny varies in shaping microbial diversity in the hindguts of lower termites. Appl Environ Microbiol . 2015;81:1059–1070.

  42. 42.

    Zhang X. Leadbetter JR. Evidence for cascades of perturbation and adaptation in the metabolic genes of higher termite gut symbionts mBio. 2012;3:1–11.

  43. 43.

    Duarte S, Nobre T, Borges PAV, Nunes L. Symbiotic flagellate protists as cryptic drivers of adaptation and invasiveness of the subterranean termite Reticulitermes grassei Clément. Ecol Evol. 2018;8:5242–53.

  44. 44.

    Liu N, Li H, Chevrette MG, Zhang L, Cao L, Zhou H, et al. Functional metagenomics reveals abundant polysaccharide-degrading gene clusters and cellobiose utilization pathways within gut microbiota of a wood-feeding higher termite. ISME J. 2019;13:104.

  45. 45.

    Abe T. Evolution of life types in termites. In: Kawano S, Connell JH, and Hidaka T, editor. Evolution and coadaptation in biotic communities. University of Tokyo Press; 1987. p. 125–146.

  46. 46.

    Korb J. Termites Curr Biol. 2007;17:R995–9.

  47. 47.

    Lenz M. Food resources, colony growth and caste development in wood-feeding termites. In: Nalepa CA, Hunt J, editors. Nourishment and evolution in insect societies. Boulder, CO, USA: Westview Press; 1994. p. 159–209.

  48. 48.

    Korb J, Hoffmann K, Hartfelder K. Molting dynamics and juvenile hormone titer profiles in the nymphal stages of a lower termite, Cryptotermes secundus (Kalotermitidae)--signatures of developmental plasticity. J Insect Physiol. 2012;58:376–83.

  49. 49.

    Lainé LV, Lainé LV, Wright DJ. The life cycle of Reticulitermes spp (Isoptera: Rhinotermitidae): what do we know? Bull Entomol Res. 2003;93:267–378.

  50. 50.

    Douglas AE. Lessons from studying insect symbioses. Cell Host Microbe. 2011;10:359–67.

  51. 51.

    Ma D, Leulier F. The importance of being persistent: the first true resident gut symbiont in Drosophila. PLoS Biol. 2018;16:e2006945.

  52. 52.

    Staubach F, Baines JF, Kuenzel S, Bik EM, Petrov DA. Host species and environmental effects on bacterial communities associated with Drosophila in the laboratory and in the natural environment. PLoS One. 2013;8:e70749.

  53. 53.

    Hongoh Y, Deevong P, Inoue T, Moriya S, Trakulnaleamsai S, Ohkuma M, et al. Intra- and interspecific comparisons of bacterial diversity and community structure support coevolution of gut microbiota and termite host intra- and interspecific comparisons of bacterial diversity and community structure support coevolution of gut microbiota and termite host. Appl Environ Microbiol. 2005;71:6590–9.

  54. 54.

    Hongoh Y, Ekpornprasit L, Inoue T, Moriya S, Trakulnaleamsai S, Ohkuma M, et al. Intracolony variation of bacterial gut microbiota among castes and ages in the fungus-growing termite Macrotermes gilvus. Mol Ecol . 2006;15:505–516.

  55. 55.

    Benjamino J, Graf J. Characterization of the core and caste-specific microbiota in the termite, Reticulitermes flavipes Front Microbiol 2016;7:171.

  56. 56.

    Benjamino J, Lincoln S, Srivastava R, Graf J. Low-abundant bacteria drive compositional changes in the gut microbiota after dietary alteration. Microbiome. 2018;6. https://doi.org/10.1186/s40168-018-0469-5.

  57. 57.

    Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27:326–49.

  58. 58.

    Segata N, Izard J, Waldron L, Gevers D, Miropolsky L, Garrett WS, et al. Metagenomic biomarker discovery and explanation. Genome Biol. 2011;12:R60.

  59. 59.

    Whittle G, Shoemaker NB, Salyers AA. The role of Bacteroides conjugative transposons in the dissemination of antibiotic resistance genes. Cell Mol Life Sci. 2002;59:2044–54.

  60. 60.

    Liu N, Zhang L, Zhou H, Zhang M, Yan X, Wang Q, et al. Metagenomic insights into metabolic capacities of the gut microbiota in a fungus-cultivating termite (Odontotermes yunnanensis). PLoS One. 2013;8:e69184.

  61. 61.

    Lombard V, Golaconda Ramulu H, Drula E, Coutinho PM, Henrissat B. The carbohydrate-active enzymes database (CAZy) in 2013. Nucleic Acids Res. 2014;42(Database issue):D490–5.

  62. 62.

    Déjean G, Blanvillain-Baufumé S, Boulanger A, Darrasse A, de Bernonville TD, Girard A-L, et al. The xylan utilization system of the plant pathogen Xanthomonas campestris pv campestris controls epiphytic life and reveals common features with oligotrophic bacteria and animal gut symbionts. New Phytol. 2013;198:899–915.

  63. 63.

    Hottes AK, Meewan M, Yang D, Arana N, Romero P, McAdams HH, et al. Transcriptional profiling of Caulobacter crescentus during growth on complex and minimal media. J Bacteriol. 2004;186:1448–61.

  64. 64.

    Burnum KE, Callister SJ, Nicora CD, Purvine SO, Hugenholtz P, Warnecke F, et al. Proteome insights into the symbiotic relationship between a captive colony of Nasutitermes corniger and its hindgut microbiome. ISME J. 2011;5:161–4.

  65. 65.

    Baum DA, Futuyma DJ, Hoekstra HE. Lenski RE. Moore AJ: Peichel CL, et al. The Princeton guide to evolution. Princeton University Press; 2013.

  66. 66.

    McMahan EA. Feeding relationships and radio-isotope techniques. In: Biology of Termites. Academic Press, New York, London; 1969. p. 387–406.

  67. 67.

    Blanvillain S, Meyer D, Boulanger A, Lautier M, Guynet C, Denancé N, et al. Plant carbohydrate scavenging through TonB-dependent receptors: a feature shared by phytopathogenic and aquatic bacteria. PLoS One. 2007;2:e224.

  68. 68.

    Schauer K, Rodionov DA, de Reuse H. New substrates for TonB-dependent transport: do we only see the ‘tip of the iceberg’? Trends Biochem Sci. 2008;33:330–8.

  69. 69.

    Koebnik R. TonB-dependent trans-envelope signalling: the exception or the rule? Trends Microbiol. 2005;13:343–7.

  70. 70.

    Martens EC, Lowe EC, Chiang H, Pudlo NA, Wu M, McNulty NP, et al. Recognition and degradation of plant cell wall polysaccharides by two human gut symbionts. PLoS Biol. 2011;9:e1001221.

  71. 71.

    Fengel D, Grosser D. Chemische Zusammensetzung von Nadel- und Laubhölzern. Holz als Roh-und Werkstoff. 1975;33:32–4.

  72. 72.

    Pettersen RC. The chemical composition of wood. Adv Chem Ser. 1984;207.

  73. 73.

    Korb J, Lenz M. Reproductive decision-making in the termite, Cryptotermes secundus (Kalotermitidae), under variable food conditions. Behav Ecol. 2004;15:390–5.

  74. 74.

    Scheffrahn R. Cuban subterranean termite - Reticulitermes simplex (Hagen). 2007. http://entnemdept.uft.edu/creatures/urban/termites/p_simplex.htm. .

  75. 75.

    Dedeine F, Dupont S, Guyot S, Matsuura K, Wang C, Habibpour B, et al. Historical biogeography of Reticulitermes termites (Isoptera: Rhinotermitidae) inferred from analyses of mitochondrial and nuclear loci. Mol Phylogenet Evol. 2016;94 Pt B:778–90.

  76. 76.

    Smythe RV, Carter FL. Feeding responses to sound wood by the eastern subterranean termite, Reticulitermes flavipes An Entomol Soc Am 1969;62:335–337.

  77. 77.

    Cárdenas AM, Gallardo P, Toledo D. Suitability of multiple Mediterranean oak species as a food resource for Reticulitermes grassei Clément (Isoptera: Rhinotermitidae). Bull Entomol Res. 2018;108:532–9.

  78. 78.

    Janzow MP, Judd TM. The termite Reticulitermes flavipes (Rhinotermitidae: Isoptera) can acquire micronutrients from soil. Environ Entomol. 2015;44:814–20.

  79. 79.

    Judd TM, Landes JR, Ohara H, Riley AW. A geometric analysis of the regulation of inorganic nutrient intake by the subterranean termite Reticulitermes flavipes Kollar. Insects. 2017;8:97.

  80. 80.

    Wertz JT, Breznak JA. Stenoxybacter acetivorans gen. Nov., sp. nov., an acetate-oxidizing obligate microaerophile among diverse O2-consuming bacteria from termite guts. Appl Environ Microbiol. 2007;73:6819–28.

  81. 81.

    Kuhnigk T, Branke J, Krekeler D, Cypionka H, König H. A feasible role of sulfate-reducing bacteria in the termite gut. Syst Appl Microbiol. 1996;19:139–49.

  82. 82.

    Breznak JA. Biochemical aspects of symbiosis between termites and their intestinal microbiota. In: Anderson J, editor. Invertebrate - microbial interaction. Cambridge University Press: Cambridge; 1984. p. 173–203.

  83. 83.

    Li H, Yelle DJ, Li C, Yang M, Ke J, Zhang R, et al. Lignocellulose pretreatment in a fungus-cultivating termite. Proc Natl Acad Sci U S A. 2017;114:4709–14.

  84. 84.

    Huang X-F, Bakker MG, Judd TM, Reardon KF, Vivanco JM. Variations in diversity and richness of gut bacterial communities of termites (Reticulitermes flavipes) fed with grassy and woody plant substrates. Microb Ecol. 2013;65:531–6.

  85. 85.

    Harrison MC, Jongepier E, Robertson HM, Arning N, Bitard-Feildel T, Chao H, et al. Hemimetabolous genomes reveal molecular basis of termite eusociality. Nat Ecol Evol. 2018;2:557.

  86. 86.

    Misof B, Liu S, Meusemann K, Peters RS, Donath A, Mayer C, et al. Phylogenomics resolves the timing and pattern of insect evolution. Science. 2014;346:763–7.

  87. 87.

    Bushnell B. BBMap. sourceforge.net/projects/bbmap/.

  88. 88.

    Huson DH, Auch AF, Qi J, Schuster SC. MEGAN analysis of metagenomic data. Genome Res. 2007;17:377–86.

  89. 89.

    Suzuki R, Shimodaira H. Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics. 2006;22:1540–2.

  90. 90.

    Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, Mcglinn D, et al. vegan: community ecology package. Ordination methods, diversity analysis and other functions for community and vegetation ecologists. Version 2.4–3. URL https://CRAN.R-project.org/package=vegan. 2017.

  91. 91.

    Asnicar F, Weingart G, Tickle TL, Huttenhower C, Segata N. Compact graphical representation of phylogenetic data and metadata with GraPhlAn. PeerJ. 2015;3:e1029.

  92. 92.

    Buchfink B, Xie C, Huson DH. Fast and sensitive protein alignment using DIAMOND. Nat Methods. 2015;12:59–60.

  93. 93.

    Shallom D, Shoham Y. Microbial hemicellulases. Curr Opin Microbiol. 2003;6:219–28.

Download references

Acknowledgments

We thank Jan Sobotnik for kindly providing P. simplex colonies, as well as Charles Darwin University (Australia), especially S. Garnett and the Horticulture and Aquaculture team for providing logistic support to collect C. secundus. The Parks and Wildlife Comission, Northern Territory, Department of the Environment, Water, Heritage and the Arts gave permission to collect (Permit number 59044) and export (Permit PWS2016-AU-001559) the termites. This work was funded by DFG grants STA1154/2-1 and KO1895/16-1. We thank Craig Michell and the KAUST Bioscience Core Lab for sequence library generation and sequencing. This study was supported by King Abdullah University of Science and Technology (KAUST) and the High Performance and Cloud Computing Group at the Zentrum fuer Datenverarbeitung of the University of Tuebingen, the state of Baden-Wuerttemberg through bwHPC and the German Research Foundation (DFG) through grant no INST 37/935-1 FUGG. We thank Karen Meusemann and the 1KITE consortium, in particular the 1KITE Blattodea group, Alexander Donath, Lars Podsiadlowski, Bernhard Misof, Xin Zhou for granting access of the transcriptome data of R. santonsensis syn. R. flavipes. We thank the anonymous reviewers for helpful comments on the manuscript.

Funding

This work was funded by DFG grants STA1154/2–1 and KO1895/16–1 and supported by the King Abdullah University of Science and Technology (KAUST) and the High Performance and Cloud Computing Group at the Zentrum fuer Datenverarbeitung of the University of Tuebingen, the state of Baden-Wuerttemberg through bwHPC and the German Research Foundation (DFG) through grant no INST 37/935–1 FUGG.

Author information

FS, JK, LW designed the experiment. JK, FD provided study organisms. LW performed the experiments, CRV generated the sequence data, LW and FS analyzed the data. LW, FS, JK, CRV, and FD wrote the manuscript. All authors read and approved the final manuscript.

Correspondence to Fabian Staubach.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Figure S2. Evolutionary genealogy of genes - non-supervised orthologous groups (eggNOG) annotation structure. Figure S3. Metagenome sequence coverage of functional categories. Figure S4. Cluster dendrograms of the functional profiles of the protist community. Figure S6. Venn diagram of COGs in the pathways of the eggNOG category “information storage and processing” in the protist metagenoms. Figure S7. Cluster dendrograms of the functional profiles of the bacterial community. Figure S8. Frequency of bacterial genera. Figure S11. Venn diagram of COGs in the pathways of the eggNOG category “metabolism” in the bacterial metagenoms. Supplement section S13. Supplementary Material and Methods. Figure S14. Analysis workflow. Supplement section S15. Metagenomic Shotgun Data Analysis Workflow.

Additional file 2: Table S1. Sample Overview. Number of sequences before and after quality filtering, as well as taxonomic and functional annotation statistics

Additional file 3: Table S2. Complete list of results of Linear Discriminant Analysis (LDA) of protist dataset.

Additional file 4: Table S3. Transposase abundance in the bacterial dataset in the pathway “replication, recombination and repair”.

Additional file 5: Table S4. Complete list of results of Linear Discriminant Analysis (LDA) of bacterial dataset.

Additional file 6: Table S5. Complete list of bacterial COGs found within the pathway ‘carbohydrate transport and metabolism’.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Waidele, L., Korb, J., Voolstra, C.R. et al. Ecological specificity of the metagenome in a set of lower termite species supports contribution of the microbiome to adaptation of the host. anim microbiome 1, 13 (2019) doi:10.1186/s42523-019-0014-2

Download citation

Keywords

  • Lower termites
  • Metagenomics
  • Ecology
  • Adaptation