Diet may be a significant determinant of insect gut microbiome composition. However, the extent to which dietary shifts shape both the composition and relevant functions of insect gut microbiomes, and ultimately impact host energy balance (i.e. metabolic phenotype), is not well understood. We investigated the impacts of diet switching on Diploptera punctata females maintained on a dog food (DF) diet relative to those fed a comparatively sub-optimal cellulose-amended dog food (CADF) diet for 4 weeks. After this period, dietary shift resulted in a significantly higher average mass-specific standard metabolic rate (SMR) in CADF-fed females compared with DF-fed females. We also uncovered significant 13C-enrichment in DF-fed insect samples relative to CADF-fed insect samples and lowered bacterial essential amino acid (EAA) provisioning in CADF-fed samples. Differences in SMR and EAA provisioning were not accompanied by significant differences in overall microbiome composition between the two groups. However, cellulolytic and nitrogen-fixing bacterial families dominant in wild omnivorous cockroaches and wood-feeding termites were significantly enriched in CADF-fed females than in DF-fed females, at the end of the study. We propose that these changes in microbiome composition after dietary shifts are associated with changes in EAA provisioning and possibly SMR. Further studies are needed to comprehensively understand the relative importance of gut microbial functions among the complexity of factors known to underscore SMR responses in insects under varying dietary conditions.

Animals strive to maintain a balance between energy inputs (sources) and outputs (sinks). Foraging, ingestion, digestion and nutrient absorption represent expenditures to acquire energy, whereas growth, maintenance of homeostasis, performance and reproduction represent expenditures dependent on acquired energy (Karasov, 1986). Several factors impact energy acquisition, but dietary quality, defined as the ratios of nutritional (protein:carbohydrate:lipid ratios) and non-nutritional (allelochemicals) components of ingested food is one of the most important (Behmer, 2008; Roeder and Behmer, 2014).

Dietary quality directly influences the minimum energy expenditure of ectotherms at a given ambient temperature under resting conditions, i.e. the standard metabolic rate (SMR) (Chown and Gaston, 1999). However, the impacts of dietary quality on SMR of ectotherms are not always straightforward (Cruz-Neto and Bozinovic, 2004; Burton et al., 2011). For example, in the harvestman, Pachylus paessleri (Opiliones), Naya et al. (2007) found significantly lower SMR and fecundity in response to feeding on a low-quality carbohydrate-rich diet relative to a high-quality protein-rich diet. In contrast, feeding on low-quality diets increased SMR and reduced growth in fifth instar Spodoptera eridania larvae (Karowe and Martin, 1989), and increased SMR in the locust Locusta migratoria (Zanotto et al., 1997) and the American cockroach Periplaneta americana (Ayayee et al., 2018). However, feeding on preys varying in lipid:protein ratios did not seem to impact SMR in the predatory wolf spider Pardosa prativa (Jensen et al., 2010), and in some cases SMR is constrained by morphological traits (e.g. size), developmental trajectories and physiological variation (Clark et al., 2016).

Indirectly, dietary quality may impact ectotherm energy balance through the disruption of host gut-associated microbiomes. For example, diet-driven differences in gut microbiomes has been demonstrated across various insect orders (Colman et al., 2012; Yun et al., 2014), as well as various wood-feeding beetle species, such as the generalist Asian long-horned beetle Anoplophora glabripennis (Geib et al., 2009; Scully et al., 2014), and specialists Monochamus alternatus (pine tree-feeding) and Psacothea hilaris (mulberry tree-feeding) beetles (Kim et al., 2017), various lignocellulolytic-fed omnivorous cockroach species (Bertino-Grimaldi et al., 2013; Pérez-Cobas et al., 2015), and various wood-feeding, soil-feeding and humus-feeding termite species (Raychoudhury et al., 2013; Rossmassler et al., 2015; Mikaelyan et al., 2016; Su et al., 2016). Diet-induced differences in gut microbiome compositions can also subsequently influence host gut microbiome functions. For example, differences in essential amino acid provisioning functions by bacterial and fungal members of host gut microbiomes (a function relevant for optimal insect host peptide and protein biosynthesis) under different dietary conditions have been verified in insects, such as A.glabripennis fed artificial diet versus host tree diets (Ayayee et al., 2016a), P.americana fed an optimal dog food diet or a sub-optimal cellulose-amended dog food diet (Ayayee et al., 2016b), and various soil mesofauna (Larsen et al., 2011; Pollierer et al., 2019), Enchytraeids (Larsen et al., 2016a), as well as earthworms (Larsen et al., 2016b; Potapov et al., 2019) feeding on various environmental dietary sources. These diet-induced disruptions of gut microbiome function are attributable to changes in gut microbiome functions, which may result in altered host metabolic phenotypes, and consequently changes in host SMR.

Although the impacts of diet type and quality on SMR are well studied, our understanding of the role of gut microbiota in mediating SMR responses in insects is comparatively nascent. An exception is a recent study that investigated the relationship between diet shifts and gut microbiome composition and metabolite production by host and microbiome across different genetic lines of Drosophila melanogaster (Jehrke et al., 2018). Diet had no significant effect on gut microbiome compositions or on investigated metabolite–microbiome associations among genetic lines or developmental stages. This observation was attributed to the relatively high quality of the diets used and the comparatively simple composition of the gut microbiome of D. melanogaster (Wong et al., 2011; Jehrke et al., 2018), which may have masked any substantial dietary impacts. To date, we still do not quite understand how shifts from optimal laboratory-based diets to near-sub-optimal diets impact gut microbiome composition and functions, and how these interactions may potentially mediate the manifestation of the host metabolic phenotype.

We hypothesize that the energy expenditure of hosts is mediated by differences in gut microbiome composition and function. To evaluate our hypothesis, we used the only known viviparous cockroach species Diploptera punctata. This cockroach was chosen because of its viviparous reproductive strategy and the scarcity of data detailing impacts of dietary quality on this strategy in sexually mature females prior to starting reproductive events. Post-mated D. punctata females carry developing embryos in a specialized structure called a brood sac while providing nutritive ‘milk’ secretions (Williford et al., 2004; Bell et al., 2007). Like most cockroach species, D. punctata harbors the bacterial endosymbiont Blattabacterium spp. in specialized cells in the fat body (Gier, 1936), and also has a diverse adult gut microbiome (Ayayee et al., 2017; Dietrich et al., 2014). We exposed cockroaches to two diets differing in quality over 4 weeks, and measured changes in SMR (a proxy for energy expenditure), body mass (a proxy for growth), the composition of the gut microbiome, and provisioning of essential amino acids, a crucial function of the microbiome. We predicted changes in the composition and function of the gut microbiomes of D. punctata insects fed a comparatively sub-optimal diet relative to those fed an optimal diet, underscored by changes in SMR.

Insect rearing and study layout

Diploptera punctata (Eschscholtz 1822) females were obtained from the insectary in the Department of Entomology at The Ohio State University. Adult laboratory-reared females used in this study were obtained from an initial colony fed an optimal dog food (DF) diet (Red Flannel Hi-Protein Formula dog food, PMI Nutrition, St Louis, MO, USA). Cockroaches were maintained in closed containers and provided with ventilation and water ad libitum under ambient laboratory conditions (ambient temperature 22±2°C, relative humidity 45%). The sub-optimal cellulose-amended dog food diet (CADF) was made by diluting DF with cellulose (Sigma-Aldrich, St Louis, MO, USA) at a 30:70 ratio (w/w), as used in a previous study (Ayayee et al., 2018). Cellulose was mixed with the DF diet because it is not a phago-deterrent for cockroaches (Bignell, 1978), but promotes increased food consumption for nutrient acquisition from the DF fraction in the CADF diet (A. Muñoz-Garcia, personal observations). It is also relatively recalcitrant to digestion by insects without gut microbial assistance (Kane and Breznak, 1991) and is shown to impact gut microbial compositions (Bignell, 1977, 1978; Gijzen et al., 1994; Bell et al., 2007) and microbial functions (Bignell, 1978; Ayayee et al., 2016b).

To investigate the impacts of diet on SMR and gut microbiome composition, we used 22 adult unmated females of varying ages from the DF-fed laboratory colony, to avoid any confounding effects due to sex. Selected females were first kept individually in containers without food for 24 h, to ensure that they would become post-absorptive before the start of the study on day 1. On day 1, we measured SMR and body mass, after which females were randomly assigned to a DF (N=11) or CADF (N=11) group for the duration of the study. On day 28 of the experiment, we measured body mass and SMR again in the same individuals.

A separate group of females was used to study essential amino acid (EAA) provisioning by the gut microbiota in D. punctata. For this purpose, we used previously collected samples from mated females (head capsules and six legs) and their associated embryos, unmated females and orphaned nymphs obtained from dead females from colonies that were part of another long-term study and had been feeding on the same DF and CADF diets used in this study for 6–8 weeks. Samples from both DF and CADF dietary groups were combined into DF-insect and CADF-insect sample groups, respectively, without any distinction between adult females and embryo samples for simplicity. Samples were lyophilized for 48 h, pulverized after lyophilization, and stored at −80°C, before compound-specific stable isotope analysis of EAA.

Measurement of standard metabolic rate

On day 1 and day 28, we measured oxygen consumption using the manual bolus integration method (Lighton, 2008). Briefly, females were incubated in a 60 ml airtight syringe containing dry air for 90 min at an ambient temperature of 30±2°C. Reference air samples (from an empty airtight syringe with dry air) were also collected using the same collection method. Subsequently, 40 ml of air was drawn from the airtight syringe and injected into an Oxzilla oxygen analyser (Sable Systems International, Las Vegas, NV, USA) at a flow rate of 260 ml min−1, after passing through a column of silica gel and ascarite to remove water vapor and CO2. The amount of oxygen consumed by each female was then calculated by determining the difference in the amount of oxygen (VO2) left in the airtight syringe used for sampling insects relative to the empty reference syringe (Withers, 1977). Oxygen consumed was subsequently converted to energy expenditure using the factor 20.08 J ml−1 O2 (Schmidt-Nielsen, 1995). For each female, body mass was measured before and after incubation. We calculated mass-specific SMR as the quotient of energy expenditure and body mass.

The temperature in the incubation chamber (30°C) was constant during the trials, but different from that of the laboratory (∼22°C). We used this set-up to be able to compare values for SMR across individuals and treatments, avoiding the fluctuations in ambient room temperature in the laboratory, which is a common practice (Dingha, 2009; DeVries et al., 2013; Harrison and Fewell, 2016). Moreover, the temperature difference between the insect colony and the incubation temperature (5–7°C) is within the ∼10°C range within which mass-specific SMR is not significantly affected by temperature changes in cockroaches (Dingha, 2009; Harrison and Fewell, 2016). Diploptera punctata females were incubated for 90 min to compensate for their smaller sizes relative to other cockroaches (Ayayee et al., 2018), and to ensure significant oxygen consumption given the discontinuous gas exchange cycle of cockroaches (Schimpf et al., 2012a,b).

DNA extraction and gut bacterial load determination

After SMR measurements on day 28, all females (N=22) were immobilized at 4°C for 5 min and surface-sterilized as described by Ayayee et al. (2017); the entire gut was removed and gut mass was determined. We extracted DNA from each insect sample using the Qiagen DNeasy Blood and Tissue DNA extraction kit (Qiagen, Valencia, CA, USA). We carried out quantification of bacterial loads out using the bacterial primers 357f and 519r (Turner et al., 1999) via qPCR using the IQ SYBR Green Supermix (Bio-Rad Laboratories, Hercules, CA, USA), cockroach gut DNA (∼16–18 ng µl−1) and water to a final 20 µl volume. Bacterial 16S rRNA plasmid standards were run in triplicate and insect DNA samples in duplicate. Reaction conditions consisted of an initial denaturation at 95°C for 10 min, followed by 40 cycles of 95°C for 15 s, 55.3°C for 15 s and 68°C for 20 s, and a final elongation step of 62°C for 20 s, using the Realplex Mastercycler (Eppendorf, Westbury, NY, USA). Bacterial load of template samples was calculated based on the generated standard curve of plasmid standards, after examination of melting curves and divided by gut mass as copy number/gram.

Amplicon generation, processing and analyses

Amplicon library preparation of the V4 region of the bacterial 16S rRNA gene was carried out using barcoded modified forward (515f) (Parada et al., 2016) and reverse (806r) (Apprill et al., 2015) primers for multiplexing and enhanced amplification of the bacterial and archaeal groups in environmental samples (Walters et al., 2016). DNA samples from individual cockroaches from each dietary treatment (N=10 from DF-fed individuals, one dropped due to low quality; N=12 from CADF-fed individuals, another individual was added from a similar CADF colony) were used for amplicon preparation. Briefly, we carried out triplicate PCR reactions out for all samples (with two technical DNA replicates each) using the Phusion polymerase and the Phusion High-Fidelity PCR Master Mix with HF buffer, water and barcoded-primers in a final reaction volume of 20 µl. Reaction conditions consisted of an initial denaturation for 30 s at 98°C, 30 cycles of 10 s at 98°C (denaturation), 10 s at 65°C (annealing) and 8 s at 72°C (extension), with a final extension phase of 5 min at 72°C. Gel-verified amplified products were pooled, purified using the Axygen AxyPrep Mag PCR Clean-Up kit (Axygen Biosciences, Union City, CA, USA) following the manufacturer's protocol, and concentration was determined using a Qubit 3.0 fluorometer (Invitrogen, Carlsbad, CA, USA). Equimolar amounts of purified PCR products (12 ng) were then pooled together (yielding a total of 44 samples), and the multiplexed sample was sequenced on the Illumina Miseq platform using V2 chemistry (2×250 PE) at the University of Minnesota Genomics Center (UMGC, St Paul, MN, USA). Raw sequence data for all 44 samples are available in the NCBI Sequence Read Archive under the following BioProject accession number: PRJNA598044.

All 44 de-multiplexed reads were first merged using the fastqmerge_pairs command in usearch (Edgar, 2010), with default quality filtering parameters, and minimum and maximum merged read lengths of 250 and 253, respectively. Primary quality checking (filtering and error learning and de-replication) of merged fastq reads was performed using DADA2 (Callahan et al., 2016) and removal of primer sequences using cutadapt (Martin, 2011). Subsequent chimera removal (removeChimeraDenovo) and taxonomy assignments of amplicon sequence variants (ASV) using the SILVA database (version Silva_v132) (Quast et al., 2012) were carried out using DADA2. The final ASV table was then formatted to include taxonomic information and exported to a biom format for downstream analysis in QIIME (version 1.9) (Caporaso et al., 2010). Technical replicates of each sample were subsequently combined (avoiding over-inflation in subsequent downstream statistical analyses), resulting in the original 22 insect samples.

Alpha (count, richness and evenness measures) (Simpson, 1949; Shannon, 1957; Chao, 1984) and beta (Bray–Curtis distance matrix) (Bray and Curtis, 1957) diversity measures were analysed on the ASV table at the ASV and family level to assess within- and between-sample diversity among dietary categories using non-parametric analyses. Dietary category comparisons were examined using MRPP (Mielke, 1984), ANOSIM (D'Argenio et al., 2014) and PERMANOVA (Anderson, 2017). An indicator species analysis at the family level, examining microbial community members most likely driving differences among cockroach gut microbiomes from both dietary groups, was done using the group_significance command in QIIME (version 1.9) followed by a Kruskal–Wallis test. Figures were generated using JMP Pro 14 (SAS, Cary, NC, USA).

Assessment of microbial EEA provisioning

Pulverized samples were analysed at the Stable Isotope Facility at the University of California, Davis (Davis, CA, USA). Briefly, samples were acid-hydrolysed in 6 mol l−1 HCl at 150°C for 70 min under an N2 headspace and subsequently derivatized via methoxycarbonylation-esterification (MOC) (Yarnes and Herszage, 2017) as modified from Walsh et al. (2014). Two technical replicates per D. punctata and dietary samples were analysed. Exogenous carbon addition, kinetic isotope effects from derivatization reagents, and normalization to the international reference for δ13C, V-PDB, was conducted following the procedures outlined in Walsh et al. (2014). EAA δ13C data were obtained for isoleucine, leucine, lysine, phenylalanine, methionine, threonine and valine from all samples. The mean precision of the technical replicates across EAAs was ±0.55‰, which is well below the established quality control value of ±1.25. Final accuracy, as determined by the mean absolute difference in the measured and known δ13C values of EAAs from a quality assessment mixture of amino acids, was within ±0.49‰.

EAA provisioning by gut microbiota in D. punctata was assessed using the established stable isotope fingerprinting technique (Larsen et al., 2009) and subsequently used to assess EAA provisioning in a variety of macroinvertebrates (Ayayee et al., 2016a; Larsen et al., 2013, 2016b). Briefly, the δ13C values of EAAs in consumers (in this case D. punctata) biomass should be similar to those of EAAs from dietary materials (typically within ±2‰ range of mean dietary δ13C EAAs), as consumers rely exclusively on dietary sources for EAAs (McMahon et al., 2010; Newsome et al., 2011). Origin of consumer EAAs and the relative contribution of bacterial, fungal and plant EAAs to consumer EAAs were subsequently determined through a linear discriminant function analysis (LDA), using δ13C EAAs of the original representative bacterial, fungal and plant dataset from Larsen et al. (2013) as training data, calibrated to δ13C EAAs data of the same representative samples from the Stable Isotope Facility.

Statistical analyses

We first tested all measured variables for normality and examined this using the Shapiro–Wilk’s W goodness-of-fit test (Shapiro and Wilk, 1965) before analyses. Non-normal data were subsequently log-transformed, and normality of transformed data was tested again using the Shapiro–Wilk test. Insect body mass, gut mass, mass-specific SMR and bacterial copy numbers failed this initial normality screening and were all log-transformed. A mixed-model analysis (for repeated measures) with diet, time and their interaction as fixed effects, and samples as random factor, was carried out with SMR and body mass as our dependent variables, followed by Tukey's honestly significant difference pairwise mean comparisons. Gut bacterial load (16S rRNA copy numbers) and gut masses by day 28 were analysed using one-way ANOVA, with diet as the main factor. All these statistical analyses were performed using JMP (SAS).

Analyses of δ13C EAA values and 13C-isotopic enrichment among insect samples and dietary groups (DF and CADF) were, respectively, carried out using ANOVA with treatment groups (DF-insects, CADF-insects, DF-diet, CADF-diet) and amino acids (all seven EAAs) as factors, first without and then with normalization to respective dietary δ13C EAAs in JMP (SAS). We carried out an interlaboratory calibration for δ13C EAA data from this study and the representative fungal, bacterial and plant δ13C EAAs from Larsen et al. (2013), before analyses. LDA in JMP was used to create a classification model based on δ13C EAA values from a training data set consisting of three classifiers, fungi, bacteria and plants to predict biosynthetic origins of insect EAAs. The LDAs were run using R package MASS (R; http://www.R-project.org).

Effects of diet on body mass, gut mass, bacterial 16S rRNA copy number and mass-specific SMR

Diet (F=0.72, P=0.41), time (day 1 and day 28) (F=0.005, P=0.94) and their interaction (F=2.12, P=0.16) did not significantly impact body mass over the course of the study (Fig. 1A). There was a ∼3.5% body mass loss in CADF-fed females, and a ∼4.2% body mass gain in DF-fed females by day 28. Similarly, gut mass (F=0.13, P=0.72) and 16S rRNA gene copy number (F=0.10, P=0.76) were not significantly different between DF and CADF-fed females by day 28.

Fig. 1.

Effects of diet on body mass and standard metabolic rate. Effects of dietary shift on (A) body mass and (B) standard metabolic rate of dog food-fed (DF-fed, N=11) and cellulose-amended dog food-fed (CADF-fed, N=11) Diploptera punctata females from the mixed model analysis. Different letters indicate significant differences. There were no significant differences in mean body mass between DF-fed and CADF-fed females on day 1 (CADF-fed females, 0.25±0.015 g, mass±s.e.m.; DF-fed females, 0.23±0.007 g) or day 28 (CADF-fed females, 0.241±0.01 g; DF-fed females, 0.24±0.01 g). Connecting lines between time points do not imply any linear changes.

Fig. 1.

Effects of diet on body mass and standard metabolic rate. Effects of dietary shift on (A) body mass and (B) standard metabolic rate of dog food-fed (DF-fed, N=11) and cellulose-amended dog food-fed (CADF-fed, N=11) Diploptera punctata females from the mixed model analysis. Different letters indicate significant differences. There were no significant differences in mean body mass between DF-fed and CADF-fed females on day 1 (CADF-fed females, 0.25±0.015 g, mass±s.e.m.; DF-fed females, 0.23±0.007 g) or day 28 (CADF-fed females, 0.241±0.01 g; DF-fed females, 0.24±0.01 g). Connecting lines between time points do not imply any linear changes.

We found an overall significant diet by time interaction on mass-specific SMR (F=5.50, P=0.029). There was no significant difference in mass-specific SMR between dietary groups on day 1 (P=0.99; Fig. 1B). By day 28, the mass-specific SMR of the CADF-fed females, 2.06±0.35 (mW g−1, mean±s.e.m.), was approximately twice that of DF-fed females: 1.03±0.09 (P=0.023;Fig. 1B). Body mass, gut mass, bacterial 16S rRNA copy number and mass-specific SMR data are provided in Table S1.

Effect of diet on microbiome composition

Quality checking (merge, multistep filtering, error learning and dereplication, removal of chimeras) via DADA2 retained 83.95% (1,160,292) of initial merged 1,382,481 reads resulting in a mean of 52,559±22,106 across samples. Filtering yielded a total of 2348 ASVs, which resulted in 165 bacterial families across dietary groups. Both family-level and ASV tables were rarefied to 17,000 reads per sample prior to all diversity analyses. None of the four alpha diversity indices measured yielded significant differences between the two dietary feeding groups at the taxon (ASV) level or family level (Table 1). We uncovered no significant differences in microbiome composition between DF-fed and CADF-fed dietary groups at the taxon (ASV) or family level (Table 2).

Table 1.

Microbial alpha diversity estimates at the taxon level (amplicon sequence variant, ASV) for DF-fed and CADF-fed Diplopterapunctata females, following rarefaction to 17,000 reads per sample

Microbial alpha diversity estimates at the taxon level (amplicon sequence variant, ASV) for DF-fed and CADF-fed Diploptera punctata females, following rarefaction to 17,000 reads per sample
Microbial alpha diversity estimates at the taxon level (amplicon sequence variant, ASV) for DF-fed and CADF-fed Diploptera punctata females, following rarefaction to 17,000 reads per sample
Table 2.

Between-categories tests of beta diversity at the taxon (amplicon sequence variant, ASV) and family level for DF-fed and CADF-fed D. punctata females, following rarefaction to 17,000 reads per sample

Between-categories tests of beta diversity at the taxon (amplicon sequence variant, ASV) and family level for DF-fed and CADF-fed D. punctata females, following rarefaction to 17,000 reads per sample
Between-categories tests of beta diversity at the taxon (amplicon sequence variant, ASV) and family level for DF-fed and CADF-fed D. punctata females, following rarefaction to 17,000 reads per sample

We found differences in the relative in abundances of 10 bacterial families between the two dietary groups (Table 3). Of these, Xanthomonadaceae, Sphingobacteriaceae, Lactobacillaceae, Caulobacteraceae, Rhizobiaceae and Methylophilaceae (in order of decreasing abundance) were significantly more abundant in DF-fed females, whereas the remaining four families [Clostridiaceae_1, Paracaedibacteraceae, Unassigned JGI_0000069-P22 and Absconditabacteriales_(SR1)] were more abundant in CADF-fed females (Table 3).

Table 3.

Mean abundances (read count) and calculated relative abundances of the bacterial families that were significantly different between DF-fed and CADF-fed D. punctata females, and their associated Kruskal–Wallis test statistics

Mean abundances (read count) and calculated relative abundances of the bacterial families that were significantly different between DF-fed and CADF-fed D. punctata females, and their associated Kruskal–Wallis test statistics
Mean abundances (read count) and calculated relative abundances of the bacterial families that were significantly different between DF-fed and CADF-fed D. punctata females, and their associated Kruskal–Wallis test statistics

Gut microbial EAA provisioning

Despite the lack of significant differences in gut microbiome composition, we uncovered significant differences in gut microbial EAA provisioning between DF-fed females and CADF-fed females. Across all samples, threonine and lysine were more enriched than the remaining five EAAs, with leucine having the lowest raw δ13C values (F=106.1, P<0.0001) (Table 4). Mean δ13C of the seven EAAs in both CADF-fed and DF-fed insects were significantly more 13C-enriched relative to respective diets (F=48.43, P<0.0001) (Table 4). However, the mean δ13C EAA of DF-fed females was significantly higher than that of CADF-fed females (P<0.0001) (Table 4). Isotopic offset following normalization of insect δ13C EAAs to dietary δ13C EAAs (Δδ13CEAA13Cinsect samples EAA−δ13Crespective diet EAA), revealed significant differences between insect samples and corresponding diets across all seven EAAs (F=32.80, P<0.0001) (Table 4). However, mean 13C-isotopic offsets across all seven EAAs between DF-fed and CADF-fed females were not significantly different (Table 4). Mean 13C-isotopic offsets of individual EEAs were comparable between DF-fed and CADF-fed females, except for lysine, which was more 13C-enriched in DF insects (11.30±0.66‰) relative to CADF-fed insects (8.53±0.56‰; F=1.76, P=0.04; Fig. 2). The amino acid threonine had the lowest δ13C value in both diets and insect samples and experienced the least isotopic offset. The δ13C EAA data for samples used in this analysis are provided in Table S2.

Table 4.

δ13C signatures of seven essential amino acids across samples and their calculated 13C-discrimination values, for DF-fed and CADF-fed D. punctata females and their respective diets

δ13C signatures of seven essential amino acids across samples and their calculated 13C-discrimination values, for DF-fed and CADF-fed D. punctata females and their respective diets
δ13C signatures of seven essential amino acids across samples and their calculated 13C-discrimination values, for DF-fed and CADF-fed D. punctata females and their respective diets
Fig. 2.

Gut microbial essential amino acid provisioning. The 13C-offsets (Δ13CEAA: enrichment or depletion) of essential amino acids in DF-fed (N=5) and CADF-fed females (N=7) were normalized to respective dog food (DF) and cellulose-amended (CADF) diets. Only lysine (Lys) was significantly enriched in DF-fed insects (11.30±0.66) than in CADF-fed insects (8.53±0.56) (*P=0.04). Values represent mean values (±s.e.m.) for duplicate insect samples and dietary samples. Ile, isoleucine; Leu, leucine; Lys, lysine; Met, methionine; Phe, phenylalanine; Thr, threonine; Val, valine.

Fig. 2.

Gut microbial essential amino acid provisioning. The 13C-offsets (Δ13CEAA: enrichment or depletion) of essential amino acids in DF-fed (N=5) and CADF-fed females (N=7) were normalized to respective dog food (DF) and cellulose-amended (CADF) diets. Only lysine (Lys) was significantly enriched in DF-fed insects (11.30±0.66) than in CADF-fed insects (8.53±0.56) (*P=0.04). Values represent mean values (±s.e.m.) for duplicate insect samples and dietary samples. Ile, isoleucine; Leu, leucine; Lys, lysine; Met, methionine; Phe, phenylalanine; Thr, threonine; Val, valine.

Given the significant 13C-isotopic offsets in EAAs between insect and dietary samples, we investigated the sources of isoleucine, leucine, lysine, phenylalanine and valine in the insect hosts using LDA. The model was validated by the successful separation of bacteria (N=12), fungi (N=9) and plant (N=12) into distinct groups (F=55, P<0.0001; Wilk's lambda=0.007, a test of appropriateness of separation of classifiers and for group membership prediction of non-classifier samples; Table S3). The validated predictive model was then used to predict possible sources of EAAs in insect and dietary samples, which was based on the location of insect and dietary samples in the LDA plot, and the proximity to classifier centroids (fungi, bacteria and plants). The resulting predictions of samples were subsequently validated based on the correct assignment of the test fungus, Fusarium solanum, to the fungal classifier group (Fig. 3). Plants were not a major source of EAAs in DF-fed and CADF-fed cockroaches (Fig. 4). The DF and CADF diets were closer to the fungal classifier than to the bacterial classifier group. This may be attributed to the high fungal protein component of the DF diet (brewer's yeast accounts for ∼26% of crude protein fraction) and the CADF diet (∼7.8% brewer's yeast in the CADF diet).

Fig. 3.

Predicted possible sources of essential amino acids in insect and dietary samples. A linear discriminant function analysis (LDA) plot based on δ13CEAA from DF-fed (N=5) and CADF-fed insects (N=7), their respective DF (N=1) and CADF (N=2) diets, three classifier groups [fungi (N=9), bacteria (N=12) and plants (N=12)], and the test fungus Fusarium solani (N=2). Variation of classification accounted for 64% along LD1 and 36% along LD2. The shaded ellipses signify the 95% confidence limits for each classifier group, and the dashed lines between classifiers represent the decision boundaries for each classifier group. The essential amino acids used were isoleucine, leucine, lysine, phenylalanine, and valine.

Fig. 3.

Predicted possible sources of essential amino acids in insect and dietary samples. A linear discriminant function analysis (LDA) plot based on δ13CEAA from DF-fed (N=5) and CADF-fed insects (N=7), their respective DF (N=1) and CADF (N=2) diets, three classifier groups [fungi (N=9), bacteria (N=12) and plants (N=12)], and the test fungus Fusarium solani (N=2). Variation of classification accounted for 64% along LD1 and 36% along LD2. The shaded ellipses signify the 95% confidence limits for each classifier group, and the dashed lines between classifiers represent the decision boundaries for each classifier group. The essential amino acids used were isoleucine, leucine, lysine, phenylalanine, and valine.

Fig. 4.

Relative contributions of bacterial and fungal essential amino acids to overall insect and dietary essential amino acids. Relative input of bacterial, fungal and plant essential amino acids (EAAs) in dog food-fed (DF, N=5) and cellulose-amended dog food-fed (CADF, N=7) D.punctata females, relative to respective diets, based on LDA predicted probabilities.

Fig. 4.

Relative contributions of bacterial and fungal essential amino acids to overall insect and dietary essential amino acids. Relative input of bacterial, fungal and plant essential amino acids (EAAs) in dog food-fed (DF, N=5) and cellulose-amended dog food-fed (CADF, N=7) D.punctata females, relative to respective diets, based on LDA predicted probabilities.

We calculated the relative contributions of bacterial and fungal EAAs to overall insect and dietary EAAs using the posterior probability scores from the validated predictive model. For DF and CADF diets, EAAs were primarily fungal in origin, with no bacterial input (Fig. 4). In contrast, both DF-fed and CADF-fed insects had bacterial and fungal EAA inputs, with a higher bacterial input in DF-insects relative to CADF-insects (Fig. 4). Bacterial EAAs may have originated from gut microbiota, whereas fungal EAAs may have originated from the diet or from gut-associated fungi.

We found that shifting D. punctata from an optimal (DF) to a sub-optimal (CADF) diet resulted in a significant increase in mass-specific SMR of the host. We also found that dietary shifts altered microbial EAA provisioning, a function perhaps associated with changes in relative abundances of specific bacterial families. We propose that these changes in the microbiome are concomitant with the altered metabolic phenotype of the hosts. Given that the underlying factors impacting mass-specific SMR in response to feeding on sub-optimal diets are not straightforward, potential explanations underscoring the results from this study may be a combination of host-mediated and microbe-mediated factors.

From a host-mediated point of view, the higher SMR may be attributed to higher physiological activity in the gastrointestinal tract; however, we did not observe significant differences in gut mass between the two dietary groups. Some insects may reduce food consumption on sub-optimal diets (Lee, 2007), resorting to post-absorptive processes, such as oxidation of stored reserves (lipids and proteins) (McCue et al., 2015; Zanotto et al., 1993), leading to increased SMR and relatively reduced body masses. In this study, CADF-fed females might have a higher mobilization of lipid reserves in the fat body, potentially explaining the relative reduction of body mass in the CADF-fed female group by day 28. The decrease in SMR and the lack of significant difference in body mass in response to the same sub-optimal diet (i.e. CADF diet) has been previously reported in P. americana cockroaches (Ayayee et al., 2018); in that experiment, animals were fed with CADF for only 7 days.

From a microbiome-mediated point of view, the higher SMR in CADF-fed females may be attributed to diet-induced changes in microbiome function, which might lead to increases in host energy expenditure. Some changes observed in response to feeding on high cellulose diets include reductions in microbe-derived metabolites for energy generation (acetate and lactate) and reduced provisioning of essential amino acids. For example, P. americana cockroaches fed a dog food diet had significantly higher amounts of microbe-derived acetate and lactate and higher lactic acid bacterial counts, relative to individuals fed a diet with high cellulose content (Kane and Breznak, 1991). These microbe-derived products are estimated to provide ∼14% of the respiratory requirements of cockroaches and termites (Breznak and Kane, 1990; Kane and Breznak, 1991). Although we did not quantify acetate or lactate in this study, we observed changes in 13C-enrichment and EAA provisioning that support the proposed microbe-mediated SMR response in CADF-fed females. This is corroborated by the lower bacterial EAA input in CADF-fed females relative to DF-fed females (Fig. 4), despite the relatively higher fungal protein component of the DF diet compared with the CADF diet. The higher fungal EAA input in CADF-fed females may be the result of a combination of CADF diet consumption and diet-induced changes in gut fungal community membership. For example, Kakumanu et al. (2018) determined that laboratory-reared German cockroaches, Blattella germanica, maintained on rodent chow (with low fiber content) had significantly higher fungal abundances, and lower ciliate and protist abundances in their guts, compared with field-collected (urban apartment) omnivorous individuals. Although we did not study the gut fungal composition in this study, and insect gut fungal–bacteria interactions are unclear, the higher fungal EAA input in CADF-fed females may have derived from the incorporation of fungal EAAs from the gut following increased gut fungal abundances and cellulolytic activity and subsequent EAA biosynthesis. Ultimately, the microbial function investigated in this study cannot be directly attributed to any bacterial or fungal taxon, regardless of the dietary group. This is because EAA biosynthesis is a feature shared by all free-living and endosymbiotic bacteria, as well as fungi. Coupling broad microbial functions, such as EAA provisioning with taxon-specific or physiology-specific microbial functions (e.g. nitrogen fixation or cellulolytic activity) using stable isotope techniques, may provide a possible approach with which to both uncover and assign microbial functions to specific gut microbiome members in similar studies going forward. Overall, results from this study indicate the potential functions of the gut microbiome in mediating host SMR responses to changing environmental conditions.

The observed differences in gut microbial EAA provisioning function were not associated with a significant difference in gut bacterial community composition in D. punctata females, despite the difference in diets and the duration of the feeding trial. This lack of significant difference between dietary groups adds to the growing body of work detailing the variable impacts of dietary quality on cockroach gut microbiomes. These have ranged from significant differences in composition after 9 days (Pérez-Cobas et al., 2015) and 14 days (Bertino-Grimaldi et al., 2013) of feeding, to a relatively stable gut microbiome following ∼90 days of incubation on diets varying in cellulosic content (Schauer et al., 2014; Tinker and Ottesen, 2016). However, we observed the emergence of differentially abundant and metabolically distinct taxa between DF-fed and CADF-fed females that corroborate the proposed microbiome-mediated SMR response framework in this study.

In CADF-fed females, we detected higher mean relative abundances of taxa more commonly abundant in wood-feeding termites (Hongoh et al., 2005; Köhler et al., 2012; Mikaelyan et al., 2016) and wood-feeding cockroaches (Berlanga et al., 2016). For example, in response to the presence of cellulose in the CADF diet, cellulolytic and anaerobic nitrogen-fixing Firmicutes (family Clostridiaceae_1), anaerobic Patescibacteria (families Unassigned JGI_0000069-P22 and Absconditabacteriales_SR1) that utilize citrate, malate and amino acids as both generators and recycling agents of reductive power (Pepe-Ranney et al., 2016; Sieber et al., 2019), and bacterial endosymbionts of protists Proteobacteria (family Paracaedibacteraceae) (Tashyreva et al., 2018; Boscaro et al., 2019) were more abundant in CADF-fed females relative to DF-fed females. These higher abundances may be indicative of different diet-induced gut conditions in CADF females, leading to the emergence of taxa that are essential to thrive on high-fiber diets. In particular, the increase in relative abundance of protist-associated endosymbiotic bacteria Paracaedibacteraceae may be reflective of increased abundances of gut protists in CADF-fed females. High-cellulose diets have been shown to increase the abundances of specific cellulolytic termite gut protists (Pyrsonympha sp. and Microjoenia sp.), which tend to be less abundant on more complex untreated wood diets (Duarte et al., 2017), as well as changes in protist and protist-associated bacterial symbiont abundances (Gijzen et al., 1994; Berlanga et al., 2018). Interestingly, this bacterial family has so far only been reported in free-living freshwater and marine protists (Tashyreva et al., 2018; Boscaro et al., 2019). The relationship between this bacterial species and other protist-associated bacteria from wood-feeding termites (Stingl et al., 2005; Sato et al., 2009) and cockroaches (Noda et al., 2006) remains to be investigated. In contrast, the abundant bacterial taxa in the DF-fed females were those typically detected in high abundances in other laboratory-reared cockroaches (Schauer et al., 2012, 2014; Bertino-Grimaldi et al., 2013; Tinker and Ottesen, 2016; Ayayee et al., 2017; Mohammed et al., 2018; Table 3). These bacterial members are similarly present in wild cockroaches, although with comparatively different abundances. For example, the phyla Bacteroidetes and Fusobacteria usually have lower representation in wild cockroaches, whereas the phyla Proteobacteria and Firmicutes have relatively higher representation. The higher relative abundances of these taxa in DF-fed females may be contributing to the higher bacterial EAA input.

Conclusions

In conclusion, we found that a dietary shift from an optimal diet (DF) to a relatively sub-optimal diet (CADF) for 28 days resulted in a higher SMR in CADF-fed D. punctata females relative to DF-fed females. We propose that the higher SMR in CADF-fed females may be attributed to diet-induced changes in microbial functions, such as lowered provisioning of microbe-derived metabolites (acetate, lactate, propionate and EAAs) for energy generation and biosynthesis. Our results support our hypothesis that such host gut microbiota functions might mediate the metabolic phenotype of insect hosts. Additionally, we found that changes in relatively low abundant taxa in the microbiome can significantly impact metabolic function, which may not necessarily be associated with significant changes in the whole microbiome composition. The results herein provide a framework within which to investigate the mediating roles of an organism's gut microbiome on several aspects of host metabolic phenotype in response to external stressors, such as diet, disease and environment.

The authors thank Jose Diaz at the Evolution, Ecology, and Organismal Biology Department, The Ohio State University, and also the University of Wyoming's Microbial Ecology Collaborative.

Author contributions

Conceptualization: P.A.A., G.K., C.Y., T.L., A.M.; Methodology: P.A.A., T.L., L.T.A.v.D., A.M.-G.; Validation: P.A.A., C.Y.; Formal analysis: P.A.A., C.Y., G.F.C., T.L., L.T.A.v.D., A.M.-G.; Resources: G.K., C.Y., L.T.A.v.D., A.M.-G.; Writing - original draft: P.A.A., T.L., A.M.-G.; Writing - review & editing: P.A.A., G.K., C.Y., G.F.C., T.L., L.T.A.v.D., A.M.-G.; Supervision: P.A.A., L.T.A.v.D., A.M.-G.; Funding acquisition: L.T.A.v.D., A.M.-G.

Funding

This research was supported by the University of Wyoming's Microbial Ecology Collaborative (National Science Foundation award EPS-1655726).

Anderson
,
M. J.
(
2017
).
Permutational multivariate analysis of variance (PERMANOVA)
.
Wiley StatsRef Stat. Ref. Online
.
Apprill
,
A.
,
McNally
,
S.
,
Parsons
,
R.
and
Weber
,
L.
(
2015
).
Minor revision to V4 region SSU rRNA 806R gene primer greatly increases detection of SAR11 bacterioplankton
.
Aquat. Microb. Ecol.
75
,
129
-
137
.
Ayayee
,
P. A.
,
Keeney
,
G.
,
Sabree
,
Z. L.
and
Muñoz-Garcia
,
A.
(
2017
).
Compositional differences among female-associated and embryo-associated microbiota of the viviparous Pacific Beetle cockroach, Diploptera punctata
.
FEMS Microbiol. Ecol.
93
,
fix052
.
Ayayee
,
P. A.
,
Larsen
,
T.
,
Rosa
,
C.
,
Felton
,
G. W.
,
Ferry
,
J. G.
and
Hoover
,
K.
(
2016a
).
Essential amino acid supplementation by gut microbes of a wood-feeding cerambycid
.
Environ. Entomol.
45
,
66
-
73
.
Ayayee
,
P. A.
,
Larsen
,
T.
and
Sabree
,
Z.
(
2016b
).
Symbiotic essential amino acids provisioning in the American cockroach, Periplaneta americana (Linnaeus) under various dietary conditions
.
PeerJ
4
,
e2046
.
Ayayee
,
P. A.
,
Ondrejech
,
A.
,
Keeney
,
G.
and
Munõz-Garcia
,
A.
(
2018
).
The role of gut microbiota in the regulation of standard metabolic rate in female Periplaneta americana
.
PeerJ
6
,
e4717
.
Behmer
,
S. T.
(
2008
).
Insect herbivore nutrient regulation
.
Annu. Rev. Entomol.
54
,
165
-
187
.
Bell
,
B.
,
Roth
,
L.
and
Nalepa
,
C.
(
2007
).
Cockroaches: Ecology, Behavior, and Natural History
.
Baltimoire, MD
,
USA
:
Johns Hopkins University Press
.
Berlanga
,
M.
,
Llorens
,
C.
,
Comas
,
J.
and
Guerrero
,
R.
(
2016
).
Gut bacterial community of the xylophagous cockroaches Cryptocercus punctulatus and Parasphaeria boleiriana
.
PLoS ONE
11
,
e0152400
.
Berlanga
,
M.
,
Palau
,
M.
and
Guerrero
,
R.
(
2018
).
Gut microbiota dynamics and functionality in Reticulitermes grassei after a 7-day dietary shift and ciprofloxacin treatment
.
PLoS ONE
13
,
e0209789
-
e0209789
.
Bertino-Grimaldi
,
D.
,
Medeiros
,
M. N.
,
Vieira
,
R. P.
,
Cardoso
,
A. M.
,
Turque
,
A. S.
,
Silveira
,
C. B.
,
Albano
,
R. M.
,
Bressan-Nascimento
,
S.
,
Garcia
,
E. S.
,
de Souza
,
W.
, et al. 
(
2013
).
Bacterial community composition shifts in the gut of Periplaneta americana fed on different lignocellulosic materials
.
Springerplus
2
,
609
.
Bignell
,
D. E.
(
1977
).
An experimental study of cellulose and hemicellulose degradation in the alimentary canal of the American cockroach
.
Can. J. Zool.
55
,
579
-
589
.
Bignell
,
D. E.
(
1978
).
Effects of cellulose in the diets of cockroaches
.
Entomol. Exp. Appl.
24
,
254
-
257
.
Boscaro
,
V.
,
Husnik
,
F.
,
Vannini
,
C.
and
Keeling
,
P. J.
(
2019
).
Symbionts of the ciliate Euplotes: diversity, patterns and potential as models for bacteria–eukaryote endosymbioses
.
Proc. R. Soc. B Biol. Sci.
286
,
20190693
.
Bray
,
J. R.
and
Curtis
,
J. T.
(
1957
).
An ordination of the upland forest communities of Southern Wisconsin
.
Ecol. Monogr.
27
,
325
-
349
.
Breznak
,
J. A.
and
Kane
,
M. D.
(
1990
).
Microbial H2/CO2 acetogenesis in animal guts: nature and nutritional significance
.
FEMS Microbiol. Lett.
87
,
309
-
314
.
Burton
,
T.
,
Killen
,
S. S.
,
Armstrong
,
J. D.
and
Metcalfe
,
N. B.
(
2011
).
What causes intraspecific variation in resting metabolic rate and what are its ecological consequences?
Proc. Biol. Sci.
278
,
3465
-
3473
.
Callahan
,
B. J.
,
McMurdie
,
P. J.
,
Rosen
,
M. J.
,
Han
,
A. W.
,
Johnson
,
A. J. A.
and
Holmes
,
S. P.
(
2016
).
DADA2: High-resolution sample inference from Illumina amplicon data
.
Nat. Methods
13
,
581
-
583
.
Caporaso
,
J. G.
,
Kuczynski
,
J.
,
Stombaugh
,
J.
,
Bittinger
,
K.
,
Bushman
,
F. D.
,
Costello
,
E. K.
,
Fierer
,
N.
,
Peña
,
A. G.
,
Goodrich
,
J. K.
,
Gordon
,
J. I.
, et al. 
(
2010
).
QIIME allows analysis of high-throughput community sequencing data
.
Nat. Methods
7
,
335
-
336
.
Chao
,
A.
(
1984
).
Non-parametric estimation of the number of classes in a population
.
Scand. J. Stat.
11
,
265
-
270
.
Chen
,
W.-P.
,
Yang
,
X.-Y.
,
Hegeman
,
A. D.
,
Gray
,
W. M.
and
Cohen
,
J. D.
(
2010
).
Microscale analysis of amino acids using gas chromatography–mass spectrometry after methyl chloroformate derivatization
.
J. Chromatogr. B
878
,
2199
-
2208
.
Chown
,
S. L.
and
Gaston
,
K. J.
(
1999
).
Exploring links between physiology and ecology at macro-scales: the role of respiratory metabolism in insects
.
Biol. Rev.
74
,
87
-
120
.
Clark
,
R. M.
,
Zera
,
A. J.
and
Behmer
,
S. T.
(
2016
).
Metabolic rate is canalized in the face of variable life history and nutritional environment
.
Funct. Ecol.
30
,
922
-
931
.
Colman
,
D. R.
,
Toolson
,
E. C.
and
Takacs-Vesbach
,
C. D.
(
2012
).
Do diet and taxonomy influence insect gut bacterial communities?
Mol. Ecol.
21
,
5124
-
5137
.
Cruz-Neto
,
A. P.
and
Bozinovic
,
F.
(
2004
).
The relationship between diet quality and basal metabolic rate in endotherms: Insights from intraspecific analysis
.
Physiol. Biochem. Zool.
77
,
877
-
889
.
D'Argenio
,
V.
,
Casaburi
,
G.
,
Precone
,
V.
and
Salvatore
,
F.
(
2014
).
Comparative metagenomic analysis of human gut microbiome composition using two different bioinformatic pipelines
.
Biomed Res. Int.
2014
,
325340
.
DeVries
,
Z. C.
,
Kells
,
S. A.
and
Appel
,
A. G.
(
2013
).
Standard metabolic rate of the bed bug, Cimex lectularius: effects of temperature, mass, and life stage
.
J. Insect Physiol.
59
,
1133
-
1139
.
Dietrich
,
C.
,
Köhler
,
T.
and
Brune
,
A.
(
2014
).
The cockroach origin of the termite gut microbiota: patterns in bacterial community structure reflect major evolutionary events
.
Appl. Environ. Microbiol.
80
,
2261
-
2269
.
Dingha
,
B.
(
2009
).
Effects of temperature on the metabolic rates of insecticide resistant and susceptible German cockroaches, Blattella germanica (L.) (Dictyoptera: Blattellidae)
.
Midsouth Entomol.
2
,
17
-
27
.
Duarte
,
S.
,
Duarte
,
M.
,
Borges
,
P. A. V.
and
Nunes
,
L.
(
2017
).
Dietary-driven variation effects on the symbiotic flagellate protist communities of the subterranean termite Reticulitermes grassei Clément
.
J. Appl. Entomol.
141
,
300
-
307
.
Edgar
,
R. C.
(
2010
).
Search and clustering orders of magnitude faster than BLAST
.
Bioinformatics
26
,
2460
-
2461
.
Geib
,
S. M.
,
Jimenez-Gasco
,
M. D. M.
,
Carlson
,
J. E.
,
Tien
,
M.
and
Hoover
,
K.
(
2009
).
Effect of host tree species on cellulase activity and bacterial community composition in the gut of larval Asian longhorned beetle
.
Environ. Entomol.
38
,
686
-
699
.
Gier
,
H.
(
1936
).
The morphology and behavior of the intracellular bacteroids of roaches
.
Biol. Bull.
71
,
433
-
452
.
Gijzen
,
H. J.
,
van der Drift
,
C.
,
Barugahare
,
M.
and
op den Camp
,
H. J. M.
(
1994
).
Effect of host diet and hindgut microbial composition on cellulolytic activity in the hindgut of the American cockroach, Periplaneta americana
.
Appl. Environ. Microbiol.
60
,
1822
-
1826
.
Harrison
,
J. F.
and
Fewell
,
J. H.
(
2016
).
Thermal effects on feeding behavior and net energy intake in a grasshopper experiencing large diurnal fluctuations in body temperature
.
Physiol. Zool.
68
,
453
-
473
.
Hongoh
,
Y.
,
Deevong
,
P.
,
Inoue
,
T.
,
Moriya
,
S.
,
Trakulnaleamsai
,
S.
,
Ohkuma
,
M.
,
Vongkaluang
,
C.
,
Noparatnaraporn
,
N.
and
Kudo
,
T.
(
2005
).
Intra- and interspecific comparisons of bacterial diversity and community structure support coevolution of gut microbiota and termite host
.
Appl. Environ. Microbiol.
71
,
6590
-
6599
.
Jehrke
,
L.
,
Stewart
,
F. A.
,
Droste
,
A.
and
Beller
,
M.
(
2018
).
The impact of genome variation and diet on the metabolic phenotype and microbiome composition of Drosophila melanogaster
.
Sci. Rep.
8
,
6215
.
Jensen
,
K.
,
Mayntz
,
D.
,
Wang
,
T.
,
Simpson
,
S. J.
and
Overgaard
,
J.
(
2010
).
Metabolic consequences of feeding and fasting on nutritionally different diets in the wolf spider Pardosa prativaga
.
J. Insect Physiol.
56
,
1095
-
1100
.
Kakumanu
,
M. L.
,
Maritz
,
J. M.
,
Carlton
,
J. M.
and
Schal
,
C.
(
2018
).
Overlapping community compositions of gut and fecal microbiomes in lab-reared and field-collected German cockroaches
.
Appl. Environ. Microbiol.
84
,
e01037-18
.
Kane
,
M. D.
and
Breznak
,
J. A.
(
1991
).
Effect of host diet on production of organic acids and methane by cockroach gut bacteria
.
Appl. Environ. Microbiol.
57
,
2628
-
2634
.
Karasov
,
W. H.
(
1986
).
Energetics, physiology and vertebrate ecology
.
Trends Ecol. Evol.
1
,
101
-
104
.
Karowe
,
D. N.
and
Martin
,
M. M.
(
1989
).
The effects of quantity and quality of diet nitrogen on the growth, efficiency of food utilization, nitrogen budget, and metabolic rate of fifth-instar Spodoptera eridania larvae (Lepidoptera: Noctuidae)
.
J. Insect Physiol.
35
,
699
-
708
.
Kim
,
J. M.
,
Choi
,
M.-Y.
,
Kim
,
J.-W.
,
Lee
,
S. A.
,
Ahn
,
J.-H.
,
Song
,
J.
,
Kim
,
S.-H.
and
Weon
,
H.-Y.
(
2017
).
Effects of diet type, developmental stage, and gut compartment in the gut bacterial communities of two Cerambycidae species (Coleoptera)
.
J. Microbiol.
55
,
21
-
30
.
Köhler
,
T.
,
Dietrich
,
C.
,
Scheffrahn
,
R. H.
and
Brune
,
A.
(
2012
).
High-resolution analysis of gut environment and bacterial microbiota reveals functional compartmentation of the gut in wood-feeding higher termites (Nasutitermes spp.)
.
Appl. Environ. Microbiol.
78
,
4691
-
4701
.
Larsen
,
T.
,
Pollierer
,
M. M.
,
Holmstrup
,
M.
,
D'Annibale
,
A.
,
Maraldo
,
K.
,
Andersen
,
N.
and
Eriksen
,
J.
(
2016b
).
Substantial nutritional contribution of bacterial amino acids to earthworms and enchytraeids: a case study from organic grasslands
.
Soil Biol. Biochem.
99
,
21
-
27
.
Larsen
,
T.
,
Taylor
,
D. L.
,
Leigh
,
M. B.
and
O'Brien
,
D. M.
(
2009
).
Stable isotope fingerprinting: a novel method for identifying plant, fungal, or bacterial origins of amino acids
.
Ecology
90
,
3526
-
3535
.
Larsen
,
T.
,
Ventura
,
M.
,
Andersen
,
N.
,
O'Brien
,
D. M.
,
Piatkowski
,
U.
and
McCarthy
,
M. D.
(
2013
).
Tracing carbon sources through aquatic and terrestrial food webs using amino acid stable isotope fingerprinting
.
PLoS ONE
8
,
e73441
.
Larsen
,
T.
,
Ventura
,
M.
,
Maraldo
,
K.
,
Triadó-Margarit
,
X.
,
Casamayor
,
E. O.
,
Wang
,
Y. V.
,
Andersen
,
N.
and
O'Brien
,
D. M.
(
2016a
).
The dominant detritus-feeding invertebrate in Arctic peat soils derives its essential amino acids from gut symbionts
.
J. Anim. Ecol.
85
,
1275
-
1285
.
Larsen
,
T.
,
Ventura
,
M.
,
O'Brien
,
D. M.
,
Magid
,
J.
,
Lomstein
,
B. A.
and
Larsen
,
J.
(
2011
).
Contrasting effects of nitrogen limitation and amino acid imbalance on carbon and nitrogen turnover in three species of Collembola
.
Soil Biol. Biochem.
43
,
749
-
759
.
Lee
,
K. P.
(
2007
).
The interactive effects of protein quality and macronutrient imbalance on nutrient balancing in an insect herbivore
.
J. Exp. Biol.
210
,
3236
-
3244
.
Lighton
,
J.
(
2008
).
Measuring Metabolic Rates: A Manual for Scientists
.
Oxford University Press
.
Martin
,
M.
(
2011
).
Cutadapt removes adapter sequences from high-throughput sequencing reads
.
EMBnet. J.
17
,
10
.
McCue
,
M. D.
,
Guzman
,
R. M.
,
Passement
,
C. A.
and
Davidowitz
,
G.
(
2015
).
How and when do insects rely on endogenous protein and lipid resources during lethal bouts of starvation? A new application for 13C-breath testing
.
PLoS ONE
10
,
e0140053
.
McMahon
,
K. W.
,
Fogel
,
M. L.
,
Elsdon
,
T. S.
and
Thorrold
,
S. R.
(
2010
).
Carbon isotope fractionation of amino acids in fish muscle reflects biosynthesis and isotopic routing from dietary protein
.
J. Anim. Ecol.
79
,
1132
-
1141
.
Mielke
,
P. W.
(
1984
).
34 Meteorological applications of permutation techniques based on distance functions
.
Handb. Stat.
4
,
813
-
830
.
Mikaelyan
,
A.
,
Meuser
,
K.
and
Brune
,
A.
(
2016
).
Microenvironmental heterogeneity of gut compartments drives bacterial community structure in wood- and humus-feeding higher termites
.
FEMS Microbiol. Ecol.
93
,
fiw210
.
Mohammed
,
W. S.
,
Ziganshina
,
E. E.
,
Shagimardanova
,
E. I.
,
Gogoleva
,
N. E.
and
Ziganshin
,
A. M.
(
2018
).
Comparison of intestinal bacterial and fungal communities across various xylophagous beetle larvae (Coleoptera: Cerambycidae)
.
Sci. Rep.
8
,
10073
.
Naya
,
D. E.
,
Lardies
,
M. A.
and
Bozinovic
,
F.
(
2007
).
The effect of diet quality on physiological and life-history traits in the harvestman Pachylus paessleri
.
J. Insect Physiol.
53
,
132
-
138
.
Newsome
,
S. D.
,
Fogel
,
M. L.
,
Kelly
,
L.
and
del Rio
,
C. M.
(
2011
).
Contributions of direct incorporation from diet and microbial amino acids to protein synthesis in Nile tilapia
.
Funct. Ecol.
25
,
1051
-
1062
.
Noda
,
S.
,
Inoue
,
T.
,
Hongoh
,
Y.
,
Kawai
,
M.
,
Nalepa
,
C. A.
,
Vongkaluang
,
C.
,
Kudo
,
T.
and
Ohkuma
,
M.
(
2006
).
Identification and characterization of ectosymbionts of distinct lineages in Bacteroidales attached to flagellated protists in the gut of termites and a wood-feeding cockroach
.
Environ. Microbiol.
8
,
11
-
20
.
Parada
,
A. E.
,
Needham
,
D. M.
and
Fuhrman
,
J. A.
(
2016
).
Every base matters: assessing small subunit rRNA primers for marine microbiomes with mock communities, time series and global field samples
.
Environ. Microbiol.
18
,
1403
-
1414
.
Pepe-Ranney
,
C.
,
Koechli
,
C.
,
Potrafka
,
R.
,
Andam
,
C.
,
Eggleston
,
E.
,
Garcia-Pichel
,
F.
and
Buckley
,
D. H.
(
2016
).
Non-cyanobacterial diazotrophs mediate dinitrogen fixation in biological soil crusts during early crust formation
.
ISME J.
10
,
287
-
298
.
Pérez-Cobas
,
A. E.
,
Maiques
,
E.
,
Angelova
,
A.
,
Carrasco
,
P.
,
Moya
,
A.
and
Latorre
,
A.
(
2015
).
Diet shapes the gut microbiota of the omnivorous cockroach Blattella germanica
.
FEMS Microbiol. Ecol.
91
,
fiv022
.
Pollierer
,
M. M.
,
Larsen
,
T.
,
Potapov
,
A.
,
Brückner
,
A.
,
Heethoff
,
M.
,
Dyckmans
,
J.
and
Scheu
,
S.
(
2019
).
Compound-specific isotope analysis of amino acids as a new tool to uncover trophic chains in soil food webs
.
Ecol. Monogr.
89
,
e01384
.
Potapov
,
A. M.
,
Tiunov
,
A. V.
,
Scheu
,
S.
,
Larsen
,
T.
and
Pollierer
,
M. M.
(
2019
).
Combining bulk and amino acid stable isotope analyses to quantify trophic level and basal resources of detritivores: a case study on earthworms
.
Oecologia
189
,
447
-
460
.
Quast
,
C.
,
Pruesse
,
E.
,
Yilmaz
,
P.
,
Gerken
,
J.
,
Schweer
,
T.
,
Yarza
,
P.
,
Peplies
,
J.
and
Glöckner
,
F. O.
(
2012
).
The SILVA ribosomal RNA gene database project: improved data processing and web-based tools
.
Nucleic Acids Res.
41
,
D590
-
D596
.
Raychoudhury
,
R.
,
Sen
,
R.
,
Cai
,
Y.
,
Sun
,
Y.
,
Lietze
,
V.-U.
,
Boucias
,
D. G.
and
Scharf
,
M. E.
(
2013
).
Comparative metatranscriptomic signatures of wood and paper feeding in the gut of the termite Reticulitermes flavipes (Isoptera: Rhinotermitidae)
.
Insect Mol. Biol.
22
,
155
-
171
.
Roeder
,
K. A.
and
Behmer
,
S. T.
(
2014
).
Lifetime consequences of food protein-carbohydrate content for an insect herbivore
.
Funct. Ecol.
28
,
1135
-
1143
.
Rossmassler
,
K.
,
Dietrich
,
C.
,
Thompson
,
C.
,
Mikaelyan
,
A.
,
Nonoh
,
J. O.
,
Scheffrahn
,
R. H.
,
Sillam-Dussès
,
D.
and
Brune
,
A.
(
2015
).
Metagenomic analysis of the microbiota in the highly compartmented hindguts of six wood- or soil-feeding higher termites
.
Microbiome
3
,
56
.
Sato
,
T.
,
Hongoh
,
Y.
,
Noda
,
S.
,
Hattori
,
S.
,
Ui
,
S.
and
Ohkuma
,
M.
(
2009
).
Candidatus desulfovibrio trichonymphae, a novel intracellular symbiont of the flagellate Trichonympha agilis in termite gut
.
Environ. Microbiol.
11
,
1007
-
1015
.
Schauer
,
C.
,
Thompson
,
C. L.
and
Brune
,
A.
(
2012
).
The bacterial community in the gut of the cockroach Shelfordella lateralis reflects the close evolutionary relatedness of cockroaches and termites
.
Appl. Environ. Microbiol.
78
,
2758
-
2767
.
Schauer
,
C.
,
Thompson
,
C.
and
Brune
,
A.
(
2014
).
Pyrotag sequencing of the gut microbiota of the cockroach Shelfordella lateralis reveals a highly dynamic core but only limited effects of diet on community structure
.
PLoS ONE
9
,
e85861
.
Schimpf
,
N. G.
,
Matthews
,
P. G. D.
and
White
,
C. R.
(
2012a
).
Cockroaches that exchange respiratory gases discontinuously survive food and water restriction
.
Evolution
66
,
597
-
604
.
Schimpf
,
N. G.
,
Matthews
,
P. G. D.
and
White
,
C. R.
(
2012b
).
Standard metabolic rate is associated with gestation duration, but not clutch size, in speckled cockroaches Nauphoeta cinerea
.
Biol. Open
1
,
1185
-
1191
.
Schmidt-Nielsen
,
K.
(
1995
).
Animal Physiology
.
New York
:
Cambridge University Press
.
Scully
,
E. D.
,
Geib
,
S. M.
,
Carlson
,
J. E.
,
Tien
,
M.
,
McKenna
,
D.
and
Hoover
,
K.
(
2014
).
Functional genomics and microbiome profiling of the Asian longhorned beetle (Anoplophora glabripennis) reveal insights into the digestive physiology and nutritional ecology of wood feeding beetles
.
BMC Genomics
15
,
1096
.
Shannon
,
C. E.
(
1957
).
A mathematical theory of communication
.
Bell Syst. Tech. J.
27
,
379
-
423
.
Shapiro
,
S. S.
and
Wilk
,
M. B.
(
1965
).
An analysis of variance test for normality (complete samples)
.
Biometrika
52
,
591
-
611
.
Sieber
,
C. M. K.
,
Paul
,
B. G.
,
Castelle
,
C. J.
,
Hu
,
P.
,
Tringe
,
S. G.
,
Valentine
,
D. L.
,
Andersen
,
G. L.
and
Banfield
,
J. F.
(
2019
).
Unusual metabolism and hypervariation in the genome of a gracilibacterium (BD1-5) from an oil-degrading community
.
MBio
10
,
e02128-19
.
Simpson
,
E. H.
(
1949
).
Measurement of diversity
.
Nature
163
,
688
.
Stingl
,
U.
,
Radek
,
R.
,
Yang
,
H.
and
Brune
,
A.
(
2005
).
‘Endomicrobia’: cytoplasmic symbionts of termite gut protozoa form a separate phylum of prokaryotes
.
Appl. Environ. Microbiol.
71
,
1473
-
1479
.
Su
,
L. J.
,
Yang
,
L. L.
,
Huang
,
S.
,
Su
,
X. Q.
,
Li
,
Y.
,
Wang
,
F. Q.
,
Wang
,
E. T.
,
Kang
,
N.
,
Xu
,
J.
and
Song
,
A. D.
(
2016
).
Comparative gut microbiomes of four species representing the higher and the lower termites
.
J. Insect Sci.
16
,
97
.
Tashyreva
,
D.
,
Prokopchuk
,
G.
,
Votýpka
,
J.
,
Yabuki
,
A.
,
Horák
,
A.
and
Lukeš
,
J.
(
2018
).
Life cycle, ultrastructure, and phylogeny of new diplonemids and their endosymbiotic bacteria
.
MBio
9
,
e02447-17
.
Tinker
,
K. A.
and
Ottesen
,
E. A.
(
2016
).
The core gut microbiome of the American cockroach, Periplaneta americana, is stable and resilient to dietary shifts
.
Appl. Environ. Microbiol.
82
,
6603
-
6610
.
Turner
,
S.
,
Pryer
,
K. M.
,
Miao
,
V. P. W.
and
Palmer
,
J. D.
(
1999
).
Investigating deep phylogenetic relationships among cyanobacteria and plastids by small subunit rRNA sequence analysis
.
J. Eukaryot. Microbiol.
46
,
327
-
338
.
Walsh
,
R. G.
,
He
,
S.
and
Yarnes
,
C. T.
(
2014
).
Compound-specific δ13C and δ15N analysis of amino acids: a rapid, chloroformate-based method for ecological studies
.
Rapid Commun. Mass Spectrom.
28
,
96
-
108
.
Walters
,
W.
,
Hyde
,
E. R.
,
Berg-Lyons
,
D.
,
Ackermann
,
G.
,
Humphrey
,
G.
,
Parada
,
A.
,
Gilbert
,
J. A.
,
Jansson
,
J. K.
,
Caporaso
,
J. G.
,
Fuhrman
,
J. A.
, et al. 
(
2016
).
Improved bacterial 16S rRNA Gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys
.
mSystems
1
,
e00009-15
.
Williford
,
A.
,
Stay
,
B.
and
Bhattacharya
,
D.
(
2004
).
Evolution of a novel function: nutritive milk in the viviparous cockroach, Diploptera punctata
.
Evol. Dev.
6
,
67
-
77
.
Withers
,
P. C.
(
1977
).
Measurement of VO2, VCO2, and evaporative water loss with a flow-through mask
.
J. Appl. Physiol.
42
,
120
-
123
.
Wong
,
C. N. A.
,
Ng
,
P.
and
Douglas
,
A. E.
(
2011
).
Low-diversity bacterial community in the gut of the fruitfly Drosophila melanogaster
.
Environ. Microbiol.
13
,
1889
-
1900
.
Yarnes
,
C. T.
and
Herszage
,
J.
(
2017
).
The relative influence of derivatization and normalization procedures on the compound-specific stable isotope analysis of nitrogen in amino acids
.
Rapid Commun. Mass Spectrom.
31
,
693
-
704
.
Yun
,
J.-H.
,
Roh
,
S. W.
,
Whon
,
T. W.
,
Jung
,
M.-J.
,
Kim
,
M.-S.
,
Park
,
D.-S.
,
Yoon
,
C.
,
Nam
,
Y.-D.
,
Kim
,
Y.-J.
,
Choi
,
J.-H.
, et al. 
(
2014
).
Insect gut bacterial diversity determined by environmental habitat, diet, developmental stage, and phylogeny of host
.
Appl. Environ. Microbiol.
80
,
5254
-
5264
.
Zanotto
,
F.
,
Simpson
,
S. J.
and
Raubenheimer
,
D.
(
1993
).
The regulation of growth by locusts through post-ingestive compensation for variation in the levels of dietary protein and carbohydrate
.
Physiol. Entomol.
18
,
425
-
434
.
Zanotto
,
F.
,
Gouveia
,
S.
,
Simpson
,
S.
and
Calder
,
D.
(
1997
).
Nutritional homeostasis in locusts: is there a mechanism for increased energy expenditure during carbohydrate overfeeding?
J. Exp. Biol.
200
,
2437
-
2448
.

Competing interests

The authors declare no competing or financial interests.

Supplementary information