Skip to main content
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
  • ASM
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Eukaryotic Cell
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
Applied and Environmental Microbiology
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Current Issue
    • Accepted Manuscripts
    • COVID-19 Special Collection
    • Archive
    • Minireviews
  • For Authors
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics Resources and Policies
  • About the Journal
    • About AEM
    • Editor in Chief
    • Editorial Board
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • Subscribe
    • Members
    • Institutions
Invertebrate Microbiology

Experimental Manipulation Shows a Greater Influence of Population than Dietary Perturbation on the Microbiome of Tyrophagus putrescentiae

Tomas Erban, Ondrej Ledvinka, Marta Nesvorna, Jan Hubert
Patrick D. Schloss, Editor
Tomas Erban
aCrop Research Institute, Prague, Czechia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Tomas Erban
Ondrej Ledvinka
bCzech Hydrometeorological Institute, Prague, Czechia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Ondrej Ledvinka
Marta Nesvorna
aCrop Research Institute, Prague, Czechia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Jan Hubert
aCrop Research Institute, Prague, Czechia
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Jan Hubert
Patrick D. Schloss
University of Michigan
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/AEM.00128-17
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

ABSTRACT

Tyrophagus putrescentiae is inhabited by bacteria that differ among mite populations (strains) and diets. Here, we investigated how the microbiome and fitness of T. putrescentiae are altered by dietary perturbations and mite populations. Four T. putrescentiae populations, referred to as dog, Koppert, laboratory, and Phillips, underwent a perturbation, i.e., a dietary switch from a rearing diet to two experimental diets. The microbiome was investigated by sequencing the V1-V3 portion of the 16S rRNA gene, and selected bacterial taxa were quantified by quantitative PCR (qPCR) using group/taxon-specific primers. The parameters observed were the changes in mite population growth and nutritional status, i.e., the total glycogen, lipid, saccharide, and protein contents in mites. The effect of diet perturbation on the variability of the microbiome composition and population growth was lower than the effect induced by mite population. In contrast, the diet perturbation showed a greater effect on nutritional status of mites than the mite population. The endosymbionts exhibited high variations among T. putrescentiae populations, including Cardinium in the laboratory population, Blattabacterium-like bacteria in the dog population, and Wolbachia in the dog and Phillips populations. Solitalea-like and Bartonella-like bacteria were present in the dog, Koppert, and Phillips populations in different proportions. The T. putrescentiae microbiome is dynamic and varies based on both the mite population and perturbation; however, the mites remain characterized by robust bacterial communities. Bacterial endosymbionts were found in all populations but represented a dominant portion of the microbiome in only some populations.

IMPORTANCE We addressed the question of whether population origin or perturbation exerts a more significant influence on the bacterial community of the stored product mite Tyrophagus putrescentiae. The microbiomes of four populations of T. putrescentiae insects subjected to diet perturbation were compared. Based on our results, the bacterial community was more affected by the mite population than by diet perturbation. This result can be interpreted as indicating high stability of the putative intracellular symbionts in response to dietary perturbation. The changes in the absolute and relative numbers of Wolbachia, Blattabacterium-like, Solitalea-like, and Cardinium bacteria in the T. putrescentiae populations can also be caused by neutral processes other than perturbation. When nutritional status is considered, the effect of population appeared less important than the perturbation. We hypothesize that differences in the proportions of the endosymbiotic bacteria result in changes in mite population growth.

INTRODUCTION

Symbiotic intracellular bacteria have been traditionally classified as primary, including ancient and obligate bacteria, and secondary, including recent and facultative bacteria (1). The primary symbionts are those adapted to the host and vice versa. The associations between secondary symbionts and hosts are facultative although secondary symbionts are also adapted to their hosts. Moreover, secondary symbionts may be absent in individuals or populations (2). Variations in symbionts affect host structure and dynamics at the population level (3). The coexistence of infected and uninfected populations or the presence of populations with incompatible strains of symbionts indicates that symbionts form a metacommunity among the populations of insects (4). The endosymbiotic bacteria Wolbachia, Rickettsia, and Cardinium are present in 52, 24, and 13%, respectively, of all arthropod species (5). It was found that Wolbachia and Cardinium maintain prevalence within many arthropod taxa by implementing a variety of reproductive manipulations (6). Although Wolbachia and Cardinium are transmitted maternally, both taxa or their genes have been found to be transmitted horizontally among spider mites (7). More recently, horizontal Wolbachia transfer between Drosophila melanogaster stocks facilitated by Tyrophagus putrescentiae has been reported; moreover, this transmission occurred through the ingestion of Wolbachia in infested corpses (8). However, Rickettsia has not yet been identified in domestic mites (9, 10).

Approximately 3 decades ago, Bowman (11) demonstrated intraspecific differences in various populations of Acarus siro, Lepidoglyphus (Glycyphagus) destructor, Tyrophagus longior, and Rhizoglyphus sp. The intrapopulation differences in the physiological status of these mite populations in terms of their digestive enzymes were higher than the differences among species (11–13). Recently, differences in population growth were observed among four tested T. putrescentiae populations: three of the four populations exhibited density-dependent growth (14). One possible explanation for this discrepancy is the presence of different symbionts or differences in their proportions in mite populations. Analyses of the Sanger sequences of 16S rRNA genes from six T. putrescentiae populations support potential intraspecific physiological differences among populations of T. putrescentiae as a result of associated bacteria (15). The differences were caused by the endosymbionts Wolbachia, Cardinium, and Blattabacterium-like, and the putative symbiont Bartonella-like bacteria. Solitalea-like bacteria were present in all tested populations of T. putrescentiae (15) and are suggested as putative symbionts in A. siro (16). The proportion of bacteriocytes localized in the T. putrescentiae fat body (17) was low in the observed individuals and increased after perturbation (diet change) (15). This effect was attributed to the perturbation associated with different nutrient requirements and subsequent acquisition of bacteria (Brevundimonas vesicularis, Stenotrophomonas maltophilia, Serratia liquefaciens, and Serratia marcescens) from the environment as well as bacteriocyte formation under certain conditions (18, 19). These endosymbionts are not primary symbionts in T. putrescentiae mites because they are not present in all individuals (18, 19).

In this study, we examined whether the T. putrescentiae bacterial community is influenced more by perturbation (diet switch) or by mite population and whether there are correlations between fitness indicators and microbiome composition in T. putrescentiae populations. To address these questions, we analyzed four populations of T. putrescentiae mites. Population growth and nutritional status were analyzed as fitness parameters in samples of adult/juvenile mites. The bacterial community was characterized by sequencing the V1-V3 fragment of the 16S rRNA gene and performing quantitative PCR (qPCR) with group/taxon-specific primers.

RESULTS

The bacterial communities of the T. putrescentiae populations differed.The bacterial community (see Table S3 in the supplemental material) was more influenced by population than perturbation (Tables 1 and 2), as indicated by both the Illumina and qPCR data sets. According to permutational multivariate analysis of variance (PERMANOVA) and redundancy analyses (RDAs), regardless of whether the Bray-Curtis or Jaccard distances were utilized (P ≤ 0.001 for all cases), the F values clearly document this finding. The two-way PERMANOVA and distance-based RDA (db-RDA) results using the Jaccard index differed from the analyses of the Bray-Curtis matrix (Tables 1 and 2). For instance, the role of interactions between populations and perturbation was slightly lower in the Jaccard space but remained significant in terms of P values (Table 1). The presence/absence of operational taxonomic units (OTUs) was more influenced by population, but certain OTUs were influenced by perturbation, as indicated by the Venn diagrams (Fig. S3 and S4).

View this table:
  • View inline
  • View popup
TABLE 1

Comparison of the effect of populations and the tested variables and their interactions on T. putrescentiae microbiome compositiona

View this table:
  • View inline
  • View popup
TABLE 2

The effect of perturbation on the microbiome composition of T. putrescentiae populations tested separately by one-way PERMANOVA and db-RDAa

The similarities of the samples and contribution of OTUs were visualized by triplots (Fig. 1A and B) of the Illumina and qPCR data (see below for environmental variables and model selection). In both data sets, the most important OTUs distinguished among the population belonged to endosymbiotic bacteria (Cardinium, Wolbachia, Solitalea-like, and Blattabacterium-like bacteria) and to Bartonella-like bacteria. For the Illumina data set, OTUs of members of Bacillaceae and Enterobacteriaceae were important, i.e., Staphylococcus cohnii (OTU 4), Bacillus cereus (OTU 9), and Xenorhabdus (OTU 14). The random forest algorithm (errors = 0) analyses indicated that the following OTUs were responsible for differences: (i) OTUS of Micrococcaceae with low similarity (90%) to Acaricomes (OTU 15), Enterobacteriaceae with low similarity (92%) to Xenorhabdus (OTU 14), Cardinium (OTU 3), and Kocuria koreensis (OTU 19) were characteristic of the laboratory population; (ii) OTUs of Staphylococcus lentus (OTU 16) and B. cereus (OTU 9) were characterized as specific to the Koppert population; (iii) OTUs of Wolbachia (OTU 1), Bartonella-like (OTU 2), and Solitalea-like (OTU 6) bacteria were characteristic of both the Phillips and dog populations; and (iv) an OTU of Blattabacterium-like bacteria (OTU 8) was associated with the dog population only (Table S5). The associations were visible in the Krona projections (Fig. S2).

FIG 1
  • Open in new tab
  • Download powerpoint
FIG 1

Triplot visualization of principal coordinates in redundancy analysis of the T. putrescentiae bacterial community. The analysis includes population, diet, and growth as environmental factors. A correlation triplot containing sample scores given by the weighted sums of OTUs was constructed. (A) Standardized Illumina data set transformed into Jaccard matrix. (B) qPCR data set transformed into a Euclidean matrix. The perturbation was coded as follows: 0, rearing diet, i.e., SPMd for all populations with the exception of the dog population (DDFd); 1, SPMd/DDFd or, for the dog population, DDFd/SPMd; 2, SPMd/OSRd or, for the dog population, DDFd/OSRd. Lab, laboratory population; Phil, Phillips population; Kop, Koppert population; P1, perturbation 1; P2, perturbation 2.

Perturbation of bacterial communities.A heat map (Fig. 2) was created to visualize the differences in the numbers of OTUs or taxa in the Illumina data set. There were differences in the log2-transformed data, indicating some contribution from low-abundance OTUs (Fig. S5 indicates cluster significance). When perturbation and population were incorporated into the random forest analyses, the error rate (0.222) was higher than for population only (Table S4). The most important OTUs were those of Micrococcaceae with low similarity to Acaricomes (OTU 15), Blattabacterium-like (OTU 8), S. cohnii (OTU 4), and Cardinium (OTU 3).

FIG 2
  • Open in new tab
  • Download powerpoint
FIG 2

Heat map representing the distribution of samples and OTUs in the T. putrescentiae bacterial community in the standardized Illumina subsample data set. Data were log2 transformed. The samples are characterized by different colors, as follows: red, Koppert; blue, laboratory; green, Phillips; violet, dog population. The legend for the sample codes is provided in Table S2, and the legend for the OTUs is provided in Table S3.

The perturbation had different effects on the microbiome depending on the mite population. The changes were random rather than systematic in the qPCR data set (Fig. 3). The number of Wolbachia bacteria was increased by perturbation in the Phillips, Koppert, and dog populations but decreased in the laboratory population. In contrast, Cardinium and Solitalea-like bacteria were not influenced by perturbation in the above populations, apart from Solitalea-like bacteria in the Koppert population (Fig. 3).

FIG 3
  • Open in new tab
  • Download powerpoint
FIG 3

Quantification of the bacterial community of T. putrescentiae based on qPCR with universal and taxon-specific primers (primers are given in Table S1). The numbers of copies were recalculated per single mite and log10 transformed for the universal and taxon-specific primers, as indicated. The perturbation was coded as follows: 0, rearing diet, i.e., SPMd for all populations with the exception of the dog population (DDFd); 1, SPMd/DDFd or, for the dog population, DDFd/SPMd; 2, SPMd/OSRd or, for the dog population, DDFd/ORSd.

Based on the Illumina data set, the laboratory population contained 13 core OTUs, and perturbation resulted in the addition of three OTUs: an OTU of Lactobacillus amylovorus (OTU 37) for perturbation 1 and two OTUs of Virgibacillus (OTUs 21 and 34) for perturbation 2. The dog population included nine core OTUs, and perturbation 1 resulted in a loss of three noncore OTUs (from L. amylovorus [OTU 37], Lactobacillus sp. [OTU 71], and Pantoea vagans [OTU 102]) and the acquisition of four OTUs (from Wolbachia sp. strain 2 [OTU 24], Prevotella paludivivens [OTU 45], Lactobacillus paracollinoides [OTU 106], and Pseudomonas putida [OTU 203]). The Phillips population included nine core species, and perturbation resulted in the disappearance of three OTUs that originally occurred with the control stored product mite diet (SPMd): P. paludivivens (OTU 45), L. paracollinoides (OTU 106), and Selenomonas lacticifex (OTU 202). Perturbation 1 provided three new OTUs: from L. amylovorus (OTU 37), Staphylococcus saprophyticus (OTU 38), and S. lentus (OTU 16). The Koppert population included 11 core OTUs, and perturbation resulted in the disappearance of OTUs 27, 30, and 31 (classified as Bacillaceae_2), whereas L. amylovorus (OTU 37), Lactobacillus sp. (OTU 71), and P. vagans (OTU 102) appeared after perturbation 1 (Fig. S4).

The observed F and P values resulting from the RDAs (Tables 3 and 4) suggest that the inverse Simpson diversity index was more influenced by perturbation than by population. The effect of the perturbation was apparent in all populations, but there were exceptions in the dog population (perturbation 1, dry dog food kernel diet switched to oilseed rape diet [DDFd/OSRd]) and Phillips population (perturbation 1, SPMd/DDFd). Alpha diversity appeared to be influenced more by perturbation, whereas the population showed a greater influence on beta diversity (Fig. 4).

View this table:
  • View inline
  • View popup
TABLE 3

Comparison of the effects of population and perturbation and their interaction on T. putrescentiae population growth and nutritional status

View this table:
  • View inline
  • View popup
TABLE 4

Comparison of the effects of population on T. putrescentiae growth and nutritional statusa

FIG 4
  • Open in new tab
  • Download powerpoint
FIG 4

Comparison of the observed parameters of the T. putrescentiae populations under the control condition and after diet switch. The points are means; the bars are standard deviations. (A) Inverse Simpson diversity index (1/D) of bacterial community characterized by the Illumina data set. (B) Population growth, expressed as the logarithms of numbers of individuals after 21 days of growth from 10 individuals. (C to F) The nutrient contents, as indicated, were recalculated per milligram of nutrient to gram of fresh weight of mite. The perturbation was coded as follows: 0, rearing diet, i.e., SPMd for all populations with the exception of the dog population (DDFd); 1, SPMd/DDFd or, for the dog population, DDFd/SPMd; 2, SPMd/OSRd or, for the dog population, DDFd/ORSd.

Fitness parameters.Based on the results of RDA applied to normalized growth, population clearly had a greater effect than diet (P < 0.001 versus P = 0.015). The Koppert population differed somewhat from the others (Fig. 4). Perturbation 2 exhibited almost no effect in the Koppert population although its influence was clearly apparent in the other populations (Tables 3 and 4). The interactions between population and perturbation were also significant (P < 0.001). RDA also suggested a substantial influence of both population and perturbation on the overall distribution of the studied nutrients (Table 3). However, perturbation contributed more heavily to the mutual effects of these two factors. Indeed, the respective F values for interaction, population, and perturbation were 4.92, 3.5, and 11.51. All P values were less than 0.005. When the effects on the status of each nutritional component (i.e., glycogen, lipids, proteins, and saccharides) were studied separately (Table 3), all influences (including their interactions) were significant at least at the 0.05 level, with the exceptions of population in the case of glycogen status (F = 2.35, P = 0.07) and the interaction between population and perturbation in the case of protein contents (F = 1.38, P = 0.233) (Fig. 4).

Bacteria associated with nutrient contents.The combination of the forward selection of explanatory variables and variance inflation factor (VIF) values greater than 10 indicated that it is not advisable to add nutritional status and growth to the basic db-RDA model including population and diet together. However, when we focused on perturbation instead of the diets themselves, the VIFs decreased remarkably. In the case of the Illumina data set, lipid status alone could be added to the basic model obtained using the Jaccard distances. In the case of the qPCR data set, glycogen content showed similar results. Therefore, triplots were constructed with these variables (population, perturbation, and lipids for the Illumina data set and population, perturbation, and glycogen status for the qPCR data set) (Fig. 1A and B). According to the best RDA model, population was responsible for the highest variability in the Illumina data set (adjusted R2, cumulative, 0.8), followed by lipid status (adjusted R2, cumulative, 0.814) and perturbation (adjusted R2, cumulative, 0.829). Additionally, the analysis of variance (ANOVA)-like permutation test applied to the terms in the RDA model confirmed the significant effects of population, lipid status, and perturbation and the interaction between population and perturbation (all with P ≤ 0.01) on the OTU distribution. All populations differed from each other in terms of bacterial composition, as shown in the triplot visualization of principal coordinates in RDA (Fig. 1A). Lipids were closely related to perturbation 2. RDA models were constructed to describe the interactions between environmental variables and absolute numbers of 16S rRNA gene copies (qPCR data). The combination of the forward selection of the environmental variables with the VIF computation and the ANOVA-like permutation tests applied to the individual terms confirmed the inclusion of all populations, glycogen content, perturbation 2, and the interaction between population and perturbation in general (for all cases, P < 0.001). Based on the triplot (Fig. 1B) constructed for the model including all these mentioned variables, predominantly negative correlations were apparent for the qPCR data set. However, certain positive associations were also identified. Specifically, Cardinium was closely associated with the laboratory population and its interactions, along with perturbation 1. The vector of glycogen status formed very small angles with the vectors of Actinobacteria and Blattabacterium-like bacteria (and to some extent also with the vector of Wolbachia). Perturbation 2 explained a certain degree of variance among the other constituents of the microbiome.

DISCUSSION

Differences in bacterial communities among the mite populations.Previous analyses of insect species possessing various eating habits, i.e., omnivores, herbivores, and predators, have revealed significant influences of both diet and taxonomic origin on the spectrum of associated bacteria (20–23). This study examined four different populations (strains) of one astigmatid mite species. Our experiments were based on the assumption that the associated bacteria are influenced by mite origin (population) or different diets (perturbation). The observation that population is associated with microbiome composition indicates that population-specific factors, including genotype and metacommunity composition, influence the microbiome. According to the results of the RDAs, the different T. putrescentiae populations exerted a greater influence than perturbation on microbiome variability. We hypothesize that this result is caused by high stability of the symbionts in response to dietary perturbation in short-term adaptation. Once the microbiome is established, the symbionts appear relatively resilient to subsequent perturbations, such as diet changes. We cannot exclude different situations in the microbiome after long-term diet adaptation when the mite physiology is changed (24). Although this factor is a selective process, we cannot rule out neutral processes, including passive dispersion (25), during perturbation.

The endosymbionts Wolbachia, Cardinium, Blattabacterium-like, and Solitalea-like bacteria contributed to the explained variability in both the Illumina and qPCR data sets. Although Cardinium and Blattabacterium-like bacteria were detected in all populations by Illumina technology, they were under the detection limits for qPCR analysis (1 copy per mite). Therefore, we considered Cardinium and Blattabacterium-like absent, similar to the findings of a previous study in which we used 16S rRNA Sanger sequences (15). The next identified core taxa included Bartonella-like, S. cohnii, B. cereus, Brevibacterium siliguriense, and Enterobacteriaceae with low similarity to Xenorhabdus, which were all resistant to perturbation. Here, we also provide evidence that the mites acquire bacteria (Lactobacillus, Pantoea, Prevotella, and Virgibacillus) after perturbation, as previously observed in cultivation experiments (17, 18, 26). This is also supported by the beta diversity data according to the Jaccard index (Tables 1 and 2). The T. putrescentiae populations also differed in terms of the proportions of B. cereus, S. cohnii, Virgibacillus, Brevibacterium siliguriense, Micrococcaceae, and Enterobacteriaceae taxa. Analysis of a wide spectrum of insect species revealed the presence of Firmicutes as the main group of bacteria found in insect guts (23), which was confirmed in this study in the astigmatid mite T. putrescentiae. The presence of Enterobacteriaceae taxa in T. putrescentiae raises the question of whether these bacterial taxa are symbiotic in the mite. Previously, we observed sequences with 99% similarity to Xenorhabduscabanillasii among clones from L. destructor and A. siro (10), whereas sequences similar to those of Erwinia bacteria comprised the majority of the microbiome of Blomia tropicalis (9). “Candidatus Erwinia dacicola” resides intracellularly in the gastric ceca of the larval midgut but extracellularly in the lumen of the foregut and the ovipositor diverticulum in adult flies (27). Moran et al. (28) described three symbiotic Enterobacteriaceae species, “Candidatus Serratia symbiotica,” “Candidatus Hamiltonella defensa,” and “Candidatus Regiella insecticola,” as sister groups to one another; moreover, they revealed a relationship to a species of Photorhabdus. These bacterial species occur in aphids and are localized to bacteriocytes, sheath cells, and the hemocoel (28). Based on these findings in the literature, the bacteria observed in this study may be located in both the gut and the parenchymal or reproductive tissues of the mite; future localization studies are necessary to confirm this assumption. The Enterobacteriaceae have been suggested to be beneficial based on their contribution to the longevity of the Mediterranean fruit fly (Ceratitis capitata) (29). The sequences obtained in the current analyses differed from those identified previously in the stored product mite Xenorhabdus (10). The Enterobacteriaceae sequences with low similarity (92%) to those of Xenorhabdus from T. putrescentiae should be identified in future studies using markers other than the 16S rRNA.

Differences in bacterial communities caused by perturbation and the effect on population growth.In insects, different nutritive factors, such as the protein-carbohydrate balance, influence the animal's life history, including reproduction (30, 31). Perturbation can lead to dysbiosis, i.e., a disruption of the proportions of gut core bacteria, and homeostasis, i.e., changes in the host physiology (32). In our experimental design, we used two initial (rearing) diets, SPMd for the majority of the T. putrescentiae populations and DDFd for the dog population, which might influence the results obtained. The perturbation factor overcame this issue, so the statistical analyses were not influenced by the difference in the initial diets.

In this study, we used population growth and increased nutrient contents in the mite samples as indirect indicators of fitness. Under suitable conditions, we assumed shorter development and a higher number of offspring than under unsuitable conditions, resulting in increased growth (33). All selected diets used in this study are natural for T. putrescentiae (34). The connection between diet perturbation and changes in the gut bacterial community is well documented (21, 23, 35). We observed that diet perturbation was associated with a rapid change in T. putrescentiae population growth, as previously reported (24), and we confirmed the strong effect of perturbation on both population growth and nutrient status. The effect of diet perturbation on the variability of the microbiome composition and population growth was lower than the effect induced by the mite population. In contrast, perturbation had a greater effect on nutritional status in mites than the mite population. Experiments with bacterial communities in Drosophila melanogaster have revealed that the bacterial community influences the relationship between energy storage and body mass and coordinates metabolism with body size (36). The RDA statistical model revealed the influence of lipid and glycogen status on the microbiome, which is correlated with Actinobacteria and Blattabacterium-like bacteria.

All populations contained Bartonella-like, Bacillus, Staphylococcus, Brenneria, or Kocuria bacteria. Perturbation resulted in changes in the numbers of these five taxa and the subsequent acquisition or disappearance of three or four taxa (Lactobacillus, Pantoea, and Virgibacillus). The Koppert population was associated with high numbers of Bacillaceae with similarity (90 to 99%) to Virgibacillus. The Sanger sequencing of the 16S rRNA previously identified sequences with similarity (94 to 99%) to Virgibacillus after a diet shift from SPMd to a fungal diet (37). These Bacillaceae include bacteria that produce enzymes facilitating the utilization of nutrients by the mites (18). Bacterial exoenzymes are active against difficult-to-digest structural polysaccharides such as cellulose and chitin (18) and are increased in the bacterial communities of T. putrescentiae on Fusarium diets (37). Bacillus cereus isolates originating from dog kernels infested by T. putrescentiae produce alkaline and neutral proteases (38). The participation of mite-associated bacteria in the digestion of skin, hair, nails, and feathers has been suggested (39) and could enhance the diet switch for T. putrescentiae.

A recent study by Serbus et al. (40) showed that the indirect effects of diet on Wolbachia titers in Drosophila ovaria include influencing cyclic AMP (cAMP)-regulated transcriptional coactivators in the TORC pathway (40). Experiments with the white fly Bemisia tabaci and the symbiont Hamiltonella showed that for insects reared on a standard-nitrogen diet, no cost or benefit was associated with Hamiltonella infection, whereas on low-nitrogen diets, Hamiltonella-infected whiteflies often grew better than uninfected whiteflies (41). In this study, diet perturbation influenced endosymbiotic bacteria, i.e., Blattabacterium-like and Wolbachia bacteria, whereas Solitalea-like and Cardinium bacteria were influenced marginally. However, the effect of perturbation on the absolute numbers of endosymbionts in the T. putrescentiae populations appeared to be random rather than systematic.

Multi-infection by Cardinium and Wolbachia is common at the individual level in spider mites (Acari: Tetranychidae) (42). At the population level, we observed different OTUs for Wolbachia in one population, in which low-abundance OTUs appeared after perturbation. In addition, the T. putrescentiae populations in our study were multi-infected with more taxa of endosymbiotic bacteria. The infection frequencies in the mite populations ranged from completely uninfected to a polymorphic state (43). Tetranychus truncates mites doubly infected with Cardinium and Wolbachia exhibit greater fecundity than uninfected individuals (44). Based on our data, we are unable to conclude whether perturbation causes dysbiosis in terms of the changes in the proportions of specific bacterial taxa in T. putrescentiae populations. Both the microbiome and population growth were more influenced by population than by diet perturbation, indicating the effect of endosymbionts or other bacteria on population growth. We hypothesize that the different proportions of Cardinium/Wolbachia or Blattabacterium-like/Wolbachia bacteria at the population level after perturbation could result in changes in population growth due to the various combinations of endosymbiotic bacteria in T. putrescentiae individuals. Similar changes have been observed in insects (2, 4) although the effects of these symbionts could also be the results of neutral processes (25). However, additional experiments are necessary to compare the effects of endosymbiotic bacteria on T. putrescentiae population growth.

MATERIALS AND METHODS

Mites and rearing diets.The bacterial communities of four populations of Tyrophagus putrescentiae (Schrank, 1781) were compared: (i) the laboratory population, which was sampled from a grain store in Bustehrad, Czechia, by E. Zdarkova in 1996; (ii) the dog population, which was sampled from dog food in the United States by J. Hubert in 2007; (iii) the Koppert population, which was obtained from a grain-derived diet by E. van Baal in the Netherlands in 2012; and (iv) the Phillips population, which was obtained from dog food by T. Phillips in Manhattan, KS, USA, in 2014. The mites were mass reared in Iwaki flasks and collected as previously described (45).

Diets utilized for perturbation.Except for the mites of the dog population, T. putrescentiae was reared on a wheat-derived diet designed for stored product mites (SPMd), which included a mixture of oat flakes, wheat germ, and Pangamin-dried yeast extract (Rapeto, Bezdruzice, Czechia) (10:10:1 wt/wt). The mites were transferred to the diets up to 6 months after collection and maintained until the experiments.

The dog population was maintained on a ground dry dog food kernel diet (DDFd) (Friskies Junior Life Plus Nutrition; Nestle Purina, Buk, Hungary) (see references 24 and 14). Both diets were heated at 70°C for 30 min and rehydrated. SPMd, DDFd, and the oilseed rape (Brassica napus L.) diet (OSRd) were compared in the dietary switch experiments. SPMd and DDFd are described in the previous section. The oilseeds were derived from a mixture of cultivars and were prepared according to a previously described protocol (46).

Treatments and samples.The experimental design included a rearing diet, SPMd, and two experimental diets, OSRd and DDFd. However, for the dog population, DDFd was the rearing diet, and OSRd and SPMd were the experimental diets. For the statistical analyses, the factor variable “perturbation” was coded in the following manner: 0, rearing diet, i.e., SPMd for all populations except the dog population (DDFd); 1, SPMd/DDFd for all populations or DDFd/SPMd for the dog population; 2, SPMd/OSRd for all populations or DDFd/OSRd for the dog population.

Unsexed mites (ca. 1,000 individuals) from approximately 1-month-old cultures were used in all experiments. Flasks containing mites and food were placed in Secador desiccator cabinets (Bel-Art Products, Wayne, NJ, USA) at 85% relative humidity and 25°C. The design included at least 10 chambers per diet and mite population. The mites that remained on the plug or the surfaces of the chambers were collected using a fine-tipped artist's paint brush. Then, the samples of mite populations were transferred into Eppendorf tubes and weighed on a microbalance (MS Mettler-Toledo, Greifensee, Switzerland); each sample had an accurate weight ranging from 0.015 to 0.025 g of wet weight. Multiple mites from each population were used for DNA extraction and the evaluation of protein, saccharide, lipid, and glycogen contents. DNA was extracted from the surface-cleaned and homogenized mites according to a previously described protocol (15). The extracted DNA was stored in a freezer at −28°C until analysis.

Growth test.The growth test began by switching 10 adult mites onto the two experimental diets and one rearing diet in Iwaki flasks. We performed 10 replicates per diet and mite population. The mites were maintained under the conditions described above. The experiment was terminated as described previously (38), and the mites were counted under a dissection microscope. To standardize the diverse levels of individual numbers associated with the diets, the vectors of the experimental diets were projected onto the vector of the rearing diet (SPMd or, in the case of the dog population, DDFd). The well-known procedure of vector projection is detailed, for instance, by Gentle (47). Then, a redundancy analysis (RDA) (48, 49) was applied to examine whether the perturbation exerted any influence or if the changes in growth were determined by populations. Using the vegan package (50) of the R statistical software (51), models with diets and perturbation were built. In addition, the interactions between population and perturbation were investigated. In particular, we aimed at changes possibly caused by perturbation within separate populations, where a special RDA was employed as well. Permutation tests (1,000 replications) were employed to determine the effect of the explanatory variables.

Nutritional status: protein, saccharide, lipid, and glycogen contents in mites.Using a method modified from Kaufmann (52), the mites were homogenized by beating in a FastPrep-24 homogenizer (MP Biomedicals, Santa Ana, CA, USA) in 200 μl of 2% sodium sulfate. Twenty-five microliters of each homogenate was used to estimate protein contents using the Bradford reagent (catalog number B6916; Sigma-Aldrich, St. Louis, MO, USA). The optical density (650 nm) was measured on a Secomam spectrophotometer (Secomam, Ales, France) according to the manufacturer's instructions. The reagent was first calibrated to a protein standard (catalog number P0834; Sigma-Aldrich). The remainder of the homogenates was precipitated in 800 μl of chloroform-methanol mixed 1:1 (vol/vol) and centrifuged (3,000 × g for 60 s). The pellets were retained for glycogen analysis. The supernatants were transferred to new Eppendorf tubes, and 600 μl of distilled water was added. After centrifugation (3,000 × g for 60 s), the top fraction (water-methanol) was collected for sugar analysis, and the bottom portion (chloroform) was used for lipid analysis. For lipid analysis, chloroform was evaporated on a heat block (110°C), and 200 μl of sulfuric acid (95 to 98%) was added for 10 min at 110°C. Next, 1 ml of vanillin-phosphoric acid reagent was added. The reagent contained 0.6 g of vanillin in 100 ml of hot distilled water and 400 ml of 85% phosphoric acid. The developed color was measured on a Secomam spectrophotometer at 450 nm. The reaction was calibrated to a lipid standard. For the sugar and glycogen assays, the tubes were placed on a heat block (110°C) to evaporate the solvent to a total volume of 100 μl, and 1 ml of anthrone reagent was added. The anthrone reagent contained 150 ml of distilled water, to which 385 ml of sulfuric acid (95 to 98%) and 0.75 g of anthrone were added. The samples containing the anthrone solution were mixed and heated for 17 min, and the optical density was then measured at 625 nm. The reaction was calibrated using a glucose standard (catalog number G6918; Sigma-Aldrich).

The data were recalculated to milligrams/gram of fresh weight. Two types of analyses were performed using R statistical software (51) and the vegan package (50). First, RDA models for individual nutrients were built similarly to those for growth as a response variable. Again, the interactions between the explanatory variables and changes within separate populations were studied as well. Second, an RDA model based on the correlation matrix computed from all nutrients together was utilized. Finally, permutation tests (1,000 replicates) were employed to determine whether population and perturbation exhibited important effects on the variability in nutritional status.

Amplification and sequencing.The design included four populations on three diets and six biological replicates (i.e., 72 samples). DNA samples were sent to MR DNA (Shallowater, TX, USA) for sequencing of the V1-V3 portions of the 16S rRNA genes using the universal primers 27Fmod and 519Rmod on an Illumina MiSeq platform, employing methods based on the bTEFAP (bacterial tag-encoded FLX amplicon pyrosequencing) process (53) in MR DNA.

Bioinformatic analyses.The multiplex sequences were processed using the mothur, version 1.36.1, software (54), according to the MiSeq standard operation procedure (SOP) (55) and UPARSE pipeline in USEARCH software (56). The mothur commands, which were used with default settings, are available online (http://www.mothur.org/wiki/MiSeq_SOP ). Read lengths, both forward and reverse, were 300 bp. The sequences were demultiplexed using mothur. The list of samples and barcodes is provided in Table S2 in the supplemental material, and primers are given in Table S1. The sequences were assembled in UPARSE, and 6,927,385 sequences were obtained. The barcodes and primers were trimmed in mothur and fastq, and renamed fasta files were prepared. The renamed fasta file contained the sample name in the beginning of the first line. Then the both fastq and renamed fasta files were processed in UPARSE. Individual OTUs were constructed by binning sequences into clusters with 97% similarity, discarding putative chimeric OTUs (2.6%) and singletons (93.4% from unique sequences) in the process. The OTUs were classified according to Ribosomal Database Project training set 15 (57). Representative sequences for OTUs (from UPARSE) were processed using the blastn program on the NCBI platform (58) to identify closely related sequences in GenBank. We identified Cardinium, Solitalea-like, and Bartonella-like OTUs by comparison of the representative sequences to previously obtained, nearly full-length 16S rRNA cloned sequences from T. putrescentiae (15) (Table S3). The taxonomical features of the samples were visualized by performing KRONA projections (59).

The shared file was constructed and processed using mothur and later in PAST, version 3.06, software (60). The obtained sequences differed in the number of reads, which ranged from 25,071 to 152,889 (Fig. S1). The data standardization was based on the subsample created in mothur on 25,071 sequences. The subsampled data set (Table S4) was used for all further statistical analyses with Illumina data. The inverse Simpson index was calculated in mothur and then subjected to another RDA very similar to the previous one (see “Growth test” and the paragraph above on nutritional status). For Venn diagrams, OTUs were selected based on the criterion of presence in all replicates per treatment. Based on the selection, the Venn diagrams were constructed per population or treatment in mothur. The effects of population and diet on OTU distribution were tested by performing two-way and one-way permutational multivariate analysis of variance (PERMANOVA) (61) in the vegan package of the R platform, function adonis (50), using the Bray-Curtis and Jaccard distance matrices with 1,000 permutations. In the two-way analysis, we explored the interactions not only between population and diets (or perturbation) but also between population and growth and individual nutrients. In the one-way analysis, we studied how the diet or its switch influenced OTUs within every single population. All computations were performed using the vegan package in R (50). Furthermore, using the same package, several db-RDA models (occasionally referred to as models of the constrained analysis of principal coordinates) based on the Bray-Curtis or Jaccard dissimilarity matrices were constructed that included population and diet (eventually perturbation) as the explanatory variables. Fitness variables (growth and nutrients) were also added to explore whether they contribute to model improvement in terms of decreasing the P value. We simultaneously controlled the values of variance inflation factors (VIFs) (62) to avoid collinearity in the models. The R package packfor (63), which can indicate significant explanatory variables via the forward selection procedure, was used. To visualize the coordinates resulting from the best db-RDA models, triplots were created. Furthermore, a heat map depicting OTU abundance and clustering in dendrograms was produced using the gplots package in R (64). The OTU abundance data were logarithmically transformed as suggested by Anderson et al. (65). We prefer the log2 transformation, which permits better visualization of less abundant species. When constructing the dendrograms, we used the Bray-Curtis transformation. The similarity profile (SIMPROF) procedure (66) was applied to determine whether sample clusters were significant using the clustsig package in R (67).

A random forest algorithm with logarithmically transformed (log2) data and METASTATS (68) were applied to describe the effects of diets on OTUs for each population separately. The analyses were performed in mothur (1,500 trees and 1,000 replicates, respectively).

qPCR analyses.Amplifications were performed on a StepOnePlus real-time PCR system (Life Technologies, Grand Island, NY, USA) in 96-well plates using GoTaq qPCR master mix (Promega). SYBR green (Bio-Rad Laboratories, The Netherlands) was used as a double-stranded DNA (dsDNA) binding dye. The protocol and standard preparation were as described previously (16), and the primers and annealing temperatures are described in Table S1. The design employed six biological replicates, each analyzed with two technical repeats. To evaluate the effects of population and perturbation, the obtained numbers of copies were first logarithmically transformed (log10). Since the data were already transformed this way, we built on the Euclidean indices in subsequent analyses using the dissimilarity or distance matrices. Again, two-way and one-way PERMANOVA were performed together with db-RDAs. Whereas the two-way PERMANOVA technique provided specific insights regarding the relationships between interactions (population and perturbation), its one-way counterpart permitted the study of the influence of diets (or their switches) within each population separately. Using selected db-RDA models and indicating significance according to their P values obtained via the permutation tests (1,000 replicates), we mainly constructed triplots and visualized the situation. Experimentally, other explanatory variables such as growth and nutrients were added to the RDA models to determine whether they contribute to improvements in the P values.

Accession number(s).Sequences were deposited in the NCBI database under BioProject number PRJNA328005 and SRA number SRP078077 .

ACKNOWLEDGMENTS

The authors are obliged to the referees for their highly useful comments on the previous version of the manuscript. We are grateful to Elmer van Baal and Tom Phillips for providing the T. putrescentiae populations, to Jan Kopecky for primer design and his advice and help with sequence processing, and to Martin Markovic for technical help.

This study was supported by project number GA15-09038S of the Czech Science Foundation.

FOOTNOTES

    • Received 17 January 2017.
    • Accepted 21 February 2017.
    • Accepted manuscript posted online 24 February 2017.
  • Supplemental material for this article may be found at https://doi.org/10.1128/AEM.00128-17 .

  • Copyright © 2017 American Society for Microbiology.

All Rights Reserved .

REFERENCES

  1. 1.↵
    1. Dale C,
    2. Moran NA
    . 2006. Molecular interactions between bacterial symbionts and their hosts. Cell126:453–465. doi:10.1016/j.cell.2006.07.014.
    OpenUrlCrossRefPubMedWeb of Science
  2. 2.↵
    1. Feldhaar H
    . 2011. Bacterial symbionts as mediators of ecologically important traits of insect hosts. Ecol Entomol36:533–543. doi:10.1111/j.1365-2311.2011.01318.x.
    OpenUrlCrossRef
  3. 3.↵
    1. McLean AHC,
    2. Parker BJ,
    3. Hrcek J,
    4. Henry LM,
    5. Godfray HCJ
    . 2016. Insect symbionts in food webs. Philos Trans R Soc Lond B Biol Sci371:20150325. doi:10.1098/rstb.2015.0325.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Ferrari J,
    2. Vavre F
    . 2011. Bacterial symbionts in insects or the story of communities affecting communities. Philos Trans R Soc Lond B Biol Sci366:1389–1400. doi:10.1098/rstb.2010.0226.
    OpenUrlCrossRefPubMed
  5. 5.↵
    1. Weinert LA,
    2. Araujo-Jnr EV,
    3. Ahmed MZ,
    4. Welch JJ
    . 2015. The incidence of bacterial endosymbionts in terrestrial arthropods. Proc Biol Sci282:20150249. doi:10.1098/rspb.2015.0249.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Werren JH,
    2. Baldo L,
    3. Clark ME
    . 2008. Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol6:741–751. doi:10.1038/nrmicro1969.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    1. Ros VID,
    2. Fleming VM,
    3. Feil EJ,
    4. Breeuwer JAJ
    . 2012. Diversity and recombination in Wolbachia and Cardinium from Bryobia spider mites. BMC Microbiol12:S13. doi:10.1186/1471-2180-12-S1-S13.
    OpenUrlCrossRefPubMed
  8. 8.↵
    1. Brown AN,
    2. Lloyd VK
    . 2015. Evidence for horizontal transfer of Wolbachia by a Drosophila mite. Exp Appl Acarol66:301–311. doi:10.1007/s10493-015-9918-z.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Hubert J,
    2. Kopecky J,
    3. Sagova-Mareckova M,
    4. Nesvorna M,
    5. Zurek L,
    6. Erban T
    . 2016. Assessment of bacterial communities in thirteen species of laboratory-cultured domestic mites (Acari: Acaridida). J Econ Entomol109:1887–1896. doi:10.1093/jee/tow089.
    OpenUrlCrossRefPubMed
  10. 10.↵
    1. Hubert J,
    2. Kopecky J,
    3. Perotti MA,
    4. Nesvorna M,
    5. Braig HR,
    6. Sagova-Mareckova M,
    7. Macovei L,
    8. Zurek L
    . 2012. Detection and identification of species-specific bacteria associated with synanthropic mites. Microb Ecol63:919–928. doi:10.1007/s00248-011-9969-6.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Bowman CE
    . 1984. Comparative enzymology of economically important astigmatid mites, p 993–1001. InGriffiths DA, Bowman CE (ed), Acarology VI, vol 2. Ellis Horwood, Chichester, United Kingdom.
    OpenUrl
  12. 12.↵
    1. Childs M,
    2. Bowman CE
    . 1981. Lysozyme activity in six species of economically important astigmatid mites. Comp Biochem Physiol B Comp Biochem70:615–617. doi:10.1016/0305-0491(81)90305-9.
    OpenUrlCrossRefWeb of Science
  13. 13.↵
    1. Bowman CE
    . 1981. Hide protease in stored product mites (Astigmata: Acaridae). Comp Biochem Physiol B Comp Biochem70:803–805. doi:10.1016/0305-0491(81)90022-5.
    OpenUrlCrossRef
  14. 14.↵
    1. Rybanska D,
    2. Hubert J,
    3. Markovic M,
    4. Erban T
    . 2016. Dry dog food integrity and mite strain influence the density-dependent growth of the stored-product mite Tyrophagus putrescentiae (Acari: Acaridida). J Econ Entomol109:454–460. doi:10.1093/jee/tov298.
    OpenUrlCrossRef
  15. 15.↵
    1. Erban T,
    2. Klimov PB,
    3. Smrz J,
    4. Phillips TW,
    5. Nesvorna M,
    6. Kopecky J,
    7. Hubert J
    . 2016. Populations of stored product mite Tyrophagus putrescentiae differ in their bacterial communities. Front Microbiol7:1046. doi:10.3389/fmicb.2016.01046.
    OpenUrlCrossRef
  16. 16.↵
    1. Hubert J,
    2. Kopecky J,
    3. Nesvorna M,
    4. Perotti MA,
    5. Erban T
    . 2016. Detection and localization of Solitalea-like and Cardinium bacteria in three Acarus siro populations (Astigmata: Acaridae). Exp Appl Acarol70:309–327. doi:10.1007/s10493-016-0080-z.
    OpenUrlCrossRef
  17. 17.↵
    1. Smrz J
    . 2003. Microanatomical and biological aspects of bacterial associations in Tyrophagus putrescentiae (Acari: Acaridida). Exp Appl Acarol31:105–113. doi:10.1023/B:APPA.0000005111.05959.d6.
    OpenUrlCrossRefPubMed
  18. 18.↵
    1. Smrz J,
    2. Catska V
    . 2010. Mycophagous mites and their internal associated bacteria cooperate to digest chitin in soil. Symbiosis52:33–40. doi:10.1007/s13199-010-0099-6.
    OpenUrlCrossRef
  19. 19.↵
    1. Smrz J
    . 2006. Types of haemocytes in saprophagous soil mites (Acari: Oribatida, Acaridida), and the correlation between their presence and certain processes within mites. Eur J Entomol103:679–686. doi:10.14411/eje.2006.088.
    OpenUrlCrossRef
  20. 20.↵
    1. Broderick NA,
    2. Raffa KF,
    3. Goodman RM,
    4. Handelsman J
    . 2004. Census of the bacterial community of the gypsy moth larval midgut by using culturing and culture-independent methods. Appl Environ Microbiol70:293–300. doi:10.1128/AEM.70.1.293-300.2004.
    OpenUrlAbstract/FREE Full Text
  21. 21.↵
    1. Robinson CJ,
    2. Schloss P,
    3. Ramos Y,
    4. Raffa K,
    5. Handelsman J
    . 2010. Robustness of the bacterial community in the cabbage white butterfly larval midgut. Microb Ecol59:199–211. doi:10.1007/s00248-009-9595-8.
    OpenUrlCrossRefPubMedWeb of Science
  22. 22.↵
    1. Colman DR,
    2. Toolson EC,
    3. Takacs-Vesbach CD
    . 2012. Do diet and taxonomy influence insect gut bacterial communities?Mol Ecol21:5124–5137. doi:10.1111/j.1365-294X.2012.05752.x.
    OpenUrlCrossRefPubMedWeb of Science
  23. 23.↵
    1. Yun J-H,
    2. Roh SW,
    3. Whon TW,
    4. Jung M-J,
    5. Kim M-S,
    6. Park D-S,
    7. Yoon C,
    8. Nam Y-D,
    9. Kim Y-J,
    10. Choi J-H,
    11. Kim J-Y,
    12. Shin N-R,
    13. Kim S-H,
    14. Lee W-J,
    15. Bae J-W
    . 2014. Insect gut bacterial diversity determined by environmental habitat, diet, developmental stage, and phylogeny of host. Appl Environ Microbiol80:5254–5264. doi:10.1128/AEM.01226-14.
    OpenUrlAbstract/FREE Full Text
  24. 24.↵
    1. Erban T,
    2. Rybanska D,
    3. Hubert J
    . 2015. Population growth of the generalist mite Tyrophagus putrescentiae (Acari: Acaridida) following adaptation to high- or low-fat and high- or low-protein diets and the effect of dietary switch. Environ Entomol44:1599–1604. doi:10.1093/ee/nvv129.
    OpenUrlCrossRefPubMed
  25. 25.↵
    1. Burns AR,
    2. Stephens WZ,
    3. Stagaman K,
    4. Wong S,
    5. Rawls JF,
    6. Guillemin K,
    7. Bohannan BJM
    . 2016. Contribution of neutral processes to the assembly of gut microbial communities in the zebrafish over host development. ISME J10:655–664. doi:10.1038/ismej.2015.142.
    OpenUrlCrossRefPubMed
  26. 26.↵
    1. Smrz J,
    2. Soukalova H,
    3. Catska V,
    4. Hubert J
    . 2016. Feeding patterns of Tyrophagus putrescentiae (Sarcoptiformes: Acaridae) indicate that mycophagy is not a single and homogeneous category of nutritional biology. J Insect Sci16:94. doi:10.1093/jisesa/iew070.
    OpenUrlCrossRefPubMed
  27. 27.↵
    1. Estes AM,
    2. Hearn DJ,
    3. Bronstein JL,
    4. Pierson EA
    . 2009. The olive fly endosymbiont, “Candidatus Erwinia dacicola,” switches from an intracellular existence to an extracellular existence during host insect development. Appl Environ Microbiol75:7097–7106. doi:10.1128/AEM.00778-09.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Moran NA,
    2. Russell JA,
    3. Koga R,
    4. Fukatsu T
    . 2005. Evolutionary relationships of three new species of Enterobacteriaceae living as symbionts of aphids and other insects. Appl Environ Microbiol71:3302–3310. doi:10.1128/AEM.71.6.3302-3310.2005.
    OpenUrlAbstract/FREE Full Text
  29. 29.↵
    1. Behar A,
    2. Yuval B,
    3. Jurkevitch E
    . 2008. Gut bacterial communities in the Mediterranean fruit fly (Ceratitis capitata) and their impact on host longevity. J Insect Physiol54:1377–1383. doi:10.1016/j.jinsphys.2008.07.011.
    OpenUrlCrossRefPubMedWeb of Science
  30. 30.↵
    1. Clark RM,
    2. Zera AJ,
    3. Behmer ST
    . 2015. Nutritional physiology of life-history trade-offs: how food protein-carbohydrate content influences life-history traits in the wing-polymorphic cricket Gryllus firmus. J Exp Biol218:298–308. doi:10.1242/jeb.112888.
    OpenUrlAbstract/FREE Full Text
  31. 31.↵
    1. Lee KP
    . 2015. Dietary protein:carbohydrate balance is a critical modulator of lifespan and reproduction in Drosophila melanogaster: a test using a chemically defined diet. J Insect Physiol75:12–19. doi:10.1016/j.jinsphys.2015.02.007.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Schwarz RS,
    2. Moran NA,
    3. Evans JD
    . 2016. Early gut colonizers shape parasite susceptibility and microbiota composition in honey bee workers. Proc Natl Acad Sci U S A113:9345–9350. doi:10.1073/pnas.1606631113.
    OpenUrlAbstract/FREE Full Text
  33. 33.↵
    1. Matsumoto K
    . 1965. Studies on environmental factors for breeding of grain mites VII. Relationship between reproduction of mites and kind of carbohydrates in the diet. Jpn J Sanit Zool16:118–122. (In Japanese.)
    OpenUrl
  34. 34.↵
    1. Hughes AM
    . 1976. The mites of stored food and houses, 2nd ed. Technical bulletin 9. Ministry of Agriculture, Fisheries and Food, London, United Kingdom.
  35. 35.↵
    1. Santo Domingo JW,
    2. Kaufman MG,
    3. Klug MJ,
    4. Holben WE,
    5. Harris D,
    6. Tiedje JM
    . 1998. Influence of diet on the structure and function of the bacterial hindgut community of crickets. Mol Ecol7:761–767. doi:10.1046/j.1365-294x.1998.00390.x.
    OpenUrlCrossRefWeb of Science
  36. 36.↵
    1. Wong A-C,
    2. Dobson AJ,
    3. Douglas AE
    . 2014. Gut microbiota dictates the metabolic response of Drosophila to diet. J Exp Biol217:1894–1901. doi:10.1242/jeb.101725.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Hubert J,
    2. Nesvorna M,
    3. Sagova-Mareckova M,
    4. Kopecky J
    . 2012. Shift of bacterial community in synanthropic mite Tyrophagus putrescentiae induced by Fusarium fungal diet. PLoS One7:e48429. doi:10.1371/journal.pone.0048429.
    OpenUrlCrossRefPubMed
  38. 38.↵
    1. Erban T,
    2. Rybanska D,
    3. Harant K,
    4. Hortova B,
    5. Hubert J
    . 2016. Feces derived allergens of Tyrophagus putrescentiae reared on dried dog food and evidence of the strong nutritional interaction between the mite and Bacillus cereus producing protease bacillolysins and exo-chitinases. Front Physiol7:53. doi:10.3389/fphys.2016.00053.
    OpenUrlCrossRefPubMed
  39. 39.↵
    1. Erban T,
    2. Hubert J
    . 2010. Comparative analyses of proteolytic activities in seven species of synanthropic acaridid mites. Arch Insect Biochem Physiol75:187–206. doi:10.1002/arch.20388.
    OpenUrlCrossRefPubMed
  40. 40.↵
    1. Serbus LR,
    2. White PM,
    3. Silva JP,
    4. Rabe A,
    5. Teixeira L,
    6. Albertson R,
    7. Sullivan W
    . 2015. The impact of host diet on Wolbachia titer in Drosophila. PLoS Pathog11:e1004777. doi:10.1371/journal.ppat.1004777.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Su Q,
    2. Xie W,
    3. Wang S,
    4. Wu Q,
    5. Liu B,
    6. Fang Y,
    7. Xu B,
    8. Zhang Y
    . 2014. The endosymbiont Hamiltonella increases the growth rate of its host Bemisia tabaci during periods of nutritional stress. PLoS One9:e89002. doi:10.1371/journal.pone.0089002.
    OpenUrlCrossRef
  42. 42.↵
    1. Zhang Y-K,
    2. Chen Y-T,
    3. Yang K,
    4. Qiao G-X,
    5. Hong X-Y
    . 2016. Screening of spider mites (Acari: Tetranychidae) for reproductive endosymbionts reveals links between co-infection and evolutionary history. Sci Rep6:27900. doi:10.1038/srep27900.
    OpenUrlCrossRef
  43. 43.↵
    1. Ros VID,
    2. Breeuwer JAJ
    . 2009. The effects of, and interactions between, Cardinium and Wolbachia in the doubly infected spider mite Bryobia sarothamni. Heredity102:413–422. doi:10.1038/hdy.2009.4.
    OpenUrlCrossRefPubMed
  44. 44.↵
    1. Zhao D-X,
    2. Zhang X-F,
    3. Hong X-Y
    . 2013. Host-symbiont interactions in spider mite Tetranychus truncates doubly infected with Wolbachia and Cardinium. Environ Entomol42:445–452. doi:10.1603/EN12354.
    OpenUrlCrossRefPubMed
  45. 45.↵
    1. Erban T,
    2. Hubert J
    . 2008. Digestive function of lysozyme in synanthropic acaridid mites enables utilization of bacteria as a food source. Exp Appl Acarol44:199–212. doi:10.1007/s10493-008-9138-x.
    OpenUrlCrossRefPubMed
  46. 46.↵
    1. Stara J,
    2. Stejskal V,
    3. Nesvorna M,
    4. Plachy J,
    5. Hubert J
    . 2011. Efficacy of selected pesticides against synanthropic mites under laboratory assay. Pest Manag Sci67:446–457. doi:10.1002/ps.2083.
    OpenUrlCrossRefPubMedWeb of Science
  47. 47.↵
    1. Gentle JE
    . 2007. Matrix algebra: theory, computations, and applications in statistics. Springer, New York, NY.
  48. 48.↵
    1. Borcard D,
    2. Gillet F,
    3. Legendre P
    . 2011. Numerical ecology with R, 1st ed.Springer, New York, NY.
  49. 49.↵
    1. Legendre P,
    2. Legendre L
    . 2012. Numerical ecology, 3rd ed.Elsevier, Amsterdam, the Netherlands.
  50. 50.↵
    1. Oksanen J,
    2. Blanchet FG,
    3. Kindt R,
    4. Legendre P,
    5. Minchin PR,
    6. O'Hara RB,
    7. Simpson GL,
    8. Solymos P,
    9. Stevens MHH,
    10. Wagner H
    . 2016. vegan: community ecology package. R Foundation for Statistical Computing, Vienna, Austria. https://CRAN.R-project.org/package=vegan.
  51. 51.↵
    R Core Team. 2016. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org.
  52. 52.↵
    1. Kaufmann C
    . 2014. Determination of lipid, glycogen and sugars in mosquitoes. InBenedict M (ed), MR4 methods in Anopheles research, 4th ed. BEI Resources, Manassas, VA.
  53. 53.↵
    1. Chiodini RJ,
    2. Dowd SE,
    3. Chamberlin WM,
    4. Galandiuk S,
    5. Davis B,
    6. Glassing A
    . 2015. Microbial population differentials between mucosal and submucosal intestinal tissues in advanced Crohn's disease of the ileum. PLoS One10:e0134382. doi:10.1371/journal.pone.0134382.
    OpenUrlCrossRef
  54. 54.↵
    1. Schloss PD,
    2. Westcott SL,
    3. Ryabin T,
    4. Hall JR,
    5. Hartmann M,
    6. Hollister EB,
    7. Lesniewski RA,
    8. Oakley BB,
    9. Parks DH,
    10. Robinson CJ,
    11. Sahl JW,
    12. Stres B,
    13. Thallinger GG,
    14. Van Horn DJ,
    15. Weber CF
    . 2009. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol75:7537–7541. doi:10.1128/AEM.01541-09.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Kozich JJ,
    2. Westcott SL,
    3. Baxter NT,
    4. Highlander SK,
    5. Schloss PD
    . 2013. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol79:5112–5120. doi:10.1128/AEM.01043-13.
    OpenUrlAbstract/FREE Full Text
  56. 56.↵
    1. Edgar RC
    . 2013. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods10:996–998. doi:10.1038/nmeth.2604.
    OpenUrlCrossRefPubMedWeb of Science
  57. 57.↵
    1. Cole JR,
    2. Wang Q,
    3. Fish JA,
    4. Chai B,
    5. McGarrell DM,
    6. Sun Y,
    7. Brown CT,
    8. Porras-Alfaro A,
    9. Kuske CR,
    10. Tiedje JM
    . 2014. Ribosomal Database Project: data and tools for high throughput rRNA analysis. Nucleic Acids Res42:D633–D642. doi:10.1093/nar/gkt1244.
    OpenUrlCrossRefPubMedWeb of Science
  58. 58.↵
    1. Altschul SF,
    2. Madden TL,
    3. Schaffer AA,
    4. Zhang J,
    5. Zhang Z,
    6. Miller W,
    7. Lipman DJ
    . 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res25:3389–3402. doi:10.1093/nar/25.17.3389.
    OpenUrlCrossRefPubMedWeb of Science
  59. 59.↵
    1. Ondov BD,
    2. Bergman NH,
    3. Phillippy AM
    . 2011. Interactive metagenomic visualization in a Web browser. BMC Bioinformatics12:385. doi:10.1186/1471-2105-12-385.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Hammer O,
    2. Harper DAT,
    3. Ryan PD
    . 2001. PAST: paleontological statistics software package for education and data analysis. Palaeontol Electron4:9–18.
    OpenUrl
  61. 61.↵
    1. Anderson MJ,
    2. Walsh DCI
    . 2013. PERMANOVA, ANOSIM, and the Mantel test in the face of heterogeneous dispersions: what null hypothesis are you testing?Ecol Monogr83:557–574. doi:10.1890/12-2010.1.
    OpenUrlCrossRef
  62. 62.↵
    1. Kutner MH,
    2. Nachtsheim CJ,
    3. Neter J,
    4. Li W
    . 2005. Applied linear statistical models, 5th ed.McGraw-Hill Irwin, Boston, MA.
  63. 63.↵
    1. Dray S,
    2. Legendre P,
    3. Blanchet G
    . 2013. packfor: Forward selection with permutation (Canoco p 46). R Project for Statistical Computing, Vienna, Austria.
  64. 64.↵
    1. Warnes GR,
    2. Bolker B,
    3. Bonebakker L,
    4. Gentleman R,
    5. Huber W,
    6. Liaw A,
    7. Lumley T,
    8. Maechler M,
    9. Magnusson A,
    10. Moeller S,
    11. Schwartz M,
    12. Venables B
    . 2016. gplots: various R programming tools for plotting data. R Project for Statistical Computing, Vienna, Austria. https://CRAN.R-project.org/package=gplots.
  65. 65.↵
    1. Anderson MJ,
    2. Ellingsen KE,
    3. McArdle BH
    . 2006. Multivariate dispersion as a measure of beta diversity. Ecol Lett9:683–693. doi:10.1111/j.1461-0248.2006.00926.x.
    OpenUrlCrossRefPubMedWeb of Science
  66. 66.↵
    1. Clarke KR,
    2. Somerfield PJ,
    3. Gorley RN
    . 2008. Testing of null hypotheses in exploratory community analyses: similarity profiles and biota-environment linkage. J Exp Mar Biol Ecol366:56–69. doi:10.1016/j.jembe.2008.07.009.
    OpenUrlCrossRefWeb of Science
  67. 67.↵
    1. Whitaker D,
    2. Christman M
    . 2014. clustsig: Significant Cluster Analysis. https://CRAN.R-project.org/package=clustsig . R Project for Statistical Computing, Vienna, Austria.
  68. 68.↵
    1. White JR,
    2. Nagarajan N,
    3. Pop M
    . 2009. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol5:e1000352. doi:10.1371/journal.pcbi.1000352.
    OpenUrlCrossRefPubMed
PreviousNext
Back to top
Download PDF
Citation Tools
Experimental Manipulation Shows a Greater Influence of Population than Dietary Perturbation on the Microbiome of Tyrophagus putrescentiae
Tomas Erban, Ondrej Ledvinka, Marta Nesvorna, Jan Hubert
Applied and Environmental Microbiology Apr 2017, 83 (9) e00128-17; DOI: 10.1128/AEM.00128-17

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print

Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this Applied and Environmental Microbiology article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Experimental Manipulation Shows a Greater Influence of Population than Dietary Perturbation on the Microbiome of Tyrophagus putrescentiae
(Your Name) has forwarded a page to you from Applied and Environmental Microbiology
(Your Name) thought you would be interested in this article in Applied and Environmental Microbiology.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Experimental Manipulation Shows a Greater Influence of Population than Dietary Perturbation on the Microbiome of Tyrophagus putrescentiae
Tomas Erban, Ondrej Ledvinka, Marta Nesvorna, Jan Hubert
Applied and Environmental Microbiology Apr 2017, 83 (9) e00128-17; DOI: 10.1128/AEM.00128-17
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

Acaridae
bacteria
DIET
Gastrointestinal Microbiome
Tyrophagus putrescentiae
bacteria
symbiont
feeding
fitness
Cardinium
Solitalea
Bartonella
Blattabacterium
Wolbachia
16S rRNA

Related Articles

Cited By...

About

  • About AEM
  • Editor in Chief
  • Editorial Board
  • Policies
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Article Types
  • Ethics
  • Contact Us

Follow #AppEnvMicro

@ASMicrobiology

       

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

 

American Society for Microbiology
1752 N St. NW
Washington, DC 20036
Phone: (202) 737-3600

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

 

Print ISSN: 0099-2240; Online ISSN: 1098-5336