Next Article in Journal
Differences in Outpatient Health Care Utilization 12 Months after COVID-19 Infection by Race/Ethnicity and Community Social Vulnerability
Next Article in Special Issue
Fast, Ungapped Reads Mapping Using Squid
Previous Article in Journal
Boredom Proneness Predicts Self-Assessed Decision Errors in Sports but Is Unrelated to Risk Taking in General
Previous Article in Special Issue
Insights into Asparaginase from Endophytic Fungus Lasiodiplodia theobromae: Purification, Characterization and Antileukemic Activity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Nanopore Sequencing-Based Evaluation of the Midgut Microbiota of the Summer Chafer (Amphimallon solstitiale L.) Associated with Possible Resistance to Entomopathogenic Nematodes

1
Laboratory of Biocontrol, Production, and Application of EPN, Centre for Interdisciplinary Research, The John Paul II Catholic University of Lublin, Konstantynów 1J, 20-708 Lublin, Poland
2
Department of Environmental Microbiology and Biotechnology, Nicolaus Copernicus University in Toruń, Lwowska 1, 87-100 Toruń, Poland
3
Department of Molecular Biology, Institute of Biological Sciences, The John Paul II Catholic University of Lublin, Konstantynów 1J, 20-708 Lublin, Poland
4
Department of Industrial and Environmental Microbiology, Institute of Biological Sciences, Maria Curie-Skłodowska University, Akademicka 19, 20-033 Lublin, Poland
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(6), 3480; https://doi.org/10.3390/ijerph19063480
Submission received: 1 February 2022 / Revised: 11 March 2022 / Accepted: 11 March 2022 / Published: 15 March 2022
(This article belongs to the Special Issue New Advances in Environmental Microbiology)

Abstract

:
Root-feeding Amphimallon solstitiale larvae and certain other scarab beetles are the main soil-dwelling pests found in Europe, while entomopathogenic nematodes (EPN) have been used as a biocontrol agent against these species. Our study provides the first detailed characterization of the bacterial community of the midgut in wild A. solstitiale larvae, based on the nanopore sequencing of the 16S rRNA gene. In the whole dataset, we detected 2586 different genera and 11,641 species, with only 83 diverse bacterial genera shared by all studied individuals, which may represent members of the core midgut microbiota of A. solstitiale larvae. Subsequently, we compared the midgut microbiota of EPN-resistant and T0 (prior to EPN exposure) individuals, hypothesizing that resistance to this parasitic infection may be linked to the altered gut community. Compared to the control, the resistant insect microbiota demonstrated lower Shannon and Evenness indices and significant differences in the community structure. Our studies confirmed that the gut microbiota alternation is associated with resistant insects; however, there are many processes involved that can affect the bacterial community. Further research on the role of gut microbiota in insect-parasitic nematode interaction may ultimately lead to the improvement of biological control strategies in insect pest management.

Graphical Abstract

1. Introduction

The larvae of some members of the Scarabaeidae family (Coleoptera), such as Amphimallon solstitiale, the summer chafer, feed on roots and other organic matter. This can be extremely destructive to agricultural and forest ecosystems during population outbreaks [1,2]. These pests (also called white grubs) live underground in a complex soil environment and are difficult to control using pesticides, which can also affect the environment. Therefore, the development of alternative management strategies has received considerable interest, with a particularly promising use of entomopathogenic nematodes (EPN) [3,4]. However, the success of EPN as biological control agents depends on many abiotic and biotic factors. Unsuccessful host infection by EPN may also be a result of the evading behaviors, morphological barriers, and physiological traits of the insect host [5,6]. Although considerable research on the innate immunity of insect hosts has been carried out [7], little is known about the contribution of the insect microbiome to the course of EPN infection. Intestinal bacteria can influence host resistance against pathogens by contributing to immune system development, providing a range of essential nutrients and protection from pathogen colonization by competing for nutrients or producing inhibitory substances [8,9]. The knowledge of the interactions between gut residents, the insect host, and pathogens has great potential for the improvement of methods for managing harmful insects. Of particular interest is the prediction of the efficiency of control methods based on the composition of the insect pest microbiota or a strategy targeting the host’s microbiota for insect pest control [10,11].
Nematodes from the genera Steinernema and Heterorhabditis have engaged in mutualistic symbiosis with entomopathogenic Xenorhabdus and Photorhabdus bacteria, respectively (family Morganellaceae). The bacteria inhabit the gut of cognate nematode infective juveniles (IJs)—this is the infective stage of these parasites, which live in soil, seeking out and invading susceptible insect larvae. After gaining the insect hemocoel, IJs release their mutualistic bacteria, which replicate fast and spread over the insect body, killing the host rapidly due to septicemia [12]. It is known that the larval midgut is a site of penetration for EPN in the process of infection. The midgut microbiota appears to interact with EPN. However, it has not been elucidated whether the antagonistic bacteria that have been detected might protect the grubs from EPN, shaping resistance to this parasite pressure [6,13]. Although the information on the intricate bacterial community housed by the gut of scarab beetle larvae is available, research has mainly focused on the bacteria that are harbored in the hindgut and their contribution to the food digestion process [14,15]. To date, there are no comprehensive studies on the midgut bacterial community of scarab larvae.
Our previous study showed a shift in the bacterial composition in the midgut microbiome, associated with the resistance of the larvae of the common cockchafer, Melolontha melolontha, to EPN [16]. To verify these findings, we focused on the microbiota of another common scarab pest—the larvae of the summer chafer, A. solstitiale, the microbiota of which has not yet been studied in detail. Hypothesizing that differences in the composition of the host microbiota may be the cause of higher resistance to EPN or whether the resistance trait can modify the gut microbiota, we characterized the midgut bacterial microbiota of resistant individuals and compared them to those that were T0 (prior to EPN exposure). To avoid confounding the effects of the application of antibiotics on host physiology, which is an alternative and more commonly used method of testing the contribution of gut microbiota to resistance against pathogens, we studied the bacterial microbiota in resistant insects following exposure to a high dose of EPN. Additionally, we compared the gut microbiota of A. solstitiale larvae among different developmental stages and sampling sites. The application of the nanopore sequencing of the 16S rRNA gene, providing many advantages such as simple sample preparation and relatively long read lengths, ensured gut community profiling and allowed us to avoid the limitations of culture-dependent methods. This work highlights the importance of bacteria inhabiting the insect gut in insect–nematode interaction.

2. Materials and Methods

2.1. Collection and Preparation of Larvae

Second and third instar (L2 and L3) A. solstitiale larvae were sampled from two localities: a forest tree nursery and an urban lawn in the Lublin region, Eastern Poland (Table 1 and Table S1 (Supplementary Materials)). The collection was performed in June 2019. The grubs were collected from the soil and placed separately in 50-mL plastic cups that were partly filled with soil taken from the sampling site. After transportation to the laboratory, the larvae were identified. In total, over 200 healthy A. solstitiale L2 and L3 larvae were selected for further study, following the protocol described previously in [16]. Shortly after, gut dissection was performed in 30 larvae, referred to as T0 (prior to EPN exposure). The extracted midguts were frozen and stored until processing. Simultaneously, the other larvae were placed in 1-L containers previously filled with soil and seeded with perennial ryegrass Lolium perenne L. (family Poaceae). The larvae were then exposed to Heterorhabditis megidis, Steinernema arenarium, Steinernema bicornutum, Steinernema carpocapse, and Steinernema kraussei at a dose of 2500 IJs per container. The experiment was carried out in controlled photoperiod and temperature conditions, according to the procedure described by Koppenhöfer and Fuzy (2004) [17]. We exposed 30 insect larvae to each EPN species (Table S1, Figure S1). After 14 days of EPN treatment, the larvae were examined and qualified as EPN-resistant in the case of healthy-looking and active specimens, in contrast to the EPN-sensitive larvae, which were dead and with the characteristic cadaver turgor, color, and smell. The survival ratio was calculated as the proportion of viable larvae to dead ones after exposure to each EPN species.

2.2. DNA Extraction and Nanopore High-Throughput Sequencing

Total bacterial DNA was extracted from the A. solstitiale midguts using a Bead-Beat Micro AX Gravity kit (A&A Biotechnology, Gdynia, Poland). The obtained DNA concentration was monitored with Qubit 2.0 (Invitrogen, Carlsbad, CA, USA) using a Qubit dsDNA HS Assay Kit (Thermo Fisher Scientific, Waltham, MA, USA). Extraction of microbial community DNA, nanopore sequencing, data processing, and taxonomic identification were performed by genXone (Złotniki, Poland). For accurate detection of the bacterial communities, we targeted the fragment of the 16S rRNA bacterial gene spanning V3-V8 hypervariable domains, which were amplified with the forward primer F338 5′-ACTCCTACGGGAGGCAGC-3′ and the reverse primer R1391 5′-GACGGGCGGTGTGTRCA-3′ [18,19]. PCR amplicon libraries were prepared using a Ligation Sequencing Kit 1D and sequenced on a GridION X5 sequencer (Oxford Nanopore Technologies, Oxford, UK).

2.3. Bioinformatic Analysis

After the sequencing run, the reads were filtered for their quality using the following criteria: an average read quality ≥ 10 Phred, an average quality of read fragments with a window size of 500 bp ≥ 6 Phred, and a minimum sequence length of 800 bp. Adapters and barcodes were removed using Porechop v.0.2.4. High-quality reads were processed for taxonomic identification by matching the NGS sequences with sequences deposited in the NCBI RefSeq for the 16S ribosomal RNA database, using a modified BLAST algorithm. For each sample, an equivalent of the OTU frequency table was built using an Excel spreadsheet. Sequencing data were deposited in NCBI under the accession number PRJNA 787447.

2.4. Exploratory Data Analyses

To assess the sequencing depth, rarefaction curves were generated for all 28 samples that were sequenced, using phyloseq (v. 1.34.0) in R (v. 4.0.3). Rarefying to an even depth (82,308 reads) was performed with the same software package. The alpha diversity metrics comprised the observed number of species (Sobs), Shannon–Wiener biodiversity index (H’), Simpson biodiversity index (1/D), and community Evenness (E). These were calculated using the vegan (v. 2.5.7) package in R. Differences in the alpha diversity metrics between the group of samples that were subjected to the nematode pressure and the T0 group (prior to the EPN exposure) were tested with Student’s t-test or the Mann–Whitney test using Past v 3.08 [20]. Principal component analysis (PCoA) and an analysis of similarities (ANOSIM) between groups were performed using the vegan package in R. Additionally, to identify the differences without losing a considerable number of sequences due to the normalization procedure, the fold-changes in bacterial genera were calculated using DESeq2 in R [21]. Microbiome networks in the same groups were built using igraph, qgraph, vegan, and MCL in R [22]. They were made using the method of sparse correlations for compositional data (SparCC) as proposed by Friedman and Alm (2012) [23]. The link-analysis method was used to detect hubs or keystone taxa in the networks [24]. The keystone genera were nodes in the network that had a significantly larger number of links compared to the other hubs in the network. Keystone taxa were sorted according to the “hubbiness” scores, and the five top taxa from each group were presented and discussed further. Clusters were identified with the use of the “walktrap” algorithm [25].

3. Results

3.1. Resistance of A. solstitiale Larvae to EPN Infection

A. solstitiale larvae (L2 and L3 developmental stages) were collected from soils in eastern Poland. Selected larvae were assigned to the T0 or the EPN-resistant groups. Thirteen individuals, which constituted the T0 group, were not exposed to EPN, while fifteen randomly selected larvae, which constituted the EPN-resistant group, were selected after EPN exposure in the laboratory, based on the lowest level of susceptibility to this parasite (Table 1). Our survival test revealed that 52%, 68%, 66%, 57%, and 68% of the A. solstitiale larvae were healthy after 14 days of exposure to S. arenarium, S. bicornutum, S. carpocapse, S. kraussei, and H. megidis, respectively.

3.2. Nanopore Sequencing Results

The fourth-generation sequencing of the 16S rRNA gene using nanopore technology was applied to characterize the midgut bacterial community. The bioinformatics data processing resulted in 5,404,166 high-quality reads. The number of bacterial reads obtained across all the samples ranged from 82,317 to 646,404 (median 177,433). The rarefaction curves indicated a high coverage of bacterial diversity in the insect midguts; however, the curves for eight samples were still far from reaching the asymptote (Figure S2). The taxonomic classification rate varied from 86.7% to 99.8% (median 93.4%) at the genus level and from 69.3 to 93.2% (median 81.9%) at the species level (Table 1). In the whole datasheet, we identified 2519 different bacterial genera and 11,641 species in the range of 362–1746 and 1090–5324 for individual samples, respectively. However, since the base calling in the nanopore sequencing is still prone to errors, the results indicating the number of genera and species must be considered with caution. The bacterial genera present in all samples were considered as a core microbiome. This included 83 genera, constituting 3.3% of all identified genera. The most abundant were Bacillus, Microbacterium, Paenibacillus, Turicibacter, Bacteroides, Anaerotruncus, and Bradyrhizobium, whose average relative abundance (RA) exceeded 2%; however, a large inter-individual variation in the proportion of these genera was detected (Tables S2 and S3. The total overall relative abundance of all the core genera in the T0 and EPN-resistant groups of insects constituted 59.6% and 56.8% of the whole midgut community, respectively.

3.3. Bacterial Community Composition in the Midgut of the T0 Group of A. solstitiale Larvae

The analysis of the RA of the microbiota composition in the T0 group of larvae showed that the midgut was inhabited by bacteria belonging mainly to the phyla Firmicutes (48.9%), Proteobacteria (22.7%), Actinobacteria (14.6%), and Bacteroidetes (7.5%). Other detected phyla were present at RA < 1% (Figure 1, Table S4). The phylum Firmicutes included mainly Clostridia, Bacilli, and Erysipelotrichia (28.8%, 11.3%, and 7.0%, respectively), while Proteobacteria were represented by α- and γ-Proteobacteria (11.2 and 7.8%, respectively). The class Clostridia was detected to have the highest RA in individual samples of the T0 group—up to 58.3%. The classes Actinobacteria and Bacteroidia were the widespread components in the T0 group samples as well, with RAs of 10.8% and 7.3%, respectively. Other detected classes displayed an RA of <2% (Table S4, Supplementary Materials).
The analysis of the T0 larval midgut community at the family level revealed 8 families: Ruminococcaceae, Erysipelotrichaceae, Lachnospiraceae, Bacillaceae, Bacteroidaceae, Bradyrhizobiaceae, Propionibacteriaceae, and Hungateiclostridiaceae with RA > 2% (Figure 1, Table S5, Supplementary Materials). The most abundant genera in the midgut of the T0 larvae were Turicibacter, Bacteroides, Bacillus, Bradyrhizobium, Anaerotruncus, and Cutibacterium (Table S3), while Turicibacter sanguinis, Anaerotruncus rubiinfantis, and Cutibacterium acnes were the most abundant species (RA > 2%) (Table S6, Supplementary Materials).
Figure 1. Composition of midgut bacteria of A. solstitiale larvae: (A) phylum level, (B) family level.
Figure 1. Composition of midgut bacteria of A. solstitiale larvae: (A) phylum level, (B) family level.
Ijerph 19 03480 g001

3.4. Bacterial Community Composition in the Midgut of the EPN-Resistant Group of A. solstitiale Larvae

We did not notice important differences between the microbiota of insects exposed to certain nematode species, either in the α-diversity indices or in the PCoA and ANOSIM analyses. Therefore, we included all the samples exposed to EPN to one group “EPN-resistant”. In the midgut of the EPN-resistant group, similar to the T0 group, Firmicutes (44.0%), Actinobacteria (23.5%), and Proteobacteria (20.0%) were the dominated phyla; however, Bacteroidetes, Tenericutes (also known as Mycoplasmatota), and Synergistetes were detected with relatively high abundance as well (5.0%, 3.4%, and 1.6%, respectively). Other phyla were present at RA < 1% (Figure 1, Table S4). The most abundant classes of Firmicutes were Bacilli (23.4%) and Clostridia (18.2%), and the Bacilli class was detected with the highest RA in the individual samples of the EPN-resistant group—up to 87.1%. Other dominant phyla were Actinobacteria, α- and γ Proteobacteria, Bacteroidia, and Mollicutes (22.2%, 13.2%, 5.5%, 4.2%, and 3.4%, respectively). Other classes were detected with an RA < 2% (Table S4, Supplementary Materials).
The analysis of the midgut community of EPN-resistant larvae at the family level revealed more families, genera, and species, compared to the T0 group, with relatively high abundance. Families with an RA of >2% were represented by Bacillaceae, Microbacteriaceae, Ruminococcaceae, Lachnospiraceae, Paenibacillaceae, Rhizobiaceae, Enterobacteriaceae, Promicromonosporaceae, Spiroplasmataceae, Bradyrhizobiaceae, Mycobacteriaceae, Mycobacteriaceae, and Bacteroidaceae (Figure 1, Table S5). The most abundant genera in the midgut of the EPN-resistant larvae were Bacillus, Microbacterium, Paenibacillus, Spiroplasma, Enterobacter, Anaerotruncus, Rhizobium, Cellulosimicrobium, Virgibacillus, and Bacteroides (RA > 2) (Table S3). At the species level, the analysis revealed 4 species with an RA of >2%, i.e., Bacillus pumilus, Spiroplasma ixodetis, Microbacterium kyungheense, and Anaerotruncus rubiinfantis (Table S6, Supplementary Materials).

3.5. Comparison of the Midgut Microbiota between the T0 and the EPN-Resistant Groups of Insects

Although the mean value of the number of the observed genera in the T0 group of A. solstitiale was higher (1055) than in the EPN-resistant group (865), the difference was not significant (p-value 0.064). By contrast, both the community Evenness (E) and Shannon–Wiener (H’) diversity indices were higher in the T0 group of samples than in the group of EPN-resistant individuals (8.79 versus 5.71, p-value 0.009, and 4.36 versus 3.53, p-value 0.013, respectively). The median value of the Simpson diversity index (D-1) did not differ between the groups (8.21 and 6.29, p-value 0.059) (Table 2).
PCoA revealed a clear separation of 10 T0 samples from the rest of the specimens, indicating a different structure of the microbiota in the guts of insects subjected to the nematode exposure, prolonged incubation, and changed diet (Figure 2). The ANOSIM results fully supported this assumption (R 0.1781, p-value 0.006). The results of the same test performed on the T0 group only showed that the community structure in the midgut of the T0 group differed significantly, also given the location of the sampling site (R 0.409, p-value 0.016). However, the developmental stage did not exert such an effect (R 0.032, p-value 0.322).
The calculated fold changes revealed that the EPN-resistant larvae were enriched in 72 bacterial genera and were depleted in 96 genera. The most remarkable changes are presented in Figure 3. Beduini, Stenotrophobacter, Neisseria, unclassified Rickettsieae, Pseudomaricurvus, Shewanella, Gemmata, Sunxiuqinia, Ignavibacterium, and Nordella had significantly lower representation over the time of the experiment, i.e., comparing the group of the EPN-resistant specimens with the T0 group. On the other hand, Cellulosimicrobium, Diaminobutyricibacter, Ochrobactrum, Laspinema, Ancylobacter, Xylanimicrobium, Rathayibacter, Plantibacter, Kaistia, and Propionicimonas had much higher representation in the EPN-resistant larvae.
The microbiome networks that were created for the 100 most abundant genera of the two studied groups of A. solstitiale larvae differed substantially. The microbiome network of the EPN-resistant insects consisted of 3 clusters, while 2 clusters were detected in the T0 group (Table 3). The most important keystone genera were Sphingomonas, Colidextribacter, Cloacibacillus, Bradyrhizobium, and Nocardioides in the T0 group of larvae (Figure S3a) and Lachnoclostridium, Bosea, Mesorhizobium, Mycolicibacterium, and Devosia in the EPN-resistant group (Figure S3b, Supplementary Materials).

4. Discussion

4.1. Natural Midgut Microbiota of A. solstitiale Larvae

Several studies have shown that the digestive tract of scarab larvae contains a rich microbial community that may participate in various processes, especially digestion, and influence the immune system of the insect host [26]. However, little research has focused on the bacterial microbiota of the midgut, as this highly alkaline and oxidative compartment was considered to contain a less diverse bacterial community than in the hindgut, which is specialized for anaerobic fermentation [27,28]. In contrast, our data document a rich and diverse bacterial community associated with the midgut of A. solstitiale larvae (Coleoptera: Scarabaeidae). The considerable α-diversity indices of the midgut community were also observed in M. melolontha, Melolontha hippocastani, and Pachnoda spp., i.e., other Scarabaeidae members [16,29,30]. However, since different molecular techniques were used, a direct comparison of the α-diversity values seems to be improper.
In this study, the phyla Firmicutes, Proteobacteria, Actinobacteria, and Bacteroidetes were predominant in the analyzed midgut samples. The association of the midgut microbiota with mainly these four phyla has been reported for other rhizophagous larvae, e.g., Melolontha spp. [16,29]. However, while the phylum Proteobacteria dominated in the M. hippocastani midgut, Firmicutes representatives were dominant in M. melolontha and A. solstitiale. This may probably be explained by the remote location of the site of M. hippocastani sampling, compared with the overlapping locations in the case of the latter ones. Ruminococcaceae and Lachnospiracee were the most abundant families in the midgut of A. solstitiale and M. melolontha, respectively, unlike M. hippocastani, which was dominated by Enterobacteriaceae. The latter bacterial family was also widespread in the whole gut of other wild coleopterans, such as Nicrophorus vespiloides and Argilus planipennis [31,32].
The developmental stage, diet, host environment, and host taxonomy were found to be key factors influencing the structure of gut bacterial communities [33,34]. Research demonstrated that the gut community of rhizophagous M. hippocastani larvae was more complex than the community of grazing adults; however, they had a subset of bacteria in common [29,32]. The differences in the microbiota of scarab larvae, observed at different instars fed an identical diet, were more discrete, and they seemed to increase proportionally to the larval age/developmental stage [29,35]. In our study, the gut bacterial communities from the L2 and L3 larvae were not different, which is in agreement with the results of our previous studies on the midgut microbiota of M. melolontha [16]; however, larva development-driven changes in the midgut microbiota composition might have been obscured by the high individual variation observed. On the other hand, we noticed a difference in the community structure of the A. solstitiale larvae between the sampling locations, which is probably linked with the diet. This agrees with earlier studies describing the bacterial gut community structure in scarab beetles, which was strictly influenced by diet; however, it was clearly distinct from the microbiota of ingested food [29,36,37]. Additionally, urbanization as an environmental factor may exert an impact on the microbial-insect host association in the soil [38,39].
The assumption that the diet and the environment habitat are the determinants of gut microbial assemblages in the root-feeding scarab larvae studied can be also supported by the large intraspecies differences observed in the gut microbiota composition. In our study, the difference in the gut bacterial composition between the samples was apparent, especially at lower taxonomic levels. In many beetles, gut bacterial communities vary among individuals within a species and appear to consist largely of bacteria that are not specifically adapted to living in the guts of their host species. Large intraspecific differences in the gut community structure were also reported for the larvae of other scarab beetles, such as the rhizophagous M. melolontha [16,28] and Holotrichia parallela [35], or the humivorous Pachnoda spp. [30].
The presence of both constant and variable components in the gut bacterial community has been shown in our study on A. solstitiale larvae. A total of 83 different bacterial genera, which constituted 3.3% of all genera detected, were common to all the midgut samples analyzed; however, this number may be underestimated, considering the possible inflation of diversity derived from nanopore sequencing data. Many genera found in the core bacterial microbiota of the A. solstitiale midgut were shared with those of M. melolontha, such as the most abundant Turicibacter, Bacteroides, and Anaerotruncus, among others [16]. This implies that these bacteria may play a key role in various digestive and protective processes in insects. For example, Turicibacter sanguis, abundantly detected in the insect gut microbiota, were formerly described as mainly fermentative bacteria, but the latest research suggests that it plays an important role in the promotion of insect fitness by linking with the host’s serotonergic system [40]. It has been suggested that the inconstant gut community in rhizogenous scarab beetle larvae that is reported in many insects, such as the Pachysoma dung beetle or Drosophila spp., fits into the “minimal core” model initially described in mammals [36,41,42]. This model assumes weak partner fidelity and/or functional equivalence in diverse bacteria. On the other hand, we found a relatively high abundance of bacterial genera assigned to the core midgut microbiota of A. solstitile larvae (up to 60%), which may argue against favoring the minimal core microbiome concept. Nevertheless, since individual gut microbiota members can have an impact on host physiology, the ecological factors that determine natural variation in gut bacteria and the effect of this phenomenon on the fitness of their insect host are important issues for future study.

4.2. Alterations in the Gut Microbial Composition Associated with EPN-Resistance

It is known that EPN can infect a wide range of insects. However, this issue was frequently studied on lepidopteran caterpillars. A study on coleopteran T. molitor larvae confirmed their susceptibility to Steinernema spp. [43]. On the other hand, coleopteran larvae such as A. solstitiale or M. melolontha are relatively resistant to both Heterorhabditis and Steinernema nematodes [13,16,44]. This is consistent with the results of the survival assay, where we observed a high rate of larval survival following Steinerema and Heterorhabditis spp. treatment. Taking the impact of different EPN species exposure on the midgut microbiota community into account, we did not observe any significant differences. This is consistent with the findings reported by Cambon et al., who showed that the identity of three Steinernema spp. used for the infestation of insect larvae had no effect on the composition of the gut bacterial community in the early stage of infection [43]. However, it is also possible that we were not able to detect the discrete changes due to the high intra-individual variability of microbiota compositions.
Insects facing EPN associated with mutualistic bacteria co-evolved with them for millions of years, allowing the development of complex resistance mechanisms [5]. Several studies have reported that the composition of insect host microbiota is closely linked to the course of nematode parasitism [45,46]. Whether these resistance mechanisms, in addition to the activation of the immune response, avoidance behaviors, or morphological barriers, might include a selection of specific protective-acting microbes is still an open question. In addition to providing nutrients and improving the immune competence of the host against nematode infection, selection for resistance to nematodes may favor gut symbionts that act more directly, e.g., through the degradation of toxins released by nematodes and their native microbiota, the synthesis of a wide range of bioactive compounds that inhibit pathogen growth, or the alteration of the intestinal environment to limit successful EPN establishment. For example, certain bacterial endosymbionts of insects, Wolbachia and Spiroplasma spp., were found to be able to impede the growth of Brugia and Howardula spp. nematodes, which infect mosquitos and fruit flies, respectively [45,46]. Insect endosymbionts with antimicrobial activity, which are widespread in beetles, can also play an important role in the regulation of insect host–nematode interaction [47,48]. Previous studies showed that some members of the midgut microbiota of A. solstitiale and M. melolontha larvae, belonging to the Pseudomonas, Serratia, Acinetobacter, Citrobacter, Bacillus, Enterococcu, and Paenibacillus genera, are able to inhibit the growth of EPN symbiotic bacteria in vitro, implying that such activity can be one of the approaches to limit EPN infection in the natural world [13,44]. In all samples examined in this study, bacteria that are recognized as antagonistic to EPN symbionts were identified in the analyzed reads, reaching above the 0.001% RA thresholds (data not shown). This indicates that they are rare but widespread taxa, which is consistent with an equivalent study on M. melolontha [16].
The studied midgut bacterial communities showed considerable α-diversity in both analyzed groups of larvae; nevertheless, it was lower in the EPN-resistant group of larvae compared to the T0 group. These results are partly consistent with the study on M. melolontha larvae, where the α-diversity of midgut bacteria was comparably high, but we did not observe any significant difference between EPN-resistant and T0 individuals [16]. However, several studies on insects, including coleopteran soil-dwelling larvae, showed that resistance to entomopathogens is often correlated with simplified bacterial microbiota [49,50,51]. This finding may suggest that resistant A. solstitiale larvae living in the rich bacterial community of the soil are more selective in their internal filtering, which is aimed at avoiding species capable of becoming pathobionts in case of infection, thereby reducing the risk of secondary infections and septicemia, i.e., the major killing mechanism of EPN. Previous evidence has revealed that certain members of the insect gut community can facilitate Steinernema nematode pathogenesis, thus shaping the susceptibility of the host to infection in Manduca sexta and Tenebrio molitor larvae [43,52].
This study demonstrates that the composition of midgut bacteria is significantly different in the group of EPN-resistant insects, compared to the T0 group. Insect host resistance to parasites, bacteria, or toxins that are linked to an altered intestinal microbiota was also evidenced in several other studies [49,50,51,53,54]. Bacterial genera whose abundance is altered in A. solstitiale are usually rare members of the midgut community and belong to diverse phyla. Considering the most important changes, Actinobacteria were represented by a majority of genera that increased their abundance in the group with the EPN-resistant larvae compared to the T0 group, whereas the abundance of a majority of Proteobacteria decreased. In addition to the determination of differences in the composition of bacterial taxa in the microbiota, we assessed the keystone genera in both analyzed groups, which differed from each other. A similar study on M. melolontha bacterial microbiota also showed that the gut microbiota of the EPN-resistant group differed significantly in terms of the abundance of many genera, compared to the T0 group [16]; however, the different genera and the keystone genera differed from those shown in the present study. Nevertheless, the obtained results do not allow us to conclude whether the observed midgut bacterial microbiota alteration is a cause of resistance or rather a consequence of the EPN exposure.
In the PCoA analysis, the separation of the T0 group may rule out the explanation that specific bacteria were present in the resistant individuals before the exposure since, in this case, T0 individuals (which were not exposed to EPN) would likely be distributed among resistant ones. Therefore, the changes in the microbiota composition observed in the EPN-resistant larvae may be, at least partially, a consequence of EPN exposure, while the observed resistance of larvae to EPN is based on other mechanisms, probably a genetically determined immune response. This is in agreement with the study conducted by Tetreau et al., (2018), who concluded on the basis of a similar pattern of NMDS (Non-Metric Multidimensional Scaling) analysis that the modified microbiota associated with a high tolerance of Bacillus thuringienis was rather a consequence of the exposure of the host to this pathogen [50]. Similar to their study, our PCoA analysis displayed a higher dispersion among the resistant probes compared to the T0 group, supporting the concept of the EPN effect. In addition, the impact of the extended incubation and diet change in the case of the EPN-resistant individuals on the gut microbiota composition, leading to the spreading of these samples in PCoA analysis, could not be excluded. As previously shown, the gut microbiota of laboratory-reared insects fed with the simple diet is different (usually less diverse) from that in insects living in the wild [43].
On the other hand, in the PCoA analysis, we can observe three individuals from the T0 group that overlap with resistant insects, which suggests that they have changed microbiota and may represent the EPN-resistant phenotype, which is not connected with EPN exposure. A similar result was obtained in the experiments on the M. melolontha midgut [16]. In fact, the relationship between the gut microbiota the resistance to EPN is probably bidirectional. Selection for resistance to EPN can add pressure affecting the gut microbiota due to the activation of genes controlling bacterial overgrowth or the host’s immune response. Constantly increased immunity toward potential pathogens was detected in the Bacillus thuringiensis-resistant caterpillar Galleria mellonella [49].
Taken together, it is most probable that the changes observed in the midgut composition, apart from the extended incubation of larvae, are a result of the EPN treatment. It has been determined that EPN IJs highly transcribe genes coding for venom proteins already in the host-seeking stage, to be ready for their activation in the insect host; however, there are no data on the effects of these products on the host microbiota [55]. It is also known that other parasitic nematodes exhibit broad antimicrobial activity by secreting antibacterial factors such as lysozymes, cecropins, or C-type lectins; therefore, they can limit the growth of certain microbes and enrich the gut environment with nematode-protective microbiota members, supporting the suppression of host immunity and facilitating infection [45,56,57,58]. In addition, EPN mutualists are able to colonize the midgut, substantially affecting its microbiota, since they not only utilize gut nutrients but are also able to secrete a broad spectrum of antibiotics [59,60]. In turn, the role of non-mutualistic bacteria, isolated from different strains of Steinernema or Heterorhabditis nematodes, in their life cycle is unclear, but there is evidence that different secretively active bacteria can enter the insect gut with nematodes, such as Alcaligenes or Paenibacillus spp. [43,61].
Modification of the gut microbiome can also be an effect of tissue damage and changes in gut physiology, especially through the stimulation of the host’s innate immune system by the parasite, leading to increasing antimicrobial peptide production, which can keep gut microbiota under control [62,63]. Transcriptome analysis of Holotrichia pararella and Drosophila melanogaster treated with Steinernema and Heterorhabditis nematodes, respectively, indicated that, in the early stage of EPN infection, the activated genes were mostly associated with the immune system, especially the humoral response, leading to the secretion of antimicrobial peptides in gut epithelia [64,65]. Research has demonstrated that while Gram-positive bacteria are more sensitive to antibacterial defensins, antibacterial cecropins secreted into the gut lumen are active primarily against Gram-negative bacteria, such as Xenorhabdus and Photorhabdus bacteria [66,67,68]. This may support the suggestion that the observed decrease in the abundance of Proteobacteria is related to the insect immune response to the EPN mutualists.
On the other hand, direct toxin-mediated protection against the entomopathogens provided by defensive gut bacteria can not only mitigate EPN activity but also strongly affect the host and inhibit the members of its microbiota, particularly if the mechanisms are not specific, e.g., the production of bacteriocins or reactive oxygen species [9,69].

5. Conclusions

Despite the advent of increasing characterization of the insect pest microbiome, our understanding of the role of gut microbiota in determining infection outcomes is still insufficient. The application of more alternative tools to study the association between host microbiota and EPN calls for a determination of the role of specific gut residents, which is necessary for the elucidation of EPN resistance mechanisms. Our data confirm the high potential of 16S rRNA gene nanopore sequencing for profiling insect gut microbiota. Using this approach, we shed light on the midgut microbiota in wild coleopteran A. solstitiale larvae, highlighting compositional shifts in EPN infection-resistant individuals compared to T0 ones. Nevertheless, this work has some limitations, particularly with respect to the error-prone nanopore sequencing, inflating the biodiversity of the samples. Another limitation is the fact that we studied T0 samples (prior to the EPN exposure) as a point of reference; therefore, the observed effect can be biased by the lack of the same treatment of individuals of this group, compared to the EPN-resistant ones. Our study, however, is preliminary and we believe that the results are still relevant. Further studies on gut microbiota mediating the effect of host resistance to EPN are necessary, primarily to check whether insect resistance is dependent on nematode treatment. The high individual variation of the bacterial microbiota composition observed in this study may determine the possible effect of the microbiota on larval resistance to EPN, which is not uniform among individuals. Nevertheless, insects and their nematode parasites can be an excellent system for deciphering parasitic nematode strategies that result in pathogenic outcomes. Given the significance of scarab beetle pests, elucidation of insect-microbiota-pathogen interactions could potentially lead to the development of novel methods of pest biocontrol.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph19063480/s1, Table S1: List of sampling sites, Table S2: Proportions of the top twenty most abundant 100% prevailed bacterial genera detected in the midgut of A. solstitiale larvae, Table S3: Composition of midgut-associated bacteria in the A. solstitiale larvae at the genus level, Table S4: Composition of midgut-associated bacteria in the A. solstitiale larvae at the phylum and class level, Table S5: Composition of midgut-associated bacteria in the A. solstitiale larvae at the family level, Table S6: Composition of midgut-associated bacteria in the A. solstitiale larvae at the species level, Figure S1: Steps of sample collection and preparation—schematic diagram, Figure S2: Rarefaction curves from nanopore sequencing of the 16S rRNA gene for all samples, Figure S3: Microbiome networks of A. solstitiale midgut bacterial community. A—T0 group, B—EPN-resistant group.

Author Contributions

Conceptualization, M.S., E.S. and A.W.; data curation, A.K., M.S. and E.S.; formal analysis, A.K., M.S. and E.S.; funding acquisition, M.S., A.W. and E.S.; investigation, M.S. and E.S.; methodology, M.S., E.S., A.K. and W.K.; project administration, M.S., A.W. and E.S.; resources, M.S., M.L. and M.E.J.; supervision, E.S. and M.S., validation, A.K. and E.S.; visualization, A.K. and E.S.; writing—original draft, E.S.; writing—review and editing, E.S., A.K., M.S., W.K. and A.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the National Science Centre, Poland (grant 2016/21/B/NZ9/01865) and The John Paul II Catholic University of Lublin (grant 1/6-20-20-10-0603-0002-0503).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The authors confirm that the data supporting the findings of this study are available within the article and its Supplementary materials. Raw data that support the findings have been deposited in Sequenced Read Archive (SRA) under accession number PRJNA 787447.

Acknowledgments

The authors gratefully acknowledge the use of the services and facilities of the Center for Interdisciplinary Research of The John Paul II Catholic University of Lublin, Poland, co-funded by the EU from the European Regional Development Fund in the frame of the Operational Program Development of Eastern Poland 2007–2013 (POPW.01.03.00-06-003/09-00).

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ANOSIManalysis of similarities
BLASTbasic local alignment Search Tool
EPNentomopathogenic nematodes
IJsinfective juvenile
NCBINational Center for Biotechnology Information
OTUoperative taxonomic unit
PCoAprincipal component analysis
PCRpolymerase chain reaction
RAaverage relative abundance

References

  1. Jackson, T.A.; Klein, M.G. Scarabs as pests: Continuing problem. Coleopt. Bull. 2006, 60, 102–119. [Google Scholar] [CrossRef]
  2. Harrison Jdu, G.; Wingfield, M.J. A taxonomic review of white grubs and leaf chafers (Coleoptera: Scarabaeidae: Melolonthinae) recorded from forestry and agricultural crops in Sub-Saharan Africa. Bull. Entomol. Res. 2016, 6, 141–153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Guo, W.; Yan, X.; Zhao, C.; Han, R. Efficacy of entomopathogenic Steinernema and Heterorhabditis nematodes against white grubs (Coleoptera: Scarabaeidae) in peanut fields. J. Econ. Entomol. 2013, 106, 1112–1117. [Google Scholar] [CrossRef]
  4. Patil, J.; Vijayumar, R.; Nagesh, M.; Holajjer, P. Biocontrol potential of entomopathogenic Heterorhabditis and Steinernema nematodes against Phyllognathus dionysius Fabricius (Coleoptera: Scarabaeidae). Biol. Control 2019, 140, 104103. [Google Scholar] [CrossRef]
  5. Raja, R.K.; Arun, A.; Touray, M.; Gulsen, S.H.; Cimen, H.; Gulcu, B.; Hazir, C.; Aiswarya, D.; Ulug, D.; Cakmak, I.; et al. Antagonists and defense mechanisms of entomopathogenic nematodes and their mutualistic bacteria. Biol. Control 2020, 152, 104452. [Google Scholar] [CrossRef]
  6. Allahverdipour, H.; Karimi, J. Nematodes versus white grubs: Long but challenging association. Ann. Entomol. Soc. Am. 2021, 114, 448–458. [Google Scholar] [CrossRef]
  7. Eleftherianos, I.; Shokal, U.; Yadav, S.; Kenney, E.; Maldonado, T. Insect immunity to entomopathogenic nematodes and their mutualistic bacteria. Curr. Top. Microbiol. Immunol. 2017, 402, 123–156. [Google Scholar] [CrossRef] [PubMed]
  8. Rolhion, N.; Chassaing, B. When pathogenic bacteria meet the intestinal microbiota. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2016, 71, 20150504. [Google Scholar] [CrossRef] [Green Version]
  9. Oliver, K.M.; Perlman, S.J. Toxin-mediated protection against natural enemies by insect defensive symbionts. Adv. Insect Physiol. 2020, 58, 277–303. [Google Scholar] [CrossRef]
  10. Arora, A.K.; Douglas, A.E. Hype or opportunity? Using microbial symbionts in novel strategies for insect pest control. J. Insect Physiol. 2017, 103, 10–17. [Google Scholar] [CrossRef]
  11. Qadri, M.; Short, S.; Gast, K.; Hernandez, J.; Wong, A.C.-N. Microbiome innovation in agriculture: Development of microbial based tools for insect pest management. Front. Sustain. Food Syst. 2020, 4, 547751. [Google Scholar] [CrossRef]
  12. Stock, S.P. Diversity, biology and evolutionary relationships. In Nematode Pathogenesis of Insects and Other Pests: Ecology and Applied Technologies for Sustainable Plant and Crop Protection; Campos-Herrera, R., Ed.; Springer International Publishing: Neuchâtel, Switzerland, 2015; pp. 3–27. [Google Scholar]
  13. Skowronek, M.; Sajnaga, E.; Kazimierczak, W.; Lis, M.; Wiater, A. Screening and molecular identification of bacteria from the midgut of Amphimallon solstitiale larvae exhibiting antagonistic activity against bacterial symbionts of entomopathogenic nematodes. Int. J. Mol. Sci. 2021, 22, 12005. [Google Scholar] [CrossRef]
  14. Ziganshina, E.; Mohammed, W.; Shagimardanova, E.; Vankov, P.; Gogoleva, N.; Ziganshin, A. Fungal, bacterial, and archaeal diversity in the digestive tract of several beetle larvae (Coleoptera). Biomed Res. Int. 2018, 2018, 6765438. [Google Scholar] [CrossRef]
  15. Alonso-Pernas, P.; Arias-Cordero, E.; Novoselov, A.; Ebert, C.; Rybak, J.; Kaltenpoth, M.; Westermann, M.; Neugebauer, U.; Boland, W. Bacterial community and PHB-accumulating bacteria associated with the wall and specialized niches of the hindgut of the forest cockchafer (Melolontha hippocastani). Front. Microbiol. 2017, 28, 291. [Google Scholar] [CrossRef] [Green Version]
  16. Sajnaga, E.; Skowronek, M.; Kalwasińska, A.; Kazimierczak, W.; Ferenc, K.; Lis, M.; Wiater, A. Nanopore-sequencing characterization of the gut microbiota of Melolontha melolontha larvae: Contribution to protection against entomopathogenic nematodes? Pathogens 2021, 10, 396. [Google Scholar] [CrossRef]
  17. Koppenhöfer, A.M.; Fuzy, E. Effect of white grub developmental stage on susceptibility to entomopathogenic nematodes. J. Econ. Entomol. 2004, 97, 1842–1849. [Google Scholar] [CrossRef]
  18. Suzuki, M.T.; Giovannoni, S.J. Bias caused by template annealing in the amplification of mixtures of 16S rRNA genes by PCR. Appl. Environ. Microbiol. 1996, 62, 625–630. [Google Scholar] [CrossRef] [Green Version]
  19. Lane, D.J.; Pace, B.; Olsen, G.J.; Stahl, D.A.; Sogin, M.L.; Pace, N.R. Rapid determination of 16S ribosomal RNA sequences for phylogenetic analyses. Proc. Natl. Acad. Sci. USA 1985, 82, 6955–6959. [Google Scholar] [CrossRef] [Green Version]
  20. Hammer, Ø.; Harper, D.A.T.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. Paleontol. Electronica 2001, 4, 1–9. [Google Scholar]
  21. McMurdie, P.J.; Holmes, S. Waste not, want not: Why rarefying microbiome data is inadmissible. PLoS Comput. Biol. 2014, 10, 1003531. [Google Scholar] [CrossRef] [Green Version]
  22. Layeghifard, M.; Hwang, D.M.; Guttman, D.S. Constructing and analyzing microbiome networks in R. Methods Mol. Biol. 2018, 1849, 243–266. [Google Scholar] [CrossRef] [PubMed]
  23. Friedman, J.; Alm, E.J. Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 2012, 8, e1002687. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Brandes, U. A faster algorithm for betweenness centrality. J. Math. Sociol. 2001, 25, 163–177. [Google Scholar] [CrossRef]
  25. Pons, P.; Latapy, M. Computing communities in large networks using random walks. Lect. Notes Comput. Sci. 2005, 3733, 284–293. [Google Scholar] [CrossRef] [Green Version]
  26. Thiyonila, B.; Reneeta, N.P.; Kannan, M.; Shantkriti, S.; Krishnan., M. Dung beetle gut microbes: Diversity, metabolic and immunity related roles in host system. Int. J. Sci. Innovs. 2018, 4, 77–83. [Google Scholar] [CrossRef]
  27. Lemke, T.; Stingl, U.; Egert, M.; Friedrich, M.W.; Brune, A. Physicochemical conditions and microbial activities in the highly alkaline gut of the humus-feeding larva of Pachnoda ephippiata (Coleoptera: Scarabaeidae). Appl. Environ. Microbiol. 2003, 69, 6650–6658. [Google Scholar] [CrossRef] [Green Version]
  28. Egert, M.; Stingl, U.; Bruun, L.D.; Brune, A.; Friedrich, M.W.; Pommerenke, B. Structure and topology of microbial communities in the major gut compartments of Melolontha melolontha larvae (Coleoptera:Scarabaeidae). Appl. Environ. Microbiol. 2005, 71, 4556–4566. [Google Scholar] [CrossRef] [Green Version]
  29. Arias-Cordero, E.; Ping, L.; Reichwald, K.; Delb, H.; Platzer, M.; Boland, W. Comparative evaluation of the gut microbiota associated with the below- and above-ground life stages (larvae and beetles) of the forest cockchafer, Melolontha hippocastani. PLoS ONE 2012, 7, e51557. [Google Scholar] [CrossRef] [Green Version]
  30. Andert, J.; Marten, A.; Brandl, R.; Brune, A. Inter- and intraspecific comparison of the bacterial assemblages in the hindgut of humivorous scarab beetle larvae (Pachnoda spp.). FEMS Microbiol. Ecol. 2010, 74, 439–449. [Google Scholar] [CrossRef]
  31. Wang, Y.; Rozen, D.E. Gut microbiota colonization and transmission in the burying beetle Nicrophorus vespilloides throughout development. Appl. Environ. Microbiol. 2017, 83, e03250–e03316. [Google Scholar] [CrossRef] [Green Version]
  32. Vasanthakumar, A.; Handelsman, J.; Schloss, P.D.; Bauer, L.S.; Raffa, K.F. Gut microbiota of an invasive subcortical beetle, Agrilus planipennis Fairmaire, across various life stages. Environ. Entomol. 2008, 37, 1344–1353. [Google Scholar] [CrossRef]
  33. Colman, D.R.; Toolson, E.C.; Takacs-Vesbach, C.D. Do diet and taxonomy influence insect gut bacterial communities? Mol. Ecol. 2012, 21, 5124–5137. [Google Scholar] [CrossRef]
  34. 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. Insect Gut Bacterial Diversity Determined by Environmental Habitat, Diet, Developmental Stage, and Phylogeny of Host. Appl. Environ. Microbiol. 2014, 80, 5254–5264. [Google Scholar] [CrossRef] [Green Version]
  35. Huang, S.; Zhang, H. The impact of environmental heterogeneity and life stage on the hindgut microbiota of Holotrichia parallela larvae (Coleoptera: Scarabaeidae). PLoS ONE 2013, 8, e57169. [Google Scholar] [CrossRef]
  36. Franzini, P.Z.; Ramond, J.B.; Scholtz, C.H.; Sole, C.L.; Ronca, S.; Cowan, D.A. The gut microbiomes of two Pachysoma MacLeay desert dung beetle species (Coleoptera: Scarabaeidae: Scarabaeinae) feeding on different diets. PLoS ONE 2016, 11, e0161118. [Google Scholar] [CrossRef]
  37. Egert, M.; Wagner, B.; Lemke, T.; Brune, A.; Friedrich, M.W. Microbial community structure in midgut and hindgut of the humus-feeding larva of Pachnoda ephippiata (Coleoptera: Scarabaeidae). Appl. Environ. Microbiol. 2003, 69, 6659–6668. [Google Scholar] [CrossRef] [Green Version]
  38. Bosmans, L.; Pozo, M.I.; Verreth, C.; Crauwels, S.; Wilberts, L.; Sobhy, S.; Wäckers, F.; Jacquemyn, H.; Lievens, B. Habitat-specific variation in gut microbial communities and pathogen prevalence in bumblebee queens (Bombus terrestris). PLoS ONE 2018, 13, e0204612. [Google Scholar] [CrossRef] [Green Version]
  39. Park, R.; Dzialo, M.C.; Spaepen, S.; Nsabimana, D.; Gielens, K.; Devriese, H.; Crauwels, S.; Tito, R.Y.; Raes, J.; Lievens, B.; et al. Microbial communities of the house fly Musca domestica vary with geographical location and habitat. Microbiome 2019, 7, 147. [Google Scholar] [CrossRef]
  40. Fung, T.C.; Vuong, H.E.; Luna, C.D.G.; Pronovost, G.N.; Aleksandrova, A.; Riley, N.G.; Vavilina, A.; McGinn, J.; Rendon, T.; Forrest, L.R.; et al. Intestinal serotonin and fluoxetine exposure modulate bacterial colonization in the gut. Nat. Microbiol. 2019, 4, 2064–2073. [Google Scholar] [CrossRef]
  41. Hamady, M.; Knight, R. Microbial community profiling for human microbiome projects: Tools, techniques, and challenges. Genome Res. 2009, 19, 1141–1152. [Google Scholar] [CrossRef] [Green Version]
  42. Wong, A.N.; Chaston, J.; Douglas, A. The inconstant gut microbiota of Drosophila species revealed by 16S rRNA gene analysis. ISME J. 2013, 7, 1922–1932. [Google Scholar] [CrossRef] [Green Version]
  43. Cambon, M.C.; Lafont, P.; Frayssinet, M.; Lanois, A.; Ogier, J.C.; Pagès, S.; Parthuisot, N.; Ferdy, J.B.; Gaudriault, S. Bacterial community profile after the lethal infection of Steinernema-Xenorhabdus pairs into soil-reared Tenebrio molitor larvae. FEMS Microbiol. Ecol. 2020, 96, fiaa009. [Google Scholar] [CrossRef]
  44. Skowronek, M.; Sajnaga, E.; Pleszczyńska, M.; Kazimierczak, W.; Lis, M.; Wiater, A. Bacteria from the midgut of common cockchafer (Melolontha melolontha L.) larvae exhibiting antagonistic activity against bacterial symbionts of entomopathogenic nematodes: Isolation and molecular identification. Int. J. Mol. Sci. 2020, 21, 580. [Google Scholar] [CrossRef] [Green Version]
  45. Topalović, O.; Vestergård, M. Can microorganisms assist the survival and parasitism of plant-parasitic nematodes? Trends Parasitol. 2021, 37, 947–958. [Google Scholar] [CrossRef]
  46. Eleftherianos, I.; Yadav, S.; Kenney, E.; Cooper, D.; Ozakman, Y.; Patrnogic, J. Role of endosymbionts in insect-parasitic nematode interactions. Trends Parasitol. 2018, 34, 430–444. [Google Scholar] [CrossRef]
  47. Heise, P.; Liu, Y.; Degenkolb, T.; Vogel, H.; Schäberle, T.F.; Vilcinskas, A. Antibiotic-producing beneficial bacteria in the gut of the burying beetle Nicrophorus vespilloides. Front. Microbiol. 2019, 10, 1178. [Google Scholar] [CrossRef] [Green Version]
  48. Flórez, L.V.; Scherlach, K.; Miller, I.J.; Rodrigues, A.; Kwan, J.C.; Hertweck, C.; Kaltenpoth, M. An antifungal polyketide associated with horizontally acquired genes supports symbiont-mediated defense in Lagria villosa beetles. Nat. Commun. 2018, 9, 2478. [Google Scholar] [CrossRef] [Green Version]
  49. Dubovskiy, I.M.; Grizanova, E.V.; Whitten, M.M.; Mukherjee, K.; Greig, C.; Alikina, T.; Kabilov, M.; Vilcinskas, A.; Glupov, V.V.; Butt, T.M. Immuno-physiological adaptations confer wax moth Galleria mellonella resistance to Bacillus thuringiensis. Virulence 2016, 7, 860–870. [Google Scholar] [CrossRef] [Green Version]
  50. Tetreau, G.; Grizard, S.; Patil, C.D.; Tran, F.H.; Tran Van, V.; Stalinski, R.; Laporte, F.; Mavingui, P.; Després, L.; Valiente Moro, C. Bacterial microbiota of Aedes aegypti mosquito larvae is altered by intoxication with Bacillus thuringiensis israelensis. Parasit. Vectors 2018, 11, 121. [Google Scholar] [CrossRef] [Green Version]
  51. Paddock, K.J.; Pereira, A.E.; Finke, D.L.; Ericsson, A.C.; Hibbard, B.E.; Shelby, K.S. Host resistance to Bacillus thuringiensis is linked to altered bacterial community within a specialist insect herbivore. Mol. Ecol. 2021, 30, 1–16. [Google Scholar] [CrossRef]
  52. Singh, S.; Reese, J.M.; Casanova-Torres, A.M.; Goodrich-Blair, H.; Forst, S. Microbial population dynamics in the hemolymph of Manduca sexta infected with Xenorhabdus nematophila and the entomopathogenic nematode Steinernema carpocapsae. Appl. Environ. Microbiol. 2014, 80, 4277–4285. [Google Scholar] [CrossRef] [Green Version]
  53. Xia, X.; Zheng, D.; Zhong, H.; Qin, B.; Gurr, G.M.; Vasseur, L.; Lin, H.; Bai, J.; He, W.; You, M. DNA sequencing reveals the midgut microbiota of diamondback moth, Plutella xylostella (L.) and a possible relationship with insecticide resistance. PLoS ONE 2013, 8, e68852. [Google Scholar] [CrossRef] [PubMed]
  54. Boissière, A.; Tchioffo, M.T.; Bachar, D.; Abate, L.; Marie, A.; Nsango, S.E.; Shahbazkia, H.R.; Awono-Ambene, P.H.; Levashina, E.; Christen, R.; et al. Midgut Microbiota of the Malaria Mosquito Vector Anopheles gambiae and Interactions with Plasmodium falciparum Infection. PLoS Pathog. 2012, 8, e1002742. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Lu, D.; Macchietto, M.; Chang, D.; Barros, M.M.; Baldwin, J.; Mortazavi, A.; Dillman, A.R. Activated entomopathogenic nematode infective juveniles release lethal venom proteins. PLoS Pathog. 2017, 13, e1006302. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Dierking, K.; Yang, W.; Schulenburg, H. Antimicrobial effectors in the nematode Caenorhabditis elegans: An outgroup to the Arthropoda. Philo. Trans. R. Soc. Lond. B Biol. Sci. 2016, 371, 20150299. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Midha, A.; Janek, K.; Niewienda, A.; Henklein, P.; Guenther, S.; Serra, D.O.; Schlosser, J.; Hengge, R.; Hartmann, S. The Intestinal roundworm Ascaris suum releases antimicrobial factors which interfere with bacterial growth and biofilm formation. Front. Cell Infect. Microbiol. 2018, 8, 271. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Rausch, S.; Midha, A.; Kuhring, M.; Affinass, N.; Radonic, A.; Kühl, A.A.; Bleich, A.; Renard, B.Y.; Hartmann, S. Parasitic nematodes exert antimicrobial activity and benefit from microbiota-driven support for host immune regulation. Front. Immunol. 2018, 9, 2282. [Google Scholar] [CrossRef] [Green Version]
  59. Silva, C.P.; Waterfield, N.R.; Daborn, P.J.; Dean, P.; Chilver, T.; Au, C.P.; Sharma, S.; Potter, U.; Reynolds, S.E.; French-Constant, R.H. Bacterial infection of a model insect: Photorhabdus luminescens and Manduca sexta. Cell Microbiol. 2002, 4, 329–339. [Google Scholar] [CrossRef] [Green Version]
  60. Carneiro, C.N.; DaMatta, R.A.; Samuels, R.I.; Silva, C.P. Effects of entomopathogenic bacterium Photorhabdus temperata infection on the intestinal microbiota of the sugarcane stalk borer Diatraea saccharalis (Lepidoptera: Crambidae). J. Invertebr. Pathol. 2008, 99, 87–91. [Google Scholar] [CrossRef]
  61. El-Borai, F.E.; Duncan, L.W.; Preston, J.F. Bionomics of a phoretic association between Paenibacillus sp. and the entomopathogenic nematode Steinernema diaprepesi. J. Nematol. 2005, 37, 18–25, Erratum in: J. Nematol. 2008, 40, 55. [Google Scholar]
  62. Login, F.H.; Balmand, S.; Vallier, A.; Vincent-Monégat, C.; Vigneron, A.; Weiss-Gayet, M.; Rochat, D.; Heddi, A. Antimicrobial peptides keep insect endosymbionts under control. Science 2011, 334, 362–365. [Google Scholar] [CrossRef]
  63. Dahan, D.; Preston, G.; Sealey, J.; King, K.C. Impacts of a novel defensive symbiosis on the nematode host microbiome. BMC Microbiol. 2020, 20, 159. [Google Scholar] [CrossRef]
  64. Li, E.; Qin, J.; Feng, H.; Li, J.; Li, X.; Nyamwasa, I.; Cao, Y.; Ruan, W.; Li, K.; Yin, J. Immune-related genes of the larval Holotrichia parallela in response to entomopathogenic nematodes Heterorhabditis beicherriana LF. BMC Genom. 2021, 22, 192. [Google Scholar] [CrossRef]
  65. Yadav, S.; Daugherty, S.; Shetty, A.C.; Eleftherianos, I. RNAseq analysis of the Drosophila response to the entomopathogenic nematode Steinernema. G3 2017, 7, 1955–1967. [Google Scholar] [CrossRef] [Green Version]
  66. Xu, X.X.; Zhang, Y.Q.; Freed, S.; Yu, J.; Gao, Y.F.; Wang, S.; Ouyang, L.N.; Ju, W.Y.; Jin, F.L. An anionic defensin from Plutella xylostella with potential activity against Bacillus thuringiensis. Bull. Entomol. Res. 2016, 106, 790–800. [Google Scholar] [CrossRef]
  67. Tonk, M.; Knorr, E.; Cabezas-Cruz, A.; Valdés, J.J.; Kollewe, C.; Vilcinskas, A. Tribolium castaneum defensins are primarily active against Gram-positive bacteria. J. Invertebr. Pathol. 2015, 132, 208–215. [Google Scholar] [CrossRef]
  68. Moore, A.J.; Beazley, W.D.; Bibby, M.C.; Devine, D.A. Antimicrobial activity of cecropins. J. Antimicrob Chemother. 1996, 37, 1077–1089. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Koch, H.; Schmid-Hempel, P. Socially transmitted gut microbiota protect bumble bees against an intestinal parasite. Proc. Natl. Acad. Sci. USA 2011, 108, 19288–19292. [Google Scholar] [CrossRef] [Green Version]
Figure 2. The principal coordinate analysis (PCoA) generated from the sequencing results of the 16S rRNA gene at the genus level. G—midgut samples of individuals from the T0 group (marked in blue); GA, GB, GC, GK, GM—midgut samples of individuals exposed to S. arenarium, S. bicornutum, S. carpocapse, S. kraussei, and S. megidis, respectively (marked in black). The numbers are an integral part of the sample name. All samples are further described in Table 1.
Figure 2. The principal coordinate analysis (PCoA) generated from the sequencing results of the 16S rRNA gene at the genus level. G—midgut samples of individuals from the T0 group (marked in blue); GA, GB, GC, GK, GM—midgut samples of individuals exposed to S. arenarium, S. bicornutum, S. carpocapse, S. kraussei, and S. megidis, respectively (marked in black). The numbers are an integral part of the sample name. All samples are further described in Table 1.
Ijerph 19 03480 g002
Figure 3. Average fold changes in the genera, with different abundances in the midgut of both the T0 and EPN-resistant groups of A. solstitiale larvae. (a) Negative changes, (b) positive changes (p < 0.05).
Figure 3. Average fold changes in the genera, with different abundances in the midgut of both the T0 and EPN-resistant groups of A. solstitiale larvae. (a) Negative changes, (b) positive changes (p < 0.05).
Ijerph 19 03480 g003
Table 1. Sampling sites, NGS sequencing statistics, and α-diversity metrics for the studied midgut bacterial microbiota.
Table 1. Sampling sites, NGS sequencing statistics, and α-diversity metrics for the studied midgut bacterial microbiota.
Sample NameNematode ExposureSampling Site *Developmental Stage of LarvaTotal Number of ReadsClassification Rate (%)Diversity Metrics
Genus LevelSpecies LevelSobsShannonSimpsonEvenness
G1NoFNL2252.22691.879.18904.630.980.68
G2NoFNL2175.05696.183.86463.490.930.54
G3NoFNL3191.58693.786.19674.380.970.64
G3RNoFNL3130.81089.185.79144.090.950.60
G4NoFNL3161.82495.991.811574.560.970.65
G5NoFNL2666.34591.086.89774.250.970.62
G6NoFNL2202.97486.781.89693.960.910.58
G7NoFNL2193.90089.483.48624.220.950.62
G8NoULL3139.28995.181.514645.320.990.73
G9NoULL383.10993.887.712654.420.920.62
G10NoULL3112.61993.075.313784.820.980.67
G11NoULL3194.58290.779.111824.330.950.61
G12NoULL397.43993.483.710514.250.950.61
GA1YesFNL2137.31092.786.96473.590.930.55
GA2YesFNL2107.18692.681.39874.600.980.67
GA3YesFNL217311997.693.26222.850.790.44
GB1YesFNL2146.97193.178.213605.010.980.69
GB2YesFNL2181.11993.386.79313.930.950.58
GB2RYesFNL2199.58296.077.02250.830.280.15
GB3YesFNL2217.01499.769.39263.260.870.48
GC1YesFNL2216.09193.278.111614.890.980.69
GC2YesFNL2175.53393.285.88514.040.960.60
GC3YesFNL2217.95192.181.912283.350.790.47
GK1YesFNL2158.66695.974.85212.370.760.38
GK2YesFNL2257.76696.374.68483.370.880.50
GM1YesFNL2129.20996.873.37973.330.910.50
GM2YesFNL2300.87995.181.38983.710.940.55
GM3YesFNL2244.16294.785.09703.870.930.56
* FN: forest nursery at 51°23′37.0″ N 22°29′44.0″ E; UL: urban lawn at 51°14′19.7″ N 22°29′58.7″ E.
Table 2. Comparison of diversity metrics between the T0 group of A. solstitiale larvae and EPN-resistant larvae.
Table 2. Comparison of diversity metrics between the T0 group of A. solstitiale larvae and EPN-resistant larvae.
Diversity MetricsT0EPN-ResistantStatisticStatistic-Valuep-Value
N = 13N = 15
Observed1055865t1.940.064
Shannon4.363.53t2.820.013
Simpson8.216.29Mann–Whitney U560.059
Evenness8.795.71Mann–Whitney U400.009
Table 3. Basic characteristics of microbiome networks in the T0 group of A. solstitiale larvae and EPN-resistant larvae.
Table 3. Basic characteristics of microbiome networks in the T0 group of A. solstitiale larvae and EPN-resistant larvae.
T0 GroupEPN-Resistant Group
(N = 13)(N = 15)
Keystone generaSphingomonas (1.00 *)Lachnoclostridium (1.0)
Colidextribacter (0.99)Bosea (0.99)
Cloacibacillus (0.98)Mesorhizobium (0.98)
Bradyrhizobium (0.96)Mycolicibacterium (0.96)
Nocardioides (0.95)Devosia
Number of clusters 23
Modularity0.0670.126
* hubbiness scores are given in parentheses.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sajnaga, E.; Skowronek, M.; Kalwasińska, A.; Kazimierczak, W.; Lis, M.; Jach, M.E.; Wiater, A. Comparative Nanopore Sequencing-Based Evaluation of the Midgut Microbiota of the Summer Chafer (Amphimallon solstitiale L.) Associated with Possible Resistance to Entomopathogenic Nematodes. Int. J. Environ. Res. Public Health 2022, 19, 3480. https://doi.org/10.3390/ijerph19063480

AMA Style

Sajnaga E, Skowronek M, Kalwasińska A, Kazimierczak W, Lis M, Jach ME, Wiater A. Comparative Nanopore Sequencing-Based Evaluation of the Midgut Microbiota of the Summer Chafer (Amphimallon solstitiale L.) Associated with Possible Resistance to Entomopathogenic Nematodes. International Journal of Environmental Research and Public Health. 2022; 19(6):3480. https://doi.org/10.3390/ijerph19063480

Chicago/Turabian Style

Sajnaga, Ewa, Marcin Skowronek, Agnieszka Kalwasińska, Waldemar Kazimierczak, Magdalena Lis, Monika Elżbieta Jach, and Adrian Wiater. 2022. "Comparative Nanopore Sequencing-Based Evaluation of the Midgut Microbiota of the Summer Chafer (Amphimallon solstitiale L.) Associated with Possible Resistance to Entomopathogenic Nematodes" International Journal of Environmental Research and Public Health 19, no. 6: 3480. https://doi.org/10.3390/ijerph19063480

APA Style

Sajnaga, E., Skowronek, M., Kalwasińska, A., Kazimierczak, W., Lis, M., Jach, M. E., & Wiater, A. (2022). Comparative Nanopore Sequencing-Based Evaluation of the Midgut Microbiota of the Summer Chafer (Amphimallon solstitiale L.) Associated with Possible Resistance to Entomopathogenic Nematodes. International Journal of Environmental Research and Public Health, 19(6), 3480. https://doi.org/10.3390/ijerph19063480

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop